本发明涉及油气勘探开发领域，尤其涉及油气储层密度的预测方法。
背景技术：
油气储层的密度是最基本的岩石性质之一，它反映了储层岩石骨架的矿物成分、孔隙度和孔隙流体类型及饱和度，能够揭示关于烃类，尤其是气态烃的发育位置及其饱和度等重要信息，是一种重要的储层弹性参数。因此，储层密度对岩性和流体识别具有重要意义。
国内外应用实例表明，相对于速度等弹性参数，储层密度可以有效预测烃源岩的toc含量，可以更好的反映储层的质量，有效预测储层中的优质储层、储层中流体的产能。储层密度可以在储层非均质性强的地区预测优质烃源岩的发育区；在储集层均值或相对均值的地区研究烃源岩的toc含量。尤其是在描述和评价气态烃源岩发育区中具有巨大的潜力。因此，发展新方法指导实际勘探中储层密度的预测，可以为油气地震勘探提供更加可靠的技术支持。
现有关于反演储层密度的方法是基于模型的叠前地震反演方法，采用贝叶斯反演来获取储层的密度属性，该方法中密度项对反射系数的贡献率低，使得反演准确度较低，不能适用于实际需求。
技术实现要素：
本发明要解决的技术问题是：提出一种油气储层密度高可靠性预测方法，提高储层预测与流体识别的可靠性。
为了解决上述技术问题，本发明所采用的技术方案是：一种油气储层密度高可靠性预测方法，
设在界面两侧介质弹性变化较小且入射角在临界角范围内假设条件下，利用纵横波阻抗表示的反射系数线性近似方程为：
公式(1)中，ip为介质纵波阻抗，is为介质横波阻抗，ρ为介质密度，θ为入射角和透射角的平均；
所述油气储层密度高可靠性预测方法包括以下步骤：
第一步：将公式(1)写作：
公式(2)矩阵形式为
公式(3)中，rm为纵横波阻抗和密度的反射系数，a，b，c为对应的系数；
对公式(3)中待预测参数rm进行解耦，建立独立方程：
第二步：建立目标泛函，计算待定系数d，
公式(5)中，rpp为利用fatti反射系数线性近似方程计算的井旁道反射系数，
得到系数矩阵d写作为：
第三步：利用公式(4)计算得到待预测参数rm，再通过道积分计算得到待预测参数rm的绝对值。
本发明技术方案带来的有益效果是：本发明的储层密度高可靠性预测方法，能够有效准确预测储层密度信息。
附图说明
图1为反射系数随密度的变化示意图。
图2某工区二维地震数据剖面示意图。
图3使用本发明方法反演计算得到的储层密度剖面示意图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明，以使本领域的技术人员可以更好的理解本发明并能予以实施，但所举实施例不作为对本发明的限定。
现有储层密度反演方法如下：
在界面两侧介质弹性变化较小且入射角在临界角范围内假设条件下，fatti在aki-richard近似方程的基础上，发展了利用纵横波阻抗表示的反射系数线性近似方程：
式中，ip为介质纵波阻抗，is为介质横波阻抗，ρ为介质密度，θ为入射角和透射角的平均。以上式为正演方程，在贝叶斯反演框架下，对纵横波阻抗和密度进行反演预测。然而，在小角度入射情况下，密度项对反射系数的贡献率较低，因此很难通过反演得到准确的密度参数。为了解决这种问题，本案提出了一种新的高可靠性预测方法，并通过实例分析验证了新方法在预测密度的可行性。
本方法是基于fatti反射系数线性近似方程，建立独立方程计算弹性参数，其计算过程分为以下几步：
第一步：为简便表示，将式(1)写作
其矩阵形式为
其中，rm为纵横波阻抗和密度的反射系数，a，b，c为对应的系数。
对式(3)中待预测参数rm进行解耦，建立独立方程
第二步：建立目标泛函，计算待定系数dm
其中rpp为利用fatti反射系数线性近似方程计算的井旁道反射系数，
得到系数矩阵d，可写作为
第三步：利用地震反射系数rpp与系数矩阵d计算得到待预测参数，即利用式(4)计算得到待预测参数rm，再通过道积分计算得到待预测参数rm的绝对值。
图1表示密度项对反射系数的贡献率，从图中可以看出本发明提出的方法在小角度入射情况下，密度项对反射系数的贡献率有很大的提高，因此，反演结果也更加可靠。图1中，横坐标：incicentangle(degree)入射角(度)，纵坐标：rpp纵波反射系数。
图2是某工区二维地震数据，图3是使用本方法预测计算得到的储层密度，从图中可以看出本发明预测的密度剖面与测井数据较好的吻合，并且在含油气储层位置显低值异常，对于密度预测以及含油气储层预测具有良好的结果。图2中，纵坐标：time(s)时间(秒)，色标：seismic(-)地震数据(-)。图3中，纵坐标：time(s)时间(秒)，色标：density密度(千克/立方米)。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例，本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换，均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。