CN1185508C - 无趋肤效应电磁波测井方法 - Google Patents

无趋肤效应电磁波测井方法 Download PDF

Info

Publication number
CN1185508C
CN1185508C CNB00109775XA CN00109775A CN1185508C CN 1185508 C CN1185508 C CN 1185508C CN B00109775X A CNB00109775X A CN B00109775XA CN 00109775 A CN00109775 A CN 00109775A CN 1185508 C CN1185508 C CN 1185508C
Authority
CN
China
Prior art keywords
coil
conductivity
transmitting coil
omega
receiving coil
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
CNB00109775XA
Other languages
English (en)
Other versions
CN1332380A (zh
Inventor
陶果
刘国强
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.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CNB00109775XA priority Critical patent/CN1185508C/zh
Publication of CN1332380A publication Critical patent/CN1332380A/zh
Application granted granted Critical
Publication of CN1185508C publication Critical patent/CN1185508C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明涉及石油勘探中测量井周不同深度的电导率(电阻率)、介电常的方法,特点是用Lambert W函数求解其电磁场传播系数,进而提取地层电导率。当包含MAXT个发射线圈,MAXR个接收线圈时,用如下公式计算出地层视电导率,所得地层电导率信号不存在趋肤效应影响,测量结果准确,在发射频率较高时,可同时获得准确地层介电常数,可以把实现这种方法的电子线路装入井下,直接输出。

Description

无趋肤效应电磁波测井方法
所属技术领域
本发明涉及一种石油勘探领域中测量井周地层不同深度处的电导率(电阻率)、介电常数的无趋肤效应电磁波测井方法。
背景技术
石油勘探中现有的测量井周地层电导率的感应测井方法,按所测量提取的信号方式来分类,主要有以下两种:即测量R信号和同时测量R信号和X信号的感应测井方法。
前者将包含地层重要信息的X信号看成无用信号,采用硬件聚焦的方式压制X信号,然后据R信号来确定地层的电导率。当地层的电导率较大时,受趋肤效应的严重影响,将使地层电导率的测量值低于地层的真电导率,导致仪器响应具有很大误差。
后者同时测量R信号和X信号,利用X信号中包含大量的趋肤效应信息的性质,采用相量反褶积的方法对X信号的形状和幅度进行变换,进而校正实部信号中由于趋肤效应的影响而产生的误差。这样,在一定程度上改进了测量精度,但缺陷是:1.从R信号和X信号与地层电导率的近似关系出发,由于测量电导率与地层电导率成非线性关系必然导致地层电导率测量精度的下降;2.为弥补上述缺陷而采用的信号处理校正方法的合理性不能在物理机制上给出令人信服的解释;3.信号处理过程既要兼顾趋肤效应的校正和保证纵向分辨率不减小又要有益于实时处理,必须采用折衷方案,导致处理过程繁琐,带有人为的因素,其效果的好坏在很大程度上依赖于后验。
测量R信号方法的常规感应测井方法对单发单收仪器所得地层视电导率为:
σ R = Re ( V ) K = 2 L 2 ωμ { [ ( 1 + P ) sin P - P cos P ] + i [ ( 1 + P ) cos P + p sin P ] } e - P
L是发射线圈和接收线圈之间的距离;
NT和NR分别为发射线圈和接收线圈的匝数;
aT和aR分别为发射线圈和接收线圈的半径;
发射线圈的发射电流为I;
角频率为ω;地层电导率为σ;磁导率为μ
当线圈缩成点状时,发射线圈在接收线圈处的感应电动势V;
K为仪器常数: K = ω 2 μ 2 π α T α R 2 N T N R I 4 L
P为传播系数:P=L2ωμσ/2;
Re表示取实部。
由于上面的式子是在ωσ<<ε(σ为电导率,ε为介电常数)的假设条件下得到的。而且,视电导率只有当。小于0.01S·m-1时才接近于真电导率,在σ较大时,受趋肤效应的严重影响,将使地层电导率的测量值低于地层的真电导率,从而降低仪器测量精度。如果此条件不成立,则P值将改变,这时,σa将包含岩石电导率σ与介电常数ε的信息。这样导致所测视电导率与地层的真电导率相差更远。
同时测量R信号和X信号的相量双感应测井方法所得地层视电导率为:
σ P ( j ) = Σ L = 1 N h ( j - l ) σ R ( l ) + A ( σ x ) Σ l - 1 N b a ( j - l ) σ x ( l )
其中: σ R = Re ( V ) K ; σ X = Im ( v - V 0 ) K , Re、Im分别表示取实部、虚部,Vo为直耦信号。式中,σP是校正后的测量值,h是线性反褶积滤波因子,σR(1)是R信号的第1个子样,A是幅度拟合函数,ba是X信号拟合滤波因子,σx(1)是X信号的第1个子样。这里h,ba,A的确定复杂且带有人为因素。
综上所述,到目前为止,还没有能消除趋肤效应影响的准确测量不同深度处地层电导率(电阻率)的方法。
发明内容
本发明目的是提供一种不含趋肤效应影响的准确测量不同深度处地层电导率(电阻率)的无趋肤效应电磁波测井方法。
本发明目的是通过如下技术方案实现:
无趋肤效应电磁波测井方法对于任何线圈系组合的感应测井仪所取得的测井信号,用Lambert W函数求解其电磁场传播系数,进而提取地层电导率。
无趋肤效应电磁波测井方法线圈系包含MAXT个发射线圈,MAXR个接收线圈,分布的发射线圈发发出电磁信号,将接收线圈接收的电磁信号采用如下公式计算出地层视电导率σ:
σ mn = Im ( [ W ( V mn ω m μL mn 2 i 2 K mn e ) + 1 ] 2 ω m μ L mn 2 , ϵ mn = Re ( [ W ( V mn ω m μ L mn 2 i 2 K mn e ) + 1 ] 2 ω m μ L mn 2
式中:Lmn是第m个发射线圈和第n个接收线圈之间的距离; K m = ω m 2 μ 2 πa Tm 2 a Rn 2 N Tm N Rn I m 4 L mn 为常规感应测井的仪器常数;
NTm和NRn分别为第m个发射线圈和第n个接收线圈的匝数;
aTm和aRn分别为第m个发射线圈和第n个接收线圈的半径;
第m个发射线圈的发射电流为Im,角频率为ωm,磁导率为μ;
第m个发射线圈在第n个接收线圈处的感应电动势为Vmn
W为Lambert W函数,Re和Im分别表示取实部和虚部。
本发明主要特点是:
1.所得地层电导率信号不存在趋肤效应影响,测量结果准确,在发射频率较高时,可同时获得准确地层介电常数,可以把实现这种方法的电子线路装入井下,直接输出。
2.提取地层电导率(电阻率)时无需作趋肤效应校正。
3.在发射电流频率较高时,可同时输出地层介电常数曲线。
4.对发射线圈和接收线圈的位置可以不加限制。
5.没有对地层真电导率的大小和发射电流频率作任何限制,因而计算出的结果为地层电导率的准确值。
附图说明
本发明附图说明如下:
图1为本发明发射线圈和接收线圈的位置示意图;
图2为本发明效果对比图。
具体实施方式
以下结合附图说明实施例:
本发明所得地层视电导率为:
σ mn = Im ( [ W ( V mn ω m μ L mn 2 i 2 K mn e ) + 1 ] 2 ω m μL mn 2 , ϵ mn = Re ( [ W ( V mn ω m μL mn 2 i 2 K mn e ) + 1 ] 2 ω m μL mn 2
式中:W为Lambert W函数,Re和Im分别表示取实部和虚部,仪器包含MAXT个发射线圈,MAXR个接收线圈。发射线圈和接收线圈在仪器中的位置可以不加限制。Lmn是第m个发射线圈和第n个接收线圈之间的距离;NTm和NRn分别为第m个发射线圈和第n个接收线圈的匝数;aTm和aRn分别为第m个发射线圈和第n个接收线圈的半径;第m个发射线圈的发射电流为Im,角频率为ωm,由电磁感应的原理可知,在均匀介质空间中,若其电导率为σ,磁导率为μ,当线圈缩成点状时,设第m个发射线圈在第n个接收线圈处的感应电动势为Vmn
实施例1:单发单收感应仪器提取地层电导率的方法
设发射电流频率为200MHz,发射线圈和接收线圈的间距为L=0.8m。下面是多层模型提取电导率的方法比较,其中:L&T代表本发明所得电导率,MODEL表示模型电导率,R代表常规感应测井方法所得电导率,RX代表相量感应测井方法所得电导率。可以看出,本发明所提方法所得电导率最接近地层模型。
实施例2:多发多收感应仪器提取地层电导率的方法
采用多发多收感应线圈系组成,提取地层电导率的方法与实施例1类似,只是把多发多收线圈系分解为多个单发单收仪器,对所得信号进行组合,加权优化,得到具有合理纵向分辨率和探测深度的合成曲线。对多层理论模型来说,结果与上图相同。
实施例3:单发多收感应仪器提取地层复电导率的方法仪器由多个单发双收线圈系组成,同样引人Lambert W函数,与实施例2不同的是,这里利用双收线圈的幅度比和相位差来获得地层电导率和介电常数。再对所得信号进行组合,加权优化,得到具有合理纵向分辨率和探测深度的合成曲线。

Claims (1)

1.一种无趋肤效应电磁波测井方法,其特征在于:对于线圈系组合的感应测井仪所取得的测井信号,用Lambert W函数求解其电磁场传播系数,进而提取地层电导率,线圈系包含MAXT个发射线圈,MAXR个接收线圈,分布的发射线圈发发出电磁信号,将接收线圈接收的电磁信号采用如下公式计算出地层视电导率б:
σ mn = lm ( [ W ( V mn ω m μL mn 2 i 2 K mn e ) + 1 ] 2 ω m μL mn 2 , ϵ mn = Re ( [ W ( V mn ω m μL mn 2 i 2 K mn e ) + 1 ] 2 ω 2 μL mn 2
式中:
Lmn是第m个发射线圈和第n个接收线圈之间的距离;
K mn = ω m 2 μ 2 πa Tm 2 a Rn 2 N Tm N Rn I m 4 L mn 为常规感应测井的仪器常数;
NTm和NRn分别为第m个发射线圈和第n个接收线圈的匝数;
aTm和aRn分别为第m个发射线圈和第n个接收线圈的半径;
第m个发射线圈的发射电流为Im,角频率为ωm,磁导率为μ;
第m个发射线圈在第n个接收线圈处的感应电动势为Vmn
W为Lambert W函数,Re和Im分别表示取实部和虚部。
CNB00109775XA 2000-07-07 2000-07-07 无趋肤效应电磁波测井方法 Expired - Fee Related CN1185508C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB00109775XA CN1185508C (zh) 2000-07-07 2000-07-07 无趋肤效应电磁波测井方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB00109775XA CN1185508C (zh) 2000-07-07 2000-07-07 无趋肤效应电磁波测井方法

Publications (2)

Publication Number Publication Date
CN1332380A CN1332380A (zh) 2002-01-23
CN1185508C true CN1185508C (zh) 2005-01-19

Family

ID=4579857

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB00109775XA Expired - Fee Related CN1185508C (zh) 2000-07-07 2000-07-07 无趋肤效应电磁波测井方法

Country Status (1)

Country Link
CN (1) CN1185508C (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1252492C (zh) * 2003-12-25 2006-04-19 周仁安 大地电磁波电阻率测量方法及其仪器
US7394258B2 (en) * 2005-08-15 2008-07-01 Baker Hughes Incorporated High resolution resistivity earth imager
CN102678106B (zh) * 2012-05-02 2015-04-01 中国电子科技集团公司第二十二研究所 随钻电磁波电阻率测井仪器的数据处理方法
CN104453866B (zh) * 2013-09-18 2019-04-23 中国石油集团长城钻探工程有限公司 感应测井用电流测量的相位标定方法
KR101885666B1 (ko) * 2016-09-01 2018-08-06 (주) 멀티패스 Rf 신호를 이용한 비접촉 방식의 전도도 및 비전도체 유전율 특성변화 측정장치
CN111305813B (zh) * 2018-12-12 2023-06-27 天津大学青岛海洋技术研究院 一种基于套管井几何因子的电阻率处理方法
CN109779621B (zh) * 2019-01-30 2022-05-13 北京工业大学 一种感应测井仪器测井响应方法及装置

Also Published As

Publication number Publication date
CN1332380A (zh) 2002-01-23

Similar Documents

Publication Publication Date Title
CN104656157B (zh) 一种识别页岩气甜点区的方法及装置
US20030229450A1 (en) Induction well logging apparatus and method
CN100429529C (zh) 利用全波列、偶极横波测井资料确定气层的方法
Lamontagne Applications of wide band time domain EM measurements in mineral exploration
CN1222978A (zh) 分层地层的电测井
CN102678106A (zh) 随钻电磁波电阻率测井仪器的数据处理方法
CN109085653A (zh) 一种深部地质、硫化物矿体资源的探测方法
CN1185508C (zh) 无趋肤效应电磁波测井方法
CN102147482A (zh) 一种用于浅层csamt方法的感应式磁传感器的发明
WO2008031505A1 (en) An antenna for an electromagnetic probe for investigating geological formations and its applications
WO2007143039A2 (en) Oil- based mud resistivity imaging using resonant circuits
CN102400670A (zh) 一种三维感应测井仪器线圈系的设计方法
WO2004106973A2 (en) Electric power grid induced geophysical prospecting method and apparartus
US7714585B2 (en) Multi-frequency cancellation of dielectric effect
CN109915120B (zh) 一种基于环境因素的随钻电阻率测井系统校正方法
CN105891895B (zh) 一种确定天波传播特性的系统和方法
Duncan et al. SALTMAP—Airborne EM for the environment
CN105781538A (zh) 一种电磁波电阻率测井线圈系
Sørensen et al. TDEM in groundwater mapping—a continuous approach
CN112083507A (zh) 一种基于经验模态分解的瞬变电磁救援井探测方法
CN1492239A (zh) 高分辨率感应测井方法
CN111812726A (zh) 一种深层地质介电共振探测装置及其探测方法
CN110847887A (zh) 一种细粒沉积陆相页岩裂缝识别评价方法
CN104863578A (zh) 一种井中雷达的电导率测井方法
CN206071580U (zh) 随钻电阻率测量装置及其方位电阻率工具

Legal Events

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