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

全波形反演可以利用叠前地震波场的运动学和动力学信息重建地下速度结构，具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而，庞大的计算量和存储空间需求，限制了全波形反演的发展。在频率多尺度全波形反演中将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。数值模拟试验结果表明：此方法反演精度高，内存开销较小，能够显著提高计算效率，并且具有良好的抗噪能力，能够反演出具有速度扰动特性的介质。

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

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.