向長生 王智鋒 亢武臣

（1.中國石油大學(xué)（華東）石油工程學(xué)院，山東青島 266580；2.勝利石油管理局鉆井工藝研究院，山東東營 257017）

國內(nèi)外油氣勘探開發(fā)已經(jīng)進(jìn)入后期：一方面，老油區(qū)存在著大量的邊緣油氣藏、獨立小油氣藏、復(fù)雜斷塊油氣藏、超薄油氣藏等難動用儲量；另一方面，新探明儲區(qū)大多處于地層深處、海洋、灘海、沙漠等地區(qū)，特殊的油藏開發(fā)條件和復(fù)雜的油藏地質(zhì)結(jié)構(gòu)使得勘探開發(fā)成本和難度大大增加，對鉆井技術(shù)的要求也越來越苛刻［1-2］。傳統(tǒng)的鉆采技術(shù)已經(jīng)無法從根本上提高采收率和實現(xiàn)對上述難動用儲量的油氣資源開采。為實現(xiàn)穩(wěn)產(chǎn)，分支井、水平井、大位移井等新型鉆井工藝的應(yīng)用越來越多［3-4］，對鉆井施工過程中控制鉆頭穿越油藏儲層到達(dá)靶點的要求越來越高。因此，在隨鉆測控過程中，準(zhǔn)確判斷地層邊界位置，尤其是準(zhǔn)確預(yù)測儲層邊界，對于提高井眼軌跡控制水平和儲層有效鉆遇率具有十分重要的意義［5］。

筆者利用自然伽馬強(qiáng)度分布的基本原理，建立了無限均勻水平地層模型中探測點距離地層邊界的解析求解法，對于隨鉆測井?dāng)?shù)據(jù)解釋以及分辨地層具有一定參考價值。

1 隨鉆自然伽馬測量技術(shù)

1.1 基本原理

由于任何物體或多或少都含有天然放射性元素，尤其以自然伽馬射線穿透能力最強(qiáng)。通過求解地層中自然伽馬射線的強(qiáng)度情況，就可以得出所處地層的巖性，如地層的厚度、以及相鄰地層的邊界等特性，所以在常規(guī)LWD 中一般會提供自然伽馬測量數(shù)據(jù)。

1.2 自然伽馬強(qiáng)度分布

為簡化分析，考慮探測點所在地層為無限均勻，假設(shè)探測點所處地層的密度為ρ，每克中含q 克放射性物質(zhì)，每克放射性物質(zhì)平均每秒發(fā)射a 個自然伽馬光子，地層對自然伽馬射線的吸收系數(shù)為μ。

地層體積元dV 在探測點處形成的自然伽馬通量dJv為

采用球坐標(biāo)系，dV = r2·sinθ ·d θ d φd r ，則上式寫成

對式（2）進(jìn)行積分后得

式（3）表示半徑為r 的均勻球體空間在球心形成的自然伽馬通量。

根據(jù)式（3），設(shè)探測點處99%的自然伽馬通量由半徑為R 的球體造成，則有1-e-μr=0.99，由此可得μr=4.065，從而根據(jù)不同地層的μ 值算出相應(yīng)R 值。一般認(rèn)為自然伽馬射線在沉積巖中的平均穿透深度約30 cm，如果考慮套管、鉆桿、鉆井液等因素影響，R的實際探測深度一般不會超過20 cm［6-7］。

2 無限均勻水平地層中自然伽馬強(qiáng)度分布模型

根據(jù)式（3），以兩無限均勻的水平地層分界面為橫軸建立坐標(biāo)系，如圖1 所示：Ⅰ區(qū)與Ⅱ區(qū)地層參數(shù)分別為μ1、a1、q1、ρ1和μ2、a2、q2、ρ2，R1表示自然伽馬在Ⅰ區(qū)地層通過的距離，R2表示自然伽馬在Ⅱ區(qū)地層通過的距離，R=R1+R2。

圖1 無限均勻水平地層中自然伽馬通量計算模型

（1）地層Ⅰ的自然伽馬分布

（2）地層Ⅱ的自然伽馬分布

根據(jù)式（5）和式（7）可得：

（2）當(dāng)在地層邊界附近時，有：

當(dāng)m ≥0 時，

當(dāng)m=0 時，

當(dāng)m ≤0 時，

3 仿真與分析

圖2 是根據(jù)上述模型利用matlab 進(jìn)行模擬地層伽馬探測的結(jié)果，計算探測點從Ⅰ區(qū)地層進(jìn)入Ⅱ區(qū)地層，然后再進(jìn)入其中地層Ⅰ，其中假設(shè)各項地層參數(shù)均為已知，并保證Ⅱ區(qū)地層厚度足夠大，即保證伽馬探測值僅受相鄰兩地層影響；如果需要考慮較薄地層，只須將上面推導(dǎo)中相關(guān)m 的值變動即可，在程序?qū)崿F(xiàn)上也具有很好的可移植性。

圖2 無限均勻水平地層中自然伽馬通量matlab 模擬

從圖中可以看出，在地層邊界附近，伽馬值都有比較明顯的變化，并且在不同地層其變化的速率也不同，能夠很好地區(qū)分出地層邊界。

4 結(jié)論

（1）根據(jù)自然伽馬在地層中的強(qiáng)度分布原理，推導(dǎo)出無限均勻水平地層中探測點距離邊界的數(shù)學(xué)模型并給出了解析解。

（2）根據(jù)該模型可以很好地描述無限均勻水平地層中的自然伽馬強(qiáng)度分布情況，能夠識別探測點與地層邊界的距離，對地質(zhì)導(dǎo)向等技術(shù)具有一定參考價值。

（3）模型給出的解析解形式十分簡單，便于編程計算分析提高效率，并且可以單獨開發(fā)成計算模塊供其它程序調(diào)用。

（4）該模型還有較大的完善空間，比如方位性、薄地層的識別等，還有待進(jìn)一步分析研究以便盡快在實踐中得以應(yīng)用。

［1］ LEE W J, HOLDITCH SA F. Fracture evaluation with pressure transient testing in low-permeability gas reservoir［J］. Journal of Petroleum Technology, 1981, 9（6）∶ 1762-1776.

［2］ TURGAY Eetekin. Dynamic gas slippage-a unique pualmechanism approach to the flow gas slippage in tight formation［R］. SPE 12045.

［3］ 張桂林.勝利油田水平井鉆井技術(shù)現(xiàn)狀與發(fā)展趨勢［J］.石油鉆探技術(shù)，2005， 33（2）：66-70.

［4］ 孫清德.中國石化集團(tuán)鉆井技術(shù)現(xiàn)狀及展望［J］.石油鉆探技術(shù)，2006， 34（2）：1-6.

［5］ 狄富春，呂建云，趙立文，等.提高薄油層水平油層鉆遇率的技術(shù)研究［J］.石油鉆探技術(shù)，2008， 36（2）：25-27.

［6］ 龐巨豐.核測井物理基礎(chǔ)［M］.北京：石油工業(yè)出版社， 2005：79-80.

［7］ 吳官生.地質(zhì)導(dǎo)向法用于鉆頭前方地質(zhì)情況預(yù)測［J］.石油物探譯叢，1998，12（6）：76-82.