CN103883318A - 一种中低矿化度地层水储层的参数反演方法及装置 - Google Patents

Abstract

本发明提供了一种中低矿化度地层水储层的参数反演方法及装置，方法包括：获取双侧向测井资料、阵列感应测井资料及辅助测井曲线资料；根据双侧向测井资料与辅助测井曲线资料对储层段进行分层；根据双侧向测井曲线、阵列感应曲线生成对应层的初始的原状地层电阻率、侵入带电阻率、侵入半径值；生成对应层的双侧向模拟响应曲线以及阵列感应模拟响应曲线；分别判断对应层的双侧向模拟响应曲线与对应层的双侧向实测曲线及阵列感应模拟响应曲线与对应层的阵列感应原始响应曲线是否一致，一致，将原状地层电阻率、侵入带电阻率和侵入半径值作为储层段的反演结果输出；否则修改初始的原状地层电阻率、侵入带电阻率和侵入半径值。

Description

一种中低矿化度地层水储层的参数反演方法及装置

技术领域

本发明涉及一种石油测井领域，尤其涉及一种利用测井资料进行地层参数反演的技术，具体的讲是一种中低矿化度地层水储层的参数反演方法及装置。

背景技术

在储层段，由于受到泥浆侵入、井眼、围岩等的影响，使得测得的双侧向、阵列感应电阻率偏离真实电阻率，从而使得通过电阻率曲线进行油水层识别的能力降低，计算的含油饱和度误差较大。因此，需要利用反演方法还原真实地层面目。

传统的反演方法一直以来难以实用化，其原因在于它们都是纯数学方法，只要满足模拟出的测井曲线与实测曲线误差小于给定的ε即可，而不管在岩石物理、地质环境及工程条件等方面是否合理，因此，利用传统反演方法得到的结果不够合理。

此外，对于中低矿化度地层水储层，当储层为油气层时，其电阻率一般较高。由于感应测井对水层类低阻储层反应明显，对油气类高阻储层反应不太明显，这时单纯利用阵列感应测井资料进行储层参数的反演效果有限。

发明内容

为最大限度地恢复地层真实面目，改善油水层的识别能力和含油饱和度评价的精度。

本发明实施例提供了一种中低矿化度地层水储层的参数反演方法，方法包括：

步骤1，获取双侧向测井资料、阵列感应测井资料及辅助测井曲线资料，所述双侧向测井资料包括双侧向测井曲线，所述阵列感应测井资料包括从不同阵列单元测量的阵列感应原始响应曲线和由所述阵列感应原始响应曲线进行合成聚焦处理后的阵列感应测井曲线；

步骤2，根据所述双侧向测井资料与辅助测井曲线资料对储层段进行分层；

步骤3，根据所述的双侧向测井曲线生成对应层的初始的原状地层电阻率、侵入带电阻率；

步骤4，根据所述的阵列感应曲线生成对应层的侵入半径值；

步骤5，根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过对模拟算法进行加速生成对应层的双侧向模拟响应曲线；

步骤6，根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过数值模式匹配法生成对应层的阵列感应模拟响应曲线；

步骤7，分别判断所述对应层的双侧向模拟响应曲线与对应层的双侧向实测曲线以及所述阵列感应模拟响应曲线与对应层的阵列感应原始响应曲线是否一致，如果一致，转到步骤8，如果不一致，转到步骤9；

步骤8，将所述原状地层电阻率、侵入带电阻率和侵入半径值作为储层段的反演结果输出；

步骤9，根据所述双侧向实测曲线和阵列感应原始响应曲线，结合对应层的地质特征、油藏特征，修改所述初始的原状地层电阻率、侵入带电阻率和侵入半径值，重复步骤5-9，输出最终的原状地层电阻率、侵入带电阻率和侵入半径值。

此外，本发明还提供了一种中低矿化度地层水储层的参数反演装置，装置包括：

资料获取装置，用于获取双侧向测井资料、阵列感应测井资料及辅助测井曲线资料，所述双侧向测井资料包括双侧向测井曲线，所述阵列感应测井资料包括从不同阵列单元测量的阵列感应原始响应曲线和由所述阵列感应原始响应曲线进行合成聚焦处理后的阵列感应测井曲线；

分层装置，用于根据所述双侧向测井资料与辅助测井曲线资料对储层段进行分层；

电阻率确定装置，用于根据所述的双侧向测井曲线生成对应层的初始的原状地层电阻率、侵入带电阻率；

半径值确定装置，用于根据所述的阵列感应曲线生成对应层的侵入半径值；

双侧向模拟响应曲线生成装置，用于根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过对模拟算法进行加速生成对应层的双侧向模拟响应曲线；

阵列感应模拟响应曲线生成装置，用于根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过数值模式匹配法生成对应层的阵列感应模拟响应曲线；

判断装置，用于分别判断所述对应层的双侧向模拟响应曲线与对应层的双侧向实测曲线以及所述阵列感应模拟响应曲线与对应层的阵列感应原始响应曲线是否一致；

参数修改装置，用于根据所述双侧向实测曲线和阵列感应原始响应曲线，结合对应层的地质特征、油藏特征，修改所述初始的原状地层电阻率、侵入带电阻率和侵入半径值；

结果输出装置，用于将所述判断装置判断结果为一致时的原状地层电阻率、侵入带电阻率和侵入半径值作为储层段的反演结果输出。

侧向测井对油气类高阻储层反应明显，因此，联合应用双侧向测井资料、阵列感应测井资料进行中低矿化度地层水储层参数的反演会提高参数反演的精度，从而可大大改善对储层流体性质的识别能力。

利用双侧向测井资料、阵列感应测井资料进行中低矿化度地层水储层参数反演的方法及装置充分考虑了地质背景、油藏条件、工程环境等因素，主动修改地层模型，以模拟出的测井曲线与实测曲线是否吻合作约束，从而使得反演结果更为合理。因此，利用本发明的参数反演方法及装置可以最大限度地恢复地层真实面目，大大提高油水层的识别能力和含油饱和度评价的精度。

为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明公开的一种中低矿化度地层水储层的参数反演方法流程图；

图2为本发明公开的一种中低矿化度地层水储层的参数反演装置的框图；

图3为利用本发明的实施例进行地层参数反演生成的双侧向测井、阵列感应测井交互式反演成果图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

如图1所示，本发明公开了一种中低矿化度地层水储层的参数反演方法，包括：

步骤S101，获取双侧向测井资料、阵列感应测井资料及辅助测井曲线资料，所述双侧向测井资料包括双侧向测井曲线，所述阵列感应测井资料包括从不同阵列单元测量的阵列感应原始响应曲线和由所述阵列感应原始响应曲线进行合成聚焦处理后的阵列感应测井曲线；

步骤S102，根据所述双侧向测井资料与辅助测井曲线资料对储层段进行分层；

步骤S103，根据所述的双侧向测井曲线生成对应层的初始的原状地层电阻率、侵入带电阻率；

步骤S104，根据所述的阵列感应曲线生成对应层的侵入半径值；

步骤S105，根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过对模拟算法进行加速生成对应层的双侧向模拟响应曲线；

步骤S106，根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过数值模式匹配法生成对应层的阵列感应模拟响应曲线；

步骤S107，分别判断所述对应层的双侧向模拟响应曲线与对应层的双侧向实测曲线以及所述阵列感应模拟响应曲线与对应层的阵列感应原始响应曲线是否一致，如果一致，转到步骤8，如果不一致，转到步骤9；

步骤S108，将所述原状地层电阻率、侵入带电阻率和侵入半径值作为储层段的反演结果输出；

步骤S109，根据所述双侧向实测曲线和阵列感应原始响应曲线，结合对应层的地质特征、油藏特征，修改所述初始的原状地层电阻率、侵入带电阻率和侵入半径值，重复步骤S105-S109，输出最终的原状地层电阻率、侵入带电阻率和侵入半径值。

优选的，本发明的辅助测井曲线资料包括：自然伽马测井曲线、井径测井曲线、密度测井曲线以及自然电位测井曲线。

优选的，本发明的实施例中根据所述双侧向测井资料与辅助测井曲线资料对储层段进行分层，包括：

根据所述双侧向测井曲线特征依据电阻率曲线的斜率变化趋势确定层界面；依据碎屑岩地层沉积规律，以3层介质为例，其电阻率曲线的变化规律可分为六种，分别为大小大、小大小、大大大、小小小、大大小、小小大，对于像大小大和小大小有高低起伏变化的地层分层可通过对电阻率求一阶导数来解决，如大小大、小大小地层，其电阻率对深度的一阶导数的符号分别是-、+和+、-，对一阶导数的符号有变化的点称之为拐点，记录下储层所有拐点的位置，每每取两拐点中点位置即地层层界面位置。对于大大大、小小小这类单向变化的地层由于其一阶导数符号没有变化则不能通过此类方法解决，本发明是通过设定电阻率变化幅度，如电阻率每增加10%或每减少10%即可分为一个层，然后结合辅助测井曲线资料再进行细分。而对于大大小、小小大这类地层则是上述两种方法的结合。

根据所述辅助测井曲线资料和确定层界面，分别将所述双侧向测井曲线特征不明显的泥岩夹层、致密层、垮塌层分出来。

优选的，本发明的实施例中根据双侧向测井曲线生成对应层的初始的原状地层电阻率、侵入带电阻率包括：

对分层后的每一层，将双侧向测井曲线中的浅侧向在层中的均值作为对应层初始的侵入带电阻率；

将深侧向在层中的均值作为对应层初始的原状地层电阻率。

优选的，本发明的实施例中根据阵列感应曲线生成对应层的侵入半径值包括：

通过几何因子法对分层后的每一层计算出初步的侵入半径值：

$C_i=G_i^m{C}_m+G_i^{\mathrm{xo}}{C}_{\mathrm{xo}}^{\′}+(1-G_i^m-G_i^{\mathrm{xo}})C_t^{\′}i=1~N;$

式中C_i－阵列感应测井0.6m分辨率系列电导率值（电阻率倒数），S/m；

－井眼泥浆对阵列感应测井响应的贡献；

Cm－泥浆电导率值，S/m；

－侵入带对阵列感应测井响应的贡献，与侵入半径值有关；

C'_xo－由双侧向测井给出的侵入带初始电导率值，S/m；

C_t'－由双侧向测井给出的原状地层初始电导率值，S/m；

N－阵列感应测井0.6m分辨率系列曲线条数，且N=5或6；

通过阵列感应测井测得的每一层0.6m分辨率系列的N个值由上式可列出N个方程，通过最优化方法计算出每一层初步的侵入半径值；

其中，C_m根据所述阵列感应测井曲线计算得到，

根据所述井径测井曲线计算得到。

优选的，本发明的实施例中根据对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过对模拟算法进行加速生成对应层的双侧向模拟响应曲线包括：

根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，结合双侧向测井仪器的仪器模型参数，通过有限元素法生成对应层的双侧向模拟响应曲线；其中，所述双侧向测井仪器的仪器模型参数包括：仪器长度、仪器半径、双侧向各电极相对距离、各电极长度；其中，所述的对模拟算法进行加速包括：

缩短仪器模型中屏蔽电极尺寸、降低有限元素法中计算网格长度、对计算中各电极分场进行等效处理。

优选的，本发明的实施例中根据对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过数值模式匹配法生成对应层的阵列感应模拟响应曲线包括：

根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，结合阵列感应测井仪器的仪器模型参数，通过数值模式匹配法生成对应层的阵列感应模拟响应曲线；其中，

所述阵列感应测井仪器的仪器模型参数包括：仪器长度、仪器半径、各阵列单元线圈距、线圈匝数、工作频率。

优选的，本发明的实施例中判断所述对应层的双侧向模拟响应曲线与对应层的所述双侧向测井曲线以及阵列感应模拟响应曲线与对应层的所述阵列感应原始响应曲线是否一致，包括：

设定对应层的双侧向模拟响应曲线与对应层的所述双侧向测井曲线以及阵列感应模拟响应曲线与对应层的所述阵列感应原始响应曲线间的相关系数均为S：

$S=\frac{\mathrm{\ΣXY}-\frac{\mathrm{\ΣX\ΣY}}{N}}{\sqrt{(\mathrm{\Σ}{X}^2-\frac{{\left(\mathrm{\ΣX}\right)}^2}{N})(\mathrm{\Σ}{Y}^2-\frac{{\left(\mathrm{\ΣY}\right)}^2}{N})}},$ 式中：

N－目的层段内的测量点数；

X－对应层的目的层段内的数据集，其中，对高阻层，X为所述对应层的双侧向测井曲线目的层段内的数据集；对中阻层，X为所述对应层的阵列感应原始响应曲线、双侧向测井曲线目的层段内的数据集；对低阻层，X为所述对应层的阵列感应原始响应曲线目的层段内的数据集；

Y－与X相对应的所述层阵列感应模拟响应曲线和/或双侧向模拟曲线目的层段内的数据集；

当所述相关系数S的计算结果大于0.9，即认为对应层的阵列感应模拟响应曲线与对应层的所述阵列感应原始响应曲线、对应层的双侧向模拟响应曲线与对应层的所述双侧向实测曲线一致。

优选的，本发明的实施例中根据双侧向测井曲线和阵列感应原始响应曲线，结合对应层的地质特征、油藏特征，修改所述初始的原状地层电阻率、侵入带电阻率和侵入半径值，包括：

根据所述对应层所在井的邻井资料以及本井的录井资料、气测资料、取心资料判断所述对应层的流体性质以及所述流体性质下的地层电阻率范围；

根据所述对应层的泥浆性质以及地层属性判断侵入带电阻率范围；

根据地层浸泡时间长短及地层物性好坏判断侵入半径值范围；

根据所述对应层的双侧向模拟响应曲线与对应的双侧向实测曲线间以及阵列感应模拟响应曲线与对应的阵列感应原始响应曲线间的相关程度，并结合所述地层电阻率范围、侵入带电阻率范围以及侵入半径值范围，修改所述初始的原状地层电阻率、侵入带电阻率和侵入半径值。

此外，如图2所示，本发明还公开了一种中低矿化度地层水储层的参数反演装置，装置包括：

资料获取装置201，用于获取双侧向测井资料、阵列感应测井资料及辅助测井曲线资料，所述双侧向测井资料包括双侧向测井曲线，所述阵列感应测井资料包括从不同阵列单元测量的阵列感应原始响应曲线和由所述阵列感应原始响应曲线进行合成聚焦处理后的阵列感应测井曲线；

分层装置202，用于根据所述双侧向测井资料与辅助测井曲线资料对储层段进行分层；

电阻率确定装置203，用于根据所述的双侧向测井曲线生成对应层的初始的原状地层电阻率、侵入带电阻率；

半径值确定装置204，用于根据所述的阵列感应曲线生成对应层的侵入半径值；

双侧向模拟响应曲线生成装置205，用于根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过对模拟算法进行加速生成对应层的双侧向模拟响应曲线；

阵列感应模拟响应曲线生成装置206，用于根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过数值模式匹配法生成对应层的阵列感应模拟响应曲线；

判断装置207，用于分别判断所述对应层的双侧向模拟响应曲线与对应层的双侧向实测曲线以及所述阵列感应模拟响应曲线与对应层的阵列感应原始响应曲线是否一致；

参数修改装置208，用于根据所述双侧向实测曲线和阵列感应原始响应曲线，结合对应层的地质特征、油藏特征，修改所述初始的原状地层电阻率、侵入带电阻率和侵入半径值；

结果输出装置209，用于将所述判断装置判断结果为一致时的原状地层电阻率、侵入带电阻率和侵入半径值作为储层段的反演结果输出。

优选的，本发明的实施例中分层装置包括：

确定层界面模块，用于根据所述双侧向测井曲线特征依据斜率变化趋势确定层界面；

分层模块，用于根据所述辅助测井曲线资料和确定层界面，分别将所述双侧向测井曲线特征不明显的泥岩夹层、致密层、垮塌层分出来。

优选的，本发明的实施例中电阻率确定装置包括：

侵入带电阻率确定模块，用于对分层后的每一层，将双侧向测井曲线中的浅侧向在层中的均值作为对应层初始的侵入带电阻率；

原状地层电阻率确定模块，用于将深侧向在层中的均值作为对应层初始的原状地层电阻率。

本实施例的利用阵列感应测井资料、双侧向测井资料进行中低矿化度地层水储层参数反演的装置，充分考虑了地质背景、油藏条件和工程环境等因素，不断主动修改原状地层电阻率、侵入带电阻率和侵入半径值，使得模拟出的阵列感应原始响应与实测原始响应及双侧向模拟响应曲线与对应层的所述双侧向实测曲线是否吻合作为判断条件，如果吻合，则修改后的原状地层电阻率、侵入带电阻率和侵入半径值即为所求，这样使得反演结果更为合理。

图3为利用本发明的利用双侧向测井资料、阵列感应测井资料进行中低矿化度地层水储层参数反演的方法和装置的实施例，进行地层参数反演生成的双侧向-阵列感应测井交互式反演成果图。

在此具体实施例中，选用普遍应用的高分辨率阵列感应测井仪（HDIL）的测井资料作处理对象。在图3中：

第一道是井径曲线（CAL）、自然电位曲线（SP）以及自然伽马曲线（GR）。其中，曲线①表示井径曲线，曲线②表示自然电位曲线，曲线③表示自然伽马曲线。

第二道是HDIL合成聚焦后的0.6m分辨率阵列曲线，M2R1、M2R2、M2R3、M2R6、M2R9、M2RX分别表示探测深度为0.254m、0.508m、0.762m、1.522m、2.286m以及3.048m的测井曲线。其中，曲线④表示探测深度为0.254m的测井曲线，曲线⑤表示探测深度为0.508m的测井曲线，曲线⑥表示探测深度为0.762m的测井曲线，曲线⑦表示探测深度为1.522m的测井曲线，曲线⑧表示探测深度为2.286m的测井曲线，曲线⑨表示探测深度为3.048m的测井曲线。

第三道是双侧向测井道，其中，RLLS、RLLD分别表示浅侧向、深侧向实测曲线，2RLLS、2RLLD分别表示模拟的浅侧向、深侧向响应曲线。

第四道是地层模型道，表示经分层后各层的原状地层电阻率（2RT）、侵入带电阻率（2RXO）和侵入半径（2RI）值，通过该道可修改地层模型值。其中，曲线D表示经分层后各层的原状地层电阻率（2RT），曲线E表示经分层后各层的侵入带电阻率（2RXO），曲线F表示经分层后各层的侵入半径（2RI）。

第五道是HDIL经转换后的原始测量曲线（单位：ms/m），这里只选取了一种工作频率下6个阵列单元（图中从下至上，线圈距由短至长排列）的测量结果。其中，在这6组曲线中，每一组中的其中一条表示模拟出的阵列感应模拟响应曲线，另一条表示测量出的阵列感应原始响应曲线。

第六道是深度道，分别有井段的深度值以及各储层的解释结论，其中1号层解释结论为气层、2号层解释结论为干层、3号层解释结论为水层。

第七道是孔隙度曲线：其中，曲线A表示声波时差（AC）曲线，曲线B表示密度（DEN）曲线，曲线C表示补偿中子（CNL）曲线。当人机交互修改第四道的地层模型值，使得第三道的双侧向模拟响应曲线与实测的双侧向测井曲线以及第五道的阵列感应模拟响应曲线与实测的阵列感应原始响应曲线相一致时，表示此时的地层模型已接近真实的地层值，当满足相关系数S时即可输出所求的各层的原状地层电阻率、侵入带电阻率和侵入半径值。

本发明的此种方法可以大大提高油气水层的识别能力和含油饱和度的评价精度。

在图3中，从3445～3450m井段，试油证实为气层。由于测井时泥浆的深侵使得测量的地层电阻率下降很多，阵列感应测井探测曲线电阻率约为6Ω.m，通过常规反演，反演后略有改善，反演出的原状地层电阻率约为7.0Ω.m，由archie公式：

Rt/Rw=abΦ^-m(1-So)^-n

I=Rt/R0；

式中Rt为原状地层电阻率，Ω.m；

Rw为地层水电阻率，Ω.m；

Φ为地层孔隙度；

So为含油气饱和度；

a,b为系数，与岩性有关；

m,n为指数，分别与孔隙结构和油气水分布有关；

I为电阻率增大系数，它可作为油气水识别指标，其值越大，含油气的概率越大；

R0为同层组的纯水地层电阻率。

取m=2,n=2，a=1,b=1，Φ=15%，Rw=0.05Ω.m，则可计算出：

反演前：Rt=6.0Ω.m，I=1.5，So=40%；

常规方法反演后得到：Rt=7.0Ω.m，I=1.75，So=44%。

通过对比可以看出，常规方法反演后，油气水层识别能力略有提高，含油饱和度计算精度略有提高。

但是使用本方法后，由图3可以看出反演后Rt=14.0Ω.m，使用同样参数可计算出I＝3.5，So=61%，相比常规方法，大幅度提高了油气水层识别能力和含油饱和度计算精度。

本发明实施例的利用双侧向测井资料、阵列感应测井资料进行中低矿化度地层水储层参数反演的方法及装置，充分考虑了地质背景、油藏条件、工程环境等因素，主动修改地层模型，以模拟出的测井曲线与实测曲线是否吻合作约束，从而使得反演结果更为合理。因此，利用本发明的参数反演方法及装置可以最大限度的恢复地层真实面目，大大提高油水层的识别能力和含油饱和度评价的精度。

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (18)

1.一种中低矿化度地层水储层的参数反演方法，其特征在于，所述方法包括：

步骤1，获取双侧向测井资料、阵列感应测井资料及辅助测井曲线资料，所述双侧向测井资料包括双侧向测井曲线，所述阵列感应测井资料包括从不同阵列单元测量的阵列感应原始响应曲线和由所述阵列感应原始响应曲线进行合成聚焦处理后的阵列感应测井曲线；

步骤2，根据所述双侧向测井资料与辅助测井曲线资料对储层段进行分层；

步骤3，根据所述的双侧向测井曲线生成对应层的初始的原状地层电阻率、侵入带电阻率；

步骤4，根据所述的阵列感应曲线生成对应层的侵入半径值；

步骤5，根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过对模拟算法进行加速生成对应层的双侧向模拟响应曲线；

步骤6，根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过数值模式匹配法生成对应层的阵列感应模拟响应曲线；

步骤7，分别判断所述对应层的双侧向模拟响应曲线与对应层的双侧向实测曲线以及所述阵列感应模拟响应曲线与对应层的阵列感应原始响应曲线是否一致，如果一致，转到步骤8，如果不一致，转到步骤9；

步骤8，将所述原状地层电阻率、侵入带电阻率和侵入半径值作为储层段的反演结果输出；

步骤9，根据所述双侧向实测曲线和阵列感应原始响应曲线，结合对应层的地质特征、油藏特征，修改所述初始的原状地层电阻率、侵入带电阻率和侵入半径值，重复步骤5-9，输出最终的原状地层电阻率、侵入带电阻率和侵入半径值。

2.根据权利要求1所述的中低矿化度地层水储层的参数反演方法，其特征在于，所述辅助测井曲线资料包括：自然伽马测井曲线、井径测井曲线、密度测井曲线以及自然电位测井曲线。

3.根据权利要求1所述的中低矿化度地层水储层的参数反演方法，其特征在于，所述根据所述双侧向测井资料与辅助测井曲线资料对储层段进行分层，包括：

根据所述双侧向测井曲线特征依据电阻率曲线的斜率变化趋势确定层界面；

根据所述辅助测井曲线资料和确定层界面，分别将所述双侧向测井曲线特征不明显的泥岩夹层、致密层、垮塌层分出来。

4.根据权利要求1所述的中低矿化度地层水储层的参数反演方法，其特征在于，所述根据双侧向测井曲线生成对应层的初始的原状地层电阻率、侵入带电阻率包括：

对分层后的每一层，将双侧向测井曲线中的浅侧向在层中的均值作为对应层初始的侵入带电阻率；

将深侧向在层中的均值作为对应层初始的原状地层电阻率。

5.根据权利要求2所述的中低矿化度地层水储层的参数反演方法，其特征在于，所述的根据阵列感应曲线生成对应层的侵入半径值包括：

通过几何因子法对分层后的每一层计算出初步的侵入半径值：

$C_i=G_i^m{C}_m+G_i^{\mathrm{xo}}{C}_{\mathrm{xo}}^{\′}+(1-G_i^m-G_i^{\mathrm{xo}})C_t^{\′}i=1~N;$

式中C_i－阵列感应测井0.6m分辨率系列电导率值（电阻率倒数），S/m；

－井眼泥浆对阵列感应测井响应的贡献；

Cm－泥浆电导率值，S/m；

－侵入带对阵列感应测井响应的贡献，与侵入半径值有关；

C'_xo－由双侧向测井给出的侵入带初始电导率值，S/m；

C_t'－由双侧向测井给出的原状地层初始电导率值，S/m；

N－阵列感应测井0.6m分辨率系列曲线条数，且N=5或6；

通过阵列感应测井测得的每一层0.6m分辨率系列的N个值由上式可列出N个方程，通过最优化方法计算出每一层初步的侵入半径值；

其中，C_m根据所述阵列感应测井曲线计算得到，

根据所述井径测井曲线计算得到。

6.根据权利要求1所述的中低矿化度地层水储层的参数反演方法，其特征在于，所述根据对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过对模拟算法进行加速生成对应层的双侧向模拟响应曲线包括：

根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，结合双侧向测井仪器的仪器模型参数，通过有限元素法生成对应层的双侧向模拟响应曲线；其中，所述双侧向测井仪器的仪器模型参数包括：仪器长度、仪器半径、双侧向各电极相对距离、各电极长度；其中，所述的对模拟算法进行加速包括：

缩短仪器模型中屏蔽电极尺寸、降低有限元素法中计算网格长度、对计算中各电极分场进行等效处理。

7.根据权利要求1所述的中低矿化度地层水储层的参数反演方法，其特征在于，所述根据对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过数值模式匹配法生成对应层的阵列感应模拟响应曲线包括：

根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，结合阵列感应测井仪器的仪器模型参数，通过数值模式匹配法生成对应层的阵列感应模拟响应曲线；其中，

所述阵列感应测井仪器的仪器模型参数包括：仪器长度、仪器半径、各阵列单元线圈距、线圈匝数、工作频率。

8.根据权利要求1所述的中低矿化度地层水储层的参数反演方法，其特征在于，所述判断所述对应层的双侧向模拟响应曲线与对应层的所述双侧向测井曲线以及阵列感应模拟响应曲线与对应层的所述阵列感应原始响应曲线是否一致，包括：

设定对应层的双侧向模拟响应曲线与对应层的所述双侧向测井曲线以及阵列感应模拟响应曲线与对应层的所述阵列感应原始响应曲线间的相关系数均为S：

$S=\frac{\mathrm{\ΣXY}-\frac{\mathrm{\ΣX\ΣY}}{N}}{\sqrt{(\mathrm{\Σ}{X}^2-\frac{{\left(\mathrm{\ΣX}\right)}^2}{N})(\mathrm{\Σ}{Y}^2-\frac{{\left(\mathrm{\ΣY}\right)}^2}{N})}},$ 式中：

N－目的层段内的测量点数；

X－对应层的目的层段内的数据集，其中，对高阻层，X为所述对应层的双侧向测井曲线目的层段内的数据集；对中阻层，X为所述对应层的阵列感应原始响应曲线、双侧向测井曲线目的层段内的数据集；对低阻层，X为所述对应层的阵列感应原始响应曲线目的层段内的数据集；

Y－与X相对应的所述层阵列感应模拟响应曲线和/或双侧向模拟曲线目的层段内的数据集；

当所述相关系数S的计算结果大于0.9，即认为对应层的阵列感应模拟响应曲线与对应层的所述阵列感应原始响应曲线、对应层的双侧向模拟响应曲线与对应层的所述双侧向实测曲线一致。

9.根据权利要求1所述的中低矿化度地层水储层的参数反演方法，其特征在于，所述根据双侧向测井曲线和阵列感应原始响应曲线，结合对应层的地质特征、油藏特征，修改所述初始的原状地层电阻率、侵入带电阻率和侵入半径值，包括：

根据所述对应层所在井的邻井资料以及本井的录井资料、气测资料、取心资料判断所述对应层的流体性质以及所述流体性质下的地层电阻率范围；

根据所述对应层的泥浆性质以及地层属性判断侵入带电阻率范围；

根据地层浸泡时间长短及地层物性好坏判断侵入半径值范围；

根据所述对应层的双侧向模拟响应曲线与对应的双侧向实测曲线间以及阵列感应模拟响应曲线与对应的阵列感应原始响应曲线间的相关程度，并结合所述地层电阻率范围、侵入带电阻率范围以及侵入半径值范围，修改所述初始的原状地层电阻率、侵入带电阻率和侵入半径值。

10.一种中低矿化度地层水储层的参数反演装置，其特征在于，所述装置包括：

资料获取装置，用于获取双侧向测井资料、阵列感应测井资料及辅助测井曲线资料，所述双侧向测井资料包括双侧向测井曲线，所述阵列感应测井资料包括从不同阵列单元测量的阵列感应原始响应曲线和由所述阵列感应原始响应曲线进行合成聚焦处理后的阵列感应测井曲线；

分层装置，用于根据所述双侧向测井资料、辅助测井曲线资料对储层段进行分层；

电阻率确定装置，用于根据所述的双侧向测井曲线生成对应层的初始的原状地层电阻率、侵入带电阻率；

半径值确定装置，用于根据所述的阵列感应曲线生成对应层的侵入半径值；

双侧向模拟响应曲线生成装置，用于根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过对模拟算法进行加速生成对应层的双侧向模拟响应曲线；

阵列感应模拟响应曲线生成装置，用于根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过数值模式匹配法生成对应层的阵列感应模拟响应曲线；

判断装置，用于分别判断所述对应层的双侧向模拟响应曲线与对应层的双侧向实测曲线以及所述阵列感应模拟响应曲线与对应层的阵列感应原始响应曲线是否一致；

参数修改装置，用于根据所述双侧向实测曲线和阵列感应原始响应曲线，结合对应层的地质特征、油藏特征，修改所述初始的原状地层电阻率、侵入带电阻率和侵入半径值；

结果输出装置，用于将所述判断装置判断结果为一致时的原状地层电阻率、侵入带电阻率和侵入半径值作为储层段的反演结果输出。

11.根据权利要求10所述的中低矿化度地层水储层的参数反演装置，其特征在于，所述辅助测井曲线资料包括：自然伽马测井曲线、井径测井曲线、密度测井曲线以及自然电位测井曲线。

12.根据权利要求10所述的中低矿化度地层水储层的参数反演装置，其特征在于，所述分层装置包括：

确定层界面模块，用于根据所述双侧向测井曲线特征依据电阻率曲线的斜率变化趋势确定层界面；

分层模块，用于根据所述辅助测井曲线资料和确定层界面，分别将所述双侧向测井曲线特征不明显的泥岩夹层、致密层、垮塌层分出来。

13.根据权利要求10所述的中低矿化度地层水储层的参数反演装置，其特征在于，所述电阻率确定装置包括：

侵入带电阻率确定模块，用于对分层后的每一层，将双侧向测井曲线中的浅侧向在层中的均值作为对应层初始的侵入带电阻率；

原状地层电阻率确定模块，用于将深侧向在层中的均值作为对应层初始的原状地层电阻率。

14.根据权利要求11所述的中低矿化度地层水储层的参数反演装置，其特征在于，所述的半径值确定装置根据阵列感应曲线生成对应层的侵入半径值包括：

通过几何因子法对分层后的每一层计算出初步的侵入半径值：

$C_i=G_i^m{C}_m+G_i^{\mathrm{xo}}{C}_{\mathrm{xo}}^{\′}+(1-G_i^m-G_i^{\mathrm{xo}})C_t^{\′}i=1~N;$

式中C_i－阵列感应测井0.6m分辨率系列电导率值，S/m；

－井眼泥浆对阵列感应测井响应的贡献；

Cm－泥浆电导率值，S/m；

－侵入带对阵列感应测井响应的贡献，与侵入半径值有关；

C'_xo－由双侧向测井给出的侵入带初始电导率值，S/m；

C_t'－由双侧向测井给出的原状地层初始电导率值，S/m；

N－阵列感应测井0.6m分辨率系列曲线条数，且N=5或6；

通过阵列感应测井测得的每一层0.6m分辨率系列的N个值由上式可列出N个方程，通过最优化方法计算出每一层初步的侵入半径值；

其中，C_m根据所述阵列感应测井曲线计算得到，

根据所述井径测井曲线计算得到。

15.根据权利要求10所述的中低矿化度地层水储层的参数反演装置，其特征在于，所述的双侧向模拟响应曲线生成装置根据对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过对模拟算法进行加速生成对应层的双侧向模拟响应曲线包括：

根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，结合双侧向测井仪器的仪器模型参数，通过有限元素法生成对应层的双侧向模拟响应曲线；其中，所述双侧向测井仪器的仪器模型参数包括：仪器长度、仪器半径、双侧向各电极相对距离、各电极长度；其中，所述的对模拟算法进行加速包括：

缩短仪器模型中屏蔽电极尺寸、降低有限元素法中计算网格长度、对计算中各电极分场进行等效处理。

16.根据权利要求10所述的中低矿化度地层水储层的参数反演装置，其特征在于，所述的阵列感应模拟响应曲线生成装置根据对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，通过数值模式匹配法生成对应层的阵列感应模拟响应曲线包括：

根据所述的对应层的初始的原状地层电阻率、侵入带电阻率和侵入半径值，结合阵列感应测井仪器的仪器模型参数，通过数值模式匹配法生成对应层的阵列感应模拟响应曲线；其中，

所述阵列感应测井仪器的仪器模型参数包括：仪器长度、仪器半径、各阵列单元线圈距、线圈匝数、工作频率。

17.根据权利要求10所述的中低矿化度地层水储层的参数反演装置，其特征在于，所述判断装置判断对应层的双侧向模拟响应曲线与对应层的所述双侧向测井曲线以及阵列感应模拟响应曲线与对应层的所述阵列感应原始响应曲线是否一致，包括：

设定对应层的双侧向模拟响应曲线与对应层的所述双侧向测井曲线以及阵列感应模拟响应曲线与对应层的所述阵列感应原始响应曲线间的相关系数均为S：

$S=\frac{\mathrm{\ΣXY}-\frac{\mathrm{\ΣX\ΣY}}{N}}{\sqrt{(\mathrm{\Σ}{X}^2-\frac{{\left(\mathrm{\ΣX}\right)}^2}{N})(\mathrm{\Σ}{Y}^2-\frac{{\left(\mathrm{\ΣY}\right)}^2}{N})}},$ 式中：

N－目的层段内的测量点数；

X－对应层的目的层段内的数据集，其中，对高阻层，X为所述对应层的双侧向测井曲线目的层段内的数据集；对中阻层，X为所述对应层的阵列感应原始响应曲线、双侧向测井曲线目的层段内的数据集；对低阻层，X为所述对应层的阵列感应原始响应曲线目的层段内的数据集；

Y－与X相对应的所述层阵列感应模拟响应曲线和/或双侧向模拟曲线目的层段内的数据集；

当所述相关系数S的计算结果大于0.9，即认为对应层的阵列感应模拟响应曲线与对应层的所述阵列感应原始响应曲线、对应层的双侧向模拟响应曲线与对应层的所述双侧向实测曲线一致。

18.根据权利要求10所述的中低矿化度地层水储层的参数反演装置，其特征在于，所述参数修改装置根据双侧向测井曲线和阵列感应原始响应曲线，结合对应层的地质特征、油藏特征，修改所述初始的原状地层电阻率、侵入带电阻率和侵入半径值，包括：

根据所述对应层所在井的邻井资料以及本井的录井资料、气测资料、取心资料判断所述对应层的流体性质以及所述流体性质下的地层电阻率范围；

根据所述对应层的泥浆性质以及地层属性判断侵入带电阻率范围；

根据地层浸泡时间长短及地层物性好坏判断侵入半径值范围；

根据所述对应层的双侧向模拟响应曲线与对应的双侧向实测曲线间以及阵列感应模拟响应曲线与对应的阵列感应原始响应曲线间的相关程度，并结合所述地层电阻率范围、侵入带电阻率范围以及侵入半径值范围，修改所述初始的原状地层电阻率、侵入带电阻率和侵入半径值。

Images