|
|
Prestack inversion of anisotropic coefficients and its application in fracture prediction |
XIE Rui( ), YAN Jian-Guo( ), CHEN Qi |
Geophysical College,Chengdu University of Technology,Chengdu 610059, China |
|
|
Abstract The forward modeling of wave field response characteristics of wide-azimuth seismic data shows that the reflection coefficient has apparent azimuth anisotropy when the P-waves propagate in underground geologic bodies.It has become a research hotspot at home and abroad to make a fracture prediction using azimuth anisotropy.This study analyzed the response characteristics of the changes in seismic P-waves with azimuth using forward modeling.Based on the analytical results combined with the relationship between the anisotropy parameters of fluid-bearing fractures proposed by Bakulin et al.,this study approximately simplified the Ruger formula,deduced the equation between anisotropic parameter γ that can be used to characterize fracture development and reflection coefficient,and proposed a fracture prediction method based on anisotropic parameter inversion.The validity and applicability of the method have been proven using theoretical models and practical data application.Therefore,this study provides a feasible method and technology for fracture prediction using wide-azimuth prestack seismic data.
|
Received: 21 August 2021
Published: 17 August 2022
|
|
Corresponding Authors:
YAN Jian-Guo
E-mail: 2429259471@qq.com;yanjianguo@cdut.edu.cn
|
|
|
|
|
Variation of P-wave reflection coefficient with azimuth in isotropic and HTI media models
|
|
Variation of P-wave reflection coefficient of HTI medium model with anisotropic parameters
|
参数 | vp/(m·s-1) | vs/(m·s-1) | ρ | φ | θ | Layer 1 | 2000 | 1253 | 2.2 | 45 | 30 | Layer 2 | 2337 | 1437 | 2.3 | 45 | 30 |
|
Variation of P-wave reflection coefficient of HTI medium with anisotropic parameters
|
|
Flow chart of anisotropic coefficient inversion
|
vp/(m·s-1) | vs/(m·s-1) | ρ/(g·cm-3) | ξ | δ | γ | θ | φ | 4000 | 2300 | 2.8 | 0 | 0 | 0 | 0~60 | 0~180 | 5500 | 3100 | 2.65 | 0.05 | 0.25 | 0.2 | | |
|
According to the data of the study area, the model parameters are established
|
|
Comparison of isotropic calculation results of different equations with incident angle
|
|
Comparison of anisotropic terms of different equations with incident angle(a) and azimuth(b)
|
|
Prestack gathers are synthesized at different azimuth angles at the target layer of the corresponding well in the study area
|
|
Anisotropic parameter curve inversion results theoretical results and inversion results at well BZ-B based on least square method a—longitudinal wave impedance reflection coefficient;b—transverse wave impedance reflection coefficient;c—anisotropic parameter γ
|
|
Seismic well connection profile in the study area
|
|
AVAZ prestack fracture prediction well connection profile
|
|
Prestack anisotropic parameter inversion fracture prediction well connection profile
|
数据体编号 | 入射角范围/(°) | 中心入射角/(°) | 方位角范围/(°) | 中心方位角/(°) | 1 | 10~30 | 15 | 30~60 | 45 | 2 | 20~30 | 25 | 90~120 | 105 | 3 | 20~40 | 30 | 120~150 | 135 | 4 | 30~40 | 35 | 150~180 | 165 |
|
Division of incident angle and azimuth selected for actual data inversion
|
[1] |
关宝文, 郭建明, 杨燕, 等. 油气储层裂缝预测方法及发展趋势[J]. 特种油气藏, 2014, 21(1):12-17.
|
[1] |
Guan B W, Guo J M, Yang Y, et al. Methods of fracture prediction in oil & gas reservoirs and their development trend[J]. Special Oil & Gas Reservoirs, 2014, 21(1):12-17.
|
[2] |
张昕, 郑晓东. 裂缝发育带地震识别预测技术研究进展[J]. 石油地球物理勘探, 2005, 40(6):724-731.
|
[2] |
Zhang X, Zheng X D. Research progress of seismic identification and prediction technology in fracture zone[J]. Oil Geophysical Prospecting, 2005, 40(6):724-731.
|
[3] |
印兴耀, 张洪学, 宗兆云. OVT数据域五维地震资料解释技术研究现状与进展[J]. 石油物探, 2018, 57(2):155-178.
|
[3] |
Yin X Y, Zhang H X, Zong Z Y. Research status and progress of 5D seismic data interpretation in OVT domain[J]. Geophysical Prospecting for Petroleum, 2018, 57(2):155-178.
|
[4] |
Crampin S. Seismic-wave propagation through a cracked solid:polarization as a possible dilatancy diagnostic[J]. Geophysical Journal International, 1978, 53(3):467-496.
|
[5] |
Thomsen L. Elastic anisotropy due to aligned cracks in porous rock[J]. Geophysical Prospecting, 1995, 43(6):805-829.
|
[6] |
张琬璐. 基于Rüger公式的地震各向异性参数反演方法研究[D]. 北京: 中国石油大学, 2016.
|
[6] |
Zhang W L. Research on seismic anisotropic parameter inversion method based on Rüger’s formula[D]. Beijing: China University of Petroleum, 2016.
|
[7] |
任涛, 秦军, 周阳, 等. 利用叠前方位各向异性预测准噶尔盆地西北缘车476井区火山岩裂缝[J]. 石油地质与工程, 2021, 35(1):1-7.
|
[7] |
Ren T, Qin J, Zhou Y, et al. Prediction of volcanic fractures in Che 476 well block in northwestern margin of Junggar basin by using pre-stack azimuthal anisotropy[J]. Petroleum Geology and Engineering, 2021, 35(1):1-7.
|
[8] |
Rüger A. P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry[J]. Geophysics, 1997, 62(3):713-721.
|
[9] |
陈怀震, 印兴耀, 杜炳毅, 等. 裂缝型碳酸盐岩储层方位各向异性弹性阻抗反演[J]. 地球物理学进展, 2013, 28(6):3073-3079.
|
[9] |
Cheng H Z, Yin X Y, Du B Y, et al. Current situations and recent progress in different azimuths seismic exploration[J]. Progress in Geophysics, 2013, 28(6):3073-3079.
|
[10] |
Bakulin A, Vladimir G, Tsvankin Ilya.Estimation of fracture parameters from reflection seismic data—Part I:HTI model due to a single fracture set[J]. Geophysics, 2000, 65(6):1788-1802.
|
[11] |
龚明平, 张军华, 王延光, 等. 分方位地震勘探研究现状及进展[J]. 石油地球物理勘探, 2018, 53(3):642-658.
|
[11] |
Gong M P, Zhang J H, Wang Y G, et al. Current situations and recent progress in different azimuths seismic exploration[J]. Oil Geophysical Prospecting, 2018, 53(3):642-658.
|
[12] |
张保庆, 周辉, 左黄金, 等. 宽方位地震资料处理技术及应用效果[J]. 石油地球物理勘探, 2011, 46(3):396-400.
|
[12] |
Zhang B Q, Zhou H, Zuo H J, et al. Wide azimuth data processing techniques and their applications[J]. Oil Geophysical Prospecting, 2011, 46(3):396-400.
|
|
|
|