一种基于数值模拟的地层倾角对电阻率影响的校正方法技术

本发明专利技术公开一种基于数值模拟的地层倾角对电阻率影响的校正方法，包括以下步骤：首先运用有限元方法模拟不同地层倾角条件下感应测井的地层电阻率；再利用数值模拟结果建立电阻率正演模型；最后利用拟合出的正演模型，配合粒子群算法，将地层视电阻率Ra与正演模型的计算电阻率Rt比较，得出满足特定适应度函数的全局最优解时的地层水平电阻率Rh。本发明专利技术的有益效果为：1、本方法所利用的基础数据完全来自于现场，资料的获取方便，容易；2、利用拟合出的模型对地层进行倾角校正，既有一定的理论依据，又有实际资料得出的客观规律，所以本发明专利技术得到的地层水平电阻率更能反映地层的储层性质。

A correction method for the influence of formation angle on resistivity based on numerical simulation

【技术实现步骤摘要】
一种基于数值模拟的地层倾角对电阻率影响的校正方法
本专利技术属于地球物理勘探领域，涉及一种基于数值模拟的地层倾角对电阻率影响的校正方法。
技术介绍
随着石油工业的快速发展，大斜度井和水平井已经成为高效率、高效益的油气田勘探开发技术。在实际的测井环境中，地层上、下界面通常不会是平行的，在测井解释过程中，这种非平行界面地层环境对电阻率测井响应的影响是一个值得关注的问题。地层电阻率是影响油气评价的重要参数，当地层倾角或井斜角比较大时受地层各向异性影响最为明显，测井曲线解释过程中出现的问题越多。研究倾斜环境下地层的电测井响应特征，分析地层电性各向异性对电测井响应的影响，所得结论可为倾斜地层测井解释校正提供一定的理论参考，以期更好地将电测井资料应用于油气藏的评价中。在确定储层特性中测井电阻率起着至关重要的作用，测井电阻率是地层真电阻率的反映，其准确程度直接影响储层特性的准确描述。测井电阻率会受到井眼、围岩和层厚、泥浆、地层倾角(井斜角)以及地层各向异性的影响，因此测井所得的地层电阻率与真实地层的电阻率存在一定的差异，但它们之间同时也存在着某种联系。分析、归纳、总结地层倾角对测井电阻率的影响，结合前人的研究成果，能够探寻出地层倾角对测井电阻率响应的规律。测井曲线形态取决于所测量的特定的井眼和地层，近年来国内外学者采用了不同的方法来研究地层倾角对电阻率测井响应的影响，已取得了长足的发展和显著的成效。通过分析地层倾角对普通电阻率测井、双侧向电阻率测井、感应电阻率测井以及随钻电阻率测井响应特征的影响经验，已总结出了不同学者在研究中采用的方法以及取得的成果。
技术实现思路
本专利技术主要是克服现有技术中的不足之处，提出一种更能反映地层的储层性质的基于数值模拟的地层倾角对电阻率影响的校正方法。本专利技术解决上述技术问题所提供的技术方案是：一种基于数值模拟的地层倾角对电阻率影响的校正方法，包括以下步骤：S1、首先运用有限元方法模拟不同地层倾角条件下感应测井的地层电阻率；S2、再利用数值模拟结果建立电阻率正演模型：式中，Rt为正演模型的计算值，Rh为地层水平电阻率，Rv为地层垂直电阻率，θ为地层倾角，其中Rh就是需要通过倾角校正的地层电阻率，通过采用粒子群算法，将地层视电阻率Ra与正演模型的计算电阻率Rt比较，得出满足特定适应度函数的全局最优解时的地层水平电阻率Rh；S3、最后利用拟合出的正演模型，配合粒子群算法，将地层视电阻率Ra与正演模型的计算电阻率Rt比较，得出满足特定适应度函数的全局最优解时的地层水平电阻率Rh。进一步的是，所述步骤S1的具体过程为：S10、建立地层物理模型，其中地层的物理模型包括地层结构、模型参数，所述地层结构包括目的层的厚度h、地层的倾斜角度θ、井径D、以及线圈系的位置坐标；所述模型参数包括目的层电阻率为Rt、地层相对介电常数εr、地层相对磁导率μr、上下围岩的电阻率为Rs、井眼内泥浆的电阻率为Rf；S11、建立微分方程及边界条件：磁场的矢量波动方程：地层界面上的边界条件(表示在地层分界面上的磁场切向分量是连续的)：n×(H1-H2)＝0狄里克雷边界条件(表示在区域边界上棱边场为零)：H|∞＝0式中：H为磁场强度(A/m)，ω角频率，μ磁导率(H/m)，εc复介电常数，Ms为磁流密度矢量；S12、将地层物理模型体积V离散成有限个四面体单元；S13、选取插值函数；S14、建立有限元方程；S15、对于复合线圈系，由迭加原理可得到感应测井的视电阻率值；S16、最后对模拟结果进行分析本专利技术的有益效果为：1、本方法所利用的基础数据完全来自于现场，资料的获取方便，容易；2、利用拟合出的模型对地层进行倾角校正，既有一定的理论依据，又有实际资料得出的客观规律，所以本专利技术得到的地层水平电阻率更能反映地层的储层性质；3、该方法操作简单，实用性较强，且有较好的推广应用价值；4、本专利技术提供的一种结合地层电阻率的正反演模型，对倾斜地层进行电阻率倾角校正的方法，能够实现对倾斜地层的连续校正。附图说明图1为四种不同倾角地层物理及网格划分模型；图2为数据一不同倾角地层的深感应视电阻率变化曲线；图3为数据一不同倾角地层的中感应视电阻率变化曲线；图4为数据二不同倾角地层的深感应视电阻率变化曲线；图5为数据二不同倾角地层的中感应视电阻率变化曲线；图6为不同角度时ln(Rv/Rh)和ln(Ra/Rh)交会图；图7为不同各向异性系数时cos(θ)和ln(Ra/Rh)交会图；图8为粒子群算法流程图；图9为研究区A井倾角校正结果。具体实施方式下面结合实施例和附图对本专利技术做更进一步的说明。如图1所示，本专利技术的一种基于数值模拟的地层倾角对电阻率影响的校正方法，包括以下步骤：S10、采用有限元方法模拟不同倾角条件下的地层电阻率：。①模型建立(1)建立地层物理模型：地层的物理模型应包括如下内容：地层结构，具体表现为目的层的厚度h，地层的倾斜角度θ，井径D，以及线圈系的位置坐标等(如图1所示)；模型参数包括目的层电阻率为Rt，地层相对介电常数εr，地层相对磁导率μr，上下围岩的电阻率为Rs，井眼内泥浆的电阻率为Rf。(2)建立微分方程及边界条件：三维有限元的核心是求解电场或者磁场波动方程来确定场量在感兴趣的区域分布。在感应测井中，线圈尺寸以及波长同线圈间距相比是非常小的，因而可以将发射线圈等效为振荡的磁偶极子源。磁场的矢量波动方程：地层界面上的边界条件(表示在地层分界面上的磁场切向分量是连续的)：n×(H1-H2)＝0(2)狄里克雷边界条件(表示在区域边界上棱边场为零)：H|∞＝0(3)式中：H为磁场强度(A/m)，ω角频率，μ磁导率(H/m)，εc复介电常数，Ms为磁流密度矢量；②模型求解(1)离散化：首先取一组整数给所有的四面体单元编号，用另外一组整数给四面体顶点处的所有结点编号，用4×M的整型数组n(i，e)将单元编码和结点编码联系起来。其中i＝0，1，2，3表示单元内的结点编码，e＝1，2，3…M表示体单元编码，M表示体单元总数，n(i，e)表示第e个单元的第i个顶点在全局顶点内的编码。(2)插值函数选取：采用四面体单元插值函数，插值函数可以表示为：对于上面的基于结点插值函数的四面体单元，在矢量场中，尤其是在电磁场中应用时，不方便在界面和外边界强加边界条件，且可能造成伪解，因此我们将在它的基础上，使用基于棱边元的矢量有限元。在四面体单元中，我们导出了其线性插值函数，即为矢量函数为：假设e1表示从结点1指向结点2的单位矢量，因而是结点1处线性函数，是结点2处的线性函数，最终可得到棱边i的矢量基函数为：(3)建立有限元方程：矢量波动方程式(1)所对应的泛函为：利用矢量格林公式将其转变为：接下来对方程进行矢量元离散，e表示第e个单元，表示第e个单元的第i个矢量基函数，离散后，得到：其中，M表示总的单元个数，Fe表达式如下：式中各个参数都是在单元e内的参数，为磁偶极子源在单位e内的值，若磁偶极子不在单元e内，则将式代入式(8)，得其中，{He}表示单元e中每一条棱边权值组成的矢量，包含用Fe对每一条棱边分量求偏导，其中驻点条件即求泛函的极值，可得[A]{H}-ω2μ[B]{H}＝ω2μ{C}(15)式(15)即是我们所要求的有限元方程组，[A]，[B]，{本文档来自技高网...

【技术保护点】
一种基于数值模拟的地层倾角对电阻率影响的校正方法，其特征在于，包括以下步骤：S1、首先运用有限元方法模拟不同地层倾角条件下感应测井的地层电阻率；S2、再利用数值模拟结果建立电阻率正演模型：

【技术特征摘要】
1.一种基于数值模拟的地层倾角对电阻率影响的校正方法，其特征在于，包括以下步骤：S1、首先运用有限元方法模拟不同地层倾角条件下感应测井的地层电阻率；S2、再利用数值模拟结果建立电阻率正演模型：式中，Rt为正演模型的计算值，Rh为地层水平电阻率，Rv为地层垂直电阻率，θ为地层倾角，其中Rh就是需要通过倾角校正的地层电阻率，通过采用粒子群算法，将地层视电阻率Ra与正演模型的计算电阻率Rt比较，得出满足特定适应度函数的全局最优解时的地层水平电阻率Rh；S3、最后利用拟合出的正演模型，配合粒子群算法，将地层视电阻率Ra与正演模型的计算电阻率Rt比较，得出满足特定适应度函数的全局最优解时的地层水平电阻率Rh。2.根据权利要求1所述的一种基于数值模拟的地层倾角对电阻率影响的校正方法，其特...

【专利技术属性】
技术研发人员：赵军，蒲万丽，刘凯，武延亮，肖承文，向薪燃，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：四川,51