CN104536042B - 一种二维叠后地震资料振幅补偿方法及装置 - Google Patents

Abstract

本发明涉及一种二维叠后地震资料振幅补偿方法及装置，包括：在二维叠后地震资料中提取时间域地震子波和地震有效频带范围，并确定地震资料振幅减弱区域的时间域钻井反射系数；时间域地震子波与时间域钻井反射系数进行褶积，获得时间域钻井合成地震记录道；对二维叠后地震资料中时间域井旁地震道经过变换获得井旁地震道振幅谱，对时间域钻井合成地震记录道经过变换获得钻井合成地震记录道振幅谱；利用井旁地震道振幅谱和钻井合成地震记录道振幅谱获得地震有效频带范围内的振幅补偿因子；确定二维叠后地震资料中振幅减弱区域的边界；利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿。

Description

一种二维叠后地震资料振幅补偿方法及装置

技术领域

本发明涉及叠后地震资料处理技术领域，特别涉及一种二维叠后地震资料振幅补偿方法及装置。

背景技术

伴随着国内外油气勘探程度的不断深化，以地震资料做为主要解释平台的综合地质解释方法得到了广泛的应用，这对叠后地震资料的成像品质也提出了更高的要求。

在目前油气勘探过程中，地震解释工作仍以叠后地震资料二维解释为基础，受诸如火成岩、碳酸盐岩等一些高波阻抗特殊岩性体的影响，其下伏地层地震反射特征失真，尤其振幅衰减严重，成像品质变差，这对二维构造解释、层位追踪以及波阻抗反演等都产生了严重的影响。针对叠后地震资料中特殊岩性体造成的下伏地层振幅的异常衰减，虽然根据现有的叠前地震资料成像处理技术能够进行相应的处理，提高成像品质，但在某些地区，这些特殊岩性体的规模及其分布是没有规律的，鉴于叠前地震资料成像处理技术的工作量和成本，使得该技术不能针对每一个特殊岩性体进行处理，从而导致某些特殊岩性体其下伏地层振幅异常衰减仍然存在；而目前基于叠后地震资料来对特殊岩性体造成的下伏地层振幅异常衰减进行处理的方法还没有被提出。

发明内容

为解决现有技术的问题，本发明提出一种二维叠后地震资料振幅补偿方法及装置。

为实现上述目的，本发明提供了一种二维叠后地震资料振幅补偿方法，该方法包括：

在二维叠后地震资料中，从振幅不减弱区域通过钻井结合井旁地震道提取时间域地震子波和地震有效频带范围，并确定地震资料振幅减弱区域的时间域钻井反射系数；

所述时间域地震子波与所述时间域钻井反射系数进行褶积，获得时间域钻井合成地震记录道；

对二维叠后地震资料中时间域井旁地震道经过变换获得井旁地震道振幅谱，对时间域钻井合成地震记录道经过变换获得钻井合成地震记录道振幅谱；

利用所述井旁地震道振幅谱和所述钻井合成地震记录道振幅谱，通过最小二乘法获得地震有效频带范围内的振幅补偿因子；

通过最小二乘法确定二维叠后地震资料中振幅减弱区域的边界；

利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿。

优选地，所述利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿的步骤包括：

以边界处所在地震道为中心，以等间隔数、且间隔数从小到大依次向边界处所在地震道两边取得的地震道振幅谱分别运用最小二乘法确定两边地震道的第二振幅差异因子，直至前一次获得的第二振幅差异因子不小于当前获得的第二振幅差异因子为止；以当前获得第二振幅差异因子对应的两个地震道中靠近振幅减弱区域一侧的地震道为基准地震道，对这两个地震道之间的其他任一地震道振幅谱和所述基准地震道的振幅谱通过最小二乘法获得第三振幅差异因子；振幅补偿因子减去第三振幅差异因子再加1的结果值，这两个地震道之间的任一地震道振幅谱乘以所述结果值，获得的计算结果经过反变换得到边界处这两个地震道之间的最终振幅补偿结果；其中，所述间隔数为2n-1，n为自然数；

在振幅减弱区域范围内，除边界处已进行振幅补偿的地震道外其他地震道振幅谱乘以振幅补偿因子，获得的计算结果经过反变换得到最终振幅补偿结果。

优选地，所述利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿的步骤还包括：

对振幅补偿结果的边界道以及边界道两侧相邻的地震道进行平滑处理。

优选地，所述时间域钻井反射系数通过校正后的声波测井曲线和密度测井曲线确定。

优选地，所述确定二维叠后地震资料中振幅减弱区域的边界的步骤包括：

确定两地震道的间隔道数；并在地震有效频带范围内，对二维叠后地震资料中符合间隔道数的任两地震道对应的地震道振幅谱通过最小二乘法法获得第一振幅差异因子；其中，所述间隔道数为奇数，且满足在取定的间隔道数范围内第一振幅差异因子的极大值是唯一的条件；

根据第一振幅谱差异因子确定第一振幅差异因子的极大值，以通过观测二维叠后地震资料中的地震剖面获得的补偿边界位置为约束，利用所述极大值确定补偿边界处的边界。

为实现上述目的，本发明还提供了一种二维叠后地震资料振幅补偿装置，该装置包括：

预处理单元，用于在二维叠后地震资料中，从振幅不减弱区域通过钻井结合井旁地震道提取时间域地震子波和地震有效频带范围，并确定地震资料振幅减弱区域的时间域钻井反射系数；

褶积单元，用于所述时间域地震子波与所述时间域钻井反射系数进行褶积，获得时间域钻井合成地震记录道；

时频域转换单元，用于对二维叠后地震资料中时间域井旁地震道经过变换获得井旁地震道振幅谱，对时间域钻井合成地震记录道经过变换获得钻井合成地震记录道振幅谱；

振幅补偿因子获取单元，用于利用所述井旁地震道振幅谱和所述钻井合成地震记录道振幅谱，通过最小二乘法获得地震有效频带范围内的振幅补偿因子；

边界确定单元，用于通过最小二乘法确定二维叠后地震资料中振幅减弱区域的边界；

振幅补偿单元，用于利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿。

优选地，所述振幅补偿单元包括：

第一振幅补偿模块，用于以边界处所在地震道为中心，以等间隔数、且间隔数从小到大依次向边界处所在地震道两边取得的地震道振幅谱分别运用最小二乘法确定两边地震道的第二振幅差异因子，直至前一次获得的第二振幅差异因子不小于当前获得的第二振幅差异因子为止；以当前获得第二振幅差异因子对应的两个地震道中靠近振幅减弱区域一侧的地震道为基准地震道，对这两个地震道之间的其他任一地震道振幅谱和所述基准地震道的振幅谱通过最小二乘法获得第三振幅差异因子；振幅补偿因子减去第三振幅差异因子再加1的结果值，这两个地震道之间的任一地震道振幅谱乘以所述结果值，获得的计算结果经过反变换得到边界处这两个地震道之间的最终振幅补偿结果；其中，所述间隔数为2n-1，n为自然数；

第二振幅补偿模块，用于在振幅减弱区域范围内，除边界处已进行振幅补偿的地震道外其他地震道振幅谱乘以振幅补偿因子，获得的计算结果经过反变换得到最终振幅补偿结果。

优选地，所述振幅补偿单元还包括：

平滑模块，用于对振幅补偿结果的边界以及边界两侧相邻的地震道进行平滑处理。

优选地，所述预处理单元通过校正后的声波测井曲线和密度测井曲线确定时间域钻井反射系数。

优选地，所述边界确定单元包括：

第一振幅差异因子获取模块，用于确定两地震道的间隔道数；并在地震有效频带范围内，对二维叠后地震资料中符合间隔道数的任两地震道对应的地震道振幅谱通过最小二乘法法获得第一振幅差异因子；其中，所述间隔道数为奇数，且满足在取定的间隔道数范围内第一振幅差异因子的极大值是唯一的条件；

边界道确定模块，用于根据第一振幅谱差异因子确定第一振幅差异因子的极大值，以通过观测二维叠后地震资料中的地震剖面获得的补偿边界位置为约束，利用所述极大值确定补偿边界处的边界。

上述技术方案具有如下有益效果：

本技术方案解决了在叠后地震资料中高波阻抗特殊岩性体下伏地层振幅异常衰减的问题，在纵向上不影响地震反射轴相对振幅强弱关系的前提下，对振幅异常衰减部分进行了准确的补偿，在高波阻抗特殊岩性体边界处振幅补偿过渡平缓，符合地层反射特征，经实际资料验证，该发明应用效果良好。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明提出的一种二维叠后地震资料振幅补偿方法流程图；

图2为本发明提出的一种二维叠后地震资料振幅补偿装置框图；

图3为本实施例的特殊岩性体二维叠后地震剖面图；

图4为本实施例的振幅恢复所用的统一的地震子波示意图；

图5为本实施例的W1井合成地震记录道振幅谱与W1井井旁地震道振幅谱的对比图；

图6为本实施例的特殊岩性体二维叠后地震补偿剖面图；

图7为本实施例的振幅补偿后振幅谱对比图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本技术方案的工作原理是：从二维叠后地震资料入手，利用振幅衰减带的测井资料计算得出的地层真实反射系数对该振幅衰减带的地震资料振幅进行有效补偿。使得在保证振幅补偿结果准确可靠前提下，较叠前地震资料成像处理而言，减小了工作量，提高了处理的灵活性及工作效率。

基于上述工作原理，本发明提出一种二维叠后地震资料振幅补偿方法，如图1所示。该方法包括：

步骤101)：在二维叠后地震资料中，从振幅不减弱区域通过钻井结合井旁地震道提取时间域地震子波和地震有效频带范围，并确定地震资料振幅减弱区域的时间域钻井反射系数；

对于井旁地震道时间域地震子波来说，在同一区域相同层位内的特殊岩性体范围之外，即地震振幅正常区域，选择时深关系标定较好的多口钻井结合井旁地震道来提取一个确定的统一的平均地震子波，该地震子波是进行振幅补偿的基础，其最大振幅值与二维地震资料振幅值和测井计算的反射系数有关。

步骤102)：所述时间域地震子波与所述时间域钻井反射系数进行褶积，获得时间域钻井合成地震记录道；

以钻井结合井旁地震道提取的确定的统一的平均地震子波为前提运用褶积公式计算钻井合成地震记录道。设钻井反射系数为r(t)，地震子波为w(t)，则钻井合成地震记录道x(t)为前面二者的褶积，即：

x(t)＝w(t)*r(t) (1)

上式中“*”代表褶积。钻井反射系数r(t)是通过校正后的声波测井曲线和密度测井曲线计算得出的。设声波测井曲线AC单位为us/m，密度测井曲线DEN单位为g/cm³，首先用AC曲线求取速度v，即：

v＝10⁶/AC (2)

此处速度v的单位是m/s。设钻井采样点为i，DEN测井曲线值为ρ，AC表示声波测井曲线值。则反射系数R为：

$R_i=\frac{{\mathrm{\ρ}}_{i+1}{v}_{i+1}-{\mathrm{\ρ}}_i{v}_i}{{\mathrm{\ρ}}_{i+1}{v}_{i+1}+{\mathrm{\ρ}}_i{v}_i}---\left(3\right)$

为了和地震道相匹配，通过时深转换将深度域反射系数R转换为时间域反射系数r(t)。

步骤103)：对二维叠后地震资料中时间域井旁地震道经过变换获得井旁地震道振幅谱，对时间域钻井合成地震记录道经过变换获得钻井合成地震记录道振幅谱；

快速傅里叶变换/傅里叶变换方法被广泛应用于信号的时频分析当中，该技术应用快速傅丽叶变换将地震及钻井合成记录由时间域转换为频率域，以利于振幅补偿因子的计算。

步骤104)：利用所述井旁地震道振幅谱和所述钻井合成地震记录道振幅谱，通过最小二乘法获得地震有效频带范围内的振幅补偿因子；

钻井反射系数直接由声波测井和密度测井计算而来，而测井的分辨率远远高于地震分辨率(0.5m以下)，因此由测井计算的反射系数在纵向上不受特殊岩性体的影响能够代表地层真实的反射系数，在时深关系标定准确的前提下，由其计算得出的合成地震记录道的振幅谱能够真实反应地层的反射特征。在不受特殊岩性体影响的地震资料振幅正常的区域，经过统一的平均地震子波结合钻井反射系数求取的合成地震记录道与该位置在特定频带内地震道振幅谱在形状和幅值上是基本相同的；而在地震资料振幅异常衰减的区域，合成地震记录道与该位置地震道在特定频带内振幅谱形状相同，但幅值出现差异。因此利用钻井合成地震记录道与井旁地震道的对比就可以补偿振幅异常衰减的部分，这里用到的是最小二乘法。

最小二乘法是一种数学优化技术，它通过最小化误差的平方和寻找数据的最佳函数匹配。设补偿因子为Q，井旁地震道振幅谱为X(f)，钻井反射系数振幅谱为W(f)，则在地震有效频带范围内，根据最小二乘法使得下式的值要达到最小：

(W₁-Q·X₁)²+(W₂-Q·X₂)²+…+(W_n-Q·X_n)² (3)

式中n表示地震有效频带内，振幅谱中的采样点个数。若使上式达到最小值，则Q得满足使上式的导数为0，即：

X₁·(W₁-Q·X₁)+X₂·(W₂-Q·X₂)+…+X₂·(W₂-Q·X₂)＝0 (4)

求解上式(4)，即得补偿因子Q值。在二维地震资料上，特殊岩性体范围内钻井有时不只一口，计算出的Q值也不只一个，因为特殊岩性体阻抗值不变，所以认为其对下伏地层的影响是一致的，即补偿因子Q是固定的，所以Q值取多口井计算结果的平均值。

步骤105)：通过最小二乘法确定二维叠后地震资料中振幅减弱区域的边界；

若已知特殊岩性体边界，此处则无需计算。在特殊岩性体边界不确定的情况下，则运用最小二乘法来计算振幅补偿边界。在特殊岩性体边界处，下伏地层的地震振幅有明显的减弱，若以一定的间隔来计算二维地震资料两道之间振幅差异，则在边界处振幅差异值达到最大，而且远离边界处振幅差异很小，利用这点，便可经过计算准确求出振幅补偿边界。基于上述工作原理，在特殊岩性体边界不确定的情况下，确定边界道的方法为：

确定两地震道的间隔道数；并在地震有效频带范围内，对二维叠后地震资料中符合间隔道数的任两地震道对应的地震道振幅谱通过最小二乘法法获得第一振幅差异因子；其中，所述间隔道数为奇数，且满足在取定的间隔道数范围内第一振幅差异因子的极大值是唯一的条件；根据第一振幅谱差异因子确定第一振幅差异因子的极大值，以通过观测二维叠后地震资料中的地震剖面获得的补偿边界位置为约束，利用所述极大值确定补偿边界处的边界道。

步骤106)：利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿。

以边界处所在地震道为中心，以等间隔数、且间隔数从小到大依次向边界处所在地震道两边取得的地震道振幅谱分别运用最小二乘法确定两边地震道的第二振幅差异因子，直至前一次获得的第二振幅差异因子不小于当前获得的第二振幅差异因子为止；以当前获得第二振幅差异因子对应的两个地震道中靠近振幅减弱区域一侧的地震道为基准地震道，对这两个地震道之间的其他任一地震道振幅谱和所述基准地震道的振幅谱通过最小二乘法获得第三振幅差异因子；振幅补偿因子减去第三振幅差异因子再加1的结果值，这两个地震道之间的任一地震道振幅谱乘以所述结果值，获得的计算结果经过反变换得到边界处这两个地震道之间的最终振幅补偿结果；

在振幅减弱区域范围内，除边界处已进行振幅补偿的地震道外其他地震道振幅谱乘以振幅补偿因子，获得的计算结果经过反变换得到最终振幅补偿结果。

为使得补偿效果更好，对振幅补偿结果的边界道以及边界道两侧相邻的地震道进行平滑处理。

基于上述工作原理，本发明还提供了一种二维叠后地震资料振幅补偿装置，如图2所示。该装置包括：

预处理单元201，用于在二维叠后地震资料中，从振幅不减弱区域通过钻井结合井旁地震道提取时间域地震子波和地震有效频带范围，并确定地震资料振幅减弱区域的时间域钻井反射系数；

褶积单元202，用于所述时间域地震子波与所述时间域钻井反射系数进行褶积，获得时间域钻井合成地震记录道；

时频域转换单元203，用于对二维叠后地震资料中时间域井旁地震道经过变换获得井旁地震道振幅谱，对时间域钻井合成地震记录道经过变换获得钻井合成地震记录道振幅谱；

振幅补偿因子获取单元204，用于利用所述井旁地震道振幅谱和所述钻井合成地震记录道振幅谱，通过最小二乘法获得地震有效频带范围内的振幅补偿因子；

边界确定单元205，用于通过最小二乘法确定二维叠后地震资料中振幅减弱区域的边界；

振幅补偿单元206，用于利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿。

在某油田某层系存在高波阻抗火成岩体(特殊岩性体)，其下伏R地层受火成岩高波阻抗的影响，振幅值衰减明显，如图3所示。由于地震资料振幅的衰减，对R地层构造解释、砂体追踪、储层预测等都会产生影响，进而影响对R地层的油气勘探，因为对火成岩体下伏的R地层地震资料振幅进行有效的补偿就显得尤为必要了。

运用本发明技术，首先在同一区域R地层内的火成岩体范围之外，即地震振幅正常区域，选择时深关系标定较好的多口钻井结合井旁地震道来提取一个确定的统一的平均地震子波，如图4所示。

在火成岩体范围之内，选择钻穿R地层，通过时深关系标定较好的钻井(W1井)时间域反射系数，运用褶积公式结合前述平均地震子波计算W1井的合成地震记录道；然后通过快速傅里叶变换计算W1井合成地震记录道振幅谱及W1井井旁地震道振幅谱，通过振幅谱进行比较，选择5～28Hz做为计算振幅补偿因子的频带范围，如图5所示。图5中虚线表示W1井井旁地震道振幅谱，实线表示W1井合成地震记录道振幅谱，频带范围选择的原则是：在井旁地震道振幅谱有效频带范围内。

运用最小二乘法计算振幅补偿因子，结果为：1.293。

运用最小二乘法在地震资料有效频带范围内计算该地震剖面R地层地震振幅差异因子；采用分步法确定差异因子极值点，最终确定地震资料振幅补偿边界。

利用快速傅丽叶反变换及均值平滑法完成最终R地层振幅补偿结果的计算，如图6所示，在火成岩体下伏R地层地震资料振幅得到了有效的恢复，为验证振幅补偿结果的准确性，取原地震资料R地层第5～45道(正常振幅地震道)做频谱分析，取原地震资料R地层第75～115道(振幅衰减地震道)做频谱分析，取振幅补偿后地震资料R地层第75～115道(振幅补偿地震道)做频谱分析，其振幅谱对比结果如图7所示。在地震有效频带范围内，受火成岩影响的地震道振幅值明显小于未受火成岩影响的地震道振幅值，经补偿后地震道振幅值与未受火成岩影响的地震道振幅值基本相当，这说明该发明实现了对地震资料振幅衰减部分的有效补偿，补偿效果良好。

以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种二维叠后地震资料振幅补偿方法，其特征在于，该方法包括：

在二维叠后地震资料中，从振幅不减弱区域通过钻井结合井旁地震道提取时间域地震子波和地震有效频带范围，并确定地震资料振幅减弱区域的时间域钻井反射系数；

所述时间域地震子波与所述时间域钻井反射系数进行褶积，获得时间域钻井合成地震记录道；

对二维叠后地震资料中时间域井旁地震道经过变换获得井旁地震道振幅谱，对时间域钻井合成地震记录道经过变换获得钻井合成地震记录道振幅谱；

利用所述井旁地震道振幅谱和所述钻井合成地震记录道振幅谱，通过最小二乘法获得地震有效频带范围内的振幅补偿因子；

通过最小二乘法确定二维叠后地震资料中振幅减弱区域的边界；

利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿；

其中，所述利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿的步骤包括：

以边界处所在地震道为中心，以等间隔数、且间隔数从小到大依次向边界处所在地震道两边取得的地震道振幅谱分别运用最小二乘法确定两边地震道的第二振幅差异因子，直至前一次获得的第二振幅差异因子不小于当前获得的第二振幅差异因子为止；以当前获得第二振幅差异因子对应的两个地震道中靠近振幅减弱区域一侧的地震道为基准地震道，对这两个地震道之间的其他任一地震道振幅谱和所述基准地震道的振幅谱通过最小二乘法获得第三振幅差异因子；振幅补偿因子减去第三振幅差异因子再加1的结果值，这两个地震道之间的任一地震道振幅谱乘以所述结果值，获得的计算结果经过反变换得到边界处这两个地震道之间的最终振幅补偿结果；其中，所述间隔数为2n-1，n为自然数；

在振幅减弱区域范围内，除边界处已进行振幅补偿的地震道外其他地震道振幅谱乘以振幅补偿因子，获得的计算结果经过反变换得到最终振幅补偿结果。

2.如权利要求1所述的方法，其特征在于，所述利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿的步骤还包括：

对振幅补偿结果的边界道以及边界道两侧相邻的地震道进行平滑处理。

3.如权利要求1或2所述的方法，其特征在于，所述时间域钻井反射系数通过校正后的声波测井曲线和密度测井曲线确定。

4.如权利要求1或2所述的方法，其特征在于，所述确定二维叠后地震资料中振幅减弱区域的边界的步骤包括：

确定两地震道的间隔道数；并在地震有效频带范围内，对二维叠后地震资料中符合间隔道数的任两地震道对应的地震道振幅谱通过最小二乘法获得第一振幅差异因子；其中，所述间隔道数为奇数，且满足在取定的间隔道数范围内第一振幅差异因子的极大值是唯一的条件；

根据第一振幅谱差异因子确定第一振幅差异因子的极大值，以通过观测二维叠后地震资料中的地震剖面获得的补偿边界位置为约束，利用所述极大值确定补偿边界处的边界。

5.一种二维叠后地震资料振幅补偿装置，其特征在于，该装置包括：

预处理单元，用于在二维叠后地震资料中，从振幅不减弱区域通过钻井结合井旁地震道提取时间域地震子波和地震有效频带范围，并确定地震资料振幅减弱区域的时间域钻井反射系数；

褶积单元，用于所述时间域地震子波与所述时间域钻井反射系数进行褶积，获得时间域钻井合成地震记录道；

时频域转换单元，用于对二维叠后地震资料中时间域井旁地震道经过变换获得井旁地震道振幅谱，对时间域钻井合成地震记录道经过变换获得钻井合成地震记录道振幅谱；

振幅补偿因子获取单元，用于利用所述井旁地震道振幅谱和所述钻井合成地震记录道振幅谱，通过最小二乘法获得地震有效频带范围内的振幅补偿因子；

边界确定单元，用于通过最小二乘法确定二维叠后地震资料中振幅减弱区域的边界；

振幅补偿单元，用于利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿；

其中，所述振幅补偿单元包括：

第一振幅补偿模块，用于以边界处所在地震道为中心，以等间隔数、且间隔数从小到大依次向边界处所在地震道两边取得的地震道振幅谱分别运用最小二乘法确定两边地震道的第二振幅差异因子，直至前一次获得的第二振幅差异因子不小于当前获得的第二振幅差异因子为止；以当前获得第二振幅差异因子对应的两个地震道中靠近振幅减弱区域一侧的地震道为基准地震道，对这两个地震道之间的其他任一地震道振幅谱和所述基准地震道的振幅谱通过最小二乘法获得第三振幅差异因子；振幅补偿因子减去第三振幅差异因子再加1的结果值，这两个地震道之间的任一地震道振幅谱乘以所述结果值，获得的计算结果经过反变换得到边界处这两个地震道之间的最终振幅补偿结果；其中，所述间隔数为2n-1，n为自然数；

第二振幅补偿模块，用于在振幅减弱区域范围内，除边界处已进行振幅补偿的地震道外其他地震道振幅谱乘以振幅补偿因子，获得的计算结果经过反变换得到最终振幅补偿结果。

6.如权利要求5所述的装置，其特征在于，所述振幅补偿单元还包括：

平滑模块，用于对振幅补偿结果的边界以及边界两侧相邻的地震道进行平滑处理。

7.如权利要求5或6所述的装置，其特征在于，所述预处理单元通过校正后的声波测井曲线和密度测井曲线确定时间域钻井反射系数。

8.如权利要求5或6所述的装置，其特征在于，所述边界确定单元包括：

第一振幅差异因子获取模块，用于确定两地震道的间隔道数；并在地震有效频带范围内，对二维叠后地震资料中符合间隔道数的任两地震道对应的地震道振幅谱通过最小二乘法法获得第一振幅差异因子；其中，所述间隔道数为奇数，且满足在取定的间隔道数范围内第一振幅差异因子的极大值是唯一的条件；

边界道确定模块，用于根据第一振幅谱差异因子确定第一振幅差异因子的极大值，以通过观测二维叠后地震资料中的地震剖面获得的补偿边界位置为约束，利用所述极大值确定补偿边界处的边界。