鄧模

（中國石化石油勘探開發(fā)研究院無錫石油地質(zhì)研究所，江蘇無錫214126）

?

一種計算石油包裹體氣液比的新方法

鄧模

（中國石化石油勘探開發(fā)研究院無錫石油地質(zhì)研究所，江蘇無錫214126）

摘要：石油包裹體的氣液比是熱力學(xué)模擬的敏感性參數(shù)，應(yīng)用石油包裹體熱力學(xué)模擬方法可以獲得古油氣充注時的壓力，通過石油包裹體二維平面分析和三維形態(tài)重構(gòu)，可以獲取石油包裹體的氣液比。實例研究證明，在一定的石油包裹體形態(tài)范圍內(nèi)，通過二維平面分析和三維形態(tài)重構(gòu)獲取的石油包裹體氣液比與共聚焦激光-熒光掃描顯微鏡所測得的氣液比的相對誤差可以控制在±10%內(nèi)。在缺少共聚焦激光-熒光掃描顯微鏡條件下，石油包裹體二維平面分析和三維形態(tài)重構(gòu)方法比目估測量法能更準(zhǔn)確地計算石油包裹體氣液比，可以滿足石油包裹體捕獲壓力模擬的精度要求。

關(guān)鍵詞：石油包裹體；氣液比；熱力學(xué)模擬；捕獲壓力

流體包裹體是指在礦物結(jié)晶生長時，被捕獲在礦物晶格的缺陷或空穴內(nèi)的成礦液體，只要流體包裹體在形成以后未發(fā)生破裂和漏失，它就包含著成礦、成藏流體的原始信息[1]。油氣運移過程中形成的流體包裹體，一般生成于碳酸鹽巖和碎屑巖中的方解石脈、石英脈、石英次生加大邊、石英顆粒裂縫愈合處以及同期形成的螢石、硬石膏等自生礦物中[2-4]，是油氣運移聚集過程的直接標(biāo)志。流體包裹體，特別是其中的石油包裹體，在油氣成因、油氣運移以及成藏等研究方面具有不可替代的作用[5-6]。

包裹體氣液比可以表示氣泡的充填度，它不但是包裹體分類中的一個重要參數(shù)，同時也是包裹體捕獲溫度和捕獲壓力計算的基礎(chǔ)數(shù)據(jù)[1]，由特定流體體系的溫度、體積和壓力決定[7]。已知石油包裹體的氣液比、化學(xué)組成和均一溫度及其同期鹽水包裹體的均一溫度，就可以進行熱力學(xué)模擬，進而獲得石油包裹體最小捕獲壓力[8]。因此，石油包裹體的氣液比是熱力學(xué)模擬的一個重要敏感性參數(shù)，氣液比計算的準(zhǔn)確性直接影響到捕獲壓力預(yù)測的精度。較為常用的石油包裹體氣液比研究手段主要有共聚焦激光-熒光掃描顯微鏡測量法和目估測量法2種。

共聚焦激光-熒光掃描顯微鏡（CLFSM）以激光作為光源，通過激光孔及物鏡的聚焦作用，形成一個光點照射到樣品上，激發(fā)出熒光。共聚焦激光-熒光掃描顯微鏡測量原理：首先將石油包裹體薄片置于激光掃描共聚焦顯微鏡下，確定石油包裹體的頂、底界后，進行深度掃描，得到不同深度的系列切片，然后利用三維重建軟件，對這些切片進行三維圖像重構(gòu)，分別計算出整個石油包裹體和氣泡的體積，從而獲得石油包裹體氣液比，其精度可以達到95%[9-10].雖然共聚焦激光-熒光掃描顯微鏡測量方法的精度較高，但仍然存在一定局限性，表現(xiàn)在2方面：一是石油包裹體受激光照射時液體部分會發(fā)出較強熒光，使得氣泡部分被照亮，造成氣液比測試結(jié)果偏小；二是由于石油包裹體形態(tài)多不規(guī)則，利用圓臺體積公式來計算氣液體積比時，通常會產(chǎn)生一定的系統(tǒng)誤差[9]，使計算結(jié)果失真。

目估測量法是目前較為常用的氣液比計算方法，其原理是在熒光顯微鏡下，通過目測二維投影得到石油包裹體中不同相態(tài)的面積，對比已知各種石油包裹體三維形態(tài)圖的標(biāo)準(zhǔn)氣液比來綜合確定[11-12]。該方法用目估的面積比代替了體積比，雖然成本低、流程簡單，但由于目測過程中受人為因素影響較大，測試精度難以滿足較高的要求。

為了在保證測試精度的前提下降低實驗成本、提高便捷性，本文通過一系列實驗對比研究，結(jié)合計算機編程，提出了一種新的石油包裹體氣液比計算方法，兼顧了共聚焦激光-熒光掃描顯微鏡測量法和目估測量法的優(yōu)點。

1　石油包裹體氣液比計算方法基本原理

捕獲于成巖礦物晶體缺陷和裂紋中的石油包裹體除部分呈負(fù)晶形之外，多數(shù)呈極其不規(guī)則的形態(tài)，石油包裹體大小從一微米以下到數(shù)十微米不等，顯微鏡下的包裹體以二維形態(tài)呈現(xiàn)。在透射光下可利用長寬比和不規(guī)則度2個參數(shù)來表征石油包裹體的二維形態(tài)，并基于2個假設(shè)，一是認(rèn)為石油包裹體中的氣泡是一個球體，二是認(rèn)為石油包裹體由n個高為Lmax/n、底面積為πR2/4的圓柱體組成。在此基礎(chǔ)上，進行石油包裹體各相態(tài)體積三維重構(gòu)，運用CorelDraw軟件進行VisualBasic編程，獲取石油包裹體體積，結(jié)合氣泡體積計算，最終得到石油包裹體的氣液比。

2　石油包裹體體積及氣相體積計算

2.1石油包裹體體積計算

（1）二維形態(tài)分析石油包裹體形態(tài)種類豐富，其中不規(guī)則形態(tài)占絕大多數(shù)，負(fù)晶形石油包裹體相對較少，筆者采用長寬比（Lmax/Lmin）和不規(guī)則度C2/（4πA）參數(shù)來表征顯微鏡下觀測到的石油包裹體的二維形態(tài)。其中Lmax與Lmin分別代表將石油包裹體二維形態(tài)沿某一端點旋轉(zhuǎn)90°過程中測定的最大水平長度與最小水平長度，C和A分別為石油包裹體二維形態(tài)的周長和面積。圖1中D，E，F(xiàn)代表在顯微鏡下觀測的石油包裹體二維平面繞左端點旋轉(zhuǎn)不同角度（根據(jù)包裹體Lmax的位置，若Lmax方向垂直x軸，則旋轉(zhuǎn)90°，將包裹體旋轉(zhuǎn)到圖1中F位置）所獲得的石油包裹體形態(tài)，從而來確定Lmin和Lmax.

Lmax/Lmin和C2/（4πA）這2個參數(shù)共同決定了石油包裹體形態(tài)，若石油包裹體的形態(tài)為正圓形，則長寬比（Lmax/Lmin)和不規(guī)則度C2/（4πA）均為1，長寬比越大，表明石油包裹體的形態(tài)越扁，不規(guī)則度越大，表明石油包裹體棱角越分明。

（2）三維體積重構(gòu)石油包裹體二維平面形態(tài)分析是石油包裹體三維體積重構(gòu)的基礎(chǔ)，圖1中石油包裹體F表示二維形態(tài)分析后得到的石油包裹體二維平面圖。從圖1中看出二維分析后的包裹體形態(tài)并沒有變化，發(fā)生變化的是石油包裹體二維形態(tài)Lmax與水平面的夾角。在顯微鏡下測定的原始二維形態(tài)的Lmax與水平面的夾角總是大于或等于0（位置D或位置E），而石油包裹體二維形態(tài)變化后Lmax與水平面的夾角為0（位置F）。石油包裹體二維平面形態(tài)分析完成后就可以通過假定石油包裹體二維平面與三維形態(tài)之間的某種關(guān)系，進而重構(gòu)石油包裹體三維形態(tài)，計算公式為

在水平x方向上，石油包裹體由垂直于Lmax方向的n個直徑為R的圓所組成，R表示石油包裹體垂直Lmax方向的寬度（R≤Lmin）（圖1）。當(dāng)n較大時，可以近似認(rèn)為石油包裹體是由n個高為Lmax/n、底面積為πR2/ 4的圓柱體組成，因此石油包裹體的體積可以通過（1）式計算。

2.2氣相體積計算

石油包裹體中的氣泡可分為規(guī)則和不規(guī)則2類[13]，形狀規(guī)則的氣泡主要為球形或近球形，而形狀不規(guī)則的氣泡主要為紡錘形、酒桶形或其他立體形態(tài)。

當(dāng)氣泡較小時，石油包裹體在各個方向上的長度都比氣泡的直徑大，那么就可以認(rèn)為氣泡是個球體，則可以通過測定二維形態(tài)中氣泡直徑，計算得到氣泡體積；而當(dāng)氣泡較大或石油包裹體形狀極不規(guī)則時，由于受到石油包裹體壁的影響，氣泡表現(xiàn)出紡錘形、酒桶形或其他的不規(guī)則形態(tài)，氣泡體積的計算則需要通過三維重構(gòu)來進行，與石油包裹體體積計算方法類似。

3　計算實例及誤差分析

文獻［14］認(rèn)為包裹體熱力學(xué)模擬精度的主要影響因素有包裹體流體成分、包裹體氣液比精度、古地層壓力計算范圍有限等。因此，只有準(zhǔn)確獲取了石油包裹體的氣液比，才可以較為客觀地進行熱力學(xué)模擬，進而準(zhǔn)確獲取古地層壓力。

為客觀評價本文提出的氣液比計算方法的精度和可行性，選擇了13種不同的石油包裹體形態(tài)（圖2），開展了石油包裹體二維形態(tài)分析、石油包裹體體積計算、氣泡體積計算以及氣液比計算，同時對這13種不同形態(tài)的石油包裹體進行了CLFSM法氣液比測定（CLFSM法的計算結(jié)果可以近似認(rèn)為是真實的氣液比）。

圖2　石油包裹體二維形態(tài)類型示意

表1詳細(xì)列出了13種石油包裹體的形態(tài)參數(shù)（長寬比和不規(guī)則度）、新方法氣液比計算結(jié)果、CLFSM法計算結(jié)果以及2種方法計算結(jié)果的絕對誤差和相對誤差。對比新方法與CLFSM法計算得到的氣液比后不難發(fā)現(xiàn)，二者絕對誤差主體控制在±2%內(nèi)（圖3），但是從表1中看到，受石油包裹體形態(tài)影響，二者的相對誤差分布跨度較大，達到了±50%，表明本文提出的新方法具有一定的適用條件。

為進一步研究新方法的適用性，對造成新方法相對誤差較大的原因進行了分析（圖4），發(fā)現(xiàn)石油包裹體形態(tài)不規(guī)則是造成相對誤差較大的主要原因。石油包裹體形態(tài)越不規(guī)則，則棱角化越明顯，所計算的氣液比相對誤差越大。計算的13種石油包裹體的形態(tài)，主要有長條狀、尖棱狀、次尖棱狀、水滴狀、負(fù)晶型、近規(guī)則狀等6種類型，從這6種類型石油包裹體的相對誤差分布上來看（圖4），水滴狀、尖棱狀（圖2a，圖2l，圖2m）、次尖棱狀以及負(fù)晶形的石油包裹體誤差較大，而近規(guī)則狀（圖2e，圖2g）和大部分長條狀石油包裹體（圖2f，圖2i，圖2k）相對誤差則較小，主體在±10%內(nèi)，因此對于形態(tài)較規(guī)則的石油包裹體，本文介紹的新方法可以較為客觀地計算出氣液比。

表1　不同形態(tài)石油包裹體氣液比計算結(jié)果及誤差對比

圖3　石油包裹體形態(tài)參數(shù)與絕對誤差分布示意

圖4　不同石油包裹體形態(tài)相對誤差分布

基于上述2種方法計算結(jié)果的誤差分析，筆者認(rèn)為在進行石油包裹體熱力學(xué)模擬前，應(yīng)選擇形狀較規(guī)則（長寬比和不規(guī)則度接近1）的石油包裹體氣液比，盡量不要選擇尖棱狀、水滴狀、負(fù)晶形等不規(guī)則形狀石油包裹體氣液比，以避免熱力學(xué)模擬結(jié)果產(chǎn)生較大誤差。

4　結(jié)論

（1）長寬比和不規(guī)則度可以有效表征石油包裹體二維形態(tài)，長寬比越大，表明石油包裹體的形態(tài)越扁；不規(guī)則度越大，表明石油包裹體棱角越分明。

（2）對于形狀不規(guī)則的水滴狀、尖棱狀以及負(fù)晶形石油包裹體，用本文新方法計算氣液比相對誤差較大，但對于形態(tài)較規(guī)則的石油包裹體則可以獲得較為客觀的、相對誤差較小的氣液比。相比目估測量法，新方法更為便捷、實用。

參考文獻：

［1］米敬奎，戴金星，張水昌.含油氣盆地包裹體研究中存在的問題［J］.天然氣地球科學(xué)，2005，16（5）：602-636. MI Jingkui，DAI Jinxing，ZHANG Shuichang.Some problems existed in research of inclusion occurring in oil and gas bearing basin［J］. Natural Gas Geoscience，2005，16（5）：602-636.

［2］劉文斌，姚素平，胡文瑄，等.流體包裹體的研究方法及應(yīng)用［J］.新疆石油地質(zhì)，2003，24（3）：264-267. LIU Wenbin，YAO Suping，HU Wenxuan，et al.Application and re?search approaches of fluid inclusions［J］.Xinjiang Petroleum Geolo?gy，2003，24（3）：264-267.

［3］PIRONON J.Organic inclusions in salt.Part1:solid and liquid or?ganic matter carbon dioxide and nitrogen species in fluid inclusions from the Bresse basin［J］.Organic Geochemistry，1995，23（5）：391-402.

［4］PIRONON J.Organic inclusions in salt.Part2:oil?gas and ammoni?um inclusions from the Gabon margin［J］.Organic Geochemistry，1995，23（8）：739-750.

［5］李艷霞，趙靖舟，李靜紅.鄂爾多斯盆地東部二疊系流體包裹體與天然氣成藏［J］.新疆石油地質(zhì)，2011，32（5）：451-452. LI Yanxia，ZHAO Jingzhou，LI Jinghong.Fluid inclusions and gas ac?cumulation of Permian in eastern Ordos basin［J］.Xinjiang Petro?leum Geology，2011，32（5）：451-452.

［6］李偉，張枝煥，楊永才，等.流體包裹體在準(zhǔn)噶爾盆地油氣成藏研究中的應(yīng)用［J］.新疆石油地質(zhì)，2006，27（5）：600-603. LI Wei，ZHANG Zhihuan，YANG Yongcai，et al.Application of fluid inclusion to research of oil gas accumulations in Junggar basin［J］. XinjiangPetroleum Geology，2006，27（5）：600-603.

［7］王金志，楊少武，蔣森堡，等.流體包裹體熱動力學(xué)模擬技術(shù)的古壓力恢復(fù)方法及應(yīng)注意的問題［J］.中國石油勘探，2008，13（1）：44-47. WANG Jinzhi，YANG Shaowu，JIANG Senbao，et al.Methods for pa?leopressure reconstruction based on thermodynamic simulation tech?nologies of fluid inclusion and some problems needing attention［J］. ChinaPetroleum Exploration，2008，13（1）：44-47.

［8］張鼐，王招明，楊文靜，等.塔中東部奧陶系烴包裹體特征及成藏意義［J］.新疆石油地質(zhì)，2011，32（3）：267-271. ZHANG Nai，WANG Zhaoming，YANG Wenjing，et al.The hydrocar?bon inclusions and significance of hydrocarbon accumulation in Or?dovician reservoirs in eastern Tazhong area，Tarim basin［J］.Xinji?angPetroleum Geology，2011，32（3）：267-271.

［9］王存武，鄒華耀，姜麗娜，等.激光掃描共聚焦顯微鏡精確測量有機包裹體氣液比方法研究［J］.現(xiàn)代科學(xué)儀器，2008，24（1）：20-22. WANG Cunwu，ZOU Huayao，JIANG Li'na，et al.The method study on the determination of gas and liquid ratio of fluid inclusion using laser confocal scanning microscope［J］.Modern Scientific Instru?ments，2008，24（1）：20-22.

［10］PIRONON J，CANALS M，DUBESSY J，et al.Volumetric recon?struction of individual fluid inclusions by confocal scanning laser microscopy［J］.European Journal of Mineralogy，1998，10（6）：1 143-1 150.

［11］SHEPHERD T J，RANKIN A H，ALDERTON D H M.A practical guide to fluid inclusion studies［M］.U.K：Blackie Academic& Professional，1985.

［12］張文淮，陳紫英.流體包裹體地質(zhì)學(xué)［M］.武漢：中國地質(zhì)大學(xué)出版社，1993. ZHANG Wenhuai，CHEN Ziying.Fluid inclusion geology［M］.Wu?han:ChinaUniversity of Geosciences Press，1993.

［13］周振柱，周瑤琪，陳勇，等.一種獲取流體包裹體氣液比的便捷方法［J］.地質(zhì)論評，2011，57（1）：147-152. ZHOU Zhenzhu，ZHOU Yaoqi，CHEN Yong，et al.A convenient method for obtaining vapour/liquid ratios of fluid inclusions［J］. Geological Review，2011，57（1）：147-152.

［14］王存武，鄒華耀.利用流體包裹體獲取含氣盆地古地層壓力的新方法［J］.天然氣勘探與開發(fā)，2013，36（1）：28-32. WANG Cunwu，ZOU Huayao.A new method to obtain paleo forma?tion pressure of gas bearing basin using fluid inclusion［J］.Natural Gas Exploration and Development，2013，36（1）：28-32.

（編輯曹元婷）

A New Method for Calculating Ratio of Vapor to Liquid Phase Volume of Oil Inclusions

DENG Mo
(Wuxi Institute of Petroleum Geology,Research Institute of Petroleum Exploration and Development,Sinopec,Wuxi,Jiangsu 214151,China)

Abstract:Paleopressure can be obtained by oil inclusion thermodynamics modeling,and ratio of vapor to liquid phase volume of oil inclu?sion is a sensitive parameter in this modeling,which can be got by 2D plane analysis and 3D shape reconfiguration of petroleum inclusions. The case study shows that within a certain oil conclusions shape,this ratio obtained by 2D plane analysis and 3D shape reconfiguration has relative error within±10%,compared with the ratio given by confocal laser scanning?fluorescence scanning microscopy(CLFSM).Without CLFSM,the 2D plane analysis and 3D shape reconfiguration method is better than ocular estimation measurement method in calculation of the ratio,and can meet the accurate requirement for trappingpressure modelingof petroleum inclusions.

Keywords:oil inclusion;vapor?liquid ratio;thermodynamics modeling;trappingpressure

作者簡介：鄧模（1983-），男，貴州遵義人，工程師，石油地質(zhì)，（Tel）13815104906（E-mail）dengmo.syky@sinopec.com

收稿日期：2014-12-05

修訂日期：2015-12-22

文章編號：1001-3873（2016）02-0236-04

DOI：10.7657/XJPG20160221

中圖分類號：TE124.1

文獻標(biāo)識碼：A