一种面向叠前反演的道集振幅分频补偿方法
【技术领域】
[0001] 本发明涉及油气藏勘探技术领域，尤其是一种面向叠前反演的道集振幅分频补偿 方法。
【背景技术】
[0002] 近年来，随着岩性油气藏勘探需求的日益增多，与之配套的地震资料储层预测技 术，特别是叠前属性分析和叠前弹性参数反演技术得到了长足的发展。这些技术在岩性预 测和流体检测中得到广泛应用，其预测精度除了受反演方法本身影响之外，还受输入数据 质量的影响。其中，叠前道集作为叠前属性分析及叠前弹性参数反演的输入数据，其质量对 反演结果精度的影响较大。
[0003] 目前，本领域技术人员普遍认识到叠前道集的去噪、子波一致性、入射角范围、剩 余动校和远道畸变等预处理对叠前反演结果的影响很大。但是，这些预处理主要集中在改 善资料信噪比和分辨率等方面，缺乏对保真性，例如道集动力学信息---振幅、频率、相位 补偿的研宄，尤其对宽角度的研宄更是少见。然而，叠前弹性参数反演要得到稳定的密度 项，需要尽可能宽的入射角度。此外，由于地震波在地下介质中传播时会受吸收衰减影响， 因此地震波的高频能量会有损失，振幅、频率和相位会发生畸变。虽然在现有技术中，在叠 前道集的前期处理中也会对各种振幅进行补偿，但是这些补偿方法只是限于纵向上的几何 扩散及吸收补偿，而没有考虑道集横向上的补偿。熊定钰等在《保持地震记录叠前av0属性 的噪声衰减方法》一文中提出了基于道集随入射角变化满足抛物线规律的补偿方法，但是 其在计算抛物线系数时，采用了零入射角数据进行计算，由于近道去面波、偏移时近远偏移 距覆盖次数不均匀等原因，近-远道能量较弱，所有小角度原始数据不适合进行参数的求 取，因此具有一定的局限性。综上所述，现有的技术方法虽然对道集的振幅有所改善，但是 都是常规的单一整体补偿方法，没有考虑不同频率具有不同的吸收衰减特征及av0特征， 还是无法满足叠前地震反演的需求，不能为叠前弹性参数反演提供可靠的数据依据。
【发明内容】
[0004] 针对上述问题，本发明提出了一种面向叠前反演的道集振幅分频补偿方法。
[0005] 该补偿方法，包括以下步骤：
[0006] 预处理步骤，对原始的叠前道集资料进行包括去噪、剩余动校和转入射角道集的 预处理；
[0007] 分频处理步骤，对预处理后的叠前道集资料进行分频处理；
[0008] 参数提取步骤，从不同频率的叠前道集资料中提取符合正演av0特征趋势的入射 角以及对应的振幅值；
[0009] 拟合计算步骤，对提取的不同频率道集的入射角和振幅值进行av0抛物线拟合， 计算抛物线方程的各个系数；
[0010] 能量补偿步骤，根据拟合得到的抛物线方程重新确定不同频率道集中不同入射角 时的反射系数，根据重新确定的反射系数对近道、远端能量弱的道进行分频振幅补偿。
[0011] 根据本发明的实施例，上述预处理步骤中，基于同相叠加的射线束聚焦进行去噪 处理。
[0012] 根据本发明的实施例，基于同相叠加的射线束聚焦的去噪处理包括以下步骤：
[0013] 对于一个确定的偏移距，首先定义炮点和检波点的邻域范围，将此范围内的道集 抽取出来为一个射线束；
[0014] 在射线束内选取标志反射层，在标志反射层附近选取时窗进行叠加形成标准道， 与射线束内各道进行相关分析；
[0015] 如果相关时差超过1/4子波长度对应的时间，进行剔除；
[0016] 将满足同相叠加原理的各道进行叠加。
[0017] 根据本发明的实施例，上述剩余动校处理是采用反射轴自动追踪的方法，对给定 时窗内的各叠前远道与近、中道求取相关系数，获得相对于零偏移距的时移量，然后逐步滑 动时窗获得整个道集的时移量，将道集拉平。
[0018] 根据本发明的实施例，上述分频处理步骤，基于小波变换对预处理后的叠前道集 资料进行分频处理。
[0019] 根据本发明的实施例，上述参数提取步骤中，优选地从不同频率的叠前道集资料 中提取符合正演av0特征趋势的中角度的入射角以及对应的振幅值。
[0020] 根据本发明的实施例，上述中角度的入射角优选为10~25度。
[0021] 根据本发明的实施例，上述拟合计算步骤中，根据提取的入射角和振幅进行av0 抛物线拟合，通过下式确定抛物线方程的系数r、w、v;
【主权项】
1. 一种面向叠前反演的道集振幅分频补偿方法，其特征在于，包括 预处理步骤，对原始的叠前道集资料进行包括去噪、剩余动校和转入射角道集的预处 理； 分频处理步骤，对预处理后的叠前道集资料进行分频处理； 参数提取步骤，从不同频率的叠前道集资料中提取符合正演avo特征趋势的入射角以 及对应的振幅值； 拟合计算步骤，对提取的不同频率道集的入射角和振幅值进行avo抛物线拟合，计算 抛物线方程的各个系数； 能量补偿步骤，根据拟合得到的抛物线方程重新确定不同频率道集中不同入射角时的 反射系数，根据重新确定的反射系数对近道及远道能量弱的道进行分频振幅补偿。
2. 如权利要求1所述的道集振幅分频补偿方法，其特征在于，所述预处理步骤中，基于 同相叠加的射线束聚焦进行去噪处理。
3. 如权利要求2所述的道集振幅分频补偿方法，其特征在于，基于同相叠加的射线束 聚焦的去噪处理包括以下步骤： 对于一个确定的偏移距，首先定义炮点和检波点的邻域范围，将此范围内的道集抽取 出来为一个射线束； 在射线束内选取标志反射层，在标志反射层附近选取时窗进行叠加形成标准道，与射 线束内各道进行相关分析； 在时窗范围内，将相关时差超过1/4子波长度的道集予以剔除； 将满足同相叠加原理的各道进行叠加。
4. 如权利要求1~3任意一项所述的道集振幅分频补偿方法，其特征在于： 所述剩余动校处理是采用反射轴自动追踪的方法，对给定时窗内的各叠前远道与近、 中道求取相关系数，获得相对于零偏移距的时移量，然后逐步滑动时窗获得整个道集的时 移量，将道集拉平。
5. 如权利要求1~3任意一项所述的道集振幅分频补偿方法，其特征在于： 所述分频处理步骤，基于小波变换对预处理后的叠前道集资料进行分频处理。
6. 如权利要求1~3任意一项所述的道集振幅分频补偿方法，其特征在于，进一步地， 所述参数提取步骤中，从不同频率的叠前道集资料中提取符合正演avo特征趋势的中 角度的入射角以及对应的振幅值。
7. 如权利要求6所述的道集振幅分频补偿方法，其特征在于： 所述中角度的入射角的范围是10~25度。
8. 如权利要求1所述的道集振幅分频补偿方法，其特征在于，所述拟合计算步骤中，根 据提取的入射角和振幅进行avo抛物线拟合，通过下式确定抛物线方程的系数r、w、v ;
式中，xi=sin20i，0i是第i个入射角；yi是第i个入射角对应的中间变量。
9.如权利要求8所述的道集振幅分频补偿方法，其特征在于： 拟合计算步骤中，采用最小二乘法计算抛物线方程的各个系数。
【专利摘要】本发明提出了一种面向叠前反演的道集振幅分频补偿方法。该方法对原始的叠前道集资料进行包括去噪、剩余动校和转入射角道集的预处理；然后对预处理后的叠前道集资料进行基于小波变换的分频处理；然后以正演avo特征趋势为约束，从不同频率的叠前道集资料中提取中角度的入射角以及对应的振幅值；对提取的不同频率道集的入射角和振幅值进行avo抛物线拟合，计算抛物线方程的各个系数；最后根据拟合得到的抛物线方程重新确定不同频率道集中不同入射角时的反射系数，根据重新确定的反射系数对近道、远端能量弱的道进行分频振幅补偿，使振幅随入射角的变化满足抛物线的规律，满足叠前反演方法对输入数据的要求，从而提高叠前弹性参数反演的可靠性和稳定性。
【ipc分类】g01v1-36
【公开号】cn104820242
【申请号】cn201510142504
【发明人】陈祖庆, 李金磊, 屈大鹏, 李文成, 杨鸿飞, 王明飞, 肖秋红
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司勘探分公司
【公开日】2015年8月5日
【申请日】2015年3月27日