姬安召，王玉風(fēng)，張光生

（隴東學(xué)院，甘肅慶陽 745100）

天然氣偏差因子是在一定溫度壓力下天然氣體積與該氣體在相同條件下按理想氣體定律計(jì)算的氣體體積的比值。天然氣的偏差因子是天然氣性質(zhì)研究過程中一個(gè)很重要的參數(shù)，例如在密度、黏度、等溫壓縮系數(shù)以及體積系數(shù)計(jì)算過程中必須考慮的參數(shù)之一。該參數(shù)主要決定于氣體的組成、溫度以及壓力。天然氣的偏差因子雖然可以通過實(shí)驗(yàn)進(jìn)行測試，但實(shí)驗(yàn)的周期長，成本高，對實(shí)驗(yàn)儀器的精度要求高。但也可以通過相關(guān)的計(jì)算方法得到。目前，天然氣偏差因子獲取大致有三類：圖版插值法、基于狀態(tài)方程法和圖版擬合法的數(shù)值計(jì)算方法。

在基于狀態(tài)方程計(jì)算天然氣偏差因子的研究方面，對狀態(tài)方程的維里方程法中的AGA8 方法進(jìn)行了研究，給出了該方法的適應(yīng)性評價(jià)[1]。根據(jù)天然氣氣體組成[2]，研究了天然氣壓縮因子的計(jì)算，得出對低中壓含碳?xì)怏w，AGA8 方程計(jì)算誤差最小，而含硫氣體，Piper-Mahmoud 方法誤差最小。根據(jù)實(shí)驗(yàn)研究結(jié)果[3]，對于超高壓氣藏而言，現(xiàn)有的圖版法和經(jīng)驗(yàn)公式對于天然氣偏差因子計(jì)算的誤差很大，甚至有些超壓氣藏天然氣偏差因子誤差超過20%。對圖版擬合的經(jīng)驗(yàn)公式和基于狀態(tài)方程兩參數(shù)法計(jì)算天然氣偏差因子的公式進(jìn)行評價(jià)，給出適合超高壓氣藏天然氣偏差因子的計(jì)算模型。根據(jù)氣體狀態(tài)方程的BWRS 維里方程[4]，采用非線性回歸分析，對該方程指數(shù)參數(shù)進(jìn)行修正，提出了更高精度天然氣偏差因子的計(jì)算方法?；赟trobridge 方法[5]，對BWR（Benedict-Webb-Rubin）狀態(tài)方程進(jìn)行了修正，修正后的方法計(jì)算天然氣偏差因子的精度控制在1%以內(nèi)。根據(jù)Nishiumi-Saito 狀態(tài)方程[6]，采用多元回歸方程，給出了天然氣偏差因子計(jì)算的解析表達(dá)式。但這些方法都要基于不同的狀態(tài)方程，還要進(jìn)行復(fù)雜的數(shù)學(xué)計(jì)算，計(jì)算的結(jié)果在精度評定時(shí)，基本要參照Standing-Katz 圖版的實(shí)驗(yàn)數(shù)據(jù)。

基于圖版擬合方法，在天然氣偏差因子的研究方面，對DAK、HTP、DRP、LXF、HY、Sarem、AGA8-92DC 七種方法計(jì)算天然氣偏差因子的結(jié)果進(jìn)行精度評估[7-10]，給出了每種方法的適用條件。對DAK 方法進(jìn)行了修正[11]，修正后的方法提高了天然氣偏差因子的計(jì)算精度。根據(jù)DPR 方法計(jì)算天然氣偏差因子的方法[12]，采用一階中心差分方法計(jì)算了超高壓氣藏天然氣等溫?cái)M對比壓縮系數(shù)，拓展了Trube 使用范圍。利用低壓傳統(tǒng)天然氣偏差因子圖版數(shù)據(jù)和實(shí)驗(yàn)測試的高壓天然氣數(shù)據(jù)進(jìn)行擬合[13]，得到了適合于中低壓和高壓的天然氣偏差因子計(jì)算GXH 模型。這些圖版擬合結(jié)果，基本都是基于Standing-Katz 圖版數(shù)字化后的結(jié)果，或者有些參照了狀態(tài)方程的基本形式，其不同的模型都用一定的適用條件。

除了基于狀態(tài)方程法和圖版擬合法之外，還有采用BP 神經(jīng)網(wǎng)絡(luò)進(jìn)行預(yù)測方法，以測試壓力、溫度、天然氣的組成以及其臨界參數(shù)、測試的氣體的偏差因子作為輸入[14]，采用BP 神經(jīng)網(wǎng)絡(luò)方法對249 個(gè)實(shí)驗(yàn)數(shù)據(jù)點(diǎn)進(jìn)行了樣本學(xué)習(xí)，建立基于BP 神經(jīng)網(wǎng)絡(luò)的天然氣偏差因子的預(yù)測模型，但該模型對輸入數(shù)據(jù)的依賴性強(qiáng)，實(shí)用性不強(qiáng)?；跉獠匚镔|(zhì)平衡原理進(jìn)行天然氣偏差因子計(jì)算的模型[15]，根據(jù)生產(chǎn)動態(tài)資料，結(jié)合干氣氣藏的物質(zhì)平衡方程，提出了基于生產(chǎn)資料計(jì)算天然氣偏差因子的新方法，但計(jì)算的精度依賴于氣藏特性和生產(chǎn)動態(tài)資料的準(zhǔn)確性。在商業(yè)軟件進(jìn)行天然氣偏差因子的評價(jià)方面，對目前流行的四種商業(yè)軟件[16]（PIPESIM 2012、PROPSER 13.0、PVTp 10.0、Multiflash 4.1）對天然氣偏差因子計(jì)算模塊進(jìn)行了評價(jià)，其中PROPSER 13.0 計(jì)算精度最好，對比不同的狀態(tài)方程計(jì)算的結(jié)果，SRK 方程計(jì)算結(jié)果最理想。

天然氣中一般還有酸性氣體（CO2、H2S），對于還有酸性氣體的天然氣計(jì)算的偏差因子需要做矯正。應(yīng)用酸性氣體偏差因子的計(jì)算方法[17]，根據(jù)Keys、SBV 以及RK 混合規(guī)則，酸性氣體的臨界參數(shù)校正方法，優(yōu)選出了CO2氣藏的偏差因子計(jì)算的高精度RK+WA 模型。研究富含酸性氣CO2天然氣偏差因子的計(jì)算方法[18-20]，通過實(shí)驗(yàn)結(jié)果與理論模型的對比，優(yōu)選出來基于GXO校正的DAK 模型計(jì)算誤差最小，并給出了富含CO2天然氣偏差因子隨溫度與壓力變化規(guī)律以及CO2含量對偏差因子計(jì)算結(jié)果的影響。通過對比目前天然氣中含有酸性氣體（CO2、H2S）時(shí)偏差因子的矯正方法[21]，得出DPR 模型和DAK 模型對含有酸性氣體的矯正后天然氣偏差因子計(jì)算結(jié)果最好。其他模型（GXQ、CKB）結(jié)果不理想。

根據(jù)前人的研究結(jié)果，以Poettman、Katz 和Smith對Standing-Katz 圖版低壓范圍（0.2 ≤Ppr≤15.0、1.05≤Tpr≤3.00）數(shù)值化的5 940 個(gè)數(shù)據(jù)點(diǎn)[22，23]以及Poettman 對高壓范圍（15.0≤Ppr≤30.0、1.4≤Tpr≤2.8）數(shù)值化Katz 圖版的1 208 個(gè)數(shù)據(jù)點(diǎn)為基礎(chǔ)，采用樣條函數(shù)差值的原理，對同一視對比溫度Tpr下的天然氣視對比壓力與偏差因子進(jìn)行樣條函數(shù)差值，既保證了實(shí)驗(yàn)數(shù)據(jù)點(diǎn)不偏離差值函數(shù)，又保證了差值曲線的連續(xù)光滑性。若視對比溫度Tpr不在實(shí)驗(yàn)數(shù)據(jù)給定的曲線上，則根據(jù)實(shí)驗(yàn)數(shù)據(jù)給定的視對比溫度Tpr值在Standing-Katz 圖版上選擇距離被差值的視對比溫度Tpr最近的兩條曲線，進(jìn)行樣條曲線間的二次差值，則可以得到滿足Standing-Katz 圖版的任意視對比溫度Tpr和視對比壓力Ppr下的天然氣偏差因子的值。如果天然氣中含有酸性氣體，則進(jìn)行矯正。

1 基本原理

1.1 曲線差值方法

根據(jù)Standing-Katz 圖版數(shù)值化的數(shù)據(jù)點(diǎn)，選取任意一條以視對比溫度Tpr為參數(shù)數(shù)值化離散數(shù)據(jù)點(diǎn)，將視對比壓力Ppr作為自變量，Standing-Katz 圖版上偏差因子Z 為因變量。在每個(gè)相鄰的視對比壓力之間的子區(qū)間進(jìn)行三次多項(xiàng)式差值，并且要求在整個(gè)區(qū)域子區(qū)間上的差值多項(xiàng)式函數(shù)具有二階連續(xù)導(dǎo)數(shù)，這樣得到差值函數(shù)即為三次樣條差值函數(shù)。這樣既保證了區(qū)間內(nèi)端點(diǎn)的連續(xù)性，也保證了內(nèi)端點(diǎn)處的光滑性，并且被差值的數(shù)據(jù)點(diǎn)落在差值曲線上。

假設(shè)在［Pprmin，Pprmax］區(qū)間上有n+1 個(gè)離散數(shù)據(jù)點(diǎn)Ppr（0）＜Ppr（1）＜…＜Ppr（n-1）＜Ppr（n），Pprmin為差值區(qū)間視對比壓力的最小值，Pprmax為差值區(qū)間視對比壓力的最大值。天然氣偏差因子Zi（Ppr）在區(qū)間［Ppr（i-1），Ppr（i）］上是三次多項(xiàng)式函數(shù)，這個(gè)三次多項(xiàng)式函數(shù)可表示為，即有4 個(gè)待定系數(shù)。因?yàn)閚+1 個(gè)離散數(shù)據(jù)點(diǎn)有n 個(gè)子區(qū)間，每個(gè)子區(qū)間要確定一個(gè)三次多項(xiàng)式，因此共計(jì)有4n 個(gè)待定系數(shù)需要計(jì)算。根據(jù)差值的基本條件，即共計(jì)n+1 個(gè)條件。在區(qū)間內(nèi)端點(diǎn)處，Zi（Ppr）、Zi（Ppr）的一階導(dǎo)函數(shù)以及二階導(dǎo)函數(shù)為連續(xù)函數(shù)，即可得3n-3個(gè)條件。還需要兩個(gè)條件，即邊界條件，這里取差值區(qū)間的兩個(gè)端點(diǎn)處的一階導(dǎo)數(shù)即可，其值由區(qū)間起點(diǎn)與終點(diǎn)的相鄰兩點(diǎn)的直線的斜率確定。根據(jù)上述分析，在區(qū)間x∈［Ppr（i-1），Ppr（i）］上的三次樣條差值多項(xiàng)式的數(shù)學(xué)模型可由（1）式表示：

其中：i=1，2，3，…，n。通過（2）式，可以計(jì)算位于差值區(qū)間任意視對比壓力下的天然氣偏差因子。

1.2 曲線間差值方法

2 計(jì)算結(jié)果

根據(jù)（2）式的分析，在進(jìn)行曲線差值之前，需要計(jì)算Standing-Katz 圖版數(shù)值化后的每條曲線兩個(gè)端點(diǎn)的一階導(dǎo)函數(shù)。這里分別取曲線每條數(shù)值化兩個(gè)相鄰端點(diǎn)的點(diǎn)進(jìn)行一階導(dǎo)數(shù)的計(jì)算。為了分析一階導(dǎo)數(shù)對差值結(jié)果的影響以及偏差因子的計(jì)算精度，這里分別選取兩個(gè)方案，方案一：以曲線端點(diǎn)相鄰兩點(diǎn)之間的直線的斜率代替一階導(dǎo)數(shù)。方案二：以曲線端點(diǎn)相鄰三點(diǎn)進(jìn)行一階導(dǎo)數(shù)加權(quán)計(jì)算，其中方案二的一階導(dǎo)數(shù)計(jì)算公式由（4）式表示。

其中：當(dāng)k=0 時(shí)，則j=1；當(dāng)k=n 時(shí)，則j=n-1。

根據(jù)Standing-Katz 圖版數(shù)值化的結(jié)果，結(jié)合上述一階導(dǎo)數(shù)的計(jì)算方法，曲線的兩個(gè)端點(diǎn)的一階導(dǎo)數(shù)的結(jié)果（見表1）。從表1 可以看出，方案一與方案二計(jì)算的結(jié)果差值曲線兩個(gè)端點(diǎn)的一階導(dǎo)數(shù)在左端點(diǎn)和右端點(diǎn)的相對誤差很小，基本控制在5%以內(nèi)，其中視對比溫度為1.8 和1.9 時(shí)，方案一與方案二計(jì)算的左端點(diǎn)的一階導(dǎo)數(shù)值誤差略大于5%。

根據(jù)表1 中方案一和方案二一階導(dǎo)數(shù)的計(jì)算結(jié)果，這里取二者的平均值，將二者平均值代入公式（1）求解差值內(nèi)節(jié)點(diǎn)的一階導(dǎo)數(shù)值，將內(nèi)節(jié)點(diǎn)的一階導(dǎo)數(shù)值代入到（2）式，得到對應(yīng)視對比溫度下的樣條函數(shù)；然后在區(qū)間［0.2，30］上離散視對比壓力，這里視對比壓力離散的間距為0.001，然后將離散的視對比壓力代入到（2）式，得到樣條函數(shù)差值后結(jié)果（見圖1），圖1中的點(diǎn)表示Standing-Katz 圖版的數(shù)值化的結(jié)果。這里需要說明的是Standing-Katz 圖版數(shù)值化的數(shù)據(jù)太多，在繪圖時(shí)按照一定間隔進(jìn)行了抽稀處理。通過圖1 的對比，Standing-Katz 圖版的數(shù)值化的結(jié)果與樣條差值的計(jì)算結(jié)果匹配很好，并且樣條差值法可以計(jì)算滿足Standing-Katz 圖版差值范圍的任意給定視對比溫度與視對比壓力下的天然氣偏差因子的值。

圖1 Standing-Katz 圖版與樣條曲線擬合結(jié)果對比圖

表1 方案一與方案二計(jì)算擬合曲線的左右端點(diǎn)的一階導(dǎo)數(shù)值與相對誤差

根據(jù)曲線間的差值原理，這里選取了視對比溫度分別為1.08、1.55 和2.10 的三組參數(shù)。將這三組參數(shù)與Standing-Katz 圖版的視對比溫度對比可知，這三組參數(shù)均不是Standing-Katz 圖版所給定的參數(shù)，因此這里只能通過曲線間差值的方法得到不同視對比壓力下的偏差因子值。這三組視對比溫度參數(shù)分別介于視對比溫度為［1.05，1.10］、［1.50，1.60］和［2.00，2.20］三個(gè)區(qū)間，根據(jù)（3）式的計(jì)算方法，得到這三組參數(shù)差值的結(jié)果（見圖2），Standing-Katz 圖版數(shù)值化的數(shù)據(jù)見圖2的離散點(diǎn)所示。從圖2 可以看出，差值的實(shí)線的結(jié)果與Standing-Katz 圖版數(shù)值化的變化趨勢一致。

圖2 Standing-Katz 圖版與樣條曲線擬合結(jié)果對比圖

3 實(shí)例分析

在垂直管流計(jì)算過程中，天然氣偏差系數(shù)的計(jì)算模型有很多種，根據(jù)前人的研究結(jié)果，本次研究主要采用了四種計(jì)算方法，即DAK（Dranchk-Abu-Kasse）模型，HTP（Hankinson-Thomas-Phillips）模型，李相方（LXF）模型和本文提出的樣條函數(shù)差值法。

根據(jù)文獻(xiàn)提出的氣井井底靜壓的計(jì)算原理[25]，采用微元分析法的思想，將氣井的井筒從井口到井底離散為多個(gè)單元，每一個(gè)計(jì)算微元體內(nèi)的性質(zhì)是一致的，其屬性值可以看作常數(shù)處理?；谝陨峡紤]，可以把井筒從井口到產(chǎn)層中深劃分為n 個(gè)微元體，在每個(gè)微元體內(nèi)，天然氣的溫度、偏差因子、黏度、壓力以及與壓力有關(guān)的量可以看作常數(shù)（都是基于前一個(gè)節(jié)點(diǎn)的壓力與溫度計(jì)算出的值）。

由于969 井區(qū)的天然氣中含有酸性氣體，考慮到酸性氣體對天然氣偏差因子的影響，采用了GXQ 方法和Wichert-Aziz 方法對天然氣的臨界溫度和壓力分別做了校正，具體方法參見文獻(xiàn)[19-21]。根據(jù)試氣資料的整理，統(tǒng)計(jì)了969 井區(qū)共計(jì)76 口井生產(chǎn)前的壓力恢復(fù)數(shù)據(jù)，這些井也下入了井下壓力計(jì)，測量了實(shí)際的靜壓。結(jié)合垂直管流微元分析法的計(jì)算思想，根據(jù)關(guān)井恢復(fù)后的井口油壓和套壓的測量值計(jì)算了這76 口井的靜壓。將靜壓計(jì)算的結(jié)果與實(shí)測的進(jìn)行對比，得到了969井區(qū)分別采用Dranchk-Abu-Kasse 模型[6]、Hankinson-Thomas-Phillips 模型[8]、李相方（LXF）模型[7]以及樣條函數(shù)差值模型計(jì)算的天然氣偏差因子，然后根據(jù)微元分析進(jìn)行井底靜壓計(jì)算，井底靜壓的誤差分布（見圖3）。從圖3 可以看出，本文給出的樣條函數(shù)差值法計(jì)算的偏差因子在計(jì)算井底靜壓時(shí)的相對誤差最小；其次為LXF 法，因?yàn)樵摲椒ㄔ谟?jì)算偏差因子時(shí)，采用了分段函數(shù)進(jìn)行處理，適用性較好；然后為DAK 和HTP法。其中井底靜壓相對誤差大于5%的這些井，基本為氣水同產(chǎn)，甚至有些井存在著井底積液，導(dǎo)致井底靜壓計(jì)算的相對誤差較大。

圖3 969 井區(qū)井底靜壓誤差統(tǒng)計(jì)圖

4 結(jié)論

（1）根據(jù)Standing-Katz 圖版數(shù)值化的結(jié)果，采用三次樣條函數(shù)的差值原理，給出了三次樣條差值計(jì)算天然氣偏差因子的基本原理。

（2）對于視對比溫度Tpr不在Standing-Katz 圖版上的情況，在Standing-Katz 圖版上查找距離給定視對比溫度Tpr最近的兩條曲線上的視對比溫度，然后根據(jù)視對比溫度Tpr的值進(jìn)行反距離加權(quán)二次差值。根據(jù)反距離加權(quán)二次差值，可以得到任意給定視對比溫度Tpr（1.05≤Tpr≤3.00）和視對比壓力Ppr（0.2≤Ppr≤15.0）下的天然氣偏差因子的值。

（3）通過969 井區(qū)76 口氣井的試氣資料，本文提出的三次樣條函數(shù)差值計(jì)算偏差因子方法具有較高的精度，并且計(jì)算簡單，適用性強(qiáng)。