L-BFGS和SPGL1在全波形反演中的联合应用.pdfVIP

·714· 中国地球物理2011 专题二十：地震波传播与成像探查l◆ (20)地震波传播与成像探查 L．BFGS和SPGLl在全波形反演中的联合应用 L—BFGSand theCombinationof SPGL1toFullWaveformInversion Apply 邓武兵+韩立国 李翔 张生强 Li DengWubing‘HanLiguoXiangZhangShengqiang 吉林大学地球探测科学与技术学院长春130026 全波形反演是基于地震波的偏微分方程的数值求解和局部最优化的过程，通过该过程来求解地下介质 物性参数重建参数模型。从数学角度讲，FWI是依赖于地震采集数据体和高精度计算求解的数据拟合过程。 通常我们把全波形反演表示成非线性最小二乘优化问题，如下式： C(Ad)=1／2《Ad十Ad)，其中△d=沈缸一d厶 (1) 其中，比m比h酣分别为观测数据、人工合成数据和二者的数据残差。FWI的目的就是为了通过数据拟 表现得非常好：相对于其它牛顿类的方法，它的收敛速度快，反演精度高。但是FWI是～个不适定问题， 常密集，所以在试验中表现出L．BFGS全波形反演的计算效率依然不是很理想。为了解决这一问题，本文 L—BFGS搜索方向的求解问题)，然后反curvelt变换回来。 Den 节约内存的关系而被忽略掉。SPGLl是针对最小二乘L1．范数正则化问题的解算法。是基于[E．Van Berg，M． Den 阵．向量积Ax、AVy上有出色的表现，A可以是显示矩阵也可以是显示函数。fE．Van 2008]进行了一系列的数值测试，证明了SPGLl算法具有很强的优势。[Xiang 算法来解决FWI中的优化问题，并取得了很好的效果。这证明了SPGLl适用于全波形反演。因此本文根据 的结合点在L．BFGS算法搜索方向的求解问题上，在L—BFGS算法中，搜索方向是通过(2)式得到的： dk=-HI戤 (2) 夙矾=一g。，将搜索方向求解转换为以下问题： (B尼)叫nlldll．s．t．118d一(一g)ll，s盯 (3) 口为误差水平。当，仃=0时，(肥)就转化成了(胎)问题。显而易见，∥越大，问题就变得越简 单，本文取盯为一个接近0的值。另外，本文对Pareto曲线进行采样获取了两个更新之间的信息并加以利 用，减少了矩阵一向量积的计算次数，大大提高了计算速率。这使得本文对L—BFGS的改进更加适合于大 新的～次迭代时会继续用到上次迭代的信息，这样能加快收敛速度。 结果证明了这种结合方式的有效性及有用性。