0 引言
1 地质背景
莺歌海盆地为新生代伸展走滑型含油气沉积盆地［2］，位于中国南海北部大陆架西部地区，在海南岛以西与越南之间，形态上呈NNW走向的长菱形，面积约为11.3×104 km2［15］，盆地经历了古近系的断陷和新近系的坳陷两大演化阶段，平面上主要划分为莺西斜坡、河内凹陷、临高凸起、中央凹陷和莺东斜坡5个构造单元［3］［图1（a）］；盆地周缘物源水系复杂，但主体为北部红河物源、西部越南昆嵩隆起及东部海南隆起三大物源体系，为莺歌海盆地储层发育提供了充足的物质条件［7，16］，其中，海南岛方向对莺歌海盆地影响较大的是昌化江水系［3］。综合前人认识和控盆断层的活动特征，将莺东斜坡带分为北、中、南3段［17］，北段新近系沉积薄，埋藏较浅［4］。盆地自下而上发育崖城组、陵水组、三亚组、梅山组、黄流组、莺歌海组和乐东组［图1（b）］。
研究区A构造位于莺歌海盆地中央凹陷北部，莺东斜坡北段［图1（a）］。前人研究表明，莺歌海盆地发育中新世陆源海相烃源岩，为Ⅱ—Ⅲ型干酪根［19］，厚达4~6 km的中新统是莺歌海盆地主要的气源岩，同时，莺东斜坡带烃源广泛，发育众多油气苗［17］。黄流组沉积时期，盆地东南部乐东地区水深，西北部临高地区水浅［18］，岩性组合以浅灰色砂岩、灰质砂岩、生物灰岩、白垩系砂岩及灰质泥岩不等厚互层为主［3］。黄流组的高位体系域和低位体系域发育近源三角洲，黄流组一段沉积期，莺东斜坡带发育来自海南岛物源的高位三角洲［3，17］。研究区黄流组一段沉积相带以三角洲前缘亚相为主［图1（c）］，三角洲前缘亚相内部微相分异较大。各种成因类型的三角洲砂体及浊积岩砂体为莺北区主要储层［3］；上新统莺歌海组二段的浅海—半深海相泥岩为主要盖层［15，20-22］。莺东斜坡带弱垂向运移动力决定了天然气垂向运移的可能性较小，因此天然气运聚以侧向运聚为主［23］。莺东斜坡区、临高—海口区油气运移以侧向为主，运移通道为区域性海退—低水位期形成的区域不整合面和碎屑岩输导层［15］，油气运聚特征反映出该区带是莺歌海盆地内天然气大规模侧向运聚的有利指向区［4，17］。
2 海平面变化与体系域划分
2.1 单井可容空间特征
区域研究表明，黄流组沉积时期海平面震荡上升，发育3个海侵旋回，第一个海侵对应层序T31—T40，第二个海侵与第三个海侵对应层序T30—T31 ［24］；前人依据层序地层学划分方案和区域地层对比，认为黄流组一段为一个三级层序地层单元，进一步划分出低位域（LST）、海侵体系域（TST）和高位域（HST）［5，11］。综合钻井岩性组合，认为HK29-A井在一个长期旋回内，可划分出3个水进水退的中期旋回；黄流组一段下部为厚层粉砂岩夹薄层泥岩的组合，属于低位域沉积；中部为厚层泥岩夹薄层砂岩，为海侵域沉积；上部为数米厚的粉砂岩与泥岩互层，到黄流组顶部变为厚层泥岩，为高位域沉积（图2）。最大海泛面处为约10 m厚的泥岩，GR值最高，AC值低。
高频海平面变化幅度可以用可容空间平均累计厚度偏移（Fischer曲线）来表示［25］，是经线性沉降校正之后把米级旋回层序的累计厚度与时间的变化作图而得出的一种图解方法［26］。有学者用Fischer曲线来定义长周期海平面变化，分析可容空间变化［26-27］。近年来，有学者将Fischer曲线旋回分析技术应用于旋回地层对比并取得了很好的效果［28-29］。因此，可用已钻井Fischer曲线代表可容空间变化，反映相对海平面变化。研究区HK29-A井位于中央凹陷［图1（a）］，钻遇完整黄流组，可用于分析目的层段海平面变化过程。通过对HK29-A井目的层段进行自然伽马曲线去奇异值、滑动平均、一阶差分及旋回厚度统计，可得到该井反映的Fischer曲线，即可容空间变化曲线。该曲线波动表明，黄流组一段沉积早期可容空间较小，向上增大，随后减小，到黄流组沉积末期再次增大（图2），反映了海平面升高—降低—升高的过程。
通常可容空间较大时，表示海平面上升。综合HK29-A井岩性组合、海平面变化曲线、可容空间变化曲线，得到该井最大海泛面（MFS）深度为2 490 m（图2），根据这一深度，开展井震标定，可从地震剖面上开展地震同相轴追踪，得到最大海泛面上超点位置。
2.2 地震Wheeler域旋回特征
图4 过HK29-A井与HK30-1井的地震Wheeler域旋回及A构造的地震Wheeler域旋回对比（a） HK29-A井附近地震Wheeler域剖面特征［剖面位置见图1（a）中B—B’］；（b）HK30-1井附近地震Wheeler域剖面特征［剖面位置见图1（a）中B—B’］；（c）A构造附近地震Wheeler域剖面特征［剖面位置见图1（a）中A—A’］ Fig.4 Comparison of seismic Wheeler domain cycles passing through Wells HK29-A and HK30-1， as well as structure A |
通过HK29-A井与目标区的过井地震剖面，并根据地震剖面的Wheeler旋回剖面，对比目标区Wheeler域旋回与HK29-A1井的层序划分［图4（c）］可见，A构造黄流组同样发育一个水进—水退旋回，沉积旋回特征与HK30-1井附近一致。
综上所述，A构造黄流组一段沉积于高位域晚期，发育一个水进—水退沉积旋回；同时，由于研究区总体为三角洲沉积背景［11］，故水退期三角洲为更为有利的砂体发育区。
3 复合朵体识别与期次划分
在油气勘探开发的不同阶段，针对沉积体的研究尺度和研究对象都有所不同。自然界的沉积体系通常是由不同级次、不同规模的沉积单元在空间上复合叠置而成，从而在不同阶段需要根据不同层次开展沉积体研究，揭示不同尺度下的沉积演变过程。
在三角洲、海底扇研究中，不同学者提出了朵叶体不同的级次划分方案。李彦泽等［35］在对浅水三角洲研究中提出叠覆式浅水三角洲模式［图5（a）］，指出这类三角洲的基本构成单元是三角洲朵叶体，而并非传统的三角洲微相单元。即该类三角洲以单一的沉积朵体为沉积单元，多个同级别朵体叠置构成高级别复合朵体，从而形成三角洲骨架系统，且不能细分为其他微相。在海底扇研究中，DEPTUCK等［36］将朵叶体系划分为朵叶复合体、复合朵叶、朵叶单元、层—层组4个级次；PRELAT等［37］把朵叶体系划分为朵叶复合体、单一朵叶、朵叶单元、层—层组4个级次；张佳佳等［38］将海底扇朵叶体系细分为朵叶复合体组合、朵叶复合体、单一朵叶、朵叶单元、朵叶层、单一岩相6个级次［图5（b）］。而林煜等［14］采用地震相分析、地震沿层切片、露头观察描述等多种手段相结合的方法，分级次开展了浊积朵叶的构型模式研究，确定了深水浊积朵叶的7级构型单元划分方案。
从而在晚期高位域格架约束下，开展研究区黄流组复合朵体识别。有学者研究表明，对于海相沉积盆地，在90°相位地震剖面上，同相轴能对应砂岩和泥岩，沉积地质解释相对容易［39］，因此，在研究区90°地震剖面上，自目的层段底界开始，自下而上，观察三角洲特有的地震同相轴“S”形反射或叠瓦状反射特征，设置地震极大值和具有这些反射特征的一条同相轴为一个前积复合朵体。
然后将地震同相轴的上下相邻2个0°相位作为其顶底界面，当该同相轴尖灭、消失、断开或与其他同相轴相交时结束追踪，并将追踪结束的点设置为前积复合朵体的边界（图6）。
4 砂体富集区预测
通过追踪闭合后可编制复合朵体厚度图，结合复合朵体剖面地震内幕结构和外部形态特征预测复合朵体砂体发育区。复合朵体厚度图［图9（a）］显示，复合朵体中部厚度最大，向边缘减薄，表明复合朵体中部为低洼处，有利于粗碎屑堆积。从复合朵体上游至下游剖面可见，顺物源的NE—SW方向存在叠瓦状前积反射特征［图9（b）］，每一期叠瓦状反射即为一期单一朵体，表明该复合朵体是由多期单一朵体叠覆构成的复合体；上游垂直物源的方向，下切现象明显［图9（c）］，表明上游微古地貌限制性较强，导致水动力较强，是朵体的输砂通道；而下游垂直物源方向复合朵体厚度较薄，表明下游微古地貌限制性减弱，单一朵体间的切割作用减小，复合朵体中的砂体较薄，粒度较细。
有学者利用最小振幅属性分析认为，水道边部席状砂强能量，最小振幅值高；水道微相以中—细砂为主，最小振幅值中等［40］。从而，结合研究区最小属性分带特征（图10），编制了复合朵体E砂体厚薄程度分布图［图11（a）］。过复合朵体上游与下游剖面显示，上游的内幕朵体下切幅度大于下游下切幅度［图11（b）］，表明复合朵体上游水动力较强，下游水动力较弱。综合层序旋回特征分析认为，研究区黄流组为晚期高位域时期沉积的三角洲；其中，复合朵体E为海退期进积形成，砂体富集程度最高；朵体内幕叠置、下切等特征能指示砂体集中发育区；利用层间最小振幅属性结合复合朵体剖面分析认为，复合朵体E中部为砂体集中发育区［图11（a）］。通过对有利砂体富集区进行预测，可为有利储层分布预测、勘探井位部署、储量计算等提供基础支撑。
5 结论
（1）莺歌海盆地莺东斜坡带北段A构造黄流组沉积时期，海平面先上升后下降，早期可容空间大，晚期可容空间小；黄流组形成于高位域晚期，发育高位域晚期沉积的三角洲砂体，可进一步划分为一个海进—海退的旋回。
（2）研究区可识别出5期三角洲复合朵体，其中复合朵体A为最早沉积的朵体，复合朵体E为最晚沉积的朵体；由复合朵体A至复合朵体D为海平面上升期发育的三角洲退积型复合朵体，复合朵体E为海平面下降期发育的三角洲进积型复合朵体，砂体富集程度最高。
（3）复合朵体E内部发育多期单一朵体迁移叠置，朵体上游下切作用强、单一朵体相互切割叠置，下游单一朵体以横向排列为主；朵体中部单一叠置、下切特征最明显，指明复合朵体E中部为砂体集中发育区，呈南北向条带状。
甘公网安备 62010202000678号
