袁雨欣^1,2,3,4,,
胡婷^1,2,3,4,
王之洋^1,2,3,
郭鹏⁵,
刘洪^1,2,3,4,,
1. 中国科学院地质与地球物理研究所, 北京 100029
2. 中国科学院地球科学研究院, 北京 100029
3. 中国科学院油气资源研究院重点实验室, 北京 100029
4. 中国科学院大学, 北京 100049
5. 中国地震局兰州地震研究所, 兰州 730000

基金项目: 国家自然科学基金项目（41630319）、中国博士后科学基金（2016M591244）、国家重点研发计划深地专项项目（2016YFC0601101）、中国石油集团"弹性波地震成像技术合作研发课题"联合资助


详细信息

作者简介: 袁雨欣, 男, 1994年生, 在读博士研究生, 主要从事地震波数值模拟研究.E-mail:yuanyuxin15@mails.ucas.ac.cn

通讯作者: 刘洪, 男, 1959年生, 研究员, 主要从事地震波成像和油储地球物理研究.E-mail:liuhong@mail.iggcas.ac.cn

中图分类号: P631

收稿日期:2017-08-02
修回日期:2018-06-27
上线日期:2018-08-05

Solving second-order decoupled elastic wave equation using low-rank decomposition and low-rank finite differences

YUAN YuXin^1,2,3,4,,
HU Ting^1,2,3,4,
WANG ZhiYang^1,2,3,
GUO Peng⁵,
LIU Hong^1,2,3,4,,
1. Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
2. Institute of Earth Science, Chinese Academy of Sciences, Beijing 100029, China
3. Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
4. University of Chinese Academy of Sciences, Beijing 100049, China
5. Lanzhou Institute of Seismology, China Earthquake Administration, Lanzhou 730000, China


More Information

Corresponding author: LIU Hong,E-mail:liuhong@mail.iggcas.ac.cn

MSC: P631

--> Received Date: 02 August 2017
Revised Date: 27 June 2018
Available Online: 05 August 2018


摘要
摘要:时间域的波场延拓方法在本质上都可以归结为对一个空间-波数域算子的近似.本文基于一阶波数-空间混合域象征，提出一种新的方法求解解耦的二阶位移弹性波方程.该方法采用交错网格，连续使用两次一阶前向和后向拟微分算子，推导得到了解耦的二阶位移弹性波方程的波场延拓算子.由于该混合域象征在伪谱算子的基础上增加了一个依赖于速度模型的补偿项，可以补偿由于采用二阶中心差分计算时间微分项带来的误差，有效地减少模拟结果的数值频散，提高模拟精度.然而，在非均匀介质中，直接计算该二阶的波场延拓算子，每一个时间步上需要做N次快速傅里叶逆变换，其中N是总的网格点数.为了减少计算量，提出了交错网格低秩分解方法；针对常规有限差分数值频散问题，本文将交错网格低秩方法与有限差分法结合，提出了交错网格低秩有限差分法.数值结果表明，交错网格低秩方法和交错网格低秩有限差分法具有较高的精度，对于复杂介质的地震波数值模拟和偏移成像具有重要的价值.
关键词: 波场延拓/
波数-空间混合域象征/
交错网格低秩分解方法/
交错网格低秩有限差分法

Abstract:Seismic wavefield extrapolation is an important part of seismic imaging and full waveform inversion. Theoretically, the problem of wave extrapolation in the time domain can be reduced into analyzing numerical approximation to the space-wavenumber mixed-domain symbol. In this paper, we propose a novel method to solve decoupled second-order elastic wave equation. We successively employ backward and forward first-order pseudo differential operators to derive a wave extrapolation operator for second-order decoupled elastic wave propagation. The mixed-domain symbol incorporates the accurate spectral evaluation of spatial derivatives and the time-marching adjustment to ensure that the solution is exact for homogeneous wave propagation for time steps of an arbitrarily large size. Considering its straightforward implementation in heterogeneous media, it needs to do N times inverse fast Fourier transform (FFT) every time step, here N is the total size of the model grid. In order to reduce computational cost, we propose a staggered-grid low-rank(SGL) method, which can be applied to large time step wave extrapolation. We also propose a staggered-grid low-rank finite-difference (SGLFD) method by combining the staggered low-rank method and finite-difference to reduce numerical dispersion. The 2D numerical experiments demonstrate that these two methods can improve the accuracy of modeling results compared with the ordinary staggered-grid finite difference method.
Key words:Wavefield extrapolation/
Space-wavenumber mixed-domain symbol/
Staggered-grid low-rank decomposition/
Staggered-grid low-rank finite difference



PDF全文下载地址:

http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_14634