亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        傾斜各向異性地層隨鉆電磁波測井響應(yīng)模擬

        2013-12-03 10:30:06范宜仁胡云云李虎孫慶濤
        測井技術(shù) 2013年2期
        關(guān)鍵詞:磁偶極子偶極子斜角

        范宜仁,胡云云,李虎,孫慶濤

        (1.中國石油大學地球科學與技術(shù)學院,青島 山東266555;2.中國石油大學CNPC測井重點實驗室,青島 山東266555)

        0 引 言

        斯倫貝謝、哈里伯頓以及貝克休斯等公司相繼推出的隨鉆電磁波測井儀器[1-2]能較好地實現(xiàn)實時隨鉆地質(zhì)導(dǎo)向與儲層測井評價功能。隨鉆電磁波測井響應(yīng)理論大多仍基于各向同性地層模型,實際隨鉆測井過程中各向異性是影響儀器測井響應(yīng)的一個重要因素[3],尤其是在裂縫發(fā)育或者薄互層發(fā)育的層段,地層所呈現(xiàn)出的宏觀各向異性會使隨鉆電磁波測井響應(yīng)嚴重失真,不能有效反映地層信息[4],給測井儲層評價帶來了很大的困難。沈金松[5]、王昌學和周燦燦等[6]采用有限差分方法模擬各向異性地層隨鉆電磁波測井響應(yīng),但三維數(shù)值模擬計算速度過慢,不適應(yīng)現(xiàn)場應(yīng)用。因此,有必要開展基于解析方法的各向異性地層隨鉆電磁波測井響應(yīng)快速模擬的研究,為各向異性地層隨鉆電磁波響應(yīng)特征研究以及基于隨鉆電磁波測井的實時地質(zhì)導(dǎo)向提供支持。本文針對傾斜各向異性地層,進行了基于解析方法的隨鉆電磁波測井響應(yīng)數(shù)值模擬,并針對地層各向異性系數(shù)、儀器發(fā)射頻率、井斜角等對隨鉆電磁波測井響應(yīng)的影響進行了理論分析。

        1 原 理

        圖1為各向異性地層模型。其中,T為發(fā)射線圈,R1、R2分別為近接收線圈和遠接收線圈;ε1、ε2、ε3為各層的介電常數(shù);μ為磁導(dǎo)率;σ1h、σ2h、σ3h、σ1v、σ2v、σ3v分別為各層的水平電導(dǎo)率和垂直電導(dǎo)率;儀器軸與z軸夾角為α;接收線圈與發(fā)射線圈法向均沿儀器軸方向。

        電導(dǎo)率張量的Maxwell方程[7]可以表示為

        圖1 各向異性地層模型

        其中,σ′為復(fù)電導(dǎo)率張量

        定義赫茲矢量位函數(shù)

        其洛倫茲條件為Δ

        對于層狀介質(zhì)中任意方向的磁偶極子源,可以分解為垂直磁偶極子源(VMD)和水平磁偶極子源(HMD)等2個部分[8],通過赫茲矢量位滿足的波動方程分別描述VMD和HMD產(chǎn)生的場,進而基于模型將矢量方程轉(zhuǎn)化為標量方程,再由邊界條件確定各層電磁場的解析解。根據(jù)法拉第電磁感應(yīng)定律,接收線圈處的電動勢VR為

        式中,ω為發(fā)射線圈頻率;s為接收線圈面積;Hx為接收線圈處總的水平磁場強度;Hz為接收線圈處總的垂直磁場強度。

        隨鉆電磁波測井儀器線圈系的設(shè)計均以單發(fā)雙收的三線圈系為基礎(chǔ)結(jié)構(gòu),通過測量2個接收線圈電動勢間的相位差和幅度比轉(zhuǎn)化為相位差電阻率與幅度比電阻率

        式中,VR1、VR2分別為近接收線圈與遠接收線圈的電動勢;arg表示取相位角;Pd為2個接收線圈電位的相位差;Ar為2個接收線圈電位的幅度比。其中,垂直磁偶極子源位函數(shù)為

        其產(chǎn)生的磁場為

        根據(jù)邊界條件,式(12)、式(13)即可確定各層系數(shù)Fi、Gi,進而確定垂直偶極源子產(chǎn)生的磁場分布。

        水平磁偶極子源可以分解為TE波與TM波,其中TE波位函數(shù)

        TM波位函數(shù)

        產(chǎn)生的磁場為

        根據(jù)邊界條件,式(18)至式(21)即可確定各層系數(shù)Pi、Qi、Ti、Si,進而求取水平偶極子源產(chǎn)生的磁場強度。

        2 含Bessel函數(shù)廣義積分的求解

        2.1 被積函數(shù)的收斂性

        垂直偶極子源與水平偶極子源磁場[式(10)、式(11)、式(16)、式(17)]的計算均涉及到 Bessel函數(shù)的無窮積分;由式(6),接收線圈總的磁場強度為H=cosαHx+sinαHz。圖2給出了井斜角為0°、45°、75°、80°和85°時H中被積函數(shù)實部、虛部與積分變量λ的關(guān)系;在低角度時被積函數(shù)的震蕩性較弱,但隨著角度增加,被積函數(shù)的震蕩性明顯增強,夾角小于80°時,積分變量λ達到[0,80]即能滿足工程需要;夾角大于80°時,積分變量λ達到[0,200]才能滿足工程精度要求。當井斜角為90°時(水平偶極子),被積函數(shù)變?yōu)橥耆l(fā)散(見圖3)。

        井斜角小于80°或大于100°時,為了提高積分精度,采用高斯積分[10-11]進行求解;井斜角大于80°且小于100°時,由于積分區(qū)間較大,為了加快計算速度,采用Hankel變換[12]對含Bessel函數(shù)的廣義積分進行求解。

        圖2 被積函數(shù)與λ的關(guān)系

        2.2 算法驗證

        為了驗證算法的正確性,針對 Myung Jin Nam[13]給出的線圈結(jié)構(gòu)分別模擬了各向異性地層不同磁場分量Hxx(接收線圈與發(fā)射線圈法向皆沿x軸方向)與Hzz(接收線圈與發(fā)射線圈法向皆沿z軸方向)響應(yīng)(結(jié)果見圖4和圖5)。模擬結(jié)果與K.H.LeeEM1D的模擬結(jié)果[13]能很好地吻合。

        3 數(shù)值模擬結(jié)果與分析

        3.1 各向異性地層隨鉆電磁波模擬測井響應(yīng)特征

        圖6、圖7分別模擬了各向同性與各向異性三層地層模型隨鉆電磁波測井響應(yīng),其中2個接收線圈的線圈距分別為24in和30in,電阻率各向異性系數(shù)定義為各向同性地層中,目的層中離地層界面較遠處的視電阻率值與實際地層電阻率值接近,且基本不受井斜角的影響(見圖6)。而在各向異性地層中,不同井斜角的視電阻率發(fā)生分異,井斜角較小時(井斜角小于30°)視電阻率值基本反應(yīng)實際水平電阻率,隨著井斜角的增大,視電阻率值逐漸向垂直電阻率過渡,且相位差電阻率較幅度比電阻率對垂直電阻率更加敏感;在大斜度井中(井斜角大于70°),由于垂直電阻率的影響,視電阻率曲線在地層界面處產(chǎn)生“犄角”明顯大于各向同性地層(見圖7)。

        3.2 不同各向異性系數(shù)隨鉆電磁波測井響應(yīng)特征

        井斜角較小時(井斜角小于30°)各向異性對視電阻率的影響較小,視電阻率值基本反映實際水平電阻率;井斜角較大時(井斜角大于30°)各向異性(η>1時)對視電阻率的影響明顯,視電阻率值逐漸靠近垂直電阻率,且當各向異性系數(shù)η>3時大斜度井中(井斜角大于80°)相位差電阻率超過垂直電阻率(見圖8)。

        3.3 發(fā)射頻率對隨鉆電磁波測井響應(yīng)的影響

        發(fā)射頻率越大,視電阻率對地層各向異性越敏感,且相位差電阻率較幅度比電阻率更加敏感。發(fā)射頻率小于100kHz時,不同頻率的幅度比電阻率值幾乎一致,而相位差電阻率略有差異;發(fā)射頻率大于100kHz時,不同頻率的幅度比電阻率在井斜角大于50°時才出現(xiàn)差異,而相位差電阻率則在井斜角為30°時即出現(xiàn)較大差異;當發(fā)射頻率為2MHz時,大斜度井(井斜角大于75°)相位差電阻率已超過垂直電阻率(見圖9)。

        圖9 不同頻率條件下隨鉆電磁波測井響應(yīng)隨井斜角的變化

        4 結(jié) 論

        (1)將各向異性地層任意方向的磁偶極子源分解為水平方向與垂直方向磁偶極子源的疊加,采用赫茲矢量位函數(shù)對偶極子源產(chǎn)生的磁場進行了描述,基于高斯積分與Hankel變換,分別對不同角度Bessel函數(shù)進行求解,克服了被積函數(shù)的強震蕩性。

        (2)井斜角較小時,隨鉆電磁波測井視電阻率主要反映水平電阻率;隨井斜角增大,視電阻率逐漸由水平電阻率向垂直電阻率過渡,且相位電阻率較幅度比電阻率更為敏感;當?shù)貙痈飨虍愋韵禂?shù)與井斜角較大時,相位差電阻率甚至會超過垂直電阻率。

        (3)低頻條件下(發(fā)射頻率小于100kHz)隨鉆電磁波視電阻率受各向異性影響較小,不同發(fā)射頻率下的幅度比電阻率幾乎重合,而相位差電阻率略有差異;發(fā)射頻率越高,視電阻率受地層各向異性影響越大,高頻條件下視電阻率將遠大于低頻條件。

        [1] Zhang Z,Gonguet C,Rajani V.Directional LWD Resistivity Tools and Their Business Impacts[C]∥ SPWLA 49th Annual Logging Symposium,Texas:Austin,25-28May,2008.

        [2] Davydycheva S.3DModeling of New-generation (1999-2010)Resistivity Logging Tools,Society of Exploration[J].Geophysicists,2010,29(7):780-789.

        [3] Bittar M,Rodney P.The Effect of Rock Anisotropy on MWD Electromagnetic Wave Resistivity Sensors[J].The Log Analyst,1996,37(1):20-30.

        [4] Abubakar A,Habashy T M,Druskin V,et al.A 3D Parametric Inversion Algorithm for Triaxial Induction Data[J].Geophysics,2006,71(1):G1-G9.

        [5] 沈金松.用有限差分法計算各向異性介質(zhì)中多分量感應(yīng)測井的響應(yīng)[J].地球物理學進展,2004,19(1):101-107.

        [6] 王昌學,周燦燦,儲昭坦,等.電性各向異性地層頻率域電磁響應(yīng)模擬[J].地球物理學報,2006,49(6):1874-1883.

        [7] Li J,Bhardwaj A,et al.Computation of Triaxial Induction Logging Tools in Layered Anisotropic Dipping Formations[C]∥IEEE,2008,46(4):1149-1163.

        [8] 徐建華.層狀媒質(zhì)中的電磁場與電磁波[M].北京:石油工業(yè)出版社,1997:172-183.

        [9] 其木蘇榮,孫丙西.水平層狀介質(zhì)中任意方向磁偶極子的電磁場分布[J].內(nèi)蒙古民族大學學報:自然科學版,2004,19(3):248-256.

        [10] 莫平華.一階貝塞爾函數(shù)廣義積分的數(shù)值計算[J].數(shù)學理論與應(yīng)用,2007,27(1):66-67.

        [11] 張輝,李桐林,董瑞霞,等.利用高斯求積和連分式展開計算電磁張量格林函數(shù)積分[J].地球物理學進展,2005,20(3):667-670.

        [12] Walter L A.Computer Program Numerical Integration of Related Hankel Transforms of Oders 0and 1by Adaptive Digital Filtering [J].Geophysics,1979,44(7):1287-1305.

        [13] Nam M J,Pardo D,Torres-Verdín C.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 [J].Geophysics,2010,75(3):F83-F95.

        猜你喜歡
        磁偶極子偶極子斜角
        基于DDS的正交偶極子聲波測井儀快檢裝置研究
        化工管理(2021年7期)2021-05-13 00:46:12
        基于工具造斜率和井斜角的井眼軌跡設(shè)計
        云南化工(2020年11期)2021-01-14 00:51:02
        以“傾斜角與斜角”為例談概念教學
        弧形寬帶印刷偶極子5G天線的設(shè)計
        基于遞推更新卡爾曼濾波的磁偶極子目標跟蹤
        磁偶極子跟蹤的漸進貝葉斯濾波方法
        自動化學報(2017年5期)2017-05-14 06:20:53
        建立平面斜角坐標系在向量中的應(yīng)用
        基于磁偶極子的磁場梯度張量縮并的試驗驗證及相關(guān)參數(shù)確定
        水平磁偶極子電磁場特征研究
        一種新的無源偶極子天線輻射效率測量方法
        亚洲av精品一区二区三区| 精品少妇一区二区三区四区| 亚洲三级中文字幕乱码| 中文字幕人妻熟女人妻| 婷婷成人基地| 成人国产精品高清在线观看| 国产午夜精品视频观看| 久久国产成人精品国产成人亚洲| 免费看久久妇女高潮a| 国产亚洲精品hd网站| 日本韩国三级在线观看| 人妻少妇久久中文字幕| 无码专区天天躁天天躁在线| 精品一区二区三区在线观看l| 久久综合激情的五月天| 精品国产免费一区二区三区 | 日本看片一区二区三区| 变态另类人妖一区二区三区 | 久草国产手机视频在线观看| 国产av一区二区三区在线播放| 玩50岁四川熟女大白屁股直播| 456亚洲老头视频| 91久久国产精品综合| 亚洲国产精品久久艾草| 亚洲日韩精品欧美一区二区 | 亚洲国产天堂av成人在线播放| 老太婆性杂交视频| 精品福利视频一区二区三区 | 天天干夜夜操| 无码AⅤ最新av无码专区| 亚洲精品在线免费视频| 高清不卡一区二区三区| 国产成人九九精品二区三区| 免费黄网站一区二区三区| 丰满大爆乳波霸奶| 国产精品亚洲日韩欧美色窝窝色欲 | 久久久无码人妻精品一区| 青青国产成人久久91| 亚洲精品中文字幕91| 国内女人喷潮完整视频| 久草视频这里有精品|