牛庄洼陷沙三中亚段三角洲—重力流体系沉积特征与模式

周学文，姜在兴，汤望新，许文茂，胡晨林，贾超尘

1.海相储层演化与油气富集机理教育部重点实验室，北京　100083 2.中国地质大学(北京)能源学院，北京　100083 2.中国地质大学(北京)科学研究院，北京　100083

0　引言

深水重力流沉积拥有巨大的油气勘探前景，目前研究主要基于对现代海床未固结岩芯的描述、古代露头实测观察和水槽模拟实验研究[1]。不断进步的物理、数值模拟和海底实地监测等手段[2]为研究深水重力流沉积过程，流体性质及其对储层性质的影响提供了便利[3]。如今，在中国深水沉积领域，砂质碎屑流[4]，黏结、弱黏结和非黏结性碎屑流[5]等新概念与新模式已广泛运用于解释陆相湖盆深水块状砂岩搬运与沉积过程，如鄂尔多斯盆地延长组砂质碎屑流向浊流转化[6]，渤海湾盆地济阳坳陷水下重力流成因研究[7]，东营凹陷三角洲供给型重力流沉积模式[8]等。大量研究证实，我国陆相湖盆深水重力流的沉积特征与分布规律和大型湖泊三角洲密切相关。然而，针对陆相断陷湖盆三角洲—重力流体系，一方面对重力流类型，沉积过程与控制因素的研究仍不够深入；另一方面，三角洲砂体对深水沉积单元沉积特征与分布规律的影响尚缺乏深入认识；此外，不同分类方案所确定的重力流类型之间的对应关系不明确，容易造成术语冗余，从而阻碍不同地区深水沉积单元的对比研究。

本文以东营凹陷牛庄洼陷沙三中亚段三角洲—重力流体系为例，综合运用地震、测井和岩芯，在明确该区层序格架的基础上，系统研究三角洲砂体发育规律，重力流类型，沉积特征，分布规律和控制因素。旨在统一该区重力流的分类方案，完善大型湖泊三角洲—重力流沉积模式，以期为陆相断陷湖盆深水沉积研究提供借鉴。

1　区域地质概况

牛庄洼陷位于我国东部渤海湾盆地东营凹陷东部，东西长约40 km，南北宽约15 km，面积约600 km2[9]。其北面与民丰洼陷、利津洼陷以中央断裂背斜带相隔，南面以陈官庄—王家岗断裂构造带为界并与南部斜坡相接，西面与博兴洼陷纯化草桥断裂鼻状构造带相隔，向东延至广利南洼陷，形成一个南北界受东西向断裂所控制的似菱形沉积单元(图1)。

牛庄洼陷构造演化主要经历四个阶段：裂陷初期(E1-2k)、裂陷中期(E2s1-3)、裂陷晚期(E2s4～E3d)以及后裂谷期(N)[10]。

牛庄洼陷古近纪地层自下至上依次发育孔店组(E1-2k)、沙河街组(E2-3s)和东营组(E3d)。沙河街组细分为沙四段(E2s4)、沙三段(E2s3)、沙二段(E2-3s2)和沙一段(E3s1)，其中沙三段又进一步细分为沙三上(E2s3上)、沙三中(E2s3中)、沙三下(E2s3下)三个亚段。沙三中亚段(E2s3中)沉积时期，三角洲砂体多期次进积，重力流砂体分布广泛于洼陷斜坡与深陷带。

2　层序地层格架

沙三中亚段(E2s3中)时期，东营三角洲主力砂体进积到了牛庄洼陷内。为了精细刻画各期三角洲砂体的展布规律，前人对东营三角洲层序格架进行了大量研究，但至今仍没有形成统一的方案。邱桂强等[11]认为沙三中亚段高位体系域的发育程度远超低位和湖侵体系域；方勇等[12]将东营三角洲划分为7 个短期基准面旋回；操应长等[13]提出将东营凹陷沙三段内缺失稳定湖平面阶段的层序称为T—R层序；陈秀艳等[14]认为整个沙三中亚段是一个T—R层序，并将东营三角洲划分为8 个期次。刘建平等[15]则将东营凹陷沙三中亚段划分为一个完整的三级层序，其中细分出10 个准层序组。

通过分析牛庄洼陷地震反射终止样式(图2A)，68 口钻井的岩性和测井曲线特征(图2B,C)，笔者认为该区湖平面的变化确实缺乏相对稳定的阶段，所以将牛庄洼陷沙三中段划分为一个T—R层序，即发育快速湖侵和湖退两个体系域，并进一步细分为6 个准层序组，分别对应研究区6 期三角洲砂体。湖侵体系域对应准层序组PS6，湖退体系域包括准层序组PS1—PS5(图3)。

图1　东营凹陷大地构造位置及研究区位置图Fig.1　Tectonic setting of Dongying depression and location of the study area

图2　牛庄洼陷层序格架划分依据A.层序划分地震反射终端识别标志；B.准层序组界面识别标志(褐色油页岩)；C.准层序组界面识别标志(灰色泥岩向上突变为绿色泥岩)Fig.2　Identifications of dividing sequence stratigraphy in Niuzhuang sag

图3　王24井层序地层格架划分Fig.3　Column of sequence stratigraphy division of Well Wang 24

3　三角洲和重力流砂体沉积特征

沙三中亚段时期，牛庄洼陷古地势东南部较西北部高，在物源体系和该古地貌背景控制下，三角洲砂体自东南向西北多期次进积，使得洼陷内沉积了丰富的三角洲和重力流砂体。

3.1 三角洲砂体特征

基于工区地震、测井和岩芯特征，认为牛庄洼陷三角洲砂体主要属于三角洲前缘亚相，砂体类型主要为水下分流河道和河口坝砂体。水下分流河道砂体多具正粒序，岩性为分选好的中细砂岩至泥质粉砂岩，且发育块状层理(图4A)，槽状交错层理，冲刷充填构造，反映远物源、长距离搬运沉积作用，自然电位曲线上表现为钟型和箱型(图4C)。河口坝砂体多具反粒序，岩性为砂质纯、分选好的细砂岩和粉砂岩，发育槽状交错层理，浪成砂纹层理和生物钻孔等沉积构造(图4B)，反映湖水波浪冲刷改造和簸选作用，自然电位曲线上表现为漏斗形(图4D)。

3.2　重力流砂体特征

三角洲前缘砂体滑动、滑塌形成重力流沉积是三角洲建设—破坏作用的典型过程。迄今为止，在洼陷斜坡和深陷带的沙三中亚段，尤其是湖退体系域中已经发现了大量与重力流相关的岩性油气藏。早期对该地区重力流的研究主要集中在浊积扇类型[16-17]和储层特征[17-20]。

3.2.1存在的问题

近几年，在砂质碎屑流，砂质块体搬运沉积(SMTD)，黏结—弱黏结—非黏结性碎屑流等深水沉积新概念与模式的引领下，众多学者对渤海湾盆地东营凹陷重力流体系提出了新的观点和认识[8,21-23]。然而，大量重力流新类型和新模式的建立却容易造成术语的冗余，这不但会给后续的研究者增添困惑，而且一定程度上会阻碍不同区深水沉积单元的对比研究。

图4　三角洲砂体岩芯沉积构造和测井曲线特征A.牛106井，2 596.91 m，水下分流河道砂体块状层理构造；B.牛85井，2 800.8 m，河口坝砂体浪成沙纹层理和生物钻孔；C.牛106井水下分流河道箱型测井曲线；D.牛85井河口坝漏斗形测井曲线Fig.4　Characteristics of sedimentary structure and well logs of deltaic sand body

目前，国内外的学者对重力流的分类方案仍没有达成统一的共识。最具代表性的争论在于高密度浊流与砂质碎屑流的区别[4,24-25]。表1中列举了近几年国内学者对东营凹陷沙三中亚段重力流沉积的分类方案。陈秀艳等[14]的分类方案中将具滑塌构造的深水沉积单元也纳入了浊积岩范畴，这与浊流沉积典型沉积特征相矛盾；杨田等[21]和鲜本忠等[8]分类方案都比较注重流体成因，他们都接受了Shultz[26]将碎屑流细分为砂质碎屑流与泥质碎屑流的观点，但是在泥流沉积(含有漂浮砾石的块状泥岩[8])的归属上却产生了分歧，杨田等[21]认为泥质碎屑流包含泥流，鲜本忠等[8]则更倾向于将泥流沉积与泥质碎屑流和砂质碎屑流沉积区分开。然而，泥流和泥质碎屑流沉积的不同主要体现在黏土含量的细微差别上，在实际研究中往往较难区分；陈杰等[22]和王伟峰等[23]将三角洲或扇三角洲前缘斜坡带以滑动、滑塌或碎屑流的形式在斜坡带内发生再搬运形成的一种事件性砂质沉积体定义为坡移扇，这种命名的问题在于：1)运用相模式术语来描述深水沉积类型；2)无法准确反映深水砂体的成因；3)造成术语的冗余。

在实际研究中，重力流沉积的分类和命名应当注意以下几点：1)避开与相模式相关的术语；2)能够直接表征深水沉积的成因；3)在不同的地区具备可操作性。基于上述几点原则，结合牛庄洼陷沙三中亚段深水背景下钻井岩芯，测井和录井特征，划分出4 种不同成因的重力流砂体：滑动沉积、滑塌沉积、碎屑流沉积和浊流沉积(图5)。笔者不建议对滑塌和泥质碎屑流成因的深水沉积进一步细分，原因在于过于精细的分类方案在复杂多变的深水沉积领域往往不具备可操作性，在空间上也不便于对深水砂体展布特征的刻画和预测。

3.2.2滑动沉积

滑动沉积是指块体在自身重力或一定外界条件(地震、火山喷发、海啸、气旋、快速沉积作用等)的触发下，沿着滑动剪切面向下滑动形成的块体流沉积[27-28]。砂质滑动沉积自身由于弹性形变的力学特性，往往保留块体滑移前的原始沉积结构而不发生明显的内部变形，但是其底部的剪切作用会使下伏半固结的沉积物发生变形。牛庄洼陷滑动沉积顶底都为深水暗色泥岩，岩性主要为灰色或灰白色的粉砂岩和细砂岩，在岩芯上可见明显的滑移面，底部剪切带和下伏沉积物变形等现象(图6)。一些砂质滑动沉积内部可见交错层理，波状层理等浅水沉积沉积构造，在自然伽马测井曲线上特征为中高幅指形和漏斗形，推测来自三角洲前缘。

3.2.3滑塌沉积

滑塌沉积是指沉积物在斜坡沿下凹的滑移面运动时，在重力和剪切力的作用下，内部发生了明显旋转变形的黏结块体[8,27-28]。滑动沉积具备塑性形变的力学特性，常发育大量的变形和滑塌构造。牛庄洼陷滑塌沉积岩性主要为灰色粉砂岩和具棕色氧化色泥质粉砂岩；岩芯上常见包卷层理，砂岩脉，旋转火焰状构造，砂岩中变形泥砾，砂岩中的变形层，滑塌褶皱等变形构造(图7)；在自然伽马测井曲线上特征为中高幅齿化钟型。

3.2.4碎屑流沉积

碎屑流是一种塑性流变，层流态，且通过整体“冻结”的方式发生沉积的非牛顿流体[24,28]。碎屑流沉积中研究最热的属砂质碎屑流沉积，其流动特征介于黏性和非黏结性碎屑流之间[4]。“砂质碎屑流”最早由Hampton[29]基于水下碎屑流模拟实验提出，随后Shanmugam[28]对砂质碎屑流的概念和模式进行完善和推广。如今，砂质碎屑流和深水块体搬运模式已被国内众多学者接受并用于解释我国陆相湖盆深水块状砂岩成因[8,21,27,30-32]。

表1　东营凹陷深水重力流分类方案对比Table 1　Comparison of different classifications of deep-water gravity flow in Dongying depression

图5　牛庄洼陷沙三中亚段重力流类型Fig.5　Types of gravity flow in the middle Member 3 of Shahejie Formation in Niuzhuang sag

图6　牛庄洼陷滑动沉积岩芯照片A.块状细砂岩，内部高角度次级滑移面，牛33井，3 110.4 m；B.块状粉砂岩，滑动砂体底部主滑移面使下伏泥质粉砂岩褶皱变形，牛33井，3 063 m；C.波状层理粉细砂岩，底部剪切带使下伏泥质粉砂岩发生变形，牛33井，3 084.5 m；D.粉砂质泥岩，具擦痕的滑移面，牛22井，3 213.45 mFig.6　Cores of slide deposits in Niuzhuang sag

牛庄洼陷沙三中亚段也发育了一定规模的碎屑流沉积，且主要为砂质碎屑流沉积。其岩性为与暗色或棕色泥岩伴生的块状中细砂岩、细砂岩和粉砂岩；岩芯上可见大量泥岩撕裂屑，砂岩顶部常有漂浮泥砾，发育面状碎屑结构，与底面常突变接触(图8)；在自然伽马曲线上主要呈中高幅齿化箱形和齿形，也有少数呈钟形，可能是由块状细砂岩顶部富集泥砾造成。

图7　牛庄洼陷滑塌沉积岩芯照片A.泥质粉砂岩中肠状的砂岩脉，牛22井，3 206.5 m；B.粉细砂岩中的变形层，牛22井，3 228.5 m；C.具滑塌构造的泥质粉砂岩，牛22井，3 208.0 m；D.粉砂质泥岩中的包卷层理，牛22井，3 225.9 m；E.泥质粉砂岩中旋转火焰状构造，牛22井，3 207.10 m；F.粉砂岩中变形泥砾，牛33井， 3 125.7 m；G.粉砂岩中变形层，牛116井，3 101.2 m；H.粉砂岩中的滑塌褶皱，牛116井，3 101.4 mFig.7　Cores of slump deposits in Niuzhuang sag

图8　牛庄洼陷碎屑流沉积岩芯照片A.顶部具有泥岩撕裂屑块状砂岩，牛33井，3 079.8 m；B.含大量泥岩撕裂屑的灰白色块状砂岩，牛85井，3 066 m；C.含两种颜色泥岩撕裂屑的灰色块状砂岩，牛22井，3 211.2 m；D.具面状碎屑结构的块状砂岩，泥砾的长轴与层面平行，牛22井，3 215.85 mFig.8　Cores of debris flow deposits in Niuzhuang sag

3.2.5浊流沉积

浊流是一种遵循牛顿流变学的，紊态的非均质沉积物重力流，其沉积物靠紊流支撑，悬浮搬运，并且以悬浮沉降的方式沉积[27-28,33-34]。Lowe[35]基于颗粒粒度，浓度和支撑机理将浊流分为低密度浊流和高密度浊流，且认为高密度浊流以紊流、颗粒碰撞和流体逃逸混合支撑。高密度浊流最初被Kuenen[36]用于解释深水粒序层理形成，随后又被用于解释深水块状砂岩的成因[37-38]。然而，高密度浊流概念自提出以来就一直是深水沉积领域争论的焦点。最具代表性的争论在于高密度浊流和砂质碎屑流这两个概念都可用于解释深水块状砂岩[39]。Haughtonetal.[24]认为碎屑流是按整体冻结的方式沉积，而浊流是按持续加积的方式沉积，两者的沉积方式有着显著差别。Shanmugam[28]认为所谓的“高密度浊流”实际上就是砂质碎屑流，同时指出浊流不可能存在于高浓度值的情况下，因为当颗粒浓度超过30%时，浊流会转化为层流。Zouetal.[6]提出紊流态沉积物流的颗粒体积分数大于20%时，会转变为颗粒碰撞方式，从而形成弱黏结性碎屑流。在高密度浊流问题上，笔者更赞同Haughton[24]的观点，建议从成因上识别区分砂质碎屑流和高密度浊流。砂质碎屑流和高密度浊流沉积在岩芯上都可以含有漂浮的泥砾和泥岩撕裂屑，但是整体冻结方式沉积的砂质碎屑流不会形成高密度浊流沉积中的粒序层理和牵引构造。

牛庄洼陷沙三中亚段的浊流沉积是由低密度浊流形成，岩芯上未见典型高密度浊流沉积特征。低密度浊流沉积的岩性为灰色薄层粉砂岩和泥质粉砂岩；岩芯上主要发育鲍马序列和火焰状构造(图9)；在自然伽马曲线上常呈低幅的指形和齿形。

图9　牛庄洼陷浊流沉积岩芯照片A.火焰状构造，牛116井，3 102.1 m；B.鲍马序列Tc、Td和Te段，牛33井，3 071.0 m；C.鲍马序列Ta、Tc和Td段，牛106井，3 205.4 mFig.9　Cores of turbidites in Niuzhuang sag

4　砂体分布规律和控制因素

4.1　砂体分布规律

牛庄洼陷沙三中亚段时期，三角洲自东南向西北方向多期次进积。为了总结各期三角洲砂体与不同成因重力流沉积的平面分布规律，在地震资料的约束下，结合岩芯观察和不同成因砂体在测井上的响应特征，分析各准层序组的砂体展布特征(图10)。

PS6—PS1沉积时期三角洲主要发育水下分流河道和河口坝两种砂体，砂体平面分布的范围逐渐扩大，与此同时，在滨浅湖和半深湖—深湖亚相中形成了滑动—浊流沉积等一系列不同成因的重力流砂体。重力流砂体在湖退体系域即PS5—PS1准层序组中分布广泛，但在湖侵体系域PS6分布局限。平面上，滑动—浊流沉积的分布具有一定规律性，滑动沉积距离三角洲前缘近，呈带状平行于三角洲前缘砂体分布；滑塌和碎屑流沉积距离三角洲前缘较远，主要分布于洼陷斜坡和斜坡角；浊流沉积发育的距离最远，呈席状分布于洼陷的深湖相中。

4.2　控制因素

牛庄洼陷沙三中亚段重力流沉积主要来自快速进积的三角洲前缘砂体，其形成主要受沉积，构造和古气候这3 大因素的控制。

4.2.1沉积因素

三角洲在沙三中亚段湖退体系域PS5—PS1准层序组时期，从洼陷东南部的斜坡向西北方向洼陷的深陷部位依次进积。邱桂强等[11]对牛庄洼陷三角洲的沉积速率进行过研究，认为该三角洲的进积速率可达200 m/万年，且在400 万年时间内进积超过80 km。高进积速率必然造成其三角洲前缘的砂体在垂向上迅速地叠置，尚未固结的前缘砂体容易在地震、风暴等因素的诱发下，沿斜坡向半深湖—深湖发生块体搬运。

4.2.2构造因素

东营凹陷是受断层控制的伸展半地堑，其南部斜坡主要受八面河断裂带控制。牛庄洼陷北接中央断裂带，西临石村断裂带，南部边界与陈官庄—王家岗断裂带相接(图1)。在古近纪时期，整个东营凹陷经历了强烈的伸展活动[40]。通过统计计算东营凹陷新生代北西向和北东向断层的活动速率(图11)，认为牛庄洼陷周缘八面河断裂带、中央断裂带和石村断裂带等控洼断裂的活动高峰主要出现在沙三段(E2s3)时期。对洼陷南部斜坡陈官庄—王家岗断裂带不同地质时期生长断层断层落差的计算结果显示(表2)，该反向调节式断裂带在沙三中亚段时期断层落差最大，断裂活动明显加强[40,42]。如此频繁的地震和断裂活动会促使三角洲前缘砂体失稳向洼陷深处发生滑动，滑塌，形成一系列不同成因的重力流砂体。方勇等[12]对洼陷东南部的构造坡折带进行了研究，认为沙三中亚段时期，发育两种构造坡折带：与陈官庄—王家岗断裂密切相关的断裂坡折带和受基底挠曲控制的同沉积地形坡折带。顺三角洲物源方向发育的构造坡折带有利于三角洲前缘重力流砂体的搬运与再沉积作用。

图10　牛庄洼陷沙三中亚段三角洲—重力流砂体展布图Fig.10　Distribution of deltaic and gravity flow sand bodies in the middle Member 3 of Shahejie Formation in Niuzhuang sag

4.2.3古气候因素

温暖潮湿的气候条件下，三角洲平原河流水系茂盛，物源供给充足，三角洲的进积速率较高，有助于重力流沉积的形成。前人对渤海湾盆地济阳坳陷和东营凹陷进行了大量的古气候分析工作(表3)，结果一致显示沙三中亚段时期，气候温暖潮湿，主要属亚热带温湿气候类型。

图11　牛庄洼陷周边断裂带断层活动率随地质年代变化直方图(数据来自赵延江[41])Fig.11　Histogram of fault movement rate versus geologic time in the peripheral fault zone of Niuzhuang sag(data from Zhao[41])

地质时期E2s3下E2s3中E2s3上～E2⁃3s2下E2⁃3s2上～E3s1E3d下E3d上断层落差/m508050302030

表3　渤海湾盆地沙三中亚段古气候研究结果对比表Table 3　Comparison of the results of paleoclimate studies of the middle Member 3 of Shahejie Formation in Bohai Bay Basin

沉积岩中不同类别的孢粉组合可以直接反映古植物群落特征，它易于保存，常用于恢复古温度，地史时期干燥度或降水量和古地形[46]，是对古气候重建最有利的证据之一。对东营凹陷沙三中亚段孢粉组合的鉴定统计结果(图12)表明，该时期亚热带植物含量远大于山地植物和旱生植物。亚热带植物具有地质温度计作用，其含量与气温的变化密切相关，含量越高则指示古气温越高。旱生植物多为灌木或草本植物，其含量的纵向变化可反映古湿度变化[46]。

图12　东营凹陷沙三中亚段孢粉生态组合饼状图(数据来自朱宗浩等[46])Fig.12　Pie charts of sporopollen biotic associations of the middle Member 3 of Shahejie Formation in Dongying depression(data from Zhu et al.[46])

所以，沙三中亚段时期亚热带植物尤为发育，旱生植物极其匮乏的孢粉组合特征指示该时期的古气候条件温暖潮湿，有利于河流三角洲发育。

湖泊沉积物地球化学记录是恢复古环境的有效手段之一，湖泊沉积碳酸盐碳氧稳定同位素作为重要的古环境变化指标，常用于研究古湖水化学(盐度)和古气候(温度、降水及季节变化)[48-49]。东营凹陷沙河街组湖泊沉积碳酸盐的δ13CPDB和δ18OPDB值在沙三段时期均发生了明显的负偏移(图13)。δ18O值主要受控于湖水的蒸发作用，当气候湿润，降水充沛，湖水蒸发量小于水体的补给量时，流域河流与地下水会为湖泊带来大量贫18O的水体，从而造成δ18O值减小[50-51]。δ13C值与湖泊水体深度和水体与大气CO2交换程度相关[50-51]。表层水体易与大气沟通，使δ13C值升高。深层水体相对封闭，死亡有机质的沉降和降磷作用使水体相对富集12C，从而导致δ13C值降低[51]。基于上述原理，沙三段时期，东营凹陷δ13CPDB和δ18OPDB值的负漂移现象表明该时期气候条件相对湿润，降水较充足，湖水的补给量大于蒸发量，水深相对增加，这与孢粉组合分析结果相吻合。

在此温暖潮湿古气候条件的控制下，研究区沙三中亚段半深湖相的岩芯中多见含植物碎片块状细砂岩和棕褐色的泥岩(图14)，推测是季节性洪水将来自陆地或浅水区的沉积物同三角洲前缘砂体沿构造斜坡一起带入半深湖—深湖相的泥岩中。

5　沉积模式

基于地震、测井、岩芯和各准层序组三角洲—重力流砂体的平面分布特征，结合重力流沉积的控制因素，总结研究区三角洲—重力流体系的沉积模式(图15)。

图13　东营凹陷沙河街组湖泊沉积碳酸盐δ13CPDB、δ18OPDB变化曲线(数据来自蔡观强等[50])Fig.13　Variations of δ13CPDB and δ18OPDB values of carbonates for lacustrine sediments in Shahejie Formation, Dongying depression(data from Cai et al.[50])

牛庄洼陷三角洲主体位于洼陷东南缓坡带，且主要受陈官庄—王家岗断裂构造带控制。沙三中亚段时期，气候温暖潮湿，三角洲砂体从洼陷的东南缓坡不断地向洼陷内进积，在PS6—PS1准层序组时期，三角洲多期次进积使得水下分流河道和河口坝砂体不断在洼陷东南斜坡堆积。在地质历史演化的过程中，三角洲前缘砂体在自身重力和外部构造作用力的影响下，内部沉积物发生移动和垮塌，顺着前缘斜坡向洼陷内搬运，形成一系列不同成因的重力流砂体。滑动沉积主要呈带状分布在三角洲前缘砂体外侧；在自身重力、外部构造，风暴，季节性洪水等因素的影响下，沉积物会顺斜坡继续向洼陷深处搬运，一部分沉积物会在斜坡上发生变形、旋转并堆积下来，形成滑塌沉积；另一部分沉积物则会在斜坡角附近形成碎屑流沉积；而粒度最细的粉砂、泥质粉砂的沉积物最终在洼陷深陷部位悬浮沉降形成席状的浊流沉积。

图14　牛庄洼陷沙三中亚段半深湖相中含植物碎屑块状砂岩和棕褐色泥岩A.含植物碎屑块状细砂岩，牛22井，3 220.7 m；B.含植物碎屑块状中砂岩，牛22井，3 208.9 m；C.半深湖亚相中的棕褐色泥岩段，牛22井，3 221.0 mFig.14　Massive sandstones with plant detritus and brown mudstones in the middle Member 3 Shahejie Formation semi-deep lake facies, Niuzhuang sag

图15　牛庄洼陷沙三中亚段三角洲—重力流体系沉积模式Fig.15　Depositional model of delta and gravity flow system of the middle Member 3 of Shahejie Formation in Niuzhuang sag

6　结论

(1) 牛庄洼陷沙三中亚段可划分为一个T—R层序，其中包括一个快速湖侵体系域和一个湖退体系域，并进一步细分为6 个准层序组，分别对应6 期三角洲砂体。

(2) 牛庄洼陷沙三中亚段发育四种成因的重力流砂体：滑动沉积，滑塌沉积，碎屑流沉积和浊流沉积。重力流砂体在湖退体系域即PS5—PS1准层序组中分布广泛，但在湖侵体系域分布局限。平面上，滑动沉积主要呈带状临近三角洲前缘砂体分布；滑塌和碎屑流沉积分布于洼陷斜坡和斜坡角附近；浊流沉积呈席状分布于洼陷深陷部位。

(3) 三角洲高速的进积速率，频繁的地震与断裂活动，顺物源方向发育的构造坡折带和温暖潮湿的古气候条件控制着牛庄洼陷重力流砂体的形成与分布。

参考文献(References)

[1]Piper D J W, Normark W R. Processes that initiate turbidity currents and their influence on turbidites: A marine geology perspective[J]. Journal of Sedimentary Research, 2009, 79(6): 347-362.

[2]朱筱敏，李顺利，潘荣，等. 沉积学研究热点与进展：第32届国际沉积学会议综述[J]. 古地理学报，2016，18(5)：699-716. [Zhu Xiaomin, Li Shunli, Pan Rong, et al. Current hot topics and advances of sedimentology: A summary from 32nd IAS meeting of sedimentology[J]. Journal of Palaeogeography, 2016, 18(5): 699-716.]

[3]Marchand A M E, Apps G, Li W, et al. Depositional processes and impact on reservoir quality in deepwater Paleogene reservoirs, US Gulf of Mexico[J]. AAPG Bulletin, 2015, 99(9): 1635-1648.

[4]Shanmugam G. High-density turbidity currents: are they sandy debris flows?[J]. Journal of Sedimentary Research, 1996, 66(1): 2-10.

[5]Talling P J, Masson D G, Summer E J, et al. Subaqueous sediment density flows: depositional processes and deposit types[J]. Sedimentology, 2012, 59(7): 1937-2003.

[6]Zou C N, Wang L, Li Y, et al. Deep-lacustrine transformation of sandy debrites into turbidites, upper Triassic, central China[J]. Sedimentary Geology, 2012, 265-266: 143-155.

[7]Cao Y C, Yang T, Wang Y, et al. Classification of subaqueous sedimentary gravity flows and recognition of their deposit in lacustrine basin with example of Jiyang depression, East China[C]//Abstracts Book of 19th International Sedimentological Congress. Geneva: University of Geneva, 2014: 118.

[8]鲜本忠，王璐，刘建平，等. 东营凹陷东部始新世三角洲供给型重力流沉积特征与模式[J]. 中国石油大学学报(自然科学版)，2016，40(5)：10-21. [Xian Benzhong, Wang Lu, Liu Jianping, et al. Sedimentary characteristics and model of delta-fed turbidites in Eocene eastern Dongying depression[J]. Journal of China University of Petroleum(Natural Science Version), 2016, 40(5): 10-21.]

[9]蔡李梅. 东营凹陷牛庄洼陷沙三中亚段砂岩透镜体油藏成藏动力学过程分析[D]. 武汉：中国地质大学(武汉)，2009. [Cai Limei. Dynamic process of hydrocarbon migration and accumulation of lenticular sandstones of the third middle member of Shahejie Formation in Niuzhuang subsag, Dongying depression[D]. Wuhan: China University of Geosciences (Wuhan), 2009.]

[10]朱光有，金强，周建林，等. 渤海湾盆地东营断陷湖盆充填模式研究[J]. 石油实验地质，2003，25(2)：143-148，152. [Zhu Guangyou, Jin Qiang, Zhou Jianlin, et al. Study on the infilling patterns of lacustrine basins in the Dongying depression[J]. Petroleum Geology & Experiment, 2003, 25(2): 143-148, 152.]

[11]邱桂强，王居峰，李从先. 东营凹陷沙三中东营三角洲地层格架与油气勘探[J]. 同济大学学报，2001，29(10)：1195-1199. [Qiu Guiqiang, Wang Jufeng, Li Congxian. Preliminary study on stratigraphy architecture of Middle-Shasan Dongying delta and its significance to hydrocarbon exploration[J]. Journal of Tongji University, 2001, 29(10): 1195-1199.]

[12]方勇，邓宏文，王红亮，等. 东营三角洲高频层序特征与岩性圈闭分布[J]. 中国海上油气(地质)，2003，17(3)：160-163. [Fang Yong, Deng Hongwen, Wang Hongliang, et al. High frequency sequence and lithologic trap distribution in Dongying delta[J]. China Offshore Oil and Gas (Geology), 2003, 17(3): 160-163.]

[13]操应长，姜在兴，夏斌，等. 陆相断陷湖盆T-R层序的特点及其控制因素：以东营凹陷古近系沙河街组三段层序地层为例[J]. 地质科学，2004，39(1)：111-122. [Cao Yingchang, Jiang Zaixing, Xia Bin, et al. Characters and controlling factors of T-R sequence in lacustrine deposits of rift basin: An example from the Dongying depression, eastern China[J]. Chinese Journal of Geology, 2004, 39(1): 111-122.]

[14]陈秀艳，姜在兴，杜伟，等. 东营凹陷沙三中亚段东营三角洲沉积期次成因及对含油性的影响[J]. 沉积学报，2014，32(2)：344-353.[Chen Xiuyan, Jiang Zaixing, Du Wei, et al. Origin of depositional cycles and their influence on oil-bearing sandstone of Dongying delta in Mid-Es3, Dongying depression[J]. Acta Sedimentologica Sinica, 2014, 32(2): 344-353.]

[15]刘建平，鲜本忠，王璐，等. 渤海湾盆地东营凹陷始新世三角洲供给型重力流地震沉积学研究[J]. 古地理学报，2016，18(6)：961-975. [Liu Jianping, Xian Benzhong, Wang Lu, et al. Seismic sedimentology of delta-fed turbidites of the Eocene in Dongying sag, Bohai Bay Basin[J]. Journal of Palaeogeography, 2016, 18(6): 961-975.]

[16]王居峰. 济阳坳陷东营凹陷古近系沙河街组沉积相[J]. 古地理学报，2005，7(1)：45-58. [Wang Jufeng. Sedimentary Facies of the Shahejie Formation of Paleogene in Dongying sag, Jiyang depression[J]. Journal of Palaeogeography, 2005, 7(1): 45-58.]

[17]侯明才，田景春，陈洪德，等. 东营凹陷牛庄洼陷沙三中段浊积扇特征研究[J]. 成都理工学院学报，2002，29(5)：506-510. [Hou Mingcai, Tian Jingchun, Chen Hongde, et al. Turbidite fan characters of the intermediate section of Member 3 of Shahejie Formation in Niuzhuang depression of Dongying area[J]. Journal of Chengdu University of Technology, 2002, 29(5): 506-510.]

[18]侯明才，田景春，陈洪德，等. 牛庄洼陷沙三中亚段浊积砂体储集特征研究[J]. 石油勘探与开发，2003，30(1)：93-95. [Hou Mingcai, Tian Jingchun, Chen Hongde, et al. A study on reservoir characteristics of the turbidite sand body Shasan intermediate member in Niuzhuang sag in Dongying depression[J]. Petroleum Exploration and Development, 2003, 30(1): 93-95.]

[19]饶孟余，钟建华，王海侨，等. 山东东营牛庄洼陷沙三中亚段浊积砂体储层特征及影响因素[J]. 现代地质，2004，18(2)：256-262. [Rao Mengyu, Zhong Jianhua, Wang Haiqiao, et al. Reservoir characteristics and influence factors of the turbidite sandbody Shasan intermediate member in the Niuzhuang sag, Dongying depression, Shandong province[J]. Geoscience, 2004, 18(2): 256-262.]

[20]饶孟余，张遂安，李秀生. 牛庄洼陷沙三中亚段浊积岩储层成岩作用及主控因素分析[J]. 地质找矿论丛，2007，22(1)：66-70. [Rao Mengyu, Zhang Sui’an, Li Xiusheng. Diagenesis and controlling factors of the turbidite reservoir of sub-middle of the third member of Shahejie Formation in Niuzhuang sag[J]. Contributions to Geology and Mineral Resources Research, 2007, 22(1): 66-70.]

[21]杨田，操应长，王艳忠，等. 深水重力流类型、沉积特征及成因机制：以济阳坳陷沙河街组三段中亚段为例[J]. 石油学报，2015，36(9)：1048-1059. [Yang Tian, Cao Yingchang, Wang Yanzhong, et al. Types, sedimentary characteristics and genetic mechanisms of deep-water gravity flows: A case study of the middle submember in member 3 of Shahejie Formation in Jiyang depression[J]. Acta Petrolei Sinica, 2015, 36(9): 1048-1059.]

[22]陈杰，于正军，韩宏伟，等. 东营凹陷坡移扇形成机制与分布规律[J]. 特种油气藏，2015，22(6)：45-49. [Chen Jie, Yu Zhengjun, Han Hongwei, et al. Forming mechanism and distribution pattern of slope slip fan in Dongying sag[J]. Special Oil & Gas Reservoirs, 2015, 22(6): 45-49.]

[23]王伟锋，胡瑜，于正军，等. 东营三角洲前缘坡移扇储集体特征及成因研究[J]. 石油实验地质，2016，38(5)：600-608. [Wang Weifeng, Hu Yu, Yu Zhengjun, et al. Reservoir characteristics and genesis of slope fans in Dongying delta front[J]. Petroleum Geology & Experiment, 2016, 38(5): 600-608.]

[24]Haughton P, Davis C, McCaffrey W, et al. Hybrid sediment gravity flow deposits-classification, origin and significance[J]. Marine and Petroleum Geology, 2009, 26(10): 1900-1918.

[25]裴羽，何幼斌，李华，等. 高密度浊流和砂质碎屑流关系的探讨[J]. 地质论评，2015，61(6)：1281-1292. [Pei Yu, He Youbin, Li Hua, et al. Discuss about relationship between high-density turbidity current and sandy debris flow[J]. Geological Review, 2015, 61(6): 1281-1292.]

[26]Shultz A W. Subaerial debris-flow deposition in the upper Paleozoic cutler formation, western Colorado[J]. Journal of Sedimentary Research, 1984, 54(3): 759-772.

[27]刘芬，朱筱敏，李洋，等. 鄂尔多斯盆地西南部延长组重力流沉积特征及相模式[J]. 石油勘探与开发，2015，42(5)：577-588. [Liu Fen, Zhu Xiaomin, Li Yang, et al. Sedimentary characteristics and facies model of gravity flow deposits of late Triassic Yanchang Formation in southwestern Ordos Basin, NW China[J]. Petroleum Exploration and Development, 2015, 42(5): 577-588.]

[28]Shanmugam G. New Perspectives on Deep-Water Sandstones, Volume 9: Handbook of Petroleum Exploration and Production[M]. Amsterdam, Boston: Elsevier, 2012, 9: 1-542.

[29]Hampton M A. The role of subaqueous debris flow in generating turbidity currents[J]. Journal of Sedimentary Petrology, 1972, 42(4): 775-793.

[30]邹才能，赵政璋，杨华，等. 陆相湖盆深水砂质碎屑流成因机制与分布特征：以鄂尔多斯盆地为例[J]. 沉积学报，2009，12(6)：1065-1075. [Zou Caineng, Zhao Zhengzhang, Yang Hua, et al. Genetic mechanism and distribution of sandy debris flows in terrestrial lacustrine basin[J]. Acta Sedimentologica Sinica, 2009, 12(6): 1065-1075.]

[31]李相博，刘化清，完颜容，等. 鄂尔多斯盆地三叠系延长组砂质碎屑流储集体的首次发现[J]. 岩性油气藏，2009，21(4)：19-21. [Li Xiangbo, Liu Huaqing, Wanyan Rong, et al. First discovery of the sandy debris flow from the Triassic Yanchang Formation, Ordos Basin[J]. Lithologic Reservoirs, 2009, 21(4): 19-21.]

[32]李相博，付金华，陈启林，等. 砂质碎屑流概念及其在鄂尔多斯盆地延长组深水沉积研究中的应用[J]. 地球科学进展，2011，26(3)：286-294. [Li Xiangbo, Fu Jinhua, Chen Qilin, et al. The concept of sandy debris flow and its application in the Yanchang Formation deep water sedimentation of the Ordos Basin[J]. Advances in Earth Science, 2011, 26(3): 286-294.]

[33]Dott Jr R H. Dynamics of subaqueous gravity depositional processes[J]. AAPG Bulletin, 1963, 47(1): 104-128.

[34]Middleton G V, Hampton M A. Sediment gravity flows: mechanics of flow and deposition[C]//Middleton G V, Bouma A H. Turbidity and Deep water Sedimentation. Anaheim, CA: SEPM, SEPM Pacific section Short Course, 1973: 1-38.

[35]Lowe D R. Sediment gravity flows: Ⅱ. Depositional models with special reference to the deposits of high-density turbidity currents[J]. Journal of Sedimentary Petrology, 1982, 52(1): 279-297.

[36]Kuenen P H. Properties of turbidity currents of high density[C]//Hough J L. Turbidity Currents and the Transportation of Coarse Sediments to Deep Water: A Symposium. SEPM Special Publication, 1951, 2: 14-33.

[37]Middleton G V. Experiments on density and turbidity currents: Ⅲ. Deposition of sediment[J]. Canadian Journal of Earth Sciences, 1967, 4(3): 475-505.

[38]Pickering K T, Hiscott R N, Hein F J. Deep Marine Environments: Clastic Sedimentation and Tectonics[M]. London: Unwin Hyman, 1989: 416.

[39]Stow D A V, Johansson M. Deep-water massive sands: nature, origin and hydrocarbon implications[J]. Marine and Petroleum Geology, 2000, 17(2): 145-174.

[40]叶兴树，王伟锋，戴俊生，等. 东营凹陷沙三—东营期断裂活动特征[J]. 中国石油大学学报(自然科学版)，2006，30(4)：7-11. [Ye Xingshu, Wang Weifeng, Dai Junsheng, et al. Characteristics of fault activities of Sha-3 Member and Dongying periods in Dongying depression[J]. Journal of China University of Petroleum, 2006, 30(4): 7-11.]

[41]赵延江. 东营凹陷古近系盆地结构与充填特征研究[D]. 广州：中国科学院研究生院(广州地球化学研究所)，2007. [Zhao Yanjiang. The research of basin structure and filling characteristics of Palaeogene in Dongying depression[D]. Guangzhou: Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, 2007.]

[42]张英凯，谭红霞，李庆红. 东营凹陷的主要断裂活动性分析[J]. 山东国土资源，2009，25(7)：29-31. [Zhang Yingkai, Tan Hongxia, Li Qinghong. Analysis on major fault activities in Dongying depression[J]. Land and Resources in Shandong Province, 2009, 25(7): 29-31.]

[43]石油化学工业部石油勘探开发规划研究院. 渤海沿岸地区早第三纪孢粉[M]. 北京：科学出版社，1978. [Planning & Research Institute of Petroleum Exploration and Development. Early Tertiary sporopollen in Bohai coastal area[M]. Beijing: Science Press, 1978.]

[44]叶得泉，钟筱春，姚益民. 中国油气区第三系(I)：总论[M]. 北京：石油工业出版社，1993. [Ye Dequan, Zhong Xiaochun, Yao Yimin. Tertiary oil and gas zones in China (I): Pandect[M]. Beijing: China Petroleum Industry Press, 1993.]

[45]姚益民，梁鸿德，蔡治国，等. 中国油气区第三系(IV)：渤海湾盆地油气区分册[M]. 北京：石油工业出版社，1994. [Yao Yimin, Liang Hongde, Cai Zhiguo, et al. Tertiary oil and gas zones in China (IV): Volume of oil and gas area in Bohai Bay Basin[M]. Beijing: China Petroleum Industry Press, 1994.]

[46]朱宗浩，范乃敏，罗华仪. 胜利油区第三纪气候[J]. 山东地质，2000，16(1)：22-28. [Zhu Zonghao, Fan Naimin, Luo Huayi. Climate of Shengli oil field in Tertiary[J]. Land and Resources in Shandong Province, 2000, 16(1): 22-28.]

[47]李守军，王明镇，郑德顺，等. 山东济阳坳陷古近纪的气候恢复[J]. 山东科技大学学报(自然科学版)，2003，22(3)：6-9. [Li Shoujun, Wang Mingzhen, Zheng Deshun, et al. Recovery of climate of Palaeogene in Jiyang depression of Shandong[J]. Journal of Shandong University of Science and Technology: Natural Science, 2003, 22(3): 6-9.]

[48]王苏民，羊向东，马燕，等. 江苏固城湖15ka来的环境变迁与古季风关系探讨[J]. 中国科学(D辑)，1996，26(2)：137-141. [Wang Sumin, Yang Xiangdong, Ma Yan, et al. Environmental change of Gucheng Lake of Jiangsu in the past 15 ka and its relation to palaeomonsoon[J]. Science in China: Series D, 1996, 39(2): 144-151.]

[49]王苏民，张振克. 中国湖泊沉积与环境演变研究的新进展[J]. 科学通报，1999，44(6)：579-587. [Wang Sumin, Zhang Zhenke. New progress of lake sediments and environmental changes research in China[J]. Chinese Science Bulletin, 1999, 44(6): 579-587.]

[50]蔡观强，郭锋，刘显太，等. 东营凹陷沙河街组沉积岩碳氧同位素组成的古环境记录[J]. 地球与环境，2009，37(4)：347-354. [Cai Guanqiang, Guo Feng, Liu Xiantai, et al. Carbon and oxygen isotope characteristics and palaeoenvironmental implications of lacustrine carbonate rocks from the Shahejie Formation in the Dongying sag[J]. Earth and Environment, 2009, 37(4): 347-354.]

[51]孙晶，薛林福，李双林，等. 渤海湾盆地孔南地区孔二段湖盆性质研究[J]. 地球科学进展，2012，27(S)：408-410. [Sun Jing, Xue Linfu, Li Shuanglin, et al. Study on the lacustrine basin’s characteristics of second member of Kongdian Formation in Kongnan area, Bohai Bay Basin[J]. Advances in Earth Science, 2012, 27(S): 408-410.]