1998年，曾洪流等首次提出了“地震沉积学”概念^[1]，标志着这门学科的诞生。2004年，曾洪流教授提出地震沉积学是“利用地震资料来研究沉积岩及其形成过程的一门学科；地震岩石学与地震地貌学组成了地震沉积学的核心内容”^[2]。之后，曾洪流、Schlgaer、林承焰、董春梅、朱筱敏等国内外学者对地震沉积学的概念、方法及技术进行了研究，地震沉积学理论体系逐渐丰富和完善^{[3, 4, 5, 6, 7]}。近年来，国内学者针对中国陆相复杂储集层开展地震沉积学研究，在油气勘探与开发中取得了显著的效果，但多集中在三角洲和曲流河沉积体系地层中^{[8, 9, 10]}。湖相重力流体系的地震沉积学研究无论国外还是国内均较少涉及。本文基于歧口凹陷采集、连片处理的高精度三维地震资料，以黄骅坳陷歧口凹陷歧南斜坡沙一段(Es₁)的深水厚层砂岩为研究对象，应用地震沉积学的理论和方法，分析了断陷湖盆以沟道型砂砾质碎屑流为主体的重力流沉积的地震相模式，探讨了地震沉积学研究方法在重力流砂体形态刻画和储层预测中的应用，形成了一套适合本地区重力流水道储层预测的思路和方法，试图丰富和发展地震沉积学研究成果。

黄骅坳陷歧南斜坡是目前大港油田岩性地层圈闭勘探的重点地区之一，主要勘探目的层为沙河街组沙一段(Es₁)。歧南斜坡位于渤海湾盆地歧口凹陷西南部，其南部围限于埕宁隆起和羊三木凸起之间，北部夹持于沙三段和沙一段均强烈断陷的南大港断层和张北断层之间^[11]，总体上为向北东方向倾伏的单斜，面积大约900 km²。在歧口凹陷沙一段和沙三段底部构造图上，研究区均位于东西两侧正向构造之间的低洼部位，表明研究区为继承性发育的沉积凹槽，油田上俗称“歧南水道”，同时，区内赵北、张北、南大港等控盆断裂强烈活动，断层落差均在300 m以上，湖盆深陷扩张，沉积坡度陡，具备形成重力流的动力学背景(图 1)。

图 1 研究区构造平面位置图 Fig. 1 Sketch map of structural location in the study area

图选项

歧南斜坡为三维地震覆盖区，地震资料品质较好，沉积现象丰富，构造相对较简单，地层岩性以砂泥岩为主，比较适合开展地震沉积学的研究工作。因此，笔者在建立高频等时地层格架基础上，应用地震沉积学研究方法和技术，对研究区沙一段进行地震地貌恢复、属性提取、90°相位转换、地层切片及三维可视化等综合研究，并结合岩芯、钻测井资料，精细刻画了重力流水道沉积体的平面展布形态和规模，明确了砂体的分布规律，为研究区油气勘探奠定了坚实基础，并取得了良好的地质效果。

2.1 结合分频技术和时频分析，建立高频等时格架

建立等时地层格架既是层序地层学的精髓，也是地震沉积学中一系列关键技术应用的前提^[12]。

区域层序地层学研究表明，黄骅坳陷沙一段为一个完整的二级层序，内部可进一步划分为3个三级层序，分别对应沙一上亚段(Es^s₁)、沙一中亚段(Es^z₁)和沙一下亚段(Es^x₁)^{[13, 14]}。地震沉积学认为原始地震剖面上同相轴不完全是等时的，其等时性受地震频率控制，对地震剖面进行分频后，产状和位置在不同的分频剖面上表现一致的同相轴其等时性较可靠^[12]。本文运用小波变换计算不同频率地震体来识别层序界面，通过对地震反射同相轴等时性分析，揭示沙一段(Es₁)顶、底界面，沙一中亚段(Es^z₁)顶部、底部稳定泥岩具有很好的等时性，可作为地震沉积学研究的参考界面。在识别三级层序界面的基础上，由于不同级别的地质层序体对应着地震剖面上的不同频率特征，利用时频分析技术对关键井按不同频率进行扫描分析可以识别出由大到小的各级层序体，纵向上频率变化的方向性代表了岩性粗细的变化^[15]，如图 2a中显示了沙一中亚段砂泥岩频谱变化的关系、层序转换点和频率变化点的相对应，将沙一中亚段进一步划分为4期中期旋回砂组，在关键井高频层序划分的基础上，井震结合，对骨架剖面进行基准面旋回划分与对比，并在30 Hz高频地震数据体上识别追踪四级层序界面。通过时频分析的地震旋回划分方法所识别和标定的层序界面与地震资料中的高频成分和测井高频层序划分及旋回分析都吻合的较好，提高了旋回界面识别的准确性，从而建立了研究区高频等时地层格架(图 2b)。

图 2 歧南斜坡沙一中亚段高频层序地层格架(平面位置见图 1) a.qn2井时频分析；b.连井等时地层格架 Fig. 2 High frequency sequence stratigraphic framework of Esz1 in Qi′nan slope

图选项

最小等时单元是地震相、测井相和岩芯相三者结合时，井震统一对比的最小研究尺度^[12]。歧南地区沙一段地震反射波主频为20 Hz，频带范围7~35 Hz，砂岩层速度4 400 m/s左右，可以识别的砂岩厚度在25 m左右，沙一中亚段砂组厚度大约在40~70 m左右，砂组正韵律下部、反韵律上部砂组厚大于20 m可分辨，从图 2a可以看出，测井上可以划分到准层序或更小的层序单位，但由于地震分辨率较低，准层序界面在地震剖面上无法追踪，为了达到井震研究尺度的统一，将准层序组作为最小等时研究单元^[12]。

重力流水道砂体是指由重力流或浊流在湖盆内的断凹或沟槽中所形成的带状碎屑砂体，它可以堆积在浅水和深水中^{[16, 17]}。研究区内重力流水道砂体分布在赵北断层及其以北，主要见于沙一段，是来自南部埕宁隆起物源的辫状河三角洲前缘砂体沿歧南水道滑动、滑塌，在坡折带的下倾方向深水区形成的条带状粗碎屑岩体。

通过对5口取芯井进行系统岩芯观察，结合实验数据和测井相分析，认为沙一段属于典型的重力流水道沉积，具有以下沉积特征：

(1) 沙一中的泥岩以灰色、深灰色为主，平面分布稳定，具有质纯，有机质含量低的特点。而且在构造图上(图 1)，研究区两侧地势较高地区均为湖相泥岩，表明地势较低的歧南水道处于深水沉积环境。

(2) 岩芯观察沙一段岩石类型以细砂岩和含砾砂岩为主，进一步研究发现，沙一段厚层块状砂岩是由3种不同成因类型的砂岩相互叠置构成：一种为砂砾质碎屑流，以qn6井3 399~3 411 m井段为代表，粒度粗，砾石无分选、无粒序漂浮于较细的砂质基质中，砾石成分以碳酸盐为主，偶见直径2~5 cm的紫红色泥砾，沉积构造以块状构造为主(3 410.12 m)，显示在重力作用下以基质和杂基(水、泥及砂质颗粒)支撑的块体搬运机制以及以“冻结式”快速卸载堆积的沉积特点^[11]，在伽马测井曲线上表现为微齿化块状箱型，与下伏岩性呈突变接触，单砂层厚度较大，一般大于1m，主要为4~7 m，在重力流水道水道中心微相中最为常见。第二种为滑塌岩，以qn6井3 293~3 310 m取芯井段为代表，岩性主要为细砂岩、粉砂岩和粉砂质泥岩，表现为砂泥混杂并发生强烈的同沉积变形，常见包卷层理(3 294.6 m)、泥岩撕裂屑(3 295.00 m)、滑塌变形层理(3 307.86 m)及重荷构造(3 410.12 m)等，其电测曲线齿化明显，单砂层厚度主要为1~3 m，表现为多期钟形、漏斗形或指形的叠加。第三种为浊积岩，典型识别标志为递变层理，本区不太发育，常以砂、泥岩薄互层形式出现，在qn6井3 404.27 m井段岩芯可见平行层理，单砂层厚度为2~4 cm，由于层薄，测井曲线反映微弱(图 3a,b)。

(3) 根据岩性和沉积构造特征，研究区的重力流水道沉积常见6种岩相类型，宏观上重力流水道砂体的沉积剖面主要为大套暗色泥岩中发育厚层块状砂岩相和具变形层理砂泥岩相，总的剖面结构显示出中间粗、上下细的逆—正粒序结构，下部为反旋回、上部为正旋回(图 3b)。

(4) 沉积物粒度概率曲线以圆弧形为主，悬浮组分比例高，占30%~40%。据3 200.79~3 403.9 m井段的粒度统计结果显示，C值为120~1 200 μm，M值为35~330 μm，样点在C-M图上大致平行于C=M基线集中分布^[11]。

歧南斜坡沙一段沉积时期，受张北、南大港断裂构造影响，歧南斜坡发育地势相对较低的断槽或沟槽，发育于盆地缓坡边缘的辫状河三角洲沉积体系顺着缓坡向湖盆中心大幅推进；到达深湖后，在地震、火山喷发或水动力机制的诱导下，早期沉积的辫状河三角洲前缘沉积物就会顺着这些断槽进行搬运并发生沉积，从而在较深水环境中形成一系列平行盆地长轴方向的、顺断层走向分布的条带状重力流水道沉积，在水道末端地形较平坦的湖盆中心部位还会发育多期叠置的朵体(图 3c)。在湖盆内部因大断裂而产生的断凹或沟槽对重力流水道砂体的形成和分布起着重要的控制作用。

图 3 歧南斜坡重力流水道砂体沉积特征及沉积模式 a.重力流水道沉积构造(qn6井)；b.重力流水道沉积序列；c.歧南斜坡重力流水道沉积模式示意图 Fig. 3 Sedimentary characteristics and depositional sequence of gravity-flow channel in Qi′nan slope

图选项

通过对歧南斜坡沙一段重力流水道沉积特征分析及沉积模式的建立，为开展地震沉积学研究奠定了基础。

地震沉积学研究已经形成了较为规范的流程^[18]，目前应用最多的是90°相位转换、地层切片和分频技术等。本文在建立高频等时格架的基础上，应用地震沉积学的思路和方法对重力流水道砂体进行精细刻画。

在地震反射等时性分析约束下完成了沙一中亚段地震地貌学研究。古地貌恢复对沉积体系配置及储集砂体的展布的研究具有重要意义。常用的方法主要有残余厚度法、印模法、地球物理法、高分辨率层序地层法、沉积学方法等^{[19, 20, 21]}。目前，古地貌恢复的方法主要从构造恢复和地层恢复两方面入手^[21]，笔者主要针对研究区地层的恢复，由于古生物是确定古水深的可靠标志，研究区缺乏古生物资料，没有进行古水深校正，本文主要采用去压实校正和正断层断面厚度补偿相结合恢复歧南斜坡沙一中沉积前的相对古地貌，而不是绝对古地貌。由于歧南斜坡沙一段沉积时期构造相对稳定，是一较连续沉积过程，没有长时期的剥蚀，剥蚀量相对较小，研究中仅对斜坡边缘高部位剥蚀厚度明显的井点用地层横剖面对比法进行估算。主要基于压实前后地层骨架体积不变原理，首先分岩性建立压实曲线方程，采用回剥法进行单井地层厚度压实校正，并在钻井压实校正的基础之上，利用地震资料进行井间古地貌趋势补偿^[22]。同时，对正断层断面厚度进行补偿，避免因断面影响地震地貌出现假的古隆起，主要依据断层上、下盘最远端断点位置，删除断层区域，再应用克里格插值法恢复断层位置地层厚度。

古地貌恢复结果表明，沙一中亚段沉积时期，沿ch59—qn2—qn6—zh65—zh56井一线发育一个大型沟槽，沟槽内地层厚度较大，沙一中地层厚200~400 m，沟槽以外地层厚度减薄^[16]。大致在zh65—zh47井一带，相当于赵北断层的位置，发育古地形突然变陡的坡折带，在古地貌图上，该带沉积厚度等值线较密集(图 4a)。坡折带控制沉积相带分异^[23]，钻井取芯及测井资料证实，坡折带以南为辫状河三角洲沉积，坡折带之下于qn8—qn6—qn2井一带大规模发育重力流砂体。古地貌恢复可在宏观上预测重力流水道砂体的有利发育区，并且随着古地貌在沉积期上的变化，沉积物的卸载区也发生变化。沙一中沉积早期，研究区赵北断层下降盘qn8井以东的zh48—qn2井区古地貌上呈现明显的低洼，沉积物在此聚集(图 4b)，而到沙三中沉积晚期，南大港等控盆断层的强烈活动，沉积中心向西迁移，此时低洼区域主要位于qn8—qn2井区(图 4c)。

图 4 歧南斜坡沙一中亚段古地貌图和均方根振幅属性图 a.沙一中沉积前古地貌图；b.沙一中沉积早期古地貌图；c.沙一中沉积晚期古地貌图；d.沙一中均方根振幅属性图 Fig. 4 Paleogeomorphic map and RMS amplitude attribute of Esz1 in Qi′nan slope

图选项

研究区沙一中地层岩性以砂、泥岩为主，通过对地震属性提取与分析，优选出振幅类、统计类、方位角等多种属性与地层岩性、沉积体系的平面展布细节匹配较好，各属性之间也具有较好的相关性，与钻井资料吻合较好。这主要是由于沙一中以湖相为主体的沉积环境中砂岩波阻抗大于泥岩，砂体具有较高的波阻抗值，在水道两侧，均为低阻抗的湖相泥岩，在沙一中亚段的均方根振幅图上，重力流水道发育区为明显的强振幅的反射特征，大体上存在两个振幅高值条带。主水道一支位于qn8—qn3与qn6—qn2井之间，一支位于zh47—zh48井区(图 4d)。

地震地貌与宏观属性相结合显示歧南斜坡沙一中亚段以南部埕宁隆起物源体系为主，古沟槽和坡折带控制沉积体系发育。沙一中砂体分布局限，水道砂体位于赵北断层下降盘或三角洲前方沟槽内，由南向北呈条带状延伸。

对关键井水道砂体进行精细标定，从过qn8井垂直重力流水道的90°相位地震剖面可以看出，qn8井位于水道砂体侧翼，通过井震标定，清晰刻画了水道砂体的地震响应特征：横剖面上多期河道侧向迁移叠置，具有典型水道冲刷－充填反射特征，整体表现为顶平下凸的透镜体，其下伏层的湖盆泥岩遭受侵蚀冲刷，复合砂体内部结构清晰，单期水道表现为短轴状、不连续、强振幅的特征，水道间往往被较厚的泥岩隔层分开。根据地震反射外形的不同可进一步细分为“透镜状”、“梭状”和“蠕虫状”。其中，“透镜状”地震相呈底凹顶凸型，为下切水道的典型地震响应特征，通常对应于水道中心微相。“梭状”地震相呈“顶凹底凹”或“顶平底凹”的“梭状”外形特征，反映单期或多期突发性高能碎屑流对下伏地层的侵蚀冲刷。“蠕虫状”地震相中强振幅、断续、首尾相连呈叠瓦状排列，反映不同时期水道由于垂向叠加，侧向迁移形成的侧向加积式反射特征(图 5a)。

图 5 90°相位地震剖面与典型地层切片(沙一段) Fig. 5 90° phase and interpretation of stratal slices in Qi′nan slope

图选项

依据地震沉积学的基本原理，在垂向上无法识别的地质体，在平面上有可能通过横向分辨率被识别出来^[7]。通过90°相位地震资料基础上的地层切片技术可以较准确地在平面上分析沉积体的沉积特征，能更好地研究厚度比地震分辨率薄的沉积体。

利用地层切片技术对重力流水道砂体进行精细刻画，最关键的一步是选出具有地质时间界面意义的参照同相轴^[14]。因此，在该研究区内，在对地震数据滤波基础上，应用中石油西北分院研制的Geosed1.0软件的地震反射同相轴等时性分析技术，优选出沙一段的顶、底界面及沙一中亚段顶、底界面为等时层序界面，在对关键层序界面进行追踪解释的基础上，以90°相位化地震数据体为基础，在沙一段内部等比例内插，制作100张地层切片，目的是反映自沙一层序底部至顶部地震相的纵向变化特征，结合钻井岩芯校正，开展井震对比检验，明确重力流水道沉积微相平面分布与空间演化规律。

从层位上看，研究区重力流水道砂体主要位于沙一上底部与沙一中亚段^[16]。现对挑选的3张典型地层切片进行沉积相解释。由于区内沙一段90°相位化地震数据体的振幅信息与砂、泥岩具有较好对应关系^[11]。经过钻测井、岩芯资料证实了地层切片中红色与较厚层砂岩对应，黄色代表薄层砂岩或泥质砂岩，而蓝色对应的是泥岩层。在No58地层切片上(图 5b)，对应于沙一中下部沉积，此时处于初始湖泛期，来自南部埕宁隆起的物源供给充足，半深湖区重力流水道贯穿工区南北，规模大，对下伏地层侵蚀作用显著，是突发事件所致，重力流水道明显分为两支，自南向北发育，主水道位于qn8以东与qn6—zh48井区，到qn6井以北两支水道开始合并。到No38地层切片上(图 5c)，对应于沙一中上部地层切片，已达到最大湖泛期，主水道位于qn8和qn2井区，规模减小，两条水道近平行呈南北向顺直延伸，前端终于qn2井附近，延伸长度约9.5 km，与No58地层切片相比，qn8井已由水道侧翼到水道中心位置，表明水道向西迁移，而且由于大规模湖侵，砂体分布较为局限。在No17地层切片上(图 5d)，对应于沙一上中下部沉积，为高水位沉积期，水动力强，而且物源供给充分，重力流砂体发生大规模进积，发育三条辫状水道，水道较沙一中沉积期明显向东迁移，规模较大，水道砂体形成分叉，合并，连片、并行等多种组合形式。

与常规的沉积相研究结果相比，地层切片更为精细地刻画了重力流水道在不同时期的分布范围与演化规律。切片显示，随着时间的推移，水道主体经历了由东而西、再由西而东的演化规律。结合沙一中早晚期的微古地貌分析(图 4b,c)，揭示湖盆古地貌沙一中沉积早期较平坦，沉积中心在偏东部qn6—zh48井区，到沙一中沉积晚期南大港等控盆断层的强烈活动，导致湖盆最低洼处向西迁移，到沙一上亚段沉积时期，由于构造活动减弱及填平补齐作用，湖盆底形又复平坦的过程^[11]。

在反演剖面图上，重力流水道单砂体形态表现为孤立状，强振幅的特征，单河道边界清晰，垂向上叠置迁移，河道间普遍存在泥质隔层(图 6a)。在透视图上，重力流水道砂体规模巨大，为多期水道叠加复合而成，由南向北呈条带状展布(图 6b)。研究区在建立高频等时格架的基础上，应用地震沉积学技术和方法明确了沉积环境、沉积体系宏观展布特征，并对重力流水道砂体的内、外部特征进行精细刻画和描述。在此基础上，结合岩芯相、测井相的综合研究成果，得出沉积微相平面分布图。结果表明，沙一中亚段重力流水道砂岩体主要分布在张北、南大港断裂带形成的凹槽内和缓坡带辫状河三角洲砂体的前端，顺断层走向呈长条状或带状分布，单期水道宽度约1~2 km，长度为5 km以上，平均厚度3 m左右。按照重力流成因机制^[24]，歧南斜坡重力流水道可进一步划分为水道中心、水道边缘和水下漫溢等3种沉积微相(图 6c)。通过钻井检验，应用该方法进行沉积微相预测精度和符合率较高，较好的指导了沉积微相平面展布特征研究。

图 6 歧南斜坡沙一中亚段重力流水道展布特征 a.过zh47井波阻抗反演剖面；b.重力流水道立体透视图；c.沉积微相平面分布图 Fig. 6 Distribution of gravity-flow channel in mid-Es1 member in Qi′nan slope

图选项

研究黄骅坳陷大型陆相断陷湖盆发育的重力流水道砂体的分布对在湖相泥岩中寻找岩性油气藏具有一定的指导意义。

重力流水道呈条带状多期次发育，纵向上叠加、横向上连片，配合歧南斜坡背景，可在上倾方向形成岩性尖灭，或由于断层的遮挡形成岩性、构造—岩性圈闭。砂体为大段暗色泥质岩所夹持，生储盖组合等条件配置良好，又处于油气运移方向上，具备形成岩性油气藏的条件，是歧南斜坡重要的储油砂体类型。近几年研究区以沙一段为主要目的层系，采用地震沉积学的方法进行砂体预测和描述，圈定不同层系岩性圈闭的边界。在此基础上，分批部署钻探ch59、ch64x1井均获得工业油气流，砂体预测结果得到了钻井资料的进一步证实。重力流水道是一种重要的深水储层，qn6井的平均孔隙度为16.7%，ch59井的平均孔隙度为10.9%，且歧南斜坡沙一中亚段埋深处于晚成岩A亚期，次生孔隙较发育，相对优质储集层主要分布在水道中心微相^[17]。目前此类砂体并未大规模钻探，研究区是一个重要的油气勘探区域。

本文以歧南斜坡沙一段发育的带状重力流水道砂体为研究对象，在地震沉积学理论和方法的指导下，形成了以古地貌和地震属性约束物源方向，预测重力流水道的宏观展布；地层切片、反演与三维可视化相结合刻画砂体空间展布形态的断陷湖盆沟道型重力流地震沉积学识别方法和技术，较准确刻画了沙一段重力流水道储层空间展布和内部结构特征，取得了良好的地质效果，具有推广应用价值。同时，重力流水道的发现对深水重力流沉积体系研究具有重要意义，且可为研究区的油气勘探提供新的油气勘探目标。