張 平 李洪烈 李 燦 彭太翀 梁 懌

（國家管網(wǎng)集團西氣東輸公司）

隨著天然氣的貿(mào)易計量方式逐步由體積計量轉(zhuǎn)換為能量計量，天然氣的物性參數(shù)，尤其是密度、熱值、組分計算的準確性對貿(mào)易計量爭議的影響越來越明顯［1，2］。目前天然氣計量的方法是流量計算機通過獲取天然氣的工況流量、溫度及壓力等信息，結(jié)合預設(shè)的氣質(zhì)組分參數(shù)換算成標準狀態(tài)下（20℃，101.325kPa）的天然氣體積，用于和用戶的貿(mào)易交接［3，4］。供氣方提供天然氣日分輸體積總累積量和當日的氣質(zhì)組分表作為交接憑據(jù)。 流量計算機具備根據(jù)在線色譜分析儀提供的天然氣高、低位發(fā)熱量計算日分輸天然氣總能量的能力，這也將是推行天然氣貿(mào)易能量計量的主要指標。 因此，需要一種精確、可靠的算法來根據(jù)已知組分的天然氣測算密度、熱值等物性參數(shù)的方法。

1 研究方法

1.1 目前研究的普遍方法

目前針對氣體的物性參數(shù)計算，基于氣體狀態(tài)方程的算法主要有LKP狀態(tài)方程、SRK狀態(tài)方程和BWR狀態(tài)方程。有學者在關(guān)聯(lián)大量實驗數(shù)據(jù)的基礎(chǔ)上，對BWR方程的應用范圍和精確度做了進一步提高， 形成了修正的BWR 狀態(tài)方程（BWRS）， 這被認為是當前烴類計算中的最佳模型。文獻［5］提出可利用BWRS方程計算單組分烴類氣體和多組分混合氣體的密度和焦耳-湯姆遜系數(shù)，結(jié)果顯示BWRS計算方法簡單、精確度高。但該方法計算精度受比熱容關(guān)系的影響明顯，對于天然氣，需要提供可靠的比熱容關(guān)系式來進一步提高計算精度［5～9］。

其中，p為介質(zhì)絕對壓力，kPa；ρ為介質(zhì)密度，kmol/m3；R為氣體常數(shù)，kJ/（kmol·K）；T為介質(zhì)溫度，K；相關(guān)參數(shù)A0、B0、C0、D0、E0、a′、b、c、d、α′、γ必須通過大量的實驗得到［8］。

從式（1）可以看出，氣體壓力p是以密度ρ和溫度T為自變量的二元函數(shù)。 同樣，基于密度ρ和溫度T的熱力學函數(shù)吉布斯能量方程 （Gibbs Energy）和亥姆霍茲能量方程（Helmholtz Energy）也可用于氣體的物性參數(shù)計算。 但實際研究表明，吉布斯能量方程在介質(zhì)從液態(tài)轉(zhuǎn)換為氣態(tài)過程中并不連續(xù)，較亥姆霍茲能量方程氣液兩相連續(xù)特性的局限性較大，因此筆者采用亥姆霍茲能量方程作為計算基礎(chǔ)。

1.2 亥姆霍茲能量方程

對于某種介質(zhì)而言， 亥姆霍茲能量a通常表示為理想狀態(tài)下能量貢獻a0和能量貢獻修正值ar的總和［10～12］，即：

實際上，亥姆霍茲能量a除以氣體常數(shù)R和溫度T的乘積為一個無量綱的量， 因此自變量T和ρ也可以通過轉(zhuǎn)換，得出亥姆霍茲能量方程的無量綱形式，表述為：

表1 天然氣各組分臨界參數(shù)

對亥姆霍茲能量公式在溫度不變的情況下對體積求導，即可得到介質(zhì)壓力關(guān)于亥姆霍茲能量的顯性方程：

1.3 天然氣密度計算

對于天然氣密度而言， 將已知天然氣的壓力、溫度及各組分系數(shù)等參數(shù)代入式（9）即可求出混合后天然氣的密度。 為便于求解，將式（9）改寫成如下形式進行迭代求解：

2 計算結(jié)果分析

以西氣東輸某分輸站場為例，對該站場連續(xù)7日間氣體組分的體積百分比數(shù)據(jù)進行統(tǒng)計，結(jié)果見表2。

表2 某站場天然氣組分體積百分比 %

按前文方法計算出天然氣的相對密度，再與 色譜分析儀測量的結(jié)果比對，具體數(shù)據(jù)見表3。

表3 天然氣相對密度計算結(jié)果

由表3中的計算結(jié)果可知， 與色譜分析儀測量結(jié)果相比筆者所提方法計算出來的天然氣相對密度值在正常范圍之內(nèi)， 誤差在-0.000 6～0.000 3之間，誤差百分比不超過0.10%。 因此，該方法用于計算站場天然氣密度值是準確、 可行的。

3 結(jié)束語

筆者給出了一種基于氣體亥姆霍茲能量方程和氣體物性參數(shù)關(guān)系計算天然氣密度的方法。在既定的工況條件下，天然氣密度的計算結(jié)果精度較高，能夠滿足生產(chǎn)需求，可較好地作為一種輔助工具進行天然氣密度的計算，在天然氣能量計量的應用中起到了關(guān)鍵的作用，同時也為天然氣其他物性參數(shù)的軟測量提供了技術(shù)思路。