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

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

【技术实现步骤摘要】
一种二维叠后地震资料振幅补偿方法及装置
本专利技术涉及叠后地震资料处理
，特别涉及一种二维叠后地震资料振幅补偿方法及装置。
技术介绍
伴随着国内外油气勘探程度的不断深化，以地震资料做为主要解释平台的综合地质解释方法得到了广泛的应用，这对叠后地震资料的成像品质也提出了更高的要求。在目前油气勘探过程中，地震解释工作仍以叠后地震资料二维解释为基础，受诸如火成岩、碳酸盐岩等一些高波阻抗特殊岩性体的影响，其下伏地层地震反射特征失真，尤其振幅衰减严重，成像品质变差，这对二维构造解释、层位追踪以及波阻抗反演等都产生了严重的影响。针对叠后地震资料中特殊岩性体造成的下伏地层振幅的异常衰减，虽然根据现有的叠前地震资料成像处理技术能够进行相应的处理，提高成像品质，但在某些地区，这些特殊岩性体的规模及其分布是没有规律的，鉴于叠前地震资料成像处理技术的工作量和成本，使得该技术不能针对每一个特殊岩性体进行处理，从而导致某些特殊岩性体其下伏地层振幅异常衰减仍然存在；而目前基于叠后地震资料来对特殊岩性体造成的下伏地层振幅异常衰减进行处理的方法还没有被提出。
技术实现思路
为解决现有技术的问题，本专利技术提出一种二维叠后地震资料振幅补偿方法及装置。为实现上述目的，本专利技术提供了一种二维叠后地震资料振幅补偿方法，该方法包括：在二维叠后地震资料中，从振幅不减弱区域通过钻井结合井旁地震道提取时间域地震子波和地震有效频带范围，并确定地震资料振幅减弱区域的时间域钻井反射系数；所述时间域地震子波与所述时间域钻井反射系数进行褶积，获得时间域钻井合成地震记录道；对二维叠后地震资料中时间域井旁地震道经过变换获得井旁地震道振幅谱，对时间域钻井合成地震记录道经过变换获得钻井合成地震记录道振幅谱；利用所述井旁地震道振幅谱和所述钻井合成地震记录道振幅谱，通过最小二乘法获得地震有效频带范围内的振幅补偿因子；通过最小二乘法确定二维叠后地震资料中振幅减弱区域的边界；利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿。优选地，所述利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿的步骤包括：以边界处所在地震道为中心，以等间隔数、且间隔数从小到大依次向边界处所在地震道两边取得的地震道振幅谱分别运用最小二乘法确定两边地震道的第二振幅差异因子，直至前一次获得的第二振幅差异因子不小于当前获得的第二振幅差异因子为止；以当前获得第二振幅差异因子对应的两个地震道中靠近振幅减弱区域一侧的地震道为基准地震道，对这两个地震道之间的其他任一地震道振幅谱和所述基准地震道的振幅谱通过最小二乘法获得第三振幅差异因子；振幅补偿因子减去第三振幅差异因子再加1的结果值，这两个地震道之间的任一地震道振幅谱乘以所述结果值，获得的计算结果经过反变换得到边界处这两个地震道之间的最终振幅补偿结果；其中，所述间隔数为2n-1，n为自然数；在振幅减弱区域范围内，除边界处已进行振幅补偿的地震道外其他地震道振幅谱乘以振幅补偿因子，获得的计算结果经过反变换得到最终振幅补偿结果。优选地，所述利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿的步骤还包括：对振幅补偿结果的边界道以及边界道两侧相邻的地震道进行平滑处理。优选地，所述时间域钻井反射系数通过校正后的声波测井曲线和密度测井曲线确定。优选地，所述确定二维叠后地震资料中振幅减弱区域的边界的步骤包括：确定两地震道的间隔道数；并在地震有效频带范围内，对二维叠后地震资料中符合间隔道数的任两地震道对应的地震道振幅谱通过最小二乘法法获得第一振幅差异因子；其中，所述间隔道数为奇数，且满足在取定的间隔道数范围内第一振幅差异因子的极大值是唯一的条件；根据第一振幅谱差异因子确定第一振幅差异因子的极大值，以通过观测二维叠后地震资料中的地震剖面获得的补偿边界位置为约束，利用所述极大值确定补偿边界处的边界。为实现上述目的，本专利技术还提供了一种二维叠后地震资料振幅补偿装置，该装置包括：预处理单元，用于在二维叠后地震资料中，从振幅不减弱区域通过钻井结合井旁地震道提取时间域地震子波和地震有效频带范围，并确定地震资料振幅减弱区域的时间域钻井反射系数；褶积单元，用于所述时间域地震子波与所述时间域钻井反射系数进行褶积，获得时间域钻井合成地震记录道；时频域转换单元，用于对二维叠后地震资料中时间域井旁地震道经过变换获得井旁地震道振幅谱，对时间域钻井合成地震记录道经过变换获得钻井合成地震记录道振幅谱；振幅补偿因子获取单元，用于利用所述井旁地震道振幅谱和所述钻井合成地震记录道振幅谱，通过最小二乘法获得地震有效频带范围内的振幅补偿因子；边界确定单元，用于通过最小二乘法确定二维叠后地震资料中振幅减弱区域的边界；振幅补偿单元，用于利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿。优选地，所述振幅补偿单元包括：第一振幅补偿模块，用于以边界处所在地震道为中心，以等间隔数、且间隔数从小到大依次向边界处所在地震道两边取得的地震道振幅谱分别运用最小二乘法确定两边地震道的第二振幅差异因子，直至前一次获得的第二振幅差异因子不小于当前获得的第二振幅差异因子为止；以当前获得第二振幅差异因子对应的两个地震道中靠近振幅减弱区域一侧的地震道为基准地震道，对这两个地震道之间的其他任一地震道振幅谱和所述基准地震道的振幅谱通过最小二乘法获得第三振幅差异因子；振幅补偿因子减去第三振幅差异因子再加1的结果值，这两个地震道之间的任一地震道振幅谱乘以所述结果值，获得的计算结果经过反变换得到边界处这两个地震道之间的最终振幅补偿结果；其中，所述间隔数为2n-1，n为自然数；第二振幅补偿模块，用于在振幅减弱区域范围内，除边界处已进行振幅补偿的地震道外其他地震道振幅谱乘以振幅补偿因子，获得的计算结果经过反变换得到最终振幅补偿结果。优选地，所述振幅补偿单元还包括：平滑模块，用于对振幅补偿结果的边界以及边界两侧相邻的地震道进行平滑处理。优选地，所述预处理单元通过校正后的声波测井曲线和密度测井曲线确定时间域钻井反射系数。优选地，所述边界确定单元包括：第一振幅差异因子获取模块，用于确定两地震道的间隔道数；并在地震有效频带范围内，对二维叠后地震资料中符合间隔道数的任两地震道对应的地震道振幅谱通过最小二乘法法获得第一振幅差异因子；其中，所述间隔道数为奇数，且满足在取定的间隔道数范围内第一振幅差异因子的极大值是唯一的条件；边界道确定模块，用于根据第一振幅谱差异因子确定第一振幅差异因子的极大值，以通过观测二维叠后地震资料中的地震剖面获得的补偿边界位置为约束，利用所述极大值确定补偿边界处的边界。上述技术方案具有如下有益效果：本技术方案解决了在叠后地震资料中高波阻抗特殊岩性体下伏地层振幅异常衰减的问题，在纵向上不影响地震反射轴相对振幅强弱关系的前提下，对振幅异常衰减部分进行了准确的补偿，在高波阻抗特殊岩性体边界处振幅补偿过渡平缓，符合地层反射特征，经实际资料验证，该专利技术应用效果良好。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描本文档来自技高网...

【技术保护点】
一种二维叠后地震资料振幅补偿方法，其特征在于，该方法包括：在二维叠后地震资料中，从振幅不减弱区域通过钻井结合井旁地震道提取时间域地震子波和地震有效频带范围，并确定地震资料振幅减弱区域的时间域钻井反射系数；所述时间域地震子波与所述时间域钻井反射系数进行褶积，获得时间域钻井合成地震记录道；对二维叠后地震资料中时间域井旁地震道经过变换获得井旁地震道振幅谱，对时间域钻井合成地震记录道经过变换获得钻井合成地震记录道振幅谱；利用所述井旁地震道振幅谱和所述钻井合成地震记录道振幅谱，通过最小二乘法获得地震有效频带范围内的振幅补偿因子；通过最小二乘法确定二维叠后地震资料中振幅减弱区域的边界；利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿。

【技术特征摘要】
1.一种二维叠后地震资料振幅补偿方法，其特征在于，该方法包括：在二维叠后地震资料中，从振幅不减弱区域通过钻井结合井旁地震道提取时间域地震子波和地震有效频带范围，并确定地震资料振幅减弱区域的时间域钻井反射系数；所述时间域地震子波与所述时间域钻井反射系数进行褶积，获得时间域钻井合成地震记录道；对二维叠后地震资料中时间域井旁地震道经过变换获得井旁地震道振幅谱，对时间域钻井合成地震记录道经过变换获得钻井合成地震记录道振幅谱；利用所述井旁地震道振幅谱和所述钻井合成地震记录道振幅谱，通过最小二乘法获得地震有效频带范围内的振幅补偿因子；通过最小二乘法确定二维叠后地震资料中振幅减弱区域的边界；利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿；其中，所述利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿的步骤包括：以边界处所在地震道为中心，以等间隔数、且间隔数从小到大依次向边界处所在地震道两边取得的地震道振幅谱分别运用最小二乘法确定两边地震道的第二振幅差异因子，直至前一次获得的第二振幅差异因子不小于当前获得的第二振幅差异因子为止；以当前获得第二振幅差异因子对应的两个地震道中靠近振幅减弱区域一侧的地震道为基准地震道，对这两个地震道之间的其他任一地震道振幅谱和所述基准地震道的振幅谱通过最小二乘法获得第三振幅差异因子；振幅补偿因子减去第三振幅差异因子再加1的结果值，这两个地震道之间的任一地震道振幅谱乘以所述结果值，获得的计算结果经过反变换得到边界处这两个地震道之间的最终振幅补偿结果；其中，所述间隔数为2n-1，n为自然数；在振幅减弱区域范围内，除边界处已进行振幅补偿的地震道外其他地震道振幅谱乘以振幅补偿因子，获得的计算结果经过反变换得到最终振幅补偿结果。2.如权利要求1所述的方法，其特征在于，所述利用地震有效频带范围内的振幅补偿因子、二维叠后地震资料中振幅减弱区域的边界对二维叠后地震资料进行振幅补偿的步骤还包括：对振幅补偿结果的边界道以及边界道两侧相邻的地震道进行平滑处理。3.如权利要求1或2所述的方法，其特征在于，所述时间域钻井反射系数通过校正后的声波测井曲线和密度测井曲线确定。4.如权利要求1或2所述的方法，其特征在于，所述确定二维叠后地震资料中振幅减弱区域的边界的步骤包括：确定两地震道的间隔道数；并在地震有效频带范围内，对二维叠后地震资料中符合间隔道数的任两地震道对应的地震道振幅谱通过最小二乘法获得第一振幅差异因子；其中，所述间隔道数为奇数，且满足在取定的间隔道数范围内第一振幅差异因子的极大值是唯一的条件；根据第一振幅谱差异因子确定第一振幅差异因子的极大值，以通过观测二维叠后地震资料中的地震剖面获得的补偿边界位置为约束，利用所述极大值确定补偿边界处的边界...

【专利技术属性】
技术研发人员：姚军，李双文，倪长宽，张晶，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11