振幅补偿课件

欢迎各位光临！地震资料的振幅补偿一、概述二、几何扩散补偿三、吸收补偿四、地表一致性振幅补偿五、剩余振幅补偿六、其它振幅处理主要内容引起地震反射波振幅改变的原因十分复杂，但大致可分为以下几个方面：1、几何扩散一个点激发产生的振动能量以球面波的形式向外传播，因为球面面积与半径（传播距离）的平方成正比，所以波的能量与传播距离的平方成反比：

2、吸收衰减由于地下介质的非完全弹性性质，地震波在传播过程中能量被介质吸收转化为热能而引起衰减，这种衰减随传播距离的增加呈指数下降：其中，为衰减系数。

3、透射损失地震波遇到反射界面时，将发生反射和透射，反射和透射能量的总和等于入射波的能量。

下层界面的入射波是上层界面的透射波，它比上层界面入射波能量小，反射波能量自然也就小。更不说它在上传的过程中还要遭到更多的反射、透射。4、激发、接收条件的差异引起能量的不同

激发条件不同使炮记录之间能量有差别；接收条件不同使道记录之间能量有差别。

振幅补偿原则上是要对各种非地质因素对地震波能量的改变进行补偿

球面扩散引起深、浅层和远、近道的能量差异，球面扩散、吸收衰减和透射损失引起深、浅层的差异。引起远、近道的能量差异还有一个原因是入射角：入射角越大反射能量越小，这就是所谓的AVO效应。由于反射能量随入射角的变化包含有有用的地质信息，所以我们希望保留而不做校正。Omega系统的几何扩散补偿模块（GEOSPREAD_COMP）以上几何扩散补偿公式只有在介质均匀适用,水平层状介质中的地震波在不考虑吸收衰减的情况下,其传播衰减可近似地表示为,其中，V是沿传播路径的均方根速度，T为沿传播路径的旅行时。显然，在地震记录上同一炮点相同时刻接收到的反射波传播路径与炮检距X有关，因而均方根速度也是X的函数。为了补偿能量的扩散损失，需对地震数据乘以增益函数G(T)：式中，为所给定速度场中最小速度。这就是Omega系统的几何扩散补偿模块所用的基本几何扩散补偿增益函数。

一炮记录分为如下图所示的四个区域区域1区域2区域3区域4当我们使用几何扩散补偿模式和均方根速度平方内插方式时，各区域的补偿公式为：区域1－该区域使用常数因子，于是输出道的值为：

区域2－该区域之的速度由和线性内插，这时输出道的值为：

区域3－此区域内的速度由和线性内插，这时输出道的值为：区域4－该区域使用常数因子，输出道的值为：以上补偿公式看起来很复杂，特别是区域2和3，实际只是速度插值繁琐而已，其本质仍是使用补偿（或称增益）函数：该几何扩散补偿模块除了几何扩散补偿以外，还提供有吸收衰减补偿功能，其作用相当于Q补偿。要实现这个功能需要用户提供阿尔法、贝它、和嘎码等与地质有关的参数，所以在实际资料处理中并没有使用该功能。几何扩散补偿模块的主要参数1、模式选择参数组（*MODE_SELECT）（1）模式（MODE）该参数指定是做几何扩散补偿还是去几何扩散补偿。NORMAL：几何扩散补偿INVERSE：去几何扩散补偿缺省：几何扩散补偿（2）补偿类型（COREECTION_TYPE）该参数指定几何扩散补偿的类型。GEOSPRED：几何扩散补偿NEWMAN：纽曼补偿SECOND_ORDER：二次项补偿FOURTH_ORDER：四次项补偿缺省：几何扩散补偿（COMPUTED）（3）内插方式（INTERPOL_TYPE）该参数指定速度内插方式。RMS_VEL_SQUARED：均方根速度平方内插GAIN_FUNCTION：增益函数内插缺省（COMPUTED）：如果模式参数是去几何扩散补偿，缺省值是原几何扩散补偿的内插方式：如果模式参数是几何扩散补偿，则缺省值为均方根速度平方内插。

时间函数增益模块（TIME_FUNCTION_GAIN）

几何扩散效应在实际地震记录上的反映主要是地震波能量随记录时间的增大而减小，于是人们试图用随时间变化的增益函数来补偿几何扩散效应。Omega系统的时间函数增益模块就是基于这样的思想。它的做法是对地震道的各样点乘以其时间值的某次幂：其中的指数作为参数由用户指定。我们通常把这个模块叫做T补偿模块，对几何扩散补偿而言，它显然是很粗糙的。第一，增益函数与炮-检距没有关系，事实上从近道到远道也存在几何扩散的问题；第二，增益函数与速度无关，几何扩散的本质地震波随传播距离的衰减，因而必然与速度有关，即使从相对变化的意义上讲，也只有当速度为常数时，才可以仅从时间的变化来考虑。主要参数

时间函数增益模块增益从道头字START_TIME中的时间值开始起算，若该值不存在，则从第一个采样点起算。它的参数只有一个。增益指数（GAIN_EXPONENT）该参数指定时间增益指数。取值在－30到30之间，但不能等于0。吸收补偿也称指数增益补偿，其中β叫做增益系数。增益系数即衰减系数，它与介质的吸收特性有关。正确确定增益系数对做好指数增益补偿十分重要，但由于我们对介质的吸收特性并不掌握，何况它还应该是时变的，所以很难从理论上给出一个唯一确定的β值。在实际资料处理中，通常用以下两种办法来确定：1、参数试验，即用不同β值进行试验处理，那个处理结果令人满意就选它；2、将记录分时窗求出其平均能量，拟合出能量衰减的对数曲线，用对数曲线的衰减因子作为增益系数。

Omega系统的指数增益模块（EXPGAIN）该模块的功能是对地震道按时间进行指数增益补偿，其目的是补偿由于地层的吸收效应所引起的振幅衰减。补偿可以从该道的第一个采样点（TIME_FIRST_SAMP）开始，也可从该道的起始时间（START_TIME）开始。如果从第一个采样点开始，则补偿按以下公式计算：（1）（2）（3）如果从该道的起始时间开始，则补偿按以下公式计算：（4）（5）（6）正弦函数平滑增益函数如下图所示图中，S1是时间t1处的增益值；S2是时间t2处的增益值，它是按TL（或TL+TS）计算的；t2-t2即是转换区（t1和t2到TL的距离相等）;t0是正弦函数的原点时间。图中的实线是原计算的增益因子，虚线是正弦曲线。最后的增益因子使用是：t1之前和t2之后用实线的值，中间使用虚线的值。介质的吸收作用具有频率选择性，通常对低频吸收较少，高频吸收较多。指数增益仅对总体能量衰减进行补偿，所以补偿之后，仍保留有一定的地质吸收特征。指数增益模块的主要参数（1）增益指数（GAIN_EXPONENT）该参数即公式中的PAR，它的取值范围在－5到5之间。无缺省。（2）以分贝方式指定的增益指数（DB_GAIN_EXPONENT）该参数以分贝的方式中确定PAR的值：它的取值范围在－40到40之间。无缺省。（3）增益起始点时间（GAIN_BEGIN_TIME）该参数指定增益的起始点时间，有两种方式可供选择：TIME_FIRST_SAMP：从第一个采样点开始；START_TIME：从该道的起始时间开始。（4）增益长度（GAIN_LENGTH）该参数指定增益的时间长度，当时间大于本参数与增益起始点时间之和时，增益函数保持常量不变。透射损失补偿透射损失补偿实现起来非常困难，目前也没有这样的处理模块。在指数增益中，如果根据拟合出能量衰减对数曲线确定增益系数，则实际上以把透射损失考虑在内，所以目前在我们的处理中没有透射损失补偿这一步。Omega系统的地表一致性振幅补偿处理由地表一致性振幅补偿模块完成。该模块的基本思想是设法将地表因素对振幅的影响分离出来并予以归一化，因此它必须先在指定的时窗内计算各道的振幅统计量，有两种类型振幅统计量可供选择：均方根振幅和平均绝对值振幅。

均方根振幅p(i)

平均绝对值振幅p(i)其中，i表示道序号，a(j)是时窗内的样点值，N为时窗长度。设i点激发j点接收的道在指定的时窗内计算的振幅统计量为Aij，则Aij=Si·Rj·Gk·Ml

其中，Si=与炮点i有关的振幅变化分量Rj=与检波点j有关的振幅变化分量Gk=与第k个CMP点有关的振幅变化分量（k=(i+j)/2）Ml=与偏移距l有关的振幅变化分量（l=i-j）对上式两端取对数，得logAij=logSi+logRj+logGk+logMl

设法将各分量分解出来，并使分解后各分量之和与原来振幅量的误差能量

达到最小。由此可得以下高斯-赛德尔法的迭代公式：将每道的各振幅统计量分量分解出来后，再计算每道的振幅均衡因子。考虑到共中心点分量反映是地下反射点的信息，计算振幅均衡因子时不应使用。这时i点激发j点接收的道（偏移距为n=i-j）的振幅均衡因子Fij是其中，

分别是炮点分量、接收点分量和偏移距分量在整个测线上的几何平均。用各道的振幅均衡因子乘以输入数据就完成了地表一致性振幅补偿。

地表一致性振幅补偿由三个处理模块共同完成，它们分别是：地表一致性振幅补偿拾取、地表一致性振幅补偿分解和地表一致性振幅补偿应用，处理流程如下图所示：地震数据输入速度生成SCAC拾取SCAC分解输出SCAC分解数据时窗设计地震数据SCAC分解数据输入输入SCAC应用地表一致性振幅补偿输出第一步第二步地表一致性振幅补偿拾取（SCAC_PICK）的功能和主要参数功能计算振幅统计量。如前所述，有两种振幅统计量类型可供选择：均方根振幅和平均振幅。均方根振幅为常用方式，但它对存在大振幅噪音脉冲的资料不合适，因此在做地表一致性振幅补偿之前，应先做区域异常处理（ZAP）。大振幅噪音脉冲对平均振幅方式影响不大。在做地表一致性振幅补偿之前还应先做坏炮坏道去除，以免影响振幅统计量的计算。主要参数1、一般参数组（*GENERAL）

1）时窗文件名（WNDW_）该参数指定时窗文件名。该文件是时窗文件设计模块（WINDOW_DESIGN）的输出文件，它包含有各道用于计算振幅统计量的时窗时间。

2）时窗速度文件（WNDW_VEL_FILE）该参数指定速度文件名。前面所指定的时窗文件将用这个速度文件来做动校正。如果该参数缺省，则模块将使用包含在速度文件中的原始速度文件。在时窗设计时我们曾指定过这个速度文件名。

3）使用道头中的起始时间（USE_START_TIMES）该参数指定要不要用道头中的起始时间来计算本道第一个有效样点的位置。YES要用道头中的起始时间来计算这一道第一个有效样点的位置。这时第一个有效样点时间等于道头中的起始时间+起始延迟时间NO不用道头中的起始时间来计算这一道第一个有效样点的位置。这时第一个有效样点时间等于道头中的第一个样点定位时移（TIME_SHIFT_ALIGNMENT）+起始延迟时间缺省：YES

4）起始时间延迟（START_TIME_DELAY）该参数指定起始时间延迟值。它与前一个参数的值共同确定第一个有效样点的位置。这个参数值应适当大一些，以使起始时间位于直达波以下。缺省：200ms5）最小时窗百分比（MIN_WNDW_PCT）本参数指定时窗内有效样点的最小百分比。所谓有效样点指的是大于最小绝对振幅参数的样点。最小绝对值振幅参数将在稍后的振幅剔除参数组中讨论。缺省：506）去除道平衡因子（REMOVE_FACTOR）如果你在此项处理之前曾对输入数据做过道平衡处理，那么道头中存放有曾使用的道平衡因子。这里的参数确定，要不要把道平衡的作用去掉，即是否反道平衡。振幅统计量会因道平衡而改变，而道平衡因子一般不具备地表一致性，所以应去掉为好。YES在计算振幅统计量之前去除道平衡。NO在计算振幅统计量之前不去除道平衡。缺省：YES7）振幅统计量类型（*AMP_TYPE）本参数指定计算哪种振幅统计量类型。RMS_AMP均方根振幅MEAN_ABS_AMP平均方绝对值振幅缺省：RMS_AMP2、振幅值剔除限制参数组（AMP_REJECT）这一组参数设定一些限制，输入数据违反这些限制，本模块将不计算振幅统计量。

1）最小输入道平衡因子（MIN_TRC_BAL_FACTOR）这个参数指定一个最小输入道平衡因子。如果输入道的平衡因子小于该参数的值，模块在计算振幅统计量时将舍去这一道。缺省：-25.02）最大输入道平衡因子（MAX_TRC_BAL_FACTOR）这个参数指定一个最大输入道平衡因子。如果输入道的平衡因子大于该参数的值，模块在计算振幅统计量将舍去这一道。缺省：25.03）最大绝对值振幅排除（LARGE_ABS_AMP_ACTION）当输入数据的绝对值大于最大绝对值振幅参数或小于最小绝对值振幅参数时，怎么办？该参数提供下列选择：ZERO数据置为0，且不参加计算LARGEST数据的绝对值大于最大绝对值振幅参数，振幅统计量直接置为最大绝对值振幅参数值。数据的绝对值小于最小绝对值振幅参数，数据置为0，且不参加计算。NO_LIMIT不对数据的绝对值大小做检验，即无限制。COMPUTED根据振幅统计量类型而定：对均方根，该选择相当于ZERO；对平均绝对值振幅，该选择相当于NO_LIMIT缺省：COMPUTED4）最小绝对值振幅（SMALLEST_ABS_AMP）该参数指定最小绝对值振幅。绝对值小于该参数的值假定为0。取值：非负数。缺省：COMPUTED（根据振幅统计量类型而定：均方根振幅为1.0E-30；平均绝对值振幅为0)5）最大绝对值振幅（LARGE_ABS_AMP）该参数指定最大绝对值振幅。最大绝对值振幅排除将用这个值检验输入数据。取值：大于最小绝对值振幅的数。缺省：COMPUTED1.0E+30地表一致性振幅补偿分解的功能和主要参数功能

地表一致性振幅补偿分解（SCAC_DECOMPOSITION）模块是地表一致性振幅补偿处理的第二个模块，它将从地表一致性振幅补偿拾取模块输出的振幅统计量用高斯-赛德尔迭代法分解为共炮点、共接收点、共偏移距和共中心点分量。主要参数1、一般参数组（*GENERAL）1）平滑CMP分量的CMP点数（NUM_CMPS_SMOOTH）该参数指定CMP分量横向平滑滤波器的长度。CMP分量是地下介面反射性质的反映，它具有横向相似性。因此对其做横向平滑是有意义的，且可使迭代易于收敛。缺省：3

2、迭代参数组（*ITER）本组参数确定迭代次数、计算分量的名称及其次序等参数。1）最大迭代次数（MAX_NUM_ITER）最大迭代次数并非实际迭代次数。当某次迭代的均方根误差与上次迭代的均方根误差之差小于总均方根误差的确0.1%时，迭代自行终止。取值：大于0的整数缺省：62）CMP分量平滑的最后一次迭代（MAX_ITER_CMP_SMOOTH）该参数指定，从哪次迭代开始将不再做CMP分量平滑。缺省：最大迭代次数-13）估算分量次序（ORDER_OF_TERM）该参数指定计算分量的名称及其次序。名称必须是CMP、SOURCE、DETECTOR和OFFSET。如果本参数组没有提供，则隐含次序：SOURCE、DETECTOR和OFFSET，CMP分量不计算。没有指定的分量不计算，最多可指定四个分量，任何分量都不能重复指定。选择：‘CMP’共中心点分量‘SOURCE’共炮点分量‘DETECTOR’共检波点分量‘OFFSET’共偏移距分量缺省：None3、偏移距参数组（*OFFSET）本参数组定义如何按炮检距对数据道分组，以及要不要求偏移距分量时变（即各时窗求取自己的偏移距分量。如果在估算分量次序参数中选择了OFFSET，本参数组必须指定。如果在估算分量次序参数中未选择OFFSET，则本参数组将被忽略。1）偏移距分组数（NUM_OFFSET）该参数确定按偏移距大小将道分为多少个组，它的值应给得足够大，以使其能包含最大偏移距。缺省：#MAX_TRACE_NUM（ID头中的最大道号）2）偏移距分组间距（GROUP_INTERVAL）该参数确定偏移距分组的间距。通常它等于道间距。如果用缺省值（COMPUTED），则：对于2维资料，该参数取ID头中观测站之间的距离（STATION_SPACING），若观测站之间的距离不存在，就取ID头中共中心点间距（CMP_SPACING）的2倍。若共中心点间距也不存在，作业将输出错误信息后终止；对于3维资料，该参数取ID头中原始网格大小值（3D_PRIM_CELL_SIZE）。缺省：COMPUTED

3）第一个共偏移距组的偏移距（FIRST_OFFSET）该参数给定第一个共偏移距组的偏移距。以后各组的偏移距依次增加一个偏移距分组间距。缺省：偏移距分组间距（GROUP_INTERVAL）4、加权参数组（*PICK_WEIGHTING）这一组参数确定是否对各分量进行加权处理及加权处理参数。在存在噪声和坏脉冲的情况下，加权处理可提高模块的效能。加权的权因子按如下方法确定：

设上一次迭代得的的分解分量之和为：它与分解前的振幅统计量Aijh之间的距离定义为：

则加权因子F为：F=1.0当D<AR时F=1.0/(1.0+D-AR)当D>AR时其中AR为振幅比值参数。1）迭代次数序号（ITER_NUM）该参数指定从哪次迭代起开始做加权。建议至少填3。取值：大于等于2的数缺省：32）振幅比值（AMP_RADIO）用于确定加权因子F的振幅比值参数(AR)。选择：‘NO_WEIGHTING’不做加权大于等于1的数

缺省5

地表一致性振幅补偿应用模块的功能和主要参数功能

地表一致性振幅补偿应用模块（SCAC_APPLY）是地表一致性振幅补偿处理的最后一个模块，它将从地表一致性振幅补偿分解模块（SCAC_DECOMPOSITION）输出的振幅统计分量以地表一致性方式计算共炮点、共接收点、共偏移距道集的平衡因子，并应用于输入地震道。主要参数一般参数组（*GENERAL）1）平滑共中心点分量的相邻CMP点数（ADJ_CMPS_SMOOTH_GEOL）在计算比例因子之前可对共中心点分量做横向连续平滑。该参数确定平滑算子长度。0表示不做平滑；1表示做3点平滑；2表示做5平滑。如此等等。建议不做平滑。缺省：02）地震数据规格化（NORMALIZE）本参数确定要不要对来自地表一致性振幅补偿分解模块的各分量做规格化。所谓规格化就是将其比例放大或缩小到ID头中的均方根振幅（RMS_AMPLITUDE）水平。选择：‘NO’不规格化‘YES’要规格化缺省：YES2、SCAC分量项使用参数组（*TERM_APPLY）这一组参数指定哪些分量项用来计算并应用比例因子，建议四个分量都使用。1）使用地质（CMP）分量（TERM_APPLY_GEOL）该参数确定是否用地质（CMP）分量来计算并应用比例因子。‘YES’使用地质（CMP）分量校正‘NO’不使用地质（CMP）分量校正缺省：YES2）使用炮点分量（TERM_APPLY_SOURCE）该参数确定是否用炮点分量来计算并应用比例因子。‘YES’使用炮点分量校正‘NO’不使用炮点分量校正缺省：YES

3）使用检波点分量（TERM_APPLY_RECVR）该参数确定是否用检波点分量来计算并应用比例因子。‘YES’使用检波点分量校正‘NO’不使用检波点分量校正缺省：YES4）使用偏移距分量（TERM_APPLY_OFFSET）该参数确定是否用偏移距分量来计算并应用比例因子。‘YES’使用偏移距分量校正‘NO’不使用偏移距分量校正缺省：YES五剩余振幅补偿

在做完几何扩散补偿、吸收补偿和地表一致性振幅补偿之后，从理论上讲，地震数据的振幅随传播时间和距离的增大而减小以及因地表条件不一致引起的振幅变化已得到校正。但在实践中，由于没有详细的地下信息和准确的速度资料，以上校正并非完全准确，这就需要做进一步的振幅校正，即剩余振幅补偿。Omega系统的剩余振幅分析/补偿模块（RAAC）用统计平均方法来补偿剩余的振幅补偿量。

基本原理和计算方法因为在此之前，几何扩散、地层吸收、地表条件不一致及其它复杂因数引起的振幅衰减已经做了补偿，地震波的振幅无论在时间或空间上都不应有剧烈变化，如果有这种现象则可能存在剩余振幅补偿值。剩余振幅分析补偿模块用统计平均的方法计算振幅补偿因子，该因子在振幅较大的区域较小，在振幅较小的区域较大，用它们与响应区域的数据相乘，则各区域的能量水平就平衡了。剩余振幅分析/补偿模块分两步实现这一目标，第一步：计算振幅补偿因子，即剩余振幅分析；第二步：用振幅补偿因子与数据相乘，即剩余振幅补偿。第一步：剩余振幅分析剩余振幅分析的主要任务是计算振幅补偿因子。在计算振幅补偿因子之前，要先将数据集（例如炮集或叠加剖面）网格化，每个网格由若个道相同时窗的数据组成。下图是按炮－检距建立的二维维网格示意图，它将一个CMP道集的道沿空间方向分为四组，沿时间方向分为若干时窗，时窗是滑动的，后一个时窗和前一时窗有部分重叠。二按炮－检距建立的二维网格示意图而三维资料的分组则如下图所示说明：对于叠加剖面通常不分组，即每道都为一个组。

模块根据用户选择的剩余振幅分析类型参数计算各时窗的平均绝对值振幅（式1）或均方根振幅（式2）。

其中，N=非零样点个数WS=时窗起始时间＝时窗号*（N-1）X=样点值平均绝对值振幅或均方根振幅的倒数称为剩余振幅分析值（RAC）。它们将据时窗的中点时间放在适当的网格内。当一个网格内有多道时，存放哪一个RAC值，说明上未交代，可能是取均值，也可能是振幅量平均后再取倒数。这些剩余振幅分析值并不直接应用，还需做进一步处理。1、剩余振幅分析值的平滑考虑到地震波的振幅在时间或空间上都不应有剧烈变化，RAC值应当具有适度的时、空相关性。所以应当对所计算的RAC值进行适当的时、空平滑。一个点的平滑值由其周围若干点（奇数）的平均值计算：其中，j＝平滑起点号k＝平滑尾点号N＝平滑长度模块在进行时间方向的平滑时使用的独立平均（DiscreteAverage），而在进行空间方向的平滑时则使用的滚动平均（RunningAverage）。独立平均的意义是：每个参与平均的数都用其原来的值滚动平均的意义是：参与平均的数若已有平均值则使用平均后的值。如要对从1到6的6个数进行滚动平均，平滑长度为5,则平均过程如下：（1）第一个和最后一个数不变，仍是1和6;（2）1、2和3平均得第二个数的平均值：2’，（3）1、2’、3、4和5平均得第三个数的平均值：3’（4）2’、3’、4、5和6平均得第四个数的平均值：4’（5）4’、5和6平均得第五个数的平均值：5’符号，表示平滑后的数2、剩余振幅分析值的加权平均RAC值可能存在与合理的趋势不符的剧烈变化，尤其是在数据的顶部和底部。如果存在这种现象，则单用时间、空间的平滑效果并不令人满意，需要做加权平均。RAC值的加权平均的计算公式为：

其中，W（i）为加权系数。剩余振幅分析值计算完成，根据你的选择，模块将RAC输出到文本文件后终止或立即进入补偿应用步骤。在补偿应用步骤中，还需对对RAC值进行空间方向的线性内插，得到各个样点的剩余振幅因子，然后各样点与对应的剩余振幅因子相乘而完成剩余振幅补偿。主要参数1、一般参数组（*GENERAL）1）运行方式（RUN_TYPE）本参数指定运行RAAC的哪一步。选项：ANALY：剩余振幅分析APP：剩余振幅补偿BOTH：剩余振幅分析、补偿都做。缺省：ANALY2）RAAC函数文件名（RAAC_）本参数指定RAAC函数文件名，该文件是文本文件，它存放剩余振幅分析时所计算出的RAC值，模块在做剩余振幅补偿时将从中读出RAC值。取值：文件名，TEXT_FILES缺省：无缺省

2、剩余振幅分析参数组（*ANALY）1）第一级选排道头字（PRIM_SORT_LITERAL）本参数指定第一级选排的道头字，模块将按本道头字确定的道集格式进行剩余振幅分析，但输入数据不需预先按此格式做选排。缺省：无缺省2）第二级选排道头字（DECN_SORT_LITERAL）本参数指定第一级选排的道头字，模块将按本道头字确定的道第二级选排道头字，模块将按本道头字确定的道顺序行剩余振幅分析，但输入数据不需预先按此顺序做选排。缺省：不用（NOT_USED）3）使用第一级选排道头字的绝对值（ABS_PRIM_LITERAL）本参数确定在对地震道进行分组时，是否按第一级选排道头字的绝对值进行。NO：不按绝对值进行YES：按绝对值进行缺省：NO4）使用第二级选排道头字的绝对值（ABS_SECN_LITERAL）本参数确定在对地震道进行分组时，是否按第二级选排道头字的绝对值进行。NO：不按绝对值进行YES：按绝对值进行缺省：NO5）分析时窗起始（ANALY_WNDW_START）本参数指定分析时窗的起始位置。FIRST_SAMP：分析时窗从第一个非零样点开始START_TIME：分析时窗从道的起始时间开始缺省：无分析时窗起始起始点的选择与振幅补偿的应用无关，在做振幅补偿时，从道的起始时间到第一个非零样点所用的RAC值总是相当于此处选择START_TIME的结果。6）单个时窗选择（SINGLE_WNDW_OPTION）本参数指定是否只用一个分析时窗来计算RAC值。NO：分析时窗个数根据时窗长度和时窗滑动步进计算YES：只用一个分析时窗，每道只计算一个RAC值缺省：NO7）数据顶部时间平滑时窗个数（TIME_SMOOTH_BEGIN）本参数指定数据顶部参与平滑的RAC值个数。如果使用了单个时窗，则本参数忽略不计。该参数必须是奇数（如3、5、7等等）且不能大于道的长度/时窗滑动步进。如果该参数小于3，将不做平滑，建议使用3。取值：大于等于3的整数缺省：NOT_USED（不平滑）8）数据顶部时间平滑时窗个数（TIME_SMOOTH_BEGIN）本参数指定数据顶部参与平滑的RAC值个数。如果使用了单个时窗，则本参数忽略不计。该参数必须是奇数（如3、5、7等等）且不能大于道的长度/时窗滑动步进。如果该参数小于3，将不做平滑，建议使用3。取值：大于等于3的整数缺省：NOT_USED（不平滑）9）数据底部时间平滑时窗个数（TIME_SMOOTH_END）本参数指定数据底部参与平滑的RAC值个数。如果使用了单个时窗，则本参数忽略不计。该参数必须是奇数（如3、5、7等等）且不能大于道的长度/时窗滑动步进。数据中部的时间平滑时窗个数将根据数据顶部时间平滑时窗个数和数据底部时间平滑时窗个数进行线性内插。取值：大于等于3的整数缺省：NOT_USED（不平滑）10）第一级空间平滑宽度（RRIM_SMOOTH_WIDTH）本参数指定要以起进行平滑的第一级空间组数。该参数必须是奇数（如3、5、7等等）且不能大于第一级分组总数。如果该参数小于3，将不做平滑，建议使用3。取值：大于等于3的整数缺省：NOT_USED（不平滑）11）第二级空间平滑宽度（SECN_SMOOTH_WIDTH）本参数指定要以起进行平滑的第二空间组数。该参数必须是奇数（如3、5、7等等）且不能大于第二级分组总数。如果该参数小于3，将不做平滑，建议使用3。取值：大于等于3的整数缺省：NOT_USED（不平滑）12）时窗长度（WNDW_LENGTH）本参数指定分析时窗长度。如果使用了单个时窗，则本参数被忽略。一个时窗的RAC值放在该时窗的中点。如果不指定本参数，则缺省值为最大反射时间的十分之一。取值：大于0的数缺省：COMPUTED（最大反射时间的十分之一）13）时窗滑动步进（WNDW_ADVANCE）本参数指定时窗滑动步进。如果使用了单个时窗，则本参数被忽略。如果不指定本参数，则缺省值为半个时窗长度。时窗之间并非必须重叠。取值：大于0的数缺省：COMPUTED（半个时窗长度）14）期望输出改变（OTPUT_DECAY）本参数指定一个分贝数，所求的均方根振幅值将根据本参数予以改变。如果本参数是一个负数，则数据尾部的均方根振幅水平低于起始时间的均方根振幅水平，如果本参数是一个正数，则数据尾部的均方根振幅水平高于起始时间的均方根振幅水平，其变化量由本参数的值确定。如果不指定本参数，则均方根振幅保持不变。取值：正或负的数缺省：NOT_USED15）输出改变起始时间（OUTPUT_DECAY_START）本参数指定输出改变的起始时间，从该时间开始，道的振幅值将乘以内插后的期望输出改变参数。该时刻之前的振幅不变。如果该参数为零，则分贝斜坡区域从第一个非零样点开始。取值：大于等于0的数缺省：016）时间内包含有效样点的最少个数（MIN_WNDW_CONTENT）本参数指定时间内包含有效样点的最少个数，如果一个时窗内的非零样点个数小于本参数的值，则该时窗的RAC值不计算，而用其它方法估算，具体方法由后面的参数确定。取值：小于时窗长度的数缺省：COMPUTED（时窗步进的一半）17）有效样点不足的时窗外插（EXTRAP_WNDWS_）如果一个时窗内的非零样点个数太少，该时窗需要引入一个合理的RAC值，本参数确定用什麽方法引入该时窗的RAC值。COPY：拷贝先前的值。例如：如果最后一个时窗样点不足，就将该时窗之前一个时窗的RAC值拷贝到本时窗；如果第一个时窗样点不足，则将第二个时窗的RAC值拷贝到本时窗。LINEAR：线性外插。例如：如果最后一个时窗样点不足，就将该时窗之前的两个时窗的RAC值线性外插到本时窗；如果第一个时窗样点不足，则将其后两个时窗的RAC值线性外插到本时窗。缺省：COPY18）振幅分析类型（AMP_ANALY_TYPE）本参数指定振幅分析类型。RMS：均方根分析类型（建议使用）ABSOLUTE_VALUE：绝对值分析类型MEAN：算术平均分析类型缺省：RMS19）期望输出均方根水平（OUTPUT_RMS）本参数指定期望输出的均方根水平。如果缺省，模块将从道头字RMS_AMPLITUDE中寻找，如果找不到，则用2000。取值：大于0的整数缺省：COMPUTED3、分析时窗参数组（*ANALY_TIME_GATES）这一组参数以列表方式给出各控制点的分析时窗，控制点之间的时窗线性内插。<CNTRLPTBEGIN_TIMEEND_TIME6731804260428043404520456048892004280430043604540458041）控制点（CNTRLPT）该参数指定时窗控制点的位置。取值：递增的数组缺省：无2）时窗起始时间（BEGIN_TIME）该参数指定时窗起始时间，在一个控制点内必须是递增的。取值：递增的数组缺省：无3）时窗结束时间（BEGIN_TIME）该参数指定时窗结束时间，在一个控制点内必须是递增的。取值：递增的数组缺省：无4、数据顶部时间平滑权值参数组（*TIME_WEIGHTS_BEGIN）<WEIGHT121>1）加权值（WEIGHT）该参数指定数据顶部平滑的加权值。用户可根据数据顶部平滑时窗个数指定加权值的个数。如果给定N个加权值，则N个时窗使用相应的加权值。如果给出N个加权值，但只有N-1个时窗，则最后一组加权值不用。相反，如果只给一个加权值，而有N个时窗，则只有第一个时窗使用加权。取值：数组，个数大于等于顶部平滑时窗个数缺省：1.05、数据底部时间平滑权值参数组（*TIME_WEIGHTS_END）<WEIGHT313>1）加权值（WEIGHT）该参数指定数据底部的平滑的加权值。用此参数可实现时变加权。权值个数应与数据顶部时间平滑权值参数组指定的加权值个数相同。取值：数组，个数大于等于顶部平滑的加权值个数缺省：1.06、第一级空间平滑权值参数组（*PRIM_WEIGHTS）<WEIGHT121>1）加权值（WEIGHT）该参数指定第一级空间平滑的加权值。取值：数组，个数等于第一级空间平滑宽度缺省：1.07、第二级空间平滑权值参数组（*SECN_WEIGHTS）<WEIGHT121>

1）加权值（WEIGHT）该参数指定第二级空间平滑的加权值。取值：数组，个数等于第二级空间平滑宽度缺省：1.0六其它振幅处理

以上几种振幅补偿的方法都是从引起地震波能量改变的原因出发来考虑的，由于我们对其中地质因素并不十分掌握，补偿是否到位也就难以判断。通常只要远、近道，深、浅层以及炮与炮之间的能量看起来差别不大就认为可以了。这实际上包含了一种假设条件，即同一位置不同反射层位的反射能力相差无几，同一层位不同反射点的反射能力也基本一致。这一假设条件显然与实际的地质情况存在差距，但如果处理的资料仅用于构造解释也无不可。在Omega系统中，基于同一位置不同反射层位的反射能力相差无几这一条件的振幅处理模块有均方根增益；基于同一层位不同反射点反射的反射能力基本一致相差无几这一条件的振幅处理模块有时不变道均衡。

（一）均方根振幅增益（RMS_GAIN）

均方根振幅增益模块分别对各地震道独立进行处理，它将记录道分为若干时窗，并计算时窗内的均方根振幅：其中Xj是时窗内的样点值，N是样点个数。该时窗的均方根振幅增益因子定义为：式中的B为用户希望输出的均方根振幅输出水平。用A去与时窗的数据相乘就得到均方根振幅增益的结果：对一个时窗而言，如果我们计算均方根振幅增益后均方根振幅，则

由此可见，无论原始输入数据什么，输出数据的均方根振幅总是你指定的均方根振幅输出水平。我们来看两种极端的情况：1、当时窗短到只包含一个样点时，均方根振幅就是样点的振幅，增益结果是全道的振幅都变成了常数B，记录上的反射波特征也就不复存在了。这提醒我们，在做均方根增益时，时窗不能太短。2、当时窗长到包含整个地震道时，增益因子A为常数，增益结果是全道的均方根振幅为常数B，道内各点的振幅关系没有变化。处理效果就相当于时不变道均衡。当然模块实际上并不各时窗独立使用自已的增益因子，而是将各时窗的增益因子作为本时窗的中点的增益因子，时窗的中点之间的增益因子由相邻时窗中点增益因子线性内插获得。时窗定义方式时窗的划分对均方根振幅增益处理的效果影响很大。模块提供了两种确定时窗划分的方式：1、模块自动划分根据用户提供的第一个时窗长度和最大时窗长度参数，第一个时窗从道的起始时间开始，下一个时窗从上一个时窗结束点开始，长度为上一个时窗的2倍，直到等于最大时窗长度参数为止，以后的时窗长度保持不变。时窗划分和由此得到的增益曲线如下图：X2X4X4X时窗划分时窗中点H1H2H4H5增益因子增益曲线2、时窗由用户指定用户用参数表的方式指定各时窗的起始、结束时间，时窗之间可以重叠也可以留有空穴，但时窗中点时间必须严格递增，时窗划分和由此得到的增益曲线如下图：w1w2w3wn时窗划分时窗中点增益因子增益曲线H1H2H3Hn均方根增益模块的主要参数1、一般参数组（*GENERAL）1）模式（MODE）你是要做均方根增益还是要将已做的均方根增益去掉？APPLAY：做均方根增益REMOVE：去除已做的均方根增益2）期望的已做的均方根增益振幅水平（AMP）该参数应为大于0的常数，如果使用缺省值（COMPUTED），则为道头字RMS_AMPLITUDE）中的值，若道头中没有，缺省值为20003）增益因子镶边（TAPER）时窗中点之间的增益因子线性内插，最后一个时窗中点以后的增益因子保持常数。第一个时窗中点之前的增益因子则可以选择保持常数或镶边。NO：第一个时窗中点到第一个有效样点都用第一个时窗的增益因子YES：增益因子从第一个时窗中点到第一个有效样点线性降低到0。缺省：NO4）最小时窗百分比（MIN_WNDW_PCT）如果一个时窗内为0的样点太多（在做过道内切除的地方很可能发生），计算的增益因子会较大，它不仅对本时窗不利，还要影响相临时窗。本参数指定时窗内非0样点所占的最小百分比。如果时窗内非0样点所占的百分比小于本参数的值，则不计算增益因子，其增益因子由相邻有效时窗线性内插（或外插）。缺省：255）增益因子幂指数（AMP_POWER）对各时窗的增益因子做乘方放大或缩小，即将原有的增益因子A变为Ar，其中r为本参数的值。缺省：1（不变）2、用户定义时窗参数组（*USER_

DEFINED_WNDW）用下表形式给出时窗划分的参数<CNTRLPTWNDW_BEBIN_TIMEWNDW_END_TIME100600100090011001100150020050047001000140013001800>控制点道头字（LITERAL）控制点(CNTRLPT)的道头字缺省:SORT_LITERAL（分选道头字）

2）时窗时间起算道头字（BEGIN_WNDW_LITERAL）时窗的实际起始和结束时间将加上该道头字中的时间值。NOT_USED不用START_TIMESTART_TIME（起始时间）道头中的值缺省：NOT_USED3）控制点（CNTRLPT）4）时窗起始时间（

WNDW_BEBIN_TIME

）缺省：道头字EARLIEST_TIME（最早时间）中的值5）时窗结束时间（

WNDW_END_TIME

）缺省：道头字MAX_REFLECT_TIME（最大反射时间）中的值3、程序自动定义时窗参数组（*PROGRAM_

DEFINED_WNDW）用下表形式给出时窗划分的参数<CNTRLPTFIRST_ZONELONGEST_ZONE1002004100