1、莫西庄地区三维工区非线性反演与砂体识别 李琼 何建军 李正文 油气藏地质及开发工程国家重点实验室(成都理工大学) 摘 要： 准噶尔盆地中部莫西庄的侏罗系三工河组储层砂体在纵向、横向上变化大, 且物性不均匀。本文通过对钻井、测井、试油等地质资料的分析研究, 以三维地震资料为基础, 采用地震高分辨率非线性三维整体反演方法, 对多种参数进行了反演和处理, 获得了反映砂岩分布特征的速度、密度及伽马等属性参数的三维数据体, 并综合地质资料进行沉积分析、砂岩累计厚度和有效砂岩累计厚度及储层分布研究。通过物性特征分析, 建立了识别和评价三工河组砂体的物性标准, 利用岩性划分标准和地震多属性数据, 提取出三工
2、河组砂体和有利砂体的空间展布特征。三工河组砂岩和有利砂岩发育, 并以条带、砂坝和席状砂形式展布, 展布方向与区域沉积方向一致。其中三工河组优质储层段 J1s21 和J1s22 厚度较大, 具备良好的物性条件, 具有较高的产储能力和勘探开发前景。关键词： 莫西庄; 侏罗系; 储层特征; 物性参数; 地震属性; 地震沉积; 作者简介：李琼 (1968-) , 女, 博士, 教授, 研究方向:油气地球物理勘探及岩石物理, E-mail:。收稿日期：2017-04-26基金：国家自然科学基金项目 (41274129) Identification of sand body based on seism
3、ic high resolution 3D nonlinear inversion in Moxizhuang, Junggar Basin, ChinaLI Qiong HE Jianjun LI Zhengwen State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Chengdu University of Technology; Abstract： The sandstone reservoir of Jurassic Sangonghe Formation at Moxizhuang area
4、in JunggarBasin is controlled by lithology and structures, and characterized by inhomogeneous physicalconditions caused by superimposition of multiphase channel sand body, the vertical and horizontalvariation and the ambiguous distribution of sand bodies. Based on 3D seismic data, such as welldrilli
5、ng, logging and oil testing, a lot of parameters, including velocity, density, and gammareflecting distribution features of sandstone are obtained by inversion and processing with highresolution nonlinear 3D inversion method. These data are integrated with other geological data toanalyze the deposit
6、ion, grand total thickness of sandstone, effective total thickness of sandstone anddistribution of reservoir. The identification and evaluation standard of sandbody of SangongheFormation is established on the basis of characteristic analysis of physical property. Thecharacteristics of spatial distri
7、bution of sand body and effective sand body of the Sangonghe Formation are obtained according to lithological classification standard and seismic multi-attributedata. It reveals that the sand body and effective sand body of the Sangonghe Formation occur asstripe, sand bar and sheet-like sand paralle
8、l to the regional deposition direction. The sand bodies inthe J1s 21 and J1s22 Sections of Sangonghe Formation have relatively large thickness, and therefore, bear favorable physical property, high capacity of oil and gas and good prospect for exploration and development of natural resources.Keyword
9、： Moxizhuang; Jurassic; reservoirs characteristics; physical property parameter; seismic attributes; seismic sedimentation; Received： 2017-04-26莫西庄地区位于准噶尔盆地中部区块, 构造上属于盆地腹部的中央凹陷带。盆地腹部油气资源丰富, 具备形成大型地层岩性复合油气藏的地质条件。2002 年, 中部区块庄 1 井在侏罗系三工河组 (J 1s) 钻遇油流, 日产油 25.2 t, 随后又在2 个低幅度构造侏罗系分别取得油气突破, 展示了良好的勘探前景1-2
10、。现有的地震资料在该目的层的频带范围在 875 Hz, 主频为 2535 Hz, 目的层砂组厚度 1030 m, 单砂层厚度平均在 35 m, 10 m 的砂层很少, 地震资料在纵向上所能分辨的厚度25m (层速度选取 4 km/s) , 目的层砂体厚度远低于地震 /4 分辨率。所以在落实厚度小于分辨率的薄砂体还存在着一定问题, 用常规地震手段识别储层非常困难。针对砂体识别, 主要采用地震属性分析方法和直接反演方法。利用相干体、地震道波形分类、谱分解技术等对河道砂体、薄层砂体进行预测是行之有效的3-4。利用反演方法和三维可视化技术可以有效地进行河道砂体追踪以及储层横向预测5-6。而联合地震反演
11、和地震属性对储层砂体进行识别也是一种行之有效的方法7。地震非线性反演已成为一种比较有效的岩性反演方法, 它具有很强的稳定性、预测性, 特别是在复杂构造砂泥岩地区, 薄层和岩性尖灭体预测时, 能清晰反映砂体在纵横方向上的变化特征8-11。利用非线性反演方法, 结合测井、录井、测试资料, 可以较为准确地圈定有利储层发育区及储层的有效厚度, 反演的中低频分量物性剖面反映了大套地层的物性特征, 而高频分量反映局部特征, 可直接用于刻画有效砂体的展布特征12-13。基于目前的勘探现状, 利用先进的非线性方法技术, 开展储层的精细分析, 寻找有利的勘探区域, 已成为当务之急。主要目的层三工河组埋藏深度大,
12、 由多期河道砂叠置组成, 叠置的砂体间由于旋回性的变化造成物性条件的不均一性, 砂体的纵横方向上展布特征不清楚并且变化大。为此, 开展了物性特征分析, 建立砂岩和有利砂岩划分标准, 在岩性划分标准和三维地震属性体分析基础上, 对三工河组地震沉积相和砂体空间展布进行追踪和刻画, 对砂体进行识别和评价, 圈定有利砂岩的分布。结合物性参数、沉积相、砂体厚度和含油气分布面积等, 对三工河组产储能力作出预测和评价, 图 1 为莫西庄三维工区位置示意图。图 1 莫西庄三维工区位置示意图 Fig.1 Map showing three-dimensional study area of Moxizhuang
13、 下载原图1 储层地质特征三工河组砂岩分布于准噶尔盆地中央凹陷带的盆地腹部地区, 其构造为宽缓的古隆起, 并在隆起构造的背景上发育了一些局部低幅度背斜和鼻状构造。莫西庄地区位于马桥凸起和中拐凸起向盆地南倾没处。图 2 是 J1s2顶构造图, 由图2 可以看出, 三工河组砂岩储层表现为西南方向倾伏的鼻状构造, 正断层发育, 断层走向大多为近东西向或北东向。图 2 莫西庄地区 J1s2 顶构造图 Fig.2 The structural map illustrating the top surface of J1s2at Moxizhuang area 下载原图三工河组自上而下分为 3 段:第一段
14、 (J 1s1) 以灰、深色泥岩为主。第二段 (J1s2) 为一套河流相为主体的沉积物, 可分为 2 套砂岩组, 上部砂岩厚度较小, 为 816 m, 为河流相砂体沉积;下部砂岩较厚, 厚度为 3065 m, 以河流相砂体沉积为主的中细砂岩。第三段 (J 1s3) 为一套砂泥岩互层, 中部泥岩较为发育, 上部与下部则发育砂岩。2 储层物性特征2.1 地震反射特征图 3 是研究区目的层地震剖面。由图 3 可看出, J 1s1、J 1s2层多具平行反射组合, J1s2反射层多具前积结构, 反映砂体有叠置特征 (图 3 中的蓝色虚线所示) 。2.2 地震高分辨率非线性储层反演方法原理将遗传算法、AN
15、FIS 神经网络技术、禁忌搜索算法及混沌算法的优势相结合形成一种自适应优化方法 (AGACS) , 这种反演方法是与测井资料结合, 对 ANFIS网络进行混沌化, 接着对目标问题进行编码、选择、交叉和变异等操作, 在交叉点上加入禁忌搜索算法 (TS) , 使染色体不断进化, 从而求得全局最优解。图 3 研究区目的层地震反射特征 Fig.3 The characteristics of seismic reflection of objective layer in the research area 下载原图该算法可以对三维地震数据进行稳定及可靠的反演, 在层位控制下, 将工区多井的测井数据与
16、井旁地震道数据输入 ANFIS 网络, 同时进行整体优化训练, 得到整个三维工区的自适应权函数, 并建立非线性映射关系, 再根据储层在纵横方向上的地质变化特征更新这种非线性映射关系, 获得分辨率高且稳定的地震反演数据体 (速度数据体、波阻抗数据体、密度数据体等) , 可用于三工河组砂岩储层的定量分析。图 4 为储层高分辨率非线性三维反演流程图14。2.3 地震反演属性特征利用高分辨率反演方法, 获得了速度、密度和伽马属性数据体, 图 5 是地震反演属性剖面图。由图 5 可以看出, J 1s2层中, 底部以速度为 4.254.4km/s、密度 2.32.5 g/cm 和伽马值2.40 g/cm,
17、 伽马为 4090 API;对于 J1s2砂岩, 速度为4.34.4 km/s、密度为 2.42.5 g/cm, 伽马为 4090 API;对于 J1s2砂岩, 速度为 4.34.6km/s、密度为 2.42.5 g/cm, 伽马为 4080 API。对于有利砂岩, 速度为 4.34.6 m/s, 密度为 2.452.50 g/cm, 伽马值为 4070 API。通过对井地质与测井资料对比分析, 结合前期研究成果, 确定有利砂岩速度为4.34.6 km/s, 密度为 2.452.50 g/cm, 伽马值为 4070 API。3 储层地震沉积特征在地震沉积学分析中, 将地震反演提取的速度、密度和
18、伽马等参数进行融合处理, 形成表征沉积特征的综合属性参数, 进而获得沉积相图 (图 6) 。由图 6-A可看出, 物源自 NE 进入, 分为 2 支, 一支为 NE-SE 向延伸, 河道有小支流;另一支河道由 NE 进入后, 近 EW 向延伸, 然后转为 NE-SW 向延伸, 转向过程有分叉后合并现象。由图 6-B 可看出, 物源自 NE 向进入, 分为 2 支, 一支在庄 4 井南, NE-SW 向延伸;另一支河道在庄 106-庄 107 井一线, 特征明显, 由 NE 进入后呈 E-W 向延伸, 然后转为 NE-SW 向延伸。总体上河道分支较小, 曲流河特征明显, 在主河道边缘地带砂体发育
19、, 且处于砂泥交叉地带, 油气容易被遮挡形成岩性油气藏。由图 6-C 可看出, 物源自 NE 进入, 分为 4 支河道, 第 1 支值较低, 形态可见;第 2第 4 支主河道特征清楚;第 2、第 3 支主河道从 NE 进入后近E-W 向延伸, 在西部折向 NE-SW 向, 2 条主河道在庄 102 井和庄 5 井交叉;第 4支主河道是主砂带, 从 NE 进入 SW 过程中往南分化出多支相互交叉的小河道, 形成网状交织, 辫状河三角洲特征明显。表 1 物性参数统计表 Table 1 Statistics of property parameters 下载原表 图 5 地震反演属性剖面图 Fig.
20、5 The profile of seismic inversion 下载原图图 6 地震沉积相分布图 Fig.6 The distribution of seismic sedimentary facies 下载原图4 砂体与有利砂体空间展布特征通过钻井与测井对砂体的物性参数进行对比分析, 确定了砂体和有利砂体的速度、密度和伽马值 (表 1) 。根据砂体和有利砂体的速度、密度和伽马值计算, 提取了各层的砂体和有利砂体的厚度值, 并利用钻井和测试确定的厚度进行对比校正, 最终获得了储层的砂体和有利砂体空间展布特征。4.1 砂体空间展布特征由图 7-A 可看出, 砂岩厚度在 040 m 之间变化
21、, 最大厚度在 40 m 左右。砂体的分布分为 3 个条带:北部条带 NEE-SW 向延伸, 中部条带 NEE-SW 向展布, 南部条带 NEE-SW 向展布, 这 3 个条带相互交接, 与区域沉积方向一致。由图 7-B 可看出, J 1s2砂岩厚度在 1040 m 之间变化;最大厚度区域在庄 4 井南, 呈 NE-SW 向的条带, 厚度为 3040 m, 是一个大型砂坝。在工区西部和南部有厚度为 1025 m 的砂岩分布, 为远砂坝和席状砂分布区。由图 7-C 可看出, 砂岩厚度为 10110 m, 从东向西、由北向南减薄。厚度40 m 的区域在工区中北部呈片状分布, 仅在南部局部地区厚度6
22、0 m 的区域基本上呈条带状沿 NE-SW 向延伸, 为主河道区, 砂岩分布带延伸过程中产生多条分支。4.2 有利砂体空间展布特征由图 8-A 可看出, J 1s1层有 4 个分布区域: (1) 庄 4-庄 107-庄 102-庄 5 一线分布区, 呈 EW 向的条带分布, 厚度为 08 m, 是主要有利砂岩分布区, 是砂坝集中区; (2) 庄 4 井以北, 厚度为 45.5m, 为点砂坝分布区; (3) 庄 103 井西, 厚度为 37m, 为远砂坝区; (4) 庄 108 井南, 块状分布, 厚度为 25 m, 为远砂坝。由于 J1s1有利砂岩厚度不大, 并且分布不均, 多为含水层, 储层
23、品质较差。由图 8-B 可看出, J 1s2层有利砂岩厚度分布趋势与砂岩厚度分布趋势相比基本一致, J 1s2有效砂岩厚度为 015 m。由图 8-C 可看出, J 1s2层有利砂岩厚度分布趋势与区域砂岩厚度特征基本一致, 有利砂岩厚度仅为 022 m。在工区的东北角基本呈片状, 厚度较大, 达 1022 m, 在西部由 NE 向 SW 呈条带延伸;有利砂体呈 4 条带分布, 且基本平行延伸。5 油气预测三工河组储层的地质特征复杂, 属于河流-湖泊沉积体系, 岩性以砂泥岩互层为主, 岩性纵横向变化快, 储层为水下河道砂体。受岩性和构造共同制约, 为较典型的构造-岩性油气藏。三工河组储层的含油气
24、性与砂岩体空间展布、有利砂岩空间展布及物性特征等密切相关。5.1 地震沉积特征与储层油气预测J1s1层的沉积特征为水下河道、砂坝和远砂坝, 岩性多为粉砂岩-细砂岩;J 1s2层沉积特征为水下分流河道 (曲流河) 和席状砂, 砂体发育, 砂泥交替或交叉, 油气易被遮挡形成油气藏;J 1s2层有 4 条河道, 第 4 条河道是主砂带, 并分化多条交叉的小河道, 形成网状交织的辫状河三角洲, 是油气富集环境和富集场所。5.2 物性参数特征与储层油气预测表 1 可知, 三工河组储层的波速和密度高、伽马值低, 代表砂体质地更纯, 对油气聚集更有利。表 2 是物性参数与含油气统计表, 由表 2 可以看出,
25、 在上亚段 J1s2中, 速度为 4.254.4 km/s、密度为 2.32.5 g/cm、低伽马区 (85.0API) , 且有一定厚度和横向延伸度时, 产油气潜力大;在下亚段 J1s2中, 速度为 4.34.5 km/s、密度为 2.42.5g/cm、低伽马值区, 则产油可能性大。5.3 有利砂岩储集能力预测由有利砂岩空间展布特征可知, J 1s1层有利砂岩厚度不大, 分布不均, 储集能力较差, 无开采价值。J 1s2层有利砂岩的有效厚度在 015 m 之间, 有效厚度 4 m 是储层下限, 有效厚度较大, 储集能力强, 产能也高。J 1s2层有利砂岩的有效厚度在 022.0m 之间, 有
26、效厚度较大, 连片性好, 储集能力强, 产能也高。表3 是产能与有利砂岩厚度和含油分布面积统计表。由表 3 可看出, 产层累计有效厚度大的层段, 产能高, 因此, 有利砂岩有效厚度与产能密切相关。图 7 砂岩厚度分布图 Fig.7 The distribution of sandstone thickness 下载原图图 8 有利砂岩厚度分布图 Fig.8 The thickness distribution of favorable sandstone 下载原图表 2 物性参数与含油气统计 Table 2 The statistics of physical character parame
27、ter and prediction of oil and gas 下载原表 表 3 产能与厚度和面积统计 Table 3 The statistics of capacity, thickness and area 下载原表 应用研究不连续现象的突变理论与技术, 圈定了庄 103、庄 106、庄 107 和庄 5井区的含油面积 (表 3) , 含油面积表明了储集能力, 4 个井区含油面积为11.33 km, 庄 106 区块储集能力最大, 其次为庄 103 井区, 最小为庄 5 区块。J1s2和 J1s2层段是三工河组的优质储层段, 具有较高的产储能力和勘探开发前景。6 结论通过物性分析,
28、确定了三工河组储层各段砂岩和有利砂岩的物性判别依据, 并以此岩性划分依据采用地震高分辨率非线性三维整体反演方法对准噶尔盆地中部莫西庄三维工区提取各地震属性参数, 用多属性的砂泥比参数融合成综合参数, 划分出三工河组储层各段的沉积相。J 1s2以 2 条河道为主;J 1s2物源自 NE进入, 以 4 条河道为主。依据砂岩和有利砂岩划分标准, 利用多属性参数, 提取各层段的砂岩和有利砂岩的厚度。砂岩和有利砂岩发育, 砂岩和有利砂岩以条带状、砂坝和席状形式展布, 并与区域沉积方向一致。三工河组储层为构造-岩性储层。河流-湖泊沉积体系和水下河道砂体是油气富集环境和场所;三工河组各层段的速度、密度高和低
29、伽马值的纯砂岩有利于油气储集;J 1s2和 J1s2层段的有利砂岩厚度较大和有一定的含油气面积, 具有较高的产储能力, 为优质储层, 具有勘探开发前景。参考文献1贾发敬.庄 1-庄 101 井区油藏类型初探J.内蒙古石油化工, 2004, 30 (5) :106-108.Jia F J.A preliminary study on the types of reservoirs in Zhuang-1-Zhuang-101 WellsJ.Inner Mongolian Petrochemical Industry, 2004, 30 (5) :106-108. (in Chinese) 2陈晓
30、春.新疆准中地区老井试油复查可行性分析研究J.石油天然气学报, 2009, 31 (2) :309-311.Chen X C.Feasibility study on formation testing of old well oil-test in Xinjiang Zhunzhong areaJ.Journal of Oil and Gas Technology, 2009, 31 (2) :309-311. (in Chinese) 3刘伟, 贺振华, 陈学华.地震多属性在浊积扇河道砂体识别中的应用J.西安石油大学学报 (自然科学版) , 2011, 26 (5) :6-10.Liu W
31、, He Z H, Chen X H.Application of seismic multiattributes in the identification of the channel sand in turbidite bodyJ.Journal of Xian Shiyou University (Natural Science Edition) , 2011, 26 (5) :6-10. (in Chinese) 4田鑫, 王绪本, 张铭, 等.地震属性方法在油田开发阶段薄砂体识别中的应用 以印尼苏门答腊盆地 Gemah 油田 M 油层为例J.天然气地球科学, 2011, 22 (3
32、) :533-538.Tian X, Wang X B, Zhang M, et al.Application of seismic attributes identifying thin sand body:A example of M Formation in Gemah oilfield of Sumatrina BasinJ.Natural Gas Geoscience, 2011, 22 (3) :533-538. (in Chinese) 5马沛东.准噶尔盆地车排子地区上第三系砂体识别方法J.石油物探, 2006, 45 (5) :487-491.Ma P D.Identifica
33、tion method of Upper Third Series sand bodies in Chepaizi of Junggar BasinJ.Geophysical Prospecting for Petroleum, 2006, 45 (5) :487-491. (in Chinese) 6刘新辉, 阳建平, 袁蔚, 等.柯克亚气田西八段砂体及有利区地震预测J.天然气工业, 2012, 32 (10) :31-33.Liu X H, Yang J P, Yuan W, et al.Seismic prediction of sand bodies and plays in the
34、8thmember of the western Xihefu Formation in the Kekeya Gas Field, Southwest Tarim BasinJ.Natural Gas Industry, 2012, 32 (10) :31-33. (in Chinese) 7王晓龙, 石万忠, 何生, 等.万昌地区永二段水下滑塌扇内上倾尖灭砂体识别J.地球科学, 2013, 38 (3) :561-572.Wang X L, Shi W Z, He S, et al.Identifying updip pinchout sandstone in subaqueous slu
35、mp fans of Yonger Member using acoustic impedance and instantaneous phase in Wanchang area, Yitong BasinJ.Earth Science, 2013, 38 (3) :561-572. (in Chinese) 8吴建军, 杨培杰, 王长江.地震多属性非线性反演方法在东营三角洲中的应用J.油气地质与采收率, 2013, 20 (1) :52-54.Wu J J, Yang P J, Wang C J.Application of seismic multi-attributes nonline
36、ar inversion in Dongying deltaJ.Petroleum Geology and Recovery Efficiency, 2013, 20 (1) :52-54. (in Chinese) 9黄捍东, 张如伟, 魏世平.地震非线性随机反演方法在陆相薄砂岩储层预测中的应用J.石油学报, 2009, 30 (3) :386-390.Huang H D, Zhang R W, Wei S P.Research on application of seismic nonlinear random inversion to reservoir prediction in th
37、e thin sandstone of continental depositsJ.Acta Petrolei Sinica, 2009, 30 (3) :386-390. (in Chinese) 10李琼, 贺振华.地震高分辨率非线性反演在薄互储层识别中的应用J.成都理工大学学报 (自然科学版) , 2004, 31 (6) :708-712.Li Q, He Z H.Application of seismic high resolution nonlinear inversion to thin and interbedded reservoir recognitionJ.Journa
38、l of Chengdu University of Technology (Science&Technology Edition) , 2004, 31 (6) :708-712. (in Chinese) 11李琼, 李勇, 李正文, 等.基于 GA-BP 理论的储层视裂缝密度地震非线性反演方法J.地球物理学进展, 2006, 21 (2) :465-471.Li Q, Li Y, Li Z W, et al.A seismic nonlinear inversion for apparent fracture density of hydrocarbon reservoir based
39、on GA-BP theoryJ.Progress in Geophysics, 2006, 21 (2) :465-471. (in Chinese) 12张志伟, 王春生, 林雅平, 等.地震相控非线性随机反演在阿姆河盆地 A 区块碳酸盐岩储层预测中的应用J.石油地球物理勘探, 2011, 46 (2) :304-310.Zhang Z W, Wang C S, Lin Y P, et al.Applications of facies-controlled nonlinear random inversion for carbonate reservoir prediction in B
40、lock A of Amu-darya BasinJ.Oil Geophysical Prospecting, 2011, 46 (2) :304-310. (in Chinese) 13吴媚, 符力耘, 李维新.高分辨率非线性储层物性参数反演方法和应用J.地球物理学报, 2008, 51 (2) :546-557.Wu M, Fu L Y, Li W X.A high-resolution nonlinear inversion method of reservoir parameters and its application to oil gas explorationJ.Chinese Journal of Geophysics, 2008, 51 (2) :546-557. (in Chinese) 14李勇, 雷旭友, 李正文.整体反演方法研究及软件实现C/中国地球物理2009.合肥:中国科学技术大学出版社, 2009:113.Li Y, Lei X Y, Li Z W.Integrated inversion method and its implementationC/The Chinese
