高分辨率层序基准面变化(半)定量分析方法及其在松辽盆地泉四段的应用

林孝先, 侯中健

(油气藏地质及开发工程国家重点实验室(成都理工大学)，成都 610059)

从海相层序地层到陆相层序地层，再到陆相高分辨率层序地层，从最初的剖面观察、对比等的定性研究，到现在结合地球化学、地球物理等手段的(半)定量研究及探索，层序地层学经历了60多年的发展，不管是从理论、方法、技术手段，还是实践研究成果上，都取得了长足的进步。在中国，由于陆相地层的普遍性、复杂性，由T.A.Cross等(1994，1996)提出的高分辨率层序地层学理论[1,2]更受青睐，应用效果也较好,如其在油气田勘探、开发的各个阶段，基准面、湖(海)平面变化分析，沉积盆地分析，甚至(煤、铀等)沉积层控矿床等方面都取得了良好的研究成果。

基准面及其变化是高分辨率层序地层学研究的核心问题。基准面变化直观反映在可容纳空间的变化上，其控制因素主要有构造运动、沉积物供给速率、相对湖(海)平面变化、气候、湖-海沟通事件、河流平衡剖面迁移等[3],以前4个为主。在气候和沉积物供给速率相对稳定的前提下，可容纳空间的变化受构造运动和湖(海)平面变化控制，其叠加等效于相对湖(海)平面变化[3-5]。所以，基准面变化研究常通过可容纳空间变化、相对湖(海)平面变化研究得以实现，其中以可容纳空间变化的Fischer图解(又称可容纳空间图解)分析应用最多[3,4,6-20]。但是，基准面变化研究仍面临几个亟待解决的问题：①目前的基准面变化分析方法，如Fischer图解分析，多只能定性反映基准面变化，(半)定量效果较差；②长周期、低频(中期及以上级别)层序的基准面变化研究较多，而对短周期、高频(短期、超短期级别)层序的基准面变化研究较少、效果较差；③基准面变化研究中常用的Fischer图解始于对灰岩剖面的研究[8]，如何正确、有效地推广、应用于碎屑岩剖面研究？

基于这些问题，作者通过在鄂尔多斯盆地、渤海湾盆地、松辽盆地进行高分辨率层序地层学研究的经验和成果，总结、提出了一套高分辨率层序基准面变化(半)定量分析方法，并以其在松辽盆地红岗北地区下白垩统泉头组第四段(扶余油层)的应用为实例，对该方法及其技术流程进行了阐述和说明，以期达到与同行交流的目的。

1 高分辨率层序基准面变化(半)定量分析方法

上述基准面变化的前2个问题，归根结底是高分辨率层序的精细、定量划分问题。本文的高分辨率层序地层分析以T.A.Cross等(1994,1996)提出的高分辨率层序地层学理论[1,2]为指导。在层序级别划分上，以郑荣才等(2001，2010)提出的六级层序级别与时限[3,21]为标准。高分辨率层序地层的(半)定量研究还处于发展阶段，目前使用的方法主要有测井小波分析法、测井(深度域)频谱分析法、测井曲线活度法、测井最大熵谱分析法、砂泥比曲线法等[22-24]，其中以被誉为“数学显微镜”的测井小波分析法运用最多、最方便、效果最好。至于第3个问题，Fischer图解由Fischer(1964)提出[8]后，经Goldhammer等(1987)、Read等(1988)、Sadler等(1993)的研究、修正和完善[9-11]，目前在国内Fischer图解的应用领域除传统的碳酸盐岩剖面[12,16,25,26]外，已经初步推广到碳酸盐岩、碎屑岩混合沉积剖面[13,20,27,28]和碎屑岩剖面[15,16,19,29-32]。但必须指出的是，碎屑岩剖面应用Fischer图解的前提是相对湖(海)平面变化控制基准面、可容纳空间的变化，且选用的剖面层序必须是异旋回层序，而非自旋回层序，在构造相对稳定、无地层缺失的情况下，异旋回最低、最精细的级别可以达到超短期(Ⅵ级)旋回[33,34]。

在这种背景下，总结、提出了高分辨率层序基准面变化(半)定量分析方法及技术流程。

1.1 分析方法及技术流程

高分辨率层序基准面变化(半)定量分析方法，主要包含高分辨率层序地层精细、(半)定量划分，压实校正，基准面变化(Fischer图解)分析，综合对比分析等内容(图1)。

1.2 高分辨率层序地层划分

高分辨率层序地层划分方法和技术包括以定性方法为主的经典地质学对各级层序界面的识别，以及(半)定量的测井(曲线)小波分析。

1.2.1 经典地质学分析

高分辨率层序研究的经典地质学方法是对层序界面以及层序内部最大湖(洪)泛面的识别(表1)，识别标志主要体现在地震剖面、测井曲线、岩性(地表剖面+岩心)、构造作用面、地球化学特征等方面，以定性研究为主。

1.2.2 测井小波分析

高分辨率层序地层的测井小波分析的基本原理是：测井数据是目前连续性、分辨率最好的地质数据，其蕴含丰富的地质信息；但这些信息不分时间、级别地叠加在一起，造成从测井数据或曲线上很难识别出各种隐含的地质周期信息；而小波分析可以通过数学计算和变换，把这些信息分解成不同尺度、频带的信号，得到相应的地质周期信息。小波分析中的低频信号代表地质长周期，可用来识别大的地层旋回；相反，高频信号代表短周期，可用于小层的精细划分和对比；而频带信号中的奇异点(剧变点)、剧变程度可用于划分高分辨率层序界面和最大湖(洪)泛面。

图1 高分辨率层序基准面变化分析方法及技术流程图Fig.1 The method and technology flow chart of analysing the base level changes of the high-resolution sequence

表1 高分辨率层序界面(最大湖(洪)泛面)类型与特征Table 1 Types and characteristics of the high-resolution sequence boundary and maximum flooding surfaces

高分辨率层序地层测井小波分析可以通过Matlab软件中的小波分析工具箱实现[6,7,23,36-39],具体操作步骤如下。

a.选择数据，常用于小波分析的测井数据有自然伽马、自然电位、声波时差、电阻率、感应电导率等，在砂、泥岩互层的地区，以自然伽马使用的效果最好；测井数据在使用之前需要校正、标准化预处理。

b.准备数据，选择一种测井数据，确定研究层位的顶、底深度，把其范围内的数据单独作为一列存为txt格式文件。

c.载入、保存信号数据，运行Matlab软件，在“Command Window”中输入以下命令：

>>文件名=load(‘txt文件所在目录的完整路径’)

(回车后即完成txt文件的载入，检查正确与否后，可进行下面的保存操作)

>>save 文件名

(回车后即可在“Current Folder”中见保存的“文件名.mat”文件)

d.打开小波分析工具箱。有2种方法，一是在“Command Window”中输入“wavemenu”命令，二是执行“APPS”→下拉键→“Wavelet Design & Analysis”操作。

e.在弹出的“Wavelet Toolbox Main Menu”对话框中，选择“Wavelet 1-D”(一维离散小波)或“Continuous Wavelet 1-D”(一维连续小波)进行分析。短期(Ⅴ级)及以上级别的层序，两种分析都可选用，而超短期(Ⅵ级)层序选用一维离散小波效果较好。

f.在“Wavelet 1-D”(一维离散小波)窗口中，执行“File”→“Load”→“Signal”操作，选择刚保存的“文件名.mat”文件，载入mat文件成功后，再选择“Wavelet”的类型和“Level”，确定后点击“Analyze”会出现相应级别的分解结果；然后点击“De-noise”给数据去噪，在“Wavelet 1-D-De-noising”窗口中，设置默认，点击“De-noise”按钮；去噪完成(可保存去噪数据备用)后点击“Close”按钮，弹出对话框中选择“Yes”更新分析数据即可完成去噪信号的小波分析。

g.在“Continuous Wavelet 1-D”(一维连续小波)窗口中，执行“File”→“Load Signal”操作，选择、载入(去噪后的)mat文件成功后，再选择“Wavelet”的类型和“Sampling Period”，设置“Scale Settings”中的“Max.”值(其为小波系数曲线尺度的2倍)，其余默认设置，点击“Analyze”即可。

h.输出结果，执行“File”→“Export Setup...”操作，在弹出对话框中，各项参数可以默认(实用的是，“Rendering”目录下的“Resolution(dpi)”可以设置图片分辨率或大小)；点击“Export...”按钮，设置输出文件名、类型及位置等。

i.选取合适级别和尺度的测井小波曲线，对高分辨率层序进行划分。

小波分析工具箱中具有多种小波类型和级别，可试验选择最佳的设置、提取方案，或同时进行多种设置、提取方案对比研究。对测井小波分析所得的各级别、周期的小波曲线，可据其级别、周期特征进行相应级别、周期的高分辨率层序地层划分。

1.3 基准面变化分析

基准面变化分析主要包括压实校正、Fischer图解分析、综合对比分析等。

1.3.1 压实校正分析

压实校正分析以高分辨率层序地层的精细、定量划分为基础，是Fischer图解分析的前提。压实校正的理论基础主要有地层骨架体积不变(孔隙度-深度曲线)和地层骨架质量不变(密度-深度曲线)2种[40]，以前者应用最多[41-45]。

以地层骨架体积(厚度)不变为基础，对现今埋藏的地层或岩石进行原始沉积厚度还原，需要使用以下压实校正公式[40,42-44]

其中：Φ(z)表示孔隙度(Φ)-深度(z)的函数；d表示地层或特定岩性沉积时的原始厚度或其底面深度(m)；z1、z2分别表示地层或特定岩性的埋藏顶面、底面的深度(m)。

孔隙度-深度(曲线)函数Φ(z)是由孔隙度、深度数据拟合而成，有十多种经验公式或模型[46]，其中最常用的公式如下。

一元线性回归方程：Φ(z)=Φ0-cz(2)

指数回归方程：Φ(z)=Φ0e-cz(3)

其中Φ0表示沉积时的原始孔隙度，c表示压实(压缩)系数。

在具有不同Φ0和c的砂岩、泥岩互层的碎屑岩地区进行压实校正，对孔隙度-深度函数Φ(z)有2种处理方法。一是按照砂岩、泥岩的体积分数及其Φ砂(z) 、Φ泥(z)求解Φ(z)[16]；二是把砂岩、泥岩看做一个整体求解Φ(z)[42]。两种方法各有优缺点，可据实际情况选择，带入公式(1)，求解方程可得原始厚度(d)。如把公式(2)带入公式(1)，求得原始厚度为

(4)

1.3.2 Fischer图解分析

Fischer图解的横坐标常用的是层序旋回个数(最低级别旋回个数最好是>50个[11])，在可精确确定年龄且无地层缺失的情况下，也可使用具体地质年龄；纵坐标为经压实校正过的旋回厚度与其平均旋回厚度差值的累积，即平均厚度累计偏移，其起止点纵坐标值为零或相等(图2)。编制Fischer图解的步骤如下。

图2 Fischer图解示意图Fig.2 Schematic diagram of Fischer plots(据Sadler等(1993)[11]、郑荣才等(2010)[3]修编)

a.按照上述高分辨率层序划分及其压实校正结果，计算旋回的平均厚度和每个旋回的净增加或减少厚度(即旋回厚度与平均旋回厚度的差值)。

b.以0作为起点，1旋回的净增加或减少厚度为其纵坐标值，投点得A点。连接起点和A点，该连线即为旋回1的基准面(或可容纳空间)变化曲线。净增加厚度表示基准面上升，净减少厚度表示基准面下降。

c.以A点作为旋回2基准面变化的起点，旋回2的纵坐标为旋回1和旋回2的净增加或减少厚度之和，即二者的累计值，投点得B点，连接A、B即得旋回2的基准面变化曲线。

d.以B点作为旋回3基准面变化的起点，重复上面的工作，最后所得连线即是完整的基准面变化曲线。

1.3.3 综合对比分析

综合对比分析即用Fischer图解分析所得基准面变化曲线与高分辨率层序划分(最大湖(洪)泛面)，以及相对湖(海)平面、绝对海平面变化等前人研究成果进行对比，分析其相似性与差异性，完善基准面变化规律。

2 应用实例

以上述高分辨率层序基准面变化(半)定量分析方法为指导，对松辽盆地红岗北地区(图3)下白垩统泉头组第四段(简称“泉四段”)基准面变化规律进行了研究。泉四段沉积时期，研究区为近海内陆湖盆环境，正处于拗陷盆地发育早期，持续沉降，气候温暖、潮湿，沉积体系受西部的白城水系控制，沉积相以发育三角洲平原、前缘亚相为主，岩性以细砂岩、粉砂岩、泥岩为主[47-51]。

2.1 高分辨率层序地层划分

对研究区泉四段进行了经典地质学的层序界面和最大湖(洪)泛面识别。Ⅰ、Ⅱ级层序界面因泉四段与上覆青一段、下伏泉三段为整合接触故未发育；Ⅲ、Ⅳ、Ⅴ级层序界面和相应级别的最大湖(洪)泛面均有发育，且易识别；Ⅵ级层序界面识别较为困难，受干扰较大，只在部分钻井有所发现(图4)，难以用于层序划分。故依此划分出1个长期(Ⅲ级)层序、4个中期(Ⅳ级)层序、22个短期(Ⅴ级)层序(图5)。

图3 研究区位置及构造区划图Fig.3 The map of location and tectonic division of the northern Honggang area of Songliao Basin (据叶得泉等(2002)[50]修编)

在测井小波分析中，选择研究区代表性的Y1井的自然伽马(GR)曲线进行了分析。为方便分析结果的联井对比，对GR曲线进行了极差正规法[GR标准化=(GR-GRmin)/(GRmax-GRmin)][52]的简单校正和标准化。为使分析结果更可靠，同时选择一维离散小波和一维连续小波进行分析。在小波类型选择上，参考前人资料[6,7,23,37-39,53-59]，根据实际试验、分析，最后选择效果最好的db10(level9)和dmey9小波进行对比分析。据db10、dmey9一维离散小波分析中的a9曲线和一维连续小波分析中的a=128尺度的小波系数曲线，划分出1个长期(Ⅲ级)层序的最大湖(洪)泛面及相应的2个层序界面；据d7曲线和a=64曲线划分出4个中期(Ⅳ级)层序；据d5曲线和a=32曲线划分出22个短期(Ⅴ级)层序；据d3曲线(参考d2曲线)划分出69个超短期(Ⅵ级)层序。小波分析对长期、中期、短期层序的划分与经典地质学方法的分析结果很吻合，且精细到超短期层序(图5)。

划分出的长期、中期、短期层序在研究区内能相互对比，长期层序与泉四段对应，中期层序对应于4个亚段，短期层序对应于22个层；超短期层序的井间可对比性不是十分良好，有一定的差异性，3个左右超短期层序对应于1个小层，但在短期层序(或层)格架下的总体特征不变。另外，从高分辨率层序及其时限上，泉四段年龄为100～104 Ma[49,60,61]，那么划分出的长期(Ⅲ级)层序时限为4 Ma、中期(Ⅳ级)层序时限为1 Ma、短期(Ⅴ级)层序时限约为0.182 Ma、超短期(Ⅵ级)层序时限约为0.058 Ma，这与郑荣才等(2001，2010)的六级层序级别与时限[3,21]基本一致，可相互对比。

图4 研究区泉四段高分辨率层序界面Fig.4 High-resolution sequence boundary of Member 4 of Quantou Formation in the northern Honggang area of Songliao Basin(A)Ⅲ级界面,X1井; (B)Ⅳ级界面,X2井; (C)Ⅴ级界面,X3井; (D)Ⅵ级界面, X4井; (E)最大湖(洪)泛面(mfs),X5井

2.2 基准面变化分析

首先，对研究区嫩江组—泉头组139组684个(不同粒度)砂岩、泥岩孔隙度数据进行分析发现，一元线性回归方程Φ=42.288 56-0.016 755z对数据的拟合程度较高，相关系数R2=0.839 23(图6)。由于孔隙度数据来源于砂岩和泥岩，故所得回归方程反映的是地层整体的孔隙度-深度关系，避免了在确定砂岩和泥岩Φ0、c过程中产生的误差，且与姚秀云等(1989)、修洪文(2008)的分析结果[62,63]相似，可以对比。把研究区的Φ0=0.422 885 6和c=0.000 167 55(化为常数，不带单位)带入公式(4)，以Y1井泉四段69个超短期(Ⅵ级)层序为压实校正单元，可求解其沉积时的原始厚度。现在地层总厚度109.60 m，校正后的原始总厚度为172.82 m。

然后，对研究区Y1井泉四段经压实校正后的69个超短期层序进行了Fischer图解的绘制，采用旋回个数(图7)和柱状图实际旋回厚度(图8)两种表达方式。由此可以看出，地层亚段与基准面变化有很好的对应关系，Ⅳ、Ⅲ亚段分别对应于一次完整的基准面上升—下降旋回，Ⅱ亚段基准面上下摆动、相对稳定，Ⅰ亚段基准面持续下降，其中Ⅲ亚段的基准面上升和Ⅰ亚段的基准面下降规模较大、持续时间较长，其余规模相对较小、持续时间较短(图7、图8)。

图5 Y1井泉四段小波分析及高分辨率层序划分方案Fig.5 Wavelet analysis and classification scheme of the high-resolution sequence (Well Y1) of Member 4 of Quantou Formation in the northern Honggang area of Songliao Basin图中的a和d表示小波分析的分解、伸缩尺度，例如a9、d3、a=32尺度

图6 研究区嫩江组—泉头组孔隙度-深度关系图Fig.6 Porosity-depth relation diagram of Nenjiang-Quantou Formation in the northern Honggang area of Songliao Basin

2.3 综合对比分析

把Fischer图解按照旋回厚度投影到综合柱状图上，对其进行综合对比研究，可以检查、提高高分辨率层序地层划分和基准面变化分析的准确性。

研究区Y1井泉四段Fischer曲线与前文高分辨率层序地层划分，以及程日辉等(2008)[19]、魏魁生等(1996)[49]、Haq等(1987)[64]的研究具有很好的相似性(图8)。特别是在3-1小层，在Fischer图解上为较大型的基准面上升—下降转换期，这与前文高分辨率层序地层划分中的长期(Ⅲ级)层序的最大湖(洪)泛面一致；在松科1井中，也为Fischer图解上的一基准面上升—下降转换期[19]；在整个松辽盆地，也具相同的基准面变化规律[49]，甚至与该时期全球性的海侵[64]也有很好的吻合性。这主要与松辽盆地属于近海构造外流盆地(出口处大约在现今宾县地堑附近[47])密切相关。

3 结论与讨论

a.高分辨率层序基准面变化分析方法的核心是高分辨率层序地层学，前提和基础是高分辨率层序的精细、(半)定量划分，目的是以高分辨率层序为单元，通过压实校正计算，绘制Fischer曲线，总结基准面变化规律。

b.松辽盆地红岗北地区泉四段基准面变化与地层亚段有很好的对应关系，具一期基准面长时间、大规模上升—短时间、小规模下降规律，与研究区湖平面、松辽盆地基准面、全球海平面变化规律一致。

c.定性的经典地质学分析结合(半)定量的测井小波分析能实现高分辨率层序的精细划分，是高分辨率层序地层学未来发展的一个方向。

d.测井数据能提供大量的、连续的、定量的地质信息，如何把它们更好、更方便、更有效地提取出来，是未来值得重视和研究的方向。

图7 研究区Y1井泉四段Fischer图解分析Fig.7 Fischer plots of Well Y1 Member 4 of Quantou Formation in the northern Honggang area of Songliao Basin

图8 Y1井泉四段Fischer图解综合柱状图分析与对比Fig.8 Synthesis columnar section analysis and comparison of Fischer plots of Well Y1 of Member 4 of Quantou Formation in the northern Honggang area of Songliao Basin

[1] Cross T A. High-resolution stratigraphic correlation from the Perception of base-level cycles and sediment accommodation [C]//Proceeding of Northwestern European Sequence Stratigraphy Congress, 1994: 105-123.

[2] Cross T A, Lessenger M A. Sediment volume partitioning: rationale for stratigraphic model evaluation and high-resolution stratigraphic correlation [C]//Norwegian Petroleums-Forening Conference Volume, 1996: 1-24.

[3] 郑荣才,文华国,李凤杰.高分辨率层序地层学[M].北京:地质出版社,2010.

Zheng R C, Wen H G, Li F J. High-Resolution Sequence Stratigraphy [M]. Beijing: Geological Publishing House, 2010. (In Chinese)

[4] 胡受权,陈国能,王英民,等.Fischer图解及其沉积响应的计算机模拟——以泌阳断陷下第三系核三上段为例[J].石油与天然气地质,1999,20(1):70-75.

Hu S Q, Chen G N, Wang Y M,etal. Fischer diagram and its computer simulation of sedimentary response: A case study on the Upper Member of Eh3in Biyang fault-depression [J]. Oil & Gas Geology, 1999, 20(1): 70-75. (In Chinese)

[5] 王亚青,董春梅,邢焕清.可容空间的研究进展[J].海洋地质动态,2004,20(10):32-35.

Wang Y Q, Dong C M, Xing H Q. Development of the accommodation concept in sequence stratigraphic researches [J]. Marine Geology Letters, 2004, 20(10): 32-35. (In Chinese)

[6] 赵伟,邱隆伟,姜在兴,等.小波分析在高精度层序单元划分中的应用[J].中国石油大学学报:自然科学版,2009,33(2):18-22.

Zhao W, Qiu L W, Jiang Z X,etal. Application of wavelet analysis in high-resolution sequence unit division [J]. Journal of China University of Petroleum, 2009, 33(2): 18-22. (In Chinese)

[7] 赵伟,姜在兴,邱隆伟,等.小波分析划分层序单元的地质学理论基础、方法与应用[J].石油与天然气地质,2010,31(4):436-441.

Zhao W, Jiang Z X, Qiu L W,etal. Geological concept, method and application of sequence unit identification through wavelet analysis [J]. Oil & Gas Geology, 2010, 31(4): 436-441. (In Chinese)

[8] Fischer A G. The Lofer cyclothems of the Alpine Triassic[J]. Kansas Geological Survey Bulletin, 1964, 169(1): 107-149.

[9] Goldhammer R K, Dunn P A, Hardie L A. Highfrequency glacio-eustatic sea-level oscillations with Milankovitch characteristics recorded in Middle Triassic platform carbonates in Northern Italy[J]. American Journal of Science, 1987, 287(9): 853-892.

[10] Read J F, Goldhammer R K. Use of Fischer plots to define third-order sea-level curves in Ordovician peritidal cyclic carbonates, Appalachians[J]. Geology, 1988, 16: 895-899.

[11] Sadler P M, Osleger D A, Montanez L P. On the labeling, length, and objective basis of Fischer plots[J]. Journal of Sedimentary Petrology, 1993, 63(2): 360-368.

[12] 苏德辰,李庆谋,罗光文,等.Fischer图解及其在旋回层序研究中的应用——以北京西山张夏组为例[J].现代地质,1995,9(3):279-283.

Su D C, Li Q M, Luo G W,etal. Method for drawing Fischer plots and its applications on studying of cyclic sequences: Example from Middle Cambrian Zhangxia Formation, Western Hills of Beijing [J]. Geoscience, 1995, 9(3): 279-283. (In Chinese)

[13] 赵玉光,许效松,刘宝珺.上扬子台地西缘峨眉地区三叠纪高频层序与海平面振荡研究[J].岩相古地理,1996,16(1):1-18.

Zhao Y G, Xu X S, Liu B J. High-frequency sequences and sea-level oscillation in the Emei area on the western margin of the Upper Yangtze Platform [J]. Lithofacies Paleogeography, 1996, 16(1): 1-18. (In Chinese)

[14] Day P I. The Fischer diagram in the depth domain: a tool for sequence stratigraphy[J]. Journal of Sedimentary Research, 1997, 67(5): 982-984.

[15] 胡受权.泌阳断陷双河—赵凹地区下第三系核三上段陆相层序发育的可容空间机制[J].沉积学报,1998,16(3):102-108.

Hu S Q. Accommodation mechanics of terrigenous sequence development in fault-depressed lacustrine basin——a case study: the Upper Member of Eh3, Shuanghe-Zhaoao, Biyang Fault-depression [J]. Acta Sedimentologica Sinica, 1998, 16(3): 102-108. (In Chinese)

[16] 胡受权.泌阳断陷下第三系核三上段陆相层序可容空间变化机理及其图解[J].特种油气藏,1998,5(4):1-6.

Hu S Q. Mechanism of allowable space change and its diagram in continental sequence, upper member of He-3, Paleogene, Biyang depression [J]. Special Oil and Gas Reservoirs, 1998, 5(4): 1-6. (In Chinese)

[17] 翟永红.用Fischer图解研究山西临汾中奥陶世马家沟组旋回层序[J].岩石矿物学杂志,1999,18(2):128-133.

Zhai Y H. The application of Fischer plots to the analysis of cycles in Middle Ordovician Majiagou Formation in Linfen, Shanxi Province [J]. Acta Petrologica Et Mineralogica, 1999, 18(2): 128-133. (In Chinese)

[18] 李江涛,李增学,郭建斌,等.高分辨率层序地层分析中基准面变化的讨论[J].沉积学报,2005,23(2):297-302.

Li J T, Li Z X, Guo J B,etal. Discussion about base-level changes in the analysis of high-resolution sequence stratigraphy [J]. Acta Sedimentologica Sinica, 2005, 23(2): 297-302. (In Chinese)

[19] 程日辉,王国栋,王璞珺.松辽盆地白垩系泉三段—嫩二段沉积旋回与米兰科维奇周期[J].地质学报,2008,82(1):55-64.

Cheng R H, Wang G D, Wang P J. Sedimentary cycles of the Cretaceous Quantou-Nenjiang Formations and Milankovitch cycles of the south hole of the SLCORE-I in the Songliao Basin [J]. Acta Geologica Sinica, 2008, 82(1): 55-64. (In Chinese)

[20] 李凤杰,刘琪,刘殿鹤,等.四川盆地东北部二叠系层序发育的动力学分析[J].地层学杂志,2010,34(1):35-42.

Li F J, Liu Q, Liu D H,etal. Sequence dynamics of the Permian System in Northeastern Sichuan Basin [J]. Journal of Stratigraphy, 2010, 34(1): 35-42. (In Chinese)

[21] 郑荣才,彭军,吴朝容.陆相盆地基准面旋回的级次划分和研究意义[J].沉积学报,2001,19(2):249-255.

Zheng R C, Peng J, Wu C R. Grade division of base-level cycles of terrigenous basin and its implications [J]. Acta Sedimentologica Sinica, 2001, 19(2): 249-255. (In Chinese)

[22] 李霞,范宜仁,邓少贵,等.利用测井资料研究层序的方法述评[J].测井技术,2006,30(5):411-415.

Li X, Fan Y R, Deng S G,etal. Review on methods for studying strata sequence with log data [J]. Well Logging Technology, 2006, 30(5): 411-415. (In Chinese)

[23] 程道解.苏北盆地测井层序与沉积相研究[D].东营:中国石油大学档案馆,2008.

Cheng D X. The Study of Strata Sequence and Sedimentary Fades Based on Well Logging Data in Su-Bei Basin [D]. Dongying: The Archive of China University of Petroleum, 2008. (In Chinese)

[24] 朱剑兵.基于测井资料的高分辨率层序旋回划分方法探讨[J].油气地球物理,2011,9(4):6-11.

Zhu J B. The discussion on the division methods of high resolution sequence cycle based on logging data [J]. Petroleum Geophysics, 2011, 9(4): 6-11. (In Chinese)

[25] 于炳松.从露头剖面层序中求解海平面变化[J].岩相古地理,1995,15(3):47-52.

Yu B S. Reconstruction of sea-level curves based on outcrop sequences [J]. Lithofacies Paleogeography, 1995, 15(3): 47-52. (In Chinese)

[26] 于炳松.塔里木盆地北部寒武纪海平面变化研究[J].沉积学报,1996,14(1):33-39.

Yu B S. Study on sea-level change of Cambrian in Northern Tarim Basin [J]. Acta Sedimentologica Sinica, 1996, 14(1): 33-39. (In Chinese)

[27] 刘惠民.牛庄洼陷沙三段上亚段—沙二段高精度层序地层学[J].油气地质与采收率,2008,15(5):49-52.

Liu H M. Analysis of high-resolution sequence stratigraphy of the upper Es3to Es2of Niuzhuang Sag [J]. Petroleum Geology and Recovery Efficiency, 2008, 15(5): 49-52. (In Chinese)

[28] 李凤杰,刘殿鹤,郑荣才,等.四川盆地东北地区上二叠统层序地层特征研究[J].沉积学报,2009,27(6):1116-1123.

Li F J, Liu D H, Zheng R C,etal. The sequence stratigraphics characteristics of Upper Permian in the Northeastern Sichuan Basin [J]. Acta Sedimentologica Sinica, 2009, 27(6): 1116-1123. (In Chinese)

[29] 李凤杰,王多云,张庆龙,等.鄂尔多斯盆地陇东地区延长组沉积相特征与层序地层分析[J].沉积学报,2006,24(4):549-554.

Li F J, Wang D Y, Zhang Q L,etal. Sedimentary facies characteristic and sequence stratigraphy analysis of Yanchang Formation in Longdong area, Ordos Basin [J]. Acta Sedimentologica Sinica, 2006, 24(4): 549-554. (In Chinese)

[30] 郭少斌,陈成龙.利用米兰科维奇旋回划分柴达木盆地第四系层序地层[J].地质科技情报,2007,26(4):27-30.

Guo S B, Chen C L. Division of sequence stratigraphy of Quaternary formation in Qaidam Basin using the Milankvitch cycle [J]. Geological Science and Technology Information, 2007, 26(4): 27-30. (In Chinese)

[31] 郭玉新,隋风贵,林会喜,等.时频分析技术划分砂砾岩沉积期次方法探讨——以渤南洼陷北部陡坡带沙四段—沙三段为例[J].油气地质与采收率,2009,16(5):8-11.

Guo Y X, Sui F G, Lin H X,etal. Discussion on the division and correlation of glutinite sedimentary period by time-frequecgy analysis: A case study of Sha4 to Sha3 member of northern abrupt slope zone in Bonan Depression [J]. Petroleum Geology and Recovery Efficiency, 2009, 16(5): 8-11. (In Chinese)

[32] 高有峰,王璞珺,瞿雪姣,等.松辽盆地东南隆起区白垩系嫩江组一段沉积相、旋回及其与松科1井的对比[J].岩石学报,2010,26(1):99-108.

Gao Y F, Wang P J, Qu X J,etal. Sedimentary facies and cyclostratigraphy of the Cretaceous first member of Nenjiang Formation in the southeast uplift zone, Songliao Basin and its correlation with the CCSD-SK-I [J]. Acta Petrologica Sinica, 2010, 26(1): 99-108. (In Chinese)

[33] 纪友亮,周勇,吴胜和,等.河流相地层高精度地层构型界面形成机制及识别方法[J].中国石油大学学报:自然科学版,2012,36(2):8-15.

Ji Y L, Zhou Y, Wu S H,etal. Formation mechanism and recognizing method of high resolution strata architecture boundary in fluvial strata [J]. Journal of China University of Petroleum, 2012, 36(2): 8-15. (In Chinese)

[34] 史长林,纪友亮,刘灵童,等.自旋回—异旋回控制的扇三角洲相高精度层序地层对比模式[J].石油天然气学报,2013,35(3):6-11.

Shi C L, Ji Y L, Liu L T,etal. Research on high-resolution sequence stratigraphic correlation mode under the control of autocyclical-allocyclical fan delta [J]. Journal of Oil and Gas Technology, 2013, 35(3): 6-11. (In Chinese)

[35] 田景春,陈洪德,覃建雄,等.层序底界面的物质表现形式——以中国南方海相震旦系—中三叠统为例[J].沉积学报,2003,21(4):675-682.

Tian J C, Chen H D, Qin J X,etal. The apparent types of sequence bottom boundary surface: Taking marine Sinian System to Middle Triassic Series in South China as an example [J]. Acta Sedimentologica Sinica, 2003, 21(4): 675-682. (In Chinese)

[36] 郑治真,沈萍,杨选辉,等.小波变换及其MATLAB工具的应用[M].北京:地震出版社,2001.

Zheng Z Z, Shen P, Yang X H,etal. Wavelet Transform and Its MATLAB Tools Aapplication [M]. Beijing: Earthquake Press, 2001. (In Chinese)

[37] 余继峰,李增学.测井数据小波变换及其地质意义[J].中国矿业大学学报,2003,32(3):336-339.

Yu J F, Li Z X. Wavelet transform of logging data and its geological significance [J]. Journal of China University of Mining & Technology, 2003, 32(3): 336-339. (In Chinese)

[38] 王志坤,王多云,宋广寿,等.测井信号小波分析在高分辨率层序地层划分中的应用[J].大庆石油学院学报,2005,29(6):17-20.

Wang Z K, Wang D Y, Song G S,etal. Application of wavelet analysis to sequence stratigraphic division of high resolution sequence stratigraphy [J]. Journal of Daqing Petroleum Institute, 2005, 29(6): 17-20. (In Chinese)

[39] 寻知锋.小波变换在储层沉积学中的应用[D].青岛:山东科技大学档案馆,2009.

Xun Z F. Application of Wavelet Transformation in Reservoir Sedimentology [D]. Qingdao: The Archive of Shandong University of Science and Technology, 2009. (In Chinese)

[40] 李绍虎,吴冲龙,吴景富,等.一种新的压实校正法[J].石油实验地质,2000,22(2):110-114.

Li S H, Wu C L, Wu J F,etal. A new method for compaction correction [J]. Petroleum Geology & Experiment, 2000, 22(2): 110-114. (In Chinese)

[41] 武凤良.应用测井、地震资料进行地层压实校正的方法[J].石油地球物理勘探,1989,24(1):68-78.

Wu F L. The method for making stratigraphic compaction correction with the use of logging and seismic data [J]. Oil Geophysical Prospecting, 1989, 24(1): 68-78. (In Chinese)

[42] 关振良.利用压实模拟计算地层古厚度及差异压实量的尝试[J].石油实验地质,1992,14(2):152-158.

Guan Z L. An attempt to calculate the primary thickness and the amount of differential compaction in a formation with an application of compaction simulations [J]. Petroleum Geology & Experiment, 1992, 14(2): 152-158. (In Chinese)

[43] 秦玉娟,贾振远,蔡忠贤.海平面变化定量研究[J].地球科学:中国地质大学学报,1997,22(5):460-465.

Qin Y J, Jia Z Y, Cai Z X. Quantitative research on sea level change [J]. Earth Science, 1997, 22(5): 460-465. (In Chinese)

[44] 杨桥,漆家福.碎屑岩层的分层去压实校正方法[J].石油实验地质,2003,25(2):206-210.

Yang Q, Qi J F. Method of delaminated decompaction correction [J]. Petroleum Geology & Experiment, 2003, 25(2): 206-210. (In Chinese)

[45] Allen P A, Allen J R. Basin analysis: principles and applications [M]. Oxford: Blackwell Publishing, 2005.

[46] 马立祥,万静萍.前缘盆地砂岩孔隙度预测的一种新途经——孔隙度与深度关系多数学模型选择法[J].石油实验地质,1994,16(1):70-77.

Ma L X, Wan J P. A new approach to the porosity prediction in the sandstones of frontier oil/gas basins [J]. Petroleum Geology & Experiment, 1994, 16(1): 70-77. (In Chinese)

[47] 杨万里.松辽陆相盆地石油地质[M].北京:石油工业出版社,1985.

Yang W L. Petroleum Geology in Songliao Continental Basin [M]. Beijing: Pertroleum Industry Press, 1985. (In Chinese)

[48] 王东坡,刘招君,刘立.松辽盆地演化与海平面升降[M].北京:地质出版社,1994.

Wang D P, Liu Z J, Liu L. Basin Evolution and Eustasy in Songliao Basin [M]. Beijing: Geological Publishing House, 1994. (In Chinese)

[49] 魏魁生.非海相层序地层学——以松辽盆地为例[M].北京:地质出版社,1996.

Wei K S. Nonmarine Sequence Stratigraphy: A Case Study of the Songliao Basin [M]. Beijing: Geological Publishing House, 1996. (In Chinese)

[50] 叶得泉,黄清华,张莹,等.松辽盆地白垩纪介形类生物地层学[M].北京:石油工业出版社,2002.

Ye D Q, Huang Q H, Zhang Y,etal. Cretaceous Ostracoda Biostratigraphy in Songliao Basin [M]. Beijing: Pertroleum Industry Press, 2002. (In Chinese)

[51] 沈安江,康伟力,王艳清,等.松辽盆地南部白垩纪层序地层与岩性地层油气藏勘探[M].北京:石油工业出版社,2006.

Shen A J, Kang W L, Wang Y Q,etal. Sequence Stratigraphy and Lithologic Stratigraphic Reservoir Exploration in Cretaceous in Southern Songliao Basin [M]. Beijing: Pertroleum Industry Press, 2006. (In Chinese)

[52] 李汉林,赵永军,王海起.石油数学地质(第二版)[M].东营：中国石油大学出版社,2008.

Li H L, Zhao Y J, Wang H Q. Petroleum Mathematic Geology [M]. Dongying: China University of Petroleum Press, 2008. (In Chinese)

[53] 朱剑兵,纪友亮,赵培坤,等.小波变换在层序地层单元自动划分中的应用[J].石油勘探与开发,2005,32(1):84-86.

Zhu J B, Ji Y L, Zhao P K,etal. Application of wavelet transform in auto-identify units of stratigraphy sequence [J]. Petroleum Exploration and Development, 2005, 32(1): 84-86. (In Chinese)

[54] 刘文业.小波分析在高分辨率层序地层划分中的应用[J].地质力学学报,2006,12(1):64-70.

Liu W Y. Application of wavelet analysis in high-resolution sequence stratigraphic division [J]. Journal of Geomechanics, 2006, 12(1): 64-70. (In Chinese)

[55] 郑希民,郭彦如,刘化清,等.应用自然伽玛测井曲线小波分析划分陆相坳陷盆地三级层序的方法:以鄂尔多斯盆地延长组为例[J].天然气地球科学,2006,17(5):672-676.

Zheng X M, Guo Y R, Liu H Q,etal. The method of division third-degree succession of strata by gamma ray logging wavelet analysis in stationary continental facies depression basin: taking Yanchang Formation of Ordos Basin as an example [J]. Natural Gas Geoscience, 2006, 17(5): 672-676. (In Chinese)

[56] 李相博,郭彦如,刘化清,等.浅谈小波分析在鄂尔多斯盆地延长组层序地层划分中的应用[J].天然气地球科学,2006,17(6):779-782.

Li X B, Guo Y R, Liu H Q,etal. The application of wavelet analysis in sequence stratigraphic subdivision of the Yanchang Formation, Ordos Basin [J]. Natural Gas Geoscience, 2006, 17(6): 779-782. (In Chinese)

[57] 赵军龙,谭成仟,李娜,等.小波分析在高分辨率层序地层研究中的应用[J].地球科学与环境学报,2007,29(1):90-94.

Zhao J L, Tan C Q, Li N,etal. Application of wavelet analysis in high resolution sequence analysis [J]. Journal of Earth Sciences and Environment, 2007, 29(1): 90-94. (In Chinese)

[58] 杨勇强,邱隆伟,陈世悦,等.基于小波能量谱系图及小波曲线的层序地层划分[J].石油地球物理勘探,2011,46(5):783-789.

Yang Y Q, Qiu L W, Chen S Y,etal. Sequence stratigraphy identification based on wavelet energy spectrum and wavelet curve [J]. Oil Geophysical Prospecting, 2011, 46(5): 783-789. (In Chinese)

[59] 高迪,郭变青,邵龙义,等.基于MATLAB的小波变换在沉积旋回研究中的应用[J].物探化探计算技术,2012,34(4):444-448.

Gao D, Guo B Q, Shao L Y,etal. Application of wavelet transform in sedimentary cycle research based on Matlab [J]. Computing Techniques for Geophysical and Geochemical Exploration, 2012, 34(4): 444-448. (In Chinese)

[60] 王璞珺,杜小弟,王东坡.松辽盆地白垩纪湖侵沉积层序与湖海沟通事件的地球化学记录[J].岩相古地理,1995,15(4):14-20.

Wang P J, Du X D, Wang D P. Geochemical records of the Cretaceous lacustrine transgressive sequence and the marine-lake connecting events in the Songliao Basin, Northeast China [J]. Lithofacies Paleogeography, 1995, 15(4): 14-20. (In Chinese)

[61] 王璞珺,杜小弟,王俊,等.松辽盆地白垩纪年代地层研究及地层时代划分[J].地质学报,1995,69(4):372-381.

Wang P J, Du X D, Wang J,etal. The chronostratigraphy and stratigraphic classification of the Cretaceous of the Songliao Basin [J]. Acta Geologica Sinica, 1995, 69(4): 372-381. (In Chinese)

[62] 姚秀云,张凤莲,赵鸿儒.岩石物性综合测定——砂、泥岩孔隙度与深度及渗透率关系的定量研究[J].石油地球物理勘探,1989,24(5):533-541.

Yao X Y, Zhang F L, Zhao H R. Comprehensive determination of physical properties of rock formation: The quantitative analysis of the relation between porosity and depth (or permeability) for sandstone (or mudstone) [J]. Oil Geophysical Prospecting, 1989, 24(5): 533-541. (In Chinese)

[63] 修洪文.松辽盆地北部泉三四段成岩研究与储层评价[D].大庆:大庆石油学院档案馆,2008.

Xiu H W. Diagenesis Research and Reservoir Evaluation in Quan3 and Quan4 Members of Northern Songliao Basin [D]. Daqing: The Archive of Daqing Petroleum Institute, 2008. (In Chinese)

[64] Haq B U, Hardenbol J, Vail P R. Chronology of fluctuating sea levels since the Triassic [J]. Science, 1987, 235: 1156-1167.