CN102918521A - 使用时变滤波器的全波场反演 - Google Patents

Abstract

本发明涉及一种在通过局部目标函数优化（64）执行地震数据（65）多尺度反演时，减小对起始模型的精度要求的改进方法。通过引入低通滤波器到目标函数（61），然后减小从一个尺度到另一个尺度被滤除的高频数据的量，实现不同尺度的反演。而且，滤波器被设计为是时变的，其中滤波器的低通截止频率随被滤除的地震数据的传播时间的增加而减小（62）。滤波器可如下设计：使用消除局部最小值的普拉特（Pratt）标准，仅针对源和接收器位置而非传播时间，执行周期和传播时间误差的平均（或其他统计度量）（63）。

Description

使用时变滤波器的全波场反演

相关申请的交叉参考

本申请要求2010年3月29日申请的美国临时专利申请61/318561的权益，其标题为FULL WAVEFIELD INVERSION USING TIMEVARYING FILTERS，其全部内容通过引用包括在此。

技术领域

本发明一般涉及地震数据的数值转换从而推导传播介质的弹性参数的领域。更具体地，本发明是用于在反演，如地震数据反演中执行局部目标函数优化时降低对起始模型的精度要求的方法。

背景技术

反演（例如参考Tarantola，1984）试图发现最架解释观察的数据的模型。最小化测量模拟数据和观察数据之间差的目标函数值的局部反演方法通常是解决具有大量自由参数的模型的反演问题的唯一实用方法。这些局部方法要求针对要反演的模型的初始猜测。这些局部方法迭代更新模型，从而通过在基于目标函数梯度的方向上搜索当前模型的扰动使其更接近真解。遗憾的是，目标函数通常具有许多最小值，而不仅仅是一个对应于解模型的最小值。这些其他最小值被称为局部最小值，而对应于所需解的最小值被称为全局最小值。如果反演的起始模型太接近对应于一个这些局部最小值的模型，则局部反演方法在局部最小值附近卡壳，并且永远也不能从该值迭代到全局最小值。因此，无论如何努力迭代，都会产生错误解。

该局部最小值问题可通过对改变的目标函数执行第一迭代反演解决，改变的目标函数具有较少的局部最小值但在所需解的位置附近有全局最小值。对改变的目标函数迭代的结果应该产生更接近所需解的模型。然后该更精确模型被用作针对原始目标函数上的迭代的初始模型。因为这个新初始模型接近原始目标函数的全局最小值，所以对原始目标函数的迭代现在应产生精确解。在改变的目标函数上迭代的这个技术通常被称为多分辨率，或多网格，或多尺度反演，这将在下面进一步讨论。

存在大量已知的反演方法。这些方法落入两个类别，迭代反演和非迭代反演。下面是这两类中每一个通常表示的定义：

·非迭代反演——通过采用某些简单背景模型和基于输入数据更新模型实现的反演。该方法不使用更新的模型作为另一个反演步骤的输入。对于地震数据的情况，这些方法通常被称为成像、迁移、衍射层析（diffraction tomography）、线性反演或波恩（Born）反演。

·迭代反演——反演涉及重复改进地下特性模型，以便找到能满意地解释观察的数据的模型。如果反演收敛，则最终模型更好地解释观察的数据，并将更接近地逼近实际地下特性。迭代反演通常产生比非迭代反演更精确的模型，但计算成本更高。

波反演意味着任何基于波模拟器的反演，如声学或地震反演。在波反演中最常用的迭代方法是目标函数优化。相对模型M，目标函数优化涉及目标函数S(M)的值的迭代最小化，该目标函数是计算的和观察的数据之间失配的度量（这有时也称为成本函数）。计算的数据是利用被编程为使用物理学决定的当前模型表示的介质中源信号传播的计算机模拟的。该模拟计算可通过几种数值方法实现，包括但不限于有限差、有限元或光线跟踪。在Tarantola的文献【Tarantola，1984】之后，最常用的目标函数是最小二乘目标函数：

S(M)=(u(M)-d)^TC^-1(u(M)-d)，（1）

其中T表示矢量转置算符，且：

M=是N个参数的矢量[m₁,m₂,...,m_N]^T的模型，

d=测量的数据矢量（相对源、接收器和时间采样的），

u(M)=针对模型M的模拟数据矢量（相对源、接收器和时间采样的），

C=协方差矩阵。

目标函数优化方法是局部的或全局的【Fallat等人，1999】。全局方法仅涉及计算模型群体{M₁,M₂,M₃,...}的目标函数S(M)并从该模型群体中选择近似最小化S(M)的包括一个或更多模型的集合。如果需要进一步的改进，则新选择的模型集合可用作生成可再次相对目标函数S(M)测试的新群体模型的基础。全局方法比局部方法更可能对正确的解收敛，但应用于具有许多模型参数的多尺度反演问题成本太高。已知的全局反演方法包括蒙特卡洛（Montre Carlo）、模拟退火、基因和进化算法。

局部目标函数优化涉及：

算法1：用局部目标函数优化更新模型的算法。

局部反演方法比全局方法更有效，并因此是用于大尺度反演问题的唯一实用方法。通用的局部目标函数反演方法包括最陡下降、共轭梯度和牛顿方法。

应该注意，

的计算，算法1第二步要求针对N个模型参数m_i中每个参数计算S(M)的导数。当N非常大（约超过1000），如果必须为每个模型参数执行该计算，则该计算极其耗时。幸运的是，伴随法可用来一次为所有模型参数有效地执行该计算【Tarantola，1984】。用于最小二乘目标函数和网格状模型参数化（gridded modelparameterization）的伴随法由下面的算法总结：

算法2：用伴随法计算网格状模型最小二乘成本函数梯度的算法。

局部目标函数优化通常比全局目标函数优化便宜得多，但要求更精确的起始模型。要求该更精确的起始模型，是因为目标函数通常具有许多最小值，且局部优化方法通常发现这些最小值中最接近的一个。对应于真实模型的最小值是所谓的全局最小值，且所有其他最小值被称为局部最小值。如果起始模型不是最接近全局最小值，则局部优化技术可能产生不精确的反演模型，其对应于最接近的局部最小值。这在图1中示出，其中目标是针对具有两个参数m1和m2的模型M的反演。虚周线110显示目标函数的值作为参数m1和m2的函数。全局最小值120是由实线黑圆标记的，且两个局部最小值130和140由灰色填充圆示出。反演从初始模型M⁽⁰⁾(150)开始，然后由局部优化迭代一次模型M⁽¹⁾，直到模型M⁽³⁾(160)。无论再尝试多少次局部优化的迭代，反演的模型都仅更接近M⁽³⁾附近的局部最小值130，而非接近全局最小值120。

已经提出了几种试图克服该局部最小值问题的方法。如上所述，许多这类方法在反演的早期迭代中都涉及对改变的目标函数的迭代。选择这个改变的目标函数从而具有较少的局部最小值，但具有原始目标函数的全局最小值附近的全局最小值。通过该方法，早期迭代将产生这样的模型，虽然其不精确，但显著更接近原始目标函数的全局最小值。图2示出对应于图1的局部优化，但使用具有较少局部最小值的改变的目标函数。改变的目标函数具有全局最小值210（实黑圆），其靠近原始目标函数（交叉阴影圆）的全局最小值220，但与其不在相同位置。从初始模型M⁽³⁾(230)开始，该模型与图1中使用的初始模型相同，使用改变的目标函数的两次迭代导致模型M⁽²⁾(240)。该模型M⁽²⁾可用作初始模型，以便进一步迭代，但现在使用初始目标函数。这在图3中示出，其中来自图2的迭代2模型（以灰色示出），再编号为310并用作起始模型。现在迭代对靠近全局最小值220而非靠近图1中局部最小值的模型M⁽⁴⁾(320)收敛。因为起始模型比原始开始模型更精确，所以反演现在迭代到正确解。

通常在改变原始目标函数时，改变的目标函数中局部最小值的数目与原始函数的全局最小值和改变的目标函数的全局最小值之间的距离逆相关。因此，从具有较少数目的局部最小值和最不精确局部最小值的目标函数开始，然后演进到具有增加数目的局部最小值和增加精度的全局最小值的目标函数，然后以对原始目标函数迭代结束，对改变的目标函数序列进行迭代是有利的。对具有几个局部最小值的改变的目标函数执行初始迭代的方法通常被称为多尺度或多网格方法，且该技术的流程图在图4中示出。

过程从步骤410开始，选择原始目标函数的变化从而优化。该改变的目标函数取决于要拟合420的数据，并在步骤430迭代，直到改变的目标函数在步骤440被发现被充分最小化。（该值小于所选最大值或满足另一个停止条件）。当发生该情形时，在步骤450确定是否当前反演模型充分最小化原始目标函数。如果没有，则过程返回到步骤410且要么选择新改变的目标函数或原始目标函数来优化。最终，过程结束（460）。

在解决地震全波场反演（“FWI”）的局部最小值问题的文献中已经提出了两个改变的目标函数：

·高阻滤波器：Bunks的文献（Bunks等人，1995年）描述了通过对测量的数据和用于计算模拟的地震数据的源特征（source signature）都应用时间不变高阻滤波器（有时被称为低通滤波器，意味着使低于其截止频率的频率通过并摒弃截止频率以上频率的滤波器），改变最小二乘目标函数。然后，随着反演进行迭代，这些滤波器的高阻，即截止频率增加，其中没有滤波器应用于最终迭代（没有滤波器对应于原始目标函数）。已知如何设计这样的滤波器；参看Press等人的Numerical Recipes in FORTRAN，The Art of Scientifc Computing，Cambridge University Press（1992）。这也可从下面的资源获得，如Seismic Un*x（参看：http:／／www cwp mines edu/cwpcodes/）。

·剥层法：Maharramov的文献（Maharramov等人，2007）教导了反演的初始迭代应局部化在浅层，且随迭代进行而延伸到深度范围。相应地，当仅浅深度被反演时，数据中仅较短时间被反演，因为浅模型可仅预测较短时间部分的数据。

一般地，如果起始模型足够精确，从而预测任何传播模式的传播时间到该模式的半个周期内，则FWI收敛到全局最小值。这也可称为普拉特（Pratt）的标准。（“实际上，对于地震波形反演，这意味着许多波形能量必须被预测在观察波形的半个波长内（通过初始模型）；如果不是，则当预测的波形匹配观察波形内的错误循环时，则获得最小失配模型。”——Pratt，1999）。

本发明是当执行局部目标函数优化时降低对起始模型精度要求的改进方法。

发明内容

本发明方法适用于任何基于波模拟器的反演，如声学或地震反演。在其一个方面中，本发明是用于改变目标函数的特定方法，对于起始模型中的给定不精确性，该目标函数减少发现全局最小值所需的迭代数目。减少迭代数目将相应地减少成本和计算时间。改变包括引入时变滤波器到目标函数中。选择该滤波器以便测量的数据和计算的数据之间传播时间差的某统计度量小于数据主周期的某个分数（通常四分之一）的占优周期。这意味着滤波器是高阻滤波器。进一步选择该滤波器以便该滤波器的高切断频率随传播时间的增加而减小，使得其为时变滤波器。

参考图6的流程图，在一个实施例中，本发明是反演测量的地震数据（65）从而推断地下区域的物理特性模型的方法，包括通过用目标函数的局部最小值优化（64）执行测量的地震数据的迭代、多尺度反演而连续更新模型，该目标函数计算模型模拟的地震数据和测量的地震数据之间的失配（61），其中变化的低通滤波器被用于通过过滤失配计算中测量的和模拟的地震数据（62），将目标函数连续地从一个尺度变到另一个尺度，所述滤波器是时变的，其中滤波器的低通截止频率随地震数据的传播时间改变，所述地震数据在反演的一个或更多尺度上被滤波。滤波器可以用Pratt的标准设计，用于减小局部最小值的数目，针对源和接收器位置而非传播时间修改从而包括数据的统计度量（如平均值），使得滤波器是时变的（63）。多尺度反演的最后阶段优选使用未变的目标函数，且因此其更有效地被提供更精确的起始模型以帮助其收敛到全局最小值，导致优化的物理特性模型（66）。

如同任何数据反演一样，该过程在实际应用中是高度自动化的，即，借助根据本公开内容编程的计算机执行。

附图说明

本发明及其优势可通过参考下面的详细描述和附图更好地被理解，其中：

图1是收敛到局部最小值的反演的示意图；

图2是对应于图1的局部优化的示意图，但其使用具有较少局部最小值的已改变的目标函数；

图3示出利用来自图1的原始目标函数的局部优化，其使用来自图2的迭代2模型（灰色示出）作为起始模型；

图4是示出多尺度优化中的基本步骤的流程图；

图5是图解说明测量的数据和模拟的数据之间传播时间误差t_error，和地震道的瞬时周期T的示意图；以及

图6是示出本发明一个实施例中基本步骤的流程图。

下面结合示例性实施例描述本发明。然而，一定程度上，下面的详细描述是针对本发明特定实施例或特殊使用的，这仅是为了说明的目的，而不能解读为对本发明范围的限制。相反，本发明涵盖权利要求限定的本发明范围内的所有替换、修改和等价物。

具体实施方式

数学上，Pratt的标准可表述为：

$\frac{\underset{s,r,t}{\max }\left|t_{\mathrm{error}}(M,s,r,t)\right|}{\underset{s,r,t}{\max T(s,r,t)}}\≤\frac{1}{2},---\left(2\right)$

其中t_error是测量的数据和模拟的数据之间的传播时间误差，且T是测量的数据的瞬时周期，如图5中所示。在图5中，测量的数据和模拟的数据之间的传播时间误差被指示为t_error，且地震道的瞬时周期由T指示。传播时间误差是对齐测量的数据和模拟的数据所需的时移量。瞬时周期是数据的相似相位（如，峰或槽）之间的时间。（传播时间意味着从地震源发生地震波开始直到地震波在接收器被记录逝去的时间。）

传播时间误差和瞬时周期都是源s、接收器r和传播时间t的函数。此外，t_error是当前模型M的精度的函数。实际上，等式2可比必需的约束更大，以确保FWI的收敛。特别地，可使用比最大值（如，平均值）不严格的统计度量，或非1/2的目标可用在不等式的右侧。因此，实际上，等式2可由下式取代：

$\frac{\underset{s,r,t}{\max }\left|t_{\mathrm{error}}(M,s,r,t)\right|}{\underset{s,r,t}{\mathrm{statT}(s,r,t)}}\≤\frac{\mathrm{\α}}{2}---\left(3\right)$

其中stat是某个统计量，如平均值、众数或均方差，且α是近似等于一的常数。结果对统计的选择非常灵敏是没有预计到的。

通过利用高阻滤波器以增加T(s,r,t)，因此允许t_error的较大值，目标函数的Bunk的变化遵从类似于Pratt的标准的推理逻辑。可以理解，局部最小值的主要原因是跳周（cycle skipping），且较长周期导致这个可能性较低。理论上，t_error可减小，而非限制数据到较低频率；然而，实现这一点的唯一方式是具有更精确的起始模型，这非常困难并可能是不可能的。而且，FWI的目标是产生更精确的模型，因此要求非常精确的起始模型减小了FWI的值。另一方面，Maharramov的剥层法通过仅反演浅模型来避免大的传播时间误差，该浅模型仅传播具有少的传播时间的模式。因为传播时间误差通常随传播时间增加，所以限制反演到较短传播时间将t_error保持在拇指规则内。

在本发明中，提出第等式3的“拇指规则”的替换，这导致确保收敛的新策略。该拇指规则的替换规则如下：

$\frac{\underset{s,r}{\max }\left|t_{\mathrm{error}}(M,s,r,t)\right|}{\underset{s,r}{\mathrm{statT}(s,r,t)}}\≤\frac{\mathrm{\α}}{2}---\left(4\right)$

该拇指规则不同于等式3，因为统计分析不再对时间执行。在应用统计计算后，等式4左侧的分子和分母不是源位置和接收器位置的函数。等式4等价于：

$\hat{T}\left(t\right)\≥\frac{2}{\mathrm{\α}}{\hat{t}}_{\mathrm{error}}(M,t)---\left(5\right)$

其中 $\hat{T}\left(t\right)=\underset{s,r}{\mathrm{stat}}T(s,r,t),$ 且 ${\hat{t}}_{\mathrm{error}}(M,t)=\underset{s,t}{\mathrm{stat}}\left|t_{\mathrm{error}}(M,s,r,t)\right|.$

实际上，

是传播时间t的增函数，因为传播时间误差倾向于随波传播通过不精确的模型而累积。等式5表明针对多尺度反演的最优策略会改变测量的数据和模拟的数据，以便地震数据的平均瞬时周期以类似于

的方式随传播时间增加。这可通过对测量的数据应用时变高阻滤波器实现。这个时变滤波器的高截止频率应随增加的传播时间而减小。该方法相对Bunk技术的优点是更多信息（即，在少的传播时间的较高频率信息）会用于反演的早期迭代，因此更好地约束反演的模型，导致更快的收敛。我们的建议相对Maharramov的剥层的优点是更多数据（即，传播通过模型较深部分的地震模式）用于早期迭代，导致更快的收敛。

函数

取决于传播时间误差是如何测量的和什么统计度量被应用于这些传播时间误差测量两者。然而，可以预期，这个提出的多尺度反演测量将对

不够灵敏。实际上，不是努力从数据测量替换策略是为

假设简单的函数形式，如线性函数

其中M₀是初始模型。这个假定的函数形式然后被用于设计满足等式5的时变高阻滤波器，且尝试使用该滤波器的反演。如果反演不收敛，则增加β的值，并可尝试具有更保守估计的反演。继续这个测试直到发现为初始模型产生收敛的β。

在发现对于这个初始模型，M₀收敛的β之后，迭代然后将产生比初始模型更精确的当前反演的模型。这个增加的精度表明β应该随迭代的进行减小。β的这个减小意味着通过更高频率的相应时变滤波器。反演利用通过越来越多高频率的时变滤波器进行，直到数据中所有频率都用于反演的最终迭代中。

为了使反演中的时变滤波器实用，能够以与时变滤波器一致的方式计算伴随梯度（算法2）是必要的。完成这一点的数学上最直接方式是用等式1中的逆协方差矩阵C^-1实现时变滤波器。为了完成这一点，选择逆协方差矩阵C^-1为非对角线的（在时间维中），其中元素等于时变滤波器系数的时间表示。因为滤波器会是时变的，所以滤波器系数随时间改变。例子1示出示例性的C^-1的子矩阵，其对应于实现时变滤波器的特殊源和接收器。该子矩阵的第一行为零，预期对角线上为一个1。

例子1

这意味着该特殊时变滤波器在时间零不执行滤波。在子矩阵中，非对角线元素随行增加而增加，这意味着该时变滤波器的高截止频率随时间增加而减小。注意，传播时间误差可被视作源或接收器的函数。在该情形中，该方法可以以比简单时变滤波器更一般的方式被应用。例如：

1.传播时间误差通常不仅是传播时间的函数，而且通常是源和接收器之间偏移（offset）的函数。如果是该情形，则为不同源接收器偏移使用不同时变滤波器是有利的。随着迭代的进行，滤波器由于模型精度增加而得到弛豫（relaxed）。

2.传播时间误差通常是源和/或接收器域中数据的时间倾斜（timedip）的函数。这会发生是因为陡峭的时间倾斜对应于主要在水平方向上传播的波，这些波对初始模型的精度更灵敏。在该情形中，时变倾斜滤波器（如，除去具有高时间倾斜的地震事件的频率-波数滤波器）可取代时变时间滤波器使用。随着迭代的进行，倾斜滤波器将随模型精度增加而弛豫。

在任何情形中，滤波器（如空间变化和时变滤波器、时变倾斜滤波器等）可在协方差矩阵C^-1中实现，如上所述的，且然后如上所述继续梯度计算。

优选方法可能是如果可为当前模型M估计

则时变滤波器可被设计成与等式5一致。否则，使用β的初始猜值，估计是线性函数βt是合理的。再次，时变滤波器应设计被成与β的估值和等式5都一致。如果这个反演收敛，则过程完成。如果该反演不收敛，则增加β，且尝试另一反演。增加β的这个过程继续，直到实现收敛。

实际上，表示时变滤波器的矩阵C^-1是非常稀疏的，且因此其应用于算法2的步骤3中的数据残余是通过应用时变滤波器算子最佳实现的，而非通过矩阵乘法。实际上，这个反演方法不可能对用于实现时变滤波器的方法非常灵敏。因此，时变滤波器可最有效地被实现为分窗口时间不变滤波器。换句话说，数据将被分隔成时间窗口，然后不同的时间不变滤波器将被应用于不同窗口。

上面的描述针对本发明特定实施例，这是为了说明本发明。然而，对本领域技术人员来说，显然可以对这里所述实施例做出许多修改和变化。例如，本发明方法不限于地震数据，并可应用于其中多尺度反演被用来避免局部最小值的任何数据。所有这类修改和变化都在权利要求限定的本发明的范围内。

参考文献

Bunks，C．，F．M．Salcck，S．Zaleski，G．Chavcn，1995，“Multiscalc seismic waveforminversion,”Geophysics，60，pp．1457-1473．

Maharramov，M．,U．Albcrtin，2007，“Localizcd image-difference wave-equationtomography,”SEG Annual Mccting Expanded Abstracts，San Antonio,2007，pp.3009-3013．Fallat，M.R.，Dosso，S.F.，“Gcoacoustic inversion via local，global，and hybrid algorithms,”Journalofthe AcousticalSocietyofAmerica105，pp．219-3230(1999)．

Pratt,R.G.,“Seismic waveform inversion in thc frcquency domain，Part 1：Theory andvcrification in a physical scale model,”Geophysics64，pp.888-901(1999).

Tarantola,A.,”Inversion of seismic reflection data in the acoustic approximation,″Geophysics 49，pp．1259-1266(1984)．

Claims (23)

1.一种反演测量的地震数据从而推导地下区域的物理特性模型的方法，其包括通过用目标函数的局部最小值优化在计算机上对测量的地震数据执行迭代、多尺度反演而连续更新所述模型，该目标函数计算模型模拟的地震数据和所述测量的地震数据之间的失配，其中改变低通滤波器——以下称为所述滤波器，被用于通过滤波失配计算中测量的和模拟的地震数据，将所述目标函数从一个尺度到另一个尺度连续改变，所述滤波是时变的，其中所述滤波器的低通截止频率随在所述反演的一个或更多尺度被滤波的地震数据的传播时间改变。

2.根据权利要求1所述的方法，其中所述滤波器的低通截止频率随增加的传播时间减小。

3.根据权利要求1所述的方法，其中每个连续改变的目标函数对应于被改变从而摒弃较少数据的滤波器。

4.根据权利要求3所述的方法，其中所述目标函数被连续改变直到达到最终尺度，在该最终尺度处所述目标函数未被改变，且其中迭代继续直到满足收敛标准或其他停止条件。

5.根据权利要求4所述的方法，其中对应于所述最终尺度的所述滤波器通过所有数据且不摒弃数据。

6.根据权利要求3所述的方法，其中所述滤波器被改变从而通过增加所述滤波器的低通截止频率来摒弃较少数据。

7.根据权利要求1所述的方法，其中在所述多尺度反演的每个尺度，即，针对高通滤波器的每个变化和所述目标函数的改变，执行一个或更多迭代。

8.根据权利要求1所述的方法，其中所述目标函数的局部最小值优化包括针对当前物理特性模型的一个或更多参数计算所述目标函数的梯度，然后搜索更新的物理特性模型，该更新的物理特性模型是所述当前物理特性模型在基于更好解释所述测量的地震数据的梯度的方向中的干扰。

9.根据权利要求8所述的方法，其中所述梯度是通过伴随方法计算的。

10.根据权利要求1所述的方法，其中所述反演是全波场反演。

11.根据权利要求2所述的方法，其中为了确保所述迭代收敛到所述目标函数的全局最小值，下面的标准被用于设计时变低通滤波器：当迭代开始时，所述模型应足够精确从而在所述时变滤波器被应用后，预测任何传播模式到在该模型的半周期内的传播时间。

12.根据权利要求11所述的方法，其中所述收敛标准可数学表达为：

$\frac{\underset{s,r}{\mathrm{stat}}\left|t_{\mathrm{error}}(M,s,r,t)\right|}{\underset{s,r}{\mathrm{statT}}(s,r,t)}\≤\frac{\mathrm{\α}}{2},$

其中M表示所述物理特性模型，t_error是测量的地震数据和模拟的地震数据之间的传播时间误差，且T是测量的数据的瞬时周期，α是所选的常数，t是地震波传播时间，s表示地震源坐标，r表示地震接收器坐标，且stat意味着降低s和r的可变性到常数的均数或平均值或另一个度量。

13.根据权利要求12所述的方法，其中α是在1/2到1之间的数字。

14.根据权利要求12所述的方法，其中所述低通截止频率随传播时间增加而减小的速率被确定是符合所述收敛标准的。

15.根据权利要求1所述的方法，进一步包括为不同源-接收器偏移使用不同时变滤波器。

16.根据权利要求9所述的方法，其中所述滤波器是用等式1中的所述反演协方差矩阵C^-1实现的。

17.根据权利要求16所述的方法，其中所述反演协方差矩阵C^-1被选择为在时间维中是非对角的，其中的元素等于所述时变滤波器的系数的时间表示，所述系数随时间变化。

18.根据权利要求12所述的方法，其中被假定是传播时间的线性函数，其形式为

其中M₀是开始迭代时所述模型的初始选择，β是所选常数，其中如果用这个滤波器设计迭代反演不收敛，则β的值增加，并会尝试具有的更保守估计的反演，并且继续直到发现为所述初始模型产生收敛的β。

19.根据权利要求1所述的方法，其中所述反演是全波场反演。

20.根据权利要求1所述的方法，其中计算机被编程以执行所述方法的至少部分步骤。

21.一种计算机程序产品，其包括其中具有计算机可读程序代码的计算机可用介质，所述计算机可读程序代码适于实施已测量数据的全波场反演，从而推导所述波场的传播介质模型的方法，所述方法包括：

通过用目标函数的局部最小值优化，执行已测量数据的迭代、多尺度反演来连续更新所述模型，该目标函数计算模型模拟的数据和已测量数据之间的失配；

其中改变低通滤波器，以下称为所述滤波器，被用于通过滤波失配计算中所述已测量数据和模拟的数据，从一个尺度到另一个尺度连续改变所述目标函数，所述滤波器是时变的；

其中所述滤波器的低通截止频率随通过所述介质的波场传播时间改变，该传播时间对应于在所述反演的一个或更多尺度被滤波的数据。

22.根据权利要求21所述的方法，其中所述滤波器的低通截止频率随传播时间增加而减小。

23.根据权利要求21所述的方法，其中被反演的所述数据是地震数据，且所述传播介质是地球的地下区域。

Images