CN102066983A - 套管井电测井方法 - Google Patents

Abstract

本发明涉及井的地球物理研究领域，可用于测定金属套管井周围岩层的电阻率。按照本发明的套管井的电测方法，采用五个电极的探头，其中三个测量电极沿着井的轴线依次等距离放置，另外两个电流电极在中心测量电极外侧并相对其对称分布，对电极依次供给恒定电流双极性矩形脉冲，每次供给电流后测量中心测量电极的电场电位与外侧测量电极之间一阶电位差分。对所有信号的数字量以相应方法进行处理、滤波，根据这些数字量的处理滤波后的数值，计算套管井周围岩层的电阻率。

Description

套管井电测井方法

技术领域

本发明属于井的地球物理研究领域，可用于测定金属套管井周围岩层的电阻率。

背景技术

现有的套管井电测井方法是采用对称的双极性五电极探头，对电场电位及一阶、二阶电位差分进行测量(КашикА.С.，РыхлинскийН.И.等，套管井电测井方法，专利号2176802，申请日2001年2月20日，刊登在2001年第34期公报上)。在该方法中，通过在探头的电流电极中保持恰好能使测量点出现电位极值的电流来消除由于套管井单位长度的电阻不稳定而出现的测量误差。该方法的缺点是，当向探头的电流电极供给电流时，由于电力电缆从接收电极的电线旁边通过，从而在接收电路上引起诱导感应，以至使套管周围岩层的电阻率的测量动态范围大大降低至25Ω.m，误差超过10％。

使用现有的采用对称的双极性五电极探头的套井管电测井方法(КашикА.С.，РыхлинскийН.И.等，套管井电测井方法，专利号2229735，申请日2003年4月22日，刊登在2004年第15期公报上)，虽然消除了诱导感应。但是，在使用此种方法时，为了保持电位的极值，需要通过范围为微毫伏特的有效信号对位于钻井仪器上自动模拟调节器进行控制才能实现。由于热噪声、诱导感应、大地电流、接触电极电位等原因，而使有效信号远远低于干扰信号，使得自动模拟调节器无法稳定工作，从而难以达到控制的作用。因此这种测量方法在实际上无法用于套管井的电测量。

由此可见，任何一种套有整体金属套管的井电测方法均是要在毫微伏特有效信号的工作条件下进行的，如果不抑制干扰信号，则有效信号远远低于干扰信号。

同时还可以看出，上述所有测量方法都是建立在两组探头测量电极之间的二阶电位差分的测量基础上，第一方案是借助于由两个相同的电阻组成的电桥，另一种方案是对两组电位差单独测量，然后将输出端上显示的两组测量读数相减。第二种方案仍然存在着不足，因为在技术上很难制造出两个具有相同且稳定的放大倍数的放大器，该放大器适用于二阶电位差分的差动测量。

在两组测量电极间的基线保持恒定的情况下，采用这些二阶电位差分的测量方法对套管井进行电测量原则上是可行的。但是由于套管井的电检测装置的测量电极装置具有压紧结构，套管的直径是易变的，例如：由于打孔后的突起，或者由于轧材的特殊工艺，当电极压紧直径改变时，两组测量电极之间的距离可能会扩大到一厘米。由于套管的电阻率与周围岩层的电阻率相差10⁷倍甚至更多，则在电位的二阶差分的差动测量中，测量电极之间的这些距离变化可能导致超过容许误差许多倍的误差。

发明内容

在本发明的方法中，解决了在有效测量信号中去除这些干扰的问题，随之将井周围岩层的实际电阻率测定的动态范围提高100Ω.m，测量误差为5％。

这一技术效果通过以下方式获得：采用五电极探头套管井电测井方法，其中三个测量电极沿着井的轴线依次等距离放置，而在外侧相对于中心测量电极对称放置两个电流电极，向所述两个电流电极依次供给恒定电流双极性矩形脉冲，每次供给电流变极后经过一指定时间，测量中心电场的电位和两个外侧测量电极之间的一阶电位差分；在以上电子测量信号的基础上，测量井周围岩层的电阻率；这些电信号是在由于两个电极的和作用，而使最终得到的两个外侧测量电极之间的电场电位差分等于零的条件下测得的。

根据本发明，当向两个电流电极之一供给电流时，进行附加测量位于任一外侧测量电极与中心电极之间的一阶电位差分；

双极性脉冲的频率为0.25赫兹或低于0.25赫兹；以数字化频率5赫兹或高于5赫兹来数字化电位、及其一阶电位差分和两个电流电极的电流数值，另外，对一阶电位差分和电流进行数字化的开始时间不早于电流变极后的0.4秒，而对电位数字化的开始时间不早于电流变极后的1秒；

将每个位于某一标号下的电流的双极性脉冲正半周期所有信号数字化值从位于相应的标号下的负半周期的数字化值中减去；

每个获得的下面具有某一标号的电位及其一阶差分数字化的量除以具有相应的标号的电流无量纲的模。

将获得的所有数字化后的电位差分和根据电流值的模规格化后的电位及其一阶差分值的数值组通过高频滤波器进行滤波，借此使热噪声、大地电流和超标严重的量的影响降为最低。

对滤波后的电位值及其一阶电位差分值的数据组和电流值的数据组电流值分别进行求和后再取平均值，利用上述的两个平均值根据以下公式计算井周围岩层的电阻率ρ_n：

式中： ${\mathrm{\ρ}}_n=k\·{\mathrm{\Ω}}_z\·\left[\frac{U_N\left(I_{A1}\right)/\left|I_{A1}\right|+k\·U_N\left(I_{A2}\right)/\left|I_{A2}\right|}{\mathrm{\Δ}{U}_{M1N}\left(I_{A1}\right)/\left|I_{A1}\right|+k\·\mathrm{\Δ}{U}_{M1N}\left(I_{A2}\right)/\left|I_{A2}\right|}\right]$

Ω_z是探头外侧测量电极M1和M2之间的井段的电阻率；

k是会聚参数，其确定的条件是探头外侧测量电极间的规格化电位差分的合成后的话结果必须为零，其计算方程如下：

$\frac{{\mathrm{\ΔU}}_{M1M2}\left(I_{A1}\right)}{\left|I_{A1}\right|}+k\·\frac{\mathrm{\Δ}{U}_{M1M2}\left(I_{A2}\right)}{\left|I_{A2}\right|}\text{=0;}$

U_N(I_A1)、U_N(I_A2)是探头中心测量电极经过数字滤波和求平均值后得到的电位，它们分别取决于探头的第一和第二电流电极A1和A2的电流；

ΔU_M1N(I_A1)、ΔU_M2N(I_A2)是由一个测量器测出的位于两个外侧测量电极M1和中心电极N之间的经数字滤波和求平均值后的电场电位差分，它们分别取决于探头的第一和第二电流电极A1和A2的电流；

ΔU_M1M2(I_A1)、ΔU_M2M2(I_A2)是由一个测量器测出的位于两个外侧测量电极M1和M2之间的经数字滤波和求平均值后的电场电位差分，它们分别取决于探头电流电极A1和A2的电流；

I_A1、I_A2是探头电流电极A1和A2的电流；

|I_A1|、|I_A2|是电流I_A1、I_A2无量纲的模；

k是探头的几何参数。

在套管井电测井技术中，根据本发明的内容，探头的几何参数k是通过数学网状模型计算出的，从而也获得了技术效果。

附图说明

具体实施方式

图1示出了实现本发明方法的设备的框图，其中，1是套管，2是井周围的岩层，3-M1和4-N是两个电极间一阶电位差分测量传感器电极，3-M1和5-M2是两个电极间一阶电位差分测量传感器电极，6-A1和7-A2是位于测量电极边界外相对于中心测量电极N对称放置的探头电流电极，8是一阶电位差分的数字测量器，其输入端与外侧测量探头测量电极M1及M2相连接，9是一阶电位差分的数字测量器，其输入端与探头外侧测量电极M1及探头中心测量电极N相连接，10是探头中心测量电极N的电位数字测量器，也是位于套管1最上端的测量电极17-Nуд的测量器，11和13是探头电流电极6-A1和7-A2电路上的分流器，其作用是测量流经这些电极的电流强度，12和14是流经电极6A1和7-A2的电流数字测量器，15是由有阳光的地面操控的探头电流电极6-A1和7-A2的电流换向器，16是通过电缆从数字测量器8，9，10，12，14输出端向地面电子设备21远距离输送测量数据的电子钻井设备，18是恒电流双极性矩形脉冲的地面电源，以大约5安培的电流强度向探头电流电极供电，19是地面编程设备，用于对电流转换器15和向电流电极6-A1和7-A2供给电流进行控制，20是接到地面任意点的返回电流电极B，22是钻井设备控制处理器，用于处理和滤波所有测量信号并计算出井周围岩层的电阻率。模拟数字转换器在该实施例中的实际数字化频率为5赫兹，23是用于从电位测量器10的输入端连接远处电极Nуд(套管口)的电缆芯线，24是电缆电流芯线，25是用于从电子钻井设备16向地面电子装置输送数字测量信号的电缆芯线。

图2-a示出了当外侧电极M1和M2之间一阶电位差分ΔU(I_A1，I_A2)为零时，电流电极A1和A2间的电流分布图。图2-б示出了沿套管由来自电流电极A1的工作电流I_A1产生的电位分布图。图2-B示出了沿套管由来自电流电极A2的工作电流I_A2产生的电位分布图，该电位分布是在考虑到确定绘制比例的会聚参数k，并且是根据在两个外侧测量电极M1和M2间的一阶电位差分ΔU(I_A1，I_A2)等于零的条件下计算出的。

通过以上情况可以看出，任何一种整体金属套管电测井方法都是用于有效信号在毫微伏特范围内的工作环境，远远低于以下其他各种干扰：由于套管内直径的改变而使得探测器测量电极间的距离发生改变出现的干扰，由于以上电极压紧接触杆倾角的改变而产生的干扰；由于套管单位长度的可变电阻而产生的干扰；由于对探头电流电极供电不稳定所产生的干扰，其原因是供给电源无法满足对毫微伏特范围内的工作进行稳定的供给，或者是由于电路电阻不稳定；由于探头电流电极供电线路的诱导感应而产生的干扰；由于接触的电极电位而产生的干扰；热激干扰；大地干扰；随机出现的脉冲干扰。

为了避免上述干扰，待测量电位、一阶电位差分和两电流电极的电流都是根据数字化频率5赫兹或以上而确定的。

为了成功地从热激干扰、大地干扰和其他随机的干扰中过滤出有效信号，需要对探头电流电极供给高频率的电流，但由于诱导电感的影响，不得超过0.25赫兹。为了消除因电流线路与测量线路重叠而引起的诱导感应的失真影响，优选的是，通过符号变更的恒电流矩形脉冲来实现对探头电流电极的供给，所述的诱导电感干扰，在恒电流矩形脉冲变极后过一段时间会消失。在这种情况下，必须在与探头电流电极电流变极过程中有关的抽样信号停止后，再开始对测量电路上的信号进行测量和数字化(需要指出的是，数字化频率为5赫兹和交流电符号变更矩形脉冲的频率为0.25赫兹时，在每一个电流矩形脉冲里均匀分配10个数字量)。电流变极与一阶电位差分测量开始点之间的时间间隔，如试验研究所显示那样，取决于互相重叠的电流线路和测量线路的长度。测量一阶电位差分的间隔不低于0.4秒(相当于两个数字化量的频率不超过5赫兹)，由于一阶电位差分的电流线路和测量线路的重叠的距离总共只有几米，这个距离等于一阶电位差分的测量线路的长度。为避免诱导电感失真对一阶电位差分测量结果的影响，不使用来自这两个数字化量的信息。

在测量电位时，这一间隔已经不小于1秒(与同样的数字化频率为5赫兹时的两个数字量相符合)，因为在这种情况下，电流线路与测量电路在相距几千米的间隔内重叠，即连接测井仪器与地面设备的测井电缆的总长度。由此，双极性矩形脉冲的最佳频率为0.25赫兹或更低。

为了消除探头电极的接触电位，需要将电流双极性脉冲的正半周期范围内的具有某个下标编号的所有信号中的每一个数字量从负半周期范围内的具有相应下标编号的数字量中减去。

为了将供给电流电极的电流的不稳定影响降至最低，需要将每一个所得到的具有某个下标编号的电位及其一阶和一阶电位差分的数字化数值的差分值除以(规格化)具有相应下标编号的电流数字化数值的无量纲的模。

由于高导电率套管对探头的屏蔽作用，而使二阶电位差分放大器的输入端的二阶电位差分信号降至毫微伏特，因此必须设法去掉干扰中最强的一种，即热激干扰，对此，可以通过对被测量的信号数字量进行高频滤波。

将所获得的所有的数字化的差分和根据电流模规格化的电位差及其一阶和二阶电位差分值的数值组借助高频滤波器使热噪声、大地电流及其他严重超标数值的影响降到最低。高频滤波的方法如下：从每一个频率为0.25赫兹的循环中分离出差分的和按照具有相同下标编号的电流值的模进行规格化的数值，然后例如通过确定中线的方法进行滤波(Г.Корн，Т.Корн.数学参考，《科学》，莫斯科，1974,545页[3])。在经所述的滤波后，凡属于同一下标编号的电位及其一阶差分和电流始终按照同一个滤波量，所述的电位及其一阶差分和电流不取决于电流矩形脉冲供电循环的数目。

然后将留下滤波后的电位数值及其一阶电位差电流数值分别地求和，再求平均数，就是将得到的总数的和除以被求和的数值的个数。

现在研究一下如何去除套管电阻不稳定性的失真影响的原理。

在解决石油天然气勘探作业电测法中，把研究的介质近似为沿坐标轴Z和R两维变化的实体。同时井并不是一个理想的线性电极，即沿着坐标轴Z的，外侧测量电极间的线电阻Ω_Z是不固定的(Ω_Z≠const)，两段之间的线电阻可能相差几倍。

将电源放入井中的A点，通过该电源向所研究的向介质供给电流I来确定电位沿井的轴线分布。

在[1]中给出， $I_Z=-\frac{1}{{\mathrm{\Ω}}_Z}\·\frac{\mathrm{dU}\left(z\right)}{\mathrm{dz}}---\left(1\right)$

只有的当Ω_r/Ω_Z＞＞1时，J_r(z)＝U(z)/Ω_r  (2)

式中：

U(z)是坐标为Z的井内观测点上的电位

I_Z(z)是流过同一坐标Z的井的横截面的电流

J_r(z)是从井壁流向周围岩层的单位距离深度上的电流(电流线性密度为A/M)

Ω_r电流J_r(z)流过的介质所表现出的电阻

Ω_Z(如上所述那样)是外侧测量电极间的井段轴向的电流电阻，由于井的几何参数和其他参数的易变性，其值是坐标Z的函数。

在Z点从井中划分出一段中间观察点的一段高度为ΔZ的圆柱体，该Z点是圆柱的中心，并且是观测点(中心测量电极N)，对于这一段圆柱的封闭表面应用电流密度向量

的连续性方程积分的形式，

即：

S表面由圆柱的底S_p、S_q和其侧面S_b的表面组成，方程式(3)左边表示三个通量之和：

根据方程式(3)，得出：

I_Z(z+Δz/2)-I_Z(z-Δz/2)+J_r(z)·Δz＝o(Δz)    (4)

ΔI_Z(z)/Δz＝-J_r(z)+o(1)并且在范围Δz→0的条件下

从： $J_r\left(z\right)=-\frac{{\mathrm{dI}}_Z\left(z\right)}{\mathrm{dz}}---\left(5\right)$

将表达式(1)对Z求导数，考虑到Ωz是套管井体的电阻的函数，在实际的井中，该电阻随座标Z的改变而变化，即Ω_Z＝Ω_Z(z)≠const：

$\frac{{\mathrm{dI}}_Z\left(z\right)}{\mathrm{dz}}=\frac{1}{{\mathrm{\Ω}}_Z^2}\·\frac{{\mathrm{d\Ω}}_Z}{\mathrm{dz}}\·\frac{\mathrm{dU}\left(z\right)}{\mathrm{dz}}-\frac{1}{{\mathrm{\Ω}}_Z}\·\frac{d^{2}U\left(z\right)}{{\mathrm{dz}}^2}---\left(6\right)$

将等式(2)和(5)代入等式(6)中，获得电源电位沿井的轴线的分布的方程，在该方程中含有沿着井的轴线变化的电阻Ω_Z[1]

$\frac{d^{2}U\left(z\right)}{{\mathrm{dz}}^2}-\frac{1}{{\mathrm{\Ω}}_Z}\·\frac{{\mathrm{d\Ω}}_Z}{\mathrm{dz}}\·\frac{\mathrm{dU}\left(z\right)}{\mathrm{dz}}-\frac{{\mathrm{\Ω}}_Z}{{\mathrm{\Ω}}_r}\·U\left(z\right)=0---\left(7\right)$

对等式(7)分析表明，测量电位及其二阶导数由于等式中有dΩ_Z/dz项而无法确定所求的比值Ω_Z/Ω_r，dΩ_Z/dz强烈地依赖于井体电阻易变性。

在电测井方法[1]中，单位套管体电阻的不固定性实际上对测量结果不会产生影响，该方法的特点是，由于采用相应技术措施和手段，使得沿井的轴线线的电位分布曲线在测量电极区内获得极值(在坐标Z＝Z_N区域)，即dU(Z_N)/dz＝0。随之，从方程式中消除了包含有不确定大小的成分dΩ_Z/dz，这一方程式在点Z＝Z_N处采用以下的形式：

$U\left(z_N\right)/\frac{d^{2}U\left(z_N\right)}{{\mathrm{dz}}^2}={\mathrm{\Ω}}_r/{\mathrm{\Ω}}_Z---\left(8\right)$

由此求出：

${\mathrm{\Ω}}_r={\mathrm{\Ω}}_Z\·U\left(z_N\right)/\frac{d^{2}U\left(z_N\right)}{{\mathrm{dz}}^2}---\left(9\right)$

根据方程式(9)，通过测量外侧测量电极间圆柱体段中的坐标为Z_N上的电阻，电位、及其二阶导数来计算井周围电阻Ω_r，所述的坐标为Z_N的点的U(Z_N)具有极值。

为了找到测量电极处电位极值的位置而要对位于中心电极N(测量点)相同距离的两个方向上的电源A1和A2，选取合适的电流，使这两个电流值相对于N的2个对称电极M1和M2的电位差为零，即：

在测量点z＝z_N处获得极值就意味着消除了井中电流I_Z(z_N)的轴向分量，在用单极性电源激发所研究的介质条件下，这个轴向分量远高于径向分量J_r(z_N)。在实际中，为了测量电阻Ω_r，而使用与该电位的二阶导数成正比的电位有限的二阶差分代替方程式(9)中的二阶电位导数

${\mathrm{\Δ}}^{2}U\left(z_N\right)=U_{M_1}+U_{M_2}-2U_N~\frac{d^{2}U\left(z_N\right)}{{\mathrm{dz}}^2}---\left(10\right)$

这种电测套管井的方法可以在电流会聚在信号的接收点上的条件下，测量电阻Ω_r，即如果在探头中心电极N点上能保证电场电位U(z)(dU_N/dz＝0)维持极值的条件下。根据欧姆定律，在这一点上，沿井的轴线线的电流密度的轴向分量为零

根据公式(9)和(10)在确定套管井周围岩层的电阻率的基础上实现本发明的方法。

${\mathrm{\ρ}}_n=k\·{\mathrm{\Ω}}_Z\·\frac{U_N\left(z_N\right)}{{\mathrm{\Δ}}^2{U}_N\left(z_N\right)}---\left(11\right)$

上述的公式(11)必需在保证实现下述条件下才成立：保证由这对电流电极在外侧测量电极M1和M2之间产生的一阶电位的差分ΔU(I_A1，I_A2)的合成后结果等于零伏特。

式中U_N(z_N)和Δ²U_N(z_N)分别是电极N的电场电位和外侧测量电极M1和M2之间导电圆柱的区段上的电位二阶差分，所述的电位和所述的电位的二阶差分是在这两个电极间的合成过的一阶电位差分等于零伏特的条件下测得的。

k是探头的几何参数，米

在数字存储时，由于探头电流电极供电电源的不稳定性，必需考虑将被测信号的数字量根据供给探头电流电极的相应的数字化量的模进行规格化，这时需将公式(11)变为：

${\mathrm{\ρ}}_n=2k\·{\mathrm{\Ω}}_Z\·\left[\frac{U_N\left(I_{A1}\right)/\left|I_{A1}\right|+k\·U_N\left(I_{A2}\right)/\left|I_{A2}\right|}{{\mathrm{\Δ}}^{2}U\left(I_{A1}\right)/\left|I_{A1}\right|+k\·{\mathrm{\Δ}}^{2}U\left(I_{A2}\right)/\left|I_{A2}\right|}\right]---\left(12\right)$

k是会聚参数，它是在必需保证探头的两个外侧测量电极间电位差的规格化后的总和等于零的条件下根据方程(13)计算出的。

$\frac{\mathrm{\Δ}{U}_{M1M2}\left(I_{A1}\right)}{{|I}_{A1}|}+k\·\frac{\mathrm{\Δ}{U}_{M1M2}\left(I_{A2}\right)}{\left|I_{A2}\right|}=0---\left(13\right)$

式中U_N(I_A1)，U_N(I_A2)分别取决于探头第一电流电极A1和第二电流电极A2的探头的中心测量电极的电场的电位，所述的两个电位是在对它们的数值滤波求和并取平均得到的。

ΔU(I_A1)，ΔU(I_A2)，Δ²U(I_A1)，Δ²U(I_A2)分别是电场一阶和二阶电位差分，它们取决于探头第一电流电极A1的电流和第二电流电极A2的电流，是经过对滤波后的数值进行求和及取平均获得的。

I_A1、I_A2是探头电流电极A1、A2的电流。

|I_A1|，|I_A2|是电流I_A1和I_A2无量纲的模，是经过对滤波后的数值进行求和及求平均获得的。

k是探头的几何参数。

Ω_Z是在探头的两个外侧测量电极之间测量的井段上的电阻。

探头外侧测量电极之间的井段上的电阻Ω_Z通常通过公式(14)计算出

${\mathrm{\Ω}}_Z=\frac{\mathrm{\ΔU}\left(I_{A1}\right)}{I_{A1}}-\frac{\mathrm{\ΔU}\left(I_{A2}\right)}{I_{A2}}---\left(14\right)$

电阻率ρ_n在具体实施中通过公式(12)计算出。正如上面已指出的那样，推导这一公式的前提条件是：流过测量电极3-M1，5-M2之间的高传导率金属套管的总的轴向电流分量等于零。因此，消除了套管电阻不稳定性对测量结果的失真影响，中央处理器在对信号进行处理后，根据公式(12)计算岩精确的电阻率，经数学模拟机的模拟确认无误。

然而，正如上面所指出的那样，当套管周围岩层电阻率与套管的电阻率的比率为107或更高时(实际中这样的比例关系总是存在)，两个测量电极间的二阶电位差分的差动测量不论是用由两个具有相同电阻组成的桥；还是用两个单独的测量器单独测量电位差，然后再将输出端的两个读数相减；这两种方法都不能达到测量所必须的精确度。因此应当从公式(12)中消去包含二阶电位差分的项Δ²U(I_A1)，Δ²U(I_A2)。为此应利用图2和公式(13)，从中得出：借助于参数k，

U_M1(I_A1)＝k·U_M2(I_A2)                        (15)

U_N(I_A1)＝k·U_N(I_A2)                          (16)

U_M2(I_A1)＝k·U_M1(I_A2)                        (17)

从公式(12)中消去含有二阶电位差分的Δ²U(I_A1)，Δ²U(I_A2)的项(为简化分析，我们将|I_A1|，|I_A2|认定为一个单位)。根据公式(10)、(15-17)得出：

Δ²U(I_A1)+k·Δ²U(I_A2)＝

U_M1(I_A1)+U_M2(I_A1)-U_N(I_A1)-U_N(I_A1)+

k·U_M1(I_A2)+k·U_M2(I_A2)-k·U_N(I_A2)-k·U_N(I_A2)＝       (18)

2{[U_M1(I_A1)-U_N(I_A1)]+k·[U_M1(I_A2)-U_N(I_A2)]}＝

2[ΔU_M1N(I_A1)＝k·ΔU_M1N(I_A2)]

为了计算电阻率ρ_n，根据公式(18)和(12)得到：

${\mathrm{\ρ}}_n=k\·{\mathrm{\Ω}}_z\·\left[\frac{U_N\left(I_{A1}\right)/\left|I_{A1}\right|+k\·U_N\left(I_{A2}\right)/\left|I_{A2}\right|}{\mathrm{\Δ}{U}_{M1N}\left(I_{A1}\right)/\left|I_{A1}\right|+k\·\mathrm{\Δ}{U}_{M1N}\left(I_{A2}\right)/\left|I_{A2}\right|}\right]---\left(19\right)$

公式(19)对于测量电阻率ρ_n在参数上与公式(12)没有区别，其本质的不同在于：，在公式(19)中差动第二测量电位差Δ²U(I_A1)，Δ²U(I_A2)被替换为由同一测量器完整测量的在外侧测量电极之一的M1和中心电极N之间的第一测量电位差分ΔU_M1N(I_A1)，ΔU_M1N(I_A2)。因此，使电阻率ρ_n测量结果的准确性大大提高。在使用稳定性较高的探头电流电极供给电源的条件下，当I_A1＝I_A2＝const时，将公式(19)简化为：

${\mathrm{\ρ}}_n=k\·{\mathrm{\Ω}}_Z\·\left[\frac{U_N\left(I_{A1}\right)+k\·U_N\left(I_{A2}\right)}{\mathrm{\Δ}{U}_{M1N}\left(I_{A1}\right)+k\·\mathrm{\Δ}{U}_{M1N}\left(I_{A2}\right)}\right]---\left(20\right)$

在本发明的方法中，探头几何参数k、电阻率ρ_n和依本方法生产的设备的读数之间的线性区域的计算是借助网状数学模型完成的(B.Друскин，Л.Книжнерман，恒流电测井与电勘探方法，苏联科学院简报，《地球物理》系列，1987年第4期，63-71页)[4]，将已知的Ω_Z、ρ_n、I_A1、I_A2

代入公式(19)中，就可以计算出电位数值U_N(I_A1)，U_N(I_A2)和电位差分ΔU_M1N(I_A1)，ΔU_M1N(I_A2)，ΔU_M1M2(I_A1)，ΔU_M1M2(I_A2)。

按照本发明的方法制造的设备现已经过了井中试验，在套管井电测井时电阻率ρ_n的误差不超过5％。

Claims (2)

1.一种五电极探头套管井电测井方法，采用五电极探头套管井电测井方法，其中三个测量电极沿着井的轴线依次等距离放置，而在外侧相对于中心测量电极对称放置两个电流电极，向所述两个电流电极依次供给恒定电流双极性矩形脉冲，每次供给电流变极后经过一指定时间，测量中心电场的电位和两个外侧测量电极之间的一阶电位差分；在以上电子测量信号的基础上，测量井周围岩层的电阻率；这些电信号是在由于两个电极的和作用，而使最终得到的两个外侧测量电极之间的电场电位差分等于零的条件下测得的；其特征在于：

当向两个电流电极之一供给电流时，进行附加测量位于任一外侧量电极与中心电极之间的一阶电位差分；

双极性脉冲的频率为0.25赫兹或低于0.25赫兹；以数字化频率5赫兹或高于5赫兹来数字化电位、及其一阶电位差分和两个电流电极的电流数值，另外，对一阶电位差分和电流进行数字化的开始时间不早于电流变极后的0.4秒，而对电位数字化的开始时间不早于电流变极后的0.4秒；

将每个位于某一标号下的电流的双极性脉冲正半周期所有信号量子数值从位于相应的标号下的负半周期的量子数值中减去；

每个获得的下面具有某一标号的电位及其一阶差分数字化的量除以具有相应的标号的电流无量纲的模；

将获得的所有数字化后的电位差分和根据电流值的模规格化后的电位及其一阶差分值的数值组通过高频滤波器进行滤波，借此使热噪声、大地电流和超标严重的量的影响降为最低；

对每一个滤波后的电位值及其一阶电位差分值的数据组和电流值的数据组分别进行求和后再取平均值，利用上述的两个平均值根据以下公式计算井周围岩层的电阻率ρ_n：

${\mathrm{\ρ}}_n=k\·{\mathrm{\Ω}}_z\·\left[\frac{U_N\left(I_{A1}\right)/\left|I_{A1}\right|+k\·U_N\left(I_{A2}\right)/\left|I_{A2}\right|}{\mathrm{\Δ}{U}_{M1N}\left(I_{A1}\right)/\left|I_{A1}\right|+k\·\mathrm{\Δ}{U}_{M1N}\left(I_{A2}\right)/\left|I_{A2}\right|}\right]$

式中：

Ω_Z是探头外侧测量电极M1和M2之间的井段的电阻率；

k是会聚参数，其确定的条件是探头外侧测量电极间的规格化电位差分的合成后的结果必须为零，其计算方程如下：

$\frac{\mathrm{\Δ}{U}_{M1M2}\left(I_{A1}\right)}{{|I}_{A1}|}+k\·\frac{\mathrm{\Δ}{U}_{M1M2}\left(I_{A2}\right)}{\left|I_{A2}\right|}=0$

U_N(I_A1)、U_N(I_A2)是探头中心测量电极经过数字滤波和求平均值后得到的电位，它们分别取决于探头电流电极A1和A2的电流；

ΔU_M1N(I_A1)、ΔU_M2N(I_A2)是由一个测量器测出的位于外侧测量电极M1和中心电极N之间的经数字滤波和求平均值后的电场电位差分，它们分别取决于探头电流电极A1和A2的电流；

ΔU_M1M2(I_A1)、ΔU_M2M2(I_A2)是由一个测量器测出的位于外侧测量电极M1和中心电极M2之间的经数字滤波和求平均值后的电场电位差分，它们分别取决于探头电流电极A1和A2的电流；

I_A1、I_A2是探头电流电极A1和A2的电流；

|I_A1|、|I_A2|是电流I_A1、I_A2无量纲的模；

k是探头的几何参数。

2.如权利要求1所述的套管井电测井方法，其特征在于：通过数学网状模型来计算探头的几何参数k。

Images