祝 鵬，林承焰，李智強，趙文積，張華蓮

1.中國石油大學(華東)地球科學與技術學院，山東 青島 266580 2.中國電波傳播研究所，河南 新鄉(xiāng) 453003 3.大慶鉆探工程公司測井公司，黑龍江 大慶 163412 4.重慶地質礦產研究院， 重慶 400042

?

水平井和大斜度井中陣列側向測井響應數值模擬

祝 鵬1，林承焰1，李智強2，趙文積3，張華蓮4

1.中國石油大學(華東)地球科學與技術學院，山東 青島 266580 2.中國電波傳播研究所，河南 新鄉(xiāng) 453003 3.大慶鉆探工程公司測井公司，黑龍江 大慶 163412 4.重慶地質礦產研究院， 重慶 400042

針對水平井和大斜度井中陣列側向電極系的工作原理，利用多電場疊加方式進行電場合成，采用三維有限元方法模擬仿真各個分電場的場分布，進而利用電場線性疊加原理得到陣列側向測井響應。在基于計算機仿真的基礎上，得到陣列側向五條測井曲線的徑向探測深度，陣列側向徑向探測深度要小于深側向探測深度?？疾炝巳S地層模型下井斜和侵入深度變化對陣列側向測井響應的影響，分析了水平井和大斜度井中陣列側向測井響應特征。模擬結果表明，在井斜小于15°時，陣列側向測井響應受井斜影響小，可以不進行井斜校正；井斜超過60°的大斜度井以及水平井中，陣列側向測井響應視地層厚度逐漸增大，測井響應值與直井條件下響應值差別較大，必須進行井斜校正。

水平井/大斜度井；陣列側向；有限元；測井響應

0 引言

水平井/大斜度井技術是油氣勘探開發(fā)的重要技術，已經成為各大石油公司的重大戰(zhàn)略措施之一。在水平井/大斜度井中，地層不具有旋轉軸對稱性，不能將問題簡化到子午面上；原有混合法、解析方法等已經不適用于三維地層中，三維有限元方法可以有效地求取三維地層中的陣列側向儀器響應[1-5]。對于陣列側向測井而言，水平井/大斜度井中陣列側向測井響應特征與直井下的陣列側向測井響應差異較大。因此，對陣列側向測井儀進行數值模擬不僅可以認識儀器的探測特征，為儀器研發(fā)提供一定的理論依據,而且可以為測井分析者更加清楚地認識測井曲線提供理論依據。筆者利用三維有限元方法計算陣列側向測井響應，并對水平井/大斜度井中的陣列側向測井響應特征進行分析。

1 陣列側向電極系結構

斯倫貝謝公司于1998年推出高分辨率新型陣列側向測井儀(圖1)。儀器中間是主電流發(fā)射電極A0，兩側分別布置12對對稱電極。在A0和屏蔽電極A1、A’1之間，有2對監(jiān)督電極M1、M’1，M2、M’2；在屏蔽電極A1、A’1和A2、A’2之間，有2對監(jiān)督電極M3、M’3，M4、M’4；在屏蔽電極A2、A’2和A3、A’3之間，有2對監(jiān)督電極M5、M’5，M6、M’6；向外依次為A4、A’4，A5、A’5，A6、A’6。同命名號的電極等電位，例如M1與M’1等電位，以下不再贅述。

當A0發(fā)射電流，其他屏蔽電極為回路電極時，測量得到泥漿電阻率RA0；從A0向兩側每增加1對屏蔽電極為發(fā)射電極時，得到探測深度不同的響應RA1、RA2、RA3、RA4、RA5。陣列側向可以測量得到6條不同探測深度的電阻率曲線，其中探測最淺的主要反映泥餅電阻率的影響[6-8]。

圖1中給出了儀器與地層之間的關系。儀器與地層之間的夾角為θ，井眼大小為dh，泥漿電阻率為Rm，上下圍巖電阻率為Rs，儲層沖洗帶電阻率為Rxo，目的層電阻率為Rt。

圖1 陣列側向儀器及地層模型Fig.1 Array lateral instruments and stratigraphic model

2 陣列側向測井響應三維數值計算方法

2.1 泛函構造

直流電測井采用低頻交流電作為供電電源，由于其頻率較低，陣列側向測井響應可以歸結為穩(wěn)定電流場計算。陣列側向測井的電場可由微分方程表示[9-13]：

(1)

式中：σ為地層電導率；U為測量電位。在柱坐標系(r,φ,z)下，可以表示為

(2)

為求取方程(2)，需將偏微分方程問題轉化成求泛函的極值問題[14-15]：

(3)

(4)

(5)

式中：IE，UE是各電極上電流和電位；E為電極的邊界；Ω為所求取地層的邊界。通過有限元方法求得電位后，陣列側向測井響應值滿足

(6)

式中：K為儀器電極系常數；UM1為M1電極上的電位值；IA0為儀器A0電極發(fā)射電流強度；Ra為不同條件下的視電阻率。

2.2 電場疊加

在求解陣列側向電位分布時，采用電場疊加原理，將其化為多個相應場的電位函數迭加。

第一個分場只有A0電極發(fā)射單位電流I=1，其他電極均不發(fā)射電流；從第二個分場開始，第n個分場只有An-1、An-1’電極發(fā)射單位電流I=1，其他電極均不發(fā)射電流，An-1、An-1’電極等電位(2≤n≤7)。依次計算得到第一到第七分場電位U0(x,y,z)、U1(x,y,z)、U2(x,y,z)、U3(x,y,z)、U4(x,y,z)、U5(x,y,z)、U6(x,y,z)，則陣列側向電位為

U=U0+C1U1+C2U2+

(7)

式中：C1、C2、C3、C4、C5、C6為聚焦合成系數。

下面以RA1為例對陣列側向工作模式的約束條件進行論述，其約束條件為

(8)

(9)

(10)

式中：UM1，UM2為電極M1、M2的電位值；UA2，UA3，UA4，UA5，UA6為電極A2、A3、A4、A5、A6的電位值；IA0，IA1，IA2，IA3，IA4，IA5，IA6為電極A1、A2、A3、A4、A5、A6的電流值。

利用式(8)--(10)的約束條件，可以求出各分場的聚焦合成系數：

(11)

(12)

(13)

(14)

(15)

(16)

(17)

式中：UiM1、UiM2、UiA2、UiA3、UiA4、UiA5、UiA6為第i--1分場下各個電極的電位值。

將公式(6)代入公式(11)中可以得到視電阻率RA1的電阻率值，通過不同的聚焦方法可以得到RA2、RA3、RA4、RA5的電阻率值。

2.3 均勻地層電場分布

模擬條件：均勻地層電阻率為1 Ω·m，儀器半徑為0.045 m。由模擬結果(圖2)可知，隨著屏蔽電極數的增多，在相同的位置處地層的電位值逐漸增大，儀器的探測范圍也逐漸增大。

a.RA5響應；b.RA4響應；c.RA3響應；d.RA2響應；d.RA1響應。圖2 均勻地層下不同測量方式的電位分布Fig.2 Electric potential distribution of different measurement of homogeneous formation

3 水平井/大斜度井中陣列側向測井響應

3.1 偽幾何因子計算

根據表1所示的地層模型，井徑dh=20.32 cm，計算陣列側向測井的偽幾何因子，并確定其探測深度。偽幾何因子為

(18)

以50%偽幾何因子來定義儀器的探測深度，如圖3所示。陣列側向測井在高侵模擬條件下5種探測方式的探測深度分別為：RA50.70 m，RA40.45 m，RA30.37 m，RA20.29 m，RA10.23 m。陣列側向測井在低侵模擬條件下5種探測方式的探測深度分別為：RA50.75 m，RA40.54 m，RA30.45 m，RA20.38 m，RA10.30 m。

表1 地層模型

注：Rt為目的層電阻率；Rxo為侵入帶電阻率；Rm為泥漿電阻率。

圖3 低侵地層(a，地層1)和高侵地層(b，地層2)陣列側向偽幾何因子Fig.3 Array lateral pseudo geometrical factor of low invasion formation (a,Formation 1) and high invasion formation (b,Formation 2)

陣列側向測井探測深度要小于斯倫貝謝CSU(cyber service unit)深側向測井的探測深度(1.27 m)，但是陣列側向測井優(yōu)勢在于反映地層電阻率徑向的漸變，可以在徑向進行電阻率剖面成像。

3.2 井斜影響

根據表2所示的地層模型，模擬兩種地層模型在不同井斜條件下的陣列側向測井響應。井斜依次為0°、15°、30°、45°、60°、75°、90°，模擬結果見圖4、圖5。

表2 地層模型

注：H為目的層厚度；Rb為圍巖電阻率。

圖4 厚1 m的高電阻率地層(地層3)陣列側向測井響應Fig.4 Array lateral log response of 1 m thick and high resistivity formation(Formation 3)

圖5 厚1 m的低電阻率地層(地層4)陣列側向測井響應Fig.5 Array lateral log response of 1 m thick and low resistivity formation(Formation 4)

圖6 厚度1 m的低侵地層(地層5)陣列側向測井響應Fig.6 Array lateral log response of 1 m thick and low invasion formation(Formation 5)

圖7 厚度1 m的高侵地層(地層6)陣列側向測井響應Fig.7 Array lateral log response of 1 m thick and high invasion formation(Formation 6)

從圖4中可以看出，對于1 m厚的地層，在井斜為0°時為直井，由于陣列側向測井具有分辨率高的優(yōu)勢(分辨率為0.4 m)，可以很好地反映地層信息，探測深度最深的RA5響應值為18.1 Ω·m。隨著儀器探測深度的降低，探測深度淺的儀器受井眼和圍巖的影響較大，從而使RA4、RA3、RA2、RA1測量值依次下降。當井斜小于15°時，井斜對陣列側向測井響應的影響可以忽略；隨著井斜角度的逐漸增大，地層中點的陣列側向電阻率明顯降低；在井斜超過60°時，陣列側向測井響應值下降為直井中的陣列側向響應的70%，運用此時的電阻率信息代入到阿爾奇公式計算含水飽和度，必定會帶來極大的誤差。井斜的增大對探測深度最深的RA5影響最大，對探測深度最淺的RA1影響最小。

從圖5中可以看出，對于低電阻率泥巖層，當地層厚度為1 m、井斜為0°時，探測深度最深的RA5響應曲線測量值接近目的層的電阻率，而隨著探測深度的減小，RA4、RA3、RA2、RA1測量值受圍巖影響較大，響應值依次增大。當井斜小于15°時，井斜對于陣列側向響應的影響可以忽略；隨著井斜角度的增大，地層中點的陣列側向電阻率也逐漸增大，陣列側向測井響應受圍巖的影響變大；井斜同樣對于探測深度最大的RA5影響最大，對于探測深度最小的RA1影響較小。

3.3 侵入影響

對于斜井中，測井資料校正一般僅是對地層視厚度進行校正，對于井斜與侵入縱向影響研究較少。下面主要考慮不同傾斜井中陣列側向測井響應值和侵入深度及井斜角的關系，地層模型如表3所示。

井斜依次為0°、15°、30°、45°、60°、75°、90°，侵入半徑0～1.9 m，地層5、6模擬結果如圖6、圖7所示。

表3 地層模型

由圖6可見，對于低侵地層，在井斜小于15°范圍以內，即使存在侵入，井斜角的影響可以忽略。隨著侵入深度的增大，視電阻率與地層電阻率的比值逐漸降低；在同樣的侵入深度，隨著井斜角的增大，視電阻率與地層電阻率的比值同樣出現單調下降的趨勢。在低侵地層中的陣列側向5條響應曲線，RA5受井眼與井斜的影響最大，RA1受井眼井斜的影響最小。由圖7可見，對于高侵地層，在井斜小于15°范圍以內，井斜角的影響同樣可以忽略。隨著侵入深度的增大，視電阻率與地層電阻率的比值逐漸增大；在同樣的侵入深度，隨著井斜角的增大，視電阻率與地層電阻率的比值出現單調增大的趨勢。在高侵地層中的陣列側向5條響應曲線，RA5受井眼與井斜的影響最大，RA1受井眼井斜的影響最小。

4 結論

1)三維有限元方法可以計算水平井/大斜度井中的陣列側向測井響應，可以進行三維地層環(huán)境中的影響因素分析。

2)陣列側向探測深度要小于CSU深側向探測深度，但是陣列側向可以在徑向進行電阻率剖面成像，從而更加直觀地反映地層電阻率變化過程。

3)井斜對于陣列側向測井響應的影響主要體現在地層的分界面附近，隨著井斜角的增大，圍巖影響也逐漸增大，井斜對探測深度最深的響應曲線影響要大于探測深度最小的響應曲線。

4)在井斜小于15°存在侵入時，水平井/大斜度井中井斜影響可以忽略。

5)在水平井/大斜度井中，必須進行井斜校正，否則會影響到計算含水飽和度的準確性。

[1] 李善軍，肖永文，汪涵明，等.裂縫的雙側向測井響應的數學模型及裂縫孔隙度的定量計算[J].地球物理學報，1996，39(6):845-852. Li Shanjun，Xiao Yongwen，Wang Hanming，et al. Mathematical Model of Dual Laterolog Response to Fracture and Quantitative Interpretation of Fracture Porosity[J]. Acta Geophysica Sinica, 1996, 39(6): 845-852.

[2] 李善軍.用三維有限元素法反演單一灰?guī)r裂縫性儲層的裂縫孔隙度和裂縫傾角[J].測井技術，1998，22(6):412-415. Li Shanjun.Inverting Fracture Porosity and Dip of Limestone Fractured Reservior Using 3-D FEM[J].Well Logging Technology, 1998, 22(6): 412-415.

[3] 劉福平，李善軍，張庚驥.利用積分方程法計算雙側向測井儀的井眼校正曲線[J].地球物理學報，1997，40(6):857-866. Liu Fuping，Li Shanjun，Zhang Gengji.Computation of Borehole Correction Curves of Dual Laterolog Tool by Integral Equation Method[J]. Acta Geophysica Sinica, 1997, 40(6): 857-866.

[4] 劉振華，胡啟，張建華.用有限元法計算侵入條件下雙側向測井儀器響應[J].石油儀器，1994，8(3):149-152. Liu Zhenhua，Hu Qi，Zhang Jianhua.Calculating Intrusion by the Finite Element Method Under the Condition of Dual Laterolog Tool Response[J]. Petroleum Instruments, 1994, 8(3): 149-152.

[5] Chemali R, Gianzero S, Su S M.The Dual Laterolog in Common Complex Situations[C]// SPWLA 29th Annual Logging Symposium. San Antonio: [s.n.], 1988: 5-8.

[6] Chen Y H, Chew W C, Zhang H J. A Novel Array Latero Log Method[J]. Log Analyst, 1998, 39(1): 23-30.

[7] Smits J W, Dubourg I, Luling M G, et al. Improved Resistivity Interpretation Utilizing a New Array Laterolog Tool and Associated Inversion Processing[C]//Annual Technical Conference and Exhibition of the Society of Petroleum Engineers. New Orleans: SPE, 1998: 831-844.

[8] 鄧少貴，李智強，陳華.煤層氣儲層裂縫陣列側向測井響應數值模擬與分析[J].煤田地質與勘探，2010，38(3):55-60. Deng Shaogui，Li Zhiqiang，Chen Hua.The Simulation and Analysis of Lateral Log Response of Fracture in Coalbed Methane Reservior[J]. Coal Geology & Exploration, 2010, 38(3): 55-60.

[9] 李智強，范宜仁，鄧少貴，等.基于改進差分進化算法的陣列側向測井反演[J].吉林大學學報:地球科學版，2010，40(5)：1199-1204. Li Zhiqiang，Fan Yiren，Deng Shaogui，et al. Inversion of Array Laterolog by Improved Difference Evolution[J]. Journal of Jilin University: Earth Science Edition, 2010, 40 (5): 1199-1204.

[10] 朱大偉，鄧少貴，范宜仁，等.大斜度井雙側向測井響應數值模擬研究[J].測井技術，2005，29(3):208-211. Zhu Dawei, Deng Shaogui, Fan Yiren, et al. Numerical Simulation of Dual Laterolog Response in Highly Deviated Wells[J]. Well Log Technology, 2005, 29(3): 208-211.

[11] 李大潛.有限元素法在電法測井中的應用[M].北京:石油工業(yè)出版社，1980:72-129. Li Daqian. Finite Element Method in the Application of the Electrical Logging[M]. Beijing: Petroleum Industry Press, 1980: 72-129.

[12] Clavier C.The Challenge of Logging Horizontal Wells[J]. Log Analyst, 1991, 32(2): 63-84.

[13] Sibbit A M, Faivre O. The Dual Laterolog Response in Fractured Rocks[C]// SPWLA 26th Annual Logging Symposium. Dallas: [s.n.], 1985: 17-20.

[14] 張庚驥.電法測井[M].東營:石油大學出版社，1996. Zhang Gengji. Electrical Logging[M]. Dongying: University of Petroleum Press, 1996.

[15] 高杰，謝然紅.大斜度井側向測井三維正演數值模擬及曲線快速校正方法研究[J].石油勘探與開發(fā)，2000，27(2):69-71. Gao Jie, Xie Ranhong. 3D Numerical Forward Modeling and Fast Correction of Dual-Laterolog for High Angle Deviated Wells[J]. Petroleum Exploration and Development, 2000, 27(2): 69-71.

Numerical Simulation of Array Laterolog Response in Horizontal and Highly Deviated Wells

Zhu Peng1，Lin Chengyan1，Li Zhiqiang2，Zhao Wenji3，Zhang Hualian4

1.SchoolofGeosciences,ChinaUniversityofPetroleum,Qingdao266580,Shandong,China2.ChinaResearchInstituteofRadiowavePropagation,Xinxiang453003,Henan,China3.WirelineLoggingCompany，DaqingDrillingEngineeringCompany,Daqing163412,Heilongjiang,China4.ChongqingInstituteofGeology&MineralResources,Chongqing400042,China

According to the operating principles of array lateral electrode arrays in horizontal wells and highly deviated wells, electric field is studied by using multiple electric field superposition method, and the 3D finite element method is employed to simulate electric field distribution of various points, and then electric field linear superposition principles are used to get the array lateral logging response. Based on the computer simulation, five array lateral logging curve radial investigation depths are obtained. The array lateral radial detection depth is smaller than the deep lateral detection depth. The effect of well deviation and variation of depth of invasion on log response in the 3D formation model is examined; and the characteristics of array lateral log response in horizontal wells and highly deviated wells is analyzed. The simulation results are as follows: well deviation has a small impact on array lateral log response when the well deviation is less than 15°, and correction of well deflection is not necessary; while the well deviation must be corrected when the horizontal well deviation is more than 60° in a highly deviated well due to a greater difference of array lateral log response between highly deviated wells and vertical wells.

horizontal/highly deviated well; array lateral; finite element; log response

10.13278/j.cnki.jjuese.201506304.

2015-02-05

國家科技重大專項(2011ZX05009-003)

祝鵬(1983--)，男，博士研究生，主要從事測井方法、解釋及測井地質學方面的科研工作，E-mail:pursuit_joe@163.com。

10.13278/j.cnki.jjuese.201506304

P631.8

A

祝鵬，林承焰，李智強，等. 水平井和大斜度井中陣列側向測井響應數值模擬.吉林大學學報:地球科學版,2015,45(6):1862-1869.

Zhu Peng，Lin Chengyan，Li Zhiqiang, et al. Numerical Simulation of Array Laterolog Response in Horizontal and Highly Deviated Wells.Journal of Jilin University:Earth Science Edition,2015,45(6):1862-1869.doi:10.13278/j.cnki.jjuese.201506304.