matlab 非线性插值 基于非线性相关地震波场插值



     在地震数据处理中,除因数据缺失而直接丢失一部分信息外,还可能导致在后续处理流程中产生噪声,使得地震道中出现形态各异的脉冲。这样,基于多道处理算法的处理过程将受到缺失道的影响,其中受影响最严重的包括波动方程偏移,基于波动方程的抑制表面多次波的方法,与地表有关的多次波的消除和谱估计等处理过程。

  道内插是地震资料常规处理中不可缺少的步骤,其目的是在偏移之前加密空间方向(特别是三维横线方向)的采样,防止因采样过疏、在偏移过程中产生空间假频,从而提高小构造的成像精度。

  1、非线性相关技术

  把非线性相关技术用于地震道插值以填充地震勘探中的空白段,取代噪声道,产生均匀间隔的排列。首先进行初始对正(对叠前资料进行NMO校正,对叠后资料则人工插入相关线),然后把相邻地震道的相应特征一一相应对正,并确定出地震道之间的旅行时差。用这些资料构成的合成道(内插道)无论相关道之间的距离有多大,都可保持单个反射面的倾角和振幅的连续变化。

  按照Martinson等人(1982)的设想,假设两个待相关的地震道和(=接收点,t=时间),它们的波形相似,但其中的一个地震道相对于另一个地震道经历了时间轴的拉伸或压缩。这种拉伸或压缩的“畸变”可由未知形状的映射函数来描述,它使一个地震道的时间的特征与另一个地震道的时间t=T的对应特征相联系。

  Martinson等人(1982)把映射函数参数化为一条线性趋势线,此线是用含有未知系数的截断傅里叶正弦级数改形而得到的。此法试图寻找产生映射函数的系数值,该映射函数的型式能使描述两个地震道之间相似程度的相关系数(或相干系数)C最大。这样,相关系数C值视映射函数系数值而定。映射函数使地震道之间相关最好时C达到它的最大值,梯度向量中的所有为零。

 matlab 非线性插值 基于非线性相关地震波场插值
  用线性反演法以迭代方式求取系数使相关性最好。给系数使一个小的增量,则C也获得一个增量,所以。这对值就有了一个约束。然后要求必须处在使增加最快的方向上,即必须平行于,或,这样其他约束条件也可实现了,这些约束便于求出。

  2 非线性相关技术道间插值模式

  特殊情况下,如果是两个相邻相关道的映射函数,对于地震波场u(x,t)中的接收点和,有

  这里在最小平方意义上使得误差函数极小。

  由任何两个采样道之间所确定的u(x,t)的内插由下式得出

  式中:;

  为采样道间距。

  时间由稳定方程得出

  式中表示给定地震同相轴的采样道之间的相移。

  其中,相关系数C是用来定量描述两个波形的相似程度的。假设有两个波和,取某个适当的数a,使得波和()相接近,用误差能量来定量衡量相接近的程度

  和的能量和往往是确定的,所以的大小就由

  确定。被称为和的相关系数,这样便可得相关系数C。

  在讨论两个波形相似性时,常常还要分析两个波形相对时移后的相似性,即计算与的相关系数(相对于时移为)

  当变化时,就是的函数,是相对于的延迟时间,即时移。

  3、实际资料结果

  模型抽出一道地震道(该道位于断层上盘)的插值前、待插值以及插值后。对比分析可以看出经过非线性相关技术对实际断层构造模型进行道间插值,输出的结果同插值前相比较,地震道波形自然,地震剖面基本一致,能与原始剖面较好的吻合。

  

爱华网本文地址 » http://www.413yy.cn/a/9101032201/33794.html

更多阅读

插值法内插法、试误法 简要说明 技术权益的简要说明

学财务管理的时候,很多地方都会用到插值法,基本每本财务管理的书都会有,比较常用的有求适当的折现率使某项目的现金净流量为零、求债券的到期收益率、内含报酬率的计算等。其基本原理为在抛物线f(x)在很短的一段里为近似直线y=a+bx,再

ArcGIS插值图分类面积统计_阳光温热 arcgis反距离权重插值

在用ArcGIS进行地统计分析时,遇到一种情况,在插值完成后,想根据插值等级进行面积的统计分析,通过查资料解决了这个问题,现在分享一下。1.根据已有数据进行克里金插值。具体的探索性数据分析、数据变换、交叉验证之类的步骤,可以查找网上其

matlab中的最大值和最小值_待_雪 函数最大值和最小值

matlab中的的最大值和最小值MATLAB提供的求数据序列的最大值和最小值的函数分别为max和min,两个函数的调用格式和操作过程类似。1.求向量的最大值和最小值求一个向量X的最大值的函数有两种调用格式,分别是:(1) y=max(X):返回向量X的最

声明:《matlab 非线性插值 基于非线性相关地震波场插值》为网友温柔在手心分享!如侵犯到您的合法权益请联系我们删除