地震子波及其提取的方法的研究.docVIP

地震子波及其提取的方法的研究 　　摘要： 地震子波的提取和研究是反演和正演模型的关键工作，论文重点研究了复赛谱域子波提取方法的原理和所存在的问题，并用实际资料进行试算，从而得出一些有益结论。并将提取子波运用于合成地震记录和反褶积。 　　Abstract： How to extract seismic wavelet accurately is a key problem of seismic inversion. The paper focuses on the principles and problems in spectrum domain wavelet extraction method. In order to draw some useful conclusions， it uses the actual data to a spreadsheet. And extract the wavelet used in the synthetic seismograms and deconvolution. 　　关键词： 子波；提取方法；复赛谱域 　　Key words： wavelet；extraction method；spectrum 　　中图分类号：P631.4+43 文献标识码：A 文章编号：1006-4311（2015）01-0315-02 　　1 子波基本理论 　　子波包括最小相位子波、最大相位子波、混合相位子波。子波的Z变换是一个多项式。若此多项式的全部零点均在单位圆外，则为最小相位子波；在单位圆内为最大相位子波，零点在单位圆的内外都有，则为混合相位子波。 　　2 复赛谱域提取子波法 　　复赛谱域提取子波法研究的是一种非线性系统，或者说是一种线性滤波的推广。对于一个线性系统，它可以满足以下迭加原理：T[ax1+bx2]=aT[x1]+bT[x2]（1） 　　式中a、b为常数，T为表示函数。一个信号是由多个信号相加合成的，利用线性滤波的方法，就能有效地把它们分开。如果一个信号不是由多个信号相加，而是由多个信号中褶积合成的，那么，用这种线性滤波的方法就不能将其分开。对于这种信号的处理，我们仍可利用（2）式的基本思想，把它推广到褶积合成的信号。为此我们需要一个称为同态系统的变换系统D，它具有把一个褶积关系转变为相加关系的特征系统。 　　从（2）式可见，线性系统特别适用于分开相加合成的信号；也就是说，一个系统可以这样处理： 　　D[（a）x1*（b）x2]=aD[x1]+bD[x2]（2） 　　上式中（a）表示标量乘，（a）x1表示x1与自身褶积a次，（b）也是如此。 　　系统L是一个一般的线性系统，满足： 　　L[ax1+bx2]=aL[x1]+bL[x2] L[x?？啄（t±？子）]=L[x]？啄（t±？子）（3） 　　系统D-1为D的逆变换特征系统，是把经过处理的相加关系转化为褶积关系输出，即： 　　D-1[ax1+bx2]=（a）D-1[x1]*（b）D[x2]（4） 　　系统定义为这种D→L→D-1的标准格式，其突出的优点是，只要确定了特征系统D之后，它在整个反褶积过程中保持不变，并将一个非线性过程转化为线性运算。 　　我们以地震记录的数学模型为例，该模型的物理意义表示：野外地震信息是由地震子波bt和地层响应序列（反射系数序列）？孜t按褶积模型关系组成的。当地层之间的厚度与子波的视波长相当，层间波的旅行时间比子波的延续时间短时，则层与层之间的反射将发生干涉，这使得地震子波和反射序列不能在记录上清楚地表现出来，要想直接把两者分开是困难的。但是我们可以这样来处理，把bt和？孜t在时间域的褶积关系，经过Z变换后变成相乘的关系，然后再取对数（频谱取对数之后称为对数谱，与之对应的时间函数称为对数谱序列、或者音译为复赛谱），把相乘的关系变为相加系统。这样就可以把褶积分量以复赛谱的形式分开，如果两序列的复赛谱有差异的话，则可用一般线性滤波的方法，很容易实现反褶积。其具体过程是通过以下定义的几个关系式和步骤来完成： 　　此处arg表示幅角。其中B（Z）和E（Z）分别为子波bt和地层响应序列？孜t的Z变换，X（Z）为地震记录xt的Z变换。而xn、、bn为相应的离散时间序列。 　　积分回路c为下式所确定的圆周： 　　n称为对数谱序列，但必须注意的是，它并不是谱，没有谱的含意。用以上关系式的变换，就实现了将褶积运算变换为相加运算。 　　3 实验结果及其分析 　　合成理论地震记录由反射系数与理论子波褶积形成理论记录：x（t）=R（t）*b（t）（8） 　　其中R（t）反射系数，b（t）地震子波，x（t）地震记录。分别用理论最小相位子波和混合相位子波合成地震记录，进行子波提取