CN108547611A - 水平井复杂环境随钻电磁波电阻率测井快速仿真方法 - Google Patents
水平井复杂环境随钻电磁波电阻率测井快速仿真方法 Download PDFInfo
- Publication number
- CN108547611A CN108547611A CN201810171477.4A CN201810171477A CN108547611A CN 108547611 A CN108547611 A CN 108547611A CN 201810171477 A CN201810171477 A CN 201810171477A CN 108547611 A CN108547611 A CN 108547611A
- Authority
- CN
- China
- Prior art keywords
- resistivity
- electric field
- component
- longitudinal
- model
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000004088 simulation Methods 0.000 title claims abstract description 33
- 241001074085 Scophthalmus aquosus Species 0.000 title claims abstract description 17
- 238000012937 correction Methods 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 7
- 230000007613 environmental effect Effects 0.000 claims abstract description 5
- 230000005684 electric field Effects 0.000 claims description 51
- 230000014509 gene expression Effects 0.000 claims description 24
- 230000015572 biosynthetic process Effects 0.000 claims description 16
- 230000004044 response Effects 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 230000005672 electromagnetic field Effects 0.000 claims description 8
- 238000005562 fading Methods 0.000 claims description 7
- 230000010363 phase shift Effects 0.000 claims description 7
- 230000009545 invasion Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 2
- 238000005520 cutting process Methods 0.000 claims 1
- 238000004519 manufacturing process Methods 0.000 claims 1
- 238000013517 stratification Methods 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000005553 drilling Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 230000006698 induction Effects 0.000 description 3
- 230000035699 permeability Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000013589 supplement Substances 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 229920002472 Starch Polymers 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 235000019698 starch Nutrition 0.000 description 1
- 239000008107 starch Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种水平井复杂环境随钻电磁波电阻率测井快速仿真方法,将不存在解析解的复杂的三维地层模型按矢量叠加原理进行拆分成存在解析解的纵向成层、径向成层两套简单结构模型,利用解析解分别计算,并将纵向成层模型仿真计算结果叠加到径向成层模型中,最终对径向成层模型数值仿真耦合了三维地层模型全参数,对最终径向成层介质仿真结果进行环境校正处理,得到最终的数值仿真结果。本方法解决了传统随钻测井有限元三维仿真计算速度慢、计算量大从而导致在工程应用中受限,同时也解决了三维复杂地层环境中随钻电磁波仿真不存在解析解,导致解析解无法应用在三维空间中随钻电磁仿真的问题。
Description
技术领域
本发明涉及油田开发技术领域,特别是涉及到一种水平井复杂环境随钻电磁波电阻率测井快速仿真方法。
背景技术
随着水平井和大斜度井钻井技术广泛应用,电缆测井应用受到一定制约,测井数据采集广泛采用随钻测量方式。随钻电阻率由于钻开时间短,测量点受泥浆侵入影响较小,且能够提供不同探测深度的测井曲线。随钻电阻率数据能反映地层真实的参数,因此随钻测井仪器有着广阔的应用前景。随着复杂油气田勘探开发的不断深入,大斜度井、水平井等复杂工艺井的广泛应用,随钻测井技术研究与随钻测井仪器研发得到足够的重视,并得到了快速的发展。
随钻电磁波电阻率类仪器电阻率测量与电缆侧向类仪器有很大不同,侧向类仪器利用直流电模式进行测量,通过一定的仪器常数进行刻度,即可把测量采集的信号转换成地层电阻率信号;随钻电磁波电阻率仪器是利用距发射源不同距离的两个接收线圈获得感应电磁波的相位移和幅度衰减变化来描述地层电阻率。相位移和幅度衰减与地层电阻率是一种非线性关系,没有固定的刻度系数,通常利用相位移-电阻率转换关系将相位移转换成相位移电阻率,利用幅度衰减-电阻率转换关系将幅度衰减转换成幅度衰减电阻率。不同仪器、不同源距、不同工作频率所对应的电阻率转换关系不同。
传统的随钻电磁波测井仪器发射线圈与接收线圈共轴,测量得到的地层信号为地层信息的平均值;近年来三大石油测井服务公司相继推出了具有方位探测能力的随钻方位电磁波测井仪器。2005年斯伦贝谢公司推出了PeriScope方位电阻率测量仪,2006年贝克休斯公司推出了随钻方位电磁波测井仪APR,2007年哈里伯顿公司推出了方位深电阻率测量仪ADR。方位电磁波测量仪器均采用轴向,倾斜或横向线圈混合,能够更好地提供有关地层方位的信息,指示地层的各向异性,识别地层边界。
目前,随钻电磁波电阻率仪器正演仿真主要有有限元方法、模式匹配法、伪解析解方法、积分方程方法等,以上的所有方法在二维模型中应用较为成熟,而对于复杂的三维地层模型,伪解析解、积分方程法、模式匹配法等已经不再适用,而有限元方法虽然可以模拟复杂三维结构,但是面临计算速度慢、计算要求计算机内存空间大等问题,严重制约了三维数值仿真的发展和应用。
为此,我们从随钻电磁波电阻率测量信号源头出发,将不存在解析解的复杂的三维地层模型按矢量叠加原理进行拆分,获得纵向成层、径向成层两套简单结构模型,纵向成层、径向成层地层模型电磁波数值仿真存在解析解,利用解析解分别计算,并将纵向成层模型仿真计算结果叠加到径向成层模型中,最终对径向成层模型数值仿真耦合了三维地层模型全参数,对最终径向成层介质仿真结果进行环境校正处理,得到最终的数值仿真结果。本方法解决了传统随钻测井有限元三维仿真计算速度慢、计算量大从而导致在工程应用中受限,同时也解决了三维复杂地层环境中随钻电磁波仿真不存在解析解,导致解析解无法应用在三维空间中随钻电磁仿真的问题。本技术将传统有限元仿真一个计算单元耗时大于30分钟,提高到仿真一个计算单元耗时小于10秒,大大提高了计算速度,解决了三维仿真工程应用技术瓶颈问题。
发明内容
本发明的目的是提供一种水平井复杂环境随钻电磁波电阻率测井快速仿真方法。
本发明的目的可通过如下技术措施来实现:
步骤1,依照随钻电磁波电阻率仪器工作原理,设定工作频率和收发线圈源距,将发射线圈等价为点源,根据单位磁偶极子发射源随时间变化关系,建立单位磁偶极子电磁场分量关系式;
步骤2,将水平井地层模型等价拆分包含井斜、层界面、原状地层水平电阻率、各向异性系数、地层介电常数、边界距参数的纵向成层介质模型和包含仪器结构、井眼直径、泥浆电阻率、井眼偏心距、偏心方位、侵入深度、侵入带电阻率、原状地层电阻率的径向成层介质模型;
步骤3,将步骤2纵向成层介质模型利用解析方法进行包括数值模拟、信号补偿、电阻率转换的操作,获得仪器纵向成层介质中的测井响应;
步骤4,将步骤3中的测井响应,替代步骤2中径向成层介质原状地层电阻率,调用径向成层解析解数值仿真方法,计算测井响应,对计算得到的信号进行包括对称补偿、电阻率工程转换的操作获得电阻率曲线;
步骤5,对计算后得到的电阻率曲线进行包括井眼环境校正、偏心校正的环境因素校正操作,最后获得输出电阻率仿真结果。
在步骤1中包含:
步骤11:单位磁偶极子源随时间的变化关系为eiωt,其中ω为角频率,得到方向单位磁偶极子在均匀各向异性介质中产生的Hertz势;
步骤12:将所述Hertz势经过转化处理后表示为Sommerfeld积分形式;
步骤13:由电磁场与Hertz势之间的关系式可以得到沿三个方向单位磁偶极子产生的电场和磁场各分量的解析式,并将解析式表示为Sommerfeld积分形式,其中电场和磁场z分量的Sommerfeld积分形式分别表示成波模积分的形式;
步骤14:根据Maxwell方程组,电场和磁场波模的切向分量与纵向分量之间的关系以及所述电场和磁场z分量的Sommerfeld积分形式,分别计算得到电场和磁场x、y分量的Sommerfeld积分表达式。
在步骤2中包含:
步骤21:按照三维空间地层电性分布和结构参数,充分考虑层、侵入、倾斜、井眼变化,建立三维原始模型;
步骤22:按照原始模型的层数、层界面位置、层内水平电阻率、垂直电阻率、地层倾角、地层方位角、井斜角、采样深度、采样间隔建立新纵向成层介质模型,纵向成层介质模型不考虑模型参数径向上变化;
步骤23:按照原始模型在径向上变化,逐采样点建立径向成层介质模型,模型参数与原始模型在对应采样点处井径、泥浆电阻率、侵入深度、侵入带电阻率、仪器偏心距一致,原状地层电阻率不予赋值。
步骤3包括:
步骤31:纵向共有n+1层介质中,各层层界面位置、层厚和地层参数不尽相同,各层编号为k=0,1,…,n,磁偶极子源在第j层,方向单位磁偶极子产生TE波,方向单位磁偶极子既产生TM波又产生TE波,方向单位磁偶极子也是既产生TM波又产生TE波。
步骤32:根据所述电场和磁场z分量的Sommerfeld积分形式以及层界面处电场和磁场的切向量连续,得到z分量的电场和磁场强度表达式。
步骤33:根据所述z分量的电场和磁场强度表达式以及Maxwell方程组电场和磁场波模的切向分量与纵向分量之间的关系,得到方向单位磁偶极子在x,y分量上的电场和磁场强度表达式。
步骤34:对推导得到的电场和磁场强度表达式进行积分计算,即可得到各个方向磁偶极子x,y,z分量的电场和磁场值。
步骤35:根据接收线圈处的电场强度或磁场强度计算得到接收线圈处的感应电动势,分别计算出两个接收线圈的相位和幅度。
步骤36:用两个接收线圈得到的相位做差,幅度做比值,得到相位移和幅度衰减。
步骤4包括:
步骤41:径向成层介质模型采用逐源距和逐个频率仿真,利用步骤3中纵向成层介质正演结果对应源距和频率下的电阻率值,赋值径向成层模型原状地层电阻率,将径向成层介质模型补充完整;
步骤42:根据环形电流源随时间的变化关系为exp(iωt),其中ω为角频率,有ω=2πf,f为交变电流频率;并假设在地层直角坐标系中,水平面为xy面,源点的位置坐标为rt=(xt,yt,zt),场点的位置坐标为r=(x,y,z),则方向单位磁偶极子在均匀各向异性介质中产生的Hertz势可表示为:
式中,μb为均匀介质磁导率,σhb为均匀各向异性介质的水平复电导率;
步骤43:建立电磁场与Hertz势之间的关系式:
其中E为电场强度,Π为Hertz势,i为虚部单位,为均匀各向异性介质的电导率张量;
步骤44:根据Maxwell方程组,建立电场和磁场波模的切向分量与纵向分量之间方程组:
式中表示磁场波模切向分量,表示磁场波模纵向分量,表示电场切向分量,表示电场纵向分量,μ为磁导率,εh为介电常数,λ为积分变量,为各向异性系数,σvb为均匀各向异性介质的垂向复电导率;
步骤45:求解上式得到线圈处的电场强度或磁场强度计算得到接收线圈处的感应电动势,进而计算得到两个接收线圈之间的相位差和幅度比值。
步骤46:将相位差值和幅度比值通过电阻率转换链表转换成电阻率值。
本发明有益效果:随钻电磁波电阻率仪器结构复杂,三维正演仿真难度大,在水平井复杂地层环境中,现有随钻电磁波电阻率正演仿真要么不适用,要么计算速度慢,制约数值仿真在仪器开发领域的应用效率,同时更严重制约正演仿真在实时地质导向、井场资料处理等领域的应用。本发明可有效提高正演仿真计算效率,将三维地质模型按照矢量合成方式拆分,通过纵向-径向数值仿真算法联合,解决了三维复杂地层环境中随钻电磁波仿真不存在解析解,导致解析法无法应用在三维空间中随钻电磁仿真的问题,且开发出新的水平井复杂环境下三维快速正演仿真算法,将传统有限元仿真一个计算单元耗时大于30分钟,提高到仿真一个计算单元耗时小于10秒,大大提高了计算效率,解决传统三维数值仿真在井场资料处理应用中效率低、速度慢的技术瓶颈问题。
附图说明
图1为本发明的一种水平井复杂环境随钻电磁波电阻率测井快速仿真方法具体实施例的流程图;
图2为纵向成层介质地层模型;
图3为本发明的幅度比-电阻率转换链表;
图4为本发明的相位差-电阻率转换链表;
图5为本发明的复杂模型正演结果示例。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明的水平井复杂环境随钻电磁波电阻率测井快速仿真方法,包括:
步骤1,依照随钻电磁波电阻率仪器工作原理,设定工作频率和收发线圈源距,将发射线圈等价为点源,根据单位磁偶极子发射源随时间变化关系,建立单位磁偶极子电磁场分量关系式,包含几个步骤:
1)设单位磁偶极子源随时间的变化关系为exp(iωt),其中ω为角频率,并假设在地层直角坐标系中(水平面为xy面)源点的位置坐标为rt=(xt,yt,zt),场点的位置坐标为r=(x,y,z),则方向单位磁偶极子在均匀各向异性介质中产生的Hertz势可表示为:
式中,μb为均匀介质磁导率,σhb为均匀各向异性介质的水平复电导率。
2)经过转化处理,(11)式可以表示为如下Sommerfeld积分形式:
式中,Jv为v阶Bessel函数,λ为积分变量。
方向单位磁偶极子在均匀各向异性介质中产生的Hertz势可以表示为:
式中,为各向异性系数,i为复数单位,σhb、σvb分别为均匀各向异性介质的水平、垂向复电导率。经推导(13)、(14)式可以分别表示为如下Sommerfeld积分形式:
式中,
3)由电磁场与Hertz势之间的关系式
其中为均匀各向异性介质的电导率张量。将(2)、(5)、(6)式带入(7)式,可以得到沿三个方向单位磁偶极子产生的电场和磁场各分量的解析式及Sommerfeld积分形式,其中电场和磁场z分量的Sommerfeld积分形式分别表示为:
上述各分量均被表示成了波模积分的形式。例如
式中为某一λ对应的波模。
4)根据Maxwell方程组,电场和磁场波模的切向分量与纵向分量之间的关系可表示为:
步骤2,如图2所示,将水平井地层模型等价拆分包含井斜、层界面、原状地层水平电阻率、各向异性系数、地层介电常数、边界距参数的纵向成层介质模型和包含仪器结构、井眼直径、泥浆电阻率、井眼偏心距、偏心方位、侵入深度、侵入带电阻率、原状地层电阻率的径向成层介质模型,包含几个步骤:
1)按照三维空间地层电性分布和结构参数,充分考虑层、侵入、倾斜、井眼等变化,建立三维初始模型;
2)按照原始三维模型的层数、层界面位置、层内水平电阻率、垂直电阻率、地层倾角、地层方位角、井斜角、采样深度、采样间隔建立新纵向成层介质模型,纵向成层介质模型不考虑模型参数径向上变化;
3)按照原始模型在径向上变化,逐采样点建立径向成层介质模型,模型参数与原始模型在对应采样点处井径、泥浆电阻率、侵入深度、侵入带电阻率、仪器偏心距一致,原状地层电阻率未赋值。
步骤3,将步骤2纵向成层介质模型利用解析解方法进行数值模拟、信号补偿、电阻率转换等操作,获得仪器纵向成层介质中的测井响应,包括几个步骤:
1)设纵向成层各向异性介质共有n+1层,各层编号为k=0,1,…,n,磁偶极子源在第j层,各层层界面位置、层厚和地层参数如图2所示。方向单位磁偶极只产生TE波,方向单位磁偶极子既产生TM波又产生TE波,方向单位磁偶极子也是既产生TM波又产生TE波。
2)根据式(18a)-(18f)可以得到z分量的电场和磁场强度表达式为:
上式中l=0,1,…,n, μl表示第l层磁导率,σhl表示第l层水平电导率,σvl表示第l层垂直电导率,当l=j时δlj=1,否则δlj=0;Al(λ),Bl(λ),Cl(λ),Dl(λ),El(λ),Fl(λ)为待定系数,由层界面处,电场和磁场的切向量连续得到。
3)将式(20a)-(23d)中对应的变量代入式(19a)-(19d)中,即可以计算得到 方向单位磁偶极子在x,y分量上的电场和磁场强度。
4)对推导得到的电场和磁场强度表达式进行积分计算,即可得到各个方向磁偶极子x,y,z分量的电场和磁场值。
5)然后根据接收线圈处的电场强度或磁场强度计算得到接收线圈处的感应电动势,利用距发射线圈不同距离接受线圈感应电动势,可以计算出两接收线圈相位(式23)和幅度(式24):
式中PS表示相位,AT表示幅度,βi为方位角,为方位角为βi视接收线圈电压。
6)用两个接收线圈得到的相位做差,两个接收线圈得到的幅度做比值,即可得到相位移和幅度衰减。
步骤4,将步骤3中的测井响应,替代步骤2中径向成层介质原状地层电阻率,调用径向成层解析解数值仿真方法,计算测井响应,对计算得到的信号进行对称补偿、电阻率工程转换等获得电阻率曲线,包括以下步骤:
1)径向成层介质模型采用逐源距和频率仿真,利用步骤3中纵向成层介质正演结果对应源距和频率下的电阻率值,赋值径向成层模型原状地层电阻率,将径向成层介质模型补充完整;
2)根据环形电流源随时间的变化关系为exp(iωt),其中ω为角频率,有ω=2πf,f为交变电流频率;并假设在地层直角坐标系中,水平面为xy面,源点的位置坐标为rt=(xt,yt,zt),场点的位置坐标为r=(x,y,z),则方向单位磁偶极子在均匀各向异性介质中产生的Hertz势可表示为:
式中,μb为均匀介质磁导率,σhb为均匀各向异性介质的水平复电导率;
3)建立电磁场与Hertz势之间的关系式:
其中E为电场强度,Π为Hertz势,i为虚部单位,为均匀各向异性介质的电导率张量;
4)根据Maxwell方程组,建立电场和磁场波模的切向分量与纵向分量之间方程组:
式中表示磁场波模切向分量,表示磁场波模纵向分量,表示电场切向分量,表示电场纵向分量,μ为磁导率,εh为介电常数,λ为积分变量,为各向异性系数,σvb为均匀各向异性介质的垂向复电导率;
5)求解上式线圈处的电场强度或磁场强度,计算得到接收线圈处的感应电动势,进而计算得到两个接收线圈之间的相位差和幅度比值。
6)将相位差值和幅度比值通过电阻率转换链表如图3和图4所示,转换成电阻率值。
步骤5,对计算后得到的电阻率曲线进行井眼环境校正、偏心校正等环境因素校正,最后获得输出电阻率仿真结果,如图5所示,图中第一道为电阻率剖面道,第二道为侵入深度反演结果道,第三道为高频相位差电阻率道,第四道为高频幅度比电阻率道。
本发明可以解决现有数值仿真方法在复杂地层环境下随钻电磁波电阻率仪器数值仿真计算速度慢的问题,大大提高计算效率,节约计算时间。
以上,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本方法领域的方法人员在本发明揭露的方法范围内,根据本发明的方法方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (5)
1.水平井复杂环境随钻电磁波电阻率测井快速仿真方法,其特征在于,该电阻率测井快速仿真方法包括:
步骤1,依照随钻电磁波电阻率仪器工作原理,设定工作频率和收发线圈源距,将发射线圈等价为点源,根据单位磁偶极子发射源随时间变化关系,建立单位磁偶极子电磁场分量关系式;
步骤2,将水平井地层模型等价拆分包含井斜、层界面、原状地层水平电阻率、各向异性系数、地层介电常数、边界距参数的纵向成层介质模型和包含仪器结构、井眼直径、泥浆电阻率、井眼偏心距、偏心方位、侵入深度、侵入带电阻率、原状地层电阻率的径向成层介质模型;
步骤3,将步骤2纵向成层介质模型利用解析解方法进行包括数值模拟、信号补偿、电阻率转换的操作,获得仪器纵向成层介质中的测井响应;
步骤4,将步骤3中的测井响应,替代步骤2中径向成层介质原状地层电阻率,调用径向成层解析解数值仿真方法,计算测井响应,对计算得到的信号进行包括对称补偿、电阻率工程转换的操作获得电阻率曲线;
步骤5,对计算后得到的电阻率曲线进行包括井眼环境校正、偏心校正的环境因素校正操作,最后获得输出电阻率仿真结果。
2.根据权利要求1所述的水平井复杂环境随钻电磁波电阻率测井快速仿真方法,其特征在于,在步骤1中,依照随钻电磁波电阻率仪器工作原理,根据单位磁偶极子发射源随时间变化关系,建立单位磁偶极子电磁场分量关系式,包含以下几个步骤:
1)单位磁偶极子源随时间的变化关系为eiωt,其中ω为角频率,得到 方向单位磁偶极子在均匀各向异性介质中产生的Hertz势;
2)将所述Hertz势经过转化处理后表示为Sommerfeld积分形式;
3)由电磁场与Hertz势之间的关系式可以得到沿三个方向单位磁偶极子产生的电场和磁场各分量的解析式,并将解析式表示为Sommerfeld积分形式,其中电场和磁场z分量的Sommerfeld积分形式分别表示成波模积分的形式;
4)根据Maxwell方程组,电场和磁场波模的切向分量与纵向分量之间的关系以及所述电场和磁场z分量的Sommerfeld积分形式,分别计算得到电场和磁场x、y分量的Sommerfeld积分表达式。
3.根据权利要求1所述的水平井复杂环境随钻电磁波电阻率测井快速仿真方法,其特征在于,在步骤2中,将水平井地层模型等价拆分纵向成层介质模型和径向成层介质模型,包含以下几个步骤:
1)按照三维空间地层电性分布和结构参数,充分考虑层、侵入、倾斜、井眼变化,建立三维原始模型;
2)按照原始模型的层数、层界面位置、层内水平电阻率、垂直电阻率、地层倾角、地层方位角、井斜角、采样深度、采样间隔建立新纵向成层介质模型,纵向成层介质模型不考虑模型参数径向上变化;
3)按照原始模型参数在径向上变化,逐采样点建立径向成层介质模型,模型参数与原始模型在对应采样点处井径、泥浆电阻率、侵入深度、侵入带电阻率、仪器偏心距一致,原状地层电阻率不予赋值。
4.根据权利要求2所述的水平井复杂环境随钻电磁波电阻率测井快速仿真方法,其特征在于,在步骤3中,将步骤2纵向成层介质模型利用解析解方法进行数值模拟,获得仪器纵向成层介质中的测井响应,包含以下几个步骤:
1)纵向共有n+1层介质中,各层层界面位置、层厚和地层参数不尽相同,各层编号为k=0,1,…,n,磁偶极子源在第j层,方向单位磁偶极子产生TE波,方向单位磁偶极子既产生TM波又产生TE波,方向单位磁偶极子也是既产生TM波又产生TE波;
2)根据所述电场和磁场z分量的Sommerfeld积分形式以及层界面处电场和磁场的切向量连续,得到z分量的电场和磁场强度表达式;
3)根据所述z分量的电场和磁场强度表达式以及Maxwell方程组电场和磁场波模的切向分量与纵向分量之间的关系,得到方向单位磁偶极子在x,y分量上的电场和磁场强度表达式;
4)对推导得到的电场和磁场强度表达式进行积分计算,即可得到各个方向磁偶极子x,y,z分量的电场和磁场值;
5)根据接收线圈处的电场强度或磁场强度计算得到接收线圈处的感应电动势,分别计算出两个接收线圈的相位和幅度;
6)用两个接收线圈得到的相位做差,幅度做比值,得到相位移和幅度衰减。
5.根据权利要求1所述的水平井复杂环境随钻电磁波电阻率测井快速仿真方法,其特征在于,在步骤4中,将步骤2径向成层介质模型利用解析解方法进行数值模拟,获得仪器径向成层介质中的测井响应,包含以下几个步骤:
1)径向成层介质模型采用逐源距和频率仿真,利用步骤3中纵向成层介质正演结果对应源距和频率下的电阻率值,赋值径向成层模型原状地层电阻率,将径向成层介质模型补充完整;
2)根据环形电流源随时间的变化关系为exp(iωt),其中ω为角频率,有ω=2πf,f为交变电流频率;并假设在地层直角坐标系中,水平面为xy面,源点的位置坐标为rt=(xt,yt,zt),场点的位置坐标为r=(x,y,z),则方向单位磁偶极子在均匀各向异性介质中产生的Hertz势可表示为:
式中,μb为均匀介质磁导率,σhb为均匀各向异性介质的水平复电导率;
3)建立电磁场与Hertz势之间的关系式:
及
其中E为电场强度,∏为Hertz势,i为虚部单位,为均匀各向异性介质的电导率张量;
4)根据Maxwell方程组,建立电场和磁场波模的切向分量与纵向分量之间方程组:
式中表示磁场波模切向分量,表示磁场波模纵向分量,表示电场切向分量,表示电场纵向分量,μ为磁导率,εh为介电常数,λ为积分变量,为各向异性系数,σvb为均匀各向异性介质的垂向复电导率;
5)求解上式线圈处的电场强度或磁场强度,计算得到接收线圈处的感应电动势,进而计算得到两个接收线圈之间的相位差和幅度比值。
6)将相位差值和幅度比值通过电阻率转换链表转换成电阻率值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810171477.4A CN108547611B (zh) | 2018-03-01 | 2018-03-01 | 水平井复杂环境随钻电磁波电阻率测井快速仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810171477.4A CN108547611B (zh) | 2018-03-01 | 2018-03-01 | 水平井复杂环境随钻电磁波电阻率测井快速仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108547611A true CN108547611A (zh) | 2018-09-18 |
CN108547611B CN108547611B (zh) | 2022-02-11 |
Family
ID=63516077
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810171477.4A Active CN108547611B (zh) | 2018-03-01 | 2018-03-01 | 水平井复杂环境随钻电磁波电阻率测井快速仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108547611B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110454154A (zh) * | 2019-08-22 | 2019-11-15 | 中国科学院地质与地球物理研究所 | 一种校正井眼偏心对方位信号影响的方法及装置 |
CN110989010A (zh) * | 2019-12-02 | 2020-04-10 | 吉林大学 | 一种基于偏心线圈补偿的瞬变电磁勘探装置及勘探方法 |
CN111175834A (zh) * | 2020-01-16 | 2020-05-19 | 吉林大学 | 一种基于全程采集技术的拖曳式瞬变电磁探测装置及方法 |
CN111581791A (zh) * | 2020-04-22 | 2020-08-25 | 中国海洋石油集团有限公司 | 一种电阻率正演仿真方法、装置及计算机可读存储介质 |
CN113076508A (zh) * | 2021-04-02 | 2021-07-06 | 北京环境特性研究所 | 基于垂直磁偶极子在半空间下的低频近场快速计算方法 |
CN115062474A (zh) * | 2022-06-17 | 2022-09-16 | 上海玫克生储能科技有限公司 | 一种锂电池的工况仿真方法、系统、终端设备及存储介质 |
CN116398126A (zh) * | 2023-06-07 | 2023-07-07 | 山东万洋石油科技有限公司 | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 |
CN111581791B (zh) * | 2020-04-22 | 2024-05-31 | 中国海洋石油集团有限公司 | 一种电阻率正演仿真方法、装置及计算机可读存储介质 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1378652A (zh) * | 1999-09-15 | 2002-11-06 | 埃克森美孚上游研究公司 | 远程储层电阻率映像 |
US20110231098A1 (en) * | 2009-10-05 | 2011-09-22 | Dzevat Omeragic | Multilevel workflow method to extract resistivity anisotropy data from 3d induction measurements |
CN102619504A (zh) * | 2012-04-17 | 2012-08-01 | 中国电子科技集团公司第二十二研究所 | 随钻电磁波电阻率仪器的径向探测深度指标的确定方法 |
CN102635347A (zh) * | 2012-03-30 | 2012-08-15 | 中国电子科技集团公司第二十二研究所 | 一种定量将薄互地层等效成水平和垂直电阻率的方法 |
CN102678106A (zh) * | 2012-05-02 | 2012-09-19 | 中国电子科技集团公司第二十二研究所 | 随钻电磁波电阻率测井仪器的数据处理方法 |
CN103323883A (zh) * | 2013-06-26 | 2013-09-25 | 中南大学 | 一种近场磁偶源高密度的几何电阻率测深法 |
CN104265279A (zh) * | 2014-07-30 | 2015-01-07 | 中国石油集团川庆钻探工程有限公司 | 断层条件下随钻测井曲线预测方法 |
CN106407574A (zh) * | 2016-09-23 | 2017-02-15 | 上海神开石油设备有限公司 | 一种多分量随钻方位电磁波仪器的快速正反演处理方法 |
CN106446408A (zh) * | 2016-09-23 | 2017-02-22 | 上海神开石油设备有限公司 | 一种随钻补偿电磁波仪器的快速正反演处理方法 |
CN107561244A (zh) * | 2017-04-20 | 2018-01-09 | 中国石油大学(华东) | 一种裂缝三维尺寸可控的缝网‐基质双重介质模型及其制备方法 |
-
2018
- 2018-03-01 CN CN201810171477.4A patent/CN108547611B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1378652A (zh) * | 1999-09-15 | 2002-11-06 | 埃克森美孚上游研究公司 | 远程储层电阻率映像 |
US20110231098A1 (en) * | 2009-10-05 | 2011-09-22 | Dzevat Omeragic | Multilevel workflow method to extract resistivity anisotropy data from 3d induction measurements |
CN102635347A (zh) * | 2012-03-30 | 2012-08-15 | 中国电子科技集团公司第二十二研究所 | 一种定量将薄互地层等效成水平和垂直电阻率的方法 |
CN102619504A (zh) * | 2012-04-17 | 2012-08-01 | 中国电子科技集团公司第二十二研究所 | 随钻电磁波电阻率仪器的径向探测深度指标的确定方法 |
CN102678106A (zh) * | 2012-05-02 | 2012-09-19 | 中国电子科技集团公司第二十二研究所 | 随钻电磁波电阻率测井仪器的数据处理方法 |
CN103323883A (zh) * | 2013-06-26 | 2013-09-25 | 中南大学 | 一种近场磁偶源高密度的几何电阻率测深法 |
CN104265279A (zh) * | 2014-07-30 | 2015-01-07 | 中国石油集团川庆钻探工程有限公司 | 断层条件下随钻测井曲线预测方法 |
CN106407574A (zh) * | 2016-09-23 | 2017-02-15 | 上海神开石油设备有限公司 | 一种多分量随钻方位电磁波仪器的快速正反演处理方法 |
CN106446408A (zh) * | 2016-09-23 | 2017-02-22 | 上海神开石油设备有限公司 | 一种随钻补偿电磁波仪器的快速正反演处理方法 |
CN107561244A (zh) * | 2017-04-20 | 2018-01-09 | 中国石油大学(华东) | 一种裂缝三维尺寸可控的缝网‐基质双重介质模型及其制备方法 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110454154B (zh) * | 2019-08-22 | 2020-08-25 | 中国科学院地质与地球物理研究所 | 一种校正井眼偏心对方位信号影响的方法及装置 |
CN110454154A (zh) * | 2019-08-22 | 2019-11-15 | 中国科学院地质与地球物理研究所 | 一种校正井眼偏心对方位信号影响的方法及装置 |
CN110989010A (zh) * | 2019-12-02 | 2020-04-10 | 吉林大学 | 一种基于偏心线圈补偿的瞬变电磁勘探装置及勘探方法 |
CN110989010B (zh) * | 2019-12-02 | 2022-02-15 | 吉林大学 | 一种基于偏心线圈补偿的瞬变电磁勘探装置及勘探方法 |
CN111175834A (zh) * | 2020-01-16 | 2020-05-19 | 吉林大学 | 一种基于全程采集技术的拖曳式瞬变电磁探测装置及方法 |
CN111175834B (zh) * | 2020-01-16 | 2021-06-22 | 吉林大学 | 一种基于全程采集技术的拖曳式瞬变电磁探测装置及方法 |
CN111581791B (zh) * | 2020-04-22 | 2024-05-31 | 中国海洋石油集团有限公司 | 一种电阻率正演仿真方法、装置及计算机可读存储介质 |
CN111581791A (zh) * | 2020-04-22 | 2020-08-25 | 中国海洋石油集团有限公司 | 一种电阻率正演仿真方法、装置及计算机可读存储介质 |
CN113076508A (zh) * | 2021-04-02 | 2021-07-06 | 北京环境特性研究所 | 基于垂直磁偶极子在半空间下的低频近场快速计算方法 |
CN115062474B (zh) * | 2022-06-17 | 2023-01-06 | 上海玫克生储能科技有限公司 | 一种锂电池的工况仿真方法、系统、终端设备及存储介质 |
CN115062474A (zh) * | 2022-06-17 | 2022-09-16 | 上海玫克生储能科技有限公司 | 一种锂电池的工况仿真方法、系统、终端设备及存储介质 |
CN116398126A (zh) * | 2023-06-07 | 2023-07-07 | 山东万洋石油科技有限公司 | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 |
CN116398126B (zh) * | 2023-06-07 | 2023-08-25 | 山东万洋石油科技有限公司 | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108547611B (zh) | 2022-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108547611A (zh) | 水平井复杂环境随钻电磁波电阻率测井快速仿真方法 | |
CN100429531C (zh) | 目标最小化的三维电磁快速反演方法 | |
CN108875122B (zh) | 利用随钻测井数据计算地质参数的人工智能方法和系统 | |
Wilt et al. | Crosshole electromagnetic tomography: A new technology for oil field characterization | |
CN104656157B (zh) | 一种识别页岩气甜点区的方法及装置 | |
CN104360403B (zh) | 一种多辐射场源瞬变电磁法多分量全域视电阻率定义方法 | |
CN104863574B (zh) | 一种适用于致密砂岩储层的流体识别方法 | |
CN111058834B (zh) | 基于瞬变多分量感应测井的各向异性地层倾角确定方法 | |
CN107191181B (zh) | 一种基于电磁散射的井周界面探测方法 | |
CN106407574A (zh) | 一种多分量随钻方位电磁波仪器的快速正反演处理方法 | |
CN108984818A (zh) | 固定翼时间域航空电磁数据拟三维空间约束整体反演方法 | |
CN109597136A (zh) | 一种矿井全空间瞬变电磁数据处理方法 | |
CN108804728B (zh) | 水平井地层储层分级分析方法及计算机可读存储介质 | |
CN105074505A (zh) | 真岩层电阻率的确定 | |
CN108875211A (zh) | 一种基于fetd与fdtd耦合的二维模型的探地雷达正演方法 | |
CN103912270A (zh) | 井间测井深度与速度同步控制定位系统及方法 | |
CN107944159B (zh) | 一种随钻电磁波正演仿真数据库高精度压缩方法 | |
Dutta et al. | New borehole transient electromagnetic system for reservoir monitoring | |
CN108169802A (zh) | 一种粗糙介质模型的时域电磁数据慢扩散成像方法 | |
Wilt et al. | Oil field reservoir characterization and monitoring using electromagnetic geophysical techniques | |
CN103615239A (zh) | 一种测井地层成像系统及其方法 | |
CN108535783B (zh) | 一种对电阻率断面进行层位校正的方法及装置 | |
CN108376204B (zh) | 一种基于地下粗糙介质模型的电磁广义趋肤深度计算方法 | |
Yu et al. | Application of tensor CSAMT with high-power orthogonal signal sources in Jiama porphyry copper deposit, South Tibet | |
Li et al. | Application of the wide field electromagnetic method for oil and gas exploration in a red-bed basin of South China |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |