CN100374879C - 确定横向或三轴阵列感应或传播测井仪的探头误差的方法 - Google Patents

确定横向或三轴阵列感应或传播测井仪的探头误差的方法 Download PDF

Info

Publication number
CN100374879C
CN100374879C CNB2004101007593A CN200410100759A CN100374879C CN 100374879 C CN100374879 C CN 100374879C CN B2004101007593 A CNB2004101007593 A CN B2004101007593A CN 200410100759 A CN200410100759 A CN 200410100759A CN 100374879 C CN100374879 C CN 100374879C
Authority
CN
China
Prior art keywords
measurement
array
height
difference
tool
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CNB2004101007593A
Other languages
English (en)
Other versions
CN1648691A (zh
Inventor
T·D·巴伯
H·王
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Schlumberger Overseas SA
Original Assignee
Schlumberger Overseas SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Schlumberger Overseas SA filed Critical Schlumberger Overseas SA
Publication of CN1648691A publication Critical patent/CN1648691A/zh
Application granted granted Critical
Publication of CN100374879C publication Critical patent/CN100374879C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G01V3/26—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
    • G01V3/28—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device using induction coils
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G01V3/30—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electromagnetic waves

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geophysics (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Electromagnetism (AREA)
  • Manufacturing & Machinery (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Current Or Voltage (AREA)
  • Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)

Abstract

一种用于确定具有横向阵列的测井仪内的探头误差的方法，其包括：当测井仪处于以所选择的方向距地面第一高度时，使用横向阵列获得第一测量值；当测井仪处于以所选择的方向距地面第二高度时，使用横向阵列获得第二测量值；由第一测量值和第二测量值获得差值测量值；利用差值测量值和预设函数确定背景信号，其中所述的预设函数在一组地层电导率处将在第二高度处所得的一组信号内的地信号与在第一高度和第二高度处所得的信号之间的一组差值信号联系起来；以及通过从第二测量值中减去背景信号来确定探头误差。

Description

确定横向或三轴阵列感应或传播测井仪的探头误差的方法
技术领域
本发明通常涉及电测井领域。尤其是，本发明涉及在电磁感应和传播测井仪中测定探头误差的方法。
背景技术
电磁(EM)感应和传播测井仪已经应用了多年，用来测量井筒周围地层的电导率。EM测井仪通过响应于AC发射器信号而在地层内感应出涡流从而测量地层的电阻率(或反之，电导率)。涡流感应出第二磁场，第二磁场反过来又在接收器天线感应出电压。由于涡流的强度取决于地层的电导率，因此所接收信号的强度反映出地层的电导率。
理想的是，感应或传播测井仪在零电导率介质中的读数应该为零。然而，EM测井仪包括传导材料(如探头配线、四轴(quadraxes)、隔板、安装在感应探头上用于的球形聚焦测井的电极，等等)，其可响应于发射器所产生的磁场。在这些金属部件内感应出的涡流在接收器内产生稳定的信号。因此，感应或传播测井仪在零电导率介质中的读数将不会为零，这一自生信号被称为探头误差。通常，探头误差能够通过在远离任何外部传导性材料的自由空间内悬挂所述的仪器而被测定。然而，大多数的感应或传播测井仪用于测量地层深处(即深的探测深度)，以测定侵入带以外的地层电阻率。这些“深层读数”的测井仪使自由空间的探头误差测定不再适用。由于这个原因，用于探头误差校正的方法已经在美国专利US 4,800,496中公开了，该专利授权给Barber等人(“Barber专利”)并转让给本发明的受让人。这一专利整体引入作为参考。
Barber专利中所公开的方法通过距地面两种不同的距离进行测量来对常规的感应测井仪进行探头误差校正。之后，一种算法将两个高度处的电压的差异与地信号相关联。最后，将地信号从探头读数中减去，以获得真实的探头误差。然后，将真实的探头误差用于校正测井测量值。
常规的感应测井仪具有以一种方式设置的发射器和接收器，以便它们的磁矩与测井仪的纵向轴线对齐。这些纵向感应阵列测井仪在垂直于测井仪纵向轴线的线圈内感应出涡流。常规的感应测井仪不能提供精确的各向异性地层的电阻率估测值。地层各向异性是由地层自然沉积的方式所决定的。含油气的地层通常在地层电阻率方面显示各向异性。在这样的地层中，平行于岩层平面的方向的水平电导率σh(或电阻率，Rh)不同于垂直于岩层平面的方向的垂向电导率σv(或电阻率，Rv)。
为了测量各向异性地层的电导率，新的EM感应或传播测井仪通常包括设置有发射器和/或接收器天线的横向阵列，以便它们的磁矩基本上垂直于仪器的轴线。比如参加Moran和Gianzero的“Effect ofFormation Anisotropy on resistivity Logging Measurements”，Geophysics，44，1266-1286(1979)。横向阵列测井仪包括具有三轴阵列装置。这些装置内的每一个三轴阵列包括三个正交发射器线圈和同样正交方向的三个接收器线圈。在作业中，三轴发射器在三个正交方向内被激活。之后，在同样的三个正交方向排列的每个接收器线圈测量由周围地层内流动的涡流感应出的电压。三轴测井仪的示例可在授权给Huston的美国专利US 3,510,757、授权给Forgang等人的美国专利US 5,781,436、授权给Desbrandes的美国专利US 3,609,521、授权给Segesman的美国专利US 4,360,777、授权给Kriegshuser等人的美国专利US 6,553,314中发现。
与常规的感应出在垂直于仪器纵向轴线的平面内流动的涡流的感应测井仪形成对比的是，横向阵列具有在垂直于仪器纵向轴线的平面内流动的涡流。三轴阵列具有发射器和接收器，每个发射器和接收器都具有三个以正交方向设置的线圈。因此，在三轴阵列内的发射器与接收器之间具有9个耦合。每一耦合对涡流流动的不同方向敏感。此外，每一EM感应或传播测井仪通常包括多个阵列。因而，对具有横向或三轴阵列的EM测井仪的探头误差校正比对常规的感应测井仪的探头误差校正更加复杂，并且需要具有对具有横向或三轴阵列的EM测井仪的探头误差进行校正的方法。
发明内容
一方面，本发明的实施例涉及用于确定具有横向阵列的测井仪内的探头误差的方法。根据本发明实施例的用于确定具有横向阵列的测井仪内的探头误差的方法包括：当测井仪处于以所选择的方向距地面第一高度时，使用横向阵列获得第一测量值；当测井仪处于以所选择的方向距地面第二高度时，使用横向阵列获得第二测量值；由第一测量值和第二测量值获得差值测量值；利用差值测量值和预设函数确定背景信号，其中所述的预设函数在一组地层电导率处将在第二高度处所得的一组信号内的地信号与在第一高度和第二高度处所得的信号之间的一组差值信号联系起来；以及通过从第二测量值中减去背景信号来确定探头误差。
一方面，本发明的实施例涉及用于确定具有三轴阵列的测井仪内的探头误差的方法。根据本发明实施例的用于确定具有三轴阵列的测井仪内的探头误差的方法包括：当测井仪处于以所选择的方向距地面第一高度时，获得用于三轴阵列的至少一个耦合的第一测量值；当测井仪处于以所选择的方向距地面第二高度时，获得用于三轴阵列的至少一个耦合的第二测量值；由用于三轴阵列的至少一个耦合的第一测量值和第二测量值获得差值测量值；利用差值测量值和用于三轴阵列的至少一个耦合的预设函数确定背景信号，其中所述的预设函数在一组地层电导率处将在第二高度处所得的一组信号内的地信号与在第一高度和第二高度处所得的信号之间的一组差值信号联系起来；以及通过从用于三轴阵列的至少一个耦合的第二测量值中减去背景信号来确定探头误差。
一方面，本发明的实施例涉及用于确定具有横向阵列的测井仪内的探头误差的系统。根据本发明实施例的用于确定具有横向阵列的测井仪内的探头误差的系统，该系统包括存储具有指令的程序的存储器，其用于：当测井仪处于以所选择的方向距地面第一高度时，使用横向阵列获得第一测量值；当测井仪处于以所选择的方向距地面第二高度时，使用横向阵列获得第二测量值；由第一测量值和第二测量值获得差值测量值；利用差值测量值和预设函数确定背景信号，其中所述的预设函数在一组地层电导率处将在第二高度处所得的一组信号内的地信号与在第一高度和第二高度处所得的信号之间的一组差值信号联系起来；以及通过从第二测量值中减去背景信号来确定探头误差。
一方面，本发明的实施例涉及用于确定具有三轴阵列的测井仪内的探头误差的系统。根据本发明实施例的用于确定具有三轴阵列的测井仪内的探头误差的系统，该系统包括存储具有指令程序的存储器，其用于：当测井仪处于以所选择的方向距地面第一高度时，获得用于三轴阵列的至少一个耦合的第一测量值；当测井仪处于以所选择的方向距地面第二高度时，获得用于三轴阵列的至少一个耦合的第二测量值；由用于三轴阵列的至少一个耦合的第一测量值和第二测量值获得差值测量值；利用差值测量值和用于三轴阵列的至少一个耦合的预设函数确定背景信号，其中所述的预设函数在一组地层电导率处将在第二高度处所得的一组信号内的地信号与在第一高度和第二高度处所得的信号之间的一组差值信号联系起来；以及通过从用于三轴阵列的至少一个耦合的第二测量值中减去背景信号来确定探头误差。
从下面的描述和所附的技术方案中本发明的其它方面和优点将变得很明显。
附图说明
图1A和1B分别示出了现有技术的水平和垂向探头误差校正的装置。
图2示出了用于常规感应测井仪的作为处于两个高度处的差值测量值函数的地信号的交会图。
图3示出了三轴阵列的示意图，其中示出了每一线圈的磁矩。
图4示出了根据本发明的实施例的探头误差校正的装置。
图5示出了根据本发明的实施例的探头误差校正的方法。
图6示出了根据本发明的实施例作为处于两个高度处的差值测量值函数的三轴阵列的zz耦合内的地信号的交会图。
图7示出了根据本发明的实施例作为处于两个高度处的差值测量值函数的三轴阵列的zx耦合内的地信号的交会图。
图8示出了根据本发明的实施例作为处于两个高度处的差值测量值函数的三轴阵列的yy耦合内的地信号的交会图。
图9示出了根据本发明的实施例作为处于两个高度处的差值测量值函数的三轴阵列的xx耦合内的地信号的交会图。
图10示出了根据本发明的实施例作为处于两个高度处的差值测量值函数的三轴阵列的xz耦合内的地信号的交会图。
图11示出了现有技术的计算机工作站的典型示图。
具体实施方式
本发明的实施例涉及测定具有横向或三轴阵列的感应或传播测井仪内地信号的方法，以便如果测井仪设置在自由空间内，测井仪内的探头误差(或自信号)能够被精确测定。该方法能够从距地面两个距离处测量横向或者三轴阵列内的所要耦合产生的信号并使用算法使两个距离处的测量值与地信号相关。
图1A示出了水平或“平行”探头误差装置10，图1B示出了“垂直”探头误差装置，如Barber专利所公开的。相同或相似的装置也可以与本发明的实施例一起使用。本领域的普通技术人员可以理解，虽然“垂直”或“水平”的装置是优选的实施例，但是其它装置也可以用于本发明的实施例。如图1A所示，探头误差装置10可以相对地面30基本水平放置。探头误差装置10优选具有最小量的金属或传导性材料。例如，其可以由木质杆和木销钉构成。使用绳索轮系统26将测井仪20、测井探头22和盒24带至第一高度h1和随后的高度h2，反之亦然。为了示出简便，盒24包括常用的感应或传播测井仪的电子盒和遥测电子线路短节。移动测试和控制装置40设置于探头误差装置10的附近。测试装置42包括处理器44、电源46，记录器48为测井领域的标准设备。测井电缆52通过测井绞车50将感应或传播测井仪20连接到测试装置40。
根据Barber专利所公开的方法，在两个高度处测量探头响应。两个高度处探头读数差为地信号的函数。图2示出了用于常规的感应或传播测井仪的两个高度处的电压差相对于地信号的校正曲线图。曲线为作为上水平(12英尺)和下水平(4英尺)处的测量值(δσ)之差的函数在上水平(12英尺)处的交会图。通过使用这张图和两个高度处的测量值之间的差值，可以测定上部位置处的地信号(σe)。之后，将地信号(σe)从上部位置(如12英尺)的探头读数中减去，从而得到固有的探头读数(探头误差)。然后，探头误差用于校正测井测量值。
虽然上述的方法已经用于常规的感应测井仪具有一段时间了，具有横向阵列(如三轴阵列)的测井仪建立了不同的涡流模式。例如。横向发射器感应出在平行于测井仪纵向轴线的平面内流动的涡流。相比之下，常规的纵向磁偶极子天线在垂直于测井仪纵向轴线的环圈内感应出涡流。三轴阵列包括发射器和接收器，发射器和接收器具有三个通常在正交方向的线圈。也就是，三轴阵列除具有标准的z轴线圈之外，还具有定向于正交方向的附加配置线圈(也就是横向线圈)。
图3示出了典型三轴阵列的示意图。发射器天线内的三个线圈具有排列在x，y和z向的磁矩，即MT x，MT y和MT z。接收器具有三个排列在相同正交方向的线圈，从而它们具有磁距MR x，MR y和MR z。在作业中，每一发射器线圈可被激活并且所述的信号由每一接收器线圈探测到。因此，对于每一三轴阵列具有9个可能的测量值。EM感应或传播测井仪通常包括多个三轴阵列。因而，对具有多个三轴阵列的EM感应或传播测井仪的探头误差校正比常规的感应测井仪明显地要复杂。
在三轴阵列中，所要测量的9个分量(耦合)可表示为：
xx yx zx xy yy zy xz yz zz - - - ( 1 )
这一矩阵内的每一元素表示从发射器a到接收器b的耦合(ab)，比如yx指从y发射器到x接收器的耦合。
为了校正三轴阵列的探头误差，发射器和接收器线圈的方向应该被定义。与探头误差校正装置的在图1A中示出，校正三轴阵列的装置在图4中示出。如图所示，三轴仪器设置在具有指向朝下的x向和平行于地面的y向的提升站(lift station)上，由于这一单一位置，所有的探头误差可以在不被提升站上的旋转测井仪干扰的情况下被计算出。利用图4所示的装置，xy，yx，yz和zy耦合信号应该为零，并且对这些12英尺处的耦合的测井仪的读数能够用作近似的探头误差以从各自的测井测量值中减去。
本领域的技术人员可以理解，图4所示的装置仅仅用于说明，并且其它变化可在不脱离本发明范围的情况下进行使用。例如，两个位置的高度可不同于图示的4英尺和12英尺。此外，y轴可指向地面并且x轴可平行于地面。利用这一可替换的装置，xy，yx，xz和zx耦合信号应该为零，12英尺处的测井仪读数为要从相应的测井测量值减去的近似的探头误差。在每一高度处完成两个测量值是可能的，一个是朝下的x轴的测量值和一个朝下的y轴测量值。在这一情况下，近似的探头误差可由xy，yx，yz，zy，xz和zx耦合的探头读数获得。因而，仅仅三个耦合(xx，yy和zz)需要进行校正。此外，虽然不是优选的，具有偏离向下和平行于地面位置的x轴和y轴是可能的。在这一情况下，xy，yx，yz，zy，xz和zx耦合信号应该不为零，并且用于这些耦合的探头误差就像其它耦合一样需要进行校正。
虽然图4所示的装置具有位于平行于地面方向(即水平模式)的仪器，但是仍需利用具有图1B所示的悬于两个高度处的垂向模式的仪器进行校正。在垂向模式中，z轴垂直于地面，而在垂向模式中x轴和y轴平行于地面。因而在垂向装置中，对于xy和yx的耦合应该为零。
图5示出了用于对三轴阵列的任何耦合进行探头误差校正的方法500的流程图。该方法500示出了测定每一耦合，即xx，xy和yz等的真实探头误差的步骤。首先，对上部和下部的专用仪器进行信号(测量值)的仿真。所述的模式考虑了仪器的结构(如阵列间距)、作业频率、装置的结构，以及其它因素。例如，适于这种仿真的合适的编码在Anderson等发表的在Transaction of the SPWLA 36th AnnualLogging Symposium，Paris，France，June26-29，1995，paper D上公开的文献“The Response of Induction Tools to DippingAnisotropic Formations，”中所公开。这篇文献所公开的编码能够模拟任何线圈类型的测井仪在一维水平层TI介质中的响应，所述的线圈类型的测井仪包括感应测井仪和传播测井仪。张量格林函数以谱积分的形式被表示出来并被分解成径向方位角波和轴向平面波。这一分解将三维积分降低为一维积分，其中，径向方位角方法对整个系统是共有的。在接触面上，轴向平面波的传播和反射用于描述分层介质中的EM波。边界层处的阻抗边界情况用于确定递归函数的反射系数和传播系数。超过径向波数目的一维积分要精心选择以确保精度。在一定的地层传导率范围内模拟测量值以提供一系列的测量值(步骤502)。
信号可独立地模拟实数(R)分量和正交(X)分量。就X分量来说，由发射器产生的不平衡的互感(mutual)也是产生探头误差的主要因素。之后，R和X信号可独立应用于下述步骤中。
不像常规的感应测井仪，其测量值通常代表传导率，三轴测井仪的原始数据以阻抗为单位(也就是，由于大地耦合而产生的发射器和计数器之间的阻抗)。阻抗单位保持了简单的旋转矩阵所必要的张量特征以用于坐标变换(从仪器坐标系统到大地坐标系统，反之亦然)。因此，优选计算出探头误差，并且探头误差拟合为阻抗。为了生成对于浮点计算合理的单位，阻抗单位优选为微欧姆。
接下来，对于每一地层电导率，通过将上部信号(U)从相应的下部信号(L)中减去而计算出差值测量值(8)(步骤504)，即δ＝L-U。
需要指出的是，上部信号(U)和下部信号(L)(测量值)各自包括来自真实探头误差(TSE)的分量(信号)和来自大地电导率(σe)的分量(信号)。例如，上部测量值(U)和下部测量值(L)可表示为：
U＝TSE+ZU    (2)
以及L＝TSE+ZL    (3)
其中，ZU分别表示上部和下部测量值内的大地电导率产生的信号。因此，差值(δ)可表示为：
δ＝ZL-ZU＝f’(σe)    (4)
等式(4)示出了差值信号(δ)为大地电导率的函数，因为大地背景信号ZU和ZL都是大地电导率(σe)的函数。
如等式(2)所示，上部测量值(U)包括来自大地电导率的信号(ZU)。因此，上部测量值(U)也是大地电导率(σe)的函数。由于差值测量值(δ)和上部测量值(U)都取决于大地电导率(σe)，所以它们之间存在一种关系。因而，接下来的步骤得出这一相关函数(步骤506)。用于推导出这一相关函数的方法在下面示出。
差值测量值(δ)和上部测量值(U)之间的关系较好地限定在这些测量值的自然对数内，因为接收器所接收的信号依赖于电导率的指数函数。例如，这一关系可定义为ln[abs(U)＝f{ln[abs(δ)]}，其中abs(δ)和abs(U)分别表示差值信号(δ)和上部测量值(U)的绝对值。使用绝对值是因为大多数的耦合产生负的阻抗值。确定函数f的精确形式是不必要的。然而，函数f可近似于多项式函数P，除了xz耦合是正数之外，对于所有的耦合，其可定义为：
P = - Σ n = 1 N a n [ ln ( - δ ) ] n - - - ( 5 )
对于xz耦合，
P = - Σ n = 1 N a n [ ln ( - δ ) ] n - - - ( 6 )
多项式系数a0，a1，...an取决于仪器结构的特征(比如阵列间距)、计算出的阻抗分量(R或者X信号)、以及测量值频率。这些系数可由[lnabs(δ)]到[lnabs(U)]的最小二乘方拟合来进行限定。虽然上部利用上部测量值(U)作为示例，但是本领域的普通技术人员可以理解该处理也可利用下部测量值(L)。
一旦定义了相关函数，也就是相关系数由最小二乘方拟合进行了定义，每一大地电导率(σe)处的上部信号测量值内的大地背景信号分量就能够被计算出(步骤508)。上部测量值ZU内的背景信号可定义为eP，即ZU＝eP，其中P为等式(5)或(6)中所示的多项式函数，并且系数可按照步骤506所述的最小二乘方拟合进行限定。作为地层电导率的函数或作为两个高度处的差值测量值的函数的背景信号可存储起来以在后面使用，要么作为查找表要么表示为后面要进行描述的图6所示的交会图。
一旦限定了差值测量值与上部信号(或者上部信号内的背景信号)之间的相关函数，就能够确定专用仪器内的真实探头误差。首先，使用同一校正装置内的仪器获得一个在第一高度处和一个在第二高度处的两个测量值(步骤501)。然后，由两个测量值获得差值测量值(步骤503)。之后，这一差值测量值与步骤506所获得的相关函数或者步骤508所获得的背景信号表或曲线一起使用以得到上部测量值内的背景信号(步骤505)。接着，将估算出的上部测量值内的背景信号从上部测量值内减去从而得到真实的探头误差(TSE)(步骤507)。
前述的将差值信号与上部测量值联系起来的过程和多项式函数的拟合最好用交会图表示。图6示出了大地(背景)信号的交会图60，所述的地信号为使用电磁传播测井仪的39英寸的阵列的4英尺与12英尺zz测量值之间差值(δ)的函数。所述的测井仪在图4中示出，并且R信号在从1mS/m至2000mS/m之间的各个电导率处以26.8074KHz的频率进行测量。在这张图中，x轴表示差值测量值(δ)，其为大地电导率(σe)的函数，y轴表示计算出的上部测量值内的背景信号。曲线60基于由差值测量值多项式拟合到上部测量值所推导出的系数。
利用这张曲线图和两个高度处的zz测量值之间的差值(δ)，由在上部位置处的大地电导率(Z12)产生的信号分量可以被确定。例如，线61可在两个高度处的测量值之间的差值处画出。线61在点62处与交会图曲线60相交。之后，水平线63从点62引出以获得由大地电导率产生的背景信号的值。然后，将背景信号从上部位置(12英尺)处的探头读数中减去以得到真实探头误差。
虽然上面的描述利用上部测量值来获得相关函数并确定真实探头误差，本领域的普通技术人员可以理解，利用下部测量值也可进行同样的步骤。因此，本发明的实施例并不局限于使用上部测量值。校正步骤可与垂向模式或水平模式的测井仪一起使用。此外，本发明的实施例也可用于感应测井仪或者传播测井仪。
应该指出的是，真实探头误差和用于上述校正函数的系数取决于阵列结构(如阵列间距和线圈的方向)。因此，对于每一耦合(xx，yy，zz，xz，zx的等)真实探头误差应该单独计算。图7-10示出了类似于图6所示的其它耦合的交会图。此外，由于探头误差依赖于作业频率，这些校正应该以测井作业应用的每一频率来进行。
本发明的一些实施例涉及用于完成上述方法的系统。根据本发明实施例的系统可以为独立的单位以进行本发明的方法或者可安装在钻具中(如图1A所示的电子盒24或者测试装置42)。根据本发明的系统通常包括处理器和存储器。在一些实施例中，系统可具有含有处理器、存储器的一般用途的计算机，并可优选包括其它硬件。例如，如图11所示，一般的计算机(150)包括处理器(152)、随机存取存储器(154)和存储装置(如永久存储器或硬盘)(156)。计算机(150)还可包括输入设备，如键盘(158)和鼠标(160)，和输出设备，如监视器(162)。应该指出的是，一般用途的计算机仅用于示例性的示出，本发明的实施例可采用其它形式(如结合在测井仪内)。
根据本发明的系统，存储器存储了可由处理器读取的程序。比如，该程序可包括完成上述方法的指令：获得两个高度处的测量值；由两个高度处的测量值得出差值测量值；利用将上部测量值与地层电导率(或者一组地层电导率处的差值测量值)联系起来的预设函数由差值测量值得出背景信号；通过将背景信号从上部信号中减去而得到真实探头误差。预设函数(相关函数)可从前地层模型和从上部测量值到上述的差值测量值的运算拟合而得出。
程序设计可通过使用一个或多个可由计算机处理器读取的程序存储设备和编译一个或多个可由计算机执行的指令程序来完成，以进行这里描述的操作。例如，程序存储设备可采用如下形式：一个或多个软盘；CD-ROM或者其它光盘；磁带；只读存储器片(ROM)；以及其它本领域公知的形式。指令程序可为“目标码”，即可直接由计算机执行的二进制形式；可为在执行之前需要编译或数据整理分析的“源代码”；或者为如部分编译码的一些中间形式。程序存储设备和指令解码的精确形式在这里是不重要的。
本发明的优点包括下述的优点。用于探头误差校正的方法通常可应用于感应测井仪或传播测井仪。此外，这些方法可与单线圈阵列或者完整的三轴阵列一起使用。即使所述的阵列并不是正交阵列，即偏斜阵列，也可以使用该方法。本发明的实施例不需要旋转测井仪来对不同方向的线圈进行校正，仅仅需要两个或者多个高度处的测量值。由于不需要旋转测井仪，所以不会产生由不精确的旋转而产生的误差。一旦进行了校正，要么探头误差要么校正函数和/或系数可存储起来以在后面进行测井测量值内探头误差的校正。
虽然结合有限数量的实施例对本发明进行了描述，受益于这一公开的本领域的技术人员可以理解，可以在不脱离这里对本发明所公开的范围而设计出其它实施例。因此，本发明的范围应该仅仅由所附的技术方案进行限制。

Claims (19)

1.一种用于确定具有横向阵列的测井仪内的探头误差的方法，其包括：
当测井仪处于以所选择的方向距地面第一高度时，使用横向阵列获得第一测量值；
当测井仪处于以所选择的方向距地面第二高度时，使用横向阵列获得第二测量值；
由第一测量值和第二测量值获得差值测量值；
利用差值测量值和预设函数确定背景信号，其中所述的预设函数在一组地层电导率处将在第二高度处所得的一组信号内的地信号与在第一高度和第二高度处所得的信号之间的一组差值信号联系起来；以及
通过从第二测量值中减去背景信号来确定探头误差。
2.如权利要求1所述的方法，其中所述的第一高度小于第二高度。
3.如权利要求1所述的方法，其中所述的预设函数通过将该组差值信号拟合到在第二高度处所得的该组信号而获得。
4.如权利要求1所述的方法，其中所述在第二高度处所得的一组信号和一组差值信号由地层模型的仿真而获得。
5.如权利要求1所述的方法，其中所述的测井仪水平设置在地面上方。
6.一种用于确定具有横向阵列的测井仪内的探头误差的方法，其包括：
当测井仪处于以所选择的方向距地面第一高度时，获得用于三轴阵列的至少一个耦合的第一测量值；
当测井仪处于以所选择的方向距地面第二高度时，获得用于三轴阵列的至少一个耦合的第二测量值；
由用于三轴阵列的至少一个耦合的第一测量值和第二测量值获得差值测量值；
利用差值测量值和用于三轴阵列的至少一个耦合的预设函数确定背景信号，其中所述的预设函数在一组地层电导率处将在第二高度处所得的一组信号内的地信号与在第一高度和第二高度处所得的信号之间的一组差值信号联系起来；以及
通过从用于三轴阵列的至少一个耦合的第二测量值中减去背景信号来确定探头误差。
7.如权利要求6所述的方法，其中所述的第一高度小于第二高度。
8.如权利要求6所述的方法，其中所述的预设函数通过将该组差值信号拟合到在第二高度处所得的该组信号而获得。
9.如权利要求6所述的方法，其中所述在第二高度处所得的一组信号和一组差值信号由地层模型的仿真而获得。
10.如权利要求6所述的方法，其中所述的测井仪水平设置在地面上方。
11.一种用于确定具有横向阵列的测井仪内的探头误差的系统，该系统包括存储具有指令的程序的存储器，其用于：
当测井仪处于以所选择的方向距地面第一高度时，使用横向阵列获得第一测量值；
当测井仪处于以所选择的方向距地面第二高度时，使用横向阵列获得第二测量值；
由第一测量值和第二测量值获得差值测量值；
利用差值测量值和预设函数确定背景信号，其中所述的预设函数在一组地层电导率处将在第二高度处所得的一组信号内的地信号与在第一高度和第二高度处所得的信号之间的一组差值信号联系起来；以及
通过从第二测量值中减去背景信号来确定探头误差。
12.如权利要求11所述的系统，其中所述的第一高度小于第二高度。
13.如权利要求11所述的系统，其中所述的预设函数通过将该组差值信号拟合到在第二高度处所得的该组信号而获得。
14.如权利要求11所述的系统，其中所述在第二高度处所得的一组信号和一组差值信号由地层模型的仿真而获得。
15.一种用于确定具有三轴阵列的测井仪内的探头误差的系统，该系统包括存储具有指令的程序的存储器，其用于：
当测井仪处于以所选择的方向距地面第一高度时，获得用于三轴阵列的至少一个耦合的第一测量值；
当测井仪处于以所选择的方向距地面第二高度时，获得用于三轴阵列的至少一个耦合的第二测量值；
由用于三轴阵列的至少一个耦合的第一测量值和第二测量值获得差值测量值；
利用差值测量值和用于三轴阵列的至少一个耦合的预设函数确定背景信号，其中所述的预设函数在一组地层电导率处将在第二高度处所得的一组信号内的地信号与在第一高度和第二高度处所得的信号之间的一组差值信号联系起来；以及
通过从用于三轴阵列的至少一个耦合的第二测量值中减去背景信号来确定探头误差。
16.如权利要求15所述的系统，其中所述的第一高度小于第二高度。
17.如权利要求15所述的系统，其中所述的预设函数通过将该组差值信号拟合到在第二高度处所得的该组信号而获得。
18.如权利要求15所述的系统，其中所述在第二高度处所得的一组信号和一组差值信号由地层模型的仿真而获得。
19.如权利要求15所述的系统，其中所述的测井仪水平设置在地面上方。
CNB2004101007593A 2003-12-12 2004-12-13 确定横向或三轴阵列感应或传播测井仪的探头误差的方法 Expired - Fee Related CN100374879C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US10/707427 2003-12-12
US10/707,427 US7027923B2 (en) 2003-12-12 2003-12-12 Method for determining sonde error for an induction or propagation tool with transverse or triaxial arrays

Publications (2)

Publication Number Publication Date
CN1648691A CN1648691A (zh) 2005-08-03
CN100374879C true CN100374879C (zh) 2008-03-12

Family

ID=33565366

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2004101007593A Expired - Fee Related CN100374879C (zh) 2003-12-12 2004-12-13 确定横向或三轴阵列感应或传播测井仪的探头误差的方法

Country Status (7)

Country Link
US (1) US7027923B2 (zh)
CN (1) CN100374879C (zh)
CA (1) CA2488144C (zh)
GB (1) GB2409044B (zh)
MX (1) MXPA04012088A (zh)
NO (1) NO335636B1 (zh)
RU (1) RU2401442C2 (zh)

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7375530B2 (en) * 2002-03-04 2008-05-20 Baker Hughes Incorporated Method for signal enhancement in azimuthal propagation resistivity while drilling
US7932723B2 (en) * 2004-05-07 2011-04-26 Baker Hughes Incorporated Borehole conductivity simulator verification and transverse coil balancing
US7319331B2 (en) * 2004-05-07 2008-01-15 Baker Hughes Incorporated Two loop calibrator
US7205770B2 (en) * 2004-05-07 2007-04-17 Baker Hughes Incorporated Borehole conductivity simulator verification and transverse coil balancing
US7652478B2 (en) * 2004-05-07 2010-01-26 Baker Hughes Incorporated Cross-component alignment measurement and calibration
US7408355B1 (en) * 2004-05-07 2008-08-05 Baker Hughes Incorporated Borehole conductivity simulator verification and transverse coil balancing
US7969153B2 (en) * 2004-05-07 2011-06-28 Baker Hughes Incorporated Borehole conductivity simulator verification and transverse antenna balancing
US8736270B2 (en) 2004-07-14 2014-05-27 Schlumberger Technology Corporation Look ahead logging system
US20060208737A1 (en) * 2005-03-16 2006-09-21 Baker Hughes Incorporated Calibration of xx, yy and zz induction tool measurements
EA012740B1 (ru) * 2005-08-03 2009-12-30 Шелл Интернэшнл Рисерч Маатсхаппий Б.В. Способ и система для определения электромагнитного отклика из формации земной коры, и способ бурения буровой скважины, и способ добычи углеводородного флюида
WO2007089438A2 (en) * 2006-01-26 2007-08-09 Baker Hugues Incorporated Borehole conductivity simulator verification and transverse coil balancing
US7268555B1 (en) 2006-04-06 2007-09-11 Baker Hughes Incorporated Correction of cross-component induction measurements for misalignment using comparison of the XY formation response
US8931335B2 (en) * 2006-04-07 2015-01-13 Baker Hughes Incorporation Processing of multi-component induction data in the presence of borehole abnormalities
US7545145B2 (en) * 2006-05-03 2009-06-09 Baker Hughes Incorporated Method and apparatus for tensorial micro-resistivity imaging in oil-based muds
US7457707B2 (en) * 2006-11-30 2008-11-25 Schlumberger Technology Corporation Method for determining an effective formation conductivity for induction log borehole correction
US20080314582A1 (en) * 2007-06-21 2008-12-25 Schlumberger Technology Corporation Targeted measurements for formation evaluation and reservoir characterization
US7915895B2 (en) * 2007-06-22 2011-03-29 Baker Hughes Incorporated Method of calibrating an azimuthal inductive cross-coil or tilted coil instrument
AU2008272905B2 (en) * 2007-07-03 2011-05-19 Shell Internationale Research Maatschappij B.V. System and method for measuring a time-varying magnetic field and method for production of a hydrocarbon fluid
CN101932955B (zh) * 2007-08-27 2014-12-03 普拉德研究及开发股份有限公司 先行测井系统
US9250352B2 (en) 2008-04-17 2016-02-02 Richard H. Hardman Methods for producing a log of material properties
US8441261B2 (en) * 2010-06-16 2013-05-14 Schlumberger Technology Corporation Determination of conductive formation orientation by making wellbore sonde error correction
US8762107B2 (en) * 2010-09-27 2014-06-24 Baker Hughes Incorporated Triaxial induction calibration without prior knowledge of the calibration area's ground conductivity
US10371852B2 (en) 2011-12-21 2019-08-06 Schlumberger Technology Corporation Formation properties from conductivity tensor
EP2815070B1 (en) 2012-02-17 2016-08-24 Services Pétroliers Schlumberger Inversion-based calibration of downhole electromagnetic tools
RU2641628C2 (ru) * 2013-07-26 2018-01-18 Хэллибертон Энерджи Сервисиз Инк. Компьютерная программа для калибровки прибора для проведения каротажа сопротивления в скважине
CA2920169A1 (en) * 2013-09-10 2015-03-19 Halliburton Energy Services, Inc Surface calibration of a wellbore resistivity logging tool
CN106837299B (zh) * 2017-01-19 2019-11-19 中国科学院声学研究所 一种井眼校正的系统及方法
CA3064876C (en) 2017-08-07 2021-12-28 Halliburton Energy Services, Inc. Component-based look-up table calibration for modularized resistivity tool
EP3492952B1 (en) 2017-12-01 2022-01-26 Services Pétroliers Schlumberger Calibration of electromagnetic measurement tool
US11782180B2 (en) 2021-09-02 2023-10-10 Halliburton Energy Services, Inc. Mixed air-hang response for correcting a formation signal of a mixed set measuring device
CN115016034B (zh) * 2022-06-01 2022-11-25 中国科学院地质与地球物理研究所 一种随钻测量装置的标定方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4800496A (en) * 1987-09-28 1989-01-24 Schlumberger Technology Corporation Method for determining induction sonde error
US6594584B1 (en) * 1999-10-21 2003-07-15 Schlumberger Technology Corporation Method for calculating a distance between a well logging instrument and a formation boundary by inversion processing measurements from the logging instrument
US20030155923A1 (en) * 1999-12-01 2003-08-21 Schlumberger Technology Corporation Shielding method and apparatus using flexible strip
US20030155924A1 (en) * 2000-12-15 2003-08-21 Schlumberger Technology Corporation Borehole current correction methods and apparatus
US20030184488A1 (en) * 2002-03-29 2003-10-02 Smith David L. Simplified antenna structures for logging tools

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3510757A (en) 1966-09-01 1970-05-05 Schlumberger Technology Corp Formation dip measuring methods and apparatus using induction coils
FR1578713A (zh) 1967-07-12 1969-08-22
US4360777A (en) 1979-12-31 1982-11-23 Schlumberger Technology Corporation Induction dipmeter apparatus and method
DE4123760C2 (de) * 1991-07-18 2000-01-20 Dade Behring Marburg Gmbh Seroreaktive Bereiche auf den HPV 16 Proteinen E1 und E2
US5757997A (en) * 1995-12-22 1998-05-26 Minnesota Mining And Manufacturing Company Optical fiber connector using fiber spring force alignment groove
US5781436A (en) 1996-07-26 1998-07-14 Western Atlas International, Inc. Method and apparatus for transverse electromagnetic induction well logging
US6068957A (en) 1998-06-29 2000-05-30 Eastman Kodak Company Lubricating layer in photographic elements
US6304086B1 (en) 1999-09-07 2001-10-16 Schlumberger Technology Corporation Method and apparatus for evaluating the resistivity of formations with high dip angles or high-contrast thin layers
US6466872B1 (en) 1999-11-08 2002-10-15 Baker Hughes Incorporated Method for determination of apparent resistivities of anisotropic reservoirs
US6297639B1 (en) 1999-12-01 2001-10-02 Schlumberger Technology Corporation Method and apparatus for directional well logging with a shield having sloped slots
US6351127B1 (en) 1999-12-01 2002-02-26 Schlumberger Technology Corporation Shielding method and apparatus for selective attenuation of an electromagnetic energy field component
US6727705B2 (en) 2000-03-27 2004-04-27 Schlumberger Technology Corporation Subsurface monitoring and borehole placement using a modified tubular equipped with tilted or transverse magnetic dipoles
US6836218B2 (en) 2000-05-22 2004-12-28 Schlumberger Technology Corporation Modified tubular equipped with a tilted or transverse magnetic dipole for downhole logging
US6541979B2 (en) 2000-12-19 2003-04-01 Schlumberger Technology Corporation Multi-coil electromagnetic focusing methods and apparatus to reduce borehole eccentricity effects
US6584408B2 (en) 2001-06-26 2003-06-24 Schlumberger Technology Corporation Subsurface formation parameters from tri-axial measurements
US6556016B2 (en) * 2001-08-10 2003-04-29 Halliburton Energy Services, Inc. Induction method for determining dip angle in subterranean earth formations
US6969994B2 (en) 2001-09-26 2005-11-29 Schlumberger Technology Corporation Directional electromagnetic measurements insensitive to dip and anisotropy
US6556015B1 (en) 2001-10-11 2003-04-29 Schlumberger Technology Corporation Method and system for determining formation anisotropic resistivity with reduced borehole effects from tilted or transverse magnetic dipoles
US6690170B2 (en) 2002-03-29 2004-02-10 Schlumberger Technology Corporation Antenna structures for electromagnetic well logging tools
US6667620B2 (en) 2002-03-29 2003-12-23 Schlumberger Technology Corporation Current-directing shield apparatus for use with transverse magnetic dipole antennas
US6998844B2 (en) 2002-04-19 2006-02-14 Schlumberger Technology Corporation Propagation based electromagnetic measurement of anisotropy using transverse or tilted magnetic dipoles
US6794875B2 (en) 2002-05-20 2004-09-21 Halliburton Energy Services, Inc. Induction well logging apparatus and method
US7414391B2 (en) * 2002-07-30 2008-08-19 Schlumberger Technology Corporation Electromagnetic logging tool calibration system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4800496A (en) * 1987-09-28 1989-01-24 Schlumberger Technology Corporation Method for determining induction sonde error
US6594584B1 (en) * 1999-10-21 2003-07-15 Schlumberger Technology Corporation Method for calculating a distance between a well logging instrument and a formation boundary by inversion processing measurements from the logging instrument
US20030155923A1 (en) * 1999-12-01 2003-08-21 Schlumberger Technology Corporation Shielding method and apparatus using flexible strip
US20030155924A1 (en) * 2000-12-15 2003-08-21 Schlumberger Technology Corporation Borehole current correction methods and apparatus
US20030184488A1 (en) * 2002-03-29 2003-10-02 Smith David L. Simplified antenna structures for logging tools

Also Published As

Publication number Publication date
CN1648691A (zh) 2005-08-03
CA2488144C (en) 2008-08-19
GB2409044A (en) 2005-06-15
NO335636B1 (no) 2015-01-12
MXPA04012088A (es) 2005-09-21
US7027923B2 (en) 2006-04-11
GB0426185D0 (en) 2004-12-29
RU2401442C2 (ru) 2010-10-10
RU2004136270A (ru) 2006-05-20
GB2409044B (en) 2006-04-19
CA2488144A1 (en) 2005-06-12
NO20044894L (no) 2005-06-13
US20050143920A1 (en) 2005-06-30

Similar Documents

Publication Publication Date Title
CN100374879C (zh) 确定横向或三轴阵列感应或传播测井仪的探头误差的方法
RU2378669C2 (ru) Способ внесения поправок на влияние скважины в данные трехмерных индукционных зондов
US9547102B2 (en) Resistivity logging systems and methods employing ratio signal set for inversion
US10358911B2 (en) Tilted antenna logging systems and methods yielding robust measurement signals
US8841913B2 (en) Determining formation parameters using electromagnetic coupling components
CA2010398C (en) Method to improve directional survey accuracy
CN104169524B (zh) 井下电磁工具的基于反演的校正
US8762107B2 (en) Triaxial induction calibration without prior knowledge of the calibration area's ground conductivity
AU2011380527B2 (en) Multi-component induction logging systems and methods using real-time OBM borehole correction
US20070236221A1 (en) Method and Apparatus for the Use of Multicomponent Induction Tool for Geosteering and Formation Resistivity Data Interpretation in Horizontal Wells
US20110254552A1 (en) Method and apparatus for determining geological structural dip using multiaxial induction measurements
US20150177412A1 (en) Calibration method for multi-component induction tools
US7952358B2 (en) Methods and apparatus for determination of a formation resistivity property using eccentered measurements
CN111502648A (zh) 一种随钻远探测电磁波电阻率测井仪器刻度方法与装置
US20050083061A1 (en) Methods and systems for estimating formation resistivity that are less sensitive to skin effects, shoulder-bed effects and formation dips
US10365395B2 (en) Multi-component induction logging systems and methods using blended-model inversion
CA2432684A1 (en) Processing well logging data with neural network
CA2447179C (en) Method for asymptotic dipping correction
US5987386A (en) Interpolation of induction tool response using geometrical factors as basis functions
Zhou et al. A calibration method of ground-air frequency domain electromagnetic receving system based on magnetic source excitation
Zhu et al. Forward Modeling of the Triaxial Induction Logging Instrument in Layered Anisotropic Medium
Anderson et al. Synthetic Deep Propagation Tool Response by the Finite Element Method
Center Triaxial Induction Logging

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20080312

Termination date: 20171213

CF01 Termination of patent right due to non-payment of annual fee