一种基于频谱融合的地震波阻抗反演方法
【专利摘要】本发明涉及一种基于频谱融合的地震波阻抗反演方法，从测井数据提取背景场因子来控制反演波阻抗的整体趋势；将测井数据的低频成分与地震数据的有效频段成分进行频谱融合；并在地震数据的约束下将井中波阻抗曲线扩展到整个地震剖面。通过对测井数据和地震数据进行频谱融合，并引入背景场因子控制反演波阻抗的整体趋势，在地震数据的约束下将井中波阻抗曲线扩展到整个地震剖面。本发明比较完整的保留了地震反射的基本特征（断层、产状），不存在常规基于模型反演算法的多解性问题，能够明显的反映岩相、岩性的空间变化，具有计算速度快、稳定性好、计算精度高以及抗噪能力强的特点。
【专利说明】一种基于频谱融合的地震波阻抗反演方法
【技术领域】
[0001]本发明属于石油勘探【技术领域】，特别地涉及基于频谱融合的地震波阻抗反演方法。
【背景技术】
[0002]地震勘探技术是油气勘探中应用最为广泛的一种地球物理学方法，利用地震波在不同介质中传播的速度、振幅、频率、相位、波形等参数的变化来分析、预测油气储层分布范围及储层物性特征。随勘探开发程度逐渐深入，地震勘探已经从认识地下构造形态的构造勘探，逐渐发展成直接应用地震信息判断岩性、分析岩相、定量计算岩层物性参数的岩性勘探。这些复杂沉积环境下的储层通常表现为具有厚度薄、物性高度非均质、有效储层规模小、分布分散、岩石物理关系复杂、储层岩性差异小等特征。波阻抗是反映储集岩储层特征的一种重要物性参数，与振幅等反映界面性质的参数相比，波阻抗具有更高的储层横向预测能力，因而地震波阻抗反演技术是岩性油气勘探的重要技术手段。近几十年来地震波阻抗反演技术获得了长足的进展，如道积分法、模型法、随机反演、谱反演等。这些方法均存在一定适用条件，无法适应薄储层等复杂岩性储层的储层预测工作。
[0003]道积分反演方法的特点是无需钻井控制，反演结果直接反应了地层波阻抗的相对变化，但是该方法在数据处理过程中不能使用地质或测井资料进行约束控制，并且受地震数据固有频带的限制，分辨率低，无法适于薄层解释的需要。求取的相对波阻抗和绝对波阻抗不能用于储层参数的定量计算。
[0004]递推反演是利用反摺积从地震剖面求取反射系数剖面，进而将反射系数剖面转换为波阻抗剖面，在岩性相对稳定条件下能较好的反应储层的物性变化。在利用反褶积从地震数据求取反射系数时，选择有效的反摺积方法是递推反演的关键，如频域反褶积、LI模反褶积、最大似然反褶积、最小熵反褶积等。但是反褶积后得到的地震道振幅实际上还不是反射系数，而是相当于反射系数再褶积一个剩余子波，这个剩余子波在浅部主频较高，在深部较低。主频低的在积分得到的波阻抗数值偏大，会导致深层产生偏大的波阻抗。
[0005]基于模型的波阻抗反演技术将测井数据与地震数据有机结合，反演的波阻抗分辨率高，适于薄层油气藏的的精细描述。主要取决于初始模型与实际地质情况的符合程度，在钻井资料较多的情况下应用效果较好，多解性是基于模型反演的固有特性，在初始模型建立不准的情况下通常导致反演结果的错误。
[0006]在勘探开发前期，通常使用道积分和递推反演来定性的进行岩相分析确定地层的沉积体系，结合钻井资料揭示的储层特征进行横向预测，确定评价井位。在油田开发阶段，由于钻井资料较为丰富，可以结合地震数据建立用于反演的初始模型，利用基于模型的波阻抗反演来进行定量的反演，实现储层精细描述。但是常规叠后波阻抗反演技术存在以下缺点:分辨率低，反演结果具有多解性，抗噪能力差，不适于储层参数定量计算及薄储层预测，计算速度慢不适于实时解释，应用受到限制等缺点。
【发明内容】
[0007]为了解决上述问题，本发明的目的就是针对现有地震反演技术在储层预测中存在的不足，围绕地震数据的有限带宽及地震反演中存在的多解性、反演计算速度慢、无法处理大数据体的实时解释等问题，提供一种基于频谱融合的地震波阻抗反演方法，充分利用地震数据和测井数据的有效频带，提高反演方法的抗噪能力，将得到的波阻抗信息用于地震资料的岩性参数解释，便于精细储层描述，为寻找岩性(和地层)圈闭油气藏和非常规油气藏(煤层气、页岩气等)提供可靠的资料。
[0008]本发明是通过以下技术方案实现的:
[0009]一种基于频谱融合的地震波阻抗反演方法，所述方法从测井数据提取背景场因子来控制反演波阻抗的整体趋势，将测井数据的低频成分与地震数据的有效频段成分进行频谱融合。并在地震数据的约束下将井中波阻抗曲线扩展到整个地震剖面。
[0010]所述方法包括以下步骤:
[0011]8)对井中波阻抗曲线Al在最小二乘意义下进行线性拟合，得到背景场因子LAI ；
[0012]所述的井中波阻抗曲线Al通过下面公式得到:
[0013]AI=Vp.P(I)
[0014]上式中，vp为纵波速度，P为密度。vp和P由常规声波测井和密度测井资料得到。
[0015]设变量y与变量t之间的函数关系满足下式:
[0016]y = a+b.t(2)
[0017]求出a，b使得以下公式达到最小值:
[0018]
【权利要求】
1.一种基于频谱融合的地震波阻抗反演方法，其特征在于，从测井数据提取背景场因子来控制反演波阻抗的整体趋势；将测井数据的低频成分与地震数据的有效频段成分进行频谱融合；并在地震数据的约束下将井中波阻抗曲线扩展到整个地震剖面。
2.根据权利要求1所述的一种基于频谱融合的地震波阻抗反演方法，其特征在于，包括以下步骤: 1)对井中波阻抗曲线Al在最小二乘意义下进行线性拟合，得到背景场因子LAI； 所述的井中波阻抗曲线Al通过下面公式得到: AI=Vp.P(I) 上式中，VP为纵波速度，P为密度；VP和P由常规声波测井和密度测井资料得到； 设变量I与变量t之间的函数关系满足下式: y=a+b.t(2) 求出a，b使得以下公式达到最小值:
3.根据权利要求2所述的一种基于频谱融合的地震波阻抗反演方法，其特征在于，所述的步骤I)中所述的采用常`规最优化算法，包括牛顿法、共轭梯度法。
【文档编号】G01V1/28GK103792573SQ201210418887
【公开日】2014年5月14日 申请日期:2012年10月26日 优先权日:2012年10月26日
【发明者】白俊雨, 岳承琪, 宋志翔, 肖云飞 申请人:中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院