梯形网格伪谱法地震波场模拟  

Seismic wave simulation using a trapezoid grid pseudo-spectral method

作　　者：谭文卓 吴帮玉 李博[2] 雷军 TAN Wenzhuo;WU Bangyu;LI Bo;LEI Jun(School of Mathematics and Statistics,Xi’an JiaotongUniversity,Xi’an,Shaanxi 710049,China;Sinopec INOPEC Geophysical Research Institute,Nanjing,Jiangsu 211103,China;The Sixth Gas Production Plant,Changqing OilfieldCompany,PetroChina,Xi’an,Shaanxi 710018,China)

机构地区：[1]西安交通大学数学与统计学院,陕西西安710049 [2]中国石化石油物探技术研究院,江苏南京211103 [3]中国石油长庆油田分公司第六采气厂,陕西西安710018

出　　处：《石油地球物理勘探》2020年第6期1282-1291,1163-1164,共12页Oil Geophysical Prospecting

基　　金：国家自然科学基金项目“致密气储层弱反射特征的地震信号高分辨率相空间表征及其智能化解释”(41974137);陕西省科技计划项目“非正交坐标系伪谱法地震波模拟方法研究”(2020JM-018)联合资助

摘　　要：波动方程的数值求解是地震勘探高精度成像和反演方法的核心。为保证精度,常规方法往往根据模型最小速度设置固定的空间采样间隔,易在高速区形成"过采样",计算存在一定冗余性。由于重力引起的岩石压实效应,波的传播速度由浅入深整体增大,梯形坐标变换更契合这种速度变化趋势。为此,提出梯形网格伪谱法地震波模拟方法。利用基于梯形坐标变换的梯形网格法对介质进行剖分,在浅层低速区采用细网格、深部高速区采用粗网格,可有效减少网格点数;同时,通过伪谱法求解坐标变换后的变系数声波方程,并引入完全匹配层消除人工边界虚假反射,以兼顾波场模拟的效率与精度。Marmousi模型的测试结果表明,相较于常规固定网格剖分,梯形网格剖分的网格数可减少约69%;相较于常规网格伪谱法和高阶有限差分法,梯形网格伪谱法的计算时间分别减少约58%和60%,同时具有较小的数值频散。可见,所提方法是一种高效、高精度地震波模拟方法。The numerical solution to wave equation is the computational engine of many high-precision imaging and inversion methods in seismic exploration.In order to ensure the accuracy,regular methods generally set a fixed spatial sampling interval according to the minimum velocity of the model.This always causes over sampling in high velocity layers,and the calculation is redundant to some extent.Under the compaction of rocks cause by gravity,the wave propagation speed usually increases along with depth.The trapezoid coordinate transformation can incorporate the general increasing trend of wave velocity.So,we propose seismic wave simulation using a trapezoid grid pseudo-spectral method.We use the trapezoid grid based on trapezoid coordinate transformation to divide medium.The fine grid is used in the shallow low velocity region and the coarse grid is used in the deep high velocity region,which can effectively reduce the number of sampling points.Meanwhile,considering the accuracy of wave field simulation,we use pseudo-spectral method to solve the acoustic wave equation with variable coefficients after coordinate transformation and use the perfectly matched layer to eliminate the fictitious return wave caused by artificial boundary.The test on Marmousi model shows that comparing with the general method,the number of grids generated by trapezoidal grid can reduce 69%,and comparing with the regular grid pseudo-spectral method and the high order finite difference method,the calculation time of trapezoidal grid pseudo-spectral method can respectively reduce 58%and 60%.Therefore,this method is an efficient and high-precision seismic wave simulation method.

关 键 词：梯形坐标变换 地震波模拟 伪谱法 声波方程 

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