Wavefield separation of P- and S-waves based on modified pseudo-Helmholtz decomposition operator in transverse isotropic media

Authors

LIU Zhaoshun, Shandong Provincial Key Laboratory of Deep Oil & Gas, China University of Petroleum (East China), Qingdao 266580, ChinaFollow
DU Qizhen, Shandong Provincial Key Laboratory of Deep Oil & Gas, China University of Petroleum (East China), Qingdao 266580, China; Qingdao National Laboratory for Marine Science and Technology, Laboratory for Marine Mineral Resources, Qingdao 266237, ChinaFollow
LYU Wenhao, Shandong Provincial Key Laboratory of Deep Oil & Gas, China University of Petroleum (East China), Qingdao 266580, China
ZHOU Shihao, Shandong Provincial Key Laboratory of Deep Oil & Gas, China University of Petroleum (East China), Qingdao 266580, China
FU Liyun, Shandong Provincial Key Laboratory of Deep Oil & Gas, China University of Petroleum (East China), Qingdao 266580, China; Qingdao National Laboratory for Marine Science and Technology, Laboratory for Marine Mineral Resources, Qingdao 266237, China

Abstract

The separation of P- and S-wave fields is a significant step in elastic reverse-time migration (ERTM), which can effectively eliminate the waveform crosstalk and image artifacts, thus improving the precision of imaging. In anisotropic media, nonstationary spatial filters or low-rank approximation methods are commonly used to separate the wavefields. However, the application of multiple Fourier transforms results in the high computation costs of the above wavefield separation methods. Based on the idea of constructing the isotropic decoupling equations of P- and S- waves with the Helmholtz operator, a modified pseudo-Helmholtz decomposition operator was proposed to eliminate the amplitude distortion. Then, the expressions of shear and dilatational waves in transverse isotropic media were solved with the method of undetermined coefficients. Meanwhile, these expressions were converted into the first-order decoupled pseudo-elastic wave equations, so as to realize the wavefield decoupling of P- and S-waves in time and space. Through the wavefield separation test of the simple model, the separated P and S-wave fields were obtained, which verified the effectiveness of the method for wavefield separation. Furthermore, the vector P- and S-waves obtained by decoupling were applied to the elastic reverse time migration, and the clear imaging results of elastic reverse-time migration were obtained with the vector dot product cross-correlation imaging conditions. Thereby, it is indicated that this method is well applicable to the complex media. Meanwhile, it could be effectively applied to the elastic reverse-time migration of VTI media.

Keywords

VTI media, elastic wave, wavefield separation, reverse time migration, Helmholtz decomposition

DOI

10.12363/issn.1001-1986.22.06.0477

Recommended Citation

LIU Zhaoshun, DU Qizhen, LYU Wenhao, et al. (2023) "Wavefield separation of P- and S-waves based on modified pseudo-Helmholtz decomposition operator in transverse isotropic media," Coal Geology & Exploration: Vol. 51: Iss. 3, Article 58.
DOI: 10.12363/issn.1001-1986.22.06.0477
Available at: https://cge.researchcommons.org/journal/vol51/iss3/58

Reference

[1] BAYSAL E,KOSLOFF D D,SHERWOOD J W C. Reverse time migration[J]. Geophysics,1983,48(11):1514−1524.

[2] MCMECHAN G A. Migration by extrapolation of time–dependent boundary values[J]. Geophysical Prospecting,1983,31(3):413−420.

[3] WHITMORE N D. Iterative depth migration by backward time propagation[C]//SEG Technical Program Expanded Abstracts 1983. Society of Exploration Geophysicists，1983：382–385.

[4] AKI K，RICHARDS P G. Quantitative seismology：Theory and methods[M]. San Francisco，CA：Freeman，1980.

[5] SUN R,MCMECHAN G A. Scalar reverse–time depth migration of prestack elastic seismic data[J]. Geophysics,2001,66(5):1519−1527.

[6] YAN Jia,SAVA P. Isotropic angle–domain elastic reverse–time migration[J]. Geophysics,2008,73(6):S229−S239.

[7] DU Qizhen,ZHU Yitong,BA Jing. Polarity reversal correction for elastic reverse time migration[J]. Geophysics,2012,77(2):S31−S41.

[8] DUAN Yuting,SAVA P. Scalar imaging condition for elastic reverse time migration[J]. Geophysics,2015,80(4):S127−S136.

[9] ZHU Hejun. Elastic wavefield separation based on the Helmholtz decomposition[J]. Geophysics,2017,82(2):S173−S183.

[10] YANG Jidong,ZHU Hejun,WANG Wenlong,et al. Isotropic elastic reverse time migration using the phase−and amplitude–corrected vector P– and S–wavefields[J]. Geophysics,2018,83(6):S489−S503.

[11] 马德堂,朱光明. 弹性波波场 P 波和 S 波分解的数值模拟[J]. 石油地球物理勘探,2003,38(5):482−486.

MA Detang,ZHU Guangming. Numerical modeling of P–wave and S–wave separation in elastic wavefield[J]. Oil Geophysical Prospecting,2003,38(5):482−486.

[12] 李振春,张华,刘庆敏,等. 弹性波交错网格高阶有限差分法波场分离数值模拟[J]. 石油地球物理勘探,2007,42(5):510−515.

LI Zhenchun,ZHANG Hua,LIU Qingmin,et al. Numeric simulation of elastic wavefield separation by staggering grid high–order finite–difference algorithm[J]. Oil Geophysical Prospecting,2007,42(5):510−515.

[13] XIAO Xiang,LEANEY W S. Local vertical seismic profiling (VSP) elastic reverse–time migration and migration resolution:Salt–flank imaging with transmitted P–to–S waves[J]. Geophysics,2010,75(2):S35−S49.

[14] DU Qizhen,GUO Chengfeng,ZHAO Qiang,et al. Vector–based elastic reverse time migration based on scalar imaging condition[J]. Geophysics,2017,82(2):S111−S127.

[15] TANG Chen,MCMECHAN G A. Multidirectional–vector–based elastic reverse time migration and angle–domain common–image gathers with approximate wavefield decomposition of P– and S–waves[J]. Geophysics,2018,83(1):S57−S79.

[16] DELLINGER J,ETGEN J. Wave–field separation in two–dimensional anisotropic media[J]. Geophysics,1990,55(7):914−919.

[17] YAN Jia,SAVA P. Elastic wave–mode separation for VTI media[J]. Geophysics,2009,74(5):WB19−WB32.

[18] ZHANG Qunshan,MCMECHAN G A. 2D and 3D elastic wavefield vector decomposition in the wavenumber domain for VTI media[J]. Geophysics,2010,75(3):D13−D26.

[19] CHENG Jiubing,FOMEL S. Fast algorithms for elastic−wave−mode separation and vector decomposition using low−rank approximation for anisotropic media[J]. Geophysics,2014,79(4):C97−C110.

[20] 郭成锋. 横向各向同性介质声波与弹性波逆时偏移[D]. 青岛：中国石油大学(华东)，2016.

GUO Chengfeng. Study on acoustic and elastic reverse time migration in transversely isotropic media[D]. Qingdao：China University of Petroleum (East China)，2016.

[21] ZHOU Yang,WANG Huazhong. Efficient wave–mode separation in vertical transversely isotropic media[J]. Geophysics,2017,82(2):C35−C47.

[22] LU Yongming,LIU Qiancheng,ZHANG Jianfeng,et al. Poynting and polarization vectors based wavefield decomposition and their application on elastic reverse time migration in 2D transversely isotropic media[J]. Geophysical Prospecting,2019,67(5):1296−1311.

[23] WANG Wenlong,HUA Biaolong,MCMECHAN G A,et al. P– and S−decomposition in anisotropic media with localized low–rank approximations[J]. Geophysics,2018,83(1):C13−C26.

[24] YANG Jidong,ZHANG Houzhu,ZHAO Yang,et al. Elastic wavefield separation in anisotropic media based on eigenform analysis and its application in reverse–time migration[J]. Geophysical Journal International,2019,217(2):1290−1313.

[25] 周进举,王德利,李博文,等. 基于解耦传播的波场分解方法在VTI介质弹性波逆时偏移中的应用[J]. 吉林大学学报(地球科学版),2018,48(3):909−921.

ZHOU Jinju,WANG Deli,LI Bowen,et al. Application of wavefield decomposition based on decoupled propagation in elastic RTM for VTI media[J]. Journal of Jilin University (Earth Science Edition),2018,48(3):909−921.

[26] 孔丽云,胡志方,王一博. 基于“方程解耦”算法的各向异性介质纵横波分离方法[J]. 地球物理学报,2022,65(7):2719−2728.

KONG Liyun,HU Zhifang,WANG Yibo. Elastic wave–mode separation for anisotropic media based on equation–decoupling method[J]. Chinese Journal of Geophysics,2022,65(7):2719−2728.

[27] DU Qizhen,GONG Xufei,ZHANG Mingqiang,et al. 3D PS–wave imaging with elastic reverse–time migration[J]. Geophysics,2014,79(5):S173−S184.

[28] WANG Wenlong,MCMECHAN G A. Vector–based elastic reverse time migration[J]. Geophysics,2015,80(6):S245−S258.

[29] THOMSEN L. Weak elastic anisotropy[J]. Geophysics,1986,51(10):1954−1966.

[30] 李宾. 横向各向同性介质有限差分法波场模拟方法研究[D]. 青岛：中国石油大学(华东)，2009.

LI Bin. Finite difference numerical modeling of seismic wavefields in tansversely isotropic media with a vertical axis[D]. Qingdao：China University of Petroleum (East China)，2009.