2. 西南大学地理科学学院, 重庆 400715;
3. 中国石油化工股份有限公司西南油气分公司, 四川成都 610041
2. School of Geography Sciences, Southwest University, Chongqing 400715, China;
3. Southwest Petroleum Branch Company of SINOPEC, Chengdu 610041, China
利用以地震技术为主的地球物理勘探技术实现地下地层成像是目前石油天然气勘探开发重要技术手段, 一方面可实现地下构造、断层等成像; 另一方面, 可实现不同岩层三维空间形态成像, 目前已衍生了地震地层学、地震沉积学等多个分支学科[1]。就陆相地层中岩性油气藏勘探开发而言, 河流及三角洲相带所沉积的砂岩是主要的储集层, 识别和精细刻画地下埋深数千米的古河道砂岩储层一直是石油天然气勘探开发的重要研究内容。
川西坳陷侏罗系岩性气藏埋深1000~3000m, 主力层系蓬莱镇组(JP)、沙溪庙组(JS)砂岩储层主要受三角洲平原-前缘分流河道控制。三维地震储层预测技术在该区的应用始于上世纪九十年代, 是前期气藏勘探开发关键技术[2-3]。2010年以前, 主要在成都凹陷周边正向构造上建成新场、马井、洛带等气田。2011年以来, 勘探开发主要目标逐渐转向斜坡区及凹陷区河道砂岩气藏[4-5], 气藏具有“深源浅聚、断砂输导、相带控砂、河道控储、物性控产”的总体特征, 复杂致密薄窄河道砂岩储层成为主要开发建产目标, 储层主要特征为:储层致密, 沙溪庙组及蓬莱镇组孔隙度平均值分别为8.6%和9.3%, 渗透率平均值分别为0.21mD(1mD≈0.987×10-3μm2)和0.55mD; 河道宽度一般低于1000m, 多为200~500m, 较早期新场气田沙溪庙组毯状河道砂宽度明显变窄[2]; 河道纵向多期次发育, 横向多河道交织, 其中, 蓬莱镇组纵向划分出29套砂组, 沙溪庙组纵向划分出11套砂组, 各砂组由1~3层小层砂体组成, 砂体厚度薄, 一般低于20m, 多为5~15m。河道相带精细刻画是困扰川西侏罗系气藏勘探开发主要技术瓶颈之一, 岩石物理分析表明, 受储层致密影响, 储层波阻抗高, 与围岩波阻抗差异小, 储层与围岩波阻抗叠置较为严重, 隐蔽性强, 以前采用常规振幅、相干、分频解释等河道相带技术刻画的结果难以满足生产需要[5-10]; 利用岩性反演方法识别河道相带不理想, 如:在已钻井较少的情况下, 目前普遍应用的地质统计学反演方法精度受限[11]。此外, 沉积相带准确识别亦是困扰复杂河道砂岩气藏地质认识与评价的难题, 尤其是钻井分布密度小的斜坡及凹陷区, 地质认识程度低, 地质-地球物理结合不够, 对沉积相带的划分不准确, 制约了储层预测的精度。
针对复杂致密河道砂岩储层预测难点, 形成了“相控找砂、砂中选优”的研究思路和“宏观-微观、定性-定量、逐级约束、逐步逼近”的技术策略, 以地震相带及储层刻画与描述为技术创新研发为重点, 通过持续攻关与勘探开发实践, 形成了波形分类技术、体分频像素成像技术、多河道分期剥离及三维空间刻画等关键技术, 较为系统地解决了河道边界及内幕准确刻画、多河道分期、三维空间多河道描述等难题, 河道边界及内幕刻画效果较传统单一属性明显提升, 同一地震相位内多条河道分期次剥离, 为复杂河道砂岩气藏单河道评价、多河道立体开发发挥了极为重要的支撑作用。在此基础上, 采用将今论古的研究思路, 开展地震成像结果与现今地表典型地貌对比研究, 结合沉积相划分原则及沉积学规律, 综合利用少数钻井资料即可准确识别和划分地下河道沉积相带特征, 有效圈定受优势河道相带控制的优质储层的分布, 提升了有利沉积相带预测精度。所部署的评价及开发井河道砂岩储层钻遇率明显提高, 应用效果较好。
1 河道相带刻画关键技术 1.1 波形分类技术陆相沉积地层地震反射波形信息(振幅、频率、相位及连续性等)与砂泥岩组合样式、岩性、厚度、物性以及含油气性等多种因素相关。地震模型正演及统计分析表明, 储层沉积微相、岩性组合的变化是引起地震反射波形变化的重要原因之一, 这为利用地震波形变化信息预测沉积相带、刻画相带内地层的非均质性提供了理论依据。三维地震波形分类技术是沉积相带预测普遍采用的技术, 尤其对于地层厚度变化较小, 地震反射结构以平行-亚平行为主的地层应用效果较好。该技术采用神经网络方法对地震波形进行分类处理, 将不同的波形以相类形式反映出来, 结合地质、钻井、测井分析, 把代表同一类沉积相的地震反射波分为一类, 以此来代表同一沉积相带或微相, 针对叠置河道发育区, 该技术还可描述河道相带内砂体叠置关系[12-13]。
图 1是成都凹陷什邡气田侏罗系蓬三气藏JP310砂组波形分类图。采用自组织波形分类方法, 识别出北东—南西走向的条带状地震波形带, 在该条带状波形带中部分布有垂直于走向的扇状或舌状地震异常。根据地震相带外形、该波形条带内外钻井河道砂岩储层钻遇情况, 结合河道相带发育特点, 识别出北东—南西向主河道。根据区域沉积背景分析结果, 认为主物源方向为北东—南西向, 中部的河道外缘可能存在天然堤或决口扇沉积, 在河道南部远离物源区的地方可能存在河口坝沉积。同一条河道内部亦存在地震波形特征差异, 表明河道内部砂岩储层特征亦有差异, 这为河道内部开发单元的划分和分类部署提供了重要依据。
分频属性融合技术对沉积相带及储层的刻画效果较常规单一属性分析的效果改善明显, 该分频解释技术对三维地震资料进行特殊处理, 产生具有单一频率的一系列振幅调谐体和离散频率能量体, 人工筛选低、中、高单频振幅属性进行三原色平面融合[9], 但刻画结果差异性较大, 影响预测精度, 且操作不太方便。近年来发展的体分频像素成像技术是频变能量融合处理与体分频处理两项技术的联合运用[12], 针对河道砂岩厚度多变、陆相沉积围岩多变、砂岩与泥岩组合方式多变、物性多变等非均质性特点, 在高保真去噪目标处理基础上, 利用小时窗内不同频段对不同厚度地质异常体边界识别的敏感性, 通过高低频融合实现相带的清晰成像。该方法是在像素去噪处理后的地震数据体上, 开展基于小波变换的体分频处理, 即采用一系列不同的频率和带宽Gabor子波, 对地震道进行分频处理, 得到各个不同频率带的振幅能量体和相位体, 赋予每一单频振幅能量体一种颜色(通常是红、绿、蓝色中的一种), 单频振幅能量体内不同采样点处依据能量差异赋予同种颜色不同的亮度和强度, 最后将各单频地震振幅能量体同一点处不同颜色的亮度和强度按照RGB混合方式融合形成一个频率能量体。对地震数据进行全频段扫描, 用不同频率的调谐振幅变化来研究河道砂体的横向变化特征, 调谐振幅的强弱与亮度对应, 可观察到厚度不同的同期沉积地质体的分布范围, 实现地质体外形的整体刻画, 使相带连续性、相带边界刻画能力提高, 并以此来识别地质体边界、叠置关系及发育期次, 其识别结果较常规地震振幅属性纵横向分辨率和清晰度明显提升。
川西侏罗系三维地震资料主频25~30Hz, 频宽5~80Hz, 相带刻画一般选取25, 35, 45Hz频段对应的调谐振幅分别用红、绿、蓝三色叠合成像。图 2对比了成都凹陷什邡气田侏罗系蓬二气藏JP23砂组振幅属性(图 2a)与体分频像素成像结果(图 2b), 振幅属性图上75井~105井钻遇的北东向河道表征效果不理想, 强振幅的黄色区域隐约可见河道特征, 但河道边缘不清晰, 连续性差, 而体分频像素成像图则非常好地实现了河道边界刻画, 河道连续性亦明显改善。该技术对中江气田沙溪庙组河道外形及边界刻画较常规振幅属性亦取得了更为明显的效果。图 3是沙一气藏JS14砂组河道砂体体分频像素成像刻画结果, 该技术刻画河道边界轮廓、内幕结构等方面能力明显提升, 北部可见多个新月形边滩砂坝连接而成的相对较宽的主河道, 在新月形边滩内部明显可以看到3~4期侧积砂体。此外, 就河道外部形态而言, 南部河道宽度明显窄、河道弯度也明显降低, 与北部存在明显差异, 对后续沉积相带准确识别及划分具有极为重要的支撑作用。
多河道交织是河流三角洲相的普遍特征之一, 而复合相位或单相位内叠置河道分期次刻画是气藏开发必须要解决的技术难点。在借鉴前人单河道识别方法的基础上[14], 基于不同埋深、走向、厚度、物性等河道地质特征差异对地震反射时间、能量及频率的影响, 综合应用地震剖面河道沉积位置比较、三维可视化子体追踪技术、体分频像素成像技术, 通过井震一体化综合分析, 可实现同一相位内叠置河道期次识别与划分, 在中江气田沙溪庙组气藏地震同一复合相位内, 实现了3~4期河道分期剥离。
三维可视化技术、体分频像素成像技术是河流三角洲沉积的多河道分期剥离与三维空间展布刻画最为关键的地球物理技术。在地震剖面小层追踪对比线的控制下, 首先利用三维可视化子体追踪技术自动进行河道砂岩的三维空间识别和追踪, 不同河道赋予不同颜色; 基于沉积时限差异, 还可实现不同埋深河道三维地震体中时间域差异识别。由于不同河道地质特征的差异, 可利用体分频像素成像技术, 通过颜色与亮度的变化表征出各条河道的差异。两种技术相互补充、结果相互印证, 在多期多河道纵横交织的情况下, 结合横切多河道的地震剖面上河道的发育位置差异, 可以有效区分河道走向, 从而划分出河道期次, 实现河道三维空间刻画, 这是常规单一属性分析难以实现的。
图 4为川西中江气田沙一气藏JS11砂组多河道剖面特征及刻画结果, 图中可见在1~2个相位内多条河道纵向叠置、横向频繁迁移的现象。图 4a所示的地震剖面上, 河道所对应的透镜状强反射纵向跳跃性分布, 常规的地震层位追踪难度大, 闭合困难。针对该河道地震反射时窗, 采用三维可视化联合体分频像素成像技术, 解决河道分期次刻画和剥离这一难题。在体分频像素成像处理基础上, 根据不同颜色、宽度、走向等, 锁定不同河道的宏观分布及走向特征(图 4d); 利用三维可视化解释软件, 针对各单河道实施种子点控制与自动追踪, 不同河道用不同颜色表征(图 4b), 从而将原本交叉叠置的多期河道在空间上有效分开。在此基础上开展沉积期次识别, 由于不同河道地震剖面上反射时间不同, 从沉积序次上来说, 时间大的河道沉积时间最早。图 4c可以看出①号河道埋深最大, 其次为②、③号河道, ④、⑤号河道埋深相对较浅。结合平面河道走向识别和三维单河道子体追踪结果、各条河道间的叠置关系及发育时间可以确定河道沉积序次:河道发育时间按编号序次发育, ①号河道沉积最早, ⑤号河道沉积时间最晚。
常规三维地震解释常常是在地震层位拾取基础上, 沿层位开时窗或制作沿层切片, 提取振幅、相干等属性, 识别并刻画河道。对于纵横交织的复杂河道而言, 受层位拾取精度、地震属性横向变化及河道纵向下切深度差异等多因素影响, 河道边界清晰度差异大, 连续性不高, 相带刻画精度受限。经过多年探索, 创立了三维单河道或河道组系三维地震解释方法, 该方法不再针对三维数据开展全三维区层位拾取, 而是依托三维可视化及体分频像素成像技术, 以单河道或河道沉积体系三维空间精细刻画为切入点, 沿河道逐一开展层位拾取, 最终将一个地震相位所对应的同一时期或相近时期沉积的河道在三维空间立体成像, 该方法不仅使复杂河道层位解释精度得到大幅度提升, 还可明显提高三维地震解释效率, 这在大区域三维地震河道识别及刻画工作效率的提升上效果尤为明显。采用单河道或河道组系三维地震解释方法, 对中江气田沙溪庙组气藏、什邡气田蓬莱镇组气藏各砂组进行了全三维空间精细刻画, 以中江气田沙三气藏JS31砂组为例, 该砂组沉积时期发育9个河道沉积体系, 每一个沉积体系, 发育1~4条单河道(图 5)。
河道砂岩气藏优质储层往往发育于河道、河口坝等沉积相带, 将今论古是地质研究的核心思路, 地表地貌及露头剖面结构研究已成为识别沉积相特征极为重要的研究方法。但钻井资料较少的区域, 沉积相带划分精度是有限的, 因此, 基于地球物理河道相带刻画结果的沉积相解释显得尤为重要[15]。
根据区域沉积背景、钻井沉积微相的划分及河道刻画结果, 结合对国内外河流地貌的系统调研, 找到与地下河道特征极为相似的现今河流, 可以准确识别沉积相特征, 编制沉积相图, 从而准确圈定出有利于天然气富集的河道砂岩储层的分布。
利用体分频像素成像技术得到的河道刻画结果以及单河道或河道组系三维地震解释结果, 结合露头和少量钻井岩心、测井等资料可更为直观和准确地进行古地貌分析和沉积相划分。根据陆相河流三角洲沉积亚相划分原则, 河道分叉或多河道的出现往往作为河流相演变为三角洲相的标志, 在三角洲相中, 三角洲平原与前缘亚相的划分依据主要就是古地貌及不同类型河道特有的相带外形差异, 即曲流河是三角洲平原的典型标志, 其前端的顺直河是三角洲前缘的标志。
根据区域沉积相带沉积规律研究成果, 中江气田沙溪庙组河道物源来自北部米仓山地区[3], 以远距离搬运的长轴物源为典型特征, 图 6a中可见沙一气藏JS11砂组沉积时期北部主要为曲流河道沉积, 河道变宽处发育边滩沉积, 南部河道变窄、河道弯度降低, 演变为三角洲前缘。图 6b是在国内外河流地貌的系统调研基础上, 针对中江气田JS11河道特征优选出的相似地貌-加拿大中南部河流萨斯喀彻温河, 萨斯喀彻温河源出落基山东麓, 向东南在马尼托巴省注入温尼伯湖。对比萨斯喀彻温河地貌, 可非常明确地判断出中江气田沙一气藏JS11砂组北部曲流河发育部位上以三角洲平原亚相为主(图 7a), 而中南部地区基本处于三角洲前缘内侧环境, 三角洲前缘入湖部位尚位于中江气田南部。
中江气田沙一气藏JS14砂组沉积时期, 研究区中部发育一条北东—南西走向的复合河道, 为多期曲流河河道及边滩沉积叠加而成(图 7b), 三角洲平原沉积范围较JS11砂组明显扩大, 三角洲平原与三角洲前缘分界线南移到区块中部, 向南部迁移约20km, 预测三角洲入湖部位距离研究区更远(图 7b)。针对中江气田沙溪庙组气藏开发, 在沉积相特征识别基础上, 后续针对不同沉积相带开展砂体构型、储层特征研究及精细描述, 评价落实天然气储量, 编制开发评价及产能建设方案, 并实施钻井部署, 所部署的评价井及开发井河道砂岩钻遇率在98%以上, 成功率较早期明显提高。2012年以来, 先后在成都凹陷什邡气田蓬莱镇组气藏及川西东坡中江气田沙溪庙组气藏新建天然气产能13×108m3/a, 技术应用效果显著。
3 结束语1) 以三维地震资料为基础, 波形分析技术、体分频像素成像技术、多河道分期剥离及三维空间刻画技术等的综合应用, 在刻画不同期次单河道外形及内幕等方面的预测效果均优于常规属性, 刻画效果得到明显提高。
2) 针对多期叠置复杂河道, 以单河道或河道沉积体系为目标的三维地震解释方法较传统三维层位拾取与属性提取方法更为适用, 效果明显, 并可明显提高三维地震资料解释效率。
3) 地震成像结果与现今地表相似地貌对比研究, 是预测和分析沉积相带的有效方法, 可准确识别和划分沉积相单元, 大大提升有利沉积相带预测精度。
[1] |
曾洪流. 地震沉积学在中国:回顾和展望[J]. 沉积学报, 2011, 29(3): 417-426. ZENG H L. Seismic sedimentology in China:A review[J]. Acta Sedimentologica Sinica, 2011, 29(3): 417-426. |
[2] |
吴朝容, 段文燊. 川西坳陷XS地区河道砂体预测研究[J]. 石油天然气学报, 2011, 33(11): 81-84. WU C R, DUAN W S. Prediction of channel sandstone in XS area of western Sichuan Depression[J]. Journal of Oil and Gas Technology, 2011, 33(11): 81-84. DOI:10.3969/j.issn.1000-9752.2011.11.017 |
[3] |
唐建明. 川西坳陷致密非均质气藏储层空间展布刻划[J]. 天然气工业, 2002, 22(4): 19-23. TANG J M. Describing the space distribution of the heterogeneous tight gas reservoirs in west Sichuan Depression[J]. Natural Gas Industry, 2002, 22(4): 19-23. DOI:10.3321/j.issn:1000-0976.2002.04.005 |
[4] |
杨克明, 朱宏权. 川西叠覆型致密砂岩气区地质特征[J]. 石油实验地质, 2013, 35(1): 1-8. YANG K M, ZHU H Q. Geological characteristics of superposed tight sandstone gas-bearing areas in western Sichuan[J]. Petroleum Geology & Experiment, 2013, 35(1): 1-8. |
[5] |
谢刚平, 叶素娟, 田苗. 川西坳陷致密砂岩气藏勘探开发实践新认识[J]. 天然气工业, 2014, 34(1): 6-15. XIE G P, YE S J, TIAN M. New understandings of exploration and development practices in tight sandstones gas reservoirs in Western Sichuan Depression[J]. Natural Gas Industry, 2014, 34(1): 6-15. |
[6] |
刘伟, 尹成, 王敏, 等. 河流相砂泥岩薄互层基本地震属性特征研究[J]. 石油物探, 2014, 53(4): 468-476. LIU W, YIN C, WANG M, et al. Study on the characteristics of basic seismic attributes in fluvial sand-shale interbedded layers[J]. Geophysical Prospecting for Petroleum, 2014, 53(4): 468-476. DOI:10.3969/j.issn.1000-1441.2014.04.013 |
[7] |
陈恭洋, 陈玲, 朱洁琼, 等. 地震属性分析在河流相储层预测中的应用[J]. 西南石油大学学报(自然科学版), 2012, 34(3): 1-8. CHEN G Y, CHEN L, ZHU J Q, et al. Application of seismic attributes analysis fluvial reservoir to prediction[J]. Journal of Southwest Petroleum University(Science & Technology Edition), 2012, 34(3): 1-8. DOI:10.3863/j.issn.1674-5086.2012.03.001 |
[8] |
王世瑞, 王树平, 狄帮让, 等. 基于地震属性特征的河道砂体预测方法[J]. 石油地球物理勘探, 2009, 44(3): 304-313. WANG S R, WANG S P, DI B R, et al. Prediction of channel sand body based on seismic attributes[J]. Oil Geophysical Prospecting, 2009, 44(3): 304-313. DOI:10.3321/j.issn:1000-7210.2009.03.011 |
[9] |
叶泰然, 苏锦义, 刘兴艳. 分频解释技术在川西砂岩储层预测中的应用[J]. 石油物探, 2008, 47(1): 72-76. YE T R, SU J Y, LIU X Y. Application of seismic frequency division interpretation technology in predicting continental sandstone reservoir in the west of Sichuan province[J]. Geophysical Prospecting for Petroleum, 2008, 47(1): 72-76. DOI:10.3969/j.issn.1000-1441.2008.01.012 |
[10] |
毕俊凤, 顾汉明, 刘书会, 等. 河道砂体地震响应特征及影响因素分析——以垦东1地区馆陶组上段河道砂为例[J]. 石油物探, 2013, 52(1): 97-103. BI J F, GU H M, LIU S H, et al. Analysis on seismic response characteristics of channel sands and its influence factors:Case study of Upper Guantao Formation in KD1 area[J]. Geophysical Prospecting for Petroleum, 2013, 52(1): 97-103. DOI:10.3969/j.issn.1000-1441.2013.01.015 |
[11] |
刘俊州, 周小鹰, 王箭波, 等. 地质统计学反演在非均质薄砂岩气藏开发中的实践[J]. 内蒙古石油化工, 2009(7): 16-19. LIU J Z, ZHOU X Y, WANG J B, et al. Practice of opening heterogeneous thin sand rock gases deposit in geometry statistic inversion[J]. Inner Mongolia Petrochemical Industry, 2009(7): 16-19. DOI:10.3969/j.issn.1006-7981.2009.07.006 |
[12] |
武恒志, 叶泰然, 赵迪, 等. 川西坳陷陆相致密气藏河道砂岩储层精细刻画技术及其应用[J]. 石油与天然气地质, 2015, 36(2): 230-239. WU H Z, YE T R, ZHAO D, et al. Fine-scale characterization technique and its application for channel sand reservoir of tight gas clastic reservoirs in west Sichuan[J]. Oil & Gas Geology, 2015, 36(2): 230-239. |
[13] |
徐海, 都小芳, 高君, 等. 基于波形聚类的沉积微相定量解释技术研究——以中东地区X油田为例[J]. 石油物探, 2018, 57(5): 744-755. XU H, DU X F, GAO J, et al. Quantitative interpretation of sedimentary microfacies based on waveform clustering:A case study of X oilfield, Middle East[J]. Geophysical Prospecting for Petroleum, 2018, 57(5): 744-755. DOI:10.3969/j.issn.1000-1441.2018.05.014 |
[14] |
金振奎, 时晓章, 何苗. 单河道砂体的识别方法[J]. 新疆石油地质, 2010, 31(6): 572-575. JIN Z K, SHI X Z, HE M. Identification methods for single-channel sand body[J]. Xinjiang Petroleum Geology, 2010, 31(6): 572-575. |
[15] |
夏竹, 李中超, 贾瑞忠, 等. 井震联合薄储层沉积微相表征实例研究[J]. 石油地球物理勘探, 2016, 51(5): 1002-1011. XIA Z, LI Z C, JIA R Z, et al. Thin reservoir sedimentary microfacies characterization based on well logging and seismic data:a case study[J]. Oil Geophysical Petroleum, 2016, 51(5): 1002-1011. |