一种基于岩石物理意义的曲线校正方法技术

本发明专利技术公开了一种基于岩石物理意义的曲线校正方法，包括：根据实际测试资料，建立岩石物理模型；正演获得初始纵波速度与密度曲线交汇模板，配合井径曲线判断获得需要校正的纵波速度与密度曲线；将需要校正的曲线分为多个时窗段；设定一校正系数，利用加权平均Raymer岩石物理校正法和经验公式法，分别对多个时窗段内的曲线进行校正，获得变系数加权校正结果；针对变系数加权校正结果，采用井震相关系数为依据，并参考异常点的收敛情况，来判断校正结果是否合理；当井震相关系数低于设定的阈值，或异常点未收敛，修正校正系数重新进行变系数加权校正结果的计算；当井震相关系数等于或高于设定的阈值且异常点收敛，校正结果合理，输出校正后的曲线。

【技术实现步骤摘要】
一种基于岩石物理意义的曲线校正方法
本专利技术涉及石油地球物理勘探领域，尤指一种基于岩石物理意义的曲线校正方法。
技术介绍
在石油勘探的地震解释和储层预测过程中，需要高质量的声波时差、密度测井曲线，用来制作高精度的合成地震记录或者建立较为准确的低频模型。但很多情况下是实际的声波时差和密度测井曲线均存在严重失真，影响了井震标定和反演的准确性。尤其是浅探测测井系列(声波、密度等)易受到井筒环境的影响，扩径对测井曲线质量的影响较大，易使得测井数据失真。基于这些问题，特别需要一种合理的曲线校正方法，进行测井曲线质量控制并校正，使之能反映真实的地层信息。大体上，有两类测井曲线校正方法：一种是经验公式法，利用较少受井径垮塌影响的电阻率曲线，采用Faust公式，计算声波曲线；然后采用Gardner公式输入声波曲线计算密度曲线；另一种是岩石物理法，通过建立地区性的岩石物理模型，进而正演声波曲线和密度曲线，从而实现曲线校正。上述第一种方法，当井径垮塌非常严重的时候，深电阻率曲线也会失真，这样通过经验公式进行曲线校正也存在误差；而第二种方法，受控于岩石建模的输入参数较难取得，有时候精度也较低，误差也较大。综上所述，现有曲线校正方法各有优缺点，一种方法往往难以完整的进行曲线校正。
技术实现思路
为解决单一曲线校正方法难以满足需求，且校正误差较大的问题，本专利技术提供一种岩石物理意义下的综合曲线校正方法，以提高曲线校正精度和稳定度。本专利技术提出的基于岩石物理意义的曲线校正方法包括：步骤1，根据实际测试资料，建立岩石物理模型；步骤2，利用所述岩石物理模型正演获得初始纵波速度与密度曲线交汇模板，基于该初始纵波速度与密度曲线交汇模板配合井径曲线判断获得需要校正的纵波速度与密度曲线；步骤3，将所述需要校正的纵波速度与密度曲线分为多个时窗段；步骤4，设定一校正系数，利用加权平均Raymer岩石物理校正法和经验公式法，分别对所述多个时窗段内的需要校正的纵波速度与密度曲线进行校正，获得变系数加权校正结果；步骤5，对变系数加权校正结果做地震合成记录，获得井震相关系数，将井震相关系数与一设定的阈值进行比较，并参考所述初始纵波速度与密度曲线交汇模板上的异常点的收敛情况，判断校正结果是否合理；当井震相关系数低于所述设定的阈值，或异常点未收敛，判断校正结果不合理，执行步骤4，修正所述校正系数重新进行变系数加权校正结果的计算；当井震相关系数等于或高于所述设定的阈值，且异常点收敛，判断校正结果合理，输出校正后的曲线。进一步的，步骤1包括：利用改进的Raymer方程构建低-中孔隙度固结岩石物理模型，其中，纵波速度及密度的计算公式为：Vp＝(1-φ)kVma+φVfl；(1)ρ＝(1-φ)lρma+φρfl；(2)其中，Vp、Vma、Vfl分别为流体饱和岩石、岩石骨架和流体的纵波速度，Vma通过Voigt-Reuss-Hill模型确定的砂泥岩骨架体积模量、剪切模量和密度来计算获得，Vfl通过Wood公式计算的混合流体弹性分量体积模量和密度来计算获得；ρ、ρma和ρfl分别为流体饱和岩石、岩石骨架和流体的密度；φ为岩石的孔隙度；k、l取值范围分别为1.9-2.1和0.9-1.1。进一步的，步骤4包括：根据步骤1获取的岩石物理模型，正演整条曲线的纵波速度与密度，计算得到改进型Raymer方程岩石物理校正结果Vp_Raymer和ρ_Raymer；根据步骤1获取的岩石物理模型，基于Faust公式，由电阻率曲线重构声波时差曲线计算得到校正结果为Vp_Faust，基于Gardner公式，由声波曲线重构密度曲线计算得到校正结果ρ_Gardner；设定时窗Wi内的校正系数mWi、nWi(其中：mWi为Vp_Raymer和ρ_Raymer的校正系数；nWi为Vp_Faust和ρ_Gardner的校正系数)，mWi+nWi＝1；利用加权平均Raymer岩石物理校正法和经验公式法计算变系数加权校正结果，公式如下：[Vp]Wi＝mWiVp_Raymer+nWiVp_Faust；(3)[ρ]Wi＝mWiρ_Raymer+nWiρ_Gardner；(4)其中，[Vp]Wi、[ρ]Wi为时窗Wi内曲线的综合校正结果。本专利技术提出的基于岩石物理意义的曲线校正方法，在利用区域性Raymer声波曲线和密度曲线交会模板和井径曲线判断曲线失真的基础上，对曲线进行分时窗校正，并针对不同时窗，创新性地采用不同系数加权平均改进型Raymer岩石物理校正法和经验公式法结果，得到一个综合校正结果。克服了以往单一曲线校正方法校正误差大的问题，大大提高了曲线校正的精度和稳定度，是一种岩石物理意义控制下的综合曲线校正方法。附图说明此处所说明的附图用来提供对本专利技术的进一步理解，构成本申请的一部分，并不构成对本专利技术的限定。在附图中：图1为本专利技术一实施例的基于岩石物理意义的曲线校正方法流程图；图2为本专利技术一具体实施例的QT1井盒8段纵波速度与密度校正前交会图；图3为本专利技术一具体实施例的QT1井需要校正井段分时窗结果示意图；图4为本专利技术一具体实施例的QT1井盒8段纵波速度与密度校正后交会图；图5A及图5B分别为本专利技术一具体实施例的QT1井校正前、后曲线合成记录示意图。具体实施方式以下配合图示及本专利技术的较佳实施例，进一步阐述本专利技术为达成预定专利技术目的所采取的技术手段。图1为本专利技术一实施例的基于岩石物理意义的曲线校正方法流程图。如图1所示，该方法包括：步骤1，根据实际测试资料，建立岩石物理模型。其中，步骤1的具体步骤包括：从实际测试资料中获取温压参数、流体参数等资料，具体包括：岩矿分析资料、岩心分析资料及流体化验分析资料、地区温度压力资料、地层水类型、矿化度、气比重等。岩石骨架参数较难取得，本专利技术通过选取多口已知井井径好且测井质量好的井段计算目的层段岩石矿物模量值从而确定建模岩石骨架参数，主要包括粘土和石英的体积模量、剪切模量和密度。进一步的，利用改进的Raymer方程构建低-中孔隙度固结岩石物理模型，其中，纵波速度及密度的计算公式为：Vp＝(1-φ)kVma+φVfl；(1)ρ＝(1-φ)lρma+φρfl；(2)其中，Vp、Vma、Vfl分别为流体饱和岩石、岩石骨架和流体的纵波速度，Vma通过Voigt-Reuss-Hill模型确定的砂泥岩骨架体积模量、剪切模量和密度来计算获得，Vfl通过Wood公式计算的混合流体弹性分量体积模量和密度来计算获得；ρ、ρma和ρfl分别为流体饱和岩石、岩石骨架和流体的密度；φ为岩石的孔隙度；k、l取值范围分别为1.9-2.1和0.9-1.1，本专利技术将原始Raymer方程的常指数项改进为变指数k、l，可以控制校正结果的稳定性。步骤2，利用岩石物理模型正演获得初始纵波速度与密度曲线交汇模板，基于该初始纵波速度与密度曲线交汇模板配合井径曲线判断获得需要校正的纵波速度与密度曲线。图2为本专利技术一具体实施例的QT1井盒8段纵波速度与密度校正前交会图，如图2所示，结合井径曲线，可以判定纯砂岩线和纯泥岩线之外的圆圈区域为异常点，需要校正，并将异常点反投回测井曲线标记，可以得到需要校正的曲线段。步骤3，将需要校正的纵波速度与密度曲线分为多个时窗段，以便于下一步针对不同时窗段进行校正。图3为本专利技术一具体实施例本文档来自技高网...

【技术保护点】
一种基于岩石物理意义的曲线校正方法，其特征在于，该方法包括：步骤1，根据实际测试资料，建立岩石物理模型；步骤2，利用所述岩石物理模型正演获得初始纵波速度与密度曲线交汇模板，基于该初始纵波速度与密度曲线交汇模板配合井径曲线判断获得需要校正的纵波速度与密度曲线；步骤3，将所述需要校正的纵波速度与密度曲线分为多个时窗段；步骤4，设定一校正系数，利用加权平均Raymer岩石物理校正法和经验公式法，分别对所述多个时窗段内的需要校正的纵波速度与密度曲线进行校正，获得变系数加权校正结果；步骤5，对变系数加权校正结果做地震合成记录，获得井震相关系数，将井震相关系数与一设定的阈值进行比较，并参考所述初始纵波速度与密度曲线交汇模板上的异常点的收敛情况，判断校正结果是否合理；当井震相关系数低于所述设定的阈值，或异常点未收敛，判断校正结果不合理，执行步骤4，修正所述校正系数重新进行变系数加权校正结果的计算；当井震相关系数等于或高于所述设定的阈值，且异常点收敛，判断校正结果合理，输出校正后的曲线。

【技术特征摘要】
1.一种基于岩石物理意义的曲线校正方法，其特征在于，该方法包括：步骤1，根据实际测试资料，建立岩石物理模型；步骤2，利用所述岩石物理模型正演获得初始纵波速度与密度曲线交汇模板，基于该初始纵波速度与密度曲线交汇模板配合井径曲线判断获得需要校正的纵波速度与密度曲线；步骤3，将所述需要校正的纵波速度与密度曲线分为多个时窗段；步骤4，设定一校正系数，利用加权平均Raymer岩石物理校正法和经验公式法，分别对所述多个时窗段内的需要校正的纵波速度与密度曲线进行校正，获得变系数加权校正结果；步骤5，对变系数加权校正结果做地震合成记录，获得井震相关系数，将井震相关系数与一设定的阈值进行比较，并参考所述初始纵波速度与密度曲线交汇模板上的异常点的收敛情况，判断校正结果是否合理；当井震相关系数低于所述设定的阈值，或异常点未收敛，判断校正结果不合理，执行步骤4，修正所述校正系数重新进行变系数加权校正结果的计算；当井震相关系数等于或高于所述设定的阈值，且异常点收敛，判断校正结果合理，输出校正后的曲线；其中，步骤4，设定一校正系数，利用加权平均Raymer岩石物理校正法和经验公式法，分别对所述多个时窗段内的需要校正的纵波速度与密度曲线进行校正，获得变系数加权校正结果，包括：设定时窗Wi内的校正系数mWi、nWi，mWi+nWi＝1，其中：mWi为Vp_Raymer和ρ_Raymer的校正系数；nWi为Vp_Faust和ρ_Gardner的校正系数；根据步骤1获取的岩石物理模型，正演整条曲线的纵波速度与密度，计算得到改进型...

【专利技术属性】
技术研发人员：窦玉坛，张素磊，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11