基于BISQ方程的高频极限,导出了适用于任意方位角的双相HTI介质波动方程,推导了该方程在交错网格空间中求解的高阶有限差分格式和对应的完全匹配层(PML)吸收边界条件,实现了该类介质的地震波场正演模拟。模拟结果表明,该方法能准确模拟地震波在双相HTI介质中的传播过程,得到高精度的波场快照和合成记录。分析得出:在双相各向异性介质中,弹性波的传播仍表现出各向异性特征;孔隙度、粘滞系数和固液耦合密度主要影响慢纵波的振幅和形态。
Based on the high-frequency limit of BISQ equation,this paper derived the dual-phase HTI media elastic wave equation that could be used in the arbitrary fracture azimuth angle and deduced the high-order finite difference scheme of the wave equation in the space of staggered-grid and the corresponding perfectly matched layer(PML) absorbing boundary conditions,realizing the simulation of seismic wave field for that kind of media.The simulation results showed that the method could accurately simulate the seismic wave propagation in the dual-phase HTI media and could get accurate snapshots and synthetic records.Additionally,the conclusion by analysis was that in the dual-phase anisotropic media,the propagation of elastic wave still showed the anisotropic feature and the porosity,the coefficient of viscosity and the density of solid-liquid coupling could mainly affect the amplitude and the morphology of slow longitudinal wave.