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

doi:10.6038/pg20170523

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

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

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

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-ping^1,2,3, LIU He-xiu¹, DAI Shi-xin¹, ZHAO Yan-lin¹, ANDREW Hursthouse¹

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

摘要/Abstract

摘要：

本文将二维时间空间域和频率空间域声波全波形速度反演方法分别应用到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.

中图分类号:

P631

廖建平, 刘和秀, 戴世鑫, 赵延林, ANDREW Hursthouse. 2017. 二维时间空间域和频率空间域声波全波形速度反演方法的对比研究. 地球物理学进展, 32(5): 2029-2034. doi: 10.6038/pg20170523.

LIAO Jian-ping, LIU He-xiu, DAI Shi-xin, ZHAO Yan-lin, ANDREW Hursthouse. 2017. Research on comparisons of 2D acoustic wave full waveform velocity inversion in time-space domain and frequency-space domain. Progress in Geophysics. 32(5): 2029-2034. doi: 10.6038/pg20170523.

参考文献

Bednar J B, Shin C, Pyun S. 2007. Comparison of waveform inversion, part 2:Phase approach[J]. Geophysical Prospecting, 55(4):465-475.
Bian A F, Yu W H, Zhou H W. 2010. Progress in the frequency-domain full waveform inversion method[J]. Progress in Geophys. (in Chinese), 25(3):982-993, doi:10.3969/j.issn.1004-2903.2010.03.037.
Bunks C, Saleck F M, Zaleski S, et al. 1995. Multiscale seismic waveform inversion[J]. Geophysics, 60(5):1457-1473.
Charl-Hyun J, Shin C, Suh J H. 1996. An optimal 9-point, finite-difference, frequency-space, 2-D scalar wave extrapolator[J]. Geophysics, 61(2):529-537.
Chen S C, Chen G X. 2016. Full waveform inversion of the second-order time integral wavefield[J]. Chinese Journal of Geophysics (in Chinese), 59(10):3765-3776, doi:10.6038/cjg20161021.
Cui Y F, Peng G X, Wu G C, et al. 2016. Application of full waveform inversion velocity model-building technology for the fractured-vuggy reservoir[J]. Chinese Journal of Geophysics (in Chinese), 59(7):2713-2725, doi:10.6038/cjg20160734.
Gauthier O, Virieux J, Tarantola A. 1986. Two-dimensional nonlinear inversion of seismic waveforms:Numerical results[J]. Geophysics, 51(7):1387-1403.
Gui S, Liu H, Zhang Y J. 2017. The acceleration strategy of a simplified hybrid domain full waveform inversion on multi-GPUs[J]. Chinese Journal of Geophysics (in Chinese), 60(2):665-677, doi:10.6038/cjg20170220.
Guo X B, Liu H, Shi Y. 2016. Time domain full waveform inversion based on frequency attenuation[J]. Chinese Journal of Geophysics (in Chinese), 59(10):3777-3787, doi:10.6038/cjg20161022.
Hu Y, Han L G, Xu Z, et al. 2017. Demodulation envelope multi-scale full waveform inversion based on precise seismic source function[J]. Chinese Journal of Geophysics (in Chinese), 60(3):1088-1105, doi:10.6038/cjg20170321.
Hustedt B, Operto S, Virieux J. 2004. Mixed-grid and staggered-grid finite-difference methods for frequency-domain acoustic wave modelling[J]. Geophys. J. Int., 157(3):1269-1296.
Liao J P, Liu H X, Wang H Z, et al. 2011a. Study on rapid highly accurate acoustic wave numerical simulation in frequency space domain[J]. Progress in Geophys. (in Chinese), 26(4):1359-1363, doi:10.3969/j.issn.1004-2903.2011.04.029.
Liao J P, Liu H X, Wang H Z, et al. 2011b. High resolution acoustic wave full waveform velocity inversion in frequency space domain-theoretical model[J]. Progress in Geophys. (in Chinese), 26(5):1690-1695, doi:10.3969/j.issn.1004-2903.2011.05.023.
Liu L, Liu H, Zhang H, et al. 2013. Full waveform inversion based on modified quasi-Newton equation[J]. Chinese Journal of Geophysics (in Chinese), 56(7):2447-2451, doi:10.6038/cjg20130730.
Liu Y Z, Wang G Y, Yang J Z, et al. 2015. Multi-parameter full-waveform inversion for VTI media based on Born sensitivity kernels[J]. Chinese Journal of Geophysics (in Chinese), 58(4):1305-1316, doi:10.6038/cjg20150418.
Luo J R, Wu R S, Gao J H. 2016. Local minima reduction of seismic envelope inversion[J]. Chinese Journal of Geophysics (in Chinese), 59(7):2510-2518, doi:10.6038/cjg20160716.
Mora P. 1987. Nonlinear two-dimensional elastic inversion of multioffset seismic data[J]. Geophysics, 52(9):1211-1228.
Mora P. 1989. Inversion=migration+tomography[J]. Geophysics, 54(12):1575-1586.
Operto S, Virieux J, Dessa J X, et al. 2006. Crustal seismic imaging from multifold ocean bottom seismometer data by frequency domain full waveform tomography:Application to the eastern Nankai trough[J]. Journal of Geophysical Research, 111(B9):B09306, doi:10.1029/2005JB003835.
Pratt R G. 1990. Inverse theory applied to multi-source cross-hole tomography. Part 2:Elastic wave-equation method[J]. Geophysical Prospecting, 38(3):311-329.
Pratt R G. 1999. Seismic waveform inversion in the frequency domain, Part 1:Theory and verification in a physical scale model[J]. Geophysics, 64(3):888-901.
Pratt R G, Shin C, Hick G J. 1998. Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion[J]. Geophysical Journal International, 133(2):341-362.
Pratt R G, Worthington M H. 1990. Inverse theory applied to multi-source cross-hole tomography. Part 1:Acoustic wave-equation method[J]. Geophysical Prospecting, 38(3):287-310.
Shin C, Min D J. 2006. Waveform inversion using a logarithmic wavefield[J]. Geophysics, 71(3):R31-R42.
Sirgue L, Pratt R G. 2004. Efficient waveform inversion and imaging:A strategy for selecting temporal frequencies[J]. Geophysics, 69(1):231-248.
Tarantola A. 1984. Inversion of seismic reflection data in the acoustic approximation[J]. Geophysics, 49(8):1259-1266.
Virieux J, Asnaashari A, Brossier R, et al. 2017. An introduction to full waveform inversion[A].//Grechka V, Wapenaar K eds. Encyclopedia of Exploration Geophysics[M]. Tulsa, OK:SEG, R1-1-R1-40.
Virieux J, Operto S. 2009. An overview of full-waveform inversion in exploration geophysics[J]. Geophysics, 74(6):WCC1-WCC26.
Wang B F, Wu R S, Chen X H, et al. 2016. Non-linear parameter estimation method based on T-matrix[J]. Chinese Journal of Geophysics (in Chinese), 59(6):2257-2265, doi:10.6038/cjg20160628.
Wei Z F, Gao H W, Zhang J F. 2014. Time-domain full waveform inversion based on an irregular-grid acoustic modeling method[J]. Chinese Journal of Geophysics (in Chinese), 57(2):586-594, doi:10.6038/cjg20140222.
Xu K, Wang M Y. 2001. Finite element inversion of the coefficients of acoustic equation in frequency domain[J]. Chinese Journal of Geophysics (in Chinese), 44(6):852-864, doi:10.3321/j.issn:0001-5733.2001.06.015.
Zhang W S, Luo J, Teng J W. 2015. Frequency multiscale full-waveform velocity inversion[J]. Chinese Journal of Geophysics (in Chinese), 58(1):216-228, doi:10.6038/cjg20150119.
卞爱飞, 於文辉, 周华伟. 2010. 频率域全波形反演方法研究进展[J]. 地球物理学进展, 25(3):982-993, doi:10.3969/j.issn.1004-2903.2010.03.037.
陈生昌, 陈国新. 2016. 时间二阶积分波场的全波形反演[J]. 地球物理学报, 59(10):3765-3776, doi:10.6038/cjg20161021.
崔永福, 彭更新, 吴国忱,等. 2016. 全波形反演在缝洞型储层速度建模中的应用[J]. 地球物理学报, 59(7):2713-2725, doi:10.6038/cjg20160734.
桂生, 刘洪, 张玉洁. 2017. 简化混合域全波形反演多GPU加速策略[J]. 地球物理学报, 60(2):665-677, doi:10.6038/cjg20170220.
郭雪豹, 刘洪, 石颖. 2016. 基于频域衰减的时域全波形反演[J]. 地球物理学报,59(10):3777-3787,doi:10.6038/cjg20161022.
胡勇, 韩立国, 许卓,等. 2017. 基于精确震源函数的解调包络多尺度全波形反演[J]. 地球物理学报, 60(3):1088-1105, doi:10.6038/cjg20170321.
廖建平, 刘和秀, 王华忠,等. 2011a. 快速高精度的频率空间域声波数值模拟方法研究[J]. 地球物理学进展, 26(4):1359-1363, doi:10.3969/j.issn.1004-2903.2011.04.029.
廖建平, 刘和秀, 王华忠,等. 2011b. 高分辨率的频率空间域声波全波形速度反演-理论模型[J]. 地球物理学进展, 26(5):1690-1695, doi:10.3969/j.issn.1004-2903.2011.05.023.
刘玉柱, 王光银, 杨积忠,等. 2015. 基于Born敏感核函数的VTI介质多参数全波形反演[J]. 地球物理学报, 58(4):1305-1316, doi:10.6038/cjg20150418.
刘璐, 刘洪, 张衡,等. 2013. 基于修正拟牛顿公式的全波形反演[J]. 地球物理学报, 56(7):2447-2451, doi:10.6038/cjg20130730.
罗静蕊, 吴如山, 高静怀. 2016. 地震包络反演对局部极小值的抑制特性[J]. 地球物理学报, 59(7):2510-2518, doi:10.6038/cjg20160716.
王本锋, 吴如山, 陈小宏,等. 2016. 基于T-matrix的非线性参数估计方法[J]. 地球物理学报, 59(6):2257-2265, doi:10.6038/cjg20160628.
魏哲枫, 高红伟, 张剑锋. 2014. 基于非规则网格声波正演的时间域全波形反演[J]. 地球物理学报, 57(2):586-594, doi:10.6038/cjg20140222.
许琨, 王妙月. 2001. 声波方程频率域有限元参数反演[J]. 地球物理学报, 44(6):852-864, doi:10.3321/j.issn:0001-5733.2001.06.015.
张文生, 罗嘉, 滕吉文. 2015. 频率多尺度全波形速度反演[J]. 地球物理学报, 58(1):216-228, doi:10.6038/cjg20150119.

[1]	屈念念,李家斌,张西君,刘俊. 基于重磁资料研究贵州省深部构造特征[J]. 地球物理学进展, 2019, 34(5): 1785-1793.
[2]	安鹏,曹丹平,赵宝银,杨晓利,张明. 基于LSTM循环神经网络的储层物性参数预测方法研究[J]. 地球物理学进展, 2019, 34(5): 1849-1858.
[3]	陈稳,薛国强,陈卫营,周楠楠,卢云飞. SOTEM多分量激电响应特性分析[J]. 地球物理学进展, 2019, 34(5): 1859-1865.
[4]	豆辉,徐逸鹤. 黏性逆时偏移成像研究进展[J]. 地球物理学进展, 2019, 34(5): 1866-1877.
[5]	毛庆辉,王鹏,曾隽. 水力压裂微地震事件定位方法综述[J]. 地球物理学进展, 2019, 34(5): 1878-1886.
[6]	尚新民,王延光,崔庆辉,金昌昆,赵翠霞,滕厚华,赵胜天. 准南缘山前带精细近地表建模技术研究[J]. 地球物理学进展, 2019, 34(5): 1887-1892.
[7]	石战战,夏艳晴,王元君,周怀来,庞溯,池跃龙. 基于L₁-L₁范数稀疏表示的共偏移距道集地震数据重建方法[J]. 地球物理学进展, 2019, 34(5): 1893-1899.
[8]	王迪,陈敏,王健伟,盛志超,严曙梅,马美媛,李东哲,郑欣. 基于T₂和毛管压力分类转换的渗透率计算新方法[J]. 地球物理学进展, 2019, 34(5): 1900-1909.
[9]	王海波,张伟,张宏,肖虎,高云路,顾剑华,唐美斌,韩峰. 高精度可控震源在深反射地震采集中的应用[J]. 地球物理学进展, 2019, 34(5): 1910-1916.
[10]	王静, 张军华, 王延光, 杨勇, 冯德永, 武刚, 黄德峰. 特征值相干理论诠释及效果比较[J]. 地球物理学进展, 2019, 34(5): 1917-1923.
[11]	王维红,包培楠,陈国飞,林弘喆. 基于最小平方的横波速度拟合及应用[J]. 地球物理学进展, 2019, 34(5): 1924-1929.
[12]	王亚,王孟华,田建章,王成泉,彭玲丽,冯小英,周兴海,张浩. 一种适用于深潜山及内幕的有效储层地震预测方法——以廊固凹陷杨税务潜山为例[J]. 地球物理学进展, 2019, 34(5): 1930-1937.
[13]	吴曲波,潘自强,陈金勇,Mubarak,Hani. 利用瞬态瑞雷面波法探测浅层玄武岩三维分布——以沙特Sabkhah Ad Dumathah地区为例[J]. 地球物理学进展, 2019, 34(5): 1938-1944.
[14]	严小丽,康慧敏,王光杰,何国丽,李红领. AMT方法在鳌山卫花岗岩地区深部地热构造勘探中的应用[J]. 地球物理学进展, 2019, 34(5): 1945-1953.
[15]	尹奇峰,潘冬明,郭全仕,何情. 基于快速模拟退火算法的井中微地震事件定位反演[J]. 地球物理学进展, 2019, 34(5): 1954-1961.

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

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

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

多维度评价

本文评价