沉积地质学发展新航程——地震沉积学*
朱筱敏1, 董艳蕾1, 曾洪流2, 黄捍东3, 刘强虎4, 秦祎1, 叶蕾1
1 中国石油大学(北京)地球科学学院,北京 102249
2 The University of Texas at Austin,Bureau of Economic Geology,Jackson School of Geosciences,Austin,Texas,USA 78713-8924
3 中国石油大学(北京)非常规研究院,北京 102249
4 中国地质大学(武汉)资源学院,湖北武汉 430074

第一作者简介 朱筱敏,男,1960年生,中国石油大学(北京)地球科学学院教授,主要从事沉积地质学和层序地层学教学和科研工作。E-mail: xmzhu@cup.edu.cn

摘要

地震沉积学是在地震地层学和层序地层学基础上发展起来的、沉积地质学与地球物理学相互交叉的新兴学科,在刻画薄层砂体、表征储集层和预测岩性圈闭等方面发挥了重要作用。在表述地震沉积学发展简史的基础上,介绍了地震沉积学的基本概念、基本原理和主要研究流程;在建立精细层序地层格架基础上,重点阐述了 90°相移、分频处理与频谱分解、属性优选、地层切片、地球物理反演、地震岩性解释和 RGB(红色、绿色、蓝色)地震属性融合等地震沉积学关键技术以及地震沉积相实例分析。认为中国地震沉积学应该围绕不同类型沉积盆地开展油气勘探阶段地震地貌学(地震沉积相模式)和复杂(多样)地震岩性学综合研究,在油气开发阶段,充分利用 RGB地震属性融合等多种地球物理方法技术,开展多种成因类型砂体构型和储集层定量表征研究。综合研究表明,地震沉积学是沉积地质学发展的新航程。

关键词: 地震沉积学; 90°相移; 地层切片; 研究流程; 发展趋势
中图分类号:P512.2 文献标志码:A 文章编号:1671-1505(2019)02-0189-13
New development trend of sedimentary geology: Seismic sedimentology
Zhu Xiao-Min1, Dong Yan-Lei1, Zeng Hong-Liu2, Huang Han-Dong3, Liu Qiang-Hu4, Qin Yi1, Ye Lei1
1 College of Geosciences,China University of Petroleum(Beijing),Beijing 102249,China
2 The University of Texas at Austin,Bureau of Economic Geology,Jackson School of Geosciences,Austin,Texas,78713-8924 USA
3 Institute of Unconventional Research,China University of Petroleum(Beijing),Beijing 102249,China
4 College of Resources,China University of Geosciences(Wuhan),Wuhan 430074,China;

About the first author Zhu Xiao-Min,born in 1960,is a professor and Ph.D. supervisor in China University of Petroleum(Beijing). He is mainly engaged in sedimentology and sequence stratigraphy. E-mail: xmzhu@cup.edu.cn.

Abstract

Seismic sedimentology is an emerging interdisciplinary subject between sedimentary geology and geophysics developed on the basis of seismic stratigraphy and sequence stratigraphy,which exerts a great influence on delineation of thin-bed sand-bodies,reservoir characterization,prediction of lithographic traps. The fundamental concepts,principles and main work-flows of seismic sedimentology have been well introduced based on its brief development history,several critical technologies regarding seismic sedimentology namely as 90°phase shift,frequency division processing,spectrum decomposition,attribute optimization,stratal slicing,geophysical inversion,seismic lithological interpretation,and RGB blending of seismic attributes as well as case study of seismic sedimentary facies have been well clarified within detailed sequence stratigraphic frameworks. Discussion shows that seismic geomorphological(seismic sedimentary facies model)and complicated seismic lithological study during the stage of hydrocarbon exploration should be conducted in sedimentary basins with various type when concerning seismic sedimentology in China,and the architecture of sand-bodies with different genesis and quantitative reservoir characterization are well implemented with RGB blending of seismic attributes and many other geophysical technologies during the stage of hydrocarbon exploitation,which will be the new development trend of sedimentary geology.

Key words: seismic sedimentology; 90° phase shift; stratal slicing; work-flows; development trend
1 概述

自从20世纪70年代美国学者Peter Vail 和他的同事们开创地震地层学以来, 地震地层学在区域油气勘探中发挥了重要作用(Vail et al., 1977)。它主要根据海平面周期性变化形成沉积层序的理论观点, 利用地震相的“ 相面” 或图形分析手段来研究地震同相轴组合特征(主要是地震反射内部结构和外部形态), 通过地震相的识别和组合特征去恢复盆地和区带规模的沉积体系时空展布以及可能的生储盖组合。地震地层学解释是基于沉积地质模型驱动的, 它主要在盆地或区带尺度确定不同类型的沉积体系, 但是受到有限的垂向分辨率限制, 只适用于厚层层序划分、盆地分析以及早中期油气资源勘探评价, 但对区带和油藏尺度的沉积单元和砂体研究精度尚存不足。随着当今三维地震技术开发创新和地震资料品质的明显提高(主要是地震分辨率提高), 人们对较小尺度地质体的地震资料解释精度和高分辨率地下沉积地质学研究的需求愈来愈高, 因而在少井区地震资料越来越多地被用于薄层砂体刻画、有利储集层预测、储集层精细表征和含油气区块的综合评价, 在这种背景下, 结合层序地层学基本原理和方法, 即诞生了新的学科— — 地震沉积学(Zeng et al., 1998a, 1998b; Zeng, 2018)。

地震沉积学的形成和发展得益于20世纪70年代三维地震技术的快速发展和工业化应用。Dahm和Graebner(1979)首次在地震时间切片上看到了反映曲流河道古地貌特征的高分辨率振幅影像, 展示了用地震资料水平成像技术直接显示沉积体系古地形和古地貌的先进性。随后, Brown等(1981)首先发现使用沿层地震成像方法可以从三维地震资料中提取高分辨率、具有类似卫星图片效果的沉积图像。大量研究表明, 地震地貌学能够强有力地展示地下沉积成像(Posamentier, 2001)。

20世纪80年代, 地球物理资料解释软件公司(例如Landmark, GeoQuest等)开发出层位(顺层)切片(Horizon Slice)技术, 改善了地震沉积相成像水平, 被工业界广泛采用。1998年, 美国得州大学(Austin)曾洪流、Backus和Henry教授等在《Geophysics》上发文首次提出“ 地震沉积学” 一词, 认为地震沉积学是通过地震岩性学、地震地貌学的综合分析, 研究地层岩性、沉积成因、沉积体系和盆地充填历史的学科。该新学科探讨了地震反射等时性和穿时性, 逐渐认识到地震反射并非普遍具有年代地层意义, 只有少数反射同相轴(如最大洪泛面附近的反射界面)具有可靠的等时意义, 而大多数地震反射的产状和等时性受控于地震资料频率。同时, 讨论了地震沉积学基本原理、主要研究方法和手段以及实际应用(Zeng et al., 1998a, 1998b)。

近20年来, 美国AAPG(美国石油地质学家协会)、SEG(地球物理学家协会)年会和IAS(国际沉积学家协会)国际会议均将地震沉积学作为重要论题开展学术交流。2005年和2010年先后在美国Houston(休斯顿)召开了地震沉积学国际会议, 推动了地震沉积学核心技术, 即地震岩性学和地震地貌学的发展。目前, 地震沉积学涉及到的地震地貌学和地震岩性学等理论、方法和技术在国际上得到广泛应用和不断发展, 已经掀起了地震沉积学的研究热潮, 国外已有许多学者在北美、西非和南亚等含油气盆地开展了一系列的(定量)地震沉积学研究, 并在油气勘探和开发方面取得了显著的效果(Schlager, 2000; Posamentier, 2001; Posamentier and Venkatarathnan 2003; Wood, 2007; Herná n et al., 2011; Janocko et al., 2013; Zeng, 2018)。2000年前后, 中国石油地质研究人员将先进的地震沉积学理论方法引入到中国陆相沉积盆地沉积体系和薄层沉积砂体研究中, 在陆相含油气沉积盆地沉积体系研究和薄层砂体预测、塔里木和四川盆地古生界、新元古界碳酸盐岩储集层研究以及油气精细勘探与开发方面发挥了特有的作用(董春梅等, 2006; 林承焰等, 2007; 朱筱敏等, 2009, 2011a, 2013, 2017; Zhao et al., 2011; 朱红涛等, 2011; Dong et al., 2014, 2015, 2017; Zhu et al., 2017, 2018; 岳大力, 2018a, 2018b; Yue et al., 2019)。

地震沉积学的提出使得沉积学家和油气勘探学家的研究视角产生了革命性变化, 从利用地震垂向分辨率开展常规地震相研究走向充分利用地震水平分辨率开展多种平面沉积体系研究(特别是利用地层切片进行的研究)。所以, 它不仅表现出沉积学与地球物理学学科交叉的先进性, 并且在石油工业多类储集层识别预测方面表现出实用性, 现已成为沉积地质学和石油工业的主要研究热点, 被认为是20世纪末至21世纪初石油工业重要的科技革新和沉积地质学发展的新航程。中国发育多种类型的沉积盆地, 具有多成因类型的薄层砂体和碳酸盐岩储集层, 这些薄层砂体和碳酸盐岩储集层现已成为重要的油气精细勘探和开发剩余油的对象。如何勘探开发薄层砂体或多种成因类型储集层中的油气资源, 地震沉积学将会发挥难以替代的重要作用。地震沉积学方法可识别利用地震垂向分辨率难以识别的薄层砂体(小于1/4波长, 可识别厚度3~5 m的砂体), 表述薄层砂体的形态和分布特征、重建不同类型储集层分布和沉积体系的沉积过程。显然, 在当今烃类勘探开发更多地聚焦于小尺度、薄层和地层岩性圈闭目标的阶段, 地震沉积学将在识别薄层砂体和碳酸盐岩有利储集层、岩性地层圈闭等方面发挥强有力的理论和技术支撑作用。

2 基本原理和研究程序

地震沉积学是在地震地层学和层序地层学基础上发展起来的一门新兴交叉学科。美国得州大学(Austin)曾洪流教授和Backus、Henry等在《Geophysics》上发表利用地震资料制作地层切片的论文(Zeng et al., 1998a, 1998b), 文中认为地震沉积学是利用地震资料来研究沉积岩及其形成过程的一门学科。随后, Schlager(2000)等认为, 地震沉积学是基于高精度地震资料、现代沉积环境和露头古沉积环境模式的联合反馈来识别沉积单元的三维几何形态、内部结构和沉积过程的一门学科。Posamentier(2001)提出了地震地貌学(Seismic Geomorphology)的概念。后来, Zeng 和 Hentz(2004)及其同事提出地震沉积学详细定义, 指出地震沉积学是用地震资料研究沉积岩和沉积作用的一门学科, 在当前条件下体现为地震岩性学和地震地貌学的综合。朱筱敏等(2011b)定义地震沉积学是以现代沉积学、层序地层学和地球物理学为理论基础, 利用三维地震资料及地质资料, 经过层序地层、地层切片、地震属性分析、岩心岩性和沉积相刻度研究, 确定地层岩石宏观特征、砂体成因、沉积体系发育演化、储集层质量及油气分布的地质学科(包含地震岩性学和地震地貌学)。曾洪流(2011)又对地震沉积学概念进行了深化, 认为地震沉积学是通过地震岩性学和地震地貌学的综合分析, 研究地层岩性、沉积成因、沉积体系和盆地充填历史的学科。近期, Zeng(2018)研究认为, 地震沉积学是基于地震信号响应薄层沉积单元的认识而开展的数据驱动研究, 其与沉积单元厚度、岩性— 波阻抗模型、子波相位和频率有关。地震沉积学基于地震岩性学和地震地貌学开展研究, 编制地震岩性地貌相图, 进而确定沉积环境和储集层性质。

地震沉积学研究基于下列2个基本原理: 大多数沉积体系(90%以上)具有宽度远远大于厚度的特征(Galloway and Hobday, 1983); 在地震垂向分辨率与横向分辨率相当的情况下, 可利用横向分辨率(如地层切片)将利用地震垂向分辨率难以识别的薄层砂体识别出来(Sheriff, 2002)。

三维地震勘探技术的快速发展为沉积砂体成像提供了无与伦比的基础。在有些情况下, 以地震数据层面属性为基础, 通过90° 相位调整、大量层面属性研究, 优选出振幅、方位角、相似性和方差等多种与沉积体系层面几何形态有关的地震属性。结合层面三维可视化技术、RGB融合技术和地质历史时期构造形态恢复等技术, 结合源汇系统新概念, 展现不同地质历史时期的沉积体系形态特征, 依据沉积砂体形态和沉积模式对地震平面属性资料(地层切片)进行沉积地貌分析。当地貌分析与地震平面属性分析结合后, 可以得到重要的沉积地质信息。最后, 通过岩心、钻测井和实验数据的综合研究和刻度(地震岩性学), 可识别厚3~5im的薄层砂体, 阐明沉积体系的发育演化等地质特征。

地震岩性学和地震地貌学是地震沉积学的主要支柱。应用地震岩性学可将三维地震数据体转换为测井岩性数据体。在岩性数据体中, 各井点处的岩性测井(如GR和SP)以很小的允许误差与井旁地震道建立关系, 以确保储集层段井数据与地震数据的最佳匹配。应用地震地貌学, 结合现代沉积地貌和古地形特征, 可以将地震数据进一步转换成沉积相图。地震地貌学分析的关键在于从地震数据中提取有明确地质意义的多种地震属性, 使反射界面的地貌学特征在地震属性平面图上直观成像, 进而确定沉积体系类型和时空演化。

现今, 地震沉积学已基本形成了较为规范的工业化研究流程: (1)利用层序地层学原理, 综合多种地质和地球物理资料建立高精度等时地层格架; (2)进行子波相位频率调整(常常是90° 相位调整), 以建立薄层砂体与地震反射同相轴之间的对应关系, 即测井岩性数据体; (3)根据地层产状, 优选地震数据的切片方法, 通常从时间切片、沿层切片(顺层切片)和地层切片中选用地层切片方法, 来研究非水平、非等厚地层的平面地震属性特征; (4)根据地震调谐厚度与地震频率之间的关系, 开展分频地震参数处理和统计, 优选能够反映薄层砂体的最佳频率; (5)利用岩心资料, 通过地震岩性学研究(地震属性与岩性关系、地震模型和地震反演等), 赋予地震属性的岩性和沉积体系含义, 进而开展井震对比以及多井对比检验, 建立地震属性平面特征与关键井岩性之间的良好对应关系; (6)综合利用地震属性的平面地貌学特征, 基于岩心刻度, 开展其与沉积体系关系研究, 确定沉积体系类型和砂体形态; (7)多层段高精度地震沉积学研究, 建立基于源汇系统的沉积体系和沉积砂体演化模式, 恢复沉积体系和沉积砂体演化历史; (8)开展储集层质量、成藏要素、油气富集程度与沉积体系、砂体类型之间关系研究, 结合烃类检测等方法, 分析岩性圈闭勘探有利地区, 为油气勘探开发提供精细的沉积砂体分布演化格架(曾洪流等, 2012; 朱筱敏等, 2017)。

3 主要研究方法和技术
3.1 90° 相位调整

常规地震处理的最终成果多是零相位的地震数据, 解释人员常用这类地震数据开展构造和沉积地质解释。零相位地震数据用作地质解释的优点包括子波对称、中心瓣(最大振幅)与反射界面一致以及具有较高的分辨率(Brown et al., 1981)。但是, 只有当地震反射是来自一个单一的界面(如海底、主要不整合面和厚层块状砂岩的顶面等)时, 零相位地震数据的上述优点才是真实的。对于薄层砂体来说, 使用零相位地震数据进行解释是不够准确的, 故Zeng和Backus(2005a, 2005b)提出了90° 相位子波调整方法技术来建立90° 相位地震数据体与薄层砂体之间对应关系。Zeng和Backus(2005a, 2005b)以美国Starfak油田新近系地震剖面为例, 指出在零相位地震数据体中砂体与地震同相轴之间没有直接的对应关系, 地震极性和振幅既不能很好地指示岩性, 又不是岩相位置和形态的可信参考点, 因此标准的零相位地震数据不适合做薄层岩性解释(图 1-a)。

图 1 美国Starfak油田新近系地震剖面90° 相位调整(据Zeng和Backus, 2005a, 2005b; 有修改)
a— 标准零相位地震剖面以及提取的子波; b— 90° 相位地震剖面以及提取的子波。曲线为经测井资料计算的砂泥岩含量, 曲线左偏砂岩, 曲线右偏泥岩; 井号下方黑色直线为经过相应井位的三维地震资料拼接线
Fig.1 90° phase shift on seismic profile of the Neogene in Starfak Oilfield, USA (modified from Zeng and Backus, 2005a, 2005b)

采用90° 相位子波处理的地震数据就可以克服零相位子波数据的不足, 当地震响应的主瓣(最大振幅)经过移动与薄层砂岩中心相对应时, 主要的地震同相轴与地质上限定的砂岩层对应, 解释工作(特别是地层切片解释)就变得更加准确和相对容易。图 1-b所示的地震剖面是零相位地震数据经90° 相位调整之后形成的, 在重新处理的地震剖面中, 几乎所有井的砂岩层都对应于地震波谷(红色), 不确定性极小。

3.2 地层切片

地震沉积学(地貌学)成像的一个关键条件是解释人员能够在沉积界面 (地质时间界面)上提取地震振幅信息, 以获取整个地震测网内的沉积体系图像。目前, 常用的反映沉积体平面地震属性的切片技术主要包括时间切片、沿层切片和地层切片(Zeng et al., 1998a, 1998b)或按比例切片(Proportional Slicing)(Posamentier, 2001)(图 2)。在实际研究中, 应该选取最适合于特定构造和地层产状的某一种切片方法: 假如地层是席状且平卧的, 时间切片(Time Slice)就能适用; 假如地层是席状但并非平卧状态(最好等厚), 则沿层切片(Horizon Slice)比较合适; 假如地层既不是席状的也不呈平卧状(特别是不等厚), 则必须选用地层切片(Stratal Slice)。地层切片比时间切片和沿层切片更接近于地质时间界面。

图 2 不同地震界面切片方法的对比(据Zeng和Hentz, 2004; 有修改)
A— 时间切片; B— 沿层切片; C— 地层切片
Fig.2 Slicing methods of different seismic boundaries(modified from Zeng and Hentz, 2004)

地层切片(Zeng et al., 1998a, 1998b)是从地质时代界面(等时界面)上提取的地震属性(经常是振幅属性)。该方法不仅追踪了等时地质界面, 而且能够有效地反映地震探区中常见沉积体系的地貌特征。地层切片能很容易地拾取振幅型或结构异常型沉积体系, 如现代海底峡谷和盆地扇、嵌入厚层页岩的河道、天然堤体系、硅质碎屑岩系内的薄灰岩层, 这些部位的等时地震异常容易识别和追踪。当然, RGB融合技术可有效地刻画不同厚度沉积砂体的地貌特征(见下文)。

3.3 分频处理与频谱分解

地震沉积学与地震地层学的最大不同在于它认识到地震同相轴既不简单地反映等时界面, 也不单纯反映岩性界面, 而是受到地震资料频率的控制(Zeng et al., 1998; Zeng and Kerans 2003; Zeng and Hentz, 2004)。不同频段的地震数据反映的地质信息等时性是不同的。低频资料中地震反射同相轴更多地反映复合岩性界面信息, 而高频地震资料中反射同相轴更多地反映等时沉积界面信息。

基于“ 地震资料的频率成分控制了地震反射同相轴的倾角和内部反射结构” 这一认识(Sheriff, 2002; 董春梅等, 2006; 林承焰等, 2007), 针对不同的地质目的层段使用不同频段的地震数据可有效刻画不同沉积厚度的地质体。地震资料中连续的频率变化本身蕴含了丰富的地质信息, 不同级别的地质层序体对应着地震剖面上的不同频率特征, 因此, 可采用分频解释方法充分获得不同沉积厚度的地质体信息。

近期兴起的频谱分解技术是一种新的叠后地震数据处理和解释技术, 也是进行地震沉积学研究的强有力技术手段。它是利用小波变换或傅立叶变换, 将地震数据由时间域变换到频率域, 然后, 沿层或固定时窗把地震反射波中的各频率成分对应的调谐能量识别出来而形成对应频率的能量异常图, 进而有效检测不同沉积厚度的地质体。

频谱分解技术在理论上主要是依据薄层反射的调谐原理。对于厚度小于四分之一波长(λ /4)的地震薄层而言, 在时间域, 随着薄层厚度的增加, 地震反射振幅逐渐增大; 当薄层厚度增加至λ /4时, 反射振幅达到最大值; 然后, 随着地层厚度大于λ /4, 反射振幅逐渐减小。时间域的最大反射振幅值, 对应着频率域的最大振幅能量值(图 3)。

图 3 振幅与厚度关系在不同频率段内的响应(据于建国等, 2006)Fig.3 Relation between amplitude and thickness in different frequency(from Yu et al., 2006)

由于地层对地震信号的调谐和吸收作用, 不同的地质目标对地震资料的不同频率成分的敏感程度不同。例如, 深层目标和单层厚度较大的韵律层突出低频成分, 浅层目标和单层厚度小的韵律层则突出高频成分。实际工作时, 可根据研究区资料情况将原始地震数据体分为若干个单独频率的数据体。对某一未知厚度的砂体, 先用各个频率扫描, 总有一个频率的地震波能满足其调谐条件, 利用该调谐地震波的速度和频率反算储集层的厚度, 即可达到识别薄储集层的目的。因此, 利用地震信号的特定频率或频带信息来突出地质目标的成像效果一直是地球物理勘探领域的发展方向。

3.4 RGB地震属性融合

地震资料最早为黑白双色显示, 为了更好地应用地震资料, Balch(1971)将地震数据进行彩色显示。自此, 利用彩色图像的视觉效果直接识别地震异常成为地震资料解释的有效手段(李婷婷等, 2015)。近几十年来, 勘探地球物理学家从地震数据中提取出数百种地震属性(如振幅类、频率类、相位类、波形类、构造类和叠前类属性以及谱分解类属性等)(丁峰等, 2010)。对于多个地震属性来说, 采用单个属性逐一彩色显示方法不能明显地反映整体趋势及某些隐蔽的地质特征。为了解决这一问题, 提高对局部地下地质异常体的识别能力并且更加直观地分析地震属性图, 研究人员在地震属性融合中引入基于色彩模型的多属性融合技术(RGB)(Guo et al., 2008)。

地震多属性颜色融合, 指利用不同的数字色彩模型来强化地质特征, 从而有效地实施地质目标的定性解释。对同一地质目标, 不同属性所反映的信息是有所差异的。多属性融合的关键就是充分利用多个待融合属性进行信息互补, 形成一个新的信息集合。这种技术是一种基于视觉的属性分析方法, 其原理是, 将多个地震属性通过主成分分析(PCA)技术进行降维, 并将主分量按特征值由大到小排序, 取前3个(或4个)主分量利用RGBA(Red-Green-Blue-Alpha)颜色融合原理获得1张融合图; 再结合实际地质资料, 在融合图像上依据颜色的区域性和突变异常等视觉特征, 进行地质目标和沉积厚度解释(Stark, 2006; 丁峰等, 2010)。

RGB属性融合技术的优点是能够充分利用地震属性中蕴含的构造和岩性信息, 从而克服了单一地震属性显示不足和单属性色彩不能突出区域异常的缺点(李艳芳等, 2009), 提高了从多属性中提取地质体的能力, 使图像显示更加清晰, 具有特征明显、细节丰富、信息量高和多属性联合显示的特点。但是, RGB(RGBA)融合技术也存在自身的缺陷, 该技术必须使用3种(RGB)或4种(RGBA)属性进行融合, 不适用于4种以上的地震属性融合(李婷婷等, 2015)。目前, RGB融合技术多用于河道检测(张驰等, 2013)、小断层识别、地质体异常识别和储集层预测等方面(成荣红等, 2013)。

3.5 地球物理反演

地球物理反演的发展已有百年历史, 是预测砂体分布和解决油气勘探问题最常用的方法技术(Herglotz, 1907; Backus and Gilbert, 1967)。比如, 波阻抗(道积分)直接反演方法(Lindseth, 1979)、最小平方反演方法(Draper and Smith, 1981)、广义线性方法(Cooke and Schneider, 1983)。后来, McAulay(1985)将广义线性方法用于叠前地震反演。还有适用于井数较少或无钻井资料地区的稀疏脉冲反演算法(Debeye and Riel, 1990), 出现了以非线性反演理论为基础的多种反演算法, 如遗传算法、模拟退火法和人工神经网络算法等, 将多种反演方法交叉运用如多尺度波阻抗反演和地质统计学反演等(Bortoli et al., 1992; Hass and Dubrule, 1994; Rothman, 1998)。

最近, Huang等(2016)提出地震相控非线性混沌反演方法, 在利用井资料和地震相的约束下, 对深层薄互层岩性油气藏进行了非线性混沌反演, 最终将地震道反演算法表示为:

Δzk=[Gk-1TGk-1+σsk-12(σzk-12)-1I]-1Gk-1T(s-sk-1)zk=Zk-1+Δzk

其中, Δ zk为第k次迭代的波阻抗的修改增量; G为正演算子矩阵; λ k-1= σnk-12/ σzk-12为阻尼因子, σsk-12为第k-1次迭代的合成地震记录与真实地震记录残差的方差, 表征噪音; σzk-12是第k-1次迭代前后模型残差的方差, 表征模型摄动量; I为单位矩阵; s为地震道数据; sk-1为第k-1次迭代的合成地震记录; zk为第k次迭代的波阻抗; 上标T表示转置。

使用该方法反演结果突破了常规地震分辨率的限制, 可根据迭代过程具体状态来自适应地修改λ 的值, 从而调节噪音与分辨率, 使之达到最佳平衡。图 4为该方法在塔里木盆地应用效果。图 4左侧红色曲线代表TZ4-18-8井合成地震记录轨迹, 右侧红色曲线代表TZ4-18-8井声波测井曲线。从图中黑色虚线圆圈所圈定范围可以看出, 传统反演方法不能充分利用弱地震信号, 所得反演结果分辨率较低。混沌反演能够比常规反演揭示更多细节, 尤其在薄层砂体刻画上, 展现出良好的效果, 反演结果与实测数据吻合较好, 突出了界面和岩性体, 这种综合方法还提高了沉积界面和岩性体的分辨率, 通过实践应用, 砂岩储集层预测吻合率达到了85%以上(Huang et al., 2009; Luo et al., 2018)。通过此方法, 可以针对复杂油气藏进行有效分析, 从而增加地震沉积学应用的可信程度, 使其应用效果得到进一步改善。

图 4 塔里木盆地某原始地震剖面与反演剖面对比
a— 假定λ 为常数且使用传统反演方法得到的反演结果剖面; b— 利用混沌反演方法得到的反演结果剖面(据Huang等, 2016)
Fig.4 Original and inversion seismic profiles from Tarim Basin

3.6 地震岩性解释

利用地震沉积学地层切片技术, 即沿沉积界面(地质时间界面)提取地震属性数据, 可有效反映地震工区不同沉积层序沉积体系的展布特征, 高分辨率刻画沉积砂体。在实际工作中, 要有效确定沉积相带和薄层砂体特征, 需要充分开展地震岩性学研究工作。

在地震资料预处理过程中, 已经对地震数据体进行了90° 相位调整, 所以地震反射同相轴被赋予了地质概念, 地质工作者解释的地震剖面基本近似于波阻抗剖面。但是这些同相轴究竟代表的是什么岩性, 需要对岩性与波阻抗关系进行探讨。

图 5为准噶尔盆地P601-4井伽马— 波阻抗交会图, 图中横坐标为伽马测井数值, 纵坐标为波阻抗数值, 不同颜色的数据点是按照密度值的高低进行区分的, 红色数据点代表砂岩含量高值, 蓝色数据点代表泥岩含量高值, 绿色数据点代表砂泥岩过渡类型。从图中可以直观看出不同颜色数据点在不同域值范围内的分布, 砂岩多为高阻抗, 与泥岩的分界线为20 000g/cm3· m/s左右; 泥岩多为低阻抗, 这为地层切片沉积地质解释提供了可靠依据。

图 5 准噶尔盆地P601-4井伽马— 波阻抗交会图Fig.5 Plot of GR vs. wave impedance of Well P601-4 from Junggar Basin

此外, 还需开展关键井岩心的岩性序列分析, 通过精确合成地震记录标定, 才能对地层切片进行沉积微相和砂体分布解释。下面以准噶尔盆地车排子地区白垩系为例开展研究。车排子地区白垩系发育的岩性主要为泥岩和粉砂质泥岩, 夹薄层粉砂岩。为了准确刻度砂岩和泥岩, 选取关键井P605井进行合成记录标定, 明确白垩系地震剖面中强的负振幅区(红轴)代表砂岩, 正的振幅值(黑轴)对应泥岩, 而弱的负振幅值区代表薄层砂岩或泥质砂岩(图 6), 从而实现了地震岩性的相互标定。

图 6 准噶尔盆地P605井白垩系合成地震记录标定和地层切片Fig.6 Synthetic record and stratal slicing of the Cretaceous in Well P605, Junggar Basin

3.7 地震沉积相实例解释

地震沉积学实例研究来自渤海湾盆地沾化凹陷, 凹陷面积约2800 km2。研究区三合村洼陷位于沾化凹陷西南部, 南侧临近陈家庄凸起, 研究层段为古近系沙河街组三段(划分为3个三级层序), 沉积物源来自南侧陈家庄凸起, 形成湖泊— 扇三角洲沉积体系(图7)。根据研究区三维地震资料和数十口钻井资料的分布特征, 采用下列研究步骤, 开展地震沉积学研究: (1)建立联井基干地震剖面测网。测井— 地震联合对比建立高精度层序格架。(2)地震子波相位调整。对地震意义上的薄层(单层厚度小于1/4波长), 子波相位应调整到90° 。(3)追踪闭合地震地质等时标志层。如最大洪泛面凝缩层、平行不整合和特殊岩性层(煤和薄层灰岩等)。(4)估算地震分辨率。包括用频谱分析估计地震波有效频率范围和主频; 用井资料制作正演模型以确定薄层时间分辨厚度极限以及地震切片检测最小厚度极限。(5)地震频率调整。为实现高等级层序地层格架与地震资料的最佳配合, 应尽可能将地震沉积学最小作图单元调整到高等级层序或高等级层序体系域的平均厚度范围。(6)岩石物理关系分析。用岩心实验测定数据或关键井测井曲线来统计目的层段不同岩性间的波阻抗对应关系及极性/振幅对应关系; 确定用地震参数预测岩性的可行性。(7)地震参数筛选。试验多种地震参数以确定预测岩性和沉积相最佳地震属性参数或参数组合。(8)地层切片处理。用专用软件(如Recon StatalSlice、Rockstar等), 建立地层切片和高等级层序的对应关系。(9)地震沉积相分析。根据地层切片上显示的地震岩性学和地震地貌特征, 以及岩心— 测井相标定, 辅以剖面地震相特征, 解释沉积相、沉积环境和沉积体系分布。(10)储集层和石油地质评价。综合多学科研究成果, 预测砂体厚度、分布范围、储集层质量和地层岩性圈闭等。

地震沉积学综合研究表明, 沾化凹陷三合村洼陷沙三段SQ2沉积时期, 低位体系域(LST)和高位体系域(HST)都发育扇三角洲, 且二者朵体规模较大; 湖侵体系域(TST)扇三角洲朵体展布范围局限, 含泥率较高。不同体系域的扇三角洲前缘相带富砂, 是有利勘探地区(图 7)。

图7 沾化凹陷三合村洼陷沙三段SQ2不同体系域地层切片沉积相精细刻画Fig.7 Strata slicings of different tracts systems (SQ2) of the Member 3 of Shahejie Formation from Sanhecun sub-sag Zhanhua Sag

4 讨论与结论

地震沉积学可利用地震振幅信息和属性分析技术研究沉积岩石学、古地貌学、沉积体系、沉积砂体类型、碳酸盐岩储集层以及沉积演化历史等, 其解释是以各种地震信息的平面图为基础, 从地震数据中提取有明确地质意义的多种地震平面属性, 使反射界面特征(地貌学特征)和沉积体系在地震属性上直观成像, 并进行反射界面特征的地貌学分析和进行沉积体系展布解释。与地震地层学和层序地层学相比, 地震沉积学的优势表现为: (1)以层序地层格架为基础, 保证地层对比等时性; (2)充分利用了地震水平分辨率, 可识别地震垂向分辨率难以识别的薄层砂体, 分辨能力更高; (3)将沉积学与地震属性处理结合起来, 可预测砂体空间展布及其演化; (4)可利用储集层反演等方法技术, 预测储集层和流体性质以及有利的岩性油气圈闭勘探地区。

地震沉积学目前包括地震岩性学和地震地貌学, 可有效用于薄层储集层解释。在地震沉积学研究中, 地震资料90° 相位化及地层切片是2项关键技术, 可用于大多数地震资料条件, 也经济实用。

地震沉积学经过20年的发展, 已经在石油地质学术界和石油工业界得到了广泛应用, 理论和研究方法技术不断得到完善改进。目前, 国际上地震沉积学研究主要侧重于露头地貌与地下地震地貌的类比研究、沉积体系地貌形态演变的三维模拟、地震岩性学方法(时频分析、地震参数分析等)、地层切片方法的改进以及切片结果的显示(比如RGB融合技术)、定量地震地貌学以及复杂沉积层序(河流、三角洲、浅海陆棚、重力流及碳酸盐沉积等)地震沉积学和地震地貌学的研究等(曾洪流, 2011)。

中国含油气盆地油气资源丰富, 大量油气资源富存于薄层砂岩储集层以及碳酸盐岩储集层中, 显然地震沉积学对研究复杂沉积层序中沉积储集层(尤其是薄层砂体储集层)和地层岩性油气藏将发挥不可替代的作用。针对目前的常规和非常规油气资源勘探, 不仅要充分认识到, 中国地震沉积学应该考虑中国含油气盆地地质构造、不同演化阶段和盆地类型(拗陷、断陷和前陆)等构造背景, 通过露头和现代沉积研究, 创建不同沉积体系地貌学数据库, 建立受构造作用控制的在等时地层格架中的沉积体系和砂体发育模式, 创立各类陆相盆地的地震沉积相模式(地震沉积相特指地震分辨率尺度下代表特定沉积相的地震特征组合); 还要探索地震沉积学在非常规油气勘探(致密油、页岩油等)以及深层和超深层油气资源勘探中的应用方法技术。中国陆相含油气盆地岩性— 波阻抗关系复杂(砾岩、砂岩、泥岩、灰岩和煤层混积), 地层切片的岩性解释更具多解性, 90° 相位地震剖面的岩性校定通常比较困难, 应该创新发展适合陆相盆地的地震岩性学新方法。另外, 陆相盆地湖平面升降频繁, 物源多变, 砂、泥岩波阻抗关系受沉积旋回影响明显, 为提高地层切片的岩性解释水平, 应加强岩石物理理论和实验研究, 探索新的90° 相位地震剖面的岩性校定方法与地震反演、地震参数分析、时频分析方法等之间的相互关系。

地震沉积学理论方法不仅可用于勘探阶段薄层(厚层)砂体识别和岩性圈闭预测, 而且可以用于油气开发阶段砂体定量表征。在相对稀井网条件下, 充分利用地震沉积学方法, 挖掘地震信息可有效预测表征储集层构型。通过优选频段、融合分频数据体、提取并优选地震属性的方法, 可高精度预测不同沉积砂体分布。研究表明, 振幅类属性与砂体厚度相关性最好, 以最大峰值振幅为最佳, 相比原始地震属性, 分频融合得到的最大峰值振幅属性更好地刻画了多种成因砂体边界与砂体厚度分布, 地震正演模拟和分频融合地震属性的储集层预测方法可明显提高储集层预测精度, 可定量表征砂体厚度、分布和形态等, 分析预测储集层质量。

作者声明没有竞争性利益冲突.

参考文献
[1] 成荣红, 杨斌, 黄科, 徐文圣. 2013. 地震多属性融合成像技术在储层预测中的应用. 天然气技术与经济, 7(4): 29-32.
[Cheng R H, Yang B, Huang K, Xu W S. 2013. Application of seismic attributes-imaging to reservoir prediction. Natural Gas Technology and Economy, 7(4): 29-32] [文内引用:1]
[2] 丁峰, 年永吉, 王治国, 尹成, 古发明. 2010. 地震多属性RGBA颜色融合技术的应用研究. 石油物探, 49(3): 248-252.
[Ding F, Nian Y J, Wang Z G, Yin C, Gu F M. 2010. Application of seismic multi-attributes RGBA color blending. Geophysical Prospecting for Petroleum, 49(3): 248-252] [文内引用:2]
[3] 董春梅, 张宪国, 林承焰. 2006. 地震沉积学的概念、方法和技术. 沉积学报, 24(5): 699-704.
[Dong C M, Zhang X G, Lin C Y. 2006. Concepts, methodologies, and technologies of the seismic sedimentology. Acta Sedimentologica Sinica, 24(5): 699-704] [文内引用:2]
[4] 李婷婷, 王钊, 马世忠, 王昭, 袁子龙. 2015. 地震属性融合方法综述. 地球物理学进展, 30(1): 378-385.
[Li T T, Wang Z, Ma S Z, Wang Z, Yuan Z L. 2015. Summary of seismic attributes fusion method. Progress in Geophysics, 30(1): 378-385] [文内引用:2]
[5] 李艳芳, 程建远, 朱书阶, 王成. 2009. 基于RGB渲染技术的地震多属性分析技术. 煤炭学报, 34(11): 1512-1516.
[Li Y F, Cheng J Y, Zhu S J, Wang C. 2009. Seismic multi-attribute analysis based on RGB color blending technology. Journal of China Coal Society, 34(11): 1512-1516] [文内引用:1]
[6] 林承焰, 张宪国, 董春梅. 2007. 地震沉积学及其初步应用. 石油学报, 28(2): 69-71.
[Lin C Y, Zhang X G, Dong C M. 2007. Concepts of seismic sedimentology and its preliminary applications. Acta Petrolei Sinica, 28(2): 69-71] [文内引用:2]
[7] 于建国, 韩文功, 刘力辉. 2006. 分频反演方法及应用. 石油地球物理勘探, 41(2): 193-197.
[Yu J G, Han W G, Liu L H. 2006. Frequency-divided inversion and application. Petroleum Geophysical Exploration, 41(2): 193-197] [文内引用:1]
[8] 岳大力, 李伟, 王军, 王武荣, 李健. 2018a. 基于分频融合地震属性的曲流带预测与点坝识别: 以渤海湾盆地埕岛油田馆陶组为例. 古地理学报, 20(6): 999-1008.
[Yue D L, Li W, Wang J, Wang W R, Li J. 2018a. Prediction of meand ering belt and point bar recognition based onspectral decomposed and fused seismic attributes: A case study of the Guantao Formation, Chengdao Oilfield, Bohai Bay Basin. Journal of Paleogeography(Chinese Edition), 20(6): 999-1008] [文内引用:1]
[9] 岳大力, 胡光义, 李伟, 范廷恩, 胡嘉靖, 乔慧丽. 2018b. 井震结合的曲流河储层构型表征方法及其应用: 以秦皇岛32-6油田为例. 中国海上油气, 30(1): 99-109.
[Yue D L, Hu G Y, Li W, Fan T E, Hu J J, Qiao H H. 2018b. Meand ering fluvial reservoir architecture characterization method and application by combining well logging and seismic data: A case study of QHD32-6 oilfield. China Offshore Oil and Gas, 30(1): 99-109] [文内引用:1]
[10] 曾洪流. 2011. 地震沉积学在中国: 回顾与展望. 沉积学报, 29(3): 417-426.
[Zeng H L. 2011. Seismic Sedimentology in China: A Review. Acta Sedimentologica Sinica, 29(3): 417-426] [文内引用:2]
[11] 曾洪流, 朱筱敏, 朱如凯, 张庆石. 2012. 陆相坳陷型盆地地震沉积学研究规范. 石油勘探与开发, 39(3): 295-304.
[Zeng H L, Zhu X M, Zhu R K, Zhang Q S. 2012. Guidelines for seismic sedimentologic study in non-marine postrift basins. Petroleum Exploration and Development, 39(3): 295-304] [文内引用:1]
[12] 张驰, 朱博华, 刘培金, 李宇航, 张在金. 2013. RGB多地震属性融合技术在河道检测中的应用. 见: 中国地球物理2013—第二十分会场论文集. 北京: 中国地球物理学会, 712-713.
[Zhang C, Zhu B H, Liu P J, Li Y H, Zhang Z J. 2013. The application of RGB blending of multi-seismic attributes to channel detection. In: Geophysics in China Transaction of 20th Sessions. Beijing: China Geophysical Union, 712-713] [文内引用:1]
[13] 朱红涛, 杨香华, 周心怀, 李建平, 王德英, 李敏. 2011. 基于层序地层学和地震沉积学的高精度三维沉积体系: 以渤中凹陷西斜坡BZ3-1区块东营组为例. 地球科学(中国地质大学学报), 36(6): 1073-1084.
[Zhu H T, Yang X H, Zhou X H, Li J P, Wang D Y, Li M. 2011. High resolution three-dimensional facies architecture delineation using sequence stratigraphy, seismic sedimentology: Example from Dongying Formation in BZ3-1 block of western slope of Bozhong Sag, Bohai Bay Basin. Earth Science: Journal of China University of Geosciences, 36(6): 1073-1084] [文内引用:1]
[14] 朱筱敏, 刘长利, 张义娜, 彭伟民. 2009. 地震沉积学在陆相湖盆三角洲砂体预测中的应用. 沉积学报, 27(5): 915-921.
[Zhu X M, Liu C L, Zhang Y N, Peng W M. 2009. On seismic sedimentology of lacustrine deltaic depositional systems. Acta Sedimentologica Sinica, 27(5): 915-921] [文内引用:1]
[15] 朱筱敏, 董艳蕾, 胡廷惠, 罗琪, 王敏, 杨道庆, 罗家群, 刘桂兰. 2011a. 精细层序地层格架与地震沉积学研究: 以泌阳凹陷核桃园组为例. 石油与天然气地质, 32(4): 615-624.
[Zhu X M, Dong Y L, Hu T H, Luo Q, Wang M, Yang D Q, Luo J Q, Liu G L. 2011. Seismic sedimentary study of fine sequence stratigraphy framework: A case study of the Hetaoyuan Formation in the Biyang Sag. Oil & Gas Geology, 32(4): 615-624] [文内引用:1]
[16] 朱筱敏, 曾洪流, 董艳蕾, . 译. 曾洪流, 等. 著. 2011b. 地震沉积学(译文集). 北京: 石油工业出版社.
[Zhu X M, Zeng H L, Dong Y L, et al. (Translators). Zeng H L, et al. (Authors). 2011. Collected Works of Seismic Sedimentology. Beijing: Petroleum Industry Press] [文内引用:1]
[17] 朱筱敏, 赵东娜, 曾洪流, 孙玉, 朱如凯, 黄薇, 朱世发. 2013. 松辽盆地齐家地区青山口组浅水三角洲沉积特征及其地震沉积学响应. 沉积学报, 31(5): 889-897.
[Zhu X M, Zhao D N, Zeng H L, Sun Y, Zhu R K, Huang W, Zhu S F. 2013. Sedimentary characteristics and seismic sedimentologic responses of shallow-water delta of Qingshankou Formation in Qijia area, Songliao Basin. Acta Sedimentologica Sinica, 31(5): 889-897] [文内引用:1]
[18] 朱筱敏, 曾洪流, 董艳蕾. 2017. 地震沉积学原理与应用. 北京: 石油工业出版社.
[Zhu X M, Zeng H L, Dong Y L. 2017. The principle and application of seismic sedimentology. Beijing: Petroleum Industry Press] [文内引用:2]
[19] Backus G E, Gilbert J F. 1967. Numerical applications of a formalism for geophysical inverse problems. Geophysical Journal International, 13(1-3): 247-276. [文内引用:1]
[20] Balch A H. 1971. Color sonograms: A new dimension in seismic data interpretation. Geophysics, 36(6): 1074-1098. [文内引用:1]
[21] Bortoli L J, Alabert F, Journel A G. 1992. Constraining Stochastic Images to Seismic Data. Proceedings of the International Geostatistics Congress, 33: 526-560. [文内引用:1]
[22] Brown, A R, Dahm C G, Graebner R T. 1981. A stratigraphic case history using three-dimensional seismic data in the Gulf of Thailand . Geophysical Prospecting, 29: 327-349. [文内引用:2]
[23] Cooke D A, Schneider W A. 1983. Generalized linear inversion of reflection seismic data. Geophysics, 48(6): 665-676. [文内引用:1]
[24] Dahm C G, Graebner R J. 1979. Field development with three-dimensional seismic methods in Gulf of Thailand -a case history. Offshore Technology Conference Paper 3657: 2591-2595. [文内引用:1]
[25] Debeye H W J, Riel P. 1990. Norm deconvolution. Geophysical Prospecting, 38(4): 381-403. [文内引用:1]
[26] Dong Y L, Zhu X M, Xian B Z. 2014. Mapping sediment-dispersal characteristics using seismic geomorphology: Late Paleogene to Neogene, Qinan Sag, Huanghua Depression, China. Marine and Petroleum Geology, 54: 180-197. [文内引用:1]
[27] Dong Y L, Zhu X M, Xian B Z. 2015. Seismic Geomorphology Study of the Paleogene Hetaoyuan Formation, Central-South Biyang Sag, Nanxiang Basin, China. Marine and Petroleum Geology, 64: 104-124. [文内引用:1]
[28] Dong Y L, Zhang M Y, Zhu X M, Wei M P. 2017. Seismic geomorphology and depositional system of delta and terminal fan: A case study of the Neogene Shawan Formation in the Chepaizi Uplift, Junggar Basin, China. Marine and Petroleum Geology, 83: 362-381. [文内引用:1]
[29] Draper N R, Smith H. 1981. Applied regression. Analysis, 2: 17-33. [文内引用:1]
[30] Galloway W E, Hobday D K. 1983. Terrigenous Clastic Depositional Systems: Applications to Petroleum, Coal, and Uranium Exploration. New York: Springer-Verlag, 1-423. [文内引用:1]
[31] Guo H, Lewis S, Marfurt K J. 2008. Mapping multiple attributes to three-and four-component color models-a tutorial. Geophysics, 73(3): W7-W19. [文内引用:1]
[32] Hass A, Dubrule O. 1994. Geostatistical inversion-a sequential method of stochastic reservoir modeling constrained by seismic data. Geophysics, 48(12): 561-569. [文内引用:1]
[33] Hernán M R, Posamentier H W, Bhattacharya J P. 2011. Seismic geomorphology and high-resolution seismic stratigraphy of inner-shelf fluvial, estuarine, deltaic, and marine sequences, Gulf of Thailand . AAPG Bulletin, 95(11): 1959-1990. [文内引用:2]
[34] Huang H D, Zhang R W, Luo Q, Zhao Di, Peng Y M. 2009. Subtle trap recognition based on seismic sedimentology: A case study from Shengli Oilfield. Applied Geophysics, 6(2), 175-183. [文内引用:1]
[35] Huang H D, Yuan S Y, Zhang Y T, Zeng J, Mu W T. 2016. Use of nonlinear chaos inversion in predicting deep thin lithologic hydrocarbon reservoirs: A case study from the Tazhong oil field of the Tarim Basin, China. Geophysics, 81(6): B221-B234. [文内引用:1]
[36] Janocko M, Nemeca W, Henriksen S, Warchol M. 2013. The diversity of deep-water sinuous channel belts and slope valley-fill complexes. Marine and Petroleum Geology, 41: 7-34. [文内引用:1]
[37] Lindseth R O. 1979. Synthetic sonic logs: A process for stratigraphic interpretation. Geophysics, 44(1): 3-26. [文内引用:1]
[38] Luo Y N, Huang H D, Yang Y D, Li Q X, Zhang S, Zhang J W. 2018. Deepwater reservoir prediction using broadband seismic-driven impedance inversion and seismic sedimentology in the South China Sea. Interpretation, 6(4), SO17-SO29. [文内引用:1]
[39] McAulay A D. 1985. Prestack inversion with plane-layer point source modeling. Geophysics, 50(1): 77-89. [文内引用:1]
[40] Posamentier H W. 2001. Seismic geomorphology and depositional systems of deep water environments;observations from offshore Nigeria, Gulf of Mexico, and Indonesia(abstract). AAPG Annual Convention Program, 10: 160. [文内引用:3]
[41] Posamentier H W, Venkatarathnan Kolla. 2003. Seismic Geomorphology and Stratigraphy of Depositional Elements in Deep-water Settings. Journal of Sedimentary Research, 73(3): 367-388. [文内引用:1]
[42] Rothman D H. 1998. Geostatistical inversion of 3-D seismic data for thin-sand delineation. Geophysics, 51(2): 332-346. [文内引用:1]
[43] Schlager W. 2000. The future of applied sedimentary geology. Journal of Sedimentary Research, 70: 2-9. [文内引用:2]
[44] Sheriff R E. 2002. Encyclopedic Dictionary of Applied Geophysics(4th Edition). Society of Exploration Geophysicists, 1-429. [文内引用:2]
[45] Stark T. 2006. Visualization techniques for enhancing stratigraphic inferences from 3D seismic data volumes. First Break, 24(4): 75-85. [文内引用:1]
[46] Vail P R, Mitchum Jr R M, Thompson Ⅲ S. 1977a. Stratigraphic interpretation of seismic reflection patterns in depositional sequences. In: Payton C E(ed). Seismic Stratigraphy: AAPG Memoir, 26, 63-82. [文内引用:1]
[47] Wood L J. 2007. Quantitative seismic geomorphology of Pliocene and Miocene fluvial systems in the northern Gulf of Mexico, U. S. A. Journal of Sedimentary Research, 77: 713-730. [文内引用:1]
[48] Yue D L, Li W, Wang W R, Hu G Y, Qiao H L, Hu J J, Zhang M L, Wang W F. 2019. Fused spectral-decomposition seismic attributes and forward seismic modelling to predict sand bodies in meand ering fluvial reservoirs. Marine and Petroleum Geology, 99: 27-44. [文内引用:1]
[49] Zeng H L, Backus M M, Barrow K T, Tyler N. 1998a. Stratal slicing: Part Ⅰ. Realistic 3-D seismic model. Geophysics, 63(2): 502-513. [文内引用:6]
[50] Zeng H L, Henry S C, Riola J P. 1998b. Stratal slicing: Part Ⅱ. Real 3-D seismic data. Geophysics, 63(2): 514-522. [文内引用:5]
[51] Zeng H L, Kerans C. 2003. Seismic frequency control on carbonate seismic stratigraphy: A case study of the Kingdom Abo sequence, west Texas. AAPG Bulletin, 87(2): 273-293. [文内引用:1]
[52] Zeng H L, Hentz T F. 2004. High-frequency sequence stratigraphy from seismic sedimentology: Applied to Miocene, Vermilion Block 50, Tiger Shoal area, offshore Louisiana. AAPG Bulletin, 88(2): 153-174. [文内引用:2]
[53] Zeng H L, Backus M M. 2005a. Interpretive advantages of 90°-phase wavelets: Part 1-Modeling. Geophysics, 70(3): C7-C15. [文内引用:2]
[54] Zeng H L, Backus M M. 2005b. Interpretive advantages of 90°-phase wavelets: Part 2-Seismic applications. Geophysics, 70(3): C17-C24. [文内引用:2]
[55] Zeng H L. 2018. What is seismic sedimentology?A tutorial. Interpretation, 6(2): SD 1-12. [文内引用:3]
[56] Zhao W Z, Zou C N, Chi Y L, Zeng H L. 2011. Sequence stratigraphy, seismic sedimentology, and lithostratigraphic plays: Upper Cretaceous, Sifangtuozi area, southwest Songliao Basin, China. AAPG Bulletin, 95(2): 241-265. [文内引用:1]
[57] Zhu X M, Zeng H L, Li S L, Dong Y L, Zhu S F, Zhao D N, Huang W. 2017. Sedimentary characteristics and seismic geomorphologic responses of shallow-water delta of Qingshankou Formation in Songliao Basin, China. Marine and Petroleum Geology, 79: 131-148. [文内引用:1]
[58] Zhu X M, Pan R, Li S L, Wang H B, Zhang X, Ge J W, Lu Z Y. 2018. Seismic sedimentology of sand -gravel bodies on steep slope of rift basins-a case study of Shahejie Formation, Dongying Sag, Eastern China. Interpretation, 6(2): SD13-SD27. [文内引用:1]
[59] Herglotz G. 1907. Über das Benndorfsche Problem der Fortpflanzungsgeschwindigkeit der Erdbebenstrahlen. Physikal, Z. Geophys, 8: 145-147. [文内引用:1]