張 雋，張家利，席斌斌，施偉軍，蔣 宏

(1.中國石化石油勘探開發(fā)研究院無錫石油地質(zhì)研究所，江蘇無錫 214151;2.內(nèi)蒙古自治區(qū)一一五地質(zhì)礦產(chǎn)勘查開發(fā)院，內(nèi)蒙古烏蘭浩特 137400)

儲層中CH4－H2O－NaCl體系包裹體廣泛存在，該體系包裹體的熱力學(xué)參數(shù)是了解油氣成藏期次、探討成藏過程的重要依據(jù)[1－5］。國內(nèi)外學(xué)者對該體系的熱力學(xué)性質(zhì)以及PVTx狀態(tài)方程進(jìn)行了研究，基本建立了熱力學(xué)模型并編制了相應(yīng)的熱力學(xué)計算軟件[6－12］。但是由于上述熱力學(xué)方程涉及的參數(shù)較多、計算過程復(fù)雜，且編制的軟件難以直接利用包裹體的測溫數(shù)據(jù)或激光拉曼分析數(shù)據(jù)進(jìn)行計算，因而限制了這些熱力學(xué)計算軟件在地質(zhì)應(yīng)用中的推廣。同時，上述包裹體熱力學(xué)計算軟件開發(fā)環(huán)境及開發(fā)語言多采用面向過程的程序開發(fā)模式，面對包裹體研究不斷完善的熱力學(xué)計算模型及日益增強(qiáng)的圖形功能的需求，軟件的升級往往“牽一發(fā)而動全身”，維護(hù)及移植性較差，使得軟件體現(xiàn)包裹體研究的進(jìn)展受限。基于此，本文在Visual Studio 2010平臺上，采用面向?qū)ο蟮某绦蛟O(shè)計模式，利用微軟的MFC類庫，在迭代法計算含CH4鹽水包裹體古壓力的算法基礎(chǔ)上[13］，開發(fā)出了CH4－H2O－NaCl體系包裹體PVTx計算模擬軟件，以期能更好、更方便地實現(xiàn)CH4－H2O－NaCl體系包裹體的熱力學(xué)計算與地質(zhì)應(yīng)用。

1 計算原理及步驟

迭代法計算包裹體的PVTx性質(zhì)，在CO2－H2ONaCl體系中得到了比較多的應(yīng)用[14－18］，但是由于CH4的臨界溫度(－82.6℃)遠(yuǎn)低于CO2的臨界溫度(31.0℃)，因此無法根據(jù)部分均一溫度和均一方式來計算部分均一時包裹體的一些熱力學(xué)參數(shù)，故無法將CO2－H2O－NaCl體系的計算方法直接移植到 CH4－H2O－NaCl體系中。為了解決上述問題，席斌斌等[13］提出了利用室溫條件下CH4的拉曼位移值與壓力的關(guān)系式[19－21］，來計算室溫條件下CH4－H2O－NaCl體系的一些熱力學(xué)參數(shù)，進(jìn)而建立了迭代計算CH4－H2O－NaCl體系包裹體均一壓力、組成等PVTx參數(shù)的方法。

具體的計算步驟為:(1)根據(jù)室溫下CH4的拉曼位移值計算包裹體的內(nèi)壓;(2)利用冰點溫度計算包裹體的鹽度;(3)利用Duan等[8］的狀態(tài)方程分別計算室溫條件下包裹體中氣液相的摩爾體積以及液相的組成;(4)假設(shè)室溫條件下包裹體的氣液比為F;(5)計算包裹體的總組成以及總摩爾體積;(6)根據(jù)均一溫度以及步驟(5)中計算的包裹體總組成，利用Duan等[8］的方程計算再次包裹體的總摩爾體積;(7)若步驟(5)和(6)兩次計算的總摩爾體積相等或在誤差范圍之內(nèi)，則證明設(shè)定的F值即為包裹體的氣液比，反之改變F值，重復(fù)步驟(5)、(6)直到兩者相等或在誤差范圍之內(nèi);(8)利用Duan等[8］的方程計算包裹體均一壓力。

由上述論述可知，該計算方法所需要的參數(shù)包括室溫下CH4的拉曼位移、均一溫度、冰點溫度，上述參數(shù)均可以通過顯微測微以及激光拉曼分析直接測定。除了能得到包裹體的均一壓力外，還可以得到包裹體組成、摩爾體積、氣液比等熱力學(xué)參數(shù)。

2 開發(fā)環(huán)境及核心算法

2.1 MFC應(yīng)用程序基本框架及類庫

Visual Studio 2010是一個強(qiáng)大的軟件編程平臺，在其基礎(chǔ)之上的 MFC(Microsoft Foundation Class)是微軟公司提供的類庫，以C++類的形式封裝了Windows的API函數(shù)和控件。在程序的開發(fā)過程中，使用MFC類庫可以使應(yīng)用程序開發(fā)變得簡單，同時提高了代碼的可靠性和可移植性，在不減少編程的自由和靈活性的情況下提供了更多的支持。

基于MFC類庫的應(yīng)用程序可以利用App-Wi-zard(應(yīng)用程序向?qū)?快速建立起應(yīng)用程序的基本框架，利用ClassWizard(類向?qū)?快速地在應(yīng)用程序框架上添加類(或子類)、成員變量及成員函數(shù)。在程序員編寫代碼的過程中，可以利用強(qiáng)大的MFC 類庫，完成大多數(shù)常用的功能[22－26］。

本軟件主界面程序中涉及的類主要繼承于MFC的部分基類，其中CDxCalcWnd類、CDxCurveWnd類以及CPropertiesWnd類是關(guān)鍵的3個類，它們分別定義了軟件計算面板、曲線面板以及屬性面板的顯示功能。MFC是一個龐大的類庫，目前，大約有200多個類，且其中相當(dāng)多類具不同功能。本軟件用到的只是MFC應(yīng)用程序框架的冰山一角。

2.2 迭代二分法

本軟件算法主要根據(jù)二分迭代找到符合計算精度的數(shù)值，從而求解，即:對于給定值x，從序列的中間位置開始比較，如果當(dāng)前位置值等于x，則查找成功;若x小于當(dāng)前位置值，則在數(shù)列的前半段中查找;若x大于當(dāng)前位置值，則在數(shù)列的后半段中繼續(xù)查找，直到找到并輸出為止。此算法貫穿整個軟件的關(guān)鍵數(shù)據(jù)計算始終，諸如計算摩爾體積Vr以及包裹體氣液比F。

根據(jù)室溫下包裹體中氣相摩爾體積Vm－v、室溫下包裹體中液相摩爾體積Vm－l、室溫下包裹體液相中CH4的質(zhì)量摩爾濃度mCH4以及包裹體鹽度mNaCl計算室溫下包裹體氣液比F，假設(shè)F的取值區(qū)間是[0，1］，且要求計算精度為萬分之一(以此確定迭代步長)，那么，計算F的算法流程如圖1執(zhí)行。下面列舉計算F的部分程序代碼:

圖1 迭代二分法算法流程Fig.1 Flow chart of iterative dichotomy algorithm

圖2 軟件功能模塊Fig.2 Functional module diagram of software

3 軟件開發(fā)

3.1 軟件的功能模塊

CH4－H2O－NaCl體系包裹體 PVTx計算模擬軟件涉及一個主界面程序Bleb interface，一個動態(tài)鏈接庫Calculate DLL(封裝所有算法)和一個反注冊程序(清除本軟件在安裝及運(yùn)行時在系統(tǒng)注冊表里寫入的信息)。該軟件實現(xiàn)了2個主要功能:一是精確計算包裹體的均一壓力、氣液比、CH4的摩爾含量等熱力學(xué)參數(shù);二是實現(xiàn)了該體系的PVTx動態(tài)曲線的模擬功能(圖2)。

本軟件應(yīng)用程序主界面采用常規(guī)的Windows界面，包括菜單欄、工具欄、狀態(tài)欄、操作區(qū)以及視圖區(qū)等，該功能主要利用MFC應(yīng)用程序框架向?qū)?圖3)。

動態(tài)鏈接庫是為了實現(xiàn)多個應(yīng)用程序共享數(shù)據(jù)和代碼的一種數(shù)據(jù)文件訪問方式。它獨立于可執(zhí)行文件，在程序原有函數(shù)的參數(shù)和返回值等屬性不變的前提下，僅改變動態(tài)鏈接庫中的函數(shù)或現(xiàn)有函數(shù)的功能即可完成對程序文件的修改，不需重新編譯。本軟件涉及的動態(tài)鏈接庫Calculate DLL，封裝了2個程序文件，主要是計算 CH4－H2O－NaCl體系包裹體的各類參數(shù)值，主界面在計算面板里的“計算”按鈕事件以及曲線面板“擬合”按鈕事件里，會調(diào)用動態(tài)鏈接庫。

軟件適用的范圍為:室溫T=18～30℃，均一溫度Th=0～250℃，均一壓力Ph=0～2 000 bar，鹽度mNaCl=0～6 mol/kg。

3.2 計算功能

圖3 PVTx計算模擬軟件界面Fig.3 Soft ware interface of PVTx calculating simulation

圖4 PVTx計算模擬軟件的“計算面板”界面Fig.4 Computing interface of PVTx calculating simulation software

軟件計算功能主要通過計算面板(圖4)實現(xiàn)。計算面板主要實現(xiàn)從鍵盤輸入室溫T、包裹體均一溫度Th、包裹體冰點溫度Tice以及包裹體拉曼位移Vp的具體實驗檢測數(shù)據(jù)，通過多個熱力學(xué)方程聯(lián)立求解，計算輸出CH4－H2O－NaCl體系包裹體均一壓力Ph，同時，也可以得到計算過程中多個相關(guān)參數(shù)的值，包括:室溫下包裹體中氣相的摩爾體積Vm－v、室溫下包裹體中液相的摩爾體積 Vm－l、室溫下包裹體氣液比F、室溫下包裹體內(nèi)壓P、包裹體鹽度mNaCl、室溫下液相包裹體中CH4的質(zhì)量摩爾濃度mCH4、均一溫度下包裹體摩爾體積V總以及均一溫度下包裹體中CH4的質(zhì)量摩爾濃度mCH4。計算精度可達(dá)萬分之一。

3.3 曲線擬合功能

曲線面板主要是通過曲線控件實現(xiàn)CH4－H2ONaCl體系包裹體的P－T－mCH4數(shù)值動態(tài)模擬功能。其中，P、T、mCH4參數(shù)之間的曲線擬合以計算面板中涉及到的熱力學(xué)方程為依據(jù)，其值都是通過相關(guān)的數(shù)學(xué)函數(shù)f(P，T)=mCH4計算得到。曲線面板的主要功能包括:

(1)通過函數(shù)f(P，T)=mCH4進(jìn)行二維曲線擬合，即:已知P、T、mCH4任意一個參數(shù)的數(shù)值，可以得到其他2個參數(shù)之間的二維光滑曲線關(guān)系圖。

(2)曲線圖可以疊加顯示，圖例可個性化定制，明確區(qū)分不同參數(shù)數(shù)據(jù)下的曲線特征，藉以探索P、T、mCH4參數(shù)間的關(guān)系。

(3)曲線可以任意放大縮小，且其中的參數(shù)值可以精確顯示到萬分之一。其工作原理是:計算機(jī)屏幕的像素是有限的，圖像的顯示就是把組成圖像的數(shù)據(jù)按照一定的規(guī)則分割，并依次投點到對應(yīng)的像素中。雖然一個像素可以包含若干個數(shù)據(jù)，但只能有一個被顯示出來。假設(shè)一個屏幕像素里可以“裝”X個數(shù)據(jù)，當(dāng)用戶在曲線面板界面控制“精度”工具時，放大一倍，實際上就是讓一個像素所承載的數(shù)據(jù)減小，這樣，原來一個像素“裝”了10個數(shù)據(jù)，放大一倍后，一個像素只“裝”了5個數(shù)據(jù)，在屏幕上，用戶就比原來多看到5個數(shù)據(jù)，感覺上曲線(或圖像)放大了。同理，縮小曲線(或圖像)，就是讓一個像素“裝”的數(shù)據(jù)增加。

4 計算實例

為了對軟件的使用效果進(jìn)行檢驗，對采自普光5井5 157.7 m、上二疊統(tǒng)長興組的粒屑白云巖中的包裹體進(jìn)行了應(yīng)用。包裹體主要發(fā)育于溶孔充填的方解石中，以氣液兩相包裹體和富氣包裹體為主，對其中一個氣液兩相鹽水包裹體進(jìn)行了顯微測溫以及激光拉曼分析。測試表明包裹體的均一溫度為190.3℃，冰點溫度為－7.6℃，室溫(24℃)條件下CH4的拉曼位移為2 916.86 cm－1。將上述測試結(jié)果輸入軟件的計算面板進(jìn)行計算，得到包裹體的均一壓力為655.1 bar，包裹體中CH4的摩爾質(zhì)量濃度為0.342 mol/kg，NaCl的摩爾質(zhì)量濃度為2.16 mol/kg。利用曲線面板，可以得到該組成條件下，氣—液飽和曲線，并能夠?qū)⒂嬎憬Y(jié)果直觀地表示在P－T相圖中(圖5)。

上述計算結(jié)果表明，包裹體形成于較高的溫壓條件下，包裹體捕獲于氣—水不混溶體系，該期包裹體可能代表了氣藏形成時的PVTx條件。上述結(jié)果與秦建中等[27］所提出的第三期即原油和烴源巖熱裂解生成干氣階段的特征基本一致，表明計算結(jié)果具有較高的可信性。

結(jié)果證明，微軟的MFC類庫及其應(yīng)用程序框架編制的軟件，計算快速準(zhǔn)確，界面友好，拓展性強(qiáng)，在地質(zhì)領(lǐng)域具有廣闊的應(yīng)用前景。

圖5 P－T相圖Fig.5 Phase diagram between pressure and temperature

[1]Lu Wanjun，Chou I-Ming，Burruss R C，et al.A unified equation for calculating methane vapor pressures in the CH4－H2O system with measured Raman shifts[J］.Geochimica et Cosmochimica Acta，2007，71(16):3969－3978.

[2]Xu G J.Methane-rich fluid inclusions in the Proterozoic Zn –Pb–Ag deposit at Dugald River，NW Queensland:potential as an exploration guide[J］.Applied Geochemistry，2000，15(1):1－12.

[3]殷和平，錢一雄，陳強(qiáng)路，等.流體包裹體主要成分及譜學(xué)特征研究:以塔河油區(qū)下奧陶統(tǒng)儲層為例[J］.石油實驗地質(zhì)，2009，31(3):282－286.Yin Heping，Qian Yixiong，Chen Qianglu，et al.Approaching to chemical composition and application of Raman spectrum of fluid inclusion:taking the Lower Ordovician reservoirs in Tahe as an example[J］.Petroleum Geology ＆ Experiment，2009，31(3):282－286.

[4]施偉軍，蔣宏，席斌斌，等.油氣包裹體成分及特征分析方法研究[J］.石油實驗地質(zhì)，2009，31(6):643－648.Shi Weijun，Jiang Hong，Xi Binbin，et al.Studies of analysis approaches of oil-and-gas inclusion composition and characteristics[J］.Petroleum Geology ＆ Experiment，2009，31(6):643－648.

[5]劉德漢，戴金星，肖賢明，等.普光氣田中高密度甲烷包裹體的發(fā)現(xiàn)及形成的溫度和壓力條件[J］.科學(xué)通報，2010，55(3/4):359－366.Liu Dehan，Dai Jinxing，Xiao Xianming，et al.High density methane inclusions in Puguang gasfield:discovery and a T－P genetic study[J］.Chinese Science Bulletin，2010，55(3/4):359－366.

[6]Duan Z H，M?ller N，Greenberg J，et al.The prediction of methane solubility in natural waters to high ionic strength from 0℃to 250℃ and from 0 to 1600 bar[J］.Geochimica et Cosmochimica Acta，1992，56(4):1451－1460.

[7]Duan Z H，M?ller N，Weare J H.An equation of state EOS for CH4，CO2and H2O[J］.Geochimica et Cosmochimica Acta，1992，56:2605－2618.

[8]Duan Z H，Mao S.A thermodynamic model for calculating methane solubility，density and gas phase composition of methane-bearing aqueou.s fluids from 273 to 523 K and from 1 to 2000 bar[J］.Geochimica et Cosmochimica Acta，2006，70(13):3369－3386.

[9]Aasberg－Petersen K，Stenby E，F(xiàn)redenslund A.Prediction of high-pressure gas solubilities in aqueous mixtures of electrolytes[J］.Industrial and Engineering Chemistry Research，1991，30(9):2180－2185.

[10]Amirijafari B，Campbell J M.Solubility of gaseous hydrocarbon mixtures in water[J］.Society of Petroleum Engineers Journal，1972，12(1):21－27.

[11]Ashmayan K D，Skripka V G，Namiot A Y.The solubilities of methane and nitrogen in water at high temperatures and pressures[J］.Geochem Int，1984，2(4):138－139.

[12]Avlonitis D，Danesh A，Todd A C.Prediction of VL and VLL equilibria of mixtures containing petroleum reservoir fluids and methanol with a cubic EOS[J］.Fluid Phase Equilibria，1994，94:181－216.

[13]席斌斌，施偉軍，蔣宏.迭代法計算CH4－H2O－NaCl體系包裹體古壓力及其應(yīng)用[C］//第十三屆全國有機(jī)地球化學(xué)學(xué)術(shù)會議.2011:774－775.Xi Binbin，Shi Weijun，Jiang Hong.Application of iterative method for calculating homogenization pressure of H2O－CO2－NaCl inclusion system[C］//The Thirteenth National Symposium of organic geochemistry.2011:774－775.

[14]宋玉財，胡文瑄，倪培，等.含CO2鹽水流體包裹體摩爾體積和組分求解新方法[J］.中國科學(xué):D輯:地球科學(xué)，2007，37(3):363－369.Song Yucai，Hu Wenxuan，Ni Pei，et al.New method for solving CO2containing saline fluid inclusion molar volume and composition[J］.Science in China:Series D:Earth Sciences，2007，37(3):363－369.

[15]劉斌，沈昆.流體包裹體熱力學(xué)[M］.北京:地質(zhì)出版社，1999.Liu Bin，Shen Kun.Fluid inclusion thermodynamics[M］.Beijing:Geological Publishing House，1999.

[16]劉斌，朱思林，沈昆.流體包裹體熱力學(xué)參數(shù)計算軟件及算例[M］.北京:地質(zhì)出版社，2000.Liu Bin，Zhu Silin，Shen Kun.Fluid inclusions and example calculation software[M］.Beijing:Geological Publishing House，2000.

[17]Parry W T.Estimation of XCO2，P and fluid inclusion volume from fluid inclusion temperature measurements in the system NaCl－H2O－CO2[J］.Economic Geology，1986，81(4):1009－1013.

[18]席斌斌，施偉軍，張德會，等.迭代法計算H2O－CO2－NaCl包裹體均一壓力的改進(jìn)及其應(yīng)用[J］.礦床地質(zhì)，2010，29(6):1138－1144.Xi Binbin，Shi Weijun，Zhang Dehui，et al.Improvements and application of iterative method for calculating homogenization pressure of H2O－CO2－NaCl inclusion system[J］.Mineral Deposits，2010，29(6):1138－1144.

[19]Ben－Amotz D，La Plant F，Shea D，et al.Raman studies of molecular potential energy surface changes in supercritical fluids[C］//Supercritical Fluid Technology;ACS Symposium Series 488.Washington，DC:American Chemical Society，1999:18－30.

[20]Chen Y，Zhou Y Q，Liu C Y，et al.Quantitatively analyzing the homogenization process of CH4－H2O fluid inclusion by laser Raman spectroscopy[J］.Earth Science Frontiers，2005，12(4):592－596.

[21]Chou I M，Sharma A，Burruss R C，et al.Diamond－anvil cell observations of a new methane hydrate phase in the 100－MPa pressure range[J］.Journal of Physical Chemistry A，2001，105(19):4664－4668.

[22]侯捷.深入淺出 MFC[M］.武漢:華中科技大學(xué)出版社，2000.Hou Jie.Stusy MFC[M］.Wuhan:Huazhong University of Science and Technology Press，2000.

[23]譚浩強(qiáng).C++面向?qū)ο蟪绦蛟O(shè)計[M］.北京:清華大學(xué)出版社，2006.Tan Haoqiang.Object oriented programming of C++[M］.Beijing:Tsinghua University Press，2006.

[24]李永中.基于VC的多線程編程技術(shù)及其實現(xiàn)[J］.微計算機(jī)信息，2002(7):26－28.Li Yongzhong.Multi thread programming technology based on VC and its implementation[J］.Micro Computer Information，2002(7):26－28.

[25]袁媛.基于MFC類設(shè)計ObjectARX應(yīng)用程序界面[J］.中南大學(xué)學(xué)報，2004，35(3):448－452.Yuan Yuan.MFC sequence-based design for ObjectARX application interface[J］.Journal of Central South University，2004，35(3):448－452.

[26]Scot D J K，Shepherd W G.VC++6.0 技術(shù)內(nèi)幕[M］.北京:希望電子出版社，2000.Scot D J K，Shepherd W G.VC++6.0 technology insider[M］.Beijing:Hope Electronic Press，2000.

[27]秦建中，孟慶強(qiáng)，付小東.川東北地區(qū)海相碳酸鹽巖三期成烴成藏過程[J］.石油勘探與開發(fā)，2008，35(5):548－556.Qin Jianzhong，Meng Qingqiang，F(xiàn)u Xiaodong.Three hydrocarbon generation and accumulation processes of marine carbonate rocks in northeastern Sichuan Basin，China[J］.Petroleum Exploration and Development，2008，35(5):548－556.