CN102678106B - 随钻电磁波电阻率测井仪器的数据处理方法 - Google Patents
随钻电磁波电阻率测井仪器的数据处理方法 Download PDFInfo
- Publication number
- CN102678106B CN102678106B CN201210133884.9A CN201210133884A CN102678106B CN 102678106 B CN102678106 B CN 102678106B CN 201210133884 A CN201210133884 A CN 201210133884A CN 102678106 B CN102678106 B CN 102678106B
- Authority
- CN
- China
- Prior art keywords
- resistivity
- electromotive force
- induced electromotive
- transmitting coil
- force signal
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种随钻电磁波电阻率测井仪器的数据处理方法;含有以下步骤:1、使随钻电磁波电阻率测井仪器的两个接收线圈分别接收到感应电动势信号;2、计算两个感应电动势信号的相位差和幅度比;3、根据相位差和幅度比转换出视电阻率值;4、对视电阻率值进行井眼校正;5、设置直井或斜井标志;6、将地层看作一个纵向多层的结构,反演出每层的电阻率;7、将地层看作一个径向两层的结构,进行径向三参数反演;8、根据径向三参数反演的结果合成出具有固定探测深度的几条电阻率曲线;9、结束;本发明能去除多种因素对随钻电磁波电阻率测井仪器的测量数据的影响,从而得到地层真实的电参数及结构信息,提高了该仪器在地层评价中的应用价值。
Description
(一)、技术领域:本发明涉及一种数据处理方法,特别是涉及一种随钻电磁波电阻率测井仪器的数据处理方法。
(二)、背景技术:在石油钻井行业中,随着陆上水平井和大斜度井钻井工作量的增加以及海上钻井的需求,常规电缆测井已经不能满足测井技术的需要,因此,随钻测井技术得到了非常迅速的发展。随钻测井技术可以实现钻井和测井同时进行,它是将测井仪器安装在靠近钻头的部位,在地层未受到明显侵入和污染的条件下进行参数测量,随钻测井技术和传统的电缆测井相比较,具有实时性好、测井精度高等优点。随钻电磁波电阻率测井仪器是随钻测井中最常用的仪器之一,它主要测量地层的电阻率信息,由于一般情况下油层的电阻率较高,因此,它能够有效地识别油层,并还具有指导钻头在油层中水平钻进的地质导向功能。可见,随钻电磁波电阻率测井仪器在石油钻井中具有非常重要的实际意义,它能够增强随钻测井的能力,帮助油田找到更多的油气储层,缓解油气资源紧缺的局面。
随钻电磁波电阻率测井仪器在钻井环境中的电阻率响应受很多因素的影响,如井眼、侵入、围岩、各向异性等因素,这些因素的存在使得随钻电磁波电阻率测井仪器的电阻率响应值与地层的真实电阻率值存在较大的差异,如不对随钻电磁波电阻率测井仪器的电阻率响应值进行处理而直接利用,很有可能出现错误的结果,这样会使随钻电磁波电阻率测井仪器的实际应用效果大打折扣,因此,如何从随钻电磁波电阻率测井仪器的电阻率响应值中去除这些影响,从而获得地层真实的电阻率信息,将直接影响到该仪器的实际应用价值。
(三)、发明内容:
本发明要解决的技术问题是:克服现有技术的缺陷,提供一种随钻电磁波电阻率测井仪器的数据处理方法,该方法能去除多种因素对随钻电磁波电阻率测井仪器的测量数据的影响,从而得到地层真实的电参数及结构信息,提高了随钻电磁波电阻率测井仪器在地层评价中的应用价值。
本发明的技术方案:
一种随钻电磁波电阻率测井仪器的数据处理方法,含有以下步骤:
步骤1、使随钻电磁波电阻率测井仪器的两个接收线圈分别接收到感应电动势信号,从而分别得到两个接收线圈上的感应电动势信号的相位和幅度数据;
步骤2、计算两个接收线圈上的感应电动势信号的相位差和幅度比;
步骤3、利用电阻率转换图版,根据相位差和幅度比转换出视电阻率值;
步骤4、利用井眼校正图版对视电阻率值进行井眼校正,以去除井眼对视电阻率值的影响;一般情况下,井眼越大,视电阻率值越低,井眼对仪器响应的影响越大,越需要校正;如果井眼尺寸非常小,且采用的是油基泥浆,实际中对仪器的影响非常小的情况下也可以不进行井眼校正。
步骤5、判断井眼为直井或斜井?如为直井,设置直井标志;如为斜井,则设置斜井标志;
步骤6、将地层看作一个纵向多层的结构,反演出每层的电阻率,反演时去除了各层间的相互影响,如果井眼为斜井,在反演时还去除了井斜的影响;反演时,将每层看成是均匀各向同性的地层介质;
步骤7、将地层看作一个径向两层的结构,该径向两层结构从内到外依次为侵入带、原状地层,然后进行径向三参数反演,反演出侵入带电阻率、原状地层电阻率和侵入半径;在径向三参数反演中,反演方法采用高斯牛顿快速下降反演(Levenberg-Marquardt)算法,正演方法采用径向成层介质中的格林函数方法;高斯牛顿快速下降反演(Levenberg-Marquardt)算法适用于所有已知量大于待反演量的反演计算。
步骤8、根据径向三参数反演的结果合成出具有固定探测深度的几条电阻率曲线;
步骤9、结束。
在步骤5以后的过程中,当井眼为斜井时,再进行以下步骤:
步骤10、根据现场实测的泥浆电阻率的大小以及视电阻率曲线分离的次序来判断地层是否为各向异性?如否则执行步骤12,如是则执行步骤11;
判断地层各向异性的主要依据是视电阻率曲线的分离,但视电阻率曲线分离的现象并不只是各向异性引起的,还有可能是由于仪器靠近地层边界以及侵入等因素引起的。
步骤11、对地层各向异性进行反演,得出地层的水平电阻率和垂直电阻率;在对地层各向异性进行反演时,模型选择的是均匀各向异性介质,反演方法采用高斯牛顿快速下降反演(Levenberg-Marquardt)算法,正演方法采用磁偶极子源在均匀各向异性介质中电磁场的解析表达式,该解析表达式适用于计算所有可以将发射源和接收源看成磁偶极子的随钻电磁波电阻率类仪器在均匀各向异性介质中的响应。
步骤12、结束。
步骤5以后的过程为步骤5执行结束后,或者为步骤6执行结束后,或者为步骤7执行结束后,或者为步骤8执行结束后。
步骤1和步骤2中:随钻电磁波电阻率测井仪器的测量线圈中含有N个发射线圈和两个接收线圈,N个发射线圈分时与两个接收线圈匹配组成N组测量线圈;在每组测量线圈中,两个接收线圈分别接收到感应电动势信号;N为大于等于1的自然数;
离发射线圈距离近的接收线圈上的感应电动势信号的表达式为:其中,φ1为离发射线圈距离近的接收线圈上的感应电动势信号的相位,|V1|为离发射线圈距离近的接收线圈上的感应电动势信号的幅度,i为离发射线圈距离远的接收线圈上的感应电动势信号的表达式为:其中,φ2为离发射线圈距离远的接收线圈上的感应电动势信号的相位,|V2|为离发射线圈距离远的接收线圈上的感应电动势信号的幅度,i为
相位差的表达式为:
幅度比的表达式为:
或者,步骤1和步骤2中:随钻电磁波电阻率测井仪器的测量线圈中含有N对发射线圈和两个接收线圈,N对发射线圈分时与两个接收线圈匹配组成N组测量线圈;在每组测量线圈中,每个接收线圈均接收到一对发射线圈产生的两个感应电动势信号,对每一个发射线圈来说,得到一个分相位差和一个分幅度比,对一对发射线圈来说,得到两个分相位差和两个分幅度比,将这两个分相位差和两个分幅度比分别取平均值即可得到补偿后的相位差和幅度比;N为大于等于1的自然数;
对每一个发射线圈来说,离发射线圈距离近的接收线圈上的感应电动势信号的表达式为:其中,φ1为离发射线圈距离近的接收线圈上的感应电动势信号的相位,|V1|为离发射线圈距离近的接收线圈上的感应电动势信号的幅度,i为离发射线圈距离远的接收线圈上的感应电动势信号的表达式为:其中,φ2为离发射线圈距离远的接收线圈上的感应电动势信号的相位,|V2|为离发射线圈距离远的接收线圈上的感应电动势信号的幅度,i为
分相位差的表达式为:
分幅度比的表达式为:
步骤6中反演的具体步骤如下:
步骤6.1、对实际测得的分辨率较高的视电阻率曲线通过一次求导特征识别法定出地层边界位置以及每个地层介质电阻率的初始值;
步骤6.2、在每个地层找出距离中点最近的测量点作为反演点;
步骤6.3、在每一个反演点处用电阻率的初始值计算出测井响应,计算时考虑了井的倾斜度;
步骤6.4、在每一个反演点处将计算的测井响应与实测的测井响应进行比较,根据比较后的差值来校正每一反演点处的电阻率值;
步骤6.5、判断校正后的反演点处正演响应与实测的测井响应的误差是否小于给定的极小值?如否则执行步骤6.4,如是则执行步骤6.6;
步骤6.6、根据每一反演点的校正值绘制出校正后的原始地层电阻率参数曲线,由于用反演点处的电阻率作为该层的电阻率值,所以校正后的电阻率曲线是矩形曲线。
步骤6.3中的计算方法采用了层状介质中的格林函数法,该方法适用于计算所有可以将发射源和接收源看成磁偶极子的随钻电磁波类电阻率仪器在纵向成层各向异性介质中的响应。
步骤8的具体步骤如下:
步骤8.1、伪源距计算:根据侵入带电阻率、原状地层电阻率计算出探测深度与源距的关系曲线,根据该关系曲线,通过插值法得到固定探测深度所对应的源距,将此源距称为伪源距;
步骤8.2、固定深度处视电阻率计算:根据侵入带电阻率、原状地层电阻率和侵入半径,计算不同频率和源距下的视电阻率,得到不同频率下视电阻率与源距的关系曲线;将伪源距在该关系曲线上进行插值便得到固定探测深度处的视电阻率值。
●关于层状介质中的格林函数法的介绍如下:
设层状单轴各向异性介质共有n+1层,(如图2所示),各层编号为l=0,1,…,n,源在第j层,各层参数分别为μl、εvl、εhl,则 层厚hl=zl-zl-1(l=1,…,n-1)。只要给出层状单轴各向异性介质中沿不同方向单位磁偶极子产生的电场和磁场z分量的表达式,其切向量分量可由麦克斯韦方程组得到。
层状单轴各向异性介质中的并矢Green函数GHM可表示为背景项和散射项之和,即GHM=PGHM+SGHM,对于非含源层(l≠j),则只有散射项。对于PGHM,可以根据均匀单轴各向异性介质中单位磁偶极子产生的Hertz势较容易地写出其解析式。对于任意l层内SGHM的各分量,经计算可表示为如下形式,
式中的6个Sommerfeld积分分别为:
●关于步骤11中的正演方法介绍如下:
在正演模拟中将发射源看做磁偶极子,并设其随时间的变化关系为exp(iωt),其中ω为角频率。模型如图3所示,图中θ为井眼相对倾角。MT和MR为发射和接收天线磁矩,MTx、MTz和MRx、MRz分别为发射天线水平和垂直分量磁矩及接收天线水平和垂直分量磁矩。σhb为水平电导率,σvb为垂直电导率。
设源点位置坐标r′=(x′,y′,z′)、场点位置坐标r=(x,y,z),令 其中 并令K的表达式内σhb和σvb分别为向异性介质的水平和垂直复电导率,εhb和εvb分别为各向异性介质的水平和垂直复介电系数,它们的关系为σhb=iωεhb,σvb=iωεvb。各向异性介质中方向单位磁偶极子产生的Hertz势为:
式中: μb为均匀介质磁导率。
各向异性介质中方向单位磁偶极子产生的Hertz势为:
式中:
各向异性介质中方向单位磁偶极子产生的Hertz势为:
由 及 其中 为电导率张量,可以得到均匀各向异性介质中沿不同方向单位磁偶极子产生的电场和磁场各分量的解析表达式。
由图3,在xz平面内磁矩为MT(MT=ITNTAT,IT、NT、AT分别为发射天线的电流强度、匝数和面积)的发射天线可视为磁矩为MTx的水平方向磁偶极子与磁矩为MTz的垂直方向磁偶极子的叠加。则发射天线在接收天线处产生的磁场的x分量和z分量为:
其中,分别为x方向和z方向单位磁偶极子在x方向产生的磁场, 分别为x方向和z方向单位磁偶极子在z方向产生的磁场。
接收天线处的磁场强度为:
HR=Hxsinθ+Hzcosθ (6)
由此得到接收天线的感应电动势为:
V=-iωμbHRNRAR (7)
式中NR和AR分别为接收天线的匝数和面积。对单发双收三线圈系统,设近接收天线感应电动势的幅度和相位分别为|V1|和Φ1,远接收天线感应电动势的幅度和相位分别为|V2|和Φ2,则有:
由式(8),两个接收天线之间的幅度比A和相位差ΔΦ定义为:
●关于高斯牛顿快速下降反演(Levenberg-Marquardt)算法的介绍如下:
均匀各向异性介质中的正演问题可以写成下面的非线性方程
其中 表示未知的地层参数, 表示不同的函数值,若已知实测的函数值则求解未知量的反演问题可归结为最小化与的方差问题,可写成如下形式
其中 在给定初始值的情况下得到下面的线性最小方差问题
其中为函数值的差,为雅可比矩阵,每个元素为
根据给定的初始值求解(12)式后得到值,则经过一次计算后,之后计算值,如大于给定的误差限则将此作为初始值重新计算(12)式,直至计算出的值满足给定的误差限时,将此作为方程(11)解,迭代结束。
一般情况下,只有当值较小时可以用(12)式来近似(11)式,而在实际中由于雅可比矩阵可能存在奇异性,使(12)式的解变的较大,因此为了使迭代过程稳定,在解(12)式时需要限定约束条件其中Δ为给定的边界最大值,在此范围内可以认为(12)式是(11)式的线性近似。因此反演过程就转化为求解带约束条件的线性最小方程问题
完成(14)式的求解需要先进行正演计算得出进而得出然后利用正演通过微分得到雅可比矩阵之后便可求解出满足的校正量利用上面介绍的正演计算方法,则雅可比矩阵可以通过(15)式计算得出
这样我们就得出了计算(14)式所需的基础条件。
对(14)式求解首先要对雅可比矩阵进行奇异值分解,如下式
其中和分别为M×N和N×N阶矩阵,且满足 为N×N阶的对角矩阵,包含奇异值
λ1≥λ2≥…≥λj≥…≥λp>0 (17)
和
λp+1=λp+2=…=λN=0
其中p为雅可比矩阵的秩。若不考虑约束条件(14)式的一般形式为
其中, 是对角矩阵,其元素
实际用迭代法解方程(18)时,由于非零奇异值有可能会很小,容易出现上溢,因此为了使迭代过程稳定,引进阻尼因子μ>0,将1/λj修正为根据方程(16)-(19),则方程(14)的解可以写成下式
当时,方程(14)的解为否则存在一个唯一的μ使
则相应的即为方程(14)的解。实际上并不需要非常精确的求解(21)式,只需令满足即可。这样通过迭代的方法很容易得到
其中
高斯牛顿快速下降反演(Levenberg-Marquardt)算法具有收敛速度快、迭代稳定的优点。
●关于径向成层介质中的格林函数方法的介绍如下:
设径向成层介质共有n+1层(如图4所示),由里向外各层介质编号为l=0,1,…,n,源在第j层,各层参数分别为μl、εl,则κl 2=ω2μlεl 。可将径向成层介质任意第l层内的Green函数表示为如下形式,
式(24)等号右侧第一项为源项,对于非含源层(l≠j),则没有该项。Al(λ)和Bl(λ)为待定系数,由圆柱形层界面处电场和磁场的连续性条件确定。另由r→0和r→∞的边界条件可知,若l=0,则只有A0(λ),若l=n,则只有Bn(λ)。在层界面处,电场和磁场的切向量连续,即μΓ(r,r′)和连续,由此可得到确定所有待定系数Al(λ)(l=0,1,…,n-1)和Bl(λ)(l=1,…,n)的线性方程组。经整理该方程组可表示为如下矩阵形式
AX=S, (25)
式中:A∈C2n×2n,X,S∈C2n。A的各非零元素为:
X的各元素为:
X1=A0,X2n=Bn,X2i=Ai,X2i+1=Bi,(i=1,…,n-1)
S的各非零元素为:
若j=0,则只有S2j+1和S2j+2;若j=n,则只有S2j-1和S2j;其余元素为0。
由于矩阵A各行中间元素具有呈指数衰减的特点,式(25)可采用递推矩阵方法快速求解,递推过程不会出现溢出现象。这种求解含指数衰减项带状稀疏矩阵的递推算法无需进行矩阵求逆运算,只需要一次正向递推和一次逆向递推即可。其正向递推过程仅需要n次迭代,逆向递推过程亦仅需n次迭代,大大加快了计算速度。另外由于式(24)中的和已扣除指数项,该式的被积函数呈指数衰减,计算过程不存在上溢现象。将式(25)的计算结果代入式(24)后可得到各层的Green函数,并通过下面三个式子得到任意半径的圆形载流线圈在径向成层介质中产生的电磁场。
A(r,z)=μbNTIΓ(r,r′). (26)
由于接收点在第l层内,故式(27)中的μb须换成μl。只需改变式(25)内源项S的元素的位置,就可以方便地得到当源点r′和场点r在任意层时的所有待定系数,并进而计算出各层的Green函数。本发明表达式可适用于计算具有任意半径的圆形电流在不同介质层内产生的电磁场,其结果更具有普遍意义。另外本发明采用的递推矩阵方法很给出不同情况下径向成层介质中Green函数的通用表达式,无论是从表达方式的角度还是从编程的角度看均非常简洁。
本发明的有益效果:
1、本发明实现了对随钻电磁波电阻率测井仪器的数据处理从定性到定量的质变,如同把一个大脑植入了仪器之中,仪器的测量值如同眼睛看到的表面现象,经过本发明的数据处理流程的分析,量化出更进一步的油层厚度、侵入及原状地层等信息,这些信息能更好的被油田技术人员所应用,进一步量化出准确的含油饱和度,实现准确的储能及产能评价,并以此制定出合理的开采计划。
2、本发明不论对直井还是斜井都能将原始曲线处理成具有固定探测深度的曲线,这些能让油田技术人员直观地看到不同半径的地层剖面信息,一方面方便与阵列感应等其它具有多个探测深度的仪器测井曲线进行对比,另一方面便于分析地层的侵入特性,有利于评价地层的渗透率。
3、在水平井中,随钻电磁波电阻率测井仪器具有地质导向功能,本发明能对地质导向的效果进行判断,实现对地质导向钻井效果的评价,同时也可以帮助技术人员根据评价结果制定出相应的开采计划,实现产能的最大化。本发明提升了随钻电磁波电阻率测井仪器在地层评价中的应用价值,使该仪器成为了一种在评价地层及地质导向中有重要价值的仪器,同时也为今后开发出更新的方法来实现该仪器更大的价值提供了一种借鉴。
(四)、附图说明:
图1为随钻电磁波电阻率测井仪器的数据处理方法的流程示意图;
图2为层状单轴各向异性介质地层模型示意图;
图3为仪器和各向异性层状地层模型示意图;
图4为径向成层介质的模型示意图;
图5为三参数反演物理模型示意图;
图6为2MHz相位差视电阻率探测深度与源距的关系曲线示意图;
图7为2MHz相位差视电阻率与源距的关系曲线示意图;
图8为Oklahoma地层原始响应示意图;
图9为Oklahoma地层测井响应的一次导数曲线示意图;
图10为反演结果示意图。
(五)、具体实施方式:
实施例一:参见图1、图5~图10,随钻电磁波电阻率测井仪器的数据处理方法含有以下步骤(如图1所示):
步骤1、使随钻电磁波电阻率测井仪器的两个接收线圈分别接收到感应电动势信号,从而分别得到两个接收线圈上的感应电动势信号的相位和幅度数据;
步骤2、计算两个接收线圈上的感应电动势信号的相位差和幅度比;
步骤3、利用电阻率转换图版,根据相位差和幅度比转换出视电阻率值;
步骤4、利用井眼校正图版对视电阻率值进行井眼校正,以去除井眼对视电阻率值的影响;一般情况下,井眼越大,视电阻率值越低,井眼对仪器响应的影响越大,越需要校正;如果井眼尺寸非常小,且采用的是油基泥浆,实际中对仪器的影响非常小的情况下也可以不进行井眼校正。
步骤5、判断井眼为直井或斜井?如为直井,设置直井标志;如为斜井,则设置斜井标志;
步骤6、将地层看作一个纵向多层的结构,反演出每层的电阻率,反演时去除了各层间的相互影响,如果井眼为斜井,在反演时还去除了井斜的影响;反演时,将每层看成是均匀各向同性的地层介质;
步骤7、将地层看作一个径向两层的结构(如图5所示),该径向两层结构从内到外依次为侵入带、原状地层,然后进行径向三参数反演,反演出侵入带电阻率、原状地层电阻率和侵入半径;在径向三参数反演中,反演方法采用高斯牛顿快速下降反演(Levenberg-Marquardt)算法,正演方法采用径向成层介质中的格林函数方法;高斯牛顿快速下降反演(Levenberg-Marquardt)算法适用于所有已知量大于待反演量的反演计算。
表1是在假定侵入电阻率Rxo=1.0Ohmm,原状地层电阻率Rt=100.0Ohmm,侵入半径r为8.0in的地层参数下,利用径向成层介质格林函数计算得出的随钻电磁波电阻率测井响应。
表1
选择频率为2MHz,源距分别为20in、30in、46in时的相位差视电阻率为反演中的实测值,并设定迭代初始值为Rxo=2.0Ohmm,Rt=120.0Ohmm,侵入半径r=10.0in,迭代结果如下表,从表2中看出迭代8次后的结果与原始值能比较好的吻合。
表2
步骤8、根据径向三参数反演的结果合成出具有固定探测深度的几条电阻率曲线;
步骤9、结束。
在步骤5以后的过程中,当井眼为斜井时,再进行以下步骤:
步骤10、根据现场实测的泥浆电阻率的大小以及视电阻率曲线分离的次序来判断地层是否为各向异性?如否则执行步骤12,如是则执行步骤11;
判断地层各向异性的主要依据是视电阻率曲线的分离,但视电阻率曲线分离的现象并不只是各向异性引起的,还有可能是由于仪器靠近地层边界以及侵入等因素引起的。
步骤11、对地层各向异性进行反演,得出地层的水平电阻率和垂直电阻率;在对地层各向异性进行反演时,模型选择的是均匀各向异性介质,反演方法采用高斯牛顿快速下降反演(Levenberg-Marquardt)算法,正演方法采用磁偶极子源在均匀各向异性介质中电磁场的解析表达式,该解析表达式适用于计算所有可以将发射源和接收源看成磁偶极子的随钻电磁波电阻率类仪器在均匀各向异性介质中的响应。
假定水平电阻率Rh=10(Ohmm),垂直电阻率Rv=100(Ohmm),仪器相对倾角θ=80°,经正演计算后的结果见表4:
表4
利用表4中的数据进行反演,得到的反演结果见表5:
表5
步骤12、结束。
步骤5以后的过程为步骤5执行结束后,或者为步骤6执行结束后,或者为步骤7执行结束后,或者为步骤8执行结束后。
步骤1和步骤2中:随钻电磁波电阻率测井仪器的测量线圈中含有N个发射线圈和两个接收线圈,N个发射线圈分时与两个接收线圈匹配组成N组测量线圈;在每组测量线圈中,两个接收线圈分别接收到感应电动势信号;N为大于等于1的自然数;
离发射线圈距离近的接收线圈上的感应电动势信号的表达式为:其中,φ1为离发射线圈距离近的接收线圈上的感应电动势信号的相位,|V1|为离发射线圈距离近的接收线圈上的感应电动势信号的幅度,i为离发射线圈距离远的接收线圈上的感应电动势信号的表达式为:其中,φ2为离发射线圈距离远的接收线圈上的感应电动势信号的相位,|V2|为离发射线圈距离远的接收线圈上的感应电动势信号的幅度,i为
相位差的表达式为:
幅度比的表达式为:
步骤6中反演的具体步骤如下:
步骤6.1、对实际测得的分辨率较高的视电阻率曲线通过一次求导特征识别法定出地层边界位置以及每个地层介质电阻率的初始值;
步骤6.2、在每个地层找出距离中点最近的测量点作为反演点;
步骤6.3、在每一个反演点处用电阻率的初始值计算出测井响应,计算时考虑了井的倾斜度;
步骤6.4、在每一个反演点处将计算的测井响应与实测的测井响应进行比较,根据比较后的差值来校正每一反演点处的电阻率值;
步骤6.5、判断校正后的反演点处正演响应与实测的测井响应的误差是否小于给定的极小值?如否则执行步骤6.4,如是则执行步骤6.6;
步骤6.6、根据每一反演点的校正值绘制出校正后的原始地层电阻率参数曲线,由于用反演点处的电阻率作为该层的电阻率值,所以校正后的电阻率曲线是矩形曲线。
步骤6.3中的计算方法采用了层状介质中的格林函数法,该方法适用于计算所有可以将发射源和接收源看成磁偶极子的随钻电磁波类电阻率仪器在纵向成层各向异性介质中的响应。
以Oklahoma地层为例,2MHz、20in工作模式下直井中的相位差响应见图8,一次求导特征曲线见图9。将图9中导数的极大值点的深度坐标定为层界面,极小值点对应的视电阻率值作为迭代的初始值。采用下面的迭代方法进行迭代:
其中,为第j个点经过第n次迭代后的值,σj原始为第j个点原始实际测井值,为利用第n-1次的迭代结果计算得出的第j个点的模拟测井值。σn和βn分别为第n次迭代的权重因子。
将上例的井斜变为60°,利用上述反演方法计算可得到的反演结果如图10。
步骤8的具体步骤如下:
步骤8.1、伪源距计算:根据侵入带电阻率、原状地层电阻率计算出探测深度与源距的关系曲线,根据该关系曲线,通过插值法得到固定探测深度所对应的源距,将此源距称为伪源距;
在经过步骤7得到了Rxo和Rt的情况下,计算出探测深度与源距的关系曲线,以Rxo=1.0Ohmm、Rt=10.0Ohmm为例计算的结果如图6所示。
步骤8.2、固定深度处视电阻率计算:根据侵入带电阻率、原状地层电阻率和侵入半径,计算不同频率和源距下的视电阻率,得到不同频率下视电阻率与源距的关系曲线;将伪源距在该关系曲线上进行插值便得到固定探测深度处的视电阻率值。
在步骤7反演出了Rxo、Rt和ri的情况下,计算不同频率与源距下的视电阻率,以Rxo=1.0Ohmm、Rt=10.0Ohmm、ri=15.0in为例,计算结果如图7所示,观察图7,发现视电阻率与源距的关系曲线是条平滑的曲线。
步骤8的数值模拟如下:
选取固定探测深度为10in、20in、30in、40in和60in。在给定Rxo=1.0Ohmm、Rt=10.0Ohmm、ri=15.0in的初始条件下,利用步骤8中的两次插值便可以得出固定探测深度处的视电阻率,计算结果如表3。
表3
实施例二:参见图1、图5~图10,本实施例中的随钻电磁波电阻率测井仪器的数据处理方法与实施例一中相同的部分不再重述,不同之处是:
步骤1和步骤2中:随钻电磁波电阻率测井仪器的测量线圈中含有N对发射线圈和两个接收线圈,N对发射线圈分时与两个接收线圈匹配组成N组测量线圈;在每组测量线圈中,每个接收线圈均接收到一对发射线圈产生的两个感应电动势信号,对每一个发射线圈来说,得到一个分相位差和一个分幅度比,对一对发射线圈来说,得到两个分相位差和两个分幅度比,将这两个分相位差和两个分幅度比分别取平均值即可得到补偿后的相位差和幅度比;N为大于等于1的自然数;
对每一个发射线圈来说,离发射线圈距离近的接收线圈上的感应电动势信号的表达式为:其中,φ1为离发射线圈距离近的接收线圈上的感应电动势信号的相位,|V1|为离发射线圈距离近的接收线圈上的感应电动势信号的幅度,i为离发射线圈距离远的接收线圈上的感应电动势信号的表达式为:其中,φ2为离发射线圈距离远的接收线圈上的感应电动势信号的相位,|V2|为离发射线圈距离远的接收线圈上的感应电动势信号的幅度,i为
分相位差的表达式为:
分幅度比的表达式为:
Claims (7)
1.一种随钻电磁波电阻率测井仪器的数据处理方法,其特征是:含有以下步骤:
步骤1、使随钻电磁波电阻率测井仪器的两个接收线圈分别接收到感应电动势信号,从而分别得到两个接收线圈上的感应电动势信号的相位和幅度数据;
步骤2、计算两个接收线圈上的感应电动势信号的相位差和幅度比;
步骤3、利用电阻率转换图版,根据相位差和幅度比转换出视电阻率值;
步骤4、利用井眼校正图版对视电阻率值进行井眼校正,以去除井眼对视电阻率值的影响;
步骤5、判断井眼为直井或斜井?如为直井,设置直井标志;如为斜井,则设置斜井标志;
步骤6、将地层看作一个纵向多层的结构,反演出每层的电阻率,反演时去除了各层间的相互影响,如果井眼为斜井,在反演时还去除了井斜的影响;反演时,将每层看成是均匀各向同性的地层介质;
步骤7、将地层看作一个径向两层的结构,该径向两层结构从内到外依次为侵入带、原状地层,然后进行径向三参数反演,反演出侵入带电阻率、原状地层电阻率和侵入半径;在径向三参数反演中,反演方法采用高斯牛顿快速下降反演算法,正演方法采用径向成层介质中的格林函数方法;
步骤8、根据径向三参数反演的结果合成出具有固定探测深度的几条电阻率曲线;
步骤9、结束;
在步骤5以后的过程中,当井眼为斜井时,再进行以下步骤:
步骤10、根据现场实测的泥浆电阻率的大小以及视电阻率曲线分离的次序来判断地层是否为各向异性?如否则执行步骤12,如是则执行步骤11;
步骤11、对地层各向异性进行反演,得出地层的水平电阻率和垂直电阻率;在对地层各向异性进行反演时,模型选择的是均匀各向异性介质,反演方法采用高斯牛顿快速下降反演算法,正演方法采用磁偶极子源在均匀各向异性介质中电磁场的解析表达式;
步骤12、结束。
2.根据权利要求1所述的随钻电磁波电阻率测井仪器的数据处理方法,其特征是:所述步骤5以后的过程为步骤5执行结束后,或者为步骤6执行结束后,或者为步骤7执行结束后,或者为步骤8执行结束后。
3.根据权利要求1所述的随钻电磁波电阻率测井仪器的数据处理方法,其特征是:所述步骤1和步骤2中:随钻电磁波电阻率测井仪器的测量线圈中含有N个发射线圈和两个接收线圈,N个发射线圈分时与两个接收线圈匹配组成N组测量线圈;在每组测量线圈中,两个接收线圈分别接收到感应电动势信号;N为大于等于1的自然数;
离发射线圈距离近的接收线圈上的感应电动势信号的表达式为:其中,φ1为离发射线圈距离近的接收线圈上的感应电动势信号的相位,|V1|为离发射线圈距离近的接收线圈上的感应电动势信号的幅度,i为离发射线圈距离远的接收线圈上的感应电动势信号的表达式为:其中,φ2为离发射线圈距离远的接收线圈上的感应电动势信号的相位,|V2|为离发射线圈距离远的接收线圈上的感应电动势信号的幅度,i为
相位差的表达式为:
幅度比的表达式为:
4.根据权利要求1所述的随钻电磁波电阻率测井仪器的数据处理方法,其特征是:所述步骤1和步骤2中:随钻电磁波电阻率测井仪器的测量线圈中含有N对发射线圈和两个接收线圈,N对发射线圈分时与两个接收线圈匹配组成N组测量线圈;在每组测量线圈中,每个接收线圈均接收到一对发射线圈产生的两个感应电动势信号,对每一个发射线圈来说,得到一个分相位差和一个分幅度比,对一对发射线圈来说,得到两个分相位差和两个分幅度比,将这两个分相位差和两个分幅度比分别取平均值即可得到补偿后的相位差和幅度比;N为大于等于1的自然数;
对每一个发射线圈来说,离发射线圈距离近的接收线圈上的感应电动势信号的表达式为:其中,φ1为离发射线圈距离近的接收线圈上的感应电动势信号的相位,|V1|为离发射线圈距离近的接收线圈上的感应电动势信号的幅度,i为离发射线圈距离远的接收线圈上的感应电动势信号的表达式为:其中,φ2为离发射线圈距离远的接收线圈上的感应电动势信号的相位,|V2|为离发射线圈距离远的接收线圈上的感应电动势信号的幅度,i为
分相位差的表达式为:
分幅度比的表达式为:
5.根据权利要求1所述的随钻电磁波电阻率测井仪器的数据处理方法,其特征是:所述步骤6中反演的具体步骤如下:
步骤6.1、对实际测得的分辨率较高的视电阻率曲线通过一次求导特征识别法定出地层边界位置以及每个地层介质电阻率的初始值;
步骤6.2、在每个地层找出距离中点最近的测量点作为反演点;
步骤6.3、在每一个反演点处用电阻率的初始值计算出测井响应;
步骤6.4、在每一个反演点处将计算的测井响应与实测的测井响应进行比较,根据比较后的差值来校正每一反演点处的电阻率值;
步骤6.5、判断校正后的反演点处正演响应与实测的测井响应的误差是否小于给定的极小值?如否则执行步骤6.4,如是则执行步骤6.6;
步骤6.6、根据每一反演点的校正值绘制出校正后的原始地层电阻率参数曲线,由于用反演点处的电阻率作为该层的电阻率值,所以校正后的电阻率曲线是矩形曲线。
6.根据权利要求5所述的随钻电磁波电阻率测井仪器的数据处理方法,其特征是:所述步骤6.3中的计算方法采用了层状介质中的格林函数法。
7.根据权利要求1所述的随钻电磁波电阻率测井仪器的数据处理方法,其特征是:所述步骤8的具体步骤如下:
步骤8.1、伪源距计算:根据侵入带电阻率、原状地层电阻率计算出探测深度与源距的关系曲线,根据该关系曲线,通过插值法得到固定探测深度所对应的源距,将此源距称为伪源距;
步骤8.2、固定深度处视电阻率计算:根据侵入带电阻率、原状地层电阻率和侵入半径,计算不同频率和源距下的视电阻率,得到不同频率下视电阻率与源距的关系曲线;将伪源距在该关系曲线上进行插值便得到固定探测深度处的视电阻率值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210133884.9A CN102678106B (zh) | 2012-05-02 | 2012-05-02 | 随钻电磁波电阻率测井仪器的数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210133884.9A CN102678106B (zh) | 2012-05-02 | 2012-05-02 | 随钻电磁波电阻率测井仪器的数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102678106A CN102678106A (zh) | 2012-09-19 |
CN102678106B true CN102678106B (zh) | 2015-04-01 |
Family
ID=46810610
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210133884.9A Expired - Fee Related CN102678106B (zh) | 2012-05-02 | 2012-05-02 | 随钻电磁波电阻率测井仪器的数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102678106B (zh) |
Families Citing this family (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103015995B (zh) * | 2012-12-14 | 2015-06-24 | 中国电子科技集团公司第二十二研究所 | 通过带倾斜线圈的电阻率设备的地层倾角测量方法及装置 |
CN104453866B (zh) * | 2013-09-18 | 2019-04-23 | 中国石油集团长城钻探工程有限公司 | 感应测井用电流测量的相位标定方法 |
CN103485772B (zh) * | 2013-09-30 | 2016-08-17 | 中国电子科技集团公司第二十二研究所 | 测井设备、方法及装置 |
CN103558641A (zh) * | 2013-11-04 | 2014-02-05 | 中国石油天然气集团公司 | 测井探测器 |
CN104747177B (zh) * | 2013-12-31 | 2017-12-01 | 中国石油化工集团公司 | 利用刻度天线消除随钻电磁波电阻率系统误差的方法 |
CN105317428B (zh) * | 2014-06-24 | 2019-05-07 | 中国石油天然气集团公司 | 确定电磁信道模型的方法及装置 |
CN104088630B (zh) * | 2014-07-09 | 2016-08-10 | 中国石油集团长城钻探工程有限公司 | 随钻方位电阻率测井仪器获取地层边界距离的方法 |
CN104100261B (zh) * | 2014-07-09 | 2016-08-17 | 中国石油集团长城钻探工程有限公司 | 获取随钻方位电阻率测井仪器距地层边界距离的系统 |
CN104612661B (zh) * | 2014-12-09 | 2018-09-11 | 中国科学院声学研究所 | 一种随钻电磁波测井装置和方法 |
AU2015375485B2 (en) * | 2015-01-07 | 2018-03-15 | Halliburton Energy Services, Inc. | Functional earth model parameterization for resistivity inversion |
CN105863614A (zh) * | 2015-02-10 | 2016-08-17 | 中国石油集团长城钻探工程有限公司 | 一种三维感应测井数据实时处理方法 |
CN104895561A (zh) * | 2015-04-28 | 2015-09-09 | 中国海洋石油总公司 | 基于并行计算技术的随钻电阻率测井联合反演方法 |
CN106285655B (zh) * | 2015-06-05 | 2020-03-10 | 中国石油天然气股份有限公司 | 砂岩地层测井评估方法 |
CN104949692B (zh) * | 2015-06-25 | 2018-05-01 | 江苏中海达海洋信息技术有限公司 | 多波束测深系统安装校正值的计算方法 |
CN107035364B (zh) * | 2016-02-04 | 2021-09-14 | 中石化石油工程技术服务有限公司 | 一种井间电磁刻度方法 |
CN107152275A (zh) * | 2016-03-02 | 2017-09-12 | 中国石油化工股份有限公司 | 阻抗匹配电路及随钻电磁波电阻率测量仪 |
WO2018052820A1 (en) * | 2016-09-15 | 2018-03-22 | Shanjun Li | System and methodology of cross casing resistivity tool |
CN106446408A (zh) * | 2016-09-23 | 2017-02-22 | 上海神开石油设备有限公司 | 一种随钻补偿电磁波仪器的快速正反演处理方法 |
CN106407574A (zh) * | 2016-09-23 | 2017-02-15 | 上海神开石油设备有限公司 | 一种多分量随钻方位电磁波仪器的快速正反演处理方法 |
WO2018144029A1 (en) * | 2017-02-06 | 2018-08-09 | Halliburton Energy Services, Inc. | Multi-layer distance to bed boundary (dtbb) inversion with multiple initial guesses |
CN107191181B (zh) * | 2017-05-24 | 2021-01-05 | 中国石油大学(华东) | 一种基于电磁散射的井周界面探测方法 |
CN107461191B (zh) * | 2017-08-03 | 2021-09-14 | 中石化石油工程技术服务有限公司 | 一种随钻方位电磁波边界探测仪器温度刻度方法 |
CN107630697B (zh) * | 2017-09-26 | 2020-07-10 | 长江大学 | 基于随钻电磁波电阻率测井的地层电阻率联合反演方法 |
CN108019206B (zh) * | 2017-11-16 | 2018-10-12 | 中国石油集团长城钻探工程有限公司 | 一种高介电常数下随钻电磁波电阻率仪器量程扩展方法 |
CN108547611B (zh) * | 2018-03-01 | 2022-02-11 | 杭州迅美科技有限公司 | 水平井复杂环境随钻电磁波电阻率测井快速仿真方法 |
CN108643895B (zh) * | 2018-04-03 | 2022-03-01 | 中国石油天然气股份有限公司 | 一种渐变地层的电阻率测井响应计算方法及装置 |
CN110552689A (zh) * | 2018-05-15 | 2019-12-10 | 中国石油化工股份有限公司 | 一种确定随钻仪器到地层边界距离的方法 |
CN109209354B (zh) * | 2018-10-15 | 2019-11-22 | 中国石油大学(华东) | 一种时间域瞬变电磁波测井边界远探测方法 |
CN111127243B (zh) * | 2018-10-31 | 2023-12-01 | 中国石油化工集团有限公司 | 一种修正随钻电磁波电阻率转换曲线的方法 |
CN109711109B (zh) * | 2019-02-22 | 2023-12-22 | 中国石油天然气集团有限公司 | 一种电阻率测井仪器电极系结构参数智能优化的方法及装置 |
CN110005398B (zh) * | 2019-04-04 | 2024-03-22 | 中国石油大学(北京) | 随钻电磁波电阻率测井仪器设计参数确定方法及装置 |
CN111852458B (zh) * | 2019-04-10 | 2023-06-20 | 中国石油化工股份有限公司 | 一种随钻成像测井数据传输方法及系统 |
CN109915120B (zh) * | 2019-04-16 | 2022-06-03 | 北京恒泰万博石油技术股份有限公司 | 一种基于环境因素的随钻电阻率测井系统校正方法 |
CN112034257B (zh) * | 2020-09-11 | 2023-09-29 | 斯伦贝谢油田技术(山东)有限公司 | 一种井下电阻率的计算方法 |
CN112593919B (zh) * | 2020-12-01 | 2022-12-09 | 中海油田服务股份有限公司 | 一种校正电阻率方法、装置和存储介质 |
CN113204056B (zh) * | 2021-05-06 | 2023-06-13 | 中煤科工集团西安研究院有限公司 | 一种用于确定煤岩界面剖面分布位置的反演方法 |
CN116398126B (zh) * | 2023-06-07 | 2023-08-25 | 山东万洋石油科技有限公司 | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 |
CN117189089B (zh) * | 2023-09-06 | 2024-06-07 | 中国科学院地质与地球物理研究所 | 一种随钻前视测量设备的标定方法、装置及系统 |
CN117231212B (zh) * | 2023-11-16 | 2024-02-13 | 齐鲁工业大学(山东省科学院) | 一种随钻方位电磁波测井仪器和磁场方向调节的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1332380A (zh) * | 2000-07-07 | 2002-01-23 | 石油大学(北京) | 无趋附效应电磁波测井方法 |
CN2821156Y (zh) * | 2005-06-17 | 2006-09-27 | 中国石化集团胜利石油管理局测井公司 | 阵列感应测井仪器 |
CN101609169A (zh) * | 2009-04-08 | 2009-12-23 | 中国石油集团钻井工程技术研究院 | 一种提高电磁波电阻率测量精度和扩展其测量范围的方法 |
CN102121374A (zh) * | 2011-02-17 | 2011-07-13 | 长江大学 | 一种探测金属套管外地层电阻率的井中时域脉冲电磁法 |
-
2012
- 2012-05-02 CN CN201210133884.9A patent/CN102678106B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1332380A (zh) * | 2000-07-07 | 2002-01-23 | 石油大学(北京) | 无趋附效应电磁波测井方法 |
CN2821156Y (zh) * | 2005-06-17 | 2006-09-27 | 中国石化集团胜利石油管理局测井公司 | 阵列感应测井仪器 |
CN101609169A (zh) * | 2009-04-08 | 2009-12-23 | 中国石油集团钻井工程技术研究院 | 一种提高电磁波电阻率测量精度和扩展其测量范围的方法 |
CN102121374A (zh) * | 2011-02-17 | 2011-07-13 | 长江大学 | 一种探测金属套管外地层电阻率的井中时域脉冲电磁法 |
Also Published As
Publication number | Publication date |
---|---|
CN102678106A (zh) | 2012-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102678106B (zh) | 随钻电磁波电阻率测井仪器的数据处理方法 | |
CN110685600B (zh) | 一种用于地质导向的钻头调整预测方法 | |
CN102635347B (zh) | 一种定量将薄互地层等效成水平和垂直电阻率的方法 | |
Rosthal et al. | Field test results of an experimental fully-triaxial induction tool | |
US9322910B2 (en) | Method of real time subsurface imaging using electromagnetic data acquired from moving platforms | |
US8112227B2 (en) | Processing of multi-component induction measurements in a biaxially anisotropic formation | |
US8624969B2 (en) | Methods of electromagnetic migration imaging of geologic formation | |
Wang et al. | Efficient and reliable simulation of multicomponent induction logging response in horizontally stratified inhomogeneous TI formations by numerical mode matching method | |
US6950749B2 (en) | Method for resistivity anisotropy determination in near vertical wells | |
CN106407574A (zh) | 一种多分量随钻方位电磁波仪器的快速正反演处理方法 | |
US6442488B2 (en) | Inhomogeneous background based focusing method for multiarray induction measurements in a deviated well | |
Chen et al. | Inversion-based workflow for quantitative interpretation of the new-generation oil-based-mud resistivity imager | |
EP1718999A1 (en) | Method and apparatus for gradient electromagnetic induction well logging | |
Li et al. | Using new LWD measurements to evaluate formation resistivity anisotropy at any dip angle | |
Epov et al. | Comparison of the Russian and Western resistivity logs in typical Western Siberian reservoir environments: a numerical study | |
Wang et al. | Efficient finite-volume simulation of the LWD orthogonal azimuth electromagnetic response in a three-dimensional anisotropic formation using potentials on cylindrical meshes | |
CN108169802A (zh) | 一种粗糙介质模型的时域电磁数据慢扩散成像方法 | |
Chemali et al. | The effect of shale anisotropy on focused resistivity devices | |
CN102628360B (zh) | 随钻电磁波电阻率仪器的纵向分辨率指标的确定方法 | |
US20190079209A1 (en) | Multi-Well Resistivity Anisotropy Modeling, Used to Improve the Evaluation of Thinly Bedded Oil and Gas Reservoirs | |
Hou et al. | Characterization of formation fractures with multicomponent induction logging based on biaxial anisotropy models: Method and case studies | |
Faivre et al. | Using array induction and array laterolog data to characterize resistivity anisotropy in vertical wells | |
Rabinovich et al. | Multi-component induction logging: 10 years after | |
Tabanou et al. | Thinly laminated reservoir evaluation in oil-base mud: high resolution versus bulk anisotropy measurement-a comprehensive evaluation | |
Hou et al. | Integration of 3D modeling and real-time processing for enhancement of anisotropic formation evaluation with borehole multicomponent induction measurements |
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: 20150401 Termination date: 20190502 |
|
CF01 | Termination of patent right due to non-payment of annual fee |