去噪学习笔记

ZoneAnomalyProcess（区域异常处理）分频、多域去噪手段在地震记录中，有时还会出现一些孤立的大脉冲振幅和振幅或频率异常的能量团，它们既不具备相干噪声的性质也不符合随机噪声的条件。Omega系统中有一个称为区域异常处理的方法，是专为压制这样的异常振幅或频率而设计的，即ZoneAnomalyProcess（区域异常处理）。区域异常处理是一种地表一致性处理，它的基本原理是：假定在炮点i激发，检波点j接收到来自深度h的反射振幅值为A..h，影响A..h的因素有激发条件、接收条件、偏移距大小和地下介面的反射性质。前三项ijh ijh称为地表一致性因素，一般它们应当是比较一致的。如果在某个地方出现了大的变化，说明这里有异常值。第四项是地质因素，这是我们的有用信息，称为地下一致性因素。区域异常处理就是要计算一个新的振幅量度Bijh并得到比例因子Cijh=Bijh/A..h再将这个比例因子应用于地震道数据。ijh， ijhijhijh,这样做横向上不会改变均方根平均振幅，但时窗与时窗之间的振幅得到了调整。从而达到衰减振幅异常和频率异常的目的。例如，假定一个时窗内的均方根平均振幅是完全地表地下一致性的，那么C..h=l•如果时ijh窗内存在一个脉冲噪声，且它又不是地表一致或地下一致的，则C..h<1,用该比例因子乘时窗内数据，使ijh其能量受到压制。在一个时窗内是否存在异常用以下四个振幅统计量来判定：均方根振幅、平均绝对值振幅、最大绝对值振幅和优势频率。这几个统计量可以用其中一个或多个。具体用哪些统计量根据噪声的类型而定。均方根振幅、平均绝对值振幅适用于高振福噪声团，而最大绝对值振幅和优势频率则对衰减分散的高值脉冲效果较四种振幅统计量的定义如下：A设：i=时窗序号t=时窗起始时间N=时窗长度j=时窗内样点序号（从t到t+N）a(j)=样点j的振幅值则四川石油管理局地球物理勘探公司

对这些振幅统计量还要用高斯-賽德尔代法分解成共炮点、共检波点、共偏移距和共中心点分量，再将各道的分量加在一起，最终得到各道在指定时窗内的振幅量度Bijh,它与实际反射振幅值Aijh的比值即是异常振幅压制因子。四川石油管理局地球物理勘探公司区域异常处理算法设在炮点iit发,检波点j接收到来自深度h的反射振幅值为Aijh.Aijh的大小由激发条件的彩响、接收条件的形响、偏移距的大小和地下反射介面性质（共中心点）确定.因此有：Aijh-Si*RjxGkhxMn其中,Si与炮点i有关的振幅分童Si与炮点i有关的振幅分童Rj与检波点j有关的振幅分量Gkh-与第k共中心点的反射层位h有关的振幅分量,k=（i+j）/2与保移距n有关的振搐分童，n-i-j四川石油管理局地球物理勖探公司对（1）式两端取对數，得:logAljh-logSi+logRj+logGkh+logMn如果假定偏移距项也与反射屡位h有关，则（2）式可写为:logAljh-logSl+logRJ+logGkh+logMnhMnh表示与偏移距11有关的反射屋位h的掘幅分童.达到最小。在以上（3）式戎（2）式中，变化1.j,将得到一系列方程.方程个数一般大于未知數个數.单其中独立的方程个數一般小于未知數个數.因此该方程组无确定解.为了求解该方程组，可用最小平方法，即要求得到的解使误差能量达到最小。E=给-k>gSTgRj-k>gG祓-MJ?四川石油管理局地球物理勘探公司区域异常处理的实现步骤第一步：振幅统计量计算使用模块：ZAP.P