时间域黏滞波动方程及其数值模拟新方法  

A novel time-domain viscoacoustic wave equation and its numerical simulation

作　　者：罗文山[1,2] 陈汉明 王成祥 周辉[3] 王狮虎 

机构地区：[1]成都理工大学,四川成都610059 [2]东方地球物理公司物探技术研究中心,河北涿州072751 [3]油气资源与探测国家重点实验室,CNPC物探重点实验室,中国石油大学(北京),北京102249

出　　处：《石油地球物理勘探》2016年第4期707-713,4,共7页Oil Geophysical Prospecting

基　　金：国家“973”项目(2013CB228603);国家自然科学基金项目(41174119)联合资助

摘　　要：基于常Q模型的应力—应变关系,推导了分数阶拉普拉斯算子黏滞波动方程的一阶速度—压力形式,与现有的一阶速度—压力—应变形式相比,新推导的方程形式更简单,数值模拟时耗费内存更少。针对新推导的一阶速度—压力黏滞波动方程,采用交错网格伪谱法进行数值模拟,并利用卷积型完美匹配层(CPML)作为吸收边界,压制截断边界的反射。数值模拟实验证实,新推导的黏滞波动方程能很好地描述地震波在黏滞介质中的衰减和频散,交错网格伪谱法和CPML的组合是一种高效的数值模拟方法。We develop a fractional Laplacian viscoacoustic wave equation based on the dispersion relation of the constant-Q model, and formulate the wave equation into the first-order velocity-pressure system. Compared with the existing velocity-pressure-strain formulation, our formulation is more compact, and saves memory storage after discretization. The staggered-grid pseudo-spectral (SGPS) method is adopted to numerically solve our viscoacoustic wave equation. The convolutional perfectly matched layer (CPML) is applied to suppress artificial edge reflections in numerical simulation. Test results demonstrate that numerical simulation of our viscoacoustic wave can correctly describe seismic wave attenuation and dispersion. The SGPS approach combined with the CPML is verified to be an efficient numerical simulation scheme. © 2016, Editorial Department OIL GEOPHYSICAL PROSPECTING. All right reserved.

关 键 词：时间域 黏滞波动方程 伪谱法 卷积型完美匹配层 

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