隨鉆脈沖中子測井識別天然氣的數(shù)值模擬

摘 要

摘要：由于天然氣具有密度低、黏度小等特點，利用3種孔隙度資料可以定性識別天然氣層，但定量評價存在困難。為此，在隨鉆過程中依據(jù)天然氣與油、水的含氫指數(shù)不同，利用脈沖中子測井
摘要：由于天然氣具有密度低、黏度小等特點，利用3種孔隙度資料可以定性識別天然氣層，但定量評價存在困難。為此，在隨鉆過程中依據(jù)天然氣與油、水的含氫指數(shù)不同，利用脈沖中子測井技術記錄的近、遠探測器熱中子或俘獲伽馬計數(shù)率比的相對變化量來定量確定含氣飽和度；在此基礎上，利用蒙特卡羅方法建立計算模型，模擬不同井眼和地層條件下脈沖中子測井遠近探測器記錄的熱中子或伽馬計數(shù)，研究其比值與含氣飽和度的測井響應。結(jié)果表明：油水層和氣層的計數(shù)相對變化量能反映地層的含氣飽和度，孔隙度越大，相對比值越大，對氣層的定量評價越準確；巖性、泥質(zhì)含量、地層水的礦化度、井眼流體和尺寸以及鉆井液侵入等因素都會對天然氣地層的脈沖中子測井響應產(chǎn)生影響?？傊?，利用脈沖中子測井技術可以定量評價天然氣層，對提高天然氣識別能力和氣田高效勘探開發(fā)具有重要意義。
關鍵詞：天然氣；地層；脈沖中子測井；含氣飽和度；影響；因素；蒙特卡羅法
    中國天然氣工業(yè)經(jīng)歷了近10年的快速發(fā)展，天然氣儲量進入新的增長高峰期。但是天然氣勘探具有埋深大、物性差、圈閉條件復雜的特點，其發(fā)現(xiàn)難度逐漸增大。目前在常規(guī)天然氣層的識別和評價方法很多，主要包括聲波時差與中子伽馬曲線重疊、雙孔隙度重疊、三孔隙度重疊、差值或者比值法、縱波時差差比法、核磁共振測井、成像測井、補償中子時間推移測井、巖性密度測井、熱中子衰減時間測井、碳氧比測井等[1～4]。近年來脈沖中子測井技術在氣層的識別和評價過程中發(fā)揮了重要作用，Trcka D等[5]利用儲層動態(tài)監(jiān)測儀(RPM)不同探測器的非彈伽馬計數(shù)比值建立了與含氣飽和度的響應關系，Hamada G M[6]等利用中子壽命和過套管電阻率測井進行砂巖氣層含烴物質(zhì)監(jiān)測。
    隨鉆條件下利用電阻率來判斷油氣層，但由于電阻率相差不大，很難定量確定含氣飽和度。斯倫貝謝公司在新一代隨鉆測井儀器增加了地層宏觀俘獲截面測量，Gauthier P J等[7]結(jié)合隨鉆電阻率和地層宏觀俘獲截面來確定含油飽和度，取得了較好的效果。由于天然氣含氫指數(shù)低且具有挖掘效應，其對中子的作用與油和水存在較大差別，通過脈沖中子測井技術得到的遠近探測器的熱中子或俘獲伽馬計數(shù)率比，建立其相對變化量與孔隙度的響應關系，可以定量確定含氣飽和度。利用蒙特卡羅方法模擬不同孔隙度、巖性等地層和井眼條件下的中子分布，研究熱中子或俘獲伽馬計數(shù)率比值和含氣飽和度的響應及井眼和地層因素對確定含氣飽和度的影響，為隨鉆條件下天然氣儲層含氣飽和度的定量評價提供技術支持。
1 確定含氣飽和度的方法
    新型集成隨鉆測井儀器包括方位自然伽馬測量、伽馬-伽馬密度測量和脈沖中子測量等3部分，而脈沖中子測量不僅能確定孔隙度，還可確定地層宏觀俘獲截面、元素含量以及中子伽馬密度等多種地層參數(shù)。
    天然氣和原油的元素成分比較接近，但其密度比原油和水小得多，含氫指數(shù)遠小于油和水，天然氣地層的熱中子通量分布會因氫的含量降低而有所變化，利用不同源距探測器的熱中子計數(shù)可以反映地層的含氣特性。
    根據(jù)脈沖中子測井技術得到的近、遠探測器的熱中子或俘獲伽馬計數(shù)比值為：
 
式中N1、N2分別為近探測器和遠探測器的熱中子或俘獲伽馬計數(shù)率；r1、r2分別為短源距和長源距；Lf為地層的快中子減速長度。
顯然飽含氣層的計數(shù)比值分別為：
 
式中Lfg為氣層的快中子減速長度。
參照Trcka D[5]提出反映含氣飽和度的非彈性散射伽馬計數(shù)比值，定義熱中子或俘獲伽馬計數(shù)比值的相對變化量參數(shù)(D)為：
 
    不同孔隙度和含氣飽和度地層的快中子減速長度(Lf)不同，其D也不同；利用孔隙度資料和D的響應關系可定量確定含氣飽和度。
2 蒙特卡羅模擬方法
    MCNP[8]是一種能夠模擬連續(xù)能量的中子、光子、電子、中子/光子和光子/電子粒子在任意幾何形狀介質(zhì)中的輸運過程，在核測井方法基礎和應用研究中發(fā)揮了巨大作用。
    利用蒙特卡羅方法建立隨鉆條件下的計算模型，井眼直徑為20cm；地層徑向半徑分別為10～70cm、長為120cm的柱體，把整個地層劃分成環(huán)距為7.5cm、厚度5cm的相鄰柵元，地層分別填充不同巖性、孔隙度和流體的物質(zhì)；測井儀器外徑為45mm，脈沖中子源處于儀器左端15cm處，脈沖寬度為3μs，長、短源距分別為72.5cm和42.5cm。為了簡化模擬過程，沒有考慮熱中子探測器He-3管的響應，中子源和探測器間放置理想屏蔽體(圖1)。模擬時選取的截面數(shù)據(jù)庫為ENDF/B-VI.0，采用F4柵元計數(shù)方式，記錄熱中子的時間譜道寬為30μs，記錄時間間隔為30～1800μs，得到長短源距處的熱中子時間譜，研究利用遠近探測器熱中子計數(shù)比值相對變化量的方法來定量評價氣層。

3 模擬結(jié)果及討論
3.1 含氣飽和度和孔隙度不同時的測井響應
    利用上述計算模型，井眼內(nèi)流體為淡水，地層為分別飽含油和氣的砂巖，其中油的密度為0.87g/cm3，天然氣(CH4)的密度為0.2g/cm3，在地層孔隙度為30%，改變含氣飽和度依次為0、20%、40%、60%、80%和100%，含氣飽和度分別為0、20%、50%和80%，改變孔隙度依次為10%、15%、20%、25%、30%、35%和40%時，模擬記錄遠近探測器處的熱中子計數(shù)，得到不同含氣飽和度時D與孔隙度的關系曲線(圖2)。
 

    由圖2可以得到：遠近探測器熱中子計數(shù)率比的相對變化量與孔隙度和含氣飽和度都有關；在孔隙度一定時，D隨地層含氣飽和度的增加呈線性減小，D能夠反映含氣飽和度的大?。缓瑲怙柡投纫欢〞r，D隨孔隙度的增加而增加，且含氣飽和度越小這種變化越劇烈；孔隙度越大時不同含氣飽和度引起的D的差異越大。故利用脈沖中子測井技術得到的遠近探測器熱中子計數(shù)率比的相對變化量和孔隙度資料可以確定含氣飽和度，且孔隙度越大的地層越有利于評價氣層。
3.2 巖性和泥質(zhì)含量的影響
    其他條件不變時，地層孔隙中分別飽含油和氣，骨架分別為砂巖和石灰?guī)r，改變地層的孔隙度，模擬得到的熱中子計數(shù)比相對變化量與孔隙度的變化關系見圖3a。在考慮泥質(zhì)時，泥質(zhì)由細砂、黏土和水組成，體積百分比分別為8%、72%和20%，黏土礦物包括20%的蒙脫石、21%的伊利石、24%的綠泥石、20%的高嶺石和5%黑云母，密度為2.21568g/cm3，泥質(zhì)含量分別為0、20%和40%的砂巖地層時D與孔隙度的關系如圖3-b所示。
 

    由圖3可以看出：砂巖地層和石灰?guī)r地層的相對比值(D)隨孔隙度的變化趨勢相同，在孔隙度較小時變化較大，當孔隙度增加到一定程度時開始變緩；孔隙度較小時不同巖性引起D的差異較小，隨著孔隙度的增加地層巖性不同D的差異越大；孔隙度相同時下砂巖地層的D要比石灰?guī)r地層的大，原因是石灰?guī)r對中子的減速能力比砂巖大；泥質(zhì)含量引起的D的差異隨孔隙度的增加而增加，且孔隙度相同時泥質(zhì)含量越大D越小。因此利用熱中子計數(shù)率比相對變化量確定含氣飽和度時受巖性和泥質(zhì)含量的影響，砂巖地層對氣反應靈敏，泥質(zhì)含量越大對確定地層含氣飽和度越不利。
3.3 地層水礦化度的影響
    利用同樣的計算模型，分別飽含不同礦化度地層水和氣的砂巖地層，其中地層水礦化度分別為0、30000、50000和100000μg/g，改變地層孔隙度，模擬得到熱中子計數(shù)相對比值與孔隙度的關系(圖4)。
 

    由圖4可以看出：不同礦化度地層水時熱中子計數(shù)率相對變化量(D)仍隨孔隙度的增加而增加，地層水礦化度越高變化越快；同一孔隙度條件下，地層水礦化度越大時，D越大，主要由于地層水礦化度越大，對遠近探測器處熱中子俘獲的能力差異越大。因此地層水礦化度越高越有利于含氣評價，但需對地層水礦化度的影響進行校正。
3.4 井眼尺寸和流體的影響
    其他計算條件不變，井眼內(nèi)流體分別為氣(CH4，密度為0.19g/cm3)、淡水和30000μg/g的礦化水，地層為飽含油和氣的砂巖，依次改變地層的孔隙度，模擬得到熱中子計數(shù)比的相對變化量與孔隙度的關系(圖5-a)；同樣其他條件不變，井眼內(nèi)充滿氣，井眼直徑為分別為20、22.827、27.907cm和30.447cm時的關系如圖5-b所示。

    由圖5可以看出：低孔隙度地層井眼流體不同引起的相對比值差異較小，高孔隙度地層差異較大，且井眼內(nèi)充滿氣時相對變化量(D)越小，對確定含氣飽和度越不利；井眼內(nèi)流體為水時，礦化度的高低對D的影響很小。在孔隙度一定的地層，井眼尺寸越大時D越小，這是主要是因為井眼流體為氣，井眼直徑越大挖掘效應影響越大，故井眼尺寸越大越不利于含氣飽和度的評價。
3.5 鉆井液侵入的影響
    計算模型不變，建立孔隙度為15%和30%飽含天然氣的砂巖地層，井眼內(nèi)充滿淡水鉆井液，其礦化度為5000μg/g，密度為1.33402g/cm3，向井軸周圍均勻侵入深度分別為0、3、6、9、12、15、18、21、24、27、30、33、36、39、42、45、48、51、54、57、60cm，在不考慮泥餅的理想條件下，利用前面相同的方法得到響應關系如圖6所示。
    由于天然氣的含氫指數(shù)小于鉆井液，隨著侵入深度的增加，遠近探測器熱中子計數(shù)率比的相對變化量也隨之增加，增加至30cm左右時幾乎不再發(fā)生變化，且孔隙度不同的地層受鉆井液侵入深度的影響差別不大。因此在識別氣層時要考慮到鉆井液侵入的影響，但在隨鉆測井過程中鉆井液侵入深度很淺，脈沖中子測井受其影響很小，實際資料處理時可忽略鉆井液侵入的影響。
4 結(jié)論
    1) 由于天然氣和油水的含氫指數(shù)不同，利用脈沖中子測井技術記錄的近、遠探測器熱中子或俘獲伽馬計數(shù)率比的相對變化量隨含氣飽和度的增加呈線性減小，可以定量確定含氣飽和度。
    2) 利用熱中子或俘獲伽馬計數(shù)比相對變化量確定含氣飽和度時受井眼和地層條件影響，通過利用蒙特卡羅模擬研究得到：相對變化量隨孔隙度增加而增加，孔隙度越大越有利于評價含氣飽和度；孔隙度相同時，砂巖地層的相對變化量比石灰?guī)r地層的大，而泥質(zhì)含量越小，地層水礦化度越高，相對變化量越大，對確定含氣飽和度越有利。
    3) 井眼流體為水、井眼尺寸越小時，相對變化量越大，對確定含氣飽和度越有利；在高孔隙度地層，井眼流體為氣時的相對變化量越小，而井眼內(nèi)水的礦化度對其影響很??；隨著鉆井液侵入深度的增加，相對變化量也隨之增加，增加至一定深度不再發(fā)生變化。顯然在利用脈沖中子測井技術確定含氣飽和度時需對各種影響因素進行校正，在實際應用過程中還需進一步研究。
參考文獻
[1] 黃布宙，潘保芝，李舟波.大慶徐家圍子地區(qū)深部致密砂礫巖氣層識別[J].吉林大學學報：地球科學版，2003，33(4)：490-494.
[2] 于世元.復雜氣層識別技術研究[J].天然氣工業(yè)，1999，19(1)：67-70.
[3] 張筠.川西地區(qū)致密碎屑巖儲層測井評價技術[J].測井技術，2004，28(5)：419-422.
[4] 沈慧萍，史振勤，喬國安，等.深層氣層識別方法研究[J].天然氣工業(yè)，2005，25(10)：47-48.
[5] TRCKA D，GILCHRIST A，RILEY S，et al.Field trials of a new method for the measurement of formation gas using Pulsed-Neutron Instrumentation[C]∥SPE Annual Technical Conference and Exhibition，24-27 September 2006，San Antoniu，Texas，USA：SPE，2006，SPE 102350.
[6] HAMADA G M，HEGAXY A A.Hydrocarbon potential monitoring in gas sandstone reservoirs using TDT and CHFR techniques[C]∥EUROPEC/EAGE Conference and Exhibition，11-14 June 2007，London，U K：SPE，2007，SPE 105003.
[7] GAUTHIER P J，HUSSAIN H，BOWLING J，et al.Determination of water-producing zones while underbalanced drilling horizontal wells-integration of Sigma log and realtime production data[C]∥SPE Middle East Oil and Gas Show and Conference，11-14 March 2007，Kingdom of Bahrain：SPE，2007，SPE 105166.
[8] BRIESMEISTER J F.MCNPTM-A general Monte Carlo N-Particle transport code，version 4C[CP].Los Alamos National Laboratory Report LA-13709-M，2000.
 
(本文作者：張鋒1，2 袁超1，2 侯爽1，2 王新光1，2 1.中國石油大學(華東)地球資源與信息學院；2.中石油測井重點實驗室中國石油大學(華東)研究室)