岩性油气藏  2021, Vol. 33 Issue (3): 1-12       PDF    
×
坳陷湖盆斜坡区深水重力流水道地震响应及沉积特征——以松辽盆地LHP地区嫩江组一段为例
刘化清, 冯明, 郭精义, 潘树新, 李海亮, 洪忠, 梁苏娟, 刘彩燕, 徐云泽     
中国石油勘探开发研究院 西北分院, 兰州 730020
摘要: 松辽盆地湖盆深水区重力流沉积广泛发育,由于资料分辨率低的限制,以往对其时空演化、内部结构方面的认识不清。基于高精度三维地震资料,利用90°相位转换、地层切片、光照显示等技术,研究了LHP地区下白垩统嫩一段湖盆深水斜坡上发育的重力流水道的地震特征以及水道的内部结构。结果表明:①嫩一段沉积时期,本区坡降为0.5‰~1.0‰,发育早、晚2套共9条迥然不同的水道系统。早期水道E-I处于整个重力流水道系统的下游近末端,为非限制性水道,宽度小于40 m,广泛发育决口分流水道,对下伏早期沉积侵蚀弱。晚期水道A-D处于重力流水道系统的中段,为限制性水道,其宽度为50~370 m,在曲流环发育部位宽度可达500 m,这类水道的下切深度可达20 m,剖面形态为U或V型,局部发育2~5 m高的天然堤。发育早、晚2套不同的水道系统可能与构造持续沉降导致盆地斜坡区的坡降变化有关。②坡降与水道曲率、深度与宽度之间均呈良好的正相关关系。③水道中细砂岩、粉砂岩包裹在优质烃源岩当中,具有形成岩性油气藏的良好条件,是今后重要的勘探领域。
关键词: 湖盆深水斜坡    重力流水道    地震沉积学    地震地貌学    水道结构参数    嫩江组    松辽盆地    
Seismic reflection and sedimentary characteristics of deep-water gravity flow channels on the slope of lacustrine depression basin: First member of Nenjiang Formation in LHP area, Songliao Basin
LIU Huaqing, FENG Ming, GUO Jingyi, PAN Shuxin, LI Hailiang, HONG Zhong, LIANG Sujuan, LIU Caiyan, XU Yunze     
PetroChina Research Institute of Petroleum Exploration & Development-Northwest, Lanzhou 730020, China
Abstract: Multiple deep-water slope channels have been observed within the lacustrine Songliao Basin, northeastern China. To understandthe spatial-temporal evolution and internal structures of these channels, techniques including 90° phase rotation, stratal slicing, opacity stacks and illumination were used to study the seismic characteristics and internal structure of gravity flow channels developed on the deep-water slope of lacustrine basin of Lower Cretaceous Nen 1 member in LHP area. The results show that: (1)Nine channels within the 50 km2 LHP area were developed in Nen 1 member of the Late Cretaceous Nenjiang Formation. The depositional setting was a sublacustrine slope characterized by a gradient of 0.5‰ - 1.0‰.(2) Based on morphology, the channels were grouped into two channel systems based on size and degree of associated erosion: small meandering channels characterized by sinuosity of 1.03-1.17, which developed early and were commonly less than 40 m wide. Erosional confinement was minimal and distal avulsions were common. These channels can be described as relatively unconfined, lying down-slope of large gravity-flow systems.(3)Larger meandering channels with sinuosity of 1.04-1.45 and widths up to 500 m, were characterized by greater erosion with maximum erosional valley relief of 20 m. These channels were sufficiently large so that cross profiles could be evaluated, revealing U/V shaped section morphology. Avulsions were uncommon and levees were observed locally. These channels likely lie within the central part of associated gravity-flow systems.(4)Morphological statistics suggest that channel slope and sinuosity, depth and width are positively correlated.(5)Well data from neighboring area indicate that the channel fills were characterized by predominantly silt/fine sandstones with thickness of 10-20 m. These reservoir-prone deposits were bracketed in good quality source rocks making them favorable hydrocarbon exploration targets.(6)The Late Cretaceous global anoxic event possibly caused the abnormal climate(frequent flooding events)and then the frequent gravity flows in the Nen 1 member of lacustrine basin. The continuous subsidence increased the gradient of the paleo-lake and subsequently caused the geomorphology difference between the channels A-D and E-I.
Key words: sublacustrine slope    gravity flow channels    seismic sedimentology    seismic geomorphology    channel structure parameters    Nenjiang Formation    Songliao Basin    
0 引言

近年来,石油勘探家对松辽盆地下白垩统嫩江组一段重力流水道系统的宏观分布特征进行了识别和描述,在嫩一段发现了陆相盆地中规模最大、清晰完整的深水重力流水道-湖底扇系统[1-4]。如大庆长垣的水道系统延伸的直线距离可达70 km,宽度为600 m,平均宽度为300 m,而在齐家—古龙地区的规模更大,长度超过80 km,宽度为100~ 900 m。潘树新等[2]认为这些规模宏大的深水重力流水道体系由陆上河流直接入湖而形成,为异重流成因。这些水道被细砂岩、粉砂岩充填,其含油性好,是盆地中央深水区岩性油气藏勘探潜在的新领域[1-2]。对水道的时空演化、水道的内部结构及各结构参数(曲率、坡降、深度、宽度等)及其相关性研究相对较少,而深入研究这些问题有助于弄清水道系统的结构特征以及水道砂岩的非均质性,从而为油气勘探开发部署提供重要依据。

以LHP地区高精度三维地震资料为基础,应用地震沉积学、地震地貌学技术方法,识别湖盆深水斜坡的重力流水道及其伴生的水下天然堤的分布特征,研究水道时空演化规律、确定水道的宽、深、曲率等结构参数及其相关性,探讨水道形成与演化的控制因素,以期为油气勘探开发部署提供依据,并为其他陆相盆地深水水道沉积研究提供借鉴。

1 沉积背景与研究区概况

松辽盆地的基底为变质的古生界及前古生界,中新生界沉积盖层的发育经历了伸展断陷阶段、裂后热沉降阶段和构造反转阶段[1]。其中,裂后热沉降阶段发育的下白垩统泉头组(K1q)和上白垩统青山口组(K2qn)、姚家组(K2y)和嫩江组(K2n)是盆地的主要含油层系。姚家组—嫩江组发育一套完整的基准面升降旋回沉积,嫩江组二段底部发育了全盆地范围内稳定分布的油页岩,这标志着最大湖泛面的形成期,是层序划分与对比的标志层(T07)[5]。姚家组—嫩江组一段沉积期为湖盆扩张期,形成了大型退积型三角洲及重力流水道-末端扇沉积体系;嫩江组二段—五段沉积期为湖盆萎缩期,强制性水退形成大型湖泊与进积型三角洲沉积[5]。嫩五段沉积以后,盆地抬升并遭受剥蚀,与上覆四方台组(K2s)形成不整合接触关系[6]

松辽盆地嫩江组一段沉积时期,湖盆水体快速扩张。潘树新等[2]对齐家—古龙地区嫩一段上部(S1)沉积体系进行过深入研究,发现该区发育多支大型异重流成因的水道—湖底扇系统。本文选取该水道系统的中下游的LHP地区开展嫩一段上部(S0)水道的精细解剖[图 1(a)]。LHP三维地震资料精度高,采集面元为10 m×10 m,主频为50 Hz。工区内有钻井10余口,录井及测井资料比较齐全。

下载原图 图 1 松辽盆地嫩一段沉积体系[2](a)和地层综合柱状图(b) Fig. 1 Sedimentary system(a)and comprehensive stratigraphic column(b)of Nen 1 member in Songliao Basin

LHP地区嫩一段岩性主要为厚层黑色、灰黑色泥岩,夹薄层粉砂岩、泥质粉砂岩,埋深为1 500~ 1 700 m,厚度为110 m,其顶、底界对应的地震反射界面分别为T07和T1。嫩一段可进一步细分为上下2个部分:下部(S1油层组)形成于湖盆快速扩张期,以黑色泥岩、薄层泥质粉砂岩、粉砂岩互层为主,粉砂岩单层厚度最大可达5 m,底部以厚层泥岩与下伏姚家组富砂地层(S2+3)接触,形成地震强反射界面(T1);上部(S0油层组)以厚层泥岩为主,粉砂岩、泥质粉砂岩零星分布,其厚度仅为2 m[图 1(b)],形成于湖盆发育的鼎盛时期。S0地层与上覆嫩二段底部油页岩形成显明的地震波阻抗界面(T07)。

2 研究方法与技术 2.1 地震地貌成像处理

在地震地貌成像过程中,综合利用了地震沉积学[7-8]和地震地貌学[9]的研究思路和技术手段:①首先将地震数据体转换为具有岩性意义的数据体。据12口井资料统计,嫩一段上部泥岩的波阻抗值为5 000~8 000 g/cm3·m/s,而粉砂岩的波阻抗值为7 500~11 500 g/cm3·m/s,二者的波阻抗差异明显[图 2(a)]。对于波阻抗差异明显的砂岩和泥岩,对地震资料进行90°相位转换即可赋予地震数据岩性意义[10-11][图 2(b)—(c)]。②通过地层时间模型(Stratal time model)把90°相位转换地震数据体转化为各个地震标志层均拉平的地层切片数据体(Stratal slice volume)[12-13],就制作成相对地质年代域数据体。③在地震标志层(T07和T1)控制之下,等比例内插形成地层切片层位数据体。从下到上多次浏览不同时代的地层切片,即可以弄清嫩一段水进过程中的沉积体系发育模式及时空演化规律;当地层切片不能完整显示地质模式,可能存在局部穿时现象时,采用非线性地层切片技术进行修正[14-15]。④在获得比较清晰的地震地貌图像后,应用三维显示及光源照明技术进一步突出地貌细节变化,突出沉积体的边界及内部结构。

下载原图 图 2 地层岩性界面和地震反射界面的相关性分析 Fig. 2 Correlation analysis of stratigraphic lithologic interface and seismic reflection interface
2.2 水道结构参数测量与换算

水道的深度和宽度均从地震资料测量获得。水道宽度为水道两岸之间的距离,或者从一侧天然堤到另一侧天然堤之间的距离,其测量主要基于地震地貌图(地层切片)并参考地震剖面获得;深度是从水道顶部到底部的距离,或者两侧天然堤连线的中点到水道底部的距离。在地震剖面上读取双程旅行时差,并使用嫩一段的层速度换算得到水道的深度。具体算法为

$ H = \frac{{{T_0}{V_0}}}{2} $ (1)

式中:H为深度或厚度,m;T0双程旅行时差,ms;V0为嫩一段的层速度,m/s,取值3 200。

水道坡降的测量方法是,沿水道轴线切地震剖面,沿距离水道发育层位最近的上部地震标志层T07 (嫩二段底部油页岩)拉平,测量水道两端对应的地震同相轴中心点之间的垂向双程旅行时差,然后根据式(1)换算为厚度落差,再根据如下方法计算出坡降

$ D = \frac{{\Delta H}}{L} $ (2)

式中:D为坡降,‰;ΔH为水道落差,m;L为水道首尾两端之间的直线距离,km。

曲率由水道在工区内的实际流径与水道首尾之间的直线距离的比值得到。

尽管单期水道有时候在地震地貌成像(地层切片)上能够被分辨出来,但受地震资料纵向分辨率低的限制,在地震剖面上很难识别。因此,从地震剖面上测量得到的宽度、深度等水道结构参数实际上是水道带或水道复合体的结构参数。

3 水道的地震响应识别、结构参数确定及沉积物分析 3.1 水道的地震识别

地震地层切片成像结果显示(图 3),松辽盆地LHP地区嫩一段发育A—D和E—H等2组特征截然不同的水道系统。第一组包括A—D等4条水道,其宽度较大,呈蛇曲带状穿越工区。其中,水道A发育在工区西南部,NW—SE向呈弧形延伸约3.5 km后并入水道B。水道B呈南北向贯穿整个工区,延伸10 km。水道C和D发育在工区东侧,并在东部边界中段相交而过,显然二者并未同时活动。其中,水道D在工区内最宽,在曲流环发育部位宽度可达500 m左右。对地层切片进行多光源照射发现,水道D实际上为1个水道带,由4期以上的水道复合而成,在3个发育完好的曲流环之上可以清晰地观察到因水道迁移摆动而形成的单期水道[图 3(a)]。进一步观察发现,经黄光照射的彩色振幅地层切片[图 3(a)],水道A—D的地震异常均呈不均一分布,在水道延伸方向,多段不规则分布的纺锤形强振幅异常通过各自变细的端部相互连接,呈现宽-窄交互式分布。第二组自西而东分别是E,F,G,H,I,水道窄,呈较为顺直的线状,地震成像相对模糊[图 3(d)—(e)],推测大多为单一重力流事件形成的水道。其中,水道F最为清晰,可完全追踪,它呈S型南北向贯穿整个工区,在下游(工区西南角水道A和水道B交汇部位附近)频繁决口形成多个决口水道,展现出末端扇的特征[图 3(b)]。水道G,H,I的上游部分(工区东北端)地震成像较为模糊,但在流经水道C发育区域后清晰可见,均呈NNE—SSW向顺直状延伸至工区最南端断层附近消失。其中,水道H在W1井东侧偏北2 km处出现分流,分流水道H1和H2呈约20°的夹角向西南方向延展。过断层至W5井区,水道G及H1,H2的地震成像均较为模糊,从隐约可见的线条状强振幅异常分析,它们仍然向前延展,但开始决口分叉,接近末端。水道E只在工区西北部呈S型NW—SE向延伸约4 km后消失。水道I从北向南穿过水道C所在区域之后频繁决口改道,并最终消失在水道C以南3~4 km处,形成完整发育的末端扇。该末端扇在地震振幅切片上显示为强振幅,由5条以上决口水道及其漫溢沉积组成(图 3)。从沉积过程看,这样的末端扇又可称为末端分散体系(Terminal splay)[16]

下载原图 图 3 LHP地区嫩一段S0油层组地震地貌图像及重力流水道地质解释图 (a)地层切片(黄光照射);(b)重力流水道地质解释图;(c)连井沉积剖面;(d)未解释的地震地貌图像,位置见(a)图中白色方框;(e)为(d)解释后的地震地貌,T1—T4为水道发育的先后时间顺序,黑色实线(T4)为最后一期水道位置 Fig. 3 Seismic geomorphologic image and geological interpretation of gravity flow channels of Nen 1 member in LHP area

第一组水道A,B,C,D的断面形态大多呈对称或不对称的U型,但不同水道之间、同一条水道上下游的不同位置其断面形态不同。①水道A的底界呈宽阔且两侧对称的U型,深度为15~20 m(1个地震同相轴),宽度约为100 m。下切植入下伏地层,内部反射为强振幅。水道A的上部呈现丘状反射,在横断面上呈“透镜状”[图 4(a)—(b)],推测该水道富砂且存在水下天然堤。②水道B在北部上游地区的断面呈比较对称的V型[图 5(b)],下切侵蚀现象清晰可见,时间厚度接近1个同相轴;在南部的断面演变为右陡左缓不对称的U型[图 4(a)—(b)],宽度为50~100 m。较陡的右岸为水道的凹岸(外侧),以侵蚀作用为主;较缓的左岸为水道的凸岸(内侧),以沉积为主,地震地貌成像结果显示出同样的特征[图 3(a)]。当与水道A汇合后,宽度变为150~180 m,截面特征与水道A相似[图 4 (c)—(d)]。③水道C和D的横断面同样呈现U型,对下伏地层的下切侵蚀作用强烈[图 4(e)]。其中水道C的宽度为200~270 m,水道带D在工区内上游的宽度为230 m,到下游的大型曲流环一带宽度可达340~370 m。[图 4(f)—(g)]是横穿水道D的2个曲流环的地震剖面,在下切谷下部的红色强振幅反射中隐约可见向西倾斜的叠瓦状反射,说明水道由东向西侧向迁移摆动,西侧为侵蚀岸(凹岸),东侧是沉积岸(凸岸),这与图 5中的平面特征一致。水道D两岸地形明显高于周围地区,表现出分别向岸外倾斜的反射特征,西侧和东侧的倾角约为50°和30°,显示发育有3~5 m高的水下天然堤,且保存完整[图 4(f)]。

下载原图 图 4 水道A,B,C,D的地震响应特征(沿T07拉平)(剖面位置见[图 3(a)]) (a)S2剖面,水道A和B的分别呈现出标准的U型和不规则U型或宽V型,水道A的上拱地震反射可能暗示天然堤的存在以及水道可能富砂;(b)S2剖面,水道A的地质解释;(c)S1剖面,与水道A合并之后的水道B地震剖面响应;(d)S1剖面水道B的地质解释。横断面为规则的U型,上拱的地震反射预示天然堤的存在以及富砂的水道充填;(e)S7剖面,水道C的右岸和水道D的左岸直立;(f)S5剖面,示水道带的U型底界、陡峭的水道壁,以及水道两侧比外围明显高出的正振幅(红色)异常,预示存在着天然堤;(g)S3剖面,宽阔的U型下切谷,天然堤不明显。水道底界呈宽阔而对称的U型,水道充填物的中下部波峰(红色)强振幅反映富砂特征,上部波谷(蓝色)富泥 Fig. 4 Seismic response characteristics of channel A, B, C and D

第二组水道E—I在地震剖面上大多为强振幅反射,难以分辨对下伏地层的侵蚀[图 5(a)—(b)]。只有水道F在下游决口分叉前有1 km左右的地段对下伏地层产生比较明显的侵蚀,在地震剖面上表现为尖锐的V型[图 5(b)],时间厚度约1/3地震同相轴,揭示局部地形的突然变化。

下载原图 图 5 水道E,F,G,H,I的地震剖面响应(剖面位置见[图 3(a)]) (a)S6剖面,水道对下伏地层的侵蚀很小,地震剖面上为微弱的V型,或者仅表现为强振幅能量团(红色);(b)S4剖面,B,F呈现为V型,D为不对称U型,G,H1和H2表现为强振幅能量团(红色),对下伏地层侵蚀小 Fig. 5 Seismic section response of channel E, F, G, H, I
3.2 水道结构参数及其相关性

对各水道实际流径、水道首尾直线距离、坡降和曲率等参数进行统计和相关性分析表明,水道的曲率与其坡降具有很好的正相关关系,曲率大的水道(A,C)具有较大坡降,曲率小的水道坡降较小(G,I等)[图 6(a)]。采用随机采样方式,从地震剖面上测量各水道的宽度和深度,将这些统计数据进行散点拟合,发现水道的宽度和深度之间也具有较好的正相关关系[图 6(b)]。

下载原图 图 6 水道结构参数交会图 Fig. 6 Cross plot of channel structure parameters

比较水道平均深度与坡降、曲率关系发现,水道深度与曲率、坡降之间总体呈正相关关系[图 6(c)—(d)]。总体而言,研究区水道的深度、宽度、曲率、坡降之间具有较好的相关性。这与Wood等[17]对南美洲Trinidad和Tobago地区深海水道体系的研究结果相似。Posamentier等[18]在研究印度尼西亚Kalimantan地区Borneo东部海域更新统的水道-天然堤系统时,就发现深水水道的曲率与海底斜坡的坡降呈正相关关系,海底坡度增加,水道深度增大,曲率也增大。可见海洋和湖泊的深水水道的结构特征具有一定的相似性。

3.3 水道沉积物分析

LHP地区的10口钻井没有一口真正钻遇水道,而统计分析表明,粉砂岩与泥岩具有明显的波阻抗差异(参见图 3),因此,在地震剖面上水道内的强振幅异常推测为砂质沉积。从邻区钻遇深水水道的Ta-X井来看,水道内的沉积物以细砂岩为主,厚度为25 m,发育块状层理、平行层理、波状交错层理,常见暗色泥砾沿层理面叠瓦状排列,局部见厚度可达10 cm的泥砾层(图 7)。水道位置对应的测井曲线呈箱状,地震地貌特征为条带状强振幅异常,剖面特征为宽阔的U型,宽度可达600 m[2]。LHP地区水道A—D的地震特征与Ta-X井区的情况相似,且均位于整个水道系统的中游,因此,推测水道内的沉积物为细砂岩沉积。水道外侧主要发育漫溢沉积。录井资料揭示,这些漫溢沉积以粉砂岩、泥质粉砂岩为主,其厚度通常小于1 m,夹于大套暗色泥岩当中[参见图 1(b)]。邻区资料表明,漫溢沉积中常见厘米级波状层理细砂岩与粉砂质泥岩互层[图 7(e)]。

下载原图 图 7 Ta-X井嫩一段重力流水道沉积物岩心照片 (a)块状层理细砂岩,1 259.21~1 259.43 m;(b)示细砂岩中叠瓦状排列的泥砾,1 251.40~1 251.62 m;(c)平行层理细砂岩,1 249.50~ 1249.7 m;(d)泥砾层,水道滞留沉积,1 250.50~1 258.59 m;(e)波状层理细砂岩与粉砂质泥岩互层,1 227.82~1 227.93 m Fig. 7 Core photos of channel fills of Nen 1 member in well Ta-X

水道A,D在局部地段发育有水下天然堤,天然堤与外围斜坡区漫溢沉积的地震反射振幅相近,因而预测具有相似的岩性特征,为粉砂岩或泥质粉砂岩[参见图 4(a)—(d)]。需要注意的是,水道D局部地段的天然堤较外围斜坡区漫溢沉积的地震振幅略强,揭示强水流背景下的较粗粒溢岸沉积[参见图 4(f)—(g)]。

实际上,水道C,D以及B的上游地段在地震剖面上可以分为上下2个部分(参见图 4),中下部侵蚀谷内沉积物的地震振幅(红色)明显大于水道两侧的漫溢相(浅红色),应为细砂岩—粉砂岩,厚度为10~20 m,而上部沉积物的反射呈反极性特征(蓝色),从图 2(a)的岩石物理特征分析结果判断,其岩性应为泥岩,代表了重力流水道被废弃后,残余水道谷地在安静深水环境下的沉积。水道A,B在下游地区,其横截面在水道中轴线上地震呈上拱反射[参见图 4(a)—(b)],这可能与差异压实作用有关,富砂的水道沉积物在后期成岩作用中压实率小,而泥质的天然堤及岸后溢岸沉积压实率大[16, 19]。对比发现,水道A,B在下游地区的沉积厚度更大。

水道E—I在地震剖面上显示为比围岩相对较强的亮点地震反射,在振幅切片上显示较水道A— D更弱的振幅异常。根据振幅与岩性的相关性分析认为,水道E—I的沉积物的泥质含量略高,很可能大多为粉砂及泥质粉砂岩,且厚度薄,小于1/4子波长度。

4 水道系统时空演化与成因机制 4.1 水道系统时空演化

浏览对比上下层位的地震切片是判别沉积体发育先后顺序的有效手段[20]。在地震属性平面图上可以看出,水道E—I呈线型,总体曲率小,顺直分布,侵蚀作用弱,而水道A—D呈带状,曲率较大,因侧向迁移摆动常形成宽度较大的曲流环。因而,这2组水道系统发育时限靠近,利用地层切片序列无法把E—I与A—D完全剥离,但从E被B,F被B和A,I被C切割的情况分析,E—I的发育时间要早于A—D(图 8)。这2组水道系统可能形成于不同的气候旋回,A—D显然形成于强降雨和洪水事件中,保持着强有力的下切侵蚀作用,形成深度达10~ 20 m的斜坡峡谷,而E—I已经处于其发育末端的LHP地区。

下载原图 图 8 地震振幅的地层切片序列 (a)—(d)从早到晚的地层切片,示水道C发育时期略早于A,B,D,而水道E,F,G,H,I是最早发育的水道系统;(e)地层切片在剖面上的位置 Fig. 8 Strata slice series of seismic amplitude

Posamentier等[16]将深水浊流沉积体系按照源—汇系统的思路划分为上、中、下游3段。上游靠近陆上,发育下切侵蚀的支流峡谷;中游发育以限制性水道为特征的单一水道体系和天然堤;下游发育非限制性水道,决口改道频发,是沉积物的主要卸载区,多条先后决口的水道呈分散状分布,构成末端扇。PHL地区的水道A—D单一水道复合体特征明显,局部保留天然堤,可类比于Posamentier等[16]划分的浊流沉积体系的中游,而水道E—I处于重力流系统的末端,决口分流频发,可类比于其浊流沉积体系下游末端扇部分。

就A—D等4条水道而言,在重力流水道集中发育层段自下而上制作的地层切片中,水道C的发育时间略早于A,B,D。最下部的地貌特征显示,在A,B,D出现之前,水道C表现为条带状强振幅异常,已经在工区东北部开始发育[图 8(a)],而图 8 (b)—(c)反映的是过渡层段的地震影像,水道A— D同时存在是由于地震成像过程中子波旁伴效应的影响,不能同时反映它们的真实发育位置,图 8(d)为该层段顶部的地震地貌,水道A,B,D的影像清晰,而水道C已经不见踪影。A,B,D等3条水道很可能是在同一期、或者时间上非常靠近的洪水事件中形成的。

4.2 水道成因机制 4.2.1 水道形成与演化的控制因素

松辽盆地嫩一段沉积于湖盆快速扩张期,其下部(S1)以砂泥岩互层为主,为滨浅湖沉积;上部(S0) 以厚层黑色、灰黑色泥岩夹薄层粉砂岩沉积为特色,为半深湖—深湖沉积,水道沉积主要发育在这套厚层泥岩当中。上覆嫩二段底部的油页岩,厚度约为4 m[参见图 1(b)],是湖盆水体扩张达到鼎盛时期的标志。

冯志强等[5]认为,嫩一段沉积时期大庆长垣地区从三角洲前缘演变到十分发育的深湖区重力流水道,延伸直线距离达50~100 km,在水道末端常见湖底扇。潘树新等[2]发现在齐家—古龙地区嫩一段上部(S0)发育多支大型异重流成因的水道-湖底扇系统,直线延伸距离超过80 km,宽度为100~ 900 m。

无论是限制性水道系统A—D,还是非限制性水道系统E—I,均属于异重流成因,是洪水期陆上河流携带大量泥沙直接入湖,并经斜坡区长距离奔流而成。水道A—D发育期间,包括松辽盆地在内的广大地区经历了强降雨天气,河流水位暴涨,泥沙含量剧增,携带大量泥沙的河流向盆地汹涌而来,由于河水密度远大于湖泊静水水体的密度,因而河流水体入湖后在惯性力及重力作用下,沿着湖水底部向盆地中央深水区强力推进,并对湖盆斜坡区形成强烈侵蚀,形成多个峡谷系统。后期,降雨量减小,洪水能量减弱,导致水下异重流的流速降低,以底床和悬浮负载方式搬运的细砂—粉砂质沉积物得以沉积充填于水道当中。最后,水道被废弃,在残余的水道谷地内沉积了深水泥质沉积。

黄清华等[21]认为,松辽盆地青一段和嫩一段、嫩二段发育的3套厚层暗色泥岩、油页岩、页岩分别对应于3期全球性的缺氧事件,导致了湖盆内介形虫、叶肢介等生物的3次变革。缺氧事件通常会导致大气密度降低,光照强度及蒸发量增大,进而引起温度升高和降雨量增大[22]。因此,缺氧事件可能是松辽盆地嫩一段沉积时期气候异常、洪水泛滥、盆地内部深水区广泛发育异重流水道的主要原因。

前文已经指出,与水道A—D对比,E—I发育时间略早,处在盆地整个重力流水道系统的末端,而A—D则处于稍晚发育的水道系统的中游地带。这2组水道出现显著差异的原因,可能既有构造变动的影响,也有气候变化的因素。冯志强等[1]研究表明,嫩江组沉积时期松辽盆地处于构造缓慢沉降阶段。从图 6(a)看出,水道E—I的坡降较A—D的更小,前者为0.55‰~0.69‰,后者为0.61‰~ 1.0‰。这些数据印证了嫩江组沉积时期盆地持续沉降的认识。相对较陡的湖泊底形可能是异重流水道A—D较E— I在湖泊底部能够长距离流动的原因之一。另外,气候变化也可能是导致2组重力流水道在湖泊底部延伸距离存在差异的原因。水道E—I发育时期只是气候条件剧变的开始,此时松辽盆地所在流域开始发生大面积的洪水事件,高密度河水沿着水下沟谷以异重流的方式向湖盆深水区推进。到水道A—D发育时期,气候变化达到极端条件,较前期有更大的洪水事件,也就形成了在湖盆底部延伸更远的异重流沉积。同时,强降雨导致大量的浑浊地表水流沿着河道注入湖泊,引起湖平面快速上升,也为嫩一段上部(S0)发育厚层泥岩创造了条件。

嫩一段水道系统发育之后,盆地在相当长的时期内都处在相对宁静时期,湖盆水体安静,水生藻类繁盛,发育了嫩二段底部大面积分布的油页岩,其TOC的质量分数为5.3%~11.9%(平均为7.9%),是松辽盆地一套极其优质的烃源岩[23]。也正是由于以这套油页岩为代表的安静的深水沉积环境,才使得嫩一段的重力流水道体系得以完整地保存。

4.2.2 水道沿流向“宽-窄”相间的成因

水道A,B,C在地震地貌图上沿水流方向均呈“宽-窄”间互分布的地震振幅异常[参见图 3(a)]。利用非线性地层切片技术上下浏览观察,发现这种现象反映了地下真实的沉积地貌特征。与水道D宽-窄变化原因一样,水道A,B,C的宽-窄相间变化为水道侧向迁移摆动所致,只是地震分辨率有限而无法呈现其细节。图像中出现的水道A,B,C,D实际上不是由单期水道形成,而是由一定时期内发育的多期水道侧向迁移摆动形成的水道带,其中的较宽部分在位置上可与陆上河流的点坝对比,为曲流水道沿凹岸反复侵蚀、在凸岸不断沉积所致,是曲流环(Meandering Loop)的位置之所在;较窄的部位为2个相邻、但摆动方向相反的曲流环之间的转换部位,其特征是新水道对老水道的交叉切割[参见图 3(d)—(e)]。

Posamentier等[16]将深水水道划分为有序水道序列(Organized channel successions)和无序水道序列(Disorganized channel successions)等2种类型。有序水道序列是指重力流事件结束后,水道未被完全充填,从而为后续水流留下空间,这样曲流水道的凹岸不断被后续重力流侵蚀,在凹岸靠下游位置遭受的侵蚀作用尤其强烈,导致水道在侧向摆动(Swing)的同时,还常常引起曲流环向下游方向不断迁移,形成扫动(Sweep)特征。通过曲流环扫动方向分析还可以判断古水流方向。无序水道序列是指重力流事件结束时,水道已经被沉积物完全充填,因而后续水流只能另辟蹊径,从而形成无序分布的水道模式。从形态上看,水道A,B,C,D与有序水道序列的特征更为接近,水道D体现得更为明显[图 3(d)—(e)]。

4.3 讨论

(1) 深水重力流水道与天然堤相伴生似乎是海相盆地的常见现象[20, 24-27],特别是在重力流水道的中游,受陆上供源河流洪水事件的影响,流经水道的水流时常高于水道壁,从而产生漫溢沉积,形成水下天然堤[17],而对于陆相盆地,鲜有深水区发育天然堤的报道。从部分地震剖面看(参见图 4),松辽盆地嫩一段发育的限制性水道系统(水道A,B,D)在局部地段水道两侧的沉积地貌明显高于外围同期沉积,预测存在着天然堤,但是从目前的地震地貌图上还无法确定其分布范围。这或许与其规模小有关,也可能是大部分水道天然堤高度普遍较小,经后期成岩压实作用后幅度更低,无法用地震资料分辨出来。

(2) 受沉积期后差异压实等后生成岩作用的影响及地震资料分辨率的限制,在地震反射剖面上所测量的坡降、水道宽度、深度,以及天然堤高度等数据与真实的沉积地貌之间存在一定的误差。此外,文中水道相关结构参数实际上是水道带的结构参数。尽管如此,这些数值至少可以帮助人们半定量认识湖泊中深水水道的结构特征以及储集层发育规模。

5 结论

(1) 松辽盆地LHP地区嫩一段发育2期重力流水道沉积事件,早期的水道体系E—I为非限制性水道体系,处于重力流系统的下游靠末端,末端多见决口水道,对下伏地层侵蚀弱,水道窄,水下天然堤不发育或发育程度很弱。晚期水道体系A—D为限制性水道体系,处于大型水道系统的中游,能量强,对下伏地层下切侵蚀达10~20 m,以单一水道复合体形式存在,局部发育天然堤。晚期的重力流事件还可以进一步细分为早、晚2个亚期,水道C发育时间略早于水道A,B和D。

(2) 水道沉积物主要为细砂岩和粉砂岩,厚度可达20 m,宽度50~370 m,局部宽度可达500 m以上。这些水道沉积与构造背景匹配可以形成构造-岩性圈闭,是夹持在优质烃源岩当中的高效勘探目标。

(3) 水道A—D与E—I等2套水道系统存在显著的差异,一方面是松辽盆地发生持续沉降,湖盆坡度发生了比较明显的变化。另一方面是古气候变化导致的降雨量显著增大,导致了异重流水道A—D较E—I具有更强的侵蚀能力和水下延伸距离。导致降雨量剧增的直接原因很可能是当时的全球性的缺氧事件。

致谢: 成文过程中得到Posamentier、顾家裕教授的悉心指导,黄银涛、于兴河教授提出了许多建设性的建议,在此深表谢意。

参考文献
[1]
冯志强, 张顺, 解习农, 等. 松辽盆地嫩江组大型陆相坳陷湖盆湖底水道的发现及其石油地质意义. 地质学报, 2006, 80(6): 1226-1233.
FENG Z Q, ZHANG S, XIE X N, et al. Discovery of a large-scale lacustrine subaqueous channel in the Nenjiang Formation of the Songliao Basin and its implication on petroleum geology. Acta Geologica Sinica, 2006, 80(6): 1226-1233.
[2]
潘树新, 刘化清, ZAVALA C, 等. 大型坳陷湖盆异重流成因的水道-湖底扇系统: 以松辽盆地白垩系嫩江组一段为例. 石油勘探与开发, 2017, 44(6): 911-922.
PAN S X, LIU H Q, ZAVALA C, et al. Sublacustrine hyperpycnal channel-fan system in a large depression basin: A case study of Nen 1 member, Cretaceous Nenjiang Formation in the Songliao Basin, NE China. Petroleum Exploration and Development, 2017, 44(6): 911-922.
[3]
莫午零, 吴朝东, 张顺. 松辽盆地北部嫩江组沉积物重力流形成机制. 天然气地球科学, 2013, 24(3): 555-565.
MO W L, WU C D, ZHANG S. Generation mechanism of sediment gravity flow of Nenjiang Formation in the north Songliao Basin. Natural Gas Geoscience, 2013, 24(3): 555-565.
[4]
FENG Z Q, ZHANG S, CROSSW T A, et al. Lacustrine turbidite channels and fans in the Mesozoic Songliao Basin, China. Basin Research, 2010, 22(1): 96-107. DOI:10.1111/j.1365-2117.2009.00442.x
[5]
冯志强, 张顺, 付秀丽. 松辽盆地姚家组-嫩江组沉积演化与成藏响应. 地学前缘, 2012, 19(1): 78-88.
FENG Z Q, ZHANG S, FU X L. Depositional evolution and accumulation response of Yaojia-Nenjiang Formation in Songliao Basin. Earth Science Frontiers, 2012, 19(1): 78-88.
[6]
萧德铭, 迟元林, 蒙启安, 等. 松辽盆地北部向斜区岩性油藏勘探认识与实践. 北京: 石油工业出版社, 2005: 1-65.
XIAO D M, CHI Y L, MENG Q A, et al. Understanding and practice of lithologic reservoir exploration in syncline area of northern Songliao Basin. Beijing: Petroleum Industry Press, 2005: 1-65.
[7]
ZENG H L, LOUCKS R G, BROWN L F. Mapping sedimentdispersal patterns and associated systems tracts in fourth-and fifth-order sequences using seismic sedimentology: Example from Corpus Christi Bay, Texas. AAPG Bulletin, 2007, 91(7): 981-1003. DOI:10.1306/02060706048
[8]
ZENG H L, TUCKER F H. High-frequency sequence stratigraphy from seismic sedimentology: Applied to Miocene, Vermilion Block 50, Tiger Shoal area, offshore Louisiana. AAPG Bulletin, 2004, 88(2): 153-174. DOI:10.1306/10060303018
[9]
POSAMENTIER H W, DAVIES R J, CARTWRIGHT R A, et al. Seismic geomorphology: An overview//DAVIES R J, POSAMENTIER H W, WOOD L J, et al. Seismic geomorphology: Applications to hydrocarbon exploration and production. Special Publications of The Geological Society of London, 2007, 277: 1-14.
[10]
ZENG H L, BACKUS M M. Interpretive advantages of 90 degreesphase wavelets: Part 1-Modeling. Geophysics, 2005, 70(3): C7-C15. DOI:10.1190/1.1925740
[11]
ZENG H L, BACKUS M M. Interpretive advantages of 90 degrees-phase wavelets: Part 2-Seismic applications. Geophysics, 2005, 70(3): C17-C24. DOI:10.1190/1.1925741
[12]
ZENG H L, BACKUS M M, BARROW K T, et al. Stratal slicing, Part Ⅰ: Realistic 3-D seismic model. Geophysics, 1998, 63(2): 502-513. DOI:10.1190/1.1444351
[13]
ZENG H L, HENRY S C, RIOLA J P. Stratal slicing, Part Ⅱ: Real seismic data. Geophysics, 1998, 63(2): 514-522. DOI:10.1190/1.1444352
[14]
ZENG H L. Seismic imaging for seismic geomorphology beyond the seabed: Potentials and challenges//DAVIES R J, POSAMENTIER H W, WOOD H W, et al. Seismic geomorphology: Applications to hydrocarbon exploration and production. Special Publications of The Geological Society of London, 2007, 277: 15-28.
[15]
刘化清, 苏明军, 倪长宽, 等. 薄砂体预测的地震沉积学研究方法. 岩性油气藏, 2018, 30(2): 1-11.
LIU H Q, SU M J, NI C K, et al. Thin bed prediction from interbeded background: Revised seismic sedimentological method. Lithologic Reservoirs, 2018, 30(2): 1-11.
[16]
POSAMENTIER H W, KOLLA V, 刘化清. 深水浊流沉积体系综述. 沉积学报, 2019, 37(5): 880-903.
POSAMENTIER H W, KOLLA V, LIU H Q. An overview of deep-water turbidite deposition. Acta Sedimentologica Sinica, 2019, 37(5): 880-903.
[17]
WOOD L J, MIZE-SPANSKY K L. Quantitative seismic geomorphology of a Quaternary leveed-channel system, offshore eastern Trinidad and Tobago, northeastern South America. AAPG Bulletin, 2015, 93(1): 101-125.
[18]
POSAMENTIER H W, KOLLA V. Seismic geomorphology and stratigraphy of depositional elements in deep-water settings. Journal of Sedimentary Research, 2003, 3: 367-388.
[19]
POSAMENTIER H W. Seismic geomorphology: Imaging elements of depositional systems from shelf to deep basin using 3D seismic data: Implications for exploration and development. Geological Society London Memoirs, 2004, 29(1): 11-24. DOI:10.1144/GSL.MEM.2004.029.01.02
[20]
ZENG H L. Predicting geometry and stacking pattern of thin beds by interpreting geomorphology and waveforms using sequential stratal-slices in the wheeler domain. Interpretation, 2015, 3(3): SS49-SS64. DOI:10.1190/INT-2014-0181.1
[21]
黄清华, 陈春瑞, 王平在, 等. 松辽盆地晚白垩世生物演化与古湖泊缺氧事件. 微体古生物学报, 1998, 15(4): 417-425.
HUANG H Q, CHEN C R, WANG P Z, et al. The late Cretaceous bio-evolution and anoxic events in the ancient lake in the Songliao Basin. Acta Micropalaeontologica Sinica, 1998, 15(4): 417-425.
[22]
POULSEN C J, TABOR C, WHITE J D. Long-term climate forcing by atmospheric oxygen concentrations. Science, 2015, 348(6240): 1238-1241. DOI:10.1126/science.1260670
[23]
贾建亮, 刘昭君, ACHIM B, 等. 松辽盆地嫩江组油页岩发育控制因素. 地球科学——中国地质大学学报, 2014, 39(2): 174-186.
JIA J L, LIU Z J, ACHIM B, et al. Major factors controlling formation of oil shale in Nenjiang Formation of Songliao Basin. Earth Science-Journal of China University of Geosciences, 2014, 39(2): 174-186.
[24]
POSAMENTIER H W, WALKER R G. Deep-water turbidites and submarine fans. SEPM Special Publication, 2006, 84: 397-520.
[25]
CALLOW R H T, KNELLER B, DYKSTRA M, et al. Physical, biological, geochemical and sedimentological controls on the ichnology of submarine canyon and slope channel systems. Marine and Petroleum Geology, 2014, 54: 144-166. DOI:10.1016/j.marpetgeo.2014.02.016
[26]
庞雄, 朱明, 柳保军, 等. 南海北部珠江口盆地白云凹陷深水区重力流沉积机理. 石油学报, 2014, 35(4): 646-653.
PANG X, ZHU M, LIU B J, et al. The mechanism of gravity flow deposition in Baiyun Sag deep water area of the northern South China Sea. Acta Petrolei Sinica, 2014, 35(4): 646-653.
[27]
卜范青, 张旭, 陈国宁. 尼日尔三角洲盆地重力流沉积模式及储层特征: 以AKPO油田为例. 西安石油大学学报(自然科学版), 2017, 32(1): 64-70.
BU F Q, ZHANG X, CHEN G N. Gravity flow depositional mode and reservoir characteristics of Niger Delta Basin: Taking AK-PO Oilfield as an example. Journal of Xi'an Shiyou University(Natural Science Edition), 2017, 32(1): 64-70. DOI:10.3969/j.issn.1673-064X.2017.01.010