CN102635347B - 一种定量将薄互地层等效成水平和垂直电阻率的方法 - Google Patents
一种定量将薄互地层等效成水平和垂直电阻率的方法 Download PDFInfo
- Publication number
- CN102635347B CN102635347B CN201210091428.2A CN201210091428A CN102635347B CN 102635347 B CN102635347 B CN 102635347B CN 201210091428 A CN201210091428 A CN 201210091428A CN 102635347 B CN102635347 B CN 102635347B
- Authority
- CN
- China
- Prior art keywords
- lambda
- exp
- resistivity
- centerdot
- prime
- 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
- Complex Calculations (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种电阻率测井仪数据处理方法。一种定量将电阻率各向同性薄互地层等效成水平和垂直电阻率的方法,对于一组电阻率各向同性的薄互地层,在仪器垂直于薄互地层时,选择一组发射频率和源距,计算仪器在垂直于薄互地层时的水平电阻率;在仪器倾斜于薄互地层时,选择与计算仪器垂直于薄互地层时相同的一组工作频率和源距,计算薄互地层垂直电阻率;根据获得的水平电阻率的值,通过反演方法获得垂直电阻率的值;将薄互地层等效成电阻率均匀各向异性地层,选取电阻率均匀各向异性地层作为反演地层模型,计算薄互地层等效的水平电阻率和垂直电阻率;以电阻率均匀各向异性介质中仪器的响应作为正演函数进行反演,得出薄互地层等效的垂直电阻率。
Description
技术领域
本发明涉及一种随钻电磁波电阻率测井仪数据处理方法,特别是涉及一种定量将电阻率各向同性薄互地层等效成水平和垂直电阻率的方法。
背景技术
随着陆上水平井和大斜度井钻井工作量增加以及海上钻井的需求,常规电缆测井已经不能满足测井技术的需要,因此随钻测井技术得到了非常迅速的发展。随钻电磁波电阻率测井仪器是随钻测井中最常用的仪器之一,它通过测量两个接收线圈上感应电动势的相位差和幅度比来获得地层的电阻率信息。地层中的电阻率并不全是均匀的,因为常有一些电阻率各向同性薄互层存在,这些薄互层呈现出宏观的电阻率各向异性(即水平电阻率和垂直电阻率不相等),而随钻电磁波电阻率类仪器在大斜度井中的响应受电阻率各向异性的影响很大,使得相位差和幅度比曲线出现分离的现象,但曲线的分离并不仅受电阻率各向异性的影响,泥浆的侵入及仪器在井眼中偏心也能引起曲线的分离,为了分析这种曲线分离现象产生的原因,需要弄清电阻率各向异性对随钻电磁波电阻率测井仪响应的影响规律,而电阻率各向异性通常以薄互层的形式存在,所以有必要对电阻率各向同性薄互层等效成宏观各向异性进行定量的转换。
发明内容
本发明针对现有技术不足,基于随钻电磁波电阻率类仪器在电阻率各向同性薄互地层中的响应,提出一种定量将电阻率各向同性薄互地层等效成水平和垂直电阻率的方法。
本发明所采用的技术方案:
一种定量将薄互地层等效成水平和垂直电阻率的方法,对于一组电阻率各向同性的薄互地层,采用随钻电磁波电阻率测井仪,分仪器垂直于薄互地层和倾斜于薄互地层两种情况计算该仪器在薄互地层中的响应,
1)、在仪器垂直于薄互地层时,选择一组发射频率和源距,采用纵向成层并矢格林函数法,计算出仪器通过薄互地层时的相位差视电阻曲线和幅度比视电阻率曲线,确定仪器在垂直于薄互地层时的水平电阻率;
2)、在仪器倾斜于薄互地层时,选择倾斜角度大于60°,将仪器以大于60°的倾斜角度通过薄互地层,选择与计算仪器垂直于薄互地层时相同的一组工作频率和源距,计算出仪器此时的相位差视电阻率和幅度比视电阻率响应,获得薄互地层垂直电阻率;
3)、根据获得的水平电阻率的值,将仪器倾斜于薄互地层时的两个响应值或两个响应值中质量较好的一个数据作为反演数据,通过反演方法获得垂直电阻率的值;
4)、考虑到薄互地层的总厚度足够大,使得随钻电磁波电阻率仪器在薄互地层中的响应不受围岩的影响,将薄互地层等效成电阻率均匀各向异性地层;根据薄互地层与电阻率均匀各向异性地层的等效关系,选取电阻率均匀各向异性地层作为反演地层模型,计算薄互地层等效的水平电阻率和垂直电阻率;
5)、反演时将仪器倾斜于薄互地层时计算出的相位差和幅度比视电阻率值作为已知响应,将仪器垂直于薄互地层时计算得出的视电阻率值作为已知的水平电阻率值,以电阻率均匀各向异性介质中该仪器的响应作为正演函数进行反演,得出薄互地层等效的垂直电阻率。
所述的定量将薄互地层等效成水平和垂直电阻率的方法,在步骤2)中,采用层状介质中的格林函数法计算随钻电磁波电阻率类仪器在薄互层下的响应:设层状单轴各向异性介质共有n+1层,各层编号为l=0,1,...,n,源在第j层,各层参数分别为μl、εvl、εhl,则 层厚hl=zl-zl-1(l=1,…,n-1)
层状单轴各向异性介质中的并矢格林函数GHM可表示为背景项和散射项之和,即GHM=PGHM+SGHM,对于非含源层(l≠j),则只有散射项。对于PGHM,可以根据均匀单轴各向异性介质中单位磁偶极子产生的Hertz势较容易地写出其解析式。对于任意l层内SGHM的各分量,由麦克斯韦方程组经计算可表示为如下形式,式中下脚标表示某方向的磁偶极子在某方向上产生的磁场,如表示z方向的磁偶极子在y方向上产生的散射场项。
式中的6个积分分别为:
所述的定量将薄互地层等效成水平和垂直电阻率的方法,利用仪器响应值进行反演时采用均匀各向异性介质中磁偶极子源电磁场公式:在正演模拟中将发射源看做磁偶极子,并设其随时间的变化关系为exp(iωt),其中ω为角频率,θ为井眼相对倾角,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,各向异性介质中方向单位磁偶极子产生的矢势为:
式中: μb为均匀介质磁导率,
各向异性介质中方向单位磁偶极子产生的矢势为:
式中:
各向异性介质中方向单位磁偶极子产生的矢势为:
由 及 其中 为电导率张量,可以得到均匀各向异性介质中沿不同方向单位磁偶极子产生的电场和磁场各分量的解析表达式,
在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和相位差ΔΦ定义为:
所述的定量将薄互地层等效成水平和垂直电阻率的方法,反演计算采用高斯牛顿梯度下降方法,均匀各向异性介质中的正演问题可以写成下面的非线性方程
其中表示未知的地层参数,表示不同的函数值,若已知实测的函数值则求解未知量的反演问题可归结为最小化与的方差问题,可写成如下形式
其中在给定初始值的情况下得到下面的线性最小方差问题
其中为函数值的差,为雅可比矩阵,每个元素为
根据给定的初始值求解(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)式,只需令满足即可。这样通过迭代的方法很容易得到
其中
本发明的有益效果:
1、本发明能对大斜度中出现的曲线分离现象进行解释,早先经常把曲线分离现象归结为泥浆侵入、仪器偏心以及围岩的影响,很少考虑到地层电阻率各向异性的影响,即薄互层的影响。现在利用我们的发明可以让技术人员了解薄互层与电阻率各向异性的定量转化关系,对曲线分析现象有了更深层次的认识,避免对曲线分离现象产生的原因做出错误的判断。提高了测井数据解释的质量。
2、在水平井中,随钻电磁波电阻率测井仪具有指导钻头在油层中钻进的地质导向功能,当油层为砂泥岩薄互层时,会呈现出宏观电阻率各向异性,仪器的响应会受各向异性的影响,利用本发明我们可以明确薄互层与电阻率各向异性间的定量转换关系,能够帮助我们更好的指导钻头最大限度地在油层中钻进,避免由于对薄互层与电阻率各向异性的定量关系不清楚而引起钻头钻出油层这类事件的发生,提高地质导向的质量,进而提高油井的油气产出效率,大大提高油气产量。
3、对于油层是砂泥岩薄互层的情况,在水平井中利用本发明能对该油层进行更准确的岩性和含油饱和度评价,更准确的计算出油层的储量和产出效率。另外利用邻井资料能根据仪器的响应结果对油层进行各向异性电阻率值的定量计算,实现更加准确的油层评价。
4、随钻电磁波电阻率仪器是最常用的随钻测井仪器之一,该仪器在大斜度井中的测量结果能反映出地层的各向异性,而各向异性是评价地层的一个非常重要的性质,且大多数的各向异性都是由于薄互层产生的,利用我们的发明可以更大的发挥这只仪器在储层评价中的作用,尤其是在薄互层中的作用,提高了该仪器的应用价值。
附图说明
图1:2MHz20in组合下随钻测井仪在薄互层中不同井斜下的响应;
图2:层状单轴各向异性介质地层模型;
图3:仪器和各向异性层状地层模型简图。
具体实施方式
实施例一:参见图1,本发明定量将薄互地层等效成水平和垂直电阻率的方法,对于一组电阻率各向同性的薄互地层,采用随钻电磁波电阻率测井仪,分仪器垂直于薄互地层和倾斜于薄互地层两种情况计算该仪器在薄互地层中的响应,其步骤包括:
1)、在仪器垂直于薄互地层时,选择一组发射频率和源距,采用纵向成层并矢格林函数法,计算出仪器通过薄互地层时的相位差视电阻曲线和幅度比视电阻率曲线,确定仪器在垂直于薄互地层时的水平电阻率;
2)、在仪器倾斜于薄互地层时,将仪器以大于60°的倾斜角度通过薄互地层,选择与仪器垂直于薄互地层时相同的一组发射频率和源距,计算出仪器此时的相位差视电阻率和幅度比视电阻率响应,获得薄互地层垂直电阻率;
3)、根据获得的水平电阻率的值,将仪器倾斜于薄互地层时的两个响应值或两个响应值中质量较好的一个数据作为反演数据,通过反演方法获得垂直电阻率的值;
4)、将薄互地层等效成电阻率均匀各向异性地层,根据薄互地层与电阻率均匀各向异性地层的等效关系,选取电阻率均匀各向异性地层作为反演地层模型,计算薄互地层等效的水平电阻率和垂直电阻率;
5)、反演时将仪器倾斜于薄互地层时计算出的相位差和幅度比视电阻率值作为已知响应,将仪器垂直于薄互地层时计算得出的视电阻率值作为已知的水平电阻率值,以电阻率均匀各向异性介质中该仪器的响应作为正演函数进行反演,得出薄互地层等效的垂直电阻率。
实施例二:本实施例与实施例一不同的是,在步骤2)中,采用层状介质中的格林函数法计算随钻电磁波电阻率类仪器在薄互层下的响应:
设层状单轴各向异性介质共有n+1层,各层编号为l=0,1,…,n,源在第j层,各层参数分别为μl、εvl、εhl,则 层厚hl=zl-zl-1(l=1,…,n-1)
层状单轴各向异性介质中的并矢格林函数GHM可表示为背景项和散射项之和,即GHM=PGHM+SGHM,对于非含源层(l≠j),则只有散射项;
对于任意l层内SGHM的各分量,由麦克斯韦方程组经计算可表示为如下形式,式中下脚标表示某方向的磁偶极子在某方向上产生的磁场,表示z方向的磁偶极子在y方向上产生的散射场项;
式中的6个积分分别为:
实施例三:本实施例的定量将薄互地层等效成水平和垂直电阻率的方法,与实施例二不同的是:利用仪器响应值进行反演时采用均匀各向异性介质中磁偶极子源电磁场公式:在正演模拟中将发射源看做磁偶极子,并设其随时间的变化关系为exp(iωt),其中ω为角频率,θ为井眼相对倾角,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,各向异性介质中方向单位磁偶极子产生的矢势为:
式中: μb为均匀介质磁导率,
各向异性介质中方向单位磁偶极子产生的矢势为:
式中:
各向异性介质中方向单位磁偶极子产生的矢势为:
由 及 其中 为电导率张量,可以得到均匀各向异性介质中沿不同方向单位磁偶极子产生的电场和磁场各分量的解析表达式,在xz平面内磁矩为MT(MT=ITNTAT)的发射天线可视为磁矩为MTx的水平方向磁偶极子与磁矩为MTz的垂直方向磁偶极子的叠加,其中IT、NT、AT分别为发射天线的电流强度、匝数和面积,
则发射天线在接收天线处产生的磁场的x分量和z分量为:
其中,分别为x方向和z方向单位磁偶极子在x方向产生的磁场,分别为x方向和z方向单位磁偶极子在z方向产生的磁场,
接收天线处的磁场强度为:
HR=Hxsinθ+Hzcosθ.(6)
由此得到接收天线的感应电动势为:
V=-iωμbHRNRAR,(7)
式中NR和AR分别为接收天线的匝数和面积,对单发双收三线圈系统,设近接收天线感应电动势的幅度和相位分别为|V1|和Φ1,远接收天线感应电动势的幅度和相位分别为|V2|和Φ2,则有:
由式(8),两个接收天线之间的幅度比A和相位差ΔΦ定义为:
实施例四:本实施例的定量将薄互地层等效成水平和垂直电阻率的方法,与前述各实施例不同的是:
反演计算采用高斯牛顿梯度下降方法,均匀各向异性介质中的正演问题可以写成下面的非线性方程
其中表示未知的地层参数,表示不同的函数值,若已知实测的函数值则求解未知量的反演问题可归结为最小化与的方差问题,可写成如下形式
其中在给定初始值的情况下得到下面的线性最小方差问题
其中为函数值的差,为雅可比矩阵,每个元素为
根据给定的初始值求解(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)的解,令满足通过迭代的方法得到
其中
实施例五;本实施例结合图1~图3,具体说明了本发明的一个实施方式。
选取模型的围岩电阻率为20Ohmm,目的层为层厚为2in的薄互层,薄互层电阻率分别为0.5Ohmm和10Ohmm,目的层厚5m,选择频率为2MHz,源距为20in的工作模式进行计算得出此地层下的响应曲线如图1,当倾角从0°增加到80°时,相位差和幅度比视电阻率都增大,且曲线开始出现分离,相位差视电阻率高于幅度比视电阻率,这些规律完全符合各向异性地层对曲线的影响规律,验证了薄互层地层成各向异性的特性。
从图1中得出直井中的电阻率为0.95Ohmm,将其作为水平电阻率值。井斜80°时的相位差视电阻率为3.66Ohmm,幅度比视电阻率为2.00Ohmm,将这两个值作为已知响应,利用高斯牛顿梯度下降反演方法对垂直电阻率进行反演,得出的垂直电阻率为4.71Ohmm。则利用本发明得出该薄互地层模型等效成宏观的水平电阻率为0.95Ohmm,垂直电阻率为4.71Ohmm。
计算随钻电磁波电阻率类仪器在薄互层下的响应采用纵向介质格林函数法,详述如下:
设层状单轴各向异性介质共有n+1层,见图2,各层编号为l=0,1,…,n,源在第j层,各层参数分别为μl、εvl、εhl,则 层厚hl=zl-zl-1(l=1,…,n-1)。只要给出层状单轴各向异性介质中沿不同方向单位磁偶极子产生的电场和磁场z分量的表达式,其切向量分量可由麦克斯韦方程组得到。
一、层状单轴各向异性介质中方向单位磁偶极子的轴向场
方向单位磁偶极子只产生TE波,其轴向场可表示为
(2)式最后一项为源项,对无源层,则没有该项。式中Al(λ)和Bl(λ)为待定系数,由层界面处电场和磁场的连续性条件确定,若l=0,只有A0(λ),若l=n,只有Bn(λ)(下同)。在层界面处,电场和磁场的切向量连续,即μHz和连续,由此可得到确定所有待定系数Al(λ)(l=0,1,…,n-1)和Bl(λ)(l=1,…,n)的线性方程组。经整理该方程组可表示为如下矩阵形式
AVXV=SV,(3)
式中:AV∈C2n×2n,XV,SV∈C2n。AV的各非零元素为,
i=2,…,n-1,其余元素为0。
XV的各元素为:
SV的各非零元素为:
若j=0,则只有和若j=n,则只有和其余元素为0
二、层状单轴各向异性介质中方向单位磁偶极子TM波的轴向场
方向单位磁偶极子既产生TM波又产生TE波。在层状单轴各向异性介质中TM波的轴向场可表示为:
在层界面处根据和εvEz连续可得到确定所有待定系数Cl(λ)(l=0,1,…,n-1)和Dl(λ)(l=1,…,n)的线性方程组,经整理可表示为如下矩阵形式
ATMXTM=STM.(6)
ATM的各非零元素为:
XTM的各元素为:
STM的各非零元素为:
三、层状单轴各向异性介质中方向单位磁偶极子TE波的轴向场在层状单轴各向异性介质中方向单位磁偶极子产生的TE波的轴向场可表示为:
根据层界面处电场和磁场的切向量连续性条件即和μHz连续,可得到确定所有待定系数El(λ)(l=0,1,…,n-1)和Fl(λ)(l=1,…,n)的线性方程组,经整理可表示为如下矩阵形式
ATEXTE=STE.(9)
ATE的各非零元素为:
XTE的各元素为:
STE的各非零元素为:
四、层状单轴各向异性介质中方向单位磁偶极子TM波和TE波的轴向场层状单轴各向异性介质中方向单位磁偶极子产生的TM波的轴向场可表示为:
方向单位磁偶极子产生的TE波的轴向场可表示为:
由上述各方向单位磁偶极子产生的电场和磁场的轴向分量,利用麦克斯韦方程组可计算得到电场和磁场的所有切向分量。
五、层状单轴各向异性介质中的并矢格林函数GEM和GHM
根据前面推导得到的电场和磁场的各分量,可写出层状单轴各向异性介质中的并矢格林函数GEMGHM。将GEM和GHM表示为背景项和散射项之和,即GEM=PGEM+SGEM,GHM=PGHM+SGHM,对于非含源层(l≠j),则只有散射项。对于PGEM和PGHM,可以根据均匀单轴各向异性介质中单位磁偶极子产生的矢势较容易地写出其解析式。对于任意l层内SGHM的各分量,经计算可表示为如下形式,
式中的6个积分分别为:
对于任意l层SGEM的各分量,可表示为如下形式
式中的5个积分分别为:
该计算方法不仅适用于计算仪器在薄互地层下的响应,还适用于计算每层电参数呈各向异性时仪器的响应。
在利用仪器响应值进行反演时采用的正演方法详述如下:
在正演模拟中将发射源看做磁偶极子,并设其随时间的变化关系为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。各向异性介质中方向单位磁偶极子产生的矢势为:
式中: μb为均匀介质磁导率。
各向异性介质中方向单位磁偶极子产生的矢势为:
式中:
各向异性介质中方向单位磁偶极子产生的矢势为:
由 及 其中 为电导率张量,可以得到均匀各向异性介质中沿不同方向单位磁偶极子产生的电场和磁场各分量的解析表达式。
由图3,在xz平面内磁矩为MT(MT=ITNTAT,IT、NT、AT分别为发射天线的电流强度、匝数和面积)的发射天线可视为磁矩为MTx的水平方向磁偶极子与磁矩为MTz的垂直方向磁偶极子的叠加。则发射天线在接收天线处产生的磁场的x分量和z分量为:
其中,分别为x方向和z方向单位磁偶极子在x方向产生的磁场,分别为x方向和z方向单位磁偶极子在z方向产生的磁场。
接收天线处的磁场强度为:
HR=Hxsinθ+Hzcosθ.(20)
由此得到接收天线的感应电动势为:
V=-iωμbHRNRAR,(21)
式中NR和AR分别为接收天线的匝数和面积。对单发双收三线圈系统,设近接收天线感应电动势的幅度和相位分别为|V1|和Φ1,远接收天线感应电动势的幅度和相位分别为|V2|和Φ2,则有:
由式(22),两个接收天线之间的幅度比A和相位差ΔΦ定义为:
此方法不仅适用于随钻电磁波电阻率类仪器,还适用于分析感应类及其他能将发射源和接收源看成磁偶极子的仪器。
高斯牛顿梯度下降反演方法
均匀各向异性介质中的正演问题可以写成下面的非线性方程
其牛表示未知的地层参数,表示不同的函数值,若已知实测的函数值则求解未知量的反演问题可归结为最小化与的方差问题,可写成如下形式
其中在给定初始值的情况下得到下面的线性最小方差问题
其中为函数值的差,为雅可比矩阵,每个元素为
根据给定的初始值求解(26)式后得到值,则经过一次计算后,之后计算值,如大于给定的误差限则将此作为初始值重新计算(26)式,直至计算出的值满足给定的误差限时,将此作为方程(25)解,迭代结束。
一般情况下,只有当值较小时可以用(26)式来近似(25)式,而在实际中由于雅可比矩阵可能存在奇异性,使(26)式的解变的较大,因此为了使迭代过程稳定,在解(26)式时需要限定约束条件其中Δ为给定的边界最大值,在此范围内可以认为(26)式是(25)式的线性近似。因此反演过程就转化为求解带约束条件的线性最小方程问题
完成(28)式的求解需要先进行正演计算得出进而得出然后利用正演通过微分得到雅可比矩阵之后便可求解出满足的校正量
利用上面介绍的正演计算方法,则雅可比矩阵可以通过(29)式计算得出
这样我们就得出了计算(28)式所需的基础条件。
对(28)式求解首先要对雅可比矩阵进行奇异值分解,如下式
其中和分别为M×N和N×N阶矩阵,且满足 为N×N阶的对角矩阵,包含奇异值
λ1≥λ2≥…≥λj≥…≥λp>0 (31)
和
λp+1=λp+2=…=λN=0
其中p为雅可比矩阵的秩。若不考虑约束条件(13)式的一般形式为
其中, 是对角矩阵,其元素
实际用迭代法解方程(32)时,由于非零奇异值有可能会很小,容易出现上溢,因此为了使迭代过程稳定,引进阻尼因子μ>0,将1/λj修正为根据方程(30)-(33),则方程(28)的解可以写成下式
当时,方程(13)的解为否则存在一个唯一的μ使
则相应的即为方程(28)的解。实际上并不需要非常精确的求解(35)式,只需令满足即可。这样通过迭代的方法很容易得到
其中
上面的反演方法就是高斯牛顿梯度下降方法,此方法具有收敛速度快迭代稳定的优点。
Claims (4)
1.一种定量将薄互地层等效成水平和垂直电阻率的方法,对于一组电阻率各向同性的薄互地层,采用随钻电磁波电阻率测井仪,分仪器垂直于薄互地层和倾斜于薄互地层两种情况计算该仪器在薄互地层中的响应,其特征是:
1)、在仪器垂直于薄互地层时,选择一组发射频率和源距,采用纵向成层并矢格林函数法,计算出仪器通过薄互地层时的相位差视电阻曲线和幅度比视电阻率曲线,确定仪器在垂直于薄互地层时的水平电阻率;
2)、在仪器倾斜于薄互地层时,将仪器以大于60°的倾斜角度通过薄互地层,选择与仪器垂直于薄互地层时相同的一组工作频率和源距,获得仪器此时的相位差视电阻率和幅度比视电阻率响应;
3)、根据获得的水平电阻率的值,将仪器倾斜于薄互地层时的两个视电阻率响应值或两个视电阻率响应值中质量较好的一个数据作为反演数据,通过反演方法获得垂直电阻率的值;
4)、将薄互地层等效成电阻率均匀各向异性地层,根据薄互地层与电阻率均匀各向异性地层的等效关系,选取电阻率均匀各向异性地层作为反演地层模型,计算薄互地层等效的水平电阻率和垂直电阻率;
5)、反演时将仪器倾斜于薄互地层时计算出的相位差和幅度比视电阻率值作为已知响应,将仪器垂直于薄互地层时计算得出的视电阻率值作为已知的水平电阻率值,以电阻率均匀各向异性介质中该仪器的响应作为正演函数进行反演,得出薄互地层等效的垂直电阻率。
2.根据权利要求1所述的定量将薄互地层等效成水平和垂直电阻率的方法,其特征是:在步骤2)中,采用层状介质中的格林函数法计算随钻电磁波电阻率类仪器在薄互层下的响应:设层状单轴各向异性介质共有n+1层,各层编号为l=0,1,…,n,源在第j层,各层参数分别为μl、εvl、εhl,则 层厚hl=zl-zl-1(l=1,…,n-1)
层状单轴各向异性介质中的并矢格林函数GHM可表示为背景项和散射项之和,即GHM=PGHM+SGHM,对于非含源层(l≠j),则只有散射项;
对于任意l层内SGHM的各分量,由麦克斯韦方程组经计算可表示为如下形式,式中下脚标表示某方向的磁偶极子在某方向上产生的磁场,表示z方向的磁偶极子在y方向上产生的散射场项;
式中的6个积分分别为:
3.根据权利要求1或2所述的定量将薄互地层等效成水平和垂直电阻率的方法,其特征是:利用仪器响应值进行反演时采用均匀各向异性介质中磁偶极子源电磁场公式:在正演模拟中将发射源看做磁偶极子,并设其随时间的变化关系为exp(iωt),其中ω为角频率,θ为井眼相对倾角,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,各向异性介质中方向单位磁偶极子产生的矢势为:
式中: μb为均匀介质磁导率,
各向异性介质中方向单位磁偶极子产生的矢势为:
式中:
各向异性介质中方向单位磁偶极子产生的矢势为:
由 及 其中 为电导率张量,可以得到均匀各向异性介质中沿不同方向单位磁偶极子产生的电场和磁场各分量的解析表达式,在xz平面内磁矩为MT(MT=ITNTAT)的发射天线可视为磁矩为MTx的水平方向磁偶极子与磁矩为MTz的垂直方向磁偶极子的叠加,其中IT、NT、AT分别为发射天线的电流强度、匝数和面积,
则发射天线在接收天线处产生的磁场的x分量和z分量为:
其中,分别为x方向和z方向单位磁偶极子在x方向产生的磁场,分别为x方向和z方向单位磁偶极子在z方向产生的磁场,
接收天线处的磁场强度为:
HR=Hxsinθ+Hzcosθ (6)
由此得到接收天线的感应电动势为:
V=-iωμbHRNRAR, (7)
式中NR和AR分别为接收天线的匝数和面积,对单发双收三线圈系统,设近接收天线感应电动势的幅度和相位分别为|V1|和Φ1,远接收天线感应电动势的幅度和相位分别为|V2|和Φ2,则有:
由式(8),两个接收天线之间的幅度比A和相位差ΔΦ定义为:
4.根据权利要求3所述的定量将薄互地层等效成水平和垂直电阻率的方法,其特征是:反演计算采用高斯牛顿梯度下降方法,均匀各向异性介质中的正演问题写成下面的非线性方程
其中 表示未知的地层参数, 表示不同的函数值,若已知实测的函数值则求解未知量的反演问题可归结为最小化与的方差问题,可写成如下形式
其中 在给定初始值的情况下得到下面的线性最小方差问题
其中为函数值的差,为雅可比矩阵,每个元素为
根据给定的初始值求解(12)式后得到值,则经过一次计算后,之后计算值,如大于给定的误差限则将此作为初始值重新计算(12)式,直至计算出的值满足给定的误差限时,将此作为方程(11)解,迭代结束;
为了使迭代过程稳定,在解(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)的解,令满足通过迭代的方法得到
其中
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210091428.2A CN102635347B (zh) | 2012-03-30 | 2012-03-30 | 一种定量将薄互地层等效成水平和垂直电阻率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210091428.2A CN102635347B (zh) | 2012-03-30 | 2012-03-30 | 一种定量将薄互地层等效成水平和垂直电阻率的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102635347A CN102635347A (zh) | 2012-08-15 |
CN102635347B true CN102635347B (zh) | 2014-12-10 |
Family
ID=46619908
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210091428.2A Expired - Fee Related CN102635347B (zh) | 2012-03-30 | 2012-03-30 | 一种定量将薄互地层等效成水平和垂直电阻率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102635347B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104314560B (zh) * | 2014-08-26 | 2017-10-17 | 中国石油天然气股份有限公司 | 水平井入窗前的非等比例地层对比方法和对比装置 |
CN104612661B (zh) * | 2014-12-09 | 2018-09-11 | 中国科学院声学研究所 | 一种随钻电磁波测井装置和方法 |
CN104775811B (zh) * | 2015-04-13 | 2017-09-01 | 中国海洋石油总公司 | 一种地层各向异性信息的提取及校正方法和系统 |
CN105089663B (zh) * | 2015-07-31 | 2017-11-07 | 中国石油天然气集团公司 | 一种大斜度/水平井地层电阻率各向异性校正方法 |
CN105507891B (zh) * | 2015-11-30 | 2018-06-01 | 中国石油天然气股份有限公司 | 获取大斜度井地层电阻率各向异性系数的方法及装置 |
CN106522932B (zh) * | 2016-09-28 | 2019-07-09 | 中国石油天然气股份有限公司 | 获取水平井各向异性地层电阻率的方法及装置 |
CN107121706A (zh) * | 2017-05-08 | 2017-09-01 | 厦门大学 | 基于波恩迭代法的航空瞬变电磁电导率三维反演方法 |
CN107907913A (zh) * | 2017-11-10 | 2018-04-13 | 西南石油大学 | 一种利用成像、倾角资料校正地层电阻率的方法 |
CN108169801B (zh) * | 2018-01-16 | 2020-09-15 | 陕西铁道工程勘察有限公司 | 高分辨地电阻率快速成像方法 |
CN108547611B (zh) * | 2018-03-01 | 2022-02-11 | 杭州迅美科技有限公司 | 水平井复杂环境随钻电磁波电阻率测井快速仿真方法 |
CN108398605B (zh) * | 2018-04-04 | 2023-10-20 | 中国人民解放军61489部队 | 地面核爆炸电磁脉冲复合环境模拟系统及模拟方法 |
CN109707375B (zh) * | 2018-12-26 | 2021-08-06 | 吉林大学 | 一种通过井下磁梯度张量数据辨识地下磁层信息的方法 |
CN112345989B (zh) * | 2020-11-18 | 2024-05-28 | 中国科学院电工研究所 | 一种肿瘤组织磁特性成像方法 |
CN117950065B (zh) * | 2024-03-26 | 2024-06-07 | 中国石油大学(华东) | 一种水平井阵列电阻率测井资料层厚校正方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101382599A (zh) * | 2007-09-03 | 2009-03-11 | 中国石油天然气集团公司 | 一种确定储层孔隙各向异性的瞬变电磁方法 |
CN102031964A (zh) * | 2009-10-05 | 2011-04-27 | 普拉德研究及开发股份有限公司 | 从三维感应测量值提取电阻率各向异性数据的多级工艺流程方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6819111B2 (en) * | 2002-11-22 | 2004-11-16 | Baker Hughes Incorporated | Method of determining vertical and horizontal resistivity, and relative dip in anisotropic earth formations having an arbitrary electro-magnetic antenna combination and orientation with additional rotation and position measurements |
US7617050B2 (en) * | 2007-08-09 | 2009-11-10 | Schlumberg Technology Corporation | Method for quantifying resistivity and hydrocarbon saturation in thin bed formations |
-
2012
- 2012-03-30 CN CN201210091428.2A patent/CN102635347B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101382599A (zh) * | 2007-09-03 | 2009-03-11 | 中国石油天然气集团公司 | 一种确定储层孔隙各向异性的瞬变电磁方法 |
CN102031964A (zh) * | 2009-10-05 | 2011-04-27 | 普拉德研究及开发股份有限公司 | 从三维感应测量值提取电阻率各向异性数据的多级工艺流程方法 |
Non-Patent Citations (4)
Title |
---|
径向成层介质的Green函数及其在随钻电磁波电阻率测量中的应用;杨锦舟等;《中国石油大学学报》;20090630;第33卷(第3期);第57-62页 * |
杨锦舟等.径向成层介质的Green函数及其在随钻电磁波电阻率测量中的应用.《中国石油大学学报》.2009,第33卷(第3期),第53-58页. * |
用并矢Green 函数的矢量本征函数展开式评价偏心对随钻电磁波电阻率测井响应的影响;魏宝君等;《中国石油大学学报》;20101031;第34卷(第5期);第53-58页 * |
魏宝君等.用并矢Green 函数的矢量本征函数展开式评价偏心对随钻电磁波电阻率测井响应的影响.《中国石油大学学报》.2010,第34卷(第5期),第57-62页. * |
Also Published As
Publication number | Publication date |
---|---|
CN102635347A (zh) | 2012-08-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102635347B (zh) | 一种定量将薄互地层等效成水平和垂直电阻率的方法 | |
CN102678106B (zh) | 随钻电磁波电阻率测井仪器的数据处理方法 | |
Chen et al. | Optimization of well trajectory under uncertainty for proactive geosteering | |
US8112227B2 (en) | Processing of multi-component induction measurements in a biaxially anisotropic formation | |
US9927551B2 (en) | Multifrequency processing to determine formation properties | |
EP2951619B1 (en) | Determination of true formation resistivity | |
US10416339B2 (en) | Homogeneous inversion for multi-component induction tools | |
CN103069304A (zh) | 增益校正式测量值 | |
Davydycheva et al. | Triaxial induction tool with electrode sleeve: FD modeling in 3D geometries | |
CN106446408A (zh) | 一种随钻补偿电磁波仪器的快速正反演处理方法 | |
Wilson et al. | Real-time 3D inversion of ultra-deep resistivity logging-while-drilling data | |
Liang et al. | Reservoir property mapping and monitoring from joint inversion of time-lapse seismic, electromagnetic, and production data | |
CN104903749A (zh) | 具有多极传感器的深方位系统 | |
Bhark et al. | A multiscale workflow for history matching in structured and unstructured grid geometries | |
Hong et al. | Novel and stable formulations for the response of horizontal-coil eccentric antennas in a cylindrically multilayered medium | |
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 | |
Nam et al. | Simulation of triaxial induction measurements in dipping, invaded, and anisotropic formations using a Fourier series expansion in a nonorthogonal system of coordinates and a self-adaptive hp finite-element method | |
Caudillo-Mata et al. | A framework for the upscaling of the electrical conductivity in the quasi-static Maxwell’s equations | |
CN102628360B (zh) | 随钻电磁波电阻率仪器的纵向分辨率指标的确定方法 | |
Bhark et al. | Multiscale parameterization and history matching in structured and unstructured grid geometries | |
Davydycheva | Two triaxial induction tools: sensitivity to radial invasion profile | |
Pardo et al. | Fast and automatic inversion of LWD resistivity measurements for petrophysical interpretation | |
CN107939385A (zh) | 定量计算极化值及应用的方法 | |
Yan et al. | Deep learning-assisted real-time forward modeling of electromagnetic logging in complex formations | |
Liang et al. | Electromagnetic telemetry simulation in vertical drillings |
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: 20141210 Termination date: 20190330 |
|
CF01 | Termination of patent right due to non-payment of annual fee |