2. 中国石油化工股份有限公司西北油田分公司, 新疆乌鲁木齐 830011;
3. 中国石油化工股份有限公司石油物探技术研究院, 江苏南京 211103
2. Sinopec Northwest Company, Urumqi 830011, China;
3. Sinopec Geophysical Research Institute, Nanjing 211103, China
碳酸盐岩缝洞型油气藏具有埋藏深、非均质性强等特点[1-3], 溶洞体在地震成像剖面上表现为“串珠状”特征。王小卫等[4]和胡鹏飞[5]从资料处理的角度总结了碳酸盐岩缝洞储集体成像的处理技术和手段, 指出速度建模精度是制约缝洞储集体“串珠”成像的重要因素。李鹏飞等[6]利用小尺度的网格层析速度建模策略, 提高了碳酸盐岩储层速度建模精度, 改善了缝洞体成像效果。崔永福等[7]将全波形反演技术应用于缝洞型储层探区的速度建模, 提高了火成岩建模精度, 间接改善了深部缝洞储层的成像质量。由于缝洞储层埋藏深, 常规采集资料无法满足全波形反演对排列长度的要求, 因此目前未见全波形反演应用于缝洞储层精细速度建模的公开文献。目前主要采用网格层析反演建模技术进行缝洞成像速度建模[8-10]。
利用基于成像道集拉平准则的网格层析建模技术, 只能反演得到低波数的背景速度场, 无法精细刻画实际地下介质的非均质特征。对缝洞储层速度建模而言, 质控时主要考虑大、中尺度“串珠”收敛程度, 往往忽略了对速度更敏感的小尺度缝洞体的成像质量, 最终反演得到的速度模型难以使所有缝洞储集体成像收敛效果最佳。对于缝洞精细开发而言, 大尺度缝洞不再是关注点, 收敛效果不佳或未收敛的小尺度缝洞体更为重要, 因此通过速度模型更新提高局部小尺度串珠的成像质量具有较强的实际意义。王世星等[11]提出了一种基于连续速度替换的时间域地震聚焦成像技术, 该技术采用偏移速度连续替换扫描和缝洞体反射能量聚焦分析两种方法, 调整缝洞体位置处的时间域偏移速度并完成叠前时间偏移成像处理, 提高了时间域缝洞体成像质量。
本文首先基于风化壳模型数据对比不同速度误差对深度域成像质量的影响, 分析成像剖面、成像道集与速度精度的敏感性关系, 进而提出了基于目标区速度扫描的深度域缝洞成像方法, 再对比Kirchhoff积分偏移、单程波偏移和逆时偏移的特征, 优选了一种单程波偏移成像方案, 并给出了目标区速度扫描单程波偏移成像技术流程, 最后将其应用于塔河某实际资料处理, 证明了本文方法的可行性和有效性。
1 不同速度风化壳缝洞成像分析为了分析速度模型对缝洞成像精度的影响, 设计了如图 1所示的风化壳缝洞体速度模型。该模型包括随机分布的2~10 m小尺度溶洞、20~50 m中尺度孤立溶洞和100~200 m大尺度孤立溶洞。溶洞体填充速度为4 500~5 000 m/s, 围岩速度为6 100 m/s。正演模拟参数与塔河实际采集参数基本一致, 炮间距30 m, 道间距30 m, 采样记录长度7 s, 采样间隔1 ms, 地震子波主频30 Hz。
图 2是利用准确速度模型进行逆时偏移(reverse time migration, RTM)获得的成像结果。由图 2可见, 大、中尺度溶洞均表现为明显的“串珠”状特征, 小尺度溶洞则表现为弱反射特征, 影响了风化壳顶面同相轴的下部分旁瓣, 但整体而言, 大、中尺度溶洞“串珠”收敛效果好, 小尺度溶洞成像质量高。
图 3给出了不同偏移速度模型对溶洞进行逆时偏移的成像结果, 图 3a至图 3d的偏移速度依次为准确速度的99.0%, 101.0%, 99.5%和100.5%。图中“串珠”形态圆润、无上翘或下拉的“尾巴”时表示缝洞收敛最佳; “串珠”边缘有下拉“尾巴”表示偏移不足(图 3a和图 3c); “串珠”边缘有上翘“尾巴”表示偏移过度(图 3b和图 3d)。由图 3可知, 当速度误差为1%时, “串珠”收敛质量很差, 当速度误差为0.5%时, “串珠”收敛质量与准确模型“串珠”成像结果相差较大。由图 3a和图 3c可见, 当偏移速度模型小于准确速度模型时, 小尺度溶洞表现为杂乱反射特征; 由图 3b和图 3d可见, 当偏移速度模型大于准确速度模型时, 溶洞成像表现为过偏移, 风化壳顶面同相轴的上部分旁瓣出现一系列能量异常。由此可见, “串珠”成像质量对速度模型精度非常敏感, 速度误差即便为0.5%, 依然会影响“串珠”成像的收敛效果。
速度建模基于成像道集拉平准则, 速度更新依赖于成像道集的剩余时差。我们提取图 2黄色线位置处(CDP点号为1 100)不同偏移速度的成像道集, 结果如图 4所示, 图 4中与黄色水平直线重合的同相轴对应风化壳顶面, 其下方的同相轴对应“串珠”, 图 4中两个道集的同相轴都很平, 当速度误差为0.5%时, 成像道集剩余时差非常小, 此时无法利用剩余时差更新速度模型。分析可知, 塔河缝洞体成像时若要获得精确的背景速度模型使“串珠”成像完全收敛非常困难。此外, 由于缝洞储层非均质性强, 速度横向差异大, 故很难使所有“串珠”均能准确成像并收敛。
对于缝洞地震资料而言, 奥陶系碳酸盐岩地层非均质强、速度模型复杂, 而其上覆地层层状沉积特征明显、速度相对简单。实际勘探中, 假设影响“串珠”收敛的主要因素是储层速度建模精度, 通常认为利用层析反演速度建模能够获得较高精度的上覆地层速度模型。基于上述假设, 考虑到速度误差严重制约“串珠”成像质量, 我们难以获得精确的速度模型, 为此可以在深度上选定一个包含缝洞储层的目标区, 在目标区范围内为层析反演速度模型选定一系列误差百分比, 利用不同的偏移速度进行偏移成像, 最终从多个偏移成像结果中寻找不同位置的“串珠”收敛最佳的结果。
由于扫描速度体数量通常比较多, 因此偏移方法必须计算效率高, 本文基于以下3个因素考虑选择单程波波动方程叠前深度偏移方法[12-16]: ①单程波偏移成像精度高, 能够适应横向变速介质成像, 对于构造形态复杂程度不高的缝洞探区而言, 其成像效果与逆时偏移基本相当, 且效果明显优于Kirchhoff积分偏移方法; ②单程波偏移计算效率高, 比逆时偏移速度高5倍左右, 其计算效率也优于Kirchhoff积分偏移方法; ③单程波偏移沿着深度方向进行波场延拓和成像, 在目标区深度以上速度精确的情况下, 可以将波场延拓到目标区上方, 再对目标区进行不同速度的速度扫描和偏移成像, 整体计算量较小, 而逆时偏移需要利用整体速度场进行波场传播, Kirchhoff积分偏移射线追踪时需要利用目标区以上的速度场, 二者均无法快速进行目标区速度扫描成像。
单程波波场传播方程为:
$ \begin{array}{l} \frac{{\partial U\left( {x, y, z, t} \right)}}{{\partial z}} = \\ \;\;\;\;\;\;\;\; \pm \sqrt {{{\left( {\frac{1}{{v\left( {x, y, z, t} \right)}}} \right)}^2} - {{\left( {\frac{{\partial t}}{{\partial x}}} \right)}^2} - {{\left( {\frac{{\partial t}}{{\partial y}}} \right)}^2}} \\ \;\;\;\;\;\;\;\;\frac{{\partial U\left( {x, y, z, t} \right)}}{{\partial t}} \end{array} $ | (1) |
式中: U(x, y, z; t)表示地下一点t时刻的地震波场; v(x, y, z)表示地下一点的介质速度。本文波场延拓时采用分步傅里叶算子(split-step Fourier, SSF), 每步波场延拓包含频率波数域的相移和频率空间域的时移, 表示如下:
$ \begin{array}{l} P_1^*\left( {{k_x}, {k_y}, z + \Delta z, \omega } \right) = {P^*}\left( {{k_x}, {k_y}, z, \omega } \right) \cdot \\ \;\;\;\;\;\;\;{\rm{exp}}\left\{ { \pm \sqrt {\frac{{{\omega ^2}}}{{{v_0}\left( z \right)}} - \left( {k_x^2 + k_y^2} \right)} \Delta z} \right\} \end{array} $ | (2) |
$ \begin{array}{l} P\left( {x, y, z + \Delta z, \omega } \right) = {P_1}\left( {x, y, z + \Delta z, \omega } \right) \cdot \\ \;\;\;\;\;\;\;{\rm{exp}}\left\{ { \pm {\rm{i}}\left[ {\frac{\omega }{{v\left( {x, y, z} \right)}} - \frac{\omega }{{{v_0}\left( z \right)}}} \right]\Delta z} \right\} \end{array} $ | (3) |
式中: kx、ky为水平方向的波数; ω为圆频率; v0(z)表示z深度的平均速度; Δz表示波场延拓的深度间隔; P(x, y, z, ω)表示z深度频率空间域的波场; P*(kx, ky, z, ω)为P(x, y, z, ω)变换到频率波数域的波场; P1*(kx, ky, z+Δz, ω)为频率波数域相移后深度z+Δz的波场; P1(x, y, z+Δz, ω)为P1*(kx, ky, z+Δz, ω)变换到频率空间域的波场; P(x, y, z+Δz, ω)为频率空间域时移后深度域z+Δz上的波场。
基于单程波方程的目标区速度扫描成像流程如下: ①给定目标区深度范围D1~D2, 在D1以上逐层进行震源波场、记录波场延拓和相关成像, 保存延拓到D1的两个波场; ②在目标区范围内, 将速度乘以扫描百分比, 从D1到D2逐层进行震源波场、记录波场延拓和相关成像; ③在D2深度成像后, 再取之前延拓到D1深度的两个波场, 按其它速度扫描比例再次从D1到D2逐层进行震源波场、记录波场延拓和相关成像。具体流程见图 5。
塔河油田某探区风化壳岩溶缝洞发育, 随着精细勘探开发的不断深入, 依靠大“串珠”开展井位部署的空间越来越受限, 不足以满足塔河的稳产上产需求。为了探索小尺度缝洞体成像, 选择该探区进行试验。地震数据共33 000炮, 数据量为2.1 TB, 面元大小为15 m×15 m。利用网格层析速度反演技术, 获得如图 6所示的速度模型。由图 6可见, 速度场整体较为平滑, 速度横向变化不大。根据风化壳成像要求, 将目标区深度设置为5.0~7.0 km。进行速度扫描测试, 将局部目标区速度扫描范围设定为网格层析速度的97.0%~103.0%, 即局部误差为±3%(前述数学模型分析时采用的是整体误差), 扫描个数为11, 频率范围2~70 Hz, 深度间隔10 m, 利用40个节点120个K10 GPU卡完成全部偏移计算, 耗时40.5 h。
图 7为目标区某条测线不同百分比速度扫描偏移成像剖面。由图 7a至图 7k发现“串珠”从欠偏移到过偏移的整个变化过程。图 7f是网格层析原始速度成像剖面, 其中绿色圈范围内“串珠”边缘略有上翘, 表示略微过收敛, 说明速度值略大, 而黄色圈范围内“串珠”明显欠收敛, 表明速度值偏低。利用97.0%速度扫描偏移成像时(图 7a), 绿色圈对应范围内的“串珠”收敛效果变好。当速度扫描百分比超过101.2%后(图 7h), 随着速度值的增大, 黄色圈对应范围内的“串珠”成像效果逐渐变好, 速度扫描百分比为103.0%时, 成像质量明显优于原始速度的成像效果。对比图 7f与图 7k可见, 风化壳内的小串珠成像质量明显提高, 部分原收敛质量较差或未收敛的小“串珠”已较好地收敛归位, 风化壳内的异常特征明显。由上述结果可知, 目标区速度扫描偏移成像方法对于识别小尺度缝洞体具有重要的实际应用价值。
本文方法提供了一系列成像结果, 其中某个结果中部分“串珠”聚焦性更好, 但部分“串珠”聚焦形态可能变差。如何从一系列成像结果中构建一个使得全部“串珠”都最佳收敛的速度模型, 仍有待进一步研究。
4 结束语碳酸盐岩缝洞成像效果依赖于速度模型精度, “串珠”成像质量对速度变化非常敏感。目前, 基于成像道集拉平准则的层析反演速度建模技术通常只能反演得到低波数的背景速度场, 无法精细刻画缝洞储层速度特征, 难以实现不同尺度“串珠”的完全收敛。本文提出了一种基于目标区速度扫描的单程波波动方程叠前深度偏移方法, 该方法基于目标区上覆地层速度模型精度较高的假设, 将震源波场和记录波场延拓到目标区上方, 利用不同百分比速度扫描逐次对目标区进行偏移成像, 得到一系列目标区不同速度的成像结果, 可用于提高部分“串珠”成像收敛精度。在塔河某探区的实际应用结果表明, 本文方法可使原收敛质量较差或未收敛的小溶洞形成的“串珠”较好地收敛, 在一定程度上提高了小溶洞“串珠”的识别能力。
从一系列地震成像剖面中选择合适的局部成像结果比较困难。我们下一步将研究自动化的“串珠”最佳收敛判别方法, 以及如何得到每一个串珠都能最佳收敛的速度模型, 利用该速度模型再次进行叠前深度偏移, 最终可获得整体最佳收敛的成像结果, 以进一步提升本文方法的实用性。
[1] |
韩革华, 漆立新, 李宗杰, 等. 塔河油田奥陶系碳酸盐岩缝洞型储层预测技术[J]. 石油与天然气地质, 2006, 27(6): 860-878. HAN G H, QI L X, LI Z J, et al. Prediction of the Ordovician fractured-vuggy carbonate reservoirs in Tahe oilfield[J]. Oil & Gas Geology, 2006, 27(6): 860-878. DOI:10.3321/j.issn:0253-9985.2006.06.017 |
[2] |
朱超, 刘占国, 杨少勇, 等. 利用相控分频反演预测英西湖相碳酸盐岩储层[J]. 石油地球物理勘探, 2018, 53(4): 832-841. ZHU C, LIU Z G, YANG S Y, et al. Lacustrine carbonate reservior prediction in Yingxi, Qaidam basin with the facies-constrained and segmented-frequency-band inversion[J]. Oil Geophysical Prospecting, 2018, 53(4): 832-841. |
[3] |
刘坤岩, 许杰. 塔河奥陶系隐蔽溶洞体地震精细识别[J]. 石油地球物理勘探, 2019, 54(5): 1106-1114. LIU K Y, XU J. Identification of Ordovician subtle cave reservoirs in Tahe on seismic data[J]. Oil Geophysical Prospecting, 2019, 54(5): 1106-1114. |
[4] |
王小卫, 吕磊, 刘伟方, 等. 塔里木盆地碳酸盐岩地震资料处理的几项关键技术[J]. 岩性油气藏, 2008, 20(4): 109-112. WANG X W, LV L, LIU W F, et al. Seismic data processing techniques of carbonate rocks in Tarim Basin[J]. Lithologic Reservoirs, 2008, 20(4): 109-112. DOI:10.3969/j.issn.1673-8926.2008.04.020 |
[5] |
胡鹏飞. 塔河油田碳酸盐岩缝洞型储集体成像技术研究[J]. 石油地球物理勘探, 2009, 44(2): 152-157. HU P F. Study on imaging technique of carbonate fracture-cave type's reservoir in Tahe Oilfield[J]. Oil Geophysical Prospecting, 2009, 44(2): 152-157. DOI:10.3321/j.issn:1000-7210.2009.02.006 |
[6] |
李鹏飞, 崔德育, 黄诚. 地震资料处理解释一体化技术在塔北碳酸盐岩储层识别中的应用[J]. 石油地球物理勘探, 2018, 53(S2): 306-312. LI P F, CUI D Y, HUANG C. Carbonate reservoir identification in Tabei Area, Tarim Basin with integrated seismic data processing and interpretation[J]. Oil Geophysical Prospecting, 2018, 53(S2): 306-312. |
[7] |
崔永福, 彭更新, 吴国忱, 等. 全波形反演在缝洞型储层速度建模中的应用[J]. 地球物理学报, 2016, 59(7): 2714-2725. CUI Y F, PENG G X, WU G C, et al. Application of full waveform inversion velocity model-building technology for the fractured-vuggy reservoir[J]. Chinese Journal of Geophysics, 2016, 59(7): 2714-2725. |
[8] |
闫艳琴, 费建博, 闫亮. 叠前深度偏移在跃参地区地震资料处理中的应用[J]. 工程地球物理学报, 2019, 16(6): 784-791. YAN Y Q, FEI J B, YAN L, et al. Application of pre-stack depth migration in seismic data processing in Yuecan area[J]. Chinese Journal of Engineering Geophysics, 2019, 16(6): 784-791. DOI:10.3969/j.issn.1672-7940.2019.06.002 |
[9] |
李鹏, 刘志成, 杨楠, 等. 玉北地区碳酸盐岩缝洞型储层精细成像方法应用研究[J]. 石油物探, 2015, 54(4): 443-451, 484. LI P, LIU Z C, YANG N, et al. Application of fine target reservoir imaging in fractured-vuggy carbonate reservoir Yubei area[J]. Geophysical Prospecting for Petroleum, 2015, 54(4): 443-451, 484. DOI:10.3969/j.issn.1000-1441.2015.04.011 |
[10] |
朱立华, 马灵伟, 白英哲, 等. 顺北超深断控体地震多信息约束速度建模[J]. 石油物探, 2019, 58(6): 864-873. ZHU L H, MA L W, BAI Y Z, et al. Velocity modeling of an ultra-deep fault-controlled reservoir in Shunbei area China, using a multi-information constraint[J]. Geophysical Prospecting for Petroleum, 2019, 58(6): 864-873. DOI:10.3969/j.issn.1000-1441.2019.06.009 |
[11] |
王世星, 许凯, 唐金良. 岩溶缝洞储集体地震聚焦成像与充填速度检测[J]. 石油物探, 2019, 58(6): 898-910. WANG S X, XU K, TANG J L. Seismic focused imaging and filling-velocity prediction of karst fracture-cavern reservoirs[J]. Geophysical Prospecting for Petroleum, 2019, 58(6): 898-910. DOI:10.3969/j.issn.1000-1441.2019.06.013 |
[12] |
STOFA P L, FOKKEMA J T, LUNA F, et al. Split-step Fourier migration[J]. Geophysics, 1990, 55(2): 410-421. |
[13] |
RISTOW D, RUHL T. Fourier finite-difference migration[J]. Geophysics, 1994, 59(12): 1882-1893. DOI:10.1190/1.1443575 |
[14] |
BERKHOUT A J. Wave field extrapolation techniques in seismic migration[J]. Geophysics, 1981, 46(11): 1638-1656. |
[15] |
张宇, 徐升, 张关泉, 等. 真振幅全倾角单程波方程偏移方法[J]. 石油物探, 2007, 46(6): 582-587. ZHANG Y, XU S, ZHANG G Q, et al. True amplitude turning-wave one-way wave equation migration[J]. Geophysical Prospecting for Petroleum, 2007, 46(6): 582-587. DOI:10.3969/j.issn.1000-1441.2007.06.010 |
[16] |
段心标, 王华忠, 白英哲, 等. 基于GPU的三维起伏地表单程波叠前深度偏移[J]. 石油物探, 2016, 55(2): 223-230. DUAN X B, WANG H Z, BAI Y Z, et al. 3D one-way wave equation prestack depth migration from topography based on the acceleration of GPU[J]. Geophysical Prospecting for Petroleum, 2016, 55(2): 223-230. |