用孔板流量計(jì)測量天然氣時(shí)，孔板流量計(jì)的準(zhǔn)確度除與節(jié)流裝置有關(guān)外，還與所采用的標(biāo)準(zhǔn)和天然氣物性參數(shù)的計(jì)算有關(guān)。目前，關(guān)于使用孔板流量計(jì)測量天然氣流量的標(biāo)準(zhǔn)有SYL04-83和AGA 3 ，這兩個(gè)標(biāo)準(zhǔn)在有關(guān)計(jì)算上有統(tǒng)一的地方，也存在著分歧，本文結(jié)合兩個(gè)標(biāo)準(zhǔn)討論一下天然氣參數(shù)計(jì)算和流量計(jì)算。

一 天然氣的物性參數(shù)計(jì)算孔板流量計(jì)與天然氣流量測量

使用孔板流量計(jì)測量天然氣流量時(shí)，用到的天然氣物性參數(shù)較多，SYL04-83,AGA 3均以不同的形式提供和推薦了天然氣物性參數(shù)和其計(jì)算方法，SYL04-83提供的數(shù)據(jù)和方法已不是的。為更好的使用SYL04-83，提高天然氣流量測量準(zhǔn)確度，以下討論幾個(gè)流量計(jì)計(jì)算影響較大的物性參數(shù)的計(jì)算方法。

1. 天然氣中單一組分的物性參數(shù)

SYL04-83和AGA 3 關(guān)于單一組分的物性參數(shù)與近期相關(guān)的專業(yè)發(fā)表的數(shù)據(jù)有些出入，表1是筆者根據(jù)AGA3（1985年版）和IUPAC 87 相對原子質(zhì)量整理的，建議在SYL04-83修訂之前協(xié)議的使用表1。

2. 天然氣的假臨界參數(shù)的計(jì)算

SYL04-83 以示例的形式介紹了用開氏規(guī)則計(jì)算天然氣假臨界壓力Pe和假臨界溫度Te的方法，當(dāng)天然氣含有H2S和CO2時(shí)，這種方法是有偏差的，此時(shí)應(yīng)采用Wichert—Aziz公式計(jì)算。其公式為：

式中 Te----修正后假臨界溫度 K

Pe---修正后假臨界壓力 MPa

E----溫度修正系數(shù)，E=66.667（A0.9-A1.5）+8.335（B0.5+B4）其中B為H2S的摩爾組分，A=B+X，X為CO2的摩爾組分。

3. 天然氣的密度和相對密度計(jì)算孔板流量計(jì)與天然氣流量測量

天然氣的理想氣體密度和相對密度，可以參照SYL04-83計(jì)算，空氣的摩爾質(zhì)量mair應(yīng)取為28.9625。天然氣的真實(shí)氣體的密度和相對密度應(yīng)采用GB1106-89《天然氣發(fā)熱量，密度和相對密度的計(jì)算方法》計(jì)算。

4 天然氣的等熵指數(shù)計(jì)算

SYL04-83規(guī)定一般k=1.3,這是不的，目前，關(guān)于k值的經(jīng)計(jì)算論文較少，本文結(jié)合化工行業(yè)的手冊類出版物，提出一種計(jì)算等熵指數(shù)k值的方法，供參考使用。

理想氣體的比定壓摩爾熱容一般認(rèn)為只是溫度的函數(shù)，可用下式計(jì)算：

CP=A0+A1T+A2T2+A3T2 （2）

式中： C―――理想氣體比定壓摩爾熱容 kj/mol *K;

T―――氣體溫度 K

A0 ， A1 ， A2 ， A3 為溫度系數(shù) 詳見表2

表1 天然氣中各種氣體的物性參數(shù)

天然氣的單一組分理想氣體比定壓摩爾熱容可由上式求得，混合氣的理想氣體比定壓摩爾熱容可用開氏規(guī)則計(jì)算，其真實(shí)氣體的比定壓摩爾熱容可以由下式計(jì)算：

CP=CP+CP （3）

式中 CP ――真實(shí)氣體比定壓摩爾熱容 kJ/mol*k

CP ――由開氏規(guī)則計(jì)算的理想氣體比定壓摩爾熱容 k J/mol*K

CP ―― CP的壓力修正量 KJ/mol*K

CP值可由Prry-Chilton 1973版本查出，為計(jì)算方便，筆者

用zui小二乘法將其擬合為以下的方程（0.02《 Pt《2 》

CP=APr+B （4）

式中 A, B ----系數(shù)，詳見表3

P―――天然氣的假對比壓力

天然氣的真實(shí)氣體比定壓摩爾熱容CV可由CV=CP-（CP-CV）求得 CP-CV

按下式計(jì)算：

式中： Tr------天然氣的假對比溫度

R------氣體普適常數(shù)8.314510 kj/mol *K

表2 各種氣體的溫度系數(shù)

表3 天然氣假對比溫度及系數(shù)

天然氣的k值并不嚴(yán)格等于 ，還應(yīng)對比值進(jìn)行真實(shí)氣體修正，但對于計(jì)算膨脹系數(shù)而言，可以不再修正，按下式計(jì)算：

5 天然氣的壓縮系數(shù)和超壓縮系數(shù)計(jì)算

超壓縮系數(shù)的計(jì)算偏差對流量計(jì)算影響較大。SYL04-83以附表的形式給出了壓縮系數(shù)和超壓縮系數(shù)，使用時(shí)需要內(nèi)插，認(rèn)為誤差較大;AGA 8被認(rèn)為是可行的計(jì)算方法，但計(jì)算工作量較大，不易實(shí)現(xiàn)。本文介紹了幾種計(jì)算工作量小的壓縮系數(shù)計(jì)算方法，供使用SYL04-83時(shí)參考。

1. 太平洋能源協(xié)會方程

太平洋能源協(xié)會方程是一種由相對密度計(jì)算壓縮系數(shù)的方法，這與AGA3（1985年版）推薦的壓縮系數(shù)在形式上是一致的。它的特點(diǎn)是計(jì)量工作量小，對非交接計(jì)量一般認(rèn)為是可以接收的，其方程為：

式中 Z―――天然氣壓縮系數(shù)

Pa―――天然氣的壓力 MPa

G―――天然氣的相對密度

K1 ， K2―――系數(shù)，詳見表4。

表4 天然氣相對密度及系數(shù)

2. 霍一雅方程

霍一雅方程適用于Tr 》1.05 ，Pr較大時(shí)對比狀態(tài)，此時(shí)霍一雅方程有+-0.3%準(zhǔn)確度。其方程為：

3. 雷特利奇――孔方程

對于甲烷含量較高的天然氣pt<0.6時(shí)，可用以下方程求得Z值：

對方程求導(dǎo)后，用牛頓迭代可求得Z值。

SYL04-83定義 為超壓縮系數(shù)，AGA 3定義為超壓縮系數(shù)，并定義Zf/Zb為可壓縮性，這兩個(gè)標(biāo)準(zhǔn)對超壓縮系數(shù)的定義是不一致的。天然氣的具有代表性混合物的Z值為0.9900～0.9999，并不等于1。SYL04-83推薦的壓縮系數(shù)和超壓縮系數(shù)若是在考慮Zb后得出。兩個(gè)標(biāo)準(zhǔn)就無太大差異，若沒有考慮Zb，兩者是有差異的，計(jì)算超壓縮系數(shù)時(shí)應(yīng)按AGA 3的方法計(jì)算，以便減少由Zb帶來。

由于粘度對流量系數(shù)的影響不大，粘度的計(jì)算不討論了，但應(yīng)強(qiáng)調(diào)指出計(jì)算雷諾數(shù)時(shí)應(yīng)盡量地計(jì)算粘度。

二 流量系數(shù)的計(jì)算和流量公式孔板流量計(jì)與天然氣測量

AGA 3 SYL04-83給了相互一致的流量系數(shù)公式，AGA 3（1990年版）給出了流出系數(shù)和速度接近系數(shù)的定義，并以其乘積代替了流量系數(shù)。有些數(shù)據(jù)應(yīng)用了新近研究的結(jié)果，在SYL04-83沒修訂前應(yīng)允許交接（特別是涉外交接）雙方協(xié)議的使用AGA3（1990版）中流出系數(shù)和速度系數(shù)接近系數(shù)公式。

AGA 3和SYL04-83分別給出了流量公式，由于超壓縮系數(shù)的定義不一致，采用理想氣體相對密度計(jì)量流量的公式中Zb這個(gè)因素，使用時(shí)注意這一點(diǎn)。SYL04-83在推薦導(dǎo)常數(shù)因子項(xiàng)時(shí)，使用了空氣的密度1.2046，并且由此得出常數(shù)因子Ab=0.0003619（靜壓以MPa為單位時(shí)為0.3619），推導(dǎo)常數(shù)因此，空氣的標(biāo)準(zhǔn)狀態(tài)理想氣體密度應(yīng)為1.204004，此時(shí)Ah=0.00036202（0.36202）。使用SYL04-83時(shí)注意這一問題。

使用孔板流量計(jì)測量天然氣流量時(shí)，應(yīng)地計(jì)算天然氣物性參數(shù)，盡量使用數(shù)據(jù)，提高天然氣流量測量準(zhǔn)確度。