樂東氣田非烴類氣層的測井識別

何勝林 陳 嶸 高楚橋 張海榮

1．中海石油（中國）有限公司湛江分公司 2．油氣資源與勘探技術教育部重點實驗室·長江大學

鶯歌海盆地中央泥拱帶樂東氣田的天然氣產自埋深較淺的第四系樂東組和新近系鶯歌海組一段。儲層巖性以石英砂巖為主，少量巖屑石英砂巖和長石石英砂巖。目前氣田已經投入開發，主要生產層位為樂東組三段Ⅱ下氣組和鶯歌海組一段Ⅲ氣組，主要氣體組分為烴類氣，但是在勘探階段，部分測試層位烴類氣含量較低，最少的僅為5%左右，而非烴類氣體含量較高，最高能達到90%［1］。通過天然氣組分分析資料可知，除二氧化碳外的非烴類氣含量一般較少，在所研究的儲層中總含量平均不高于11．6%，故在本次研究中近似將二氧化碳氣看作非烴類氣；烴類氣層經化驗測定其甲烷含量占主要成分，所以將甲烷氣看作烴類氣。作為非烴類氣的二氧化碳在自然界中呈氣態存在，其臨界溫度為31．1℃，臨界壓力7．23MPa，二氧化碳具有不導電高電阻率特性［2－3］。

1 烴類氣與非烴氣定性識別方法

烴類氣與二氧化碳氣層在測井曲線上的響應特征與油、水層有較明顯的區別，而它們兩者之間的差異卻較小，這是造成用測井資料區分甲烷和二氧化碳困難的根本原因。但通過對二者理論上測井響應值的研究可知，兩者還是存在一定的差別，因而兩種氣體在測井曲線上必定存在或多或少的差異，這是用測井資料定性、定量區分烴類氣與二氧化碳氣的基礎［4］。

1）電阻率曲線響應：二氧化碳氣層為高電阻率。由于二氧化碳對鉆井液侵入帶和沖洗帶的水化作用，具較強的鉆井液減阻侵入特征。與烴類氣層相比較，電阻率稍低，鉆井液減阻侵入稍強些［5］。

2）聲波時差曲線響應：二氧化碳氣層為高時差，但與烴類氣層相比時差稍小，且少見“周波跳躍”現象［6］。

3）地層密度曲線響應：二氧化碳氣層的密度較小，但與烴類氣層相比較則高些。相同狀態下，二氧化碳的體積密度是甲烷體積密度的3倍多［7］。

4）中子孔隙度曲線響應：二氧化碳氣層的中子孔隙度很低，低于烴類氣層的中子孔隙度。二氧化碳氣含氫指數為0，甲烷氣含氫指數為0．015%［8－9］。

根據上述分析，采用交會圖技術研究了兩者在測井響應上的差別。在樂東15－1／22－1研究區選取了經測試、電纜測試取樣或試生產證實含量相對比較高的甲烷或二氧化碳層段（表1），并定義甲烷含量高于50%的氣層為以甲烷為主的氣層，否則為以非烴為主（CO2為主）的氣層，然后利用這些層段的測井資料分別作出了密度—中子、密度—聲波交會圖（圖1、2）。從圖1、2中可以看出：烴類氣與二氧化碳氣體在測井響應值上是存在一定差別的，從而可以通過這種差異來區分烴類氣與二氧化碳氣層［10］。

通過分析圖1、2可以發現，中子測井與密度測井能較好地區分以甲烷為主的烴類氣層與二氧化碳氣層；聲波測井對兩類氣體的區別能力相對較差。隨鉆測井資料由于是地層剛鉆開時獲得的資料，對兩者的區分效果要比電纜測井資料好［11］。

在物性比較好，含氣飽和度高的情況下，兩類氣體的測井響應特征差別明顯；隨著泥質含量增多，物性變差，含氣飽和度降低等情況會使得以上差異不明顯。這些因素的影響將在定量識別中消除。

2 烴類氣與非烴氣定量區分方法

筆者區分烴類氣與非烴氣的基本思路：在單位體積條件下假定各組分在地層中的相對百分含量，建立雙水多礦物地層組分分析物理模型和帶約束的測井超定線性方程組，進而運用線性最小二乘法原理，將其轉換為求解極值問題的數學目標函數，最終采用最優化算法得到地面條件下烴類氣的相對含量，從而達到識別氣層類型的目的。由于除二氧化碳外的非烴氣含量一般較少，且測井響應值與二氧化碳接近，為方便起見將其歸到二氧化碳中，作為同一組分看待［12］。

2．1 物理模型

對含油氣的儲集層來說，儲集層可以看成是由具有不同性質的組分組成的，這些組分包括：不動油、可動油（或非烴氣）、天然氣、可動水、束縛水、泥質以及巖石的各種骨架礦物（表2）。

假設組分不動油、可動油（或非烴氣）、可動水、束縛水、天然氣、泥質以及巖石的各種骨架礦物在地層中的相對含量分別為：xor，xCO2，xfw，xbw，xgas，xsh，xma1，xma2，…，xmak，則有如下關系式。

孔隙度：

地層含水飽和度：

束縛水飽和度：

泥質含量：

表1 甲烷與二氧化碳測井響應值的差別表

圖1 烴類氣與非烴氣儲層中子—密度交會圖

圖2 烴類氣與非烴氣儲層中子—聲波時差交會圖

表2 雙水多礦物地層組分模型表

2．2 數學模型

根據以上物理模型，可寫出各種測井儀器的響應方程式。例如，密度測井的響應方程為：

式中ρor，ρCO2，ρfw，ρbw，ρgas，ρsh，ρma1，ρma2，…，ρmak分別表示地層中不動油氣、二氧化碳、自由水、束縛水、烴類氣、泥質、巖石骨架礦物（1～k種）的體積密度，g／cm3。

為簡便起見，將式（5）寫成：

式中n表示組成地層的組分個數；xj表示第j種組分的相對含量。

同理可寫出其他測井儀器的響應方程，用通式表示為：

式中m表示測井儀器的個數；B表示地層對測井儀器的響應值。

解以上由m個方程組成的方程組，就可以求得xj，這就是待解決的組分含量［3］。

2．3 烴類氣含量的計算

由氣體狀態方程可得到井底條件下體積為Vgf的天然氣（烴類氣）在地面條件下的體積Vgs為［7］：

式中Ts為地面溫度，K；pgf為地層壓力，MPa；Zf為烴類氣在井底條件下的壓縮因子，無因次；Tf為井底溫度，K；pg為地面壓力，MPa。

同樣可寫出井底條件下體積為CCO2f的非烴氣在地面條件下的體積VCO2s，即

式中ZCO2f為非烴類氣在井底條件下的壓縮因子，無因次。

地面條件烴類氣相對含量為：

若巖石體積為VT，則Vgf＝xgasVT，VCO2f＝將式（9）代入式（8）并整理得到：

式中xgas為烴類氣在地層中的相對含量，小數；xCO2為非烴氣在地層中的相對含量，小數。

地層流體密度為：

式中ρgas、ρCO2、ρfw分別為地層條件烴類氣、非烴氣和地層水的密度，g／cm3。

3 應用效果

圖3、4分別為 LD15－1－X井、LD22－1－Y井雙水多礦物地層組分分析程序處理的解釋成果圖。

圖3顯示LD15－1－X井1 572～1 577m 段第2道中子密度出現交會，第3道電阻率增大，第5道孔隙度差、比值識別法出現明顯的氣層“包絡”特征，第8道測井計算的含水飽和度與束縛水飽和度幾乎重合，含水飽和度大致為30%，第6道為該方法定量計算出的烴類氣相對含量，大概在40%左右，DST測試證實該層甲烷含量為34．4%，測試的結果與定量計算的結果相近。

圖4顯示LD22－1－Y井963～985m段測井定量計算烴類氣含量為75%左右，測試證實該層甲烷含量為81%。通過鶯歌海盆地15口井19個層位測井計算烴類氣含量與實測結果對比，在統計的19層中，絕對誤差在20%以內的層有15層，占79%；絕對誤差在10%以內的層有12層，占63%［13］。表明該定量方法能較為有效地計算烴類氣含量，從而達到區分烴類氣與非烴氣的目的。

圖3 LD15－1－X井處理成果圖（DST1）

圖4 LD22－1－Y井處理成果圖（DST2）

4 結論

1）樂東氣田烴類氣與二氧化碳在測井響應值上存在一定差別，通過這種差異可以定性判斷烴類與二氧化碳氣層。中子—密度交會圖能較好地區分以甲烷為主的烴類氣層與二氧化碳氣層。

2）利用地層組分分析模型和最優化理論定量計算烴類氣與非烴氣含量。計算結果表明：63%的層烴類氣含量計算誤差能控制在10%以內；79%的層烴類氣含量計算誤差能控制在20%以內。該方法可以較好地用測井資料區分烴類氣與非烴氣。

3）多口井的實際資料處理證明，常規測井非烴類氣層定量解釋確立的方法和模型適合于南海北部盆地天然氣田的測井解釋評價。

［1］劉中奇，杜春玲．非烴類氣層測井識別方法［J］．測井技術，1998，22（2）：107－110．LIU Zhongqi，DU Chunling．Identification method of nonhydrocarbon gas zones on logging［J］．Well Logging Technology，1998，22（2）：107－110．

［2］高楚橋．復雜儲層測井評價方法［M］．北京：石油工業出版社，2003．GAO Chuqiao．Well logging evaluation of complex reservoirs［M］．Beijing：Petroleum Industry Press，2003．

［3］周倫先，褚小兵．CO2氣井相態特征［J］．油氣井測試，2006，15（5）：35－37．ZHOU Lunxian，CHU Xiaobing．Phase behavior characteristic of gas well with rich－content CO2［J］．Well Testing．2006，15（5）：35－37．

［4］何勝林，張海榮，陳嶸，等．東方1－1氣田低電阻率氣層測井評價方法［J］．天然氣工業，2012，32（8）：27－30．HE Shenglin，ZHANG Hairong，CHEN Rong，et al．Logging evaluation methods for low－resistivity gas layers in the Dongfang 1－1Gas Field，Yinggehai Basin［J］．Natural Gas Industry，2012，32（8）：27－30．

［5］王德喜，曾文沖，汪志鋒．松南氣田火山巖儲層流體性質的識別［J］．西南石油大學學報：自然科學版，2009，31（6）：15－20．WANG Dexi，ZENG Wenchong，WANG Zhifeng．Reservoir fluid property identification of volcanic rock in Songnan Gas Field［J］．Journal of Southwest Petroleum University：Science ＆ Technology Edition，2009，31（6）：15－20．

［6］王磊，王學琴，吳勝，等．利用彈性模量計算含氣飽和度方法研究［J］．西南石油大學學報：自然科學版，2011，33（4）：69－72．WANG Lei，WANG Xueqin，WU Sheng，et al．Method of gas saturation prediction by elastic modules［J］．Journal of Southwest Petroleum University：Science ＆ Technology E－dition，2011，33（4）：69－72．

［7］司馬立強，趙冉，王培春，等．普光縫洞性儲層流體性質測井判別適應性［J］．西南石油大學學報：自然科學版，2010，32（1）：11－15．SIMA Liqiang，ZHAO Ran，WANG Peichun，et al．The suitability analysis of log distinguishing methods of fluid nature of fracture and cave reservoir in Puguang Field［J］．Journal of Southwest Petroleum University：Science ＆Technology Edition，2010，32（1）：11－15．

［8］李山生，趙輝．利用彈性參數識別氣、水層［J］．西南石油大學學報：自然科學版，2012，34（3）：83－88．LI Shansheng，ZHAO Hui．Identification of water and gas reservoirs with elastic parameters［J］．Journal of Southwest Petroleum University：Science＆Technology Edition，2012，34（3）：83－88．

［9］李同華．基于偶極橫波資料的火山巖裂縫及油氣識別［J］．西南石油大學學報：自然科學版，2009，31（6）：45－50．LI Tonghua．Igneous rock fractures and hydrocarbon identification based on DSL data［J］．Journal of Southwest Petroleum University：Science ＆ Technology Edition，2009，31（6）：45－50．

［10］于世元．復雜氣層識別技術研究［J］．天然氣工業，1999，19（1）：67－70．YU Shiyuan．Research on the techique of identifying the complex gas bearing formation［J］．Natural Gas Industry，1999，19（1）：67－70．

［11］VARGAFTIK N B．Tables on the thermophysical properties of liquids and gases［M］．Washington DC：Hemisphere Publishing Corporation，1975：156－166．

［12］高楚橋，袁云福，吳洪深，等．鶯歌海盆地束縛水飽和度測井評價方法研究［J］．天然氣工業，2003，23（5）：38－40．GAO Chuqiao，YUAN Yunfu，WU Hongshen，et al．Research on the log evaluation method of irreducible water saturation in Yinggehai Basin［J］．Natural Gas Industry，2003，23（5）：38－40．

［13］吳洪深，高華，林德明．南海西部海域非烴類氣層測井識別及解釋評價方法［J］．中國海上油氣，2012，24（1）：21－24．WU Hongshen，GAO Hua，LIN Deming．A methodology on logging identification and interpretation of non－hydrocarbon gas zones in the western South China Sea［J］．China Offshore Oil and Gas，2012，24（1）：21－24．