地球物理学进展 ›› 2017, Vol. 32 ›› Issue (5): 2029-2034.doi: 10.6038/pg20170523

• 应用地球物理学Ⅰ • 上一篇    下一篇

二维时间空间域和频率空间域声波全波形速度反演方法的对比研究

廖建平1,2,3, 刘和秀1, 戴世鑫1, 赵延林1, ANDREW Hursthouse1   

  1. 1. 湖南科技大学页岩气资源利用湖南省重点实验室, 湘潭 411201;
    2. 煤炭资源与安全开采国家重点实验室(中国矿业大学), 北京 100083;
    3. 湖南科技大学地质系, 湘潭 411201
  • 收稿日期:2017-02-22 修回日期:2017-08-18 出版日期:2017-10-20 发布日期:2017-10-20
  • 作者简介:廖建平,男,1976年生,博士,副教授,主要从事复杂介质中地震波传播及数值模拟、走时反演、二维和三维频率-空间域、时间-空间域声波和弹性波全波形反演方法研究.(E-mail:bobljp@163.com)
  • 基金资助:

    国家自然科学基金(41274126,51274097)、煤炭资源与安全开采国家重点实验室(中国矿业大学)开放基金(SKLCRSM11KFB01)、湖南省自然科学基金(12JJ6035)、中国博士后科学基金(2012M541116)、湖南省教育厅科研项目14C0462和页岩气资源利用湖南省重点实验室开放基金联合资助.

Research on comparisons of 2D acoustic wave full waveform velocity inversion in time-space domain and frequency-space domain

LIAO Jian-ping1,2,3, LIU He-xiu1, DAI Shi-xin1, ZHAO Yan-lin1, ANDREW Hursthouse1   

  1. 1. Hunan Provincial Key Laboratory of Shale Gas Resource Utilization, Hunan University of Science and Technology, Xiangtan 411201, China;
    2. State Key Laboratory of Coal Resources and Safe Mining, China University of Mining & Technology, Beijing 100083, China;
    3. Department of Geology of Hunan University of Science and Technology, Xiangtan 411201, China
  • Received:2017-02-22 Revised:2017-08-18 Online:2017-10-20 Published:2017-10-20

摘要:

本文将二维时间空间域和频率空间域声波全波形速度反演方法分别应用到Marmousi模型,进行数值试验.两种方法均采用相同的观测系统和其他的参数,理论模型的数值试验结果证实了:使用较多的计算集群的CPU进行二维频率空间域直接法声波全波形反演时,其加速有限(正演数值模拟的计算量主要用于稀疏矩阵的LU分解,炮点计算波场时为线性关系).二维时间空间域声波全波形反演计算时更灵活,多炮同时计算时,可以多倍提高其计算效率;二维声波全波形速度反演时,直接法求解频率空间域的计算速度远快于时间空间域,所需要的计算机内存也比时间空间域少.二维声波全波形速度反演时,相比较于时间空间域的方法,频率空间域直接法声波全波形反演具有计算速度快和节省计算机内存需求的优势.

Abstract:

The 2D time-space domain and frequency-space domain acoustic wave full waveform velocity inversion method (FWI) are applied to Marmousi model respectively. These two methods both use a same geometry and other parameters in the numerical experiments. It is confirmed that when using more computer nodes for direct FWI in frequency-space domain, its acceleration rate is limited (the forward modeling calculation is mainly used for the LU decomposition of the impedance matrix, shot point calculating the wavefield is linear relationship). 2D time-space domain FWI is more flexible, when more shots are calculated at the same time, it can be many times to improve calculation efficiency. The calculation speed of 2D frequenc-space domain direct FWI is much faster than that of the 2D time-space domain FWI, the memory requirements are much smaller than that of the 2D time-space domain FWI. Compared with the method of time-space domain FWI, frequency-space domain method of direct FWI has the advantages of high computing speed and saving computer memory requirements.

中图分类号: