波形反演中频率点选取的优化方法.doc

波形反演中选取的优化方法建征黄慧明大学科学与技术学院，湖北72） 摘 要：频域全波形反演是一种能够获得高分辨率反演结果的反演方法，但在模型较大时反演的计算用时非常巨大，这就严重阻碍了其在实际反演中的应用此问题，利用线性递增速度模型中地震波的传播规律，提出一种优化的频率选取方法，达到有效提高计算效率的目的表明，化的频率点选取在模型较大时，相比于常规选取方法保证演精度的情况下，具有更高的效率。速度模型分类号：P 631.4 标识码： An Optimization Strategy for Selecting Frequencies in Waveform Inversion Yuan Jianzheng, Huang Huiming (Wuhan University, School of Physics and Technology, Wuhan 430072, China) Abstract:Frequency domain full-waveform inversion is one inversion method that can get high resolution inversion result, but the computational cost is a major obstacle to estimating the velocity model for large-scale problems, and seriously hiders its application in the actual waveform inversion. To overcome this problem, this paper proposed an optimization strategy for selecting frequencies based on the propagation of seismic waves in linear incremental velocity model, and the strategy can achieve the purpose of improving the calculation efficiency. The experiment results show that the optimization strategy is more effective than the conventional method when the model is large. Keywords: frequencies selecting; wavenumber coverage; linear increasing velocity; waveform inversion 0 引言 波形反演在对地下介质中地震波的传播过程进行数值建模的基础上，通过使观测数据和模拟数据的残差最小来获得一个对地下介质的定量估计模型[1,2]。目前，对于时域和频域中的全波形反演方法已经有了非常深入的研究[3,4]，就频域全波形反演而言依然面临着反容易陷入局部极小和反演大模型时用时过长两个主要问题。容易陷入局部极小的问题，可以通过提供一个比较接近真实地下结构的初始速度模型多形反演选取适当的反演策略等方法来解决。地震数据在频中存在冗余5]，因此可以使用所有频率的一个子集来表示信息从低频率点到高频率点的频率迭代反演来说，其计算用时正比于演的频率点，因此如何高地去除震数据集中的冗余频率点，也即如何地选择最好频率点组合就成了决第二个问题也即此时计算效率的关键。 速度模型中震波的传播规律，而后此规律对方法进行修正， 1 线性递增速度模型结构(x, z) (其中x为水平方向坐标，z为竖直方向坐标也即深度方向坐标)可以表示为 (1) 对于从炮(0, 0)传播到接收点(s, h)处的地震波而言，≠0时，由费原理可知传播路径是光滑且一阶可导的，传播用时为T可 (2) 由费马原理(射线原理)可知[7]，真实的传播路径一定能使得式得极小值式表示的泛函解可得 (3) 其中 (4) (5) 式(3)说明了在炮点和检波器点不在同一竖直线上时地震波的传播路径的具体表达形式可以看出其路径一个圆对于0时的情形，使用类似的方法可得时的传播路径为 (6) 式(6)说明了在炮点和检波器点在同一竖直线上时地震波沿直线传播。 2 频率点选取的优化方法 .1 连续介质中反射点的波数 由前面的讨论，及线性递增速度模型的对称性可知，也即(0, 0)传播到R(s, 0)地震波，在深度为0处的反射点的坐标为0(s/2, z0)，记在点0(s/2, z0)处的地震波与竖直方向的夹角为θ由式3)可得 (