| 基于逐减随机震源采样法的频率域二维黏滞声波方程全波形反演 |
| |
| 引用本文： | 冯晅,鲁晓满,刘财,周超,金泽龙,张明贺.基于逐减随机震源采样法的频率域二维黏滞声波方程全波形反演[J].吉林大学学报(地球科学版),2016,46(6):1865-1873. |
| |
| 作者姓名： | 冯晅 鲁晓满 刘财 周超 金泽龙 张明贺 |
| |
| 作者单位： | 吉林大学地球探测科学与技术学院, 长春 130026 |
| |
| 基金项目： | 国家自然科学基金重点项目(41430322)，国家重点基础研究发展计划(“973”计划)项目(2013CB429805)Supported by the State Key Program of National Natural Science of China(41430322)，the State Key Development Program for Basic Research of China(2013CB429805) |
| |
| 摘 要： | 全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构，具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而，巨大的计算量是阻碍其发展的一个瓶颈问题。为此，研究者们提出了震源编码技术来减少计算量，但是此方法在模型更新过程中会引进随机串扰噪声，降低反演结果准确性。所以，在保证计算精度的情况下，本文提出了采用逐减随机震源采样的方法来高效计算全波形反演问题。笔者将此方法应用于频率域二维黏滞声波波动方程全波形反演，开始了在频率域进行随机震源采样类方法的研究，计算过程中共使用了依次增大的8个频率段；并应用Overthrust模型来验证此类随机震源采样法的正确性。实验结果表明：基于逐减随机震源采样法的反演结果与实际Overthrust模型的拟合误差为0.065 65，而应用基于全部震源的全波形反演方法得到的反演结果与实际Overthrust模型的拟合误差为0.064 64，两者差别不大；但计算用时由740 min减少到291.2 min，即计算效率提高了2.54倍。为了更好地确定方法的有效性，将其应用于Marmousi模型进行试算。模型试算结果表明：基于逐减随机震源和基于全部震源得到的反演结果与实际Marmousi模型的拟合误差分别为0.080 12和0.078 97，相差不大；但计算用时由1 218.9 min减少到274.4 min，计算效率提高了4.44倍。综上，在保证反演精度的情况下，基于逐减随机震源采样法的频率域全波形反演方法大大减少了计算量，具有不可替代的计算优势，并且没有引进随机串扰噪声。
|
| 关 键 词： | 逐减随机震源采样法 频率域 二维黏滞声波方程 全波形反演 无记忆拟牛顿算法 L-BFGS |
| 收稿时间： | 2016-03-06 |
Frequency-Domain Full Waveform Inversion of 2D Viscous Acoustic Wave Equation Using Decreasing Random Shot Subsampling Method |
| |
| Feng Xuan,Lu Xiaoman,Liu Cai,Zhou Chao,Jin Zelong,Zhang Minghe.Frequency-Domain Full Waveform Inversion of 2D Viscous Acoustic Wave Equation Using Decreasing Random Shot Subsampling Method[J].Journal of Jilin Unviersity:Earth Science Edition,2016,46(6):1865-1873. |
| |
| Authors: | Feng Xuan Lu Xiaoman Liu Cai Zhou Chao Jin Zelong Zhang Minghe |
| |
| Institution: | College of GeoExploration Science and Technology, Jilin University, Changchun 130026, China |
| |
| Abstract: | The full waveform inversion (FWI)method utilizes kinematic and dynamic information of pre-stack seismic data to rebuild underground velocity structure.However, the large amount of calculation is a bottleneck problem for its development.Therefore,Researchers proposed the techniques of the source-encoding in order to reduce calculation amount,but this method would introduce random crosstalk in model updates.The paper uses decreasing random shot subsampling method to invert full waveform effectively in the case of guaranteeing the calculation accuracy.The paper applied the method into frequency-domain full waveform inversion of 2D viscous acoustic wave equation,and research the type of random shot sampling method in the frequency domain.We totally use eight gradually increased frequencies in FWI process,and apply overthrust model to demonstrate the validity of such random shot sampling method.Through the experiment,we can see that the fitting error between inversion results based on decreasing random shot subsampling method and the actual overthrust model is 0.065 65,and the fitting error between inversion results based on all the shots and the actual overthrust model is 0.064 64.Namely the inversion results of the two methods have little difference.However, the calculated time reduces from 740 min to 291.2 min.Namely the computational efficiency increased 2.54 times.In order to better determine the validity of the method,we apply it to Marmousi model.The fitting error between inversion results based on decreasing random shot subsampling method and the actual Marmousi model is 0.080 12,and the fitting error between inversion results based on all the shots and the actual Marmousi model is 0.078 97.Namely the inversion results of the two methods have little difference.However, the calculated time reduces from 1 218.9 min to 274.4 min.Namely the computational efficiency increased 4.44 times.In conclusion,under the condition of the inversion accuracy,the full waveform inversion based on the decreasing random shot sampling method in frequency domain greatly reduces the amount of calculation.Namely it has irreplaceable computing advantage.Therefore there is no introduction of random crosstalk noise. |
| |
| Keywords: | decreasing random shot subsampling method frequency-domain 2D viscous acoustic wave equation full waveform inversion memory-less quasi-newton algorithm L-BFGS |
| 本文献已被 CNKI 万方数据 等数据库收录！ |
| 点击此处可从《吉林大学学报(地球科学版)》浏览原始摘要信息 |
| 点击此处可从《吉林大学学报(地球科学版)》下载免费的PDF全文 |
