[发明专利]一种渐变地层的电阻率测井响应计算方法及装置

权利要求书

1.一种渐变地层的电阻率测井响应计算方法，其特征在于，包括：

生成一过渡带指数为2的电导率指数侵入函数和一过渡带指数为2的电阻率指数侵入函数；

使用所述过渡带指数为2的电导率指数侵入函数对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率进行替换，进一步通过解析表达式计算矩阵元素；所述电导率指数侵入函数用于描述渐变地层电导率的变化；

使用所述过渡带指数为2的电阻率指数侵入函数对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率倒数进行替换，进一步通过解析表达式计算矩阵元素；所述电阻率指数侵入函数为所述电导率指数侵入函数的倒数；所述电阻率指数侵入函数用于描述渐变地层电阻率的变化；

其中，所述过渡带指数为2的电导率指数侵入函数为其中，σ(ρ)为径向半径ρ处的过渡带电导率；σ_t为原状地层电导率；σ_xo为冲洗带电导率；L_c1为过渡带中心对应的径向半径；所述过渡带中心为过渡带电导率与原状地层电导率之和的二分之一；

所述过渡带指数为2的电阻率指数侵入函数为其中，R(ρ)为径向半径ρ处的过渡带电阻率；R_t为原状地层电阻率，R_xo为冲洗带电阻率，

矩阵A表示为矩阵B表示为矩阵G表示为矩阵H表示为其中，为第m层地层第n区间的A矩阵；为第m层地层第n区间的B矩阵；和分别为第m层和第m±1层地层电磁波的传递矩阵；ρ_n为第n区间的径向半径；为基函数向量；σ(ρ)_m为第m层地层的所述电导率指数侵入函数。

2.根据权利要求1所述的渐变地层的电阻率测井响应计算方法，其特征在于，所述对地层的电阻率测井响应中的数值模式匹配法的A、B、G和H矩阵元素中包含的电导率根据所述过渡带指数为2的电导率指数侵入函数进行替换，进一步通过解析表达式计算矩阵元素，包括：

对A矩阵：B矩阵：G矩阵：和H矩阵：中形式为的积分中的σ(ρ)根据所述电导率指数侵入函数进行替换，进一步通过解析表达式计算该积分；x表示径向半径ρ的x次幂；

其中，为第m层地层第n区间的A矩阵；为第m层地层第n区间的B矩阵；和分别为第m层和第m±1层地层电磁波的传递矩阵；ρ_n为第n区间的径向半径；为基函数向量；σ(ρ)_m为第m层地层的所述电导率指数侵入函数。

3.根据权利要求1所述的渐变地层的电阻率测井响应计算方法，其特征在于，所述对地层的电阻率测井响应中的数值模式匹配法的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层地层的所述电阻率指数侵入函数。