Exact reflection coefficients of orthorhombic anisotropic media
ZHANG Xueying1, SUN Pengyuan2, MA Xuejun3, LU Jun4, LI Mengqi4
1. School of Energy Resources, China University of Geosciences(Beijing), Beijing 100083, China; 2. Geophysical Research and Development Center, BGP, CNPC, Zhuozhou, Hebei 072751, China; 3. Petroleum Exploration and Production Research Institute, Sinopec Northwest Oilfield Company, Urumqi, Xinjiang 830011, China; 4. School of Geophysics and Information Technology, China University of Geosciences(Beijing), Beijing 100083, China;
Abstract:Most of continental sedimentary formations in China are characterized by thin interbeds.Under the influence of tectonic movement,in sedimentary formations,vertical or nearly vertical high-dip fractures are developed,which result in typical orthorhombic anisotropy.This kind of orthorhombic anisotropy can be regarded as the superposition of VTI and HTI.Therefore,the study on AVO responses of orthorhombic anisotropic media is significant for fine characterization of fractured reservoirs.To the orthor-hombic anisotropy induced by a set of upright frac-tures developed in a VTI background,Tsvankin aniso- tropy parameters are used to construct the stiffness coefficient matrix.Then,the method for calculating exact reflection and transmission coefficients is presented based on the Christoffel equations and the correspon-ding boundary conditions.Model test shows that: ①if P-wave is incident from a low-impedance isotropic medium to a high-impedance orthorhombic anisotropic medium,the increase of the VTI anisotropy coefficient in the background media will lead to the increase of the PP-wave reflection coefficient and the decrease of the PS1- and PS2-wave reflection coefficients; and the increase of fracture weakness will lead to the decrease of the PP-wave reflection coefficient and the increase of the PS1- and PS2-wave reflection coefficients; ②if P-wave is incident from a high-impedance isotropic medium to a low-impedance orthorhombic anisotropic medium,the reflection coefficients of PP-,PS1- and PS2-waves present the opposite characteristics; ③the azimuthal anisotropy of PP-wave reflection coefficients are much smaller than those of PS1- and PS2-waves.By comparing with the results of Rüger approximate reflection coefficient formulas,the accuracy of our calculation method is verified.
Christoffel E B.Ueber die Fortpflanzung von Stössen durch elastische feste Körper[J].Annali Di Matematica Pura Ed Applicata,1877,8(1):193-243.
[4]
Crampin S.Seismic wave propagation through a cracked solid:Polarization as a possible dilatancy diagnostic[J].Geophysical Journal of the Royal Astronomical Society,1978,53(3):467-496.
[5]
Hudson J A.A higher order approximation to the wave propagation constants for cracked solid[J].Geophysical Journal International,1986,87(1):265-274.
[6]
Hudson J A.Crack distribution which account for a given seismic anisotropy[J].Geophysical Journal International,1991,104(3):517-521.
[7]
Schoenberg M.Reflection of elastic waves from perio-dically stratified media with interfacial slip[J].Geophysical Prospecting,1983,31(2):265-292.
[8]
Schoenberg M,Sayers C M.Seismic anisotropy of fractured rocks[J].Geophysics,1995,60(1):204-211.
Tsvankin I.Anisotropic parameters and P-wave velo-city for orthorhombic media[J].Geophysics,1997,62(4):1292-1309.
[11]
Bakulin A,Grechka V,Tsvankin I.Estimation of fracture parameters from reflection seismic data-Part Ⅰ:HTI model due to a single fracture set[J].Geophy-sics,2000,65(6):1788-1802.
[12]
Bakulin A,Grechka V,Tsvankin I.Estimation of fracture parameters from reflection seismic data-Part Ⅱ:Fractured models with orthorhombic symmetry[J].Geophysics,2000,65(6):1803-1817.
[13]
Bakulin A,Grechka V,Tsvankin I.Estimation of fracture parameters from reflection seismic data-Part Ⅲ:Fractured models with monoclinic symmetry[J].Geophysics,2000,65(6):1818-1830.
[14]
Rokhlin S I.Reflection and refraction of elastic waves on a plane interface between two generally anisotropic media[J].The Journal of the Acoustical Society of America,1986,79(4):906-918.
[15]
Rüger A.P-ware reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry[J].Geophysics,1997,62(3):713-721.
[16]
Rüger A.Variation of P-wave reflectivity with offset and azimuth in anisotropic media[J].Geophysics,1998,63(3):935-947.
[17]
Rüger A.Reflection Coefficients and Azimuthal AVO Analysis in Anisotropic Media[M].SEG,Tulsa,2002.
[18]
傅旦丹,何樵登.含垂直裂隙PTL介质长波长等效正交各向异性弹性常数[J].地球物理学报,2002,45(增刊):307-320.FU Dandan,HE Qiaodeng.Elastic constants of long-wavelength effective orthorhombic anisotropy for PTL media containing vertical cracks[J].Chinese Journal of Geophysics,2002,45(S):307-320.
[19]
张文生,宋海斌.三维正交各向异性介质三分量高精度有限差分正演模拟[J].石油地球物理勘探,2001,36(4):422-432.ZHANG Wensheng,SONG Haibin.Finite-difference forward modeling of three component records with high precision in 3-D orthorhombic anisotropic media[J].Oil Geophysical Prospecting,2001,36(4):422-432.
[20]
魏建新.正交各向异性介质波场特性的物理模拟研究[D].北京:中国地震局地球物理研究所,2003.
[21]
李娜.正交各向异性介质正演模拟方法及波场特征分析[D].山东青岛:中国石油大学(华东),2014.
[22]
许茜茹,毛伟建.正交各向异性介质中qP波入射的二阶近似反射系数与透射系数[J].地球物理学报,2019,62(1):266-281.XU Qianru,MAO Weijian.Reflection and transmi-ssion coefficients based on second-order approximation in orthorhombic anisotropic media[J].Chinese Journal of Geophysics,2019,62(1):266-281.
[23]
张繁昌,路亚威,桑凯恒,等.正交介质方位傅里叶反射系数及裂缝密度预测[J].中国石油大学学报(自然科学版),2019,43(2):37-44.ZHANG Fanchang,LU Yawei,SANG Kaiheng,et al.Azimuthal Fourier reflected coefficients in prediction of orthogonal fracture parameters[J].Journal of China University of Petroleum(Edition of Natural Science),2019,43(2):37-44.
[24]
唐杰,温雷,王浩,等.正交各向异性介质中的剪张源震源机制与矩张量特征[J].石油地球物理勘探,2018,53(6):1247-1255.TANG Jie,WEN Lei,WANG Hao,et al.Focal mecha-nisms and moment tensor in orthorhombic anisotropic media[J].Oil Geophysical Prospecting,2018,53(6):1247-1255.
[25]
郭恺,张永升,冯彦彬.正交各向异性介质多参数建模方法研究[J].石油物探,2019,58(6):874-881.GUO Kai,ZHANG Yongsheng,FENG Yanbin.Multi-parameter modeling of orthotropic media[J].Geophysical Prospecting for Petroleum,2019,58(6):874-881.
[26]
单俊臻,吴国忱,龚诚诚.HTI介质方位观测PP波反射系数一阶扰动近似[J].石油地球物理勘探,2019,54(2):371-379.SHAN Junzhen,WU Guochen,GONG Chengcheng.First-order perturbation approximation of PP wave reflection coefficient for HTI medium[J].Oil Geophysical Prospecting,2019,54(2):371-379.
[27]
秦海旭,吴国忱,刘争平.裂缝单斜介质qP波正演模拟方法[J].石油地球物理勘探,2015,50(6):1118-1126.QIN Haixu,WU Guochen,LIU Zhengping.qP wave forward modeling in fracture media with monoclinic symmetry[J].Oil Geophysical Prospecting,2015,50(6):1118-1126.
[28]
李春鹏,印兴耀,张峰.HTI介质饱和流体特性和裂缝密度对方位反射系数的影响[J].石油物探,2013,52(1):1-10.LI Chunpeng,YIN Xingyao,ZHANG Feng.Influence of HTI medium saturated fluid properties and fracture density on azimuthal reflectivity[J].Geophysical Prospecting for Petroleum,2013,52(1):1-10.
[29]
Schoenberg M,Protázio J."Zoeppritz" rationalized and generalized to anisotropy[J].Journal of the Seismic Exploration,1992,1(2):125-144.
[30]
Zhou B,Greenhalgh S.On the computation of elastic wave group velocities for a general anisotropic me-dium[J].Journal of Geophysics and Engineering,2004,1(3):205-215.
[31]
Chen H.Anisotropic Effects upon Amplitude-vs-offset Response in Realistic Earth Models[D].University of Oklahoma Graduate College,2000.
[32]
Schoenberg I,Helbig K.Orthorhombic media:mode-ling elastic wave behavior in a vertically fractured earth[J].Geophysics,1997,62(5):1954-1974.