-
辫状河油气储层在新生代陆相盆地中广泛发育,心滩作为其最重要的沉积单元成为油气储层研究的热点[1-5] 。 针对心滩内部结构,国内外学者开展了大量研究,涵盖露头、现代沉积和地下储层[6-11] 。但主要是对心滩及其内部结构的静态刻画,对其形成过程的研究甚少[8-12] ,尤其是对地下心滩沉积的形成过程,缺乏深入认识。因此选取黄骅坳陷典型辫状河沉积实例,分析心滩储层构型并利用沉积数值模拟再现心滩的形成过程,认识心滩形成机制,指导地下储层构型刻画。
-
1 研究区地质概况
-
研究区位于黄骅坳陷羊二庄油田(图 1),是黄骅坳陷赵家堡断层下降盘的一个穹隆背斜构造。 基底之上自下而上发育古近系孔店组、沙河街组和东营组、新近系馆陶组和明化镇组以及第四系平原组(图2)。 其中馆陶组和明化镇组是该油田的主要含油层系[13] 。研究区馆陶组发育辫状河沉积,心滩是最主要的含油砂体类型[14] 。经过40 余年的开发,羊二庄油田进入高含水期。心滩砂体在开发中表现出较强的沉积非均质性,开发层内矛盾突出,故认识心滩内部结构成为预测剩余油分布和提高采收率的关键。
-
图 1 研究区位置图
-
Fig. 1 Location of study area
-
图 2 黄骅坳陷地层综合柱状图
-
Fig. 2 Comprehensive stratigraphic column of Huanghua Depression
-
2 解剖区心滩构型特征
-
辫状河河道迁移频繁,具有砂体连片发育的特点,砂体以心滩沉积为主,呈“叠覆泛砂体冶的形态展布[15] 。 按照 Miall 的河流相储层构型划分方法[2,16] ,单一心滩沉积为四级构型单元,而心滩内部的加积单元为三级构型单元。
-
2.1 单一心滩识别与特征
-
地下储层辫状河沉积多是多期次叠加改造的残留沉积,其心滩沉积多以心滩复合体的形式存在,识别单一心滩是刻画心滩内部结构的前提和基础。研究区为开发后期密井网区,利用井点静态和开发动态资料可以有效识别单一心滩。
-
利用同一心滩砂体在测井相上的相似性可以初步对其进行识别。 此外心滩底部的相对高程差也是单一心滩识别的一个依据。 通常情况下,在一定范围内河床高程相对稳定,因此单一心滩的底部与标准层的相对距离具有稳定性,可以作为识别单一心滩的依据。单一心滩内部同期的落淤层沉积厚度横向变化符合自然规律,这也是识别单一心滩的一个辅助依据。
-
根据上述识别方法和原则,综合考虑心滩发育规模和平面形态,对解剖区的单一心滩进行识别。 以图 3 所示的 Ng31-3 单砂层典型心滩为例,在代表区域泛滥沉积的稳定泥岩标准层拉平之后,庄 5- 15-6 井心滩底面深度为 1885郾 29 m,与邻近的庄 5- 15-2、庄 5-15-3、庄 5-15-5、庄 5-15-10、庄 5-15- 3K 以及庄 5-15 井同期次心滩砂体底面深度相近。 而在相邻的庄 G3 井,同一单砂层心滩底界面深度为 1912.34 m,与前述的几口井深度差异大,因此庄 G3 井与前述各井钻遇的心滩不是同一个单一心滩。同时从落淤层发育情况看,庄 5-15-6 井与庄 G3 井落淤层发育位置及指示的垂积体厚度差异显著(图3(a)、(c))。
-
图 3 解剖区心滩发育特征
-
Fig. 3 Sedimentary features of braided bars in detailed research area
-
研究区 40 余年的开发积累了大量动态资料,尤其是示踪剂监测数据,可以作为判断砂体连性和单一心滩识别的辅助依据。 以庄 5 -15 -8 井组为例,注水井庄 5-15-8 井注入示踪剂后,在位于同一心滩的监测井庄 5-17K 中有示踪剂显示,指示同一心滩内部砂体的良好连通性;而由于不同心滩之间四级构型界面的渗流屏障作用,位于邻近心滩的监测井庄 5-15-2 井未见到示踪剂显示。
-
2.2 心滩内部构型特征
-
以落淤层为界面的砂岩增生体叠置心滩模式是中国油气储层研究中使用最广的心滩沉积模式,落淤层由辫状河洪泛事件末期悬浮落淤沉积形成,具有指示事件性沉积的意义[4-5,11,15] 。 研究区心滩内落淤层以夹层的形式出现,其岩性以粉砂质泥岩和泥岩为主,在测井曲线上特征显著。以图 3 中庄5- 15-5 井所在的心滩为例,测井曲线上可以识别出4期落淤层沉积,在自然电位( SP) 和自然伽马(GR) 测井曲线上与围岩幅度差明显,孔隙度和渗透率均出现低值特征(图4)。
-
在解剖区选取完整的单一心滩,在心滩的头部、中部和尾部各选一口井建立连井剖面,在落淤层单井识别的基础上,结合落淤层沉积特点进行不同井点的落淤层对比(图 3、5)。根据高程和单期垂向加积砂体的厚度变化趋势,庄5 -15 -2 井位于心滩主体部位,发育 5 期落淤层,而位于心滩头部和尾部的庄5-15-5 和庄 5-15-6 井均识别出 4 期落淤层沉积。位于心滩底部的一期落淤层规模较小,形成于心滩沉积早期,对应的其下部心滩垂积体仅覆盖庄5-15-6 井和庄5-15-2井;随着心滩发育规模的增大,在顺流加积作用下,心滩顺流加积生长,向下游推进,在 3 口井上均钻遇心滩砂体。心滩沉积过程中,在向下游顺流加积的同时,心滩头部迎水面受到水流侵蚀作用,导致心滩头部增生砂体和落淤层难以沉积或保存,因此位于心滩头部的庄5-15 -6井顶部一期落淤层不发育。
-
图 4 庄 5-15-5 井心滩内落淤层及心滩增生体
-
Fig. 4 Slit layer and accretion sand body in braided bars of well Zhuang 5-15-5
-
图 5 心滩内部结构(井点平面位置见图 3)
-
Fig. 5 Reservoir architecture of braided bar (planar location of wells is in Fig. 3)
-
3 心滩发育过程数值模拟
-
从研究区心滩构型来看,除了垂向加积外,在心滩形成过程中还存在侧向和顺流向的沉积作用。地下储层中保留的心滩形态及其内部结构是上述作用综合的结果。 研究通过沉积过程的数值模拟,重现心滩形成过程,探究不同的沉积作用如何塑造心滩砂体沉积。 基于泥沙水动力学的沉积过程数值模拟能够有效避免基于沉积物理模拟和现代沉积观测的沉积过程研究中存在的不足[17-19] ,越来越多地被应用于沉积过程和沉积机制的研究[17-23] 。
-
3.1 模型与模拟参数
-
根据实例解剖区辫状河发育规模,建立 3 km×20 km 的模拟河段,纵向上设定 3 层网格;根据实例区构型单元尺度和研究精度要求,设定平面网格尺寸为 50 m×50 m,确保单一心滩内有足够数量的网格显示其内部结构,以满足对模拟结果的观察和分析需要。根据现代沉积研究的观测尺度以及气候和降水量的周期性,同时考虑模拟运算量,将模拟时间设置为一年,模拟计算的时间步长设置为 2 min,输出观测记录的时间步长为 6 h。依据研究区古地貌恢复结果,将模型的地形坡度设置为 0.1%。
-
针对数值模拟中河流流量、含沙量等水文地质参数,采用参考沉积规模相似的现代河流相关参数的方法来确定。 根据研究区辫状河构型刻画结果, 以河道宽度、地貌(坡度)、沉积物类型和粒度、心滩规模和几何形态等为对比依据,选取与实例区具有很高相似性的黄河龙门水文站河段和潼关水文站河段,作为相似现代沉积参考河段。 根据近年来龙门水文站和潼关水文站的水文资料,设置模拟的入口边界流量为 40 m3 / s(单位网格),沉积物总质量浓度为9 kg / m3 。
-
3.2 心滩沉积演化分析
-
沉积过程模拟结果显示,在辫状河沉积开始阶段,沉积物在向下游搬运的过程中逐渐沉积下来,形成心滩的早期雏形。之后在垂向加积、顺流加积和侧向加积的作用下,心滩规模逐渐增大,同时心滩的形态也在不断的加积和侵蚀共同作用中产生变化和迁移(图 6)。
-
3.2.1 顺流加积与垂向加积
-
在心滩发育过程中,心滩的头部受到水流的冲刷侵蚀,而滩尾处则以沉积作用为主,出现顺流逐渐生长的特征,在宏观上表现为心滩的整体顺流迁移。 同时心滩顶部也在垂向加积作用下逐渐加厚(图 7)。以图7中标示的心滩为例,从模拟的 948 步到1168 步,根据平面上沉积物厚度随时间的变化可以看出,该心滩尾部顺流加积作用明显,在沉积物不断堆积的过程中形态从分叉状逐渐合并成一个完整的滩尾。心滩迎水面的侵蚀作用对心滩形态的塑造要弱于滩尾的加积作用,这一点在顺流向的纵向剖面图上更加显著。在模拟时间由948 步到 1168 步的过程中,滩头侵蚀长度为 141 m,迁移速度为 0.64m / 步,但滩尾加积范围达到 298 m,迁移速度达到1.35 m / 步(图8 纵轴正值代表沉积,负值代表侵蚀; 剖面位置见图 7 剖面线 AB;图中红色虚线框为观测的心滩位置)。 滩头侵蚀的速度和规模都要弱于滩尾的加积。
-
图 6 辫状河心滩发育过程的模拟平面结果
-
Fig. 6 Simulation result of braided river deposition process on plane
-
心滩的垂向加积在剖面上最为显著。以初始沉积厚度为基准,在上述沉积过程中,心滩最大沉积厚度增加了0. 3 m,从剖面上不同时间的砂体顶面高度可以看到(图8),垂向加积作用在心滩形成过程中具有普遍性,是主要的砂体沉积作用类型,但在不同阶段垂向加积的规模和速度存差异。
-
3.2.2 侧向侵蚀与加积
-
在心滩顺流加积和垂向加积的过程中,由于心滩两侧水道的流速和水流能量差异,流速较高的一侧水道对心滩冲刷,常出现以侵蚀为主;流速较小的一侧水道对心滩的侵蚀速度慢或者发生砂体的加积,两侧水道对心滩作用的这种差异在宏观上表现为心滩侧向迁移。 以图7中标示的心滩为例,从水流速度模拟结果的横剖面上可以看到(图9),该心滩两侧水道流速分别为 0.87 ~ 1. 09 和 0. 64 ~ 0.87 m / s。 左侧辫状水道流速大,对照该剖面上砂体厚度随时间的变化发现(图 10),心滩左侧以侵蚀作用为主,砂体厚度随时间减小,靠近心滩出现水道下切,而右侧水道的侵蚀作用较弱,从剖面上的砂体厚度变化来看,右侧砂体加积。
-
图 7 心滩顺流加积过程中的平面形态变化
-
Fig. 7 Planar evolution of braided bar with downstream accretion
-
图 8 顺流向纵剖面上心滩砂体顶面与河床底部距离随时间的变化
-
Fig. 8 Variation of distance between braided bar top and channel bottom with time at downstream section
-
图 9 河道横截面水流速度分布(剖面位置见图 7 剖面线 DC)
-
Fig. 9 Distribution of water velocity in cross profile of channel (location of the profile is marked in Fig. 7 as line DC)
-
图 10 辫状河横剖面上心滩顶面与河床底部距离随时间的变化
-
Fig. 10 Variation of distance between braided bar top and channel bottom with time at crosssection
-
从图 10 的心滩横剖面砂体厚度变化可以看到,在模拟时间从1 020 步到1140步的过程中,心滩左侧侵蚀宽度约为 26 m,右侧的侧向加积宽度约为6 m, 侵蚀速度大于加积速度,在该剖面位置心滩横向宽度变小,宏观上呈现向右侧的迁移。在侧向迁移的同时,综合图7和图8的砂体厚度变化可以看到,心滩也出现整体向下游方向的迁移。
-
4 解剖区心滩发育过程
-
通过研究心滩内部结构的刻画结果和心滩沉积过程的数值模拟,对解剖区心滩发育过程进行探讨。解剖区心滩在形成初期规模较小,心滩生长过程中经历多次洪水期,心滩规模逐渐扩大。 从保留下来的落淤层发育情况看,有明显地层记录的大规模洪水沉积事件有 5 期。在第一个大规模洪水期过后, 心滩的发育雏形已经明确,心滩在顺流加积和头部侵蚀作用下,宏观上表现为整体的不断向下游迁移。在迁移过程中,由于头部侵蚀速率小于尾部的加积速率,因此心滩整体在顺流方向上尺度变大(图5、11)。 在沉积数值模拟中也得到了相同的规律和定量化的依据。
-
心滩顶部的一期落淤层在侵蚀作用下发育范围局限,在庄5-15-6 井处未能完整保留,在测井曲线上无法明显识别(图 5);从纵剖面上看,第一期心滩规模小,这是在庄5-15 -5 井上没有看到第一期增生体顶部落淤层沉积的原因。其后在顺流加积作用下心滩顺流生长到庄5 -15-5 井位置,在该井处发育之后的几期增生砂体。从横剖面上看,心滩单侧加积作用明显,心滩头部和中部在右侧呈现出以侵蚀作用为主的特点,形成心滩两侧的不对称加积,在宏观上表现为心滩横向迁移(图 11)。
-
图 11 研究区心滩演化模式(井平面位置见图 3)
-
Fig. 11 Evolution model of braided bar in study area (well position is as in Fig. 3)
-
5 结论
-
(1)将动静态资料相结合,刻画研究区单一心滩及其内部构型,发现垂向加积、顺流加积、侧向加
-
积和侵蚀作用共同塑造心滩的形态和内部结构。
-
(2)心滩的形成是一个动态的过程,滩头的侵蚀与滩尾的顺流加积是一个同时的过程,但是二者
-
的作用强度会存在差异,侵蚀作用会造成滩头部位落淤层发育差。 心滩两侧河道的流速差异会造成心滩一侧加积而另一侧侵蚀,存在不对称的侧向加积。垂向加积伴随着心滩的整个生长过程,但是在不同阶段以及心滩的不同部位,垂向加积速率存在差异。
-
(3)辫状河水动力特征复杂,心滩的成因机制和影响因素多样,研究是对地下储层心滩发育的探
-
讨,两岸地貌差异、物源变化、气候等因素对心滩发育的影响没有单独进行详细探讨,这是下
-
一步需要关注的问题。此外基于泥沙水动力学的数值模拟多用在水利研究中,设计模拟时间尺度小于地质时间尺度,两个不同时间尺度上模拟结果的可对比性是沉积学家所关注的重点。 考虑到季节周期对河流沉积的控制,年度内的沉积模拟对认识辫状河沉积演化具有参考意义。
-
参考文献
-
[1] 刘海,林承焰,董春梅,等.基于地震沉积学的复杂曲流带沉积相研究[J].中国石油大学学报(自然科学版),2018,42(6):30-39.LIU Hai,LIN Chengyan,DONG Chunmei,et al.Sedi-mentary facies of complex meandering belt based on seis-mic sedimentology[J].Journal of China University of Pe-troleum(Edition of Natural Science),2018,42(6):30-39.
-
[2] MIALL A D.Architectural-element analysis:a new meth-od of facies analysis applied to fluvial deposits[J].Earth Science Reviews,1985,22(4):261-308.
-
[3] 张可,吴胜和,冯文杰,等.砂质辫状河心滩坝的发育演化过程探讨:沉积数值模拟与现代沉积分析启示 [J].沉积学报,2018,36(1):81-91.ZHANG Ke,WU Shenghe,FENG Wenjie,et al.Discus-sion on evolution of bar in Sandy Braided River:insights from sedi-ment numerical simulation and modern bar[J].Acta Sedimentologica Sinica,2018,36(1):81-91.
-
[4] 刘钰铭,侯加根,王连敏,等.辫状河储层构型分析 [J].中国石油大学学报(自然科学版),2009,33(1):7-11.LIU Yuming,HOU Jiagen,WANG Lianmin,et al.Ar-chitecture analysis of braided river reservoir[J].Journal of China University of Petroleum(Edition of Natural Sci-ence),2009,33(1):7-11.
-
[5] 刘海,林承焰,张宪国,等.孔店油田馆陶组辫状河储层构型及剩余油分布规律[J].吉林大学学报(地球科学版),2018,48(3):665-677.LIU Hai,LIN Chengyan,ZHANG Xianguo,et al.Reser-voir architecture and remaining oil distribution in braided river of Guantao Formation,Kongdian Oilfield[J].Jour-nal of Jilin University(Earth Science Edition),2018,48(3):665-677.
-
[6] 陈彬滔,于兴河,王天奇,等.砂质辫状河岩相与构型特征:以山西大同盆地中侏罗统云冈组露头为例[J].石油与天然气地质,2015,36(1):111-117.CHEN Bintao,YU Xinghe,WANG Tianqi,et al.Lithofa-cies and architectural characteristics of sandy braided riv-er deposits:a case from outcrops of the Middle Jurassic Yungang Formation in the Datong Basin,Shanxi Province [J].Oil &Gas Geology,2015,36(1):111-117.
-
[7] 陈仕臻,林承焰,任丽华,等.砂质辫状河沉积模式的建立:以委内瑞拉奥里诺科重油带H区块为例[J].沉积学报,2015,33(5):965-971.CHEN Shizhen,LIN Chengyan,REN Lihua,et al.Es-tablishment of the depositional model of Sandy Braided River:a case from the H Block in Orinoco heavy oil belt,Venezuela [J].Acta Sedimentologica Sinica,2015,33(5):965-971.
-
[8] 刘钰铭,侯加根,宋保全,等.辫状河厚砂层内部夹层表征:以大庆喇嘛甸油田为例[J].石油学报,2011,32(5):836-841.LIU Yuming,HOU Jiagen,SONG Baoquan,et al.Char-acterization of interlayers within braided-river thick sand-stones:a case study on the Lamadian Oilfield in Daqing [J].Acta Petrolei Sinica,2011,32(5):836-841.
-
[9] 侯加根,刘钰铭,徐芳,等.黄骅坳陷孔店油田新近系馆陶组辫状河砂体构型及含油气性差异成因[J].古地理学报,2008,10(5):459-464.HOU Jiagen,LIU Yuming,XU Fang,et al.Architecture of braided fluvial sandbody and origin for petroliferous difference of the Guantao Formation in Kongdian Oilfield of Huanghua Depression[J].Journal of Palaeogeography,2008,10(5):459-464.
-
[10] ZHANG X G,LIN C Y,ZHANG T.Seismic sedimen-tology and its application in shallow sea area,gentle slope belt of Chengning Uplift[J].Journal of Earth Sci-ence,2010,21(4):471-479.
-
[11] 孙天建,穆龙新,赵国良.砂质辫状河储集层隔夹层类型及其表征方法:以苏丹穆格莱特盆地Hegli油田为例[J].石油勘探与开发,2014,41(1):112-120.SUN Tianjian,MU Longxin,ZHAO Guoliang.Classifi-cation and characterization of barrier-intercalation in sandy braided river reservoirs:taking Hegli Oilfield of Muglad Basin in Sudan as an example [J].Petroleum Exploration and Development,2014,41(1):112-120.
-
[12] 秦国省,吴胜和,宋新民,等.远源细粒辫状河三角洲沉积特征与单砂体构型分析[J].中国石油大学学报(自然科学版),2017,41(6):9-19.QIN Guosheng,WU Shenghe,SONG Xinmin,et al.Sedimentary characteristics of distal fine-grain braided delta and architecture analysis of single sand body[J].Journal of China University of Petroleum(Edition of Natural Science),2017,41(6):9-19.
-
[13] 王松,林承焰,张宪国,等.羊二庄油田 Nm34-3 单砂层曲流河储层构型[J].大庆石油地质与发,2018,37(6):85-92.WANG Song,LIN Chengyan,ZHANG Xianguo,et al.Architectures of the Meandering River reservoir in Nm34-3 single-sandstone-layer of Yang' er Zhuang Oilfield [J].Petroleum Geology and Oilfield Development in Daqing,2018,37(6):85-92.
-
[14] 蒲秀刚,赵贤正,李勇,等.黄骅坳陷新近系古河道恢复及油气地质意义[J].石油学报,2018,39(2):16-171.PU Xiugang,ZHAO Xianzheng,LI Yong,et al.Pa-leochannel restoration and petroleum geological signifi-cance of Neogene in Huanghua Depression [J].Acta Petrolei Sinica,2018,39(2):163-171.
-
[15] 廖保方,薛培华,李列,等.辫状河现代沉积研究与相模式:中国永定河剖析[J].沉积学报,1998,16(1):34-39.LIAO Baofang,XUE Peihua,LI Lie,et al.Study on modern deposit of a braided stream and factes model:taking the Yongding river as an example[J].Acta Sedi-mentologica Sinica,1998,16(1):34-39.
-
[16] MIALL A D.Fluvial depositional systems[M].Switzer-land:Springer International Publishing,2014:25-32.
-
[17] PHILIP J A,JAMES L B,JULIE E R,et al.Morpho-logical evolution and dynamics of a large,sand braid-bar,Jamuna River,Bangladesh [J].Sedimentology,2010,47(3):533-555.
-
[18] SCHUURMAN F,KLEINHANS M G.Bar dynamics and bifurcation evolution in a modelled braided sand-bed riv-er[J].Earth Surface Processes &Landforms,2015,40(10):1318-1333.
-
[19] SCHUURMAN F,MARRA W A,KLEINHANS M G.Physics-based modeling of large braided sand-bed riv-ers:bar pattern formation,dynamics,and sensitivity [J].Journal of Geophysical Research Earth Surface,2013,118(4):2509-2527.
-
[20] 冯文杰,吴胜和,张可,等.曲流河浅水三角洲沉积过程与沉积模式探讨:沉积过程数值模拟与现代沉积分析的启示[J].地质学报,2017,91(9):2047-2064.FENG Wenjie,WU Shenghe,ZHANG Ke,et al.Depo-sitional process and sedimentary model of meandering-river shallow delta:insights from numerical simulation and modern deposition [J].Acta Geological Sinica,2017,91(9):2047-2064.
-
[21] CALDWELL R L,EDMONDS D A.The effects of sedi-ment properties on deltaic processes and morphologies:a numerical modeling study [J].Journal of Geophysical Research Earth Surface,2014,119(5):961-982.
-
[22] HAJEK E A,EDMONDS D A.Is river avulsion style controlled by floodplain morphodynamics?[J].Geolo-gy,2014,42(3):199-202.
-
[23] EDMONDS D A,SLINGERLAND R L.Significant effect of sediment cohesion on delta morphology[J].Nature Geo-science,2010,3(2):105-109.
-
摘要
选取黄骅坳陷新近系典型辫状河沉积解剖区,利用密井网区动静态资料相结合识别和刻画心滩及其内部结构,在此基础上通过沉积数值模拟探究研究区心滩形成过程,并建立研究区心滩顺流迁移和横向迁移演化模式。 结果表明:顺流加积、垂向加积、侧向加积和侵蚀作用共同塑造心滩的形态和内部结构;滩头侵蚀、滩尾顺流加积和垂向加积形成了心滩的顺流生长和垂向生长;侧向加积与侵蚀作用控制心滩的不对称侧向生长过程。
Abstract
A detailed case study is chosen in Huanghua Depression where deposited classic braided rivers in the Neogene. Braided bars and their architecture are characterized by integrated analysis of dynamic and static data in the dense well area. Furthermore, the formation of braided bars is explored by numerical simulation of deposition process,and the downstream mi- gration and lateral migration evolution model of single braided bar are discussed and built. It is concluded that downstream accretion, vertical accretion, lateral accretion and erosion shape the geometry and architecture of braided bars. Downstream and vertical growth of braided bars is controlled by bar head erosion, bar tail accretion and vertical accretion; lateral accre- tion and erosion control the lateral asymmetrical growth of braided bars.
