一种大斜度井/水平井随钻电磁波测井资料钻后反演方法技术

本发明专利技术公开了一种大斜度井/水平井随钻电磁波测井资料钻后反演方法。所述反演方法包括如下步骤：s1.获取邻近导眼井测井资料，提取储层水平电阻率序列；s2.基于大斜度井/水平井随钻成像资料，确定地层界面位置、倾角信息；s3.对步骤s2划分后的地层进行逐层处理，并建立简化的一维三层反演模型以减少反演参数个数；s4.根据步骤s1提取的储层水平电阻率序列变化范围，对反演参数施加地质约束；s5.对反演参数进行多初值自适应初始化；s6.利用自适应乘子正则化Gauss‑Newton算法对地层模型进行迭代更新；s7.优选并输出反演结果。本发明专利技术方法提高了大斜度井/水平井随钻电磁波测井钻后反演的速度与精度。

A post drilling inversion method for electromagnetic logging data of large deviated wells / horizontal wells while drilling

The invention discloses a post drilling inversion method for electromagnetic wave logging while drilling in highly deviated wells and horizontal wells. The inversion method includes the following steps: s1. obtains the logging data from adjacent guide hole, extracts the horizontal resistivity sequence of the reservoir, and s2. determines the formation interface position and dip information based on the large inclined well / horizontal well with the drilling imaging data; s3. makes the layer by layer on the strata after step S2, and establishes a simplified one dimension three layer inversion. The model is used to reduce the number of inversion parameters; s4. applies geological constraints to the inversion parameters according to the variation range of the reservoir horizontal resistivity sequence extracted by step S1; s5. adaptively initializes the initial values of the inversion parameters; s6. uses adaptive multiplier regularization Gauss Newton algorithm to update the ground model; s7. is optimized and lost. The inversion results are obtained. The method improves the speed and accuracy of drilling electromagnetic wave logging while drilling.

【技术实现步骤摘要】
一种大斜度井/水平井随钻电磁波测井资料钻后反演方法
本专利技术涉及石油勘探开发领域，属于电测井方法范畴，具体的说是涉及一种大斜度井/水平井随钻电磁波测井资料钻后反演方法。
技术介绍
大斜度井/水平井钻井过程中，随钻电磁波测井视电阻率曲线“犄角”是地层界面的良好指示，但其远高于储层真实电阻率，这对储层精确测井评价提出了巨大的挑战。同时，受地层邻近地层界面、相对倾角及各向异性等影响，不同探测深度的视电阻率曲线分离严重，地层真实电阻率获取困难。因此，地层真实电阻率信息的准确反演是储层评价首要解决的关键问题之一。截至目前，随钻电磁波测井资料快速反演在钻后储层精细评价方面主要存在以下四个问题：1)反演计算中采用的多层模型速度慢；2)多层模型反演参数多，Jacobian计算量大；3)多层反演模型假定各个层倾角一致，其在倾角快速变化的地层结构中的反演结果精度降低；4)从最优化角度考虑，目标函数存在多个局部极小值，反演结果易困于局部极小值。
技术实现思路
本专利技术的目的在于提出一种大斜度井/水平井随钻电磁波测井资料钻后反演方法，以期为储层评价提供准确地层信息。本专利技术为了实现上述目的，采用如下技术方案：一种大斜度井/水平井随钻电磁波测井资料钻后反演方法，包括如下步骤：s1.从邻井资料中提取储层水平电阻率序列，作为地质约束条件；s2.根据随钻成像资料，获取N层地层界面的位置和倾角信息；s3.对步骤s2划分的N层地层进行逐层处理，建立以i层为目的层的一维三层反演模型，其中，第i层倾角信息已知，井眼与地层位置关系明确；令i的初始值为1；上述反演模型中，待反演参数包括上/下围岩电阻率以及当前地层的各向异性电阻率；s4.根据步骤s1提取的储层水平电阻率序列，对待反演参数施加地质约束；s5.反演模型初始化，对待反演参数进行多初值自适应初始化；s6.利用自适应乘子正则化Gauss-Newton算法对i层的反演模型进行迭代更新：利用自适应乘子正则化Gauss-Newton算法，对随钻电磁波测井资料进行反演，若实测资料与模拟资料拟合差在容许精度误差允许之内，执行步骤s7；若误差较大，执行步骤s5；s7.对步骤s6得到的多个反演结果择优选取，确定最优解并输出当前地层电阻率参数；s8.令i＝i+1，判断新的i取值是否小于或等于N，若是，则建立以该新的i层为目的层的一维三层反演模型，重复执行上述步骤s4-s7；若否，则执行步骤s9；s9.输出所有地层的反演结果。优选地，所述步骤s5具体为：s5.1.基于步骤s1和s2给出的地层信息，分别确定各待反演参数初值的选取数量；各待反演参数初值选取方式，具体步骤参照s5.2-s5.4：s5.2.根据层界面位置及分层信息，确定目的层厚度，根据目的层厚度选取多组水平电阻率n1表示初值总数；s5.3.根据目的层水平电阻率初值确定垂直电阻率Rv；s5.4.三层反演模型分为三种，即高阻模型Rh>Rup,Rh>Rdn、低阻模型Rh<Rup,Rh<Rdn和阶梯状模型Rup<Rh<Rdn,Rdn<Rh<Rup；基于上述模型假设，对上/下围岩电阻率选取多组初值n2和n3分别表示上/下围岩电阻率初值总个数；s5.5.对s5.2-s5.4给出的多组初值进行随机组合，得到n1*n2*n3组初始模型；对初始模型进行自适应乘子正则化Gauss-Newton反演迭代，获取多组反演结果。优选地，所述步骤s5.2中，水平电阻率Rh的选取规则为：1)目的层层厚>2m时，选取Rh初值1个，为视电阻率曲线中间部分平均值；2)1m<目的层层厚<2m时，选取初值3个，其中一个初值为电阻率曲线中间部分平均值，另外两个初值的选取保证一个初值大于平均值，另一个小于平均值；3)目的层层厚<1m时，选取初值5个。优选地，所述步骤s5.3中，垂直电阻率Rv＝水平电阻率Rh×3。优选地，所述步骤s5.4中，若有上一层反演结果Rh_i-1，则令Rup＝Rh_i-1，n2＝1。优选地，所述步骤s6具体为：基于最小二乘方法，采用代价函数C(x),第k步迭代中可表示为：式中，代价函数的第一项表示实测数据d与正演响应S(x)的L2范数；x为待反演的参数矢量，x＝Rup,Rdn,Rh,Rv}T，其中，Rup,Rdn表示上下围岩电阻率，Rh,Rv表示当前地层的各向异性电阻率；Wd为不同测量模式数据的权重系数矩阵，Rk(x)表示自适应正则化项：式中，δ为实验值；采用Gauss-Newton最优化方法对(1)求解，另其对x导数为零，第k次迭代得：式中，J为雅克比矩阵，T表示矩阵的转置，本专利技术具有如下优点：本专利技术通过建立一维三层反演模型，结合自适应乘子正则化Gauss-Newton算法，实现了地层各向异性电阻率逐层快速提取，解决了现有反演技术速度慢、精度低的问题，满足了大斜度井/水平井随钻电磁波测井资料快速、精确评价的需求。在逐层反演过程中，待反演参数仅包括当前层水平、垂直电阻率和上、下围岩电阻率四个参数，极大的减少了Jacobian计算量，提高了反演速度。同时对围岩电阻率和目的层电阻率进行多初值选取，保证了反演结果的全局最优性和准确性。数值实验表明，利用该算法反演视深为10m的水平井井段仅需5秒钟，可以满足钻后随钻电磁波测井资料实时处理的需求。此外，本专利技术可以明确各个层与仪器的相对倾角，对倾角快速变化的地层也比较适用。附图说明图1为本专利技术中一种大斜度井/水平井随钻电磁波测井资料钻后反演方法的流程图；图2为大斜度井/水平井测井资料图；图3为地层界面的位置和倾角信息示意图；图4为一维三层反演地层模型示意图；图5为三种电阻率模型示意图；图6为反演实例中地层模型示意图；图7为反演实例中随钻电磁波测井响应示意图；图8为反演实例中随钻电磁波反演结果示意图。具体实施方式下面结合附图以及具体实施方式对本专利技术作进一步详细说明：结合图1所示，一种大斜度井/水平井随钻电磁波测井资料钻后反演方法，包括如下步骤：s1.从邻井资料中提取储层水平电阻率序列，作为地质约束条件。s2.获取大斜度井/水平井测井资料，如图2所示，根据随钻伽马成像、或中子密度成像等资料，获取N层地层界面的位置和倾角信息，如图3所示。s3.对步骤s2划分的N层地层进行逐层处理。以第i层为例，建立以i层为目的层的一维三层反演模型，如图4所示，T表示发射线圈，R1、R2分别表示近、远接收线圈。其中，第i层倾角信息已知，井眼与地层位置关系明确；令i的初始值为1。对一维三层反演模型而言，待反演参数包括上下围岩电阻率Rup、Rdn、当前层的各向异性电阻率Rh、Rv四个参数，从而减少了反演参数个数、降低反演难度。s4.根据步骤s1提取的储层水平电阻率序列，对待反演参数施加地质约束。s5.反演模型初始化，对待反演参数进行多初值自适应初始化。该步骤包括以下两个部分：基于层厚确定当前层电阻率初值个数；分别基于中间层为高阻模型、低阻模型与阶梯状模型的假设给定围岩初值。下面对该步骤s5作进一步详细描述：s5.1.基于步骤s1和s2给出的地层信息，分别确定各待反演参数初值的选取数量；各待反本文档来自技高网...

【技术保护点】
1.一种大斜度井/水平井随钻电磁波测井资料钻后反演方法，其特征在于，包括如下步骤：s1.从邻井资料中提取储层水平电阻率序列，作为地质约束条件；s2.根据随钻成像资料，获取N层地层界面的位置和倾角信息；s3.对步骤s2划分的N层地层进行逐层处理，建立以i层为目的层的一维三层反演模型，其中，第i层倾角信息已知，井眼与地层位置关系明确；令i的初始值为1；上述反演模型中，待反演参数包括上/下围岩电阻率以及当前地层的各向异性电阻率；s4.根据步骤s1提取的储层水平电阻率序列，对待反演参数施加地质约束；s5.反演模型初始化，对待反演参数进行多初值自适应初始化；s6.利用自适应乘子正则化Gauss‑Newton算法对i层的反演模型进行迭代更新：利用自适应乘子正则化Gauss‑Newton算法，对随钻电磁波测井资料进行反演，若实测资料与模拟资料拟合差在容许精度误差允许之内，执行步骤s7；若误差较大，执行步骤s5；s7.对步骤s6得到的多个反演结果择优选取，确定最优解并输出当前地层电阻率参数；s8.令i＝i+1，判断新的i取值是否小于或等于N，若是，则建立以该新的i层为目的层的一维三层反演模型，重复执行上述步骤s4‑s7；若否，则执行步骤s9；s9.输出所有地层的反演结果。...

【技术特征摘要】
1.一种大斜度井/水平井随钻电磁波测井资料钻后反演方法，其特征在于，包括如下步骤：s1.从邻井资料中提取储层水平电阻率序列，作为地质约束条件；s2.根据随钻成像资料，获取N层地层界面的位置和倾角信息；s3.对步骤s2划分的N层地层进行逐层处理，建立以i层为目的层的一维三层反演模型，其中，第i层倾角信息已知，井眼与地层位置关系明确；令i的初始值为1；上述反演模型中，待反演参数包括上/下围岩电阻率以及当前地层的各向异性电阻率；s4.根据步骤s1提取的储层水平电阻率序列，对待反演参数施加地质约束；s5.反演模型初始化，对待反演参数进行多初值自适应初始化；s6.利用自适应乘子正则化Gauss-Newton算法对i层的反演模型进行迭代更新：利用自适应乘子正则化Gauss-Newton算法，对随钻电磁波测井资料进行反演，若实测资料与模拟资料拟合差在容许精度误差允许之内，执行步骤s7；若误差较大，执行步骤s5；s7.对步骤s6得到的多个反演结果择优选取，确定最优解并输出当前地层电阻率参数；s8.令i＝i+1，判断新的i取值是否小于或等于N，若是，则建立以该新的i层为目的层的一维三层反演模型，重复执行上述步骤s4-s7；若否，则执行步骤s9；s9.输出所有地层的反演结果。2.根据权利要求1所述的一种大斜度井/水平井随钻电磁波测井资料钻后反演方法，其特征在于，所述步骤s5具体为：s5.1.基于步骤s1和s2给出的地层信息，分别确定各待反演参数初值的选取数量；各待反演参数初值选取方式，具体步骤参照s5.2-s5.4：s5.2.根据层界面位置及分层信息，确定目的层厚度，根据目的层厚度选取多组水平电阻率n1表示初值总数；s5.3.根据目的层水平电阻率初值确定垂直电阻率Rv；s5.4.三层反演模型分为三种，即高阻模型Rh>Rup,Rh>Rdn、低阻模型Rh<Rup,Rh<Rdn和阶梯状模型Rup<Rh&l...

【专利技术属性】
技术研发人员：王磊，范宜仁，霍立壮，邓少贵，邢涛，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东,37