迭代法建立的低频模型在储层反演中的应用
王伟, 李炳颖, 黄鑫, 陈易周
中海石油(中国)有限公司上海分公司,上海 200335

第一作者简介: 王伟(1988—),男,工程师,主要从事油气勘探开发中地球物理技术应用研究。Email: wangwei171@cnooc.com.cn

摘要

储层预测的精度直接制约着油气开发的经济性与有效性,而低频信息的补偿是改善储层反演效果的有利途径之一。常规的测井曲线内插低通滤波建模方法难以精细表征复杂地质背景下的低频模型。为降低研究区特殊地质体对储层反演带来的影响,采用迭代法建立低频模型,建模中分步考虑压实作用、特殊岩性等因素影响,并通过反演效果不断迭代更新低频模型,最终准确建立反映不同岩性(包含火山岩)的反演低频模型。研究表明: 该方法能够保证储层反演结果的可靠性,为储层精细描述提供了资料; 同时基于反演结果,对东海西湖凹陷W气田主力层(P1)有利储层进行重新刻画,优化了开发井水平段的位置,该井实施结果证实钻遇了优质储层。用迭代法建立的低频模型为复杂地质背景地区储层预测提供了一种更有效的建模方法。

关键词: 低频信息; 模型; 储层反演; 迭代法; 精细描述
中图分类号:P618.13;P631.4 文献标志码:A 文章编号:2095-8706(2021)01-0125-09
Application of Eeration method to build low-frequency model in reservoir inversion
WANG Wei, LI Bingying, HUANG Xin, CHEN Yizhou
Shanghai Branch of CNOOC Ltd., Shanghai 200335, China
Abstract

The accuracy of reservoir prediction has directly restricted the economy and effectiveness of oil and gas development, and the compensation of low-frequency information is one of the favorable ways to improve the inversion effect of the reservoir. However, the conventional modeling method of well interpolation low-pass filter is difficult to represent the low-frequency models in complex geological background, which seriously restricts the effect of seismic inversion. In order to reduce the impact of the presence of special geological bodies on reservoir inversion in the study area, the authors have adopted the iteration method to build low-frequency model in this paper. The factors such as compaction and special lithology have been considered step by step during modeling, and the low-frequency model was iteratively updated through the inversion effect to accurately establish the inversion low-frequency model for reflecting different lithologies (including volcanic rocks). The results show that this method can reasonably guarantee the reliability of reservoir inversion results and fully provide the data foundation for reservoir detailed description. The favorable reservoir (P1) of the main layer of W gas field in Xihu Sag was also re-characterized based on the inversion results, and the horizontal section of the development well was optimized. Besides, the good reservoirs were confirmed by the implementation results. Therefore, the iteration method provides a more effective modeling method for reservoir prediction in complex geological background.

Keyword: low-frequency information; model; reservoir inversion; iteration method; fine description
0 引言

地震记录频带范围较测井数据窄, 缺少低频和高频信息, 特别是低频段, 由于受到激发、接收条件的限制, 很多地震记录10 Hz以下的信息都是不准确的[1, 2, 3]。然而, 缺失了低频信息, 会造成速度分析质量差、深部地层成像差、反演结果不稳定等不利影响。因此要进行储层精细描述, 必须补全地震低频信息, 而低频模型的建立就成为确定性反演的关键环节之一。低频模型可以通过提高砂体横向分辨率来提高储层定量预测的精度[3, 4, 5, 6]

目前在确定性反演中建立的初始低频模型大都直接用测井曲线数据内插而成, 但一般会存在以下常见3大问题: 一是井资料质量不高, 如井眼垮塌造成的异常段多、井间一致性差等; 二是压实作用影响较大, 特别是目的层在深层, 直接导致模型质量差, 反演结果失真; 三是研究区特殊地质情况影响, 如火成岩、煤层、钙质层等特殊岩性体, 井上直接插值的初始模型不能真实反映特征岩性体的范围, 储层预测出的岩性不准。因此, 建立一个高质量的低频模型必须处理好每个环节, 而针对研究区特殊地质情况影响, 本文采用了迭代法建模, 在考虑压实和特殊岩性的基础上建立初始模型, 然后根据初始反演的结果, 补充特殊岩性并更新低频模型, 最终进行储层反演得到质量较高的储层精细描述资料, 有效指导了气田井位部署及优化。

1 研究区概况

W气田位于东海西湖凹陷中央反转构造带南部[7, 8] (图1), 构造总体表现为被一系列断层复杂化的、呈“ S” 形展布的NNE向大型挤压反转背斜, 平面上可划分为北断背斜、南断背斜和东断块3个局部构造, 其中W气田东断块为深部勘探及浅层调整的重要靶区。东断块火成岩和煤系地层较发育, 导致深部地震资料品质变差, 储层预测难度增加。W气田已钻遇地层自上而下依次为: 新生界第四系全新统— 更新统东海群, 新近系上新统三潭组, 中新统柳浪组、玉泉组、龙井组, 古近系渐新统花港组和始新统平湖组, 其中平湖组未穿。目前已发现的含油气层段主要集中在渐新统花港组和始新统平湖组[9, 10]。本次储层预测评价目的层为深层平湖组, 平湖组细分为P1— P7小层, 均属于低孔低渗甚至特低孔低渗储层, 其地震响应较浅层常规储层更加复杂, 地震属性分析存在难度, 寻找有效的储层预测手段是识别储层及含气性的关键, 而低频建模是研究区储层反演的基础。

图1 研究区地理位置[8]Fig.1 Geological location of the study area[8]

2 迭代法建模

图2为建立反演低频模型的工作流程。首先, 将测井曲线纵波阻抗、纵横波速比(VP/VS)及密度进行去压实处理, 迭代去压实趋势后的测井曲线, 再利用迭代后曲线进行插值建模得到三维体模型, 再将压实趋势补充到该三维体模型上, 得到最终初始低频模型; 其次, 利用该初始低频模型进行第一次反演, 得到纵波阻抗体, 然后基于该反演结果解释特殊异常体, 因特殊岩性体与其他岩性在纵波阻抗上存在明显差异, 可以在三维空间中雕刻出其展布范围, 进而更新原始低频模型, 进行第二次反演。以此类推, 直到建立一个稳定、可靠、准确的带有压实趋势及异常体的反演低频模型。

图2 用迭代法建立低频模型的流程Fig.2 Flow chart of the low-frequency model by iteration method

3 应用及效果
3.1 建立初始模型

因研究区构造特征复杂, 首先需建立符合构造和沉积特征、能基本反映研究区地层层序关系的地质框架模型, 以此来约束反演低频模型。储层沉积稳定地区, 地层继承关系一般选择顶底平行接触模式, 再依据地震解释层位和断层在空间上的接触关系, 搭建一个能反映研究区的地质框架模型(图3), 最后, 在地质格架内对测井数据进行空间插值。

图3 气田地质框架模型Fig.3 Geological framework model of the gas field

在建立反演低频模型过程中, 需考虑本地区强烈的压实趋势和特殊岩性体。地层压实造成断层上下盘相同地层的阻抗值不同, 而且火成岩不均匀分布, 井插值方法有其局限性, 所以在建模时利用井上波阻抗曲线拟合出研究区的压实趋势(图4), 随埋深增大, 波阻抗值逐渐增大。

图4 气田压实趋势Fig.4 Compaction interpretation of the gas field

对测井曲线进行去压实和剔除火成岩等特殊岩性处理后, 基于地质框架模型再次进行井插值并建立模型(图5), 最后补充研究区压实趋势, 建立第一轮反演所需的初始低频模型(图6)。初始低频模型充分考虑了地层压实和特殊岩性的影响。

图5 去压实后井插值模型Fig.5 Well interpolation model after removing compaction

图6 纵波阻抗的初始低频模型Fig.6 Initial low-frequency model of P-impedance

3.2 低频模型迭代优化

基于岩石物理分析[11](图7), 在初始模型得到纵波阻抗反演体上对火成岩异常体进行解释。火成岩的弹性参数特征表现为极高阻抗, 能够与其他岩性很好区分, 所以利用反演得到的纵波阻抗就可以很好地识别火山岩, 取纵波阻抗门槛值13 kg· m-3· m· s-1为限, 阻抗大于该门槛值的数据均为火成岩的反映。基于火成岩的岩石物理特征在三维空间中雕刻出相应高阻抗特征的展布范围(图8), 该套火成岩主要发育在研究区的东南部, 断层上升盘, 纵向上主要发育在平湖组P4层。基于解释结果对原始低频模型进行更新, 进行下一轮反演, 最终迭代出最佳的纵波阻抗低频模型(图9)。

图7 研究区岩石物理量图Fig.7 Physical quantities of rocks in the study area

图8 火成岩连井剖面Fig.8 Connecting-well profile of the igneous rock

图9 纵波阻抗的最终低频模型Fig.9 Final low-frequency model of P-impedance

构建的低频模型主要针对反演绝对波阻抗模型中频率为0~10 Hz部分的信息[12, 13, 14], 在研究中, 利用上述的流程建立纵波阻抗低频模型。另外, 由于叠前反演需要3个反演低频模型: 纵波阻抗、VP/VS和密度[15, 16, 17], 因而利用去压实趋势的方法建立了VP/VS及密度低频模型。图10为最终的反演低频模型(进行10 Hz低通滤波)。

图10 最终低频模型Fig.10 Final low-frequency model

3.3 模型应用效果

图11、图12为利用初始模型反演及最终低频模型叠前反演得到的纵波阻抗联井剖面。利用优化后的低频模型得出的反演剖面中火成岩的阻抗特征表现更加明显, 细节更加丰富, 与井点阻抗的变化趋势更加吻合。图13、图14为利用初始模型反演及最终低频模型叠前反演得到的VP/VS反演剖面, 优化后的低频得出的反演剖面信噪比更高, 分辨率明显提高, 而且与实钻井上岩性更吻合, 低VP/VS值表示砂岩, 与低伽马值对应。

图11 初始低频反演纵波阻抗连井剖面Fig.11 P-impedance inversion for connecting-well profile of the initial low frequency

图12 最终低频反演纵波阻抗连井剖面Fig.12 P-impedance inversion connecting-well profile of the final low frequency

图13 初始低频反演VP/VS连井剖面Fig.13 VP/VS inversion connecting-well profile of the initial low frequency

图14 最终低频反演VP/VS连井剖面Fig.14 VP/VS inversion of connecting-well profile of the final low frequency

根据储层评价结果, 对W10水平井井位进行了优化(图15)。图中背景为P1储层段VP/VS平面图。紫色线为原设计开发井轨迹。靶点1原为P1层设计位置, 蓝色实线为优化后水平井轨迹; 靶点2为优化后P1层位置。从本次反演结果分析, 原设计井位水平段的储层发育较差, 风险较大, 故根据储层预测结果, 结合地质认识, 对原水平井井轨迹进行了优化。W10井实钻表明, 优化后的水平段钻遇了较好的储层, 砂岩钻遇率高达90%以上, 从而降低了开发的风险, 实现了该气田经济有效开发的目的。

图15 W10水平井井位优化图Fig.15 Optimization diagram of W10 horizontal well

4 结论

(1)选择合适的低频建模方法是改善反演质量的关键, 常规的测井曲线内插低通滤波建模方法难以表征复杂地质背景下的低频模型, 利用迭代建模法可以建立准确反映不同岩性(包含火山岩)的反演低频模型, 剔除特殊岩性对反演质量的影响。

(2)对比常规低频建模法和迭代建模法反演的结果得出, 改进后的迭代建模方法明显提高了反演效果, 充分提供了储层精细描述的资料基础, 有效指导了气田井位部署及优化。迭代法建低频模型为复杂地质背景地区储层预测提供了一种新的有效建模方法。

(责任编辑: 常艳)

参考文献
[1] 杨立强. 测井约束地震反演综述[J]. 地球物理学进展, 2003, 18(3): 530-534.
Yang L Q. A review of well log constrained seismic inversion[J]. Prog Geophys, 2003, 18(3): 530-534. [本文引用:1]
[2] 周新桂, 杜治利. 北方新区新层系油气地质调查与勘探进展及成果[J]. 中国地质调查, 2019, 6(4): 1-10.
Zhou X G, Du Z L. Progress and achievements in oil and gas geological survey and exploration of new strata in new area of Northern China[J]. Geol Surv China, 2019, 6(4): 1-10. [本文引用:1]
[3] 刘炳杨. 低频建模技术研究及其在AVO反演中的应用[D]. 青岛: 中国石油大学(华东), 2013: 4-10.
Liu B Y. Research of Low Frequency Modeling Method and Its Application in AVO Inversion[D]. Qingdao: China University of Petroleum (East China), 2013: 4-10. [本文引用:2]
[4] 叶云飞, 崔维, 张益明, . 低频模型对波阻抗反演结果定量解释的影响[J]. 中国海上油气, 2014, 26(6): 32-36.
Ye Y F, Cui W, Zhang Y M, et al. Impacts of low-frequency models on the quantitative interpretation of acoustic impedance inversion[J]. China Offshore Oil Gas, 2014, 26(6): 32-36. [本文引用:1]
[5] 程慧. 基于地震资料的地质建模及效果评价[D]. 青岛: 中国石油大学(华东), 2009: 2-10.
Cheng H. Geological Modeling Based on Seismic Data and Its Effect Evaluation[D]. Qingdao: China University of Petroleum (East China), 2009: 2-10. [本文引用:1]
[6] Chen S Q, Wang S X, Zhang Y G, et al. Reservoir prediction using pre-stack inverted elastic parameters[J]. Appl Geophys, 2009, 6(4): 349-358. [本文引用:1]
[7] 张国华, 张建培. 东海陆架盆地构造反转特征及成因机制探讨[J]. 地学前缘, 2015, 22(1): 260-270.
Zhang G H, Zhang J P. A discussion on the tectonic inversion and its genetic mechanism in the East China Sea Shelf Basin[J]. Earth Sci Front, 2015, 22(1): 260-270. [本文引用:1]
[8] 张雷, 姜勇, 侯志强, . 西湖凹陷低孔渗储层岩石物理特征分析及叠前同步反演地震预测[J]. 中国海上油气, 2013, 25(2): 36-39.
Zhang L, Jiang Y, Hou Z Q, et al. Petrophysical analysis and simultaneous pre-stack inversion prediction of the low porous and permeable reservoirs in Xihu sag[J]. China Offshore Oil Gas, 2013, 25(2): 36-39. [本文引用:1]
[9] 张建培, 徐发, 钟韬, . 东海陆架盆地西湖凹陷平湖组-花港组层序地层模式及沉积演化[J]. 海洋地质与第四纪地质, 2012, 32(1): 35-41.
Zhang J P, Xu F, Zhong T, et al. Sequence stratigraphic models and sedimentary evolution of Pinghu and Huagang Formations in Xihu trough[J]. Mar Geol Quat Geol, 2012, 32(1): 35-41. [本文引用:1]
[10] 陈忠云, 张建培, 张涛, . 西湖凹陷层序划分及海平面变化响应[J]. 海洋地质前沿, 2013, 29(9): 15-20.
Chen Z Y, Zhang J P, Zhang T, et al. Stratigraphic sequences in Xihu sag and their reponse to sea level changes[J]. Mar Geol Front, 2013, 29(9): 15-20. [本文引用:1]
[11] 马劲风, 王学军, 谢言光, . 波阻抗反演中低频分量构建的经验与技巧[J]. 石油物探, 2000, 39(1): 27-34, 41.
Ma J F, Wang X J, Xie Y G, et al. Experience and skill of constructing low-frequency components in impedance inversion[J]. Geophys Prospect Pet, 2000, 39(1): 27-34, 41. [本文引用:1]
[12] 陈妍, 张明振, 罗霞, . 山东埕岛东部东营组三角洲砂体地震预测技术[J]. 中国地质调查, 2016, 3(2): 51-57.
Chen Y, Zhang M Z, Luo X, et al. Application of seismic technology to predicting the delta sand body of Dongying Formation in Chengdao Eastern Slope of Shand ong province[J]. Geol Surv China, 2016, 3(2): 51-57. [本文引用:1]
[13] 李爱山, 印兴耀, 张繁昌, . 叠前AVA多参数同步反演技术在含气储层预测中的应用[J]. 石油物探, 2007, 46(1): 64-68, 93.
Li A S, Yin X Y, Zhang F C, et al. Application of prestack simultaneous AVA multi-parameter inversion in gas-bearing reservoir prediction[J]. Geophys Prospect Pet, 2007, 46(1): 64-68, 93. [本文引用:1]
[14] Zoeppritz K. On the reflection and penetration of seismic waves through unstable layers[J]. Goett Nachr, 1919, 1(7B): 66-84. [本文引用:1]
[15] 李炳颖, 王伟, 黄鑫, . 叠前弹性参数反演及FFP分析技术在东海A气田开发中的应用[J]. 中国地质调查, 2020, 7(2): 101-108.
Li B Y, Wang W, Huang X, et al. Application of pre-stack elastic parameters inversion and FFP analysis techniques in the exploitation of A gas field, East China Sea[J]. Geol Surv China, 2020, 7(2): 101-108. [本文引用:1]
[16] 姚云霞, 刘江, 张雷. 叠前纵横波联合反演技术在PH油气田薄储层预测中的应用[J]. 中国海上油气, 2011, 23(6): 377-379, 386.
Yao Y X, Liu J, Zhang L. An application of pre-stack simultaneous inversion of P-and S-wave to thin reservoir prediction in PH field[J]. China Offshore Oil Gas, 2011, 23(6): 377-379, 386. [本文引用:1]
[17] 强敏, 周义军, 钟艳, . 基于部分叠加数据的叠前同时反演技术的应用[J]. 石油地球物理勘探, 2010, 45(6): 895-898.
Qiang M, Zhou Y J, Zhong Y, et al. Application of pre-stack simultaneous inversion based on partial stack data[J]. Oil Geophys Prospect, 2010, 45(6): 895-898. [本文引用:1]