② 油气资源与探测国家重点实验室, 北京 102249;
③ 中国石油长庆油田分公司勘探开发研究院, 陕西西安 710021
② State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum(Beijing), Beijing 102249, China;
③ Research Institute of Exploration and Development, PetroChina Changqing Oilfield Company, Xi'an, Shaanxi 710021, China
地层结构研究是开展地层等时对比的重要组成部分[1],亦是储层精细表征的前提,因此明确地层结构对于油气田开发具有重要的现实意义[2-4]。针对陆相油田,众多学者已经探索出一套基于井资料的“旋回对比、分级控制、逐级闭合”的地层结构分析与等时对比方法[3, 5-7]。在地震资料品质好的地区,采用“井震结合”开展油层对比的方法也得到了广泛认可。但是,受纵向分辨率的限制,地震反射波形具有一定的多解性[8-9]。比如针对前积反射,众多学者依据同相轴尖灭的特点,选取地层尖灭模式开展地层对比。然而,当地层厚度小于地震分辨率时,地震反射同相轴存在穿时的可能[8, 10-11],地层在未尖灭的情况下亦会出现地震同相轴尖灭现象[12-13]。因此,严格按照地震波形开展地层结构研究,可能会导致认知出现偏差,进而导致地层对比穿时、砂体归属混乱[14]。
地震正演模拟方法是降低地震资料多解性的有效方法之一,也是建立地震反射波形与地层结构对应关系的重要桥梁[11, 15-19]。因此,本文采用地震正演模拟方法,对鄂尔多斯盆地合水地区延长组三段(简称长三段)的地层结构开展研究,明确地层发育特征,为研究区的储层精细表征提供基础。
1 概况合水地区位于鄂尔多斯盆地陕北斜坡西南部,在长三段沉积期,其处于湖盆边缘的缓坡部位。长三段埋深为1200~1950m,发育特低渗致密砂岩储层,从上至下可分为长4+5、长6和长7油组,各油组又可进一步细分成数个砂组[20](表 1)。
地震剖面(图 1)上可见两个反射标志层:第一个标志层为长7和长8油组的分界(即长8_1砂组顶面),为连续的强波峰反射[21];第二个标志层为长6_3砂组顶界面。
研究区二维地震资料揭示在下倾方向可见前积地震反射现象。由于前期对该区地层结构认识不清,地层对比以“平对”为主[22-23],存在穿时现象,导致油气开发效果不理想。另外也有学者采用“下倾尖灭”方式进行地层对比[13],如图 1a所示,长6_3顶地震反射同相轴由z169井向z40井方向逐渐接近长8_1顶地震反射同相轴,并在z40井之前尖灭,由此可得出z40井不发育长6_3砂组(老地震层位解释方案),但实际上z40井钻遇长6_3砂组,说明该对比方式仍存在砂体归属不清的问题。因此,本文利用地震正演模拟方法针对前积地层结构展开研究,以解决地层单元细分、对比问题。
2 前积地层结构及地震波形特征 2.1 前积地层结构样式的概念模型在地震剖面上(图 1)可见,长6_1、长6_3、长7_1及长8_1砂组顶面均为波峰反射。其中,长6_3与长7_1砂组顶面反射同相轴向下倾方向逐渐尖灭在长8_1砂组顶面,形成前积反射结构。以该地震剖面(顺物源方向)前积反射结构为原型,抽象建立三种地层结构样式的概念模型(图 2)。概念模型均由相邻三套地层及其上、下围岩组成,其中1号、2号和3号地层分别对应长7_1顶面与长8_1顶、长6_3顶与长7_1顶、长6_1顶与长6_3顶所限定的地层单元。第一种地层结构样式为等厚尖灭式,1~3号地层的初始厚度(均为80m)相等,且各地层顶面基本平行,1号与2号地层依次下倾变薄、尖灭(图 2a);第二种地层结构样式为底部地层减薄式(图 2b),相邻三套地层初始厚度相等,亦为80m,1号与2号地层在底部快速减薄但未尖灭(厚度减薄至10m),3号地层厚度略有减薄(减薄至60m)。第三种地层结构样式为地层厚度等比例减薄式(图 2c),相邻三套地层初始厚度相等,均为80m,三套地层厚度整体同比例减薄,最终厚度仍近似相等。
研究区三套地层均以反旋回为主(下部以泥岩为主,上部以砂岩为主),自下而上纵波速度与密度值整体呈增大的趋势。纵波速度为3900~4100m/s,平均约为4000m/s。因此,设计概念模型中各套地层的纵波速度由底部3900m/s过渡为顶部4100m/s。为了简化模型,将围岩设置为背景相,其纵波速度设为平均速度,即4000m/s(图 2)。
2.2 概念模型的地震波形响应特征采用基于“波前速度扫描”(Wave Front Sweep Velocity)的褶积方法[24]分别对三种地层结构概念模型(图 2)进行正演模拟。研究区不同二维地震测线主频略有变化,最低约为25Hz,最高约为35Hz。在模拟过程中考虑地震资料品质,分别采用25、30、35Hz的雷克子波进行正演模拟,并与实际地震剖面进行对比、分析,探讨研究区可能的地层结构样式。
2.2.1 等厚尖灭式模型响应特征针对等厚尖灭式地层结构模型,由25Hz雷克子波激发得到的模拟结果表明,地层界面表现为强波峰反射,且反射同相轴能够较准确地反映地层界面位置及地层之间接触关系(图 3a)。伴随着2号、3号地层厚度的不断减薄,地层底部的反射同相轴逐渐与1号地层底部相交,且反射同相轴连接处无错断现象。值得注意的是,受地震子波分辨率的限制,反射同相轴相交位置略微早于该地层界面尖灭位置(图 3a)。
采用30、35Hz雷克子波激发得到的结果与25Hz相似(图 3b、图 3c),但细节上略有不同。采用30Hz雷克子波激发时,3号层底部与1号层底部反射同相轴相交位置连续性差,呈现明显的间断现象;采用35Hz雷克子波激发时,反射同相轴相交位置的连续性较25Hz的模拟结果略差,但无明显间断现象。
2.2.2 底部地层减薄式模型响应特征针对底部地层减薄式地层结构模型,由25Hz雷克子波激发得到的模拟结果(图 4a)表明,地层界面表现为强波峰反射,剖面左侧反射同相轴能够较为准确地反映地层界面的位置。但是由于剖面右侧地层减薄后的厚度(10m)远远小于地震资料可分辨厚度(25Hz条件下,1/4波长约为33.33m),1、2号地层顶部的反射同相轴逐渐与1号地层底部相交,且反射同相轴相交处早于地层厚度减至10m的位置。实际上1、2号地层并未尖灭(地层界面未相交),同相轴反映情况与地层发育情况不符,即当地层厚度小于地震资料可分辨厚度时,地震反射特征无法准确反映地层界面的位置及地层之间接触关系。
采用30Hz与35Hz雷克子波激发得到的结果与前者相似(图 4b、图 4c),且随着子波频率增大,同相轴相交位置越来越靠近地层厚度减至最薄之处。
2.2.3 地层厚度等比例减薄式模型响应特征针对地层厚度等比例减薄式地层结构模型,由25Hz雷克子波激发得到的模拟结果表明,在地层厚度减薄之前,地层界面表现为强波峰反射,同相轴能够较准确反映地层界面位置;地层厚度减至最薄(80/3m,模型右半部分)后,界面的反射振幅明显减弱,且无法准确反映地层界面位置及地层之间接触关系(图 5a)。1号地层减至最薄时,1号地层顶、底界面逐渐合并形成复合地震波形;2号地层与1号地层特征相似,地层减薄后顶、底界面形成复合地震波形;3号地层减至最薄时,其顶部反射同相轴消失,底部反射同相轴与实际地层界面相比向上偏离,地震波形无法准确反映该地层顶、底界面的位置。总之,在地层厚度等比例减薄至80/3m时,地层厚度小于地震可分辨厚度,地震波形为多个地层界面的综合响应,即地震波形无法准确识别界面的位置及地层之间接触关系。
采用30Hz与35Hz雷克子波激发得到的结果与前者相似(图 5b、图 5c)。激发频率由25Hz→30Hz→35Hz变化时,调谐厚度相应变化为40.00m→33.33m→28.57m,地震分辨率随之相应提高,模型右侧多个界面的复合地震波形逐渐呈现分离的趋势,但是受相邻界面波形的干扰,反射同相轴与实际地层界面发生偏离,波形仍难以准确反映地层界面的位置及地层之间接触关系。
与实际地震剖面(图 1)对比可见,长7_1顶、长6_3顶与长8_1砂组顶部反射同相轴依次相交,构成前积反射结构,等厚尖灭式、底部地层减薄式地层结构模型的正演模拟结果与该地震反射特征一致,即2号、3号地层底部反射同相轴依次与1号地层底部相交,组成前积反射结构(图 3、图 4);地层厚度等比例减薄式地层结构模型的正演模拟结果与实际地震剖面反射特征存在明显差异(图 5)。因此,前两种地层结构模型均为本文研究区地震前积反射的可能结构样式。
3 实际应用在恰当的约束条件下,利用测井资料可以建立一个与地下实际情况较为接近的波阻抗模型,其正演模拟结果可以接近实际地震波形。相反,当地层结构样式选取不恰当时,波阻抗模型可能会与实际情况相差甚远,其正演模拟结果与实际波形就会偏差较大。因此,可以依据地震正演模拟结果检验地层结构样式的准确性。
为进一步评价等厚尖灭式、底部地层减薄式地层结构样式在研究区的适用性,本文根据实际资料建立了与以上两种地层结构样式相对应的正演模型(图 6a、图 7a)。在等厚尖灭式正演模型中,长7、长6_3与长6_1顶界面由南(剖面左端)向北(剖面右端)逐渐下倾,长7顶界面尖灭在z169井与z18井之间,长6_3顶界面尖灭于z40井与z15井之间(图 6a);在底部地层减薄式正演模型中,长7、长6_3与长6_1顶界面由南向北逐渐下倾,各地层厚度减薄但不发生尖灭,其中长7油组厚度在z169井处减薄至10~15m,长6_3砂组厚度在z40井处同样减薄至10~15m(图 7a)。为使正演模拟结果更加合理,实际正演模型保留了长4+5_2砂组。基于声波时差与密度(没有密度资料时用Garden公式计算得到[25])得到波阻抗,在相应地层格架的限定下,由波阻抗插值建立波阻抗模型(图 6b、图 7b)。由插值结果可知,在倾斜地层中插值得到的波阻抗模型也具有倾斜特征,且每套地层单独插值,插值结果互不影响。为使模拟结果更加接近地下实际情况,采用研究区实际地震子波(目的层统计子波)激发,并加入了少量噪声,得到结果如图 6c、图 7c所示。
等厚尖灭式模型的正演模拟结果表明,各地层界面以波峰反射为主,但在细节上反射同相轴与地层界面仍存在一定差异(图 6c)。长6_3顶部反射同相轴与长8_1顶部反射同相轴相交,相交处振幅明显减弱,且相交位置略早于地层界面尖灭位置。伴随着地层厚度逐渐减薄,在地层趋于尖灭处地层厚度已小于地震可分辨厚度,地震波形为该地层顶、底界面的复合响应。
底部地层减薄式模型的模拟结果与等厚尖灭式相似,但在细节上存在差异(图 7c)。长6_3砂组在z40井处地层厚度减薄至10~15m,因此,底部反射同相轴自z40井向北表现为长8_1、长7与长6_3顶界面的复合响应,其波形分辨率变低(相位变宽)、振幅减弱。
尽管两种地层结构样式的正演模拟结果略有区别,但整体较为相似,且均与研究区实际地震剖面的前积反射结构特征相吻合,因此这两种地层结构模型均有可能与研究区地层结构特征相符。前人研究及钻井资料表明[13],研究区长7油组底部的张家滩页岩与长6_3砂组顶部的碳质泥岩在全区广泛发育,并未尖灭。因此,底部地层减薄式地层结构样式更符合研究区实际情况,采用该样式可以对本区地层进行准确对比,从而得到合理的地震解释方案(图 1b)。
4 结论本文采用地震正演模拟技术,对鄂尔多斯盆地合水地区延长组三段的地层结构特征进行研究,得出以下结论。
(1) 基于实际地震剖面的前积反射结构特征,在沉积模式的指导下,抽象建立了三种地层结构样式,即等厚尖灭式、底部地层减薄式和地层厚度等比例减薄式。
(2) 等厚尖灭式和底部地层减薄式地层模型的地震正演模拟结果均出现前积反射,表明以上两种模型均可作为前积反射的可能地层结构样式。长7油组在研究区未发生尖灭,因此底部地层减薄式结构样式更符合研究区的地层发育特征。
(3) 针对研究区的地层结构特征的认知可为后续地层等时对比、精细储层表征提供地质基础。基于地震正演的地层结构研究方法可为其他相似研究区提供借鉴。
[1] |
王国鹏, 何光玉. 双河油田厚油层内夹层分布特征[J]. 石油勘探与开发, 1995, 22(2): 55-58. WANG Guopeng, HE Guangyu. Distribution of the intercalations in thick reservoirs, Shuanghe oilfield[J]. Petroleum Exploration and Development, 1995, 22(2): 55-58. |
[2] |
邓宏文, 王红亮, 祝永军, 等. 高分辨率层序地层学——原理及应用[M]. 北京: 地质出版社, 2002.
|
[3] |
岳大力.曲流河储层构型分析与剩余油分布模式研究[D].北京: 中国石油大学(北京), 2006. YUE Dali.The Study on Architecture Analysis and Remaining Oil Distribution Patterns of Meandering River Reservoir[D]. China University of Petroleum(Beijing), Beijing, 2006. http://cdmd.cnki.com.cn/Article/CDMD-11414-2007179266.htm |
[4] |
雷蕾, 韩宏伟, 于景强. 近岸水下扇沉积样式及地震响应特征新认识[J]. 石油地球物理勘探, 2019, 54(5): 1151-1158. LEI Lei, HAN Hongwei, YU Jingqiang. New understandings of near-shore subaqueous fan sedimentary styles and its seismic responses[J]. Oil Geophysical Prospecting, 2019, 54(5): 1151-1158. |
[5] |
吴胜和, 伊振林, 许长福, 等. 新疆克拉玛依油田六中区三叠系克下组冲积扇高频基准面旋回与砂体分布型式研究[J]. 高校地质学报, 2008, 14(2): 157-163. WU Shenghe, YI Zhenlin, XU Changfu, et al. High frequency base-level cycles of alluvial fan and distribution patterns of sandbodies of Lower Karamay Formation (Triassic)in Middle 6th District of Karamayi Oilfield, Xinjiang[J]. Geological Journal of China Universities, 2008, 14(2): 157-163. DOI:10.3969/j.issn.1006-7493.2008.02.004 |
[6] |
赵翰卿. 高分辨率层序地层对比与我国的小层对比[J]. 大庆石油地质与开发, 2005, 24(1): 5-9. ZHAO Hanqing. High-resolution sequential stratigraphy correlation and Chinese subzone correlation[J]. Petroleum Geology & Oilfield Development in Daqing, 2005, 24(1): 5-9. DOI:10.3969/j.issn.1000-3754.2005.01.002 |
[7] |
伊振林, 吴胜和, 杜庆龙, 等. 克拉玛依油田六中区克下组冲积扇地层对比方法探讨[J]. 新疆石油天然气, 2009, 5(4): 1-5. YI Zhenlin, WU Shenghe, DU Qinglong, et al. A discussion of the alluvial fan stratigraphic correlation in Lower Karamay Formation in Liuzhong District Karamay Oilfield[J]. Xinjiang Petroleum Geology, 2009, 5(4): 1-5. DOI:10.3969/j.issn.1673-2677.2009.04.001 |
[8] |
Zeng H L, Charles K. Seismic frequency control on carbonate seismic stratigraphy:A case study of the Kingdom Abo sequence, west Texas[J]. AAPG Bulletin, 2003, 87(2): 273-293. DOI:10.1306/08270201023 |
[9] |
李慧琼.鄂尔多斯盆地延长组前积反射的地质意义[D].陕西西安: 西北大学, 2013. LI Huiqiong.The Interpretation of Progradational Reflection in Lacustrine Yanchang Formation in the Ordos Basin of Central China[D]. Northwest University, Xi'an, Shannxi, 2013. http://cdmd.cnki.com.cn/Article/CDMD-10697-1013253059.htm |
[10] |
邓国成. 利用正演技术识别地震解释中的假象[J]. 地球物理学进展, 2019, 34(1): 64-68. DENG Guocheng. Using forward modeling technique to identify the illusion in seismic interpretation[J]. Progress in Geophysics, 2019, 34(1): 64-68. |
[11] |
周子群. 利用正演模型识别桩西地区在沙河街组地层剥蚀点[J]. 油气地球物理, 2016, 14(4): 39-42. ZHOU Ziqun. Denudation point identification of Es in Zhuangxi area by the application of forward model[J]. Petroleum Geophysics, 2016, 14(4): 39-42. |
[12] |
赵麟天.高升探区沙三、沙四段沉积体系研究[D].黑龙江大庆: 东北石油大学, 2014. ZHAO Lintian.Three, Four Sedimentary System of Gaosheng Exploration Area[D]. Northeast Petroleum University, Daqing, Heilongjiang, 2014. http://cdmd.cnki.com.cn/Article/CDMD-10220-1014384362.htm |
[13] |
李慧琼.陇东地区延长组前积层序及有利储层预测[D].陕西西安: 西北大学, 2017. LI Huiqiong.Progradational Sequence and Reservoir Prediction of Yanchang Formation in Longdong Area[D]. Northwest University, Xi'an, Shannxi, 2017. http://cdmd.cnki.com.cn/Article/CDMD-10697-1018011735.htm |
[14] |
梁宏伟, 穆龙新, 范子菲, 等. 基于正演模拟的曲流河构型层次研究[J]. 中国矿业大学学报, 2014, 43(6): 1063-1069. LIANG Hongwei, MU Longxin, FAN Zifei, et al. The AHP study of architecture pattern of meandering re-servoir on the basis of forward simulation[J]. Journal of China University of Mining, 2014, 43(6): 1063-1069. |
[15] |
岳大力, 李伟, 王军, 等. 基于分频融合地震属性的曲流带预测与点坝识别:以渤海湾盆地埕岛油田馆陶组为例[J]. 古地理学报, 2018, 20(6): 941-950. YUE Dali, LI Wei, WANG Jun, et al. Prediction of meandering belt and point-bar recognition based on spectral-decomposed and fused seismic attributes:A case study of the Guantao Formation, Chengdao Oilfield, Bohai Bay Basin[J]. Journal of Palaeogeography, 2018, 20(6): 941-950. |
[16] |
姜岩, 程顺国, 王元波, 等. 大庆长垣油田断层阴影地震正演模拟及校正方法[J]. 石油地球物理勘探, 2019, 54(2): 320-329. JIANG Yan, CHENG Shunguo, WANG Yuanbo, et al. Seismic forward modeling for correction of fault shadow zones in Changyuan Oilfield, Daqing[J]. Oil Geophysical Prospecting, 2019, 54(2): 320-329. |
[17] |
马乃拜, 金圣林, 杨瑞召, 等. 塔里木盆地顺北地区断溶体地震反射特征与识别[J]. 石油地球物理勘探, 2019, 54(2): 398-403. MA Naibai, JIN Shenglin, YANG Ruizhao, et al. Seismic response characteristics and identification of fault-karst reservoir in Shunbei area, Tarim Basin[J]. Oil Geophysical Prospecting, 2019, 54(2): 398-403. |
[18] |
季玉新, 韩文功, 沈财余, 等. 胜利油田典型地质模型的地震正演[J]. 石油地球物理勘探, 2006, 41(5): 561-572. JI Yuxin, HAN Wengong, SHEN Caiyu, et al. Seismic forward simulation of typical geologic models in Sheng-li oilfield[J]. Oil Geophysical Prospecting, 2006, 41(5): 561-572. DOI:10.3321/j.issn:1000-7210.2006.05.015 |
[19] |
Yue D L, Li W, Wang W R, et al. Fused spectral-decomposition seismic attributes and forward seismic modelling to predict sand bodies in meandering fluvial reservoirs[J]. Marine and Petroleum Geology, 2019, 99: 27-44. DOI:10.1016/j.marpetgeo.2018.09.031 |
[20] |
Li X B, Liu H Q, Pan S X, et al. Subaqueous sandy mass-transport deposits in lacustrine facies of the Upper Triassic Yanchang Formation, Ordos Basin, Central China[J]. Marine and Petroleum Geology, 2018, 97: 66-77. DOI:10.1016/j.marpetgeo.2018.06.019 |
[21] |
康媛.鄂尔多斯盆地长6储层微观孔隙结构与渗流机理差异研究[D].陕西西安: 西北大学, 2016. KANG Yuan.The Study on the Difference of Microscopic Pore Structure and Seepage Mechanism of Chang 6 Reservior in Ordos Basin[D]. Northwest University, Xi'an, Shannxi, 2016. http://cdmd.cnki.com.cn/Article/CDMD-10697-1017027934.htm |
[22] |
曹江骏, 杨友运, 陈朝兵, 等.致密砂岩储层骨架砂体构型特征——以鄂尔多斯盆地合水地区延长组长6段砂体为例[J/OL].沉积学报, 2019, [2019-07-31].https://doi.org/10.14027/j.issn.1000-0550.2019.051. CAO Jiangjun, YANG Youyun, CHEN Chaobing, et al.Configuration characteristics analysis of skeleton sand body with tight sandstone reservior: A case study of the Triassic Chang 6 members in Heshui area, Ordos Basin, NW China[J/OL]. Acta Sedimentologica Sinica, 2019, [2019-07-31]. https://doi.org/10.14027/j.issn.1000-0550.2019.051. |
[23] |
白薷, 张金功, 李渭. 鄂尔多斯盆地合水地区延长组长6_3岩性油藏储层地质建模[J]. 吉林大学学报(地球科学版), 2012, 42(6): 1601-1609. BAI Ru, ZHANG Jingong, LI Wei. Reservoir geologi-cal model of lithological reservoir of Chang 6_3 in Heshui area, Ordos Basin[J]. Journal of Jilin University(Earth Science Edition), 2012, 42(6): 1601-1609. |
[24] |
Hilterman F J. Amplitudes of seismic waves:a quick look[J]. Geophysics, 1975, 40(5): 745-762. DOI:10.1190/1.1440565 |
[25] |
Gardner G H F, Gardner L W, Gregory A R. Formation velocity and density:the diagnostic basics for stratigraphic traps[J]. Geophysics, 1974, 39(6): 770-780. DOI:10.1190/1.1440465 |