沈 維，梁華慶，耿 敏，李 靜，竇新宇

（中國石油大學（北京）地球物理與信息工程學院，北京 102249）

雙水平井電磁測距徑向距離計算方法的ANSYS仿真研究

沈 維，梁華慶，耿 敏，李 靜，竇新宇

（中國石油大學（北京）地球物理與信息工程學院，北京 102249）

基于旋轉磁場測距（RMRS）的基本理論，借助ANSYS有限元分析軟件，研究了金屬套管條件下雙水平井電磁測距徑向距離計算方法；首先建立SAGD雙水平井RMRS井下傳播模型，通過設定單元屬性、劃分網(wǎng)格、加載邊界條件等對模型進行求解；然后研究套管厚度、直徑、相對磁導率等對雙水平電磁測距系統(tǒng)中磁場軸向分量的影響；最后利用實驗室現(xiàn)有的旋轉磁場測距模擬裝置對仿真結果進行了驗證；研究結果表明：隨著套管厚度、直徑、相對磁導率的增加，探管接收到的磁感應強度會逐漸減小，但雙水平井徑向間距計算仍可采用均勻介質中的理論測距導向計算方法；此結論可為套管的選取、磁導向儀器的研究及測量資料解釋提供理論參考。

SAGD雙水平井；ANSYS有限元分析；徑向間距計算；金屬套管

0 引言

我國剩余的石油和天然氣儲量大多屬于低品位或難動用資源，其開發(fā)難度越來越大。為了提高采收率，蒸汽輔助重力泄油（steam assisted gravity drainage，簡稱SAGD）技術已經廣泛的應用在超稠油開采中。SAGD雙水平井技術實現(xiàn)的關鍵在于兩口水平井井眼軌跡平行，徑向間距控制誤差不得超過± 0.5 m。為解決SAGD雙水平井井眼軌跡高精度控制的需求，近30年來，國外研發(fā)了一系列用于隨鉆精確監(jiān)控井間距離和方位的工具，主要有磁場定位導向工具MGT（magnetic guidance tool）［1-2］和旋轉磁場導向系統(tǒng)RMRS（rotatingmagnet ranging system）［2］，但其核心技術都被保密和壟斷。

目前，國內也在積極開展這方面的研究。其中以中國石油大學（北京）高德利院士課題組所研制的“鄰井距離隨鉆電磁探測系統(tǒng)”［3-5］、中國石油集團鉆井工程技術研究院的“多分支水平遠距離穿針工具”［6］，和克拉瑪依鉆井工藝研究院的“成對水平井磁定位系統(tǒng)”［7］為代表。SAGD雙水平井正鉆井是裸眼井，而已鉆井是金屬套管井，金屬套管對電磁信號的吸收很大，不同的套管，吸收的情況也不同。本文利用有限元方法，從數(shù)值分析的角度，建立RMRS磁場信號傳播模型，分別研究了金屬套管的厚度、相對磁導率、直徑這幾種非均勻介質因素下雙水平井徑向距離計算的方法，為實際鉆井數(shù)據(jù)的分析與計算提供理論參考。

1 SAGD雙水平井RMRS測距原理與方法

SAGD雙水平井中的RMRS工作原理［5］如圖1所示。RMRS井下儀器主要由探管、磁短節(jié)組成，其中磁短節(jié)安裝在鉆頭后方無磁鉆鋌中，放置在注入井（正鉆井）中，鉆進過程中磁短節(jié)以1～5 Hz左右的頻率隨鉆頭低速旋轉，通過旋轉產生交變磁場，作為發(fā)射磁源。探管由三軸磁場傳感器和三軸加速度傳感器組成，放置在生產井（已鉆井）中，用來探測磁短節(jié)旋轉所產生的交變磁場，通過對探管接收到的磁信號分析計算，可以求出鉆頭的位置和方向，控制鉆頭按預定軌跡鉆進，從而保證兩口井水平段的平行。

圖1 RMRS在SAGD雙水平井中的工作原理

由于探管到磁短節(jié)的距離遠大于磁短節(jié)的尺寸，因此可以將旋轉磁短節(jié)看成一個旋轉的磁偶極子。圖2所示為磁短節(jié)旋轉產生的磁場的磁感應強度計算模型，探管三軸磁場傳感器方向分別為x、y、z，磁短節(jié)與探管之間的位移矢量為r，磁偶極子的磁矩為m，與x軸的夾角為θ。由畢奧－薩伐爾定律可以推導得出磁偶極子遠場的磁感應強度可表示為

文獻［5］中給出了均勻介質中磁感應強度B在x、y、z三軸方向磁感應強度分量表達式的詳細推導過程，本文不再贅述。其中z軸方向磁感應強度為：

即當z＝±1／2r時，磁感應強度軸向分量的振幅達到最大值。因此，Bz幅值兩個最大值之間的距離即為兩口水平井水平段距離。

圖2 磁偶極子磁感應強度計算模型

2 RMRS系統(tǒng)物理模型

以文獻［5］中的實驗場景為例，設雙水平井水平段的間距為D時，將探管下入到生產井合適位置，探管到鉆頭的軸向距離大約為D，鉆頭繼續(xù)鉆進大約2D的距離，記錄在這段距離內探管接收到的磁信號，利用探管接收磁短節(jié)產生的磁信號確定SAGD雙水平井水平段空間位置。以此實驗為基礎，利用ANSYS有限元分析軟件全真還原探管和磁短節(jié)的運行軌跡。如圖3所示，現(xiàn)場實驗中鉆頭相對于探管鉆進2D的距離等效于保持鉆頭固定不動，傳感器沿著Z方向移動2D的距離，然后記錄這段距離的磁信號。根據(jù)RMRS的測距原理和方法來計算SAGD雙水平井的徑向距離。研究不同的金屬套管參數(shù)下水平井徑向距離計算的方法。

圖3 ANSYS仿真模型

3 基于有限元的數(shù)值模擬

以電磁場理論為基礎，依據(jù)有限元理論和物理模型進行網(wǎng)格劃分［8-9］，目前ANSYS三維仿真能支持四面體單元，此單元的形狀簡單、穩(wěn)定性良好。劃分方法則是根據(jù)單元的各自情況進行的自適應網(wǎng)格劃分法，并通過指定程序中有限元網(wǎng)格單元的尺寸大小，使誤差降至最低，如圖4所示為磁短節(jié)的網(wǎng)格模型。數(shù)值模擬中，永磁體材料是以釹鐵硼為代表，其電導率為6.25×10－5S／m，相對磁導率為1.35，矯頑力為50 000 N，直徑為20 mm，高為50 mm。套管厚度為8 mm，相對磁導率為100，直徑為178 mm，電導率為2×106S／m。雙水平井之間的距離以4 m為最佳，整個模型區(qū)域為20 m*20 m*20 m。由式（4）可知，計算SAGD雙水平井間距時，需獲得磁短節(jié)正交磁場的軸向分量。根據(jù)探管和磁短節(jié)的相對運動軌跡，利用有限元分析導出的ANSYS仿真模型中的y方向正4 m處磁場強度軸向分量Bz的數(shù)據(jù)，即可獲得雙水平井的間距。

圖4 磁短節(jié)的網(wǎng)格模型

3.1 套管厚度的影響

根據(jù)石油套管規(guī)格表，設置標準的套管參數(shù)。以圖3模型為基礎，設地層相對磁導率為1、套管直徑為178 mm、套管相對磁導率100，分別建立套管厚度為5.6 mm，6 mm，7 mm，8 mm，9 mm，10 mm，11.4 mm，12.4 mm，13.4 mm，15.4 mm，17.4 mm，19 mm，20.4 mm，22 mm時的RMRS井下磁場傳播的ANSYS模型。

選取其中4組代表性的套管厚度參數(shù)，仿真得到距離磁源正4 m處的Bz的變化曲線如圖5所示。其變化趨勢和均勻介質的磁場信號的軸向分量的傳播規(guī)律一樣，在8 m和12 m處達到幅值，詳細數(shù)據(jù)列于表1中。由于趨膚效應［10］的影響，隨著厚度的增加，使得金屬套管中的磁場信號衰減嚴重，Bz的值逐漸減小。而Bz兩個幅值的水平位置仍然保持在8 m和12 m處，根據(jù)均勻介質中的理論計算公式，兩者之差即為雙水平井的徑向間距，由表1可知，其與真實值4 m之間的相對誤差均小于3%。這說明在套管厚度的變化的條件下，計算雙水平井的徑向間距仍可采用理論測距導向計算方法。

表1 金屬套管的厚度對測量結果的影響

圖5 不同厚度下距磁源4 m處Bz的變化曲線

3.2 套管直徑的影響

在圖3模型中，固定參數(shù)為地層相對磁導率為1、套管相對磁導率為100、套管厚度為8 mm，不定參數(shù)為套管直徑。根據(jù)石油套管規(guī)格，分別設置套管直徑為126 mm，140 mm，168 mm，178mm，194mm，220mm，224mm，244mm，273 mm，298 mm，340 mm，逐個建立RMRS井下磁場傳播的ANSYS模型。

通過網(wǎng)格劃分、加載、后處理等步驟提取距離磁短節(jié)正4 m處的磁感應強度。選取其中3組套管直徑參數(shù)繪制Bz曲線如圖6所示，Bz曲線的變化趨勢基本上一樣。導出Bz達到幅值的水平位置，詳細數(shù)據(jù)列于表2中。隨著直徑的增大，Bz的值減小，而Bz達到幅值時的水平位置仍處于8 m和12 m，其差值與實際值4 m的相對誤差基本上在3.00%左右，所以此數(shù)據(jù)驗證了不同直徑下，SAGD雙水平井水平段的間距等于兩個磁場強度軸向分量Bz取得最大值處之間的距離。

表2 金屬套管的直徑對測量結果的影響

圖6 不同直徑下距磁源4m處Bz的變化曲線

3.3 套管相對磁導率的影響

在地層相對磁導率為1、套管直徑為178 mm、套管厚度為8 mm的條件下，由于金屬套管的相對磁導率的范圍為40～110，則分別建立套管相對磁導率為40、50、60、90、100、110時的RMRS井下磁場傳播的ANSYS模型。

通過ANSYS數(shù)據(jù)分析，選取其中3組相對磁導率參數(shù)得到Bz隨金屬相對磁導率的變化曲線如圖7所示。此圖中表明磁場信號的軸向分量分布曲線在不同的相對磁導率條件下大致相同。詳細數(shù)據(jù)如表3所示，Bz隨著套管相對磁導率的增大而減小，這是由于金屬套管中趨膚效應造成的。而其極幅值之間的間距與真實值4 m的相對誤差在0.5%以內，從而驗證了不同相對磁導率下，雙水平井水平段間距等于Bz兩個振幅最大值之間的間距。

表3 金屬套管的相對磁導率對測量結果的影響

4 室內模型實驗

4.1 空氣中模擬實驗

利用實驗室現(xiàn)有的井下微弱交變磁場信號的采集與處理模擬裝置，在地面進行模擬實驗，磁短節(jié)與探管之間徑向距離為4.45 m，探管采集到的數(shù)據(jù)被傳輸?shù)焦P記本電腦中，通過軟件系統(tǒng)對數(shù)據(jù)進行處理，信號處理方法為DFT雙譜峰的頻率重構算法，處理后的數(shù)據(jù)再通過幅值計算計算出電壓值，也就是在各個測量點接收到的電壓值，電壓值的大小與磁場強度大小對應，最后將各個測量點的三軸方向的電壓值進行曲線擬合。圖8是根據(jù)信號Bz的幅值畫出的信號幅值隨距離的變化曲線，兩個峰值位置分別為－2.3 m、2.3 m處，因此通過探管測量的磁短節(jié)與探管之間的距離為4.60 m，與真實值的誤差為0.05 m，相對誤差為1%。驗證了公式（4）推導的結論。

圖8 實測z方向響應曲線

4.2 套管中的模擬實驗

將測量探管放入直徑為178 mm，管壁厚度為12.4 mm的套管中，磁短節(jié)與探管之間的徑向距離為3.0 m，采集到Bz信號的峰峰值列于表4中，圖9是根據(jù)表4數(shù)據(jù)畫出的信號幅值隨距離的變化曲線。

表4 信號幅值隨距離的關系表

由圖9可以看出，接收信號強度最強的兩個點的之間的間隔距離為3.0 m，剛好等于磁短節(jié)與探管之間的徑向距離。該實驗結果驗證了本文的仿真結論，即金屬套管僅對磁場信號的大小產生影響，不會對雙水平井間徑向距離的計算產生影響。

圖9 信號幅值隨距離的變化曲線

5 結論

1）根據(jù)有限元理論建立了RMRS磁場信號井下傳播模型，通過數(shù)值計算表明金屬套管對探管接收的磁信號影響較大，金屬套管的厚度、直徑越大，套管對磁信號的屏蔽作用越大，金屬套管的磁導率越大，探管接收的磁信號強度就越小。為了減小金屬套管對磁性導向儀器的影響，應盡量選取套管直徑和厚度適當、非磁性或弱磁性套管。

2）以理論測距導向計算方法為基礎，研究了套管厚度、直徑、相對磁導率等因素條件下雙水平井徑向距離的計算方法。通過數(shù)值解方式驗證了在不同金屬套管參數(shù)下，雙水平井徑向間距的計算仍等于探管記錄的軸向磁信號的兩個幅值之間的距離。

3）室內模擬實驗驗證了雙水平井間徑向距離的計算，進一步需要研究不同金屬參數(shù)對計算雙水平井方位角的影響。

4）上述結論可為套管的選取、磁性導向工具的研制和鉆井資料解釋提供理論參考。

［1］Grills T L.Magnetic ranging technologies for drilling steam assisted gravity drainage well pairs and unique well geometries－a comparison of technologies［R］.SPE 79005，2002.

［2］Kuckes A F，Hay R T，Mcmahon N J，et al.New electromagnetic surveying／ranging method for drilling parallel horizontal twin wells ［R］.SPE 27466，1996.

［3］閆永維，高德利，吳志永.煤層氣連通井引導技術研究［J］.石油鉆采工藝，2010，32（2）：23-25.

［4］王德桂，高德利.管柱形磁源空間磁場矢量引導系統(tǒng)研究［J］.石油學報，2008，29（4）：608-611.

［5］刁斌斌，高德利，吳志永.雙水平井導向鉆井磁測距計算方法［J］.中國石油大學學報（自然科學版），2011，35（6）：71-75.

［6］田中蘭，喬磊，蘇義腦.鄭平01-1煤層氣多分支水平井優(yōu)化設計與實踐［J］.石油鉆采工藝，2010，32（2）：20-22.

［7］梁華慶，耿 敏，時東海.旋轉磁場井間隨鉆測距導向系統(tǒng)中微弱頻變信號的檢測方法［J］.中國石油大學學報（自然科學版），2013，37（4）：83-87.

［8］宋汐瑾，黨瑞榮，董 昭.套管井電磁測井技術及套管規(guī)律研究［J］.測井技術，2010，34（2）：143-145.

［9］魏寶君，張庚驥，梁秋錦.金屬套管對電磁測井響應的影響［J］.石油物探，1999，38（3）：67-75.

［10］呼石磊，鄢泰寧，李 曉.地層對電磁隨鉆測量信號的影響研究［J］.煤炭科學技術，2011，39（9）：114-117.

［11］李 金，鄭小林，侯文生.用于磁定位的永磁體磁場仿真與實驗研究［J］.系統(tǒng)仿真學報，2009，21（18）：5919-5922.

Double Horizontal Electromagnetic Radial Distance Calculation Method of ANSYSSimulation

Shen Wei ，Liang Huaqing，Geng Min，Li Jing，Dou Xinyu

（College of Geophysics and Information Engineering，China University of Petroleum－Beijing，Beijing 102249，China）

Based on the basic theory of rotating magnetic field ranging（RMRS），with the aid of finite element analysis software ANSYS，research influence of metal casing on Double horizontal wells electromagnetic ranging radial distance calculation.Initially，SAGD double horizontal well RMRS downhole transmission model is set up，solving the model by setting unit properties，division of grid，load boundary conditions.Then research the influence of casing thickness，diameter，and relative permeability on axial component of magnetic field of double horizontal electromagnetic distance measuring system.The results show that：with the increasing casing thickness，diameter，relative permeability，the magnetic induction probe received will gradually decrease，but the dual horizontal wells radial spacing calculated use still theoretical ranging guide calculation methods of homogeneous medium.This conclusion can be applied in selection of sleeve，the magnetic orientation of the instrument research and measurement data interpretation providing a theoretical reference.

SAGD double horizontal well；finite element；inter－well distance calculation；metal casing

1671-4598（2016）08-0248-04

10.16526／j.cnki.11－4762／tp.2016.08.068

：TE243

：A

2016-02-26；

：2016-04-08。

國家科技重大專項（2011ZX05009－005）；中國石油大學（北京）基金（KYJJ2012－05－33）。

沈 維（1991-），女，湖南岳陽人，碩士研究生，主要從事信號監(jiān)測與石油儀器方向的研究。