_设计生物化学实验动物组织DNA的提取

、STRATA程序的原理1.0前言STRATA程序的基本原理从图1.1可容易看出。图1.1STRATA基本原理从包括迭加地震数据和速度/密度信息的两个输入流（以井的测井曲线或RMS速度格式），得到一个基本的地下速度模型，然后用这个模型做整个地震数据的反演。虽然这个目的简单，但达到这个目的要产生许多困难。因此，这个程序包括做下面事的算法：合成地震记录的产生交互的测井曲线的拉伸和压缩子波提取叠后地震处理地震拾取用垂直的和横向的插值建立模型应用几种反演算法反演这些项目的相互作用在图1.2已经表明。图1.2STRATA各部分的相互关系程序设计具有高质量的图解方式和用户友好的菜单控制。对于程序运行和特定菜单项目的信息，参考这个用户手册的前面部分。在这一章里我们将叙述用于STRATA的主要算法的原理。所研究的问题包括褶积模型、反褶积、反演和属性。假设程序的其它部分，如测井曲线的拉伸和压缩、处理步骤象滤波和AGC、模型建立用户十分熟悉，或当用这个程序时是显然的。速度／密度控制建立模型地震数据地震反演测井数据地震处理合成地震记录产生地震数据拉伸/压缩测井曲线子波提取建立模型地震拾取反演2.0褶积模型地震道基本的褶积模型可写为：公式中：r(j)=以时间序列表示的地表零偏移距反射系数；W(i)=地震子波，假设固定；N(i)=附加的测量噪音。注意在这个模型里，假设多种被忽略。反演可认为是给定地震道T(i),求去反射系数r(j)的处理。在公式2-1中，反射系数与地表波阻抗的关系是公式：公式中：I(i)=ρ(j)V(j)ρ=密度V=P波速度在第5章里我们将看到，反演的目的是从地震道估算地表的速度。显然，首先需要从褶积模型中提取一个估算的反射系数，然后研究在第4章论述的相关的反褶积处理。在考虑反褶积或反演前，让我们更详细的看看褶积模型的两个主要的组成部分，子波和反射系数。2.1反射系数当震源的能量通过爆炸或一个座板突然碰撞地面释放时，这种能量以弹性波通过地层传播。人们熟悉的最简单类型的波是压缩波（纵波），象声音一种声波干扰。岩石允许声波通过的性能是通过岩石纵波的速度和密度产生的声波阻抗给定的（类似于电路里电阻的原理）。同预期的一样，岩石越强，声波阻抗越高。例如，一般压实砂岩具有比泥岩高声波阻抗（孔隙度也影响速度：孔隙度越高，速度越低）。每当地震反射出现就有声波阻抗的改变。如果我们研究两个不同声波阻抗的层之间地层走向边界的声波幅度，反射幅度可记为上述的公式2-2。公式2-2表明反射系数可正可负，它取决于I(j-1)或I(j)比较大，但它的绝对值大小不超过1。对下行波从上地层走向边界的上述分析，上行波的反射系数的确相对于下行波的反射系数是负的。显然，不是所有的入射幅度都能反射（虽然海上记录在空气-水界面这种情况，几乎反射系数接近1）。传播幅度的大小入射幅度和反射幅度之间是不同的，或：注意如果r为负，T大于1，也就大于入射幅度！解决这种明显矛盾的方法只有通过物理系统记录能量被保存得知，而不是幅度。能量与下行和上行的路径有关。对于下行路径，我们仅看r符号的变化。因此，总的透射系数可记为：我们可看到，公式2-4告诉我们总的透射幅度总是小于1。注意，随层数的增加透射损失的影响也比较大。总的影响可表示为：jinjiWjriT)()1()()(2-1)1()()1()()(jIiIjIjIjR2-2)1()()1(2)(1)(jIjIjIjrjT2-322)(1))(1))((1(jrjrjrTway2-4公式中：N=层数；r(j)=第j个界面的反射系数；rAV=平均反射系数。也就是说，总的透射损失是透射通过我们感兴趣的每层界面上面的损失产生的。公式2-5也认为，如果我们用平均反射系数替代各个反射系数，透射影响近似等于两种图景透射系数N次幂。为了获得影响的设想，假设是100层叠加（采样率2ms,相当于200ms）。如果平均反射系数是0.1，透射比例系数是0.366。然而，如果平均反射系数是0.05，总的比例仅仅是0.779。重点注意，仅在P波是垂直入射界面时公式2-2到2-5成立。非垂直入射时，波型转换发生，产生反射和透射的横波。测定岩性地震学的精确量度时，为了推测这个影响是非常重要的。对于叠后的情况，我们假设公式2-2精确地预测反射振幅。2.2地震子波地震数据记录的第一步是震源图象的产生。陆地上，两种最普遍的震源是炸药—通常比较深的埋藏，和可控震源—长时间延续振动地面震源。海上记录，空气枪是最普遍的震源。如果我们知道震源的准确类型，我们可以用反褶积处理反演。参考确定性反褶积。然而，我们常常没有这种记录信息，必须用统计反褶积技术估计震源类型。（在STRATA中包括了这两种反褶积类型）。子波通过振幅和相位谱定义。相位谱的类型我们认为是零相位、常数相位、最小相位、和非最小相位。为了了解零相位和常数相位子波的原理，不同的振幅、频率和所有零相位和常数相位（即900）正弦波的数简单地加起。然而，最小相位不是作为简单的原理去理解，有许多不同的描述方法。最简单的是一种直观的方法，时间计算起点前最小相位子波的情况不存在（数学上我们说子波是由某种原因引起的），能量的大部分集中在波前。实地上，术语“最小相位”指的是所有具有相同振幅和不同相位谱的因果子波的事实，唯一的它的相位谱最接近零的是最小相位子波。原来子波相位大量的有效信息可通过应用简单两点子波获得W1=（2，1）2-6W2=（1，2）2-7子波W1是最小相位，而子波W2是非最小相位。事实上，把子波W2称为最大相位。让我们现在推导有关这两个子波振幅和相位谱的信息。首先，设想该公式2-6的子波Z变换：W（Z）=2+Z2-8现在，作代换公式中j=√-1ω=2πf注意这是合理的代换，由于2-9复指数包含余弦和正弦函数，是傅立叶变换基本原理．现在把公式2-9代人公式2-8，我们就得到子波的傅立叶变换:)sin()cos(tjteZtj2-9ImRe)sin()cos(2)(jtjteWtj2-10NNjAVtotalrjrT)1())(1(1222-5公式中Re＝傅立叶变换的实数部分＝2+cos(ωt)Im＝傅立叶变换的虚数部分＝sin(ωt)为了从实数部分和虚数部分得到傅立叶变换的振幅和相位谱，我们以下面的方法从直角坐标变换到极坐标：公式中|W(ω)|＝Ｗ的振幅谱θ(f)＝Ｗ的相位谱用类似的方法，Ｗ２Z的变换和傅立叶变换可以用下式表示:公式中;Ｗ1和Ｗ２振幅和相位曲线图如图2.1所示注意两种信号的振幅谱是相同的,而它们的相位谱不同，图2.1没有，子波Ｗ1具有最小相位谱。子波Ｗ2是非最小相位谱。事实上，子波Ｗ2具有最大相位的特性。首先，考虑三点子波的Z变换：傅立叶变换可记为：振幅和相位谱可用前面导出的关系计算。让我们通过三种可能的组合旋转子波Ｗ1和Ｗ2产生三点子波：W3=W1*W1=（4，4，1）2-16W4=W1*W2=（2，5，2）2-17)()()(fjeWW2-112/12/122/122)]cos(45[)](sin)(cos)cos(44[)Im(Re2tttt))]cos(2/()(arctan[sinRe]arctan[Im/tt)(2222)()(21)(jeWWZZW2-122-13))]cos(21/()sin(2arctan[)()]cos(45[)(22/12tttW2210)(ZaZaaZA2-14))2sin()sin(())2cos()cos(()(212102210tatajtataaeaeaaAtjtj2-15W5=W2*W2=（1，4，4）2-18这三个子波的振幅谱如图2.2顶部所示，是相同的。当我们认识到子波Ｗ1和Ｗ2的振幅谱是相同的这有意义。相位谱如图2.2底部所示。注意子波Ｗ3应该是最小相位，因为它是两个最小相位子波的褶积。同样，子波Ｗ5应该是最大相位，因为它是两个最大相位子波的褶积。但子波Ｗ4是最小和最大相位子波褶积它怎么样呢？这称之为混和相位子波。从图2.2可以看出，混和相位子波的相位谱在最小和最大相位子波的相位谱的之间。叙述的原理仅概括出了N脉冲子波。如子波可通过N-12-点子波褶积，这些子波的每一个都具有形式：Wi=(ai,bi)2-19如果ab任何2-脉冲子波是最小相位，如果ab是最大相位（如果a=b，我们得到位置称为单位园上的原点，可能引起问题，碰巧，很少出现）。下面三种叙述可能产生：1）如果所有2-脉冲子波都是最小相位，所得到的N-脉冲子波是最小相位；2）如果所有2-脉冲子波都是最大相位，所得到的N-脉冲子波是最大相位；3）如果2-脉冲子波是最小和最大相位组合，所得到的N-脉冲子波是混合相位；。反褶积可认为是一个两步处理。第一；我们寻找子波。第二，我们通过应用其反演除去子波。再提到反褶积部分的注解之前，让我们提一个问题：任何子波都有准确的反演吗？我们将通过研究公式2-6和2-6的最小相位2-点子波试着回答的这个问题。通过公式2-8的多项式除法子波的反演可记为：Wt-1=f=(1/2,-1/4,1/8,-1/16,…)2-20也就是说，子波具有无限多的反演，为了应用子波反演，因此，必须在一些点上截取。怎样大小的井可做反演运算，我们可比较最后输出的理想结果的脉冲。让我们首先考虑截取长度为2的反演滤波器。Wtft=(2,1)*(1/2,-1/4)=(1,0,-1/4)2-21应用一个长为3的反演滤波器长，我们得到Wtft=(2,1)*(1/2,-1/4，1/8)（1，0，0，-1/8）2-22注意最后的样点总有一个错误项。虽然算子的长度越长，输出错误越小，但对于无限长的算子仅趋于零。当我们寻找时，精确的反演总不是最好的。再第三章里，我们将研究解反演的最小二乘方法。2.0子波提取为了完成反演STRATA需要有关地震子波的信息。在地震处理中子波的问题是一个复杂的问题，是现行研究的普遍领域。虽然已产生了许多子波提取方法，但下面做出了概括的叙述：1．在频率域里，我们可认为子波提取的问题由两部分组成：确定振幅谱；确定相位谱。这两部分，确定相位谱要困难的多，而且是存在于反演中的主要的错误源。2．子波提取方法分成三个主要的类型：a)完全确定性：这意味着直接应用地表接受器和其它方法测量子波：b)纯统计：这意味着仅从地震数据中确定子波。这些方法确定可靠的相位谱往往有困难。c)应用测井：这意味着除了地震数据还要用测井信息。原理上，在井位置这可提供准确的相位信息问题是这种方法决定性的取决于测井和地震之间的好的联系。特别是，把深度采样的测井转换到双程旅行时的深-时转换可能引起模糊不清降低结果。3．作为旅行时的函数子波从一道到另一道可能改变。这意味着每个地震剖面子波处理过程见确定一大组子波。实际上，设法确定可变子波比数据可以分解可引起更多的误差。实际和有用的解决办法是整个剖面提取一个单独的平均子波。下面部分叙述在STRATA中子波提取能力。3.1统计子波提取在STRATA中统计子波提取方法仅应用地震道提取子波。这个方法不计算相位谱单必须由用户提供一个独立的参数。相位谱选择上：a)等于某个数的常相位（即450）。b)最小相位。如下所示应用地震道的次自相关计算振幅谱，对每一道分析：1）提取分析时窗；2）时窗的开始和结束斜坡，斜坡长度适当小（10样点，时窗的1/4）；3）计算数据时窗的自相关。自相关长度等于期望子波长度的二分之一；4）计算自相关的振幅谱；5）获得自相关谱的平方根，这接近子波的振幅谱；6）加期望的相位；7）做反快速傅立叶变换产生子波；8）在分析时窗里加从其它道计算的子波。注意在这个过程中，在确定道振幅谱上完成的有效滤波中子波长度是关键性的参数。当增加子波长度，子波谱接近数据时窗。3.2应用全测井子波提取STRATA应用的子波提取方法的第二个类型需要利用测井。STRATA以两种方法应用。一