各向异性介质地震波场数值模拟及特征分析  

The Seismic Wave Numerical Simulation and Characteristics Analysis in Anisotropic Medium

作　　者：龚明平 张军华[1] 刘文革[2] 刘杨 吴系源 GONG Ming-ping;ZHANG Jun-hua;LIU Wen-ge;LIU Yang;WU Xi-yuan(School of Geosciences, China University of Petroleum （East China）, Qingdao 266580, China;School of Geoscience and Technology, Southwest Petroleum University, Chengdu 610500, China)

机构地区：[1]中国石油大学(华东)地球科学与技术学院,山东青岛266580 [2]西南石油大学地球科学与技术学院,成都610500

出　　处：《CT理论与应用研究（中英文）》2017年第6期661-668,共8页Computerized Tomography Theory and Applications

基　　金：国家科技重大专项(2016ZX05011-002;2017ZX05072-001)

摘　　要：由于各向异性广泛存在于地下岩石中,随着勘探精度的不断提高,对地下介质的各向同性假设越来越不能够满足于现状,因此对各向异性介质的数值模拟显得更为重要。本文推导了各向异性介质的弹性波动方程,总结了震源类型,通过PML方法处理了人工边界问题,通过快照分析验证了数值频散、稳定性条件。研究结果表明:(1)PML完全匹配层,可较好地解决人工边界问题;(2)减小空间采样间隔压制数值频散比减小时间采样间隔效果要好得多,盲目减小时间采样间隔会大大降低数值模拟的运算效率;(3)各向异性介质中弹性波场中除含有准纵波外,还含有速度较慢的准横波;(4)准纵波波前能量要比由各向异性引起的准横波能量强,准纵波和准横波的波前随着各向异性介质参数的变化而变化。Due to the anisotropy widely exists in underground rocks, and with the increasing of exploration precision, isotropic hypothesis of underground media can＇t be gradually content with the status quo, so the numerical simulation of the anisotropic medium is more important. This paper derived the elastic wave equation of anisotropic medium, summarized the source type, tackled the problem of the artificial boundary by PML method, and verified the numerical dispersion and stability conditions by snapshots. The results show that： 1） PML based on staggered grid can solve the problem of the artificial boundary; 2） reduce the spatial sampling interval is much better than reduce time sampling interval to suppress numerical dispersion and blindly reduce time sampling interval can reduce the operational efficiency of the numerical simulation; 3）there is quasi p-wave and slow quasi s-wave in the elastic wave field of anisotropic medium; 4） quasi p-wave front energy is stronger than quasi s-wave energy, quasi p- wave and quasi s-wave front are changed as parameters of anisotropic medium.

关 键 词：各向异性介质 数值模拟 有限差分 数值频散 

分 类 号：P631.4[天文地球—地质矿产勘探]