一种渐变地层的电阻率测井响应计算方法及装置制造方法及图纸

本发明专利技术提供了一种渐变地层的电阻率测井响应计算方法及装置，涉及地层电阻率测井响应技术领域。方法包括：对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率利用过渡带指数为2的电导率指数侵入函数进行替换，进一步通过解析表达式计算矩阵元素；电导率指数侵入函数用于描述渐变地层电导率的变化；对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率倒数利用过渡带指数为2的电阻率指数侵入函数进行替换，进一步通过解析表达式计算矩阵元素；电阻率指数侵入函数为电导率指数侵入函数的倒数。从而解决对于径向渐变地层，现有的数值模式匹配法的测井响应结果不准确、速度慢的问题。

【技术实现步骤摘要】
一种渐变地层的电阻率测井响应计算方法及装置
本专利技术涉及地层电阻率测井响应
，尤其涉及一种渐变地层的电阻率测井响应计算方法及装置。
技术介绍
电阻率测井是探测复杂地层油气分布的基本方法，其基本原理是通过仪器发射探头向地层发射信号，接收探头接收发射信号经地层回流的信号，由于仪器沿井眼移动时经过不同地层，会导致测量探头接收的信号不同，进一步，根据仪器特性将测量信号换算到视电阻率，再进行反演和解释，即可确定复杂地层油气分布。然而，对测井信号进行反演，需要海量的正演计算，这就要求快速高精度的正演计算方法。目前，数值模式匹配法是一种较好的正演计算方法，其将高维数值问题分解为一个方向解析解和其余方向数值解的结合，既保证精度又提高效率，是解决层状介质问题的快速有效方法，是电测井快速正演计算和反演研究的有效工具。目前数值模式匹配法的地层模型均假设地层电阻率径向分段均匀，但由于地层的孔隙性和渗透性，在钻井过程中，泥浆在压力差的作用下使地层的物性参数径向分布不均匀，这会对电阻率测井的结果产生影响，导致实际测量的地层电阻率与真电阻率相差较大，从而对准确地判断地层含油气饱和度带来很大困难。因此，要精本文档来自技高网...

【技术保护点】
1.一种渐变地层的电阻率测井响应计算方法，其特征在于，包括：生成一过渡带指数为2的电导率指数侵入函数和一过渡带指数为2的电阻率指数侵入函数；对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率根据所述过渡带指数为2的电导率指数侵入函数进行替换，进一步通过解析表达式计算矩阵元素；所述电导率指数侵入函数用于描述渐变地层电导率的变化；对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率倒数根据所述过渡带指数为2的电阻率指数侵入函数进行替换，进一步通过解析表达式计算矩阵元素；所述电阻率指数侵入函数为所述电导率指数侵入函数的倒数；所述电阻率指数侵入函数...

【技术特征摘要】
1.一种渐变地层的电阻率测井响应计算方法，其特征在于，包括：生成一过渡带指数为2的电导率指数侵入函数和一过渡带指数为2的电阻率指数侵入函数；对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率根据所述过渡带指数为2的电导率指数侵入函数进行替换，进一步通过解析表达式计算矩阵元素；所述电导率指数侵入函数用于描述渐变地层电导率的变化；对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率倒数根据所述过渡带指数为2的电阻率指数侵入函数进行替换，进一步通过解析表达式计算矩阵元素；所述电阻率指数侵入函数为所述电导率指数侵入函数的倒数；所述电阻率指数侵入函数用于描述渐变地层电阻率的变化。2.根据权利要求1所述的渐变地层的电阻率测井响应计算方法，其特征在于，所述电导率指数侵入函数为其中，σ(ρ)为径向半径ρ处的过渡带电导率；σt为原状地层电导率；σxo为冲洗带电导率；Lc1为过渡带中心对应的径向半径；所述过渡带中心为过渡带电导率与原状地层电导率之和的二分之一。3.根据权利要求2所述的渐变地层的电阻率测井响应计算方法，其特征在于，所述对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率根据所述过渡带指数为2的电导率指数侵入函数进行替换，进一步通过解析表达式计算矩阵元素，包括：对A矩阵：B矩阵：G矩阵：和H矩阵：中形式为的积分中的σ(ρ)根据所述电导率指数侵入函数进行替换，进一步通过解析表达式计算该积分；x表示径向半径ρ的x次幂；其中，为第m层地层第n区间的A矩阵；为第m层地层第n区间的B矩阵；和分别为第m层和第m±1层地层电磁波的传递矩阵；ρn为第n区间的径向半径；为基函数向量；σ(ρ)m为第m层地层的所述电导率指数侵入函数。4.根据权利要求2所述的渐变地层的电阻率测井响应计算方法，其特征在于，所述电阻率指数侵入函数为其中，R(ρ)为径向半径ρ处的过渡带电阻率；Rt为原状地层电阻率，Rxo为冲洗带电阻率，5.根据权利要求4所述的渐变地层的电阻率测井响应计算方法，其特征在于，所述对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率倒数根据所述过渡带指数为2的电阻率指数侵入函数进行替换，进一步通过解析表达式计算矩阵元素，包括：对A矩阵：B矩阵：G矩阵：和H矩阵：中形式为的积分中的根据所述电阻率指数侵入函数进行替换，进一步通过解析表达式计算该积分；x表示径向半径ρ的x次幂；其中，为第m层地层第n区间的A矩阵；为第m层地层第n区间的B矩阵；和分别为第m层和第m±1层地层电磁波的传递矩阵；ρn为第n区间的径向半径；为基函数向量；R(ρ)m为第m层地层的所述电阻率指数侵入函数。6.一种渐变地层的电阻率测井响应计算装置，其特征在于，包括：侵入函数生成单元，用于生成一过渡带指数为2的电导率指数侵入函数和一过渡带指数为2的电阻率指数侵入函数；电导率替换单元，用于对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率根据所述过渡带指数为2的电导率指数侵入函数进行替换，进一步通过解析表达式计算矩阵元素；所述电导率指数侵入函数用于描述渐变地层电导率的变化；电导率倒数替换单元，用于对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率倒数根据所述过渡带指数为2的电阻率指数侵入函数进行替换，进一步通过解析表达式计算矩阵元素；所述电...

【专利技术属性】
技术研发人员：唐章宏，程文娟，许月晨，熊焱春，代明忠，唐宇，田园诗，刘刚，宋永杨，
申请(专利权)人：中国石油天然气股份有限公司，北京工业大学，中国石油集团测井有限公司，
类型：发明
国别省市：北京,11