曲线之间的误差。在获得每一口井作为验证井时的误差后，平均验证误差最低的曲线组合以及最佳匹配函数即为T0C计算模型:
[0038]TOC = f (Log1； Log2,…，Logn) (1)
[0039]式中，Logl为第i条测井曲线或岩石弹性参数曲线，f为TOC与测井曲线或岩石弹性参数曲线之间的函数关系(即利用最小二乘法确定的最佳匹配函数)，所述函数包括线性关系、多项式关系、指数关系、对数关系等，不同地区可根据实际资料情况利用最小二乘法确定曲线数目以及与TOC相关系数最高的函数关系。
[0040](3)敏感曲线数据体预测:利用叠前地震反演技术可获得纵波速度、横波速度、密度数据体，利用这三个数据体可进一步计算T0C计算模型所需测井曲线或岩石弹性参数数据体:
[0041]叠前反演输入数据包括:①纵波速度、横波速度、密度测井曲线，用于为反演提供低频约束模型；②叠前地震道集，由于地震道集资料品质一般，需经过随机噪音衰减、同相轴拉平、AV0特征补偿等处理来提高信噪比和分辨率，恢复泥页岩储层真实AV0响应；③基于叠后数据通过构造解释获得的地震层位数据，主要用于低频模型的建立。利用上述数据通过叠前地震反演可获得纵波速度、横波速度、密度数据体，可进一步计算得到T0C计算模型所需测井曲线或岩石弹性参数数据体。
[0042](4)泥页岩T0C预测:利用式(1)所示的T0C计算模型以及步骤(3)所述反演数据体计算得到T0C数据体，通过提取目的层泥页岩段T0C数据体平面图，便可了解泥页岩T0C横向变化规律，并可进一步以某个T0C值为门槛值提取优质泥页岩厚度，完成泥页岩发育情况的描述。
[0043]本发明方法的一个实施例如下:
[0044]①地质、测井、地震数据准备:地质数据包括录井、试气以及泥页岩岩心总有机碳含量地球化学分析化验数据；测井数据包括偶极子声波测井或全波列声波测井资料，至少包括纵波时差、横波时差、密度测井曲线以及针对泥页岩的T0C测井解释数据；地震资料包括经过常规叠后/叠前偏移处理的成果数据以及经过随机噪音衰减、同相轴拉平、AV0特征补偿等面向叠前地震反演的处理手段处理后的叠前地震道集数据。
[0045]②对泥页岩岩心T0C测试数据以及测井曲线和岩石弹性参数曲线进行统计分析，优选与T0C数据相关程度高的测井曲线或岩石弹性参数曲线，建立T0C计算模型；
[0046]③利用叠前地震资料上更加丰富的能够反映地层岩性、油气的振幅、频率等信息，开展叠前地震反演获得纵波速度、横波速度、密度数据体，并进一步计算T0C计算模型所需曲线数据体。
[0047]④利用步骤②所述计算模型以及步骤③所述曲线数据体计算T0C数据体，完成对泥页岩T0C横向变化情况的预测。
[0048]图2为利用T0C计算模型计算得到的拟合T0C曲线与实测T0C曲线之间的交会图。T0C计算模型为:
[0049]T0C = -5.83 X 10 3 X μ Ρ 2+2.94 X 10 6 X Vs2-5.58 (式 2)
[0050]式中，Vs为横波速度，单位m/s ； μ Ρ为剪切模量X密度，单位GPaXg/cm3。
[0051]图3为利用式(2)及叠前地震反演结果计算得到的T0C连井剖面。其中，利用叠前地震技术反演可获得纵波速度(vp)、横波速度(vs)以及密度数据体(P);利用μ P =(P XVS)2可计算得到μ Ρ数据体；将μ Ρ和％数据体代入式(2)便可计算得到T0C数据体。
[0052]图4为利用式2及叠前地震反演结果计算并编制的T0C含量等值线图，可见由JY1井、JY2井、JY3井、JY4井所控制的构造主体部位T0C均大于2%，说明研究区富含有机质泥页岩发育，页岩气成藏地质条件优越。
[0053]图5为本发明预测的优质泥页岩厚度图，提取厚度时以T0C>1%作为门槛值，由泥页岩厚度图可知研究区JY1井、JY2井、JY井3、JY4井所控制的构造主体部位泥页岩厚度分布稳定且均大于50米，较厚的泥页岩是商业性高产的基础。
[0054]上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述【具体实施方式】所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的意义。
【主权项】
1.一种泥页岩总有机碳含量地震预测方法，其特征在于:所述方法包括: (1)地质、测井、地震数据准备； (2)对泥页岩岩心总有机碳含量地球化学分析化验数据以及测井曲线和岩石弹性参数曲线进行统计分析，优选与总有机碳含量测试数据含量相关程度高的测井曲线或岩石弹性参数曲线，建立总有机碳含量计算模型； (3)利用地震反演技术获得总有机碳含量计算模型所需的曲线反演数据体； (4)基于步骤(2)中的总有机碳含量计算模型以及步骤(3)中的曲线反演数据体计算泥页岩总有机碳含量数据体。2.根据权利要求1所述的泥页岩总有机碳含量地震预测方法，其特征在于:所述步骤(1)中的地质数据包括录井、试气以及泥页岩岩心总有机碳含量地球化学分析化验数据； 所述测井数据包括偶极子声波测井或全波列声波测井资料，至少包括纵波时差、横波时差、密度测井曲线以及经泥页岩岩心总有机碳含量地球化学分析化验数据刻度的总有机碳含量测井解释曲线； 所述地震数据包括经叠后或叠前偏移处理的成果数据以及经过面向叠前地震反演的处理手段处理后的叠前地震道集数据；所述面向叠前地震反演的处理手段至少包括随机噪音衰减、同相轴拉平、AVO特征补偿。3.根据权利要求2所述的泥页岩总有机碳含量地震预测方法，其特征在于:所述步骤(2)是这样实现的: 第一步，通过交会分析确定与总有机碳含量测井解释曲线相关系数最高的测井曲线或岩石弹性参数曲线，即第一曲线； 第二步，依次将第一曲线与其它任一测井曲线或岩石弹性参数曲线组成曲线对，针对每一曲线对利用最小二乘法确定的最佳匹配函数计算拟合总有机碳含量曲线，并计算其与总有机碳含量测井解释曲线的相关系数，相关系数最高的则为第二曲线； 第三步，按照第二步所述方法继续寻找第三曲线、第四曲线直到第N曲线； 利用盲井验证的方法确定曲线数目，即依次使用每一 口井作为验证井，利用其它井确定的曲线数目及最佳匹配函数计算验证井的拟合总有机碳含量曲线，计算该曲线与总有机碳含量测井解释曲线之间的误差；在获得每一口井作为验证井时的误差后，平均验证误差最低的曲线组合以及最佳匹配函数即为总有机碳含量计算模型:TOC = f (Log1； Log2,..., Logn) (1) 式中，TOC为总有机碳含量；Logl为第i条测井曲线或岩石弹性参数曲线；f为T0C与测井曲线或岩石弹性参数曲线之间的函数关系；所述函数关系包括线性关系、多项式关系、指数关系或对数关系；不同地区根据实际资料情况选用与总有机碳含量相关程度最高的函数关系。4.根据权利要求1所述的泥页岩总有机碳含量地震预测方法，其特征在于:所述步骤(3)是这样实现的: 输入纵波速度、横波速度、密度测井曲线、叠前地震道集和基于叠后数据通过构造解释获得的地震层位数据； 利用叠前地震反演技术获得纵波速度、横波速度、密度数据体，并进一步利用这三个数据体计算总有机碳含量计算模型所需的测井曲线和岩石弹性参数数据体，即获得所述曲线反演数据体。5.根据权利要求4所述的泥页岩总有机碳含量地震预测方法，其特征在于:所述步骤(4)是这样实现的:将步骤(3)得到的各个测井曲线和岩石弹性参数数据体代入到式(1)所示的总有机碳含量计算模型中，计算得到总有机碳含量数据体。
【专利摘要】本发明提供了一种泥页岩总有机碳含量地震预测方法，属于石油天然气勘探领域。本方法包括：(1)地质、测井、地震数据准备；(2)对泥页岩岩心总有机碳含量地球化学分析化验数据以及测井曲线和岩石弹性参数曲线进行统计分析，优选与总有机碳含量测试数据含量相关程度高的测井曲线或岩石弹性参数曲线，建立总有机碳含量计算模型；(3)利用地震反演技术获得总有机碳含量计算模型所需的曲线反演数据体；(4)基于步骤(2)中的总有机碳含量计算模型以及步骤(3)中的曲线反演数据体计算泥页岩总有机碳含量数据体。
【IPC分类】G01V1/30, G06F19/00
【公开号】CN105277982
【申请号】CN201410360448
【发明人】郭旭升, 郭彤楼, 尹正武, 李文成, 陈祖庆, 彭嫦姿, 屈大鹏, 杨鸿飞, 陈超
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司勘探南方分公司
【公开日】2016年1月27日
【申请日】2014年7月25日