基于L-BFGS算法和同时激发震源的频率多尺度全波形反演

吉林大学学报(地球科学版)

基于L-BFGS算法和同时激发震源的频率多尺度全波形反演

张生强¹，刘春成²，韩立国¹，杨小椿²

1.吉林大学地球探测科学与技术学院，长春130026;
2.中国海洋石油研究总院，北京100027

收稿日期:2012-09-30 出版日期:2013-05-26 发布日期:2013-05-26
通讯作者: 韩立国（1961-），男，教授，博士生导师，主要从事地震数据处理和解释方面研究 E-mail:hanliguo@jlu.edu.cn
作者简介:张生强(1987-)，男，博士研究生，主要从事地震波成像和反演方面的研究，E-mail:zhangsq1987@126.com
基金资助:
国家科技重大专项项目(2011ZX05025-001-04)；国家自然科学基金项目(41074075)

Frequency Multi-Scale Full Waveform Inversion Based on L-BFGS Algorithm and Simultaneous Sources Approach

Zhang Shengqiang¹, Liu Chuncheng²，Han Liguo¹, Yang Xiaochun²

1.College of GeoExploration Science and Technology, Jilin University, Changchun130026, China;
2.CNOOC Research Institute，Beijing100027, China

Received:2012-09-30 Online:2013-05-26 Published:2013-05-26

摘要/Abstract

摘要：

全波形反演可以利用叠前地震波场的运动学和动力学信息重建地下速度结构，具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而，庞大的计算量和存储空间需求，限制了全波形反演的发展。在频率多尺度全波形反演中将LBFGS数值优化算法与同时激发震源技术相结合的方法来改善这一现状。首先，对Marmousi模型进行了速度反演：在计算过程中明显发现对计算机内存的占用减少，最终反演结果与实际Marmousi模型的拟合误差为0.095 9，较小；采用10个频带单炮震源正演384炮所需时间约为32 640 s，而采用同时激发震源（384炮）正演一次所需时间仅约为700 s。然后，基于高速楔形体模型进行了抗噪能力研究：原始含噪地震记录信噪比为11.147 3 dB；对反演得到的速度模型进行正演，其地震记录信噪比为22.251 8 dB。最后，基于逆冲断层模型进行了反演速度扰动能力研究，反演得到的最终模型很清晰，与具有速度扰动特性的实际模型非常接近，拟合误差仅为0.036 0。数值模拟试验结果表明：此方法反演精度高，内存开销较小，能够显著提高计算效率，并且具有良好的抗噪能力，能够反演出具有速度扰动特性的介质。

关键词: 全波形反演, 频率域, L-BFGS算法, 同时激发震源, 多尺度

Abstract:

Full waveform inversion (FWI) can reconstruct underground velocity while utilizing the kinematic and dynamic information of pre-stack seismic data, which could reveal detail information of the structure and lithology under complex geological background. However, huge calculating amount and storage space requirements confine the development of full waveform inversion. A combined method of L-BFGS algorithm and simultaneous sources technology is applied to ameliorate this problem during frequency multi-scale full waveform inversion. First we carry out velocity inversion based on the Marmousi model. It can be obviously found that the memory spending has been improved considerably in the process of calculation, and the fitting error between the final inversion result and true Marmousi model is as small as 0.095 9. Meanwhile, single shot forward modeling 384 times takes about 32 640 seconds when using 10 frequency bands, while simultaneous sources which contain 384 shots forward modeling one time needs only about 700 seconds. Then we conduct research on the anti-noise ability based on the high-speed wedge model. The signal to noise ratio of original seismic data is 11.147 3 dB, while the signal to noise ratio of the seismic data obtained from forward modeling based on inversion velocity model is 22.251 8 dB. Finally we conduct research on inversing the medium with velocity perturbation based on the thrust fault model. The final inversion model is very clear and very satisfied with the real model with velocity perturbation characteristics, and the fitting error between them is only 0.036 0. Numerical simulation experiment result indicates that the inversion precision of this method is higher enough, and the memory spending is less, which could significantly increase calculation efficiency. What is more, this method has good anti-noise ability, and could inverse the medium with velocity perturbation characteristics.

Key words: full waveform inversion, frequency-domain, L-BFGS algorithm, simultaneous sources, multi-scale

中图分类号:

P631.4

张生强，刘春成，韩立国，杨小椿. 基于L-BFGS算法和同时激发震源的频率多尺度全波形反演[J]. 吉林大学学报(地球科学版), 2013, 43(3): 1004-1012.

Zhang Shengqiang, Liu Chuncheng，Han Liguo, Yang Xiaochun. Frequency Multi-Scale Full Waveform Inversion Based on L-BFGS Algorithm and Simultaneous Sources Approach[J]. Journal of Jilin University(Earth Science Edition), 2013, 43(3): 1004-1012.

参考文献

相关文章 15

[1]	曾昭发, 霍祉君, 李文奔, 李静, 赵雪宇, 何荣钦. 任意各向异性介质三维有限元航空电磁响应模拟[J]. 吉林大学学报(地球科学版), 2018, 48(2): 433-444.
[2]	刘新彤, 刘四新, 孟旭, 傅磊. 低频缺失下跨孔雷达包络波形反演[J]. 吉林大学学报(地球科学版), 2018, 48(2): 474-482.
[3]	崔利凯, 孙建孟, 闫伟超, 高银山, 王洪君, 宋丽媛. 基于多分辨率图像融合的多尺度多组分数字岩心构建[J]. 吉林大学学报(地球科学版), 2017, 47(6): 1904-1912.
[4]	李光, 渠晓东, 黄玲, 方广有. 基于磁偶极子的频率域电磁系统几何误差分析[J]. 吉林大学学报(地球科学版), 2017, 47(4): 1255-1267.
[5]	李振苓, 沈金松, 李曦宁, 王磊, 淡伟宁, 郭森, 朱忠民, 于仁江. 用形态学滤波从电导率图像中提取缝洞孔隙度谱[J]. 吉林大学学报(地球科学版), 2017, 47(4): 1295-1307.
[6]	冯晅, 鲁晓满, 刘财, 周超, 金泽龙, 张明贺. 基于逐减随机震源采样法的频率域二维黏滞声波方程全波形反演[J]. 吉林大学学报(地球科学版), 2016, 46(6): 1865-1873.
[7]	张广智, 孙昌路, 潘新朋, 张志明, 姜岚杰, 印兴耀. 变密度声波全波形反演中密度影响因素及反演策略[J]. 吉林大学学报(地球科学版), 2016, 46(5): 1550-1560.
[8]	刘财, 杨宝俊, 冯晅, 单玄龙, 田有, 刘洋, 鹿琪, 刘才华, 杨冬, 王世煜. 论油气资源的多元勘探[J]. 吉林大学学报(地球科学版), 2016, 46(4): 1208-1220.
[9]	赵建如, 初凤友, 金路, 杨克红, 葛倩. 珠江口西部海域表层沉积物重金属元素多尺度空间变化特征[J]. 吉林大学学报(地球科学版), 2015, 45(6): 1772-1780.
[10]	高成, 孙建国. 不同域的局部平面波分解应用与对比[J]. 吉林大学学报(地球科学版), 2015, 45(5): 1523-1529.
[11]	于德武, 龚胜平. 对迭代法位场向下延拓方法的剖析[J]. 吉林大学学报(地球科学版), 2015, 45(3): 934-940.
[12]	王卫平, 曾昭发, 李静, 吴成平. 频率域航空电磁法地形影响和校正方法[J]. 吉林大学学报(地球科学版), 2015, 45(3): 941-951.
[13]	郭振波，李振春. 起伏地表条件下频率-空间域声波介质正演模拟[J]. 吉林大学学报(地球科学版), 2014, 44(2): 683-693.
[14]	汤文武，柳建新，童孝忠. 电导率连续变化的线源FCSEM有限元正演模拟[J]. 吉林大学学报(地球科学版), 2013, 43(5): 1646-1654.
[15]	徐志锋,胡文宝. 层状大地频率域长导线源激发的电磁场[J]. 吉林大学学报(地球科学版), 2013, 43(1): 275-281.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed