层序地层分析的研究方法首先是由Exxon生产研究公司的地质学家们(Van Wagoner等,1988;Posamentier等,1988;Van Wagoner等,1990;Sangree等,1991;Scholz等,1991)提出的。自20世纪50年代开始,经历了50年代的沉积相模式,60年代的沉积体系分析和70年代的地震地层学等发展阶段(Vail等,1977;Mitchum等,1977),至80年代提出层序地层学的概念与方法体系,经历了地质学家和勘探工作者近半个世纪的努力。随着高精度地震资料的应用,结合详细的露头资料和钻孔岩芯及测井资料,Van Wag-oner等(1990)和Sangree等(1991)对层序概念及层序构成作了更精辟分析,并提出一套完整的概念体系和系统的研究方法。层序地层学是根据地震、钻孔和露头资料对地层型式作出综合解释。层序地层学成为一种划分、对比和分析沉积岩系的新方法,其研究实质就是建立沉积盆地的等时地层格架(李思田等,1992)。
1.层序地层分析进展
层序地层学研究经历了十多年的发展,不仅仅是在其应用领域取得了重要的突破,更重要的是在研究理论和方法体系上得到了完善和补充,其突出进展体现在以下三方面:
其一是强制性海退概念的提出,海进和海退分别代表岸线朝陆地和朝海洋移动的方向。岸线移动方向是沉积物注入量与陆架上沉积物所能充填空间(可容纳空间)相平衡的函数。当相对海平面上升时,在沉积物注入量较低的地区表现为海进,在沉积物注入量较高的地区仍可以表现为海退。因此,相对海平面上升与海进之间并不存在一一对应关系。但是,当相对海平面下降时,海岸线只能发生海退,并不受沉积物注入量变化的影响。相对海平面下降引起的海退响应又被称为强制性海退(Posamentier等,1992)。Posa-mentier等(1999)对经典层序地层划分进行了修正,将具有强制性海退的沉积层序划分为四个体系域,即早期低位体系域、晚期低位体系域、海侵体系域和高位体系域。其中早期低位体系域由低位三角洲、低位楔组成,有些学者将其称为强制海退楔状体(forced re-gressive wedge)、强制海退楔形体系域(forced regressive wedge systems tract,FRWST)。与高位正常海退进积体系相比,强制海退形成的早期低位体系域有着其特有的一些特征:①形成相对较粗的楔状体,出现向海的一侧,与上一个层序的高位体系域之间被一个过水区分隔;②早期低位海岸沉积以近源陡、远源坡度逐渐变缓的滨面沉积为特征;③低位进积楔的顶面由许多小的阶梯状递降序列形成,表现为一个向海倾斜的斜坡;高位正常海退的进积单元顶面或者与下伏地层平行,或者只有一个较平缓的坡度。有了强制海退的概念,一些以前认为是古陆架砂体的沉积现在则可能被重新解释为早期低位海岸沉积。尤其是在缓斜坡地貌背景下,强制海退形成的早期低位三角洲、低位楔状体粒度粗,可作为有利的储集体;垂向上存在相的突变,海侵时期细粒泥岩直接覆盖其上,如果存在好的油源,该体系域则是很好的岩性油气藏类型。
其二是高精度层序地层学(high resolution sequence stratigraphy)的发展引起了人们的广泛关注并取得了显著的进展。近年来大量的研究表明,高精度层序地层学的概念和理论可有效地应用于地下地质的研究,为精细的地层对比、沉积相和储层特征等的研究提供了有效的分析方法和预测工具。Cross T.A.(1994)倡导的高分辨率层序地层学,代表着该领域的新进展。Cross领导的成因地层组建立了一种重要的高分辨率层序地层分析方法———基准面进退旋回分析。即:层序的发育是受全球海平面变化、构造沉降、沉积物供应量及沉积地形、气候等各种因素相互作用制约的基准面(Base Level)的变化控制的。基准面是一个相对于地球表面波状起伏的、连续的、略向盆地方向下倾的抽象面(非物理面),它不一定等于海平面、海平面向陆延伸的水平面或河流平衡剖面等,其位置、运动方向及升降幅度不断随时间而变化。基准面在变化中总是向其幅度最大值或最小值单方向移动,构成一个完整的上升与下降旋回。基准面的一个上升与下降旋回称为一个基准面旋回。一个基准面旋回是等时的,在一个基准面旋回变化过程中(时间域)保存下来的岩石为一个成因地层单元,即成因层序,它以时间为界面,因而是一个时间地层单元。基准面的旋回性变化,改变了沉积的可容纳空间,导致沉积与侵蚀的发生。当基准面高于地表时,有可容纳空间,以沉积作用为主,侵蚀作用成为局部和暂时的;当基准面低于地表时无可容纳空间,侵蚀作用占优势,沉积作用成为局部且暂时的;当基准面与地表一致时,既无沉积又不侵蚀,只有沉积物路过。在基准面随时间的变化过程中,地表不同地理位置上有四种地质作用状态:沉积、侵蚀、沉积物路过(无沉积)及饥饿沉积(密集段)。
其三是不同构造背景条件下所形成盆地的层序地层格架和层序样式,特别是构造活动型盆地。构造活动性盆地包括箕状断陷盆地(又称半地堑式盆地)、双断型断陷盆地、下断上坳型盆地、前陆盆地、克拉通盆地等,这些盆地的共同特点是构造活动作用于盆地形成和充填发育的全过程。目前,对这类盆地研究比较有代表性的有:Shanley和McCabe(1994)通过对阿根廷白垩纪地层的研究,建立了半地堑盆地的层序地层样式:Use strecher的河流-湖泊环境构造地层层序。近十年来,我国在陆相伸展盆地、前陆盆地的层序地层研究中也已取得了丰硕的成果。
2.高精度层序地层分析
在碳酸盐岩中,小尺度沉积旋回长期以来被解释为米兰柯维奇(Milankovitch)气候旋回引起的全球海平面高频变化的结果。而在硅质碎屑岩地层中,过去人们倾向于用自旋回,如三角洲废弃和河流改道,来解释小型的沉积旋回,因此,类似碳酸盐岩中的Mi-lankovitch气候旋回形成的硅质碎屑岩中的高频旋回未被重视和识别。层序地层学概念提出以后,Van Wagoner(1990)和Mitchum等(1991)认为在硅质碎屑岩地层中同样可以识别出代表全球海平面高频变化的小型尺度的沉积旋回,并提出了高频层序地层学的概念。他们在Exxon全球海平面旋回图的基础上,建立了硅质碎屑岩地层中能够识别出来的全球海平面旋回变化的谱系,其二级到五级的频率变化范围在9~10Ma,1~2Ma,0.1~0.2Ma和0.01~0.02Ma,其中9~10Ma和1~2Ma为Exxon全球海平面旋回图上的二级和三级旋回,由层序组和三级层序所代表。而高频旋回0.1~0.2Ma和0.01~0.02Ma可以在测井、岩心和露头上,或者是一些地区的高分辨率地震剖面上被识别出来,这些高频旋回处于Milankovitch气候旋回的范围内,是由于地球轨道的周期性变化产生的。在过去2Ma年中,这种气候旋回与高频海平面变化的关系已经被深海钻孔氧同位素研究所证实。
三级全球性海平面变化叠加上四级(120k.y.)和五级(50k.y.)海平面旋回产生复合的海平面旋回。这种复合的海平面旋回再叠加上构造沉降就成为一个地区的相对海平面变化。Van Wagoner等(1990)定义了两种类型的四级旋回(图1-1),分别被命名为“A”和“B”。“A”型四级旋回被定义为从海平面下降到海平面下降。如果有充足的物源供给,这种四级旋回产生四级层序,其边界为陆上不整合。叠加在四级旋回上的五级旋回产生准层序,其边界为海泛面。“B”型四级旋回被定义为从海平面快速上升到海平面快速上升。这种四级旋回只产生不完整层序,其边界为海泛面。在Van Wagoner等(1990)所用的例子中,当构造沉降量较大(15cm/k.y.以上),那么只有三级海平面旋回形成层序,而四级旋回形成准层序。当构造沉降量小于15cm/k.y.时,四级旋回可以形成层序,五级旋回形成准层序。由此可见形成四级层序的条件是:稳定的大地构造背景、盆地构造沉降极低和较高的沉积物供给速率。
图1-1 不同级别的海平面变化旋回的叠加和高频层序的发育
Mitchum等(1991)进一步明确了四级层序的概念及其构成样式(图1-2),在Exxon层序体系中,三级层序是由沉积体系域构成,而沉积体系域由准层序组构成,Mitchum等(1991)将这种层序称之为基本层序(fundamental sequence)(图1-2A)。如果体系域是由四级层序组、而不是由准层序组构成,这种层序被称之为复合层序(Com-posite Sequence)(图1-2B)。与四级层序类似,三级复合层序代表了全球海平面变化的一个旋回期间的沉积。在南非海域的盆地,三级层序的时限范围是1~3Ma,平均1.3Ma,而四级层序的时限范围在0.1~0.5Ma之间。所以四级层序与基本层序类似,也是由体系域组成,但是时限更短,有时在精度高的地震剖面上可以识别出来。这种层序也被Vail等(1991)称之为简单层序(Simple Sequence),形成于可容纳空间增加的速率太低,不能够抵消高频的全球性海平面下降时,高频海平面变化旋回就可以导致四级层序的发育。
图1-2 复合层序与基本层序构成模式
综上所述,四级层序具有三级层序的基本特征,但是时限更短,属于高频层序的范畴。四级层序概念的提出虽然有很长时间,但是其应用涉及地震分辨率。随着地震数据采集与处理技术的提高,国际上著名的石油公司和产业部门开始要求在油气勘探生产中要力求分出四级层序以便更精确地进行储层预测。
高频四级层序最早发现于北美大陆晚宾夕法尼亚世碳酸盐岩地层中,其中共划分出了至少55个旋回。随后,Browm等(1995)在南非Pletomos、Breadasdorp和Orange盆地白垩纪三角洲地层中识别出了时间在0.2~0.5Ma的四级层序,其体系域构成类似于三级层序。Henriksen等(1996)在挪威大陆架下渐新统海岸三角洲中划分出了38个发育良好的四级层序。
3.构造活动型盆地层序样式
相对于稳定被动大陆边缘,许多盆地不仅在形成时期受构造活动的控制,而且在发育过程中,构造仍然表现为强烈的活动性,如控盆或控洼断裂的强烈活动(活动的结果可以形成盆地复杂的古地貌形态),以及盆地充填过程中的基底快速沉降等等都可以是构造活动型盆地的特征,此类型盆地表现形式多种多样,可以为前陆盆地,也可以表现为伸展断陷盆地;同时也可能是多种构造活动叠加而成的盆地类型,如莺歌海盆地,在其演化过程中经历了多期构造性质不同的构造活动的叠加。
国内外许多学者已经对盆地内构造活动影响和控制着盆地的沉积体系发育以及沉积地层样式进行了研究,并提出相应的理论模型。Dawers和Underhill(2000)认为北海Statif-jord地区晚侏罗世断层活动以及断层与断层之间的相互作用控制了同断陷期的层序地层发育,沉积和沉降中心均随断层活动中心变动而迁移。李思田等(1998)研究结果也表明受走滑伸展构造运动的影响,莺歌海盆地的沉积、沉降中心依次向东南方向迁移。Gupta和Cowie(2000)提出,在伸展断陷盆地内,一方面断层活动(包括其形成、生长、消亡等)影响沉积体系的发育,另一方面沉积地层响应记录了不同阶段的构造活动特征。Gawthoper和Leeder(2000)对这一过程进行了总结,并分别研究了海相和陆相断陷盆地内构造活动的沉积响应特征以及三维演化模式。
与此同时,也有不少学者注意到构造活动对层序的影响,尤其是层序内部构成、垂向叠加样式、沉积体系分布等发育的控制作用。如Lin等(2001)对二连盆地的层序构成、沉积体系综合研究后,认为不同的构造阶段发育了截然不同的沉积体系,同时对应于不同类型的层序,不同级别的层序地层界面受控于不同级别的构造活动。Ravnas和Steel(1998)通过断陷盆地沉降速率与沉积特征之间的关系,提出不同的沉降速率(代表构造活动的强与弱)控制了体系域类型的发育以及体系域的内部构成样式。Corfield等(2000)则认为在同一盆地内同一时期不同的构造部位(即构造活动强度不同的地区)沉积地层样式、沉积构成也是不同的,如断层活动性强的部位(分段活动断层之中心),在断层上盘,地层表现为楔形,沉积、沉降中心位于断层面附近,而断层活动弱的部位(分段活动断层的末端),地层呈向斜形态,沉积、沉降中心明显向盆地方向迁移;并且强调这种差异特征不仅在露头上可以识别,在高分辨率地震剖面上同样可以识别。Young等(2000)则针对Suez裂谷中发育的构造转换带的沉积特征以及层序地层进行了精细研究,认为构造活动不仅控制了盆地的可容纳空间,而且沉积物供给速率也受其影响(或控制);高沉降速率、低沉积物供给,形成以退积和加积特征为主的准层序样式;低沉降速率、高沉积物供给,会形成粗粒的进积型准层序叠加样式。
对于一个应用于构造活动盆地边缘背景的总的层序模式而言,所有的控制层序发育的3个因素(全球性海平面变化、构造运动和物源供给)必须看作为独立的变量。海底的构造运动与全球性海平面变化共同确定了相对海平面的变化,进而,相对海平面的变化又确定了沉积物聚集的可容纳空间的大小。Yoon等(2003)将全球性海平面变化、构造运动和物源供给作为独立的变量考虑,提出了9种可能的层序发育情形(图1-3),均是可容纳空间变化和物源供给变化的函数。一些重复的情形被去除以后,以准层序的叠积样式为依据,可以确定出6种层序的类型。
图1-3 构造活动型盆地层序叠置样式
A型:A型层序由一个单独的SPG准层序组构成,与稳定正常的海退有关,是相对海平面稳定上升和沉积物高速供给并超过了新的可容空间的增加(S'>A'max)的结果(图1-3a)。Emery(1996)先前提出了一个类似的层序,由一个单一的体系域构成,命名为“海退体系域”。A型层序可以发育在正断层的上盘,形成于伸展断层的初始发育阶段,并伴有迅速的沉降和断层下盘的快速的沉积物堆积。
B型:B型层序可以形成于相对海平面稳定上升并叠加在中等物源供给的部位。物源的供给速率处于最大和最小可容纳空间之间(A'min<S'>A'max,图1-3b)。该层序与Exxon研究组提出的Ⅱ型层序基本类似,其中有2个SPG准层序组被解释为陆架边缘和高位体系域,夹在中间的RG准层序组被解释为海进体系域。B型层序可以发育在迅速沉降且具有中等物源供给的伸展盆地中。
C型层序:C型层序完全由RG叠置样式组成,与相对海平面的持续上升相伴,物源供给低(S'<A'min,图1-3c),这种情形发生在海岸平原,以一系列阶梯退缩为特征,沉积物聚集不足以填充逐渐增加的可容纳空间。
D型层序:由单一的楔形准层序组构成,具有OPG叠置样式,发育在强制性海平面稳定下降的部位(A'<0),是由于海底的迅速上隆超过了海平面的上升(图1-3d-f)。OPG准层序组与Hunt等(1992)的强制性海退体系域相对应。Emery(1996)也介绍了这种类型的层序,该层序由一个单一的中位体系域(midstand)组成,在这种体系域中任何时刻的沉降都不会高过沉积物的供给,从而产生海退。Hunt等(1992)认为强制性海退OPG准层序组由丘形盆底扇/裙组成,但是不一定同时发育在一个单一的海退期间。D型层序发育在低构造沉降或上隆、高沉积物供给的盆地中,活动断层的下盘或者低坡度的前积和陆架边缘发育最好。
E型层序:该类层序由两个交替的准层序组组成,具有SPG或OPG叠置样式,且与波动的相对海平面变化相伴生,海平面变化具有上升阶段和下降阶段,而且,沉积速率大于可容纳空间产生的速率(S'>A'max)(图1-3g)。E型层序边界与Ⅰ型层序边界极其类似,因为相对海平面旋回被相对海平面下降阶段所强化。然而,存在的争议是层序边界是放在强制性海退体系域之下,还是之上。在本文的模式中,边界被置于OPG准层序组(强制性海退体系域)之上,因为这是在深部地震剖面上最容易识别和容易制图的界面。
F型层序:该层序由下到上可以划分出4个准层序组:SPG,RG,SPG,OPG。相对海平面变化具有上升和下降阶段,中等到低的物源供给,并小于最大的可容空间产生的速率(S'<A'max)(图1-3h和i)。这基本上是Exxon研究小组定义的Ⅰ型层序。相对海平面上升的早期阶段,最下部的SPG样式的准层序组对应于Exxon研究组的低位域(低位楔),或者是Hunt(1992)的低位进积楔状体。在相对海平面迅速上升时期,沉积体系满足条件S'<A',这是由于相对海平面上升加速的缘故,由此产生了RG准层序组与海进体系域相对应。在这一海侵单元之上,另一个SPG准层序组发育,这是由于可容空间产生较慢并相对地被物源供给超越。这是在强制性海平面上升的晚期形成的高位体系域,最上部的OPG样式的准层序组是相对海平面下降时期的低位域(盆底扇)上半部分的沉积。
对于应用于构造活动型盆地边缘的层序模式,构造活动和沉积物的输入量是除了全球性海平面变化之外的主要因素,这个模式提出了6种层序类型,显示了一种准层序组在数量上和叠置样式上的变化。这些准层序组反映了可容纳空间的变化(A')和沉积物(S')供给速率的复合效应(图1-3)。由单一的SPG准层序组构成的A型层序与可容纳空间的稳定产生和沉积物的高速供给(S'>A')相关。B型层序(SPG-RG-SPG)本质上与Exxon研究组的Ⅱ-型层序类似,形成于可容纳空间稳定产生叠置在中等物源供给速率的条件下(A'max>S'>A'min);C型层序(RG)与可容纳空间的稳定产生和低沉积物供给有关(A'max>S');D型层序由具有OPG叠置样式的准层序组构成,由于海底的迅速上隆超过了全球海平面上升,可容纳空间负增长(A'<0);E型层序由SPG和上覆OPG准层序组构成,它形成于具有形成和消失阶段的、波动的可容纳空间的变化与沉积速率大于可容纳空间产生的最大速率的部位(S'>A'max);F型层序与Exxon研究组的Ⅰ-型层序相对应,从下到上被划分为4个准层序组SPG、RG、SPG和OPG,其与可容纳空间产生、消失的波动变化和中低物源供给小于可容纳空间产生的最高速率的部位相关(S'<A'max)。
位于弧后的Ulleung盆地西南边缘层序研究表明,盆地层序的发育与盆地边缘的构造演化密切相关(Young等,2000)。在弧后盆地扩张晚期阶段(中新世早—中期)迅速的沉降导致了可容纳空间的持续产生和构造控制的盆地边缘的高速和最后的中等速率的物源供给,导致了A、B型层序的发育;在弧后闭合的初期阶段伴有区域性收缩变形(中到晚中新世),较慢沉降速率和来自上隆冲断带沉积物的输入的增加导致了E型层序的发育。在晚中新世,Ulleung盆地的收缩调整是平静的或不活动的,边缘沉降加速,因此上升为主的相对海平面变化与中-低沉积物供给相配合产生了F型层序,反映了全球性海平面变化的主体控制机制。