[发明专利]基于随钻电磁波测井资料的地层各向异性电阻率提取方法

说明书

本发明公开了一种基于随钻电磁波测井资料的地层各向异性电阻率提取方法，包括步骤：s1.对输入的随钻电磁波测井曲线进行环境校正；s2.获取区块地质结构信息；s3.对测井数据滑动开窗处理；s4.建立单界面反演地层模型；s5.判断是否有上一滑动窗口的反演结果，若有则执行步骤s6，否则执行步骤s8；s6.将上一滑动窗口反演结果作为当前滑动窗口的反演初始值；s7.采用非线性最优化算法对随钻电磁波测井资料反演，若实测资料与模拟资料拟合差在容许精度误差允许之内，执行步骤s9，否则执行步骤s8；s8.采用自适应多初值选取方法，并结合非线性最优化算法对随钻电磁波测井资料反演；s9.对多个反演结果进行优化选取。本发明可为实时地质导向和储层评价提供准确电阻率信息。

技术领域

本发明涉及石油勘探开发领域，属于电测井方法范畴，具体的说是涉及一种基于随钻电磁波测井资料的地层各向异性电阻率提取方法。

背景技术

随钻电磁波测井因其分辨率高、对各向异性电阻率敏感，在地质结构相对明确的油藏中得到广泛应用。但因受地层邻近地层界面、相对倾角及各向异性等影响，不同探测深度的视电阻率曲线分离严重，存在异常电阻率“犄角”，导致视电阻率值远大于地层真实值。因此，地层真实电阻率信息的准确、实时反演是地质导向与储层评价首要解决的关键问题之一。

目前，随钻电磁波测井反演主要集中在钻后储层精细评价方面，而关于随钻过程中实时电阻率曲线校正研究极少。从最优化的角度来看，随钻电磁波测井资料的实时处理主要存在以下三个问题：1)反演计算中采用的多层模型速度慢；2)多层模型反演参数多，Jacobian计算量大；3)存在多个局部极小值，且已知地层信息少，反演结果易困于局部极小值。

发明内容

本发明的目的在于提出一种基于随钻电磁波测井资料的地层各向异性电阻率提取方法，以期为实时地质导向和储层评价提供准确电阻率信息。

本发明为了实现上述目的，采用如下技术方案：

基于随钻电磁波测井资料的地层各向异性电阻率提取方法，包括如下步骤：

s1.对输入的随钻电磁波测井曲线进行井眼环境校正；

s2.获取区块的地质结构信息；

s3.对测井数据进行滑动开窗处理，将初始复杂地层模型转化为一系列开窗后的地层模型；

s4.对开窗后的地层模型进一步简化，将邻近层做等效处理；根据区块层厚，在随钻电磁波测井过程中将开窗后的地层模型简化为单界面地层模型，即单界面反演地层模型；

s5.判断是否有上一滑动窗口的反演结果，如果有，则执行步骤s6，否则执行步骤s8；

s6.将上一滑动窗口反演的结果作为当前滑动窗口的反演初始值；

s7.采用非线性最优化算法对随钻电磁波测井资料进行反演，如果实测资料与模拟资料拟合差在容许精度误差允许之内，则执行步骤s9，否则执行步骤s8；

s8.采用自适应多初值选取方法，并结合非线性最优化算法对随钻电磁波测井资料进行反演；

s9.对步骤s8中的多个反演结果进行优化选取，确定最优地层水平和垂直电阻率。

优选地，所述步骤s7和步骤s8中，非线性最优化算法具体为：

采用Taylor展开，对非线性反演问题进行线性处理，同时配合正则化的手段；

反演所用的代价函数C(m)为：

C(m)＝||W·(S(m)-d)||²+λ||m-m_ref||²+ξ||n||² (1)