3． 作为旅行时的函数子波从一道到另一道可能改变。这意味着每个地震剖面子波处理
过程见确定一大组子波。实际上，设法确定可变子波比数据可以分解可引起更多的误差。实际和有用的解决办法是整个剖面提取一个单独的平均子波。 下面部分叙述在STRATA中子波提取能力。
3.1 统计子波提取
在STRATA中统计子波提取方法仅应用地震道提取子波。这个方法不计算相位谱单必须由用户提供一个独立的参数。相位谱选择上：
a) 等于某个数的常相位（即450）。 b) 最小相位。
如下所示应用地震道的次自相关计算振幅谱，对每一道分析： 1） 提取分析时窗；
2） 时窗的开始和结束斜坡，斜坡长度适当小（10样点，时窗的1/4）； 3） 计算数据时窗的自相关。自相关长度等于期望子波长度的二分之一； 4） 计算自相关的振幅谱；
5） 获得自相关谱的平方根，这接近子波的振幅谱； 6） 加期望的相位；
7） 做反快速傅立叶变换产生子波；
8） 在分析时窗里加从其它道计算的子波。 注意在这个过程中，在确定道振幅谱上完成的有效滤波中子波长度是关键性的参数。当增加子波长度，子波谱接近数据时窗。
3.2 应用全测井子波提取
STRATA应用的子波提取方法的第二个类型需要利用测井。STRATA以两种方法应用。一种方法是应用测井确定子波的全振幅和相位谱。第二种方法利用测井应用与上面叙述的统计方法结合仅确定常相位。
首先在子波提取采单上选择选件全子波应用测井引用测井方法。这个方法需要小对每道有用的密度和声波测井分析。当然，由于测井仅是可获得的独立点，STRATA提供了漏侧测井以相同的方法通过外推和插值建立反演模型。着意味着在子波提取中把测井校正的影响（拉伸）和数据拾取引入测井应用。着允许应用井位置周围的道。
对每道分析： 1） 提取声波、密度、地震数据分析时窗； 2） 声波和密度相乘得到波阻抗。从波阻抗中计算反射系数； 3） 时窗的开始和结束反射系数序列和地震数据的斜坡，斜坡长度适当小（10样点，时
窗的1/4）； 4） 计算最小平方整形滤波器W，解决下列问题：
TRACE=W*REFLECTIVITY
由于子波的样点数一般小于道收据的样点数，也就是相当于应用最小平方体解超定的线性系统 5） 应用希尔伯特变换计算子波的振幅包络。如果包络的波峰偏离了时间零，改变测井
和地震道之间的互相关，应用第四步重新计算子波。相加之前确保测井和地震道之间随机时间偏移一道一道校正。 6） 该子波与其它子波相加计算另外一道。 7） 如下所述通过滤除高频成分稳定计算子波：对于振幅小于最大振幅1/4的原始数据
时窗谱里的每个频率值。零对应提取子波的成分。 全子波提取方法具有计算精确子波的优点，但也具有测井和地震数据之间联系非常敏感的缺点。实际上，时间调整或拉伸错可能引起迅速降低，表征为子波高频丢失、相位谱畸变、不实际的旁瓣产生。
3.3 应用测井常相位子波提取
STRATA中子波提取的第三种类型是一二两种类型的混合，该方法通过选择选件：应用常相位测井引用。在该方法中，仅应用地震数据计算振幅谱，在3.1节中已正确描述。相位谱假定为一个通过解一个自由度（平均相位）的最小平方整型滤波器而确定的常数。反射系数应用测井确定单个常相位数。因此，该方法在井联系不完善的情况下达到最稳固。
4.0反褶积
4.1维纳—莱文森整型滤波器
反褶积的目的是寻找把输入子波转换成某种期望的输出型状的滤波器。随然该期望输出常常是一个单脉冲，我们可设定为我们期望的任何。该问题典型的解法是假定在期望输出和实际输出之间最小二乘拟合。这称为维纳—莱文森整型滤波器。为了简化事情，我们在两点子波上完成我们的分析然后假设一般情况成立。
首先，让
Wt=（W0，W1） （4-1） 为子波，
ft=（f0，f1） （4-2） 为子波滤波器，和
dt=（d0，d1，d2，?） （4-3） 为期望结果。 然后，我们看
xt= ft* Wt= （f0 W0，f1 W0+ f0 W1，f1 W1） （4-4） 是实际结果。
现在，最小二乘方方法告诉我们期望和实际的结果之间差的平方和必须为最小值。为了寻找最小值，我们对滤波器值和设置结果平方差求导为零。也就是：
dI?0dfi 公式中：
222I?e0?e1?E2（4-5） et?(d0?x0,d1?x1,d2?x2) 把公式4-4和4-3代入4-5得到：
I?(d0?f0w0)2?(d1?f1w0?f0w1)2?(d2?f1w1)2 意味着：
（4-6） dI??2w0(d0?f0w0)?2w1(d1?f1w0?f0w1)?0df0dI??2w0(d1?f1w0?f0w1)?2w1(1?f1w1)?0df1
（4-7） （4-8） 注意公式4-7和4-8代表两个有两个未知数，f0和f1的线性方程，以矩阵的形式可表示为：
?w02?w12w1w0??f0??d0w0?d1w1??????dw22??fwww?w?101???0101??现在，子波Wt的自相关可写为：
（4-9） ?ww?(?ww(?1),?ww(0),?ww(?1))?(w0w1,w?w,w1w0)同样，期望的结果和子波的互相关是：
2021（4-10） ?dw?(?dw(?1)，?dw（0），?dw(?))?(d0w1,d0w0?d1w1,d1w0)公式4-10和4-11表明公式4-9可重新表示为：
（4-11） ??ww(0)?ww(?1)??f0???dw(0)??????????ww(?1)?ww(0)??f1???dw(?1)?归纳N-点子波的情况给出维纳莱文森方程： Rf=g (4-13) 公式中：R=输入的自相关矩阵； F=期望的滤波器；
g =期望输出与输入的互相关。 解该方程：
f=R-1g (4-14) 公式中： R-1 =R的逆矩阵。
注意：为了解2×2个问题，应该注意矩阵：
（4-12） ?ab?A???cd?? 通过下式给出矩阵的逆：
A??b??d1?ad?bc???ca?4.2 最小相位统计反褶积（脉冲反褶积）
给出的公式是很一般的，且有许多应用。这种应用的一种是脉冲反褶积。对于最小相位子波让我们看看脉冲反褶积。
Wmin??(2,1)（4-15） 在这种情况下，期望结果是一个在原点的尖脉冲，或 Dt=(1,0,0) (4-16) 然后我们看出：
?ww?dw?(1,2,0)?52??f0?零延迟?2 ???25??f??0??(2,5??1????,2)（4-19） （4-17） （4-18） 我们获得维纳莱—文森方程的下列格式：
解为：
应用反褶积滤波器，我们看出：
但如果我们选择了最大相位子波反褶那么：
?f0??5?2??2??10/21????1/21??25??0????4/21?（4-20） f1min4/?21)?2/21,?4/21)?,1)?(10/21?,???(20/21,??W??ft?(2 Wmax=（1，2） （4-22）
在这种情况下，自相关是：
（4-21） ?ww?(2,5,2)（4-23） 与最小相位子波自相关相同。由下式给出互相关：
gt = （1，0） （4-24）
在最小相位的情况下，是一种互相关的简单换算形式。因此最终的维纳方程可写为：
?5?2?5?5???f0??1???f??0???1??（4-25） 解为：
?f0??5???1/21??2??f1??2?5???1??5/21???0???2/21?（4-26） ???? 应用滤波器，我们得到：
Wmax?ft?(1,2)?(5/21,?2/21)?(5/21,8/21,?4/21)（4-27） 根本不是尖脉冲。
总之，如果应用了STRATA的最小相位统计反褶积，你必须认识到最好在最小位上工作。
在STRATA中应用零相位统计反褶积方法与按推导出的反子波（算子）的振幅谱的方法的最小相位选件相同。因此，在应用之前反褶积算子设置为零。
4.3 反褶积应用情况 在上述两节里，我们已看到最小相位和零相位统计反褶积的理论方面。我们应用非常简单的2—点子波的例子。因此当应用真实数据时 工作原理怎样？
对于真实的数据，我们假设公式2-1的地震道模型成立。也就是说，我们只知道地震道而不知道反射系数或单个的子波。因此，可给出统计反褶积的主要假设：
地震道的自相关相当于子波的自相关。
你可看到该假设在实际中不完全正确。同样地，由于获得自相关函数， 我们丢失了有关子波相位的所有信息。因此，我们必须做最小相位或零相位的假设。
如果我们接受子波的自相关和地震道的自相关是相同的前提，为了优化反褶积我们仍有可用的三个参数。
1． 算子长度（相当于自相关本身的长度）； 2． 控制反褶积灵敏度的预白噪系数； 3． 设置用于自相关的时窗的长度和位置。
在STRATA程序中算子长度的缺省值是60ms，应用标准的地震数据证明是合理的值。如果输入数据分辨率很低，而你希望结果的有效高频加强，你也许要附带有一个比较高的数。如果数据已经是高分辨率的，你可以去掉一比较低的数。用不同的算子长度做几个测试是可行的。
本质上，预白噪是计算逆算子之前加到数据振幅谱（时间域预白噪实际上是加到零延迟自相关系数）上噪音级。显然，如果不加预白噪可获得最好的结果。因此，某些预白噪必须做稳定化运算。在程序里，我们建议的1%的预白噪。
最后，设计时窗在反褶积结果的质量上是一个重要的系数。人们凭经验方法说指定时窗至少10的时间算子长度。把时窗定在重要的区域上是关键的，尤其除去数据的噪音部分。
必须强调反褶积充其量不过是一种技术，而不是一种科学。开始整个反褶积之前在数据体的一小块上测试不同参数的组合总是可行的。
4.4 双边的脉冲反褶积
注意脉冲反褶积的结果是时间计算起点的脉冲和在延迟时处的误差.我们可通过做比较长的算子延伸误差。为什么期望的结果在某个大于零的时刻没有脉冲呢？这具有使误差项更对称的作用。例如,让：
dt=（0，1，0） (4-28)
因此维纳-莱文森方程成为：
?5?2?2?5???f0??1??f???2????1?（4-29） 解为：
?f0??5?f??1/21??2??1??2?5???1??1/21??2???8/21?????（4-30） 应用滤波器得出反褶积的解：
?t?wt?ft?(2/21,17/21,8/21)w（4-31） 这称为双边脉冲反褶积。虽然结果不完全对称，但比第一个解更对称。因此，以这种方
法，输入子波必须已知，单边脉冲反褶积不是这种情况。
在STRATA程序中方法描述仅用于全反的确定性反褶积。如果选择仅相移选件，仅应用子波的反相位谱。
5.0 反演
叠后地震反演是分析迭加地震道和设法重建地下的速度和波阻抗结构。基于1-D褶积模型的反演主要描写为：
T(i)??r(j)W(i?j?1)n(i)j（5-1） 公式中：r(j)=以时间表示的地下零偏移距反射系数； W(i)=地震子波，假设为常数； N(i) =附加的测量噪音。
注意在该模型里假设多次波忽略不计。
反演可认为是给定地震T(i)确定反射系数r(j)的过程。由于反射系数与一系列地层的波阻抗有关：
I(j)?I(j?1)r(j)?I(j)??(Ij()j?1)I((jj))v?（5-2） 公式中：
如第一节所述，反演同样可确定地下的波阻抗。 在STRATA中应用了两种反演方法，它们主要的不同在于地震子波W的处理和与反演相关的不唯一性问题的处理。在下节里详细的描述它们。