曾冀++譚宇龍

摘 要：溫度、壓力對(duì)二氧化碳物性參數(shù)影響較大，本文采用計(jì)算精度和使用范圍較廣的Span-Wagner模型和Vesovic模型對(duì)二氧化碳物性參數(shù)進(jìn)行。結(jié)算結(jié)果表明：二氧化碳臨界點(diǎn)附近各物性參數(shù)均表現(xiàn)出較明顯的奇異性。從氣態(tài)轉(zhuǎn)變到液態(tài)（液態(tài)轉(zhuǎn)變到其他），物性參數(shù)會(huì)出現(xiàn)突變；但從氣態(tài)或液態(tài)轉(zhuǎn)變到超臨界態(tài)，物性參數(shù)呈連續(xù)變化。因此，在二氧化碳?jí)毫严嚓P(guān)模型中不能將物性參數(shù)視為定值進(jìn)行計(jì)算。

關(guān)鍵詞：二氧化碳 壓裂 物性參數(shù) 溫度 壓力

中圖分類號(hào)：TE357 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1672-3791（2014）03（c）-0088-01

隨著我國(guó)非常規(guī)儲(chǔ)層的大力開(kāi)發(fā)，傳統(tǒng)壓裂技術(shù)弊端日益凸顯，無(wú)水壓裂技術(shù)優(yōu)勢(shì)明顯，其中二氧化碳?jí)毫研Ч燃裑1～2]。與常規(guī)壓裂液流體不同，二氧化碳物性參數(shù)受溫度、壓力影響較大[3]，因此，為更好優(yōu)化二氧化碳?jí)毫咽┕?shù)，需準(zhǔn)確計(jì)算不同溫度、壓力下二氧化碳各物性參數(shù)值。

目前，使用較廣的立方型狀態(tài)方程適用范圍及計(jì)算精度都不是很盡人意[4～5]，因此，Span-Wagner基于亥姆赫茲自由能和實(shí)驗(yàn)數(shù)據(jù)修正得到了針對(duì)二氧化碳物性參數(shù)的計(jì)算模型[6]；此外，Vesovic等針對(duì)二氧化碳粘度和導(dǎo)熱系數(shù)的計(jì)算模型具有較高的精度[7]。

1 Span-Wagner模型

亥姆赫茲自由能可表示為兩個(gè)獨(dú)立的變量密度ρ和溫度T的函數(shù)，無(wú)因次亥姆赫茲自由能可由兩部分構(gòu)成：

（1）

其中：Φo為理想部分；Φr為殘余部分；δ=ρ/ρc為對(duì)比密度；τ=Tc/T為逆對(duì)比溫度；ρc為臨界密度；Tc為臨界溫度。

理想部分無(wú)因次亥姆赫茲自由能Φo：

（2）

殘余部分的無(wú)因次亥姆赫茲自由能φr表示為：

（3）

其中：

。

2 Vesovic模型

采用Vesovic模型計(jì)算二氧化碳流體粘度和導(dǎo)熱系數(shù)具有較好的精度，其通式可寫(xiě)為：

（4）

其中：Xo（T）表示僅有兩個(gè)分子相互作用的零密度極限條件下性質(zhì)；ΔcX（ρ，T）用以修正臨界點(diǎn)附近的波動(dòng)；ΔX（ρ，T）表示所以其他所有作用對(duì)性質(zhì)的貢獻(xiàn)。

3 計(jì)算結(jié)果及分析

-20 ℃、0 ℃和20 ℃三條曲線在壓力從低到高的過(guò)程中出現(xiàn)密度值的躍升，而40 ℃、60 ℃和80 ℃三條曲線所對(duì)應(yīng)的密度值則相對(duì)平緩上升，未出現(xiàn)較大的跳躍點(diǎn)。在低于臨界溫度時(shí)，對(duì)氣體二氧化碳加壓會(huì)使其從氣態(tài)轉(zhuǎn)換為液態(tài)，因此密度值發(fā)生跳躍；而在高于臨界溫度時(shí)，對(duì)氣體二氧化碳加壓，會(huì)使其從液態(tài)向超臨界態(tài)過(guò)渡，密度呈連續(xù)變化，因而無(wú)密度跳躍。

二氧化碳在不同溫度下的熱容值隨著壓力增加都是呈先增加后減小的趨勢(shì)，表現(xiàn)出明顯的峰值性，且溫度越高的曲線峰值所對(duì)應(yīng)的壓力值也越高。在低壓段，已經(jīng)進(jìn)入液態(tài)的-20 ℃、0 ℃和20 ℃三條曲線較仍處于氣態(tài)的三條曲線所對(duì)應(yīng)熱容明顯更大；但高壓段之后，六條曲線幾乎重合；而介于二者之間的區(qū)域，由于存在臨界點(diǎn)的熱容奇異點(diǎn)，此時(shí)的熱容值為溫度、壓力協(xié)同控制。

二氧化碳在不同溫度下粘度隨壓力增大而逐漸增大。-20 ℃、0 ℃和20 ℃三條曲線從氣態(tài)變?yōu)橐簯B(tài)時(shí)，二氧化碳粘度表現(xiàn)出明顯的陡升，而剩下3條曲線從氣態(tài)變?yōu)槌R界態(tài)時(shí)，其粘度升高趨勢(shì)相對(duì)平緩，且大小介于氣態(tài)和液態(tài)之間。整體看來(lái)，二氧化碳粘度主要受溫度、壓力協(xié)同控制作用。

二氧化碳導(dǎo)熱系數(shù)在不同溫度下隨壓力變化趨勢(shì)與粘度變化趨勢(shì)類似，皆是受溫度、壓力協(xié)同控制，但溫度對(duì)導(dǎo)熱系數(shù)的影響程度較其對(duì)粘度的影響程度更弱。

4 結(jié)語(yǔ)

使用S-W模型和Vesovic模型能較精確的計(jì)算二氧化碳在不同溫度、壓力下各物性參數(shù)。在整個(gè)壓裂施工過(guò)程中，熱容表現(xiàn)為壓力的函數(shù)，焓值表現(xiàn)為溫度的函數(shù)，而其他各物性參數(shù)表現(xiàn)為溫度和壓力的函數(shù)。因此，在二氧化碳?jí)毫严嚓P(guān)模型中不能將其物性參數(shù)考慮為定值計(jì)算。

參考文獻(xiàn)

[1] 唐穎，唐玄，王廣源，等.頁(yè)巖氣開(kāi)發(fā)水力壓裂技術(shù)綜述[J].地質(zhì)通報(bào)，2011，30（2）：393-399.

[2] 侯向前，盧擁軍，方波，等.非常規(guī)儲(chǔ)集層低碳烴無(wú)水壓裂液[J].石油勘探與開(kāi)發(fā)，2013，40（5）：601-605.

[3] 王海柱，沈忠厚，李根生.超臨界CO2鉆井井筒壓力溫度耦合計(jì)算[J].石油勘探與開(kāi)發(fā)，2011，38（1）：97-102.

[4] 王海柱，沈忠厚，李根生，等.CO2氣體物性參數(shù)精確計(jì)算方法研究[J].石油鉆采工藝，2011，33（5）：65-67.

[5] 里德，芝芬，怡生，等.氣體和液體性質(zhì)[M].石油工業(yè)出版社，1994.

[6] Span R，Wagner W.A new equation of state for carbon dioxide covering the fluid region from the triple‐point temperature to 1100 K at pressures up to 800 MPa[J].Journal of physical and chemical reference data，1996，25（6）：1509-1596.

[7] Vesovic V，Wakeham W A，Olchowy G A，et al.The transport properties of carbon dioxide[J].Journal of Physical and Chemical Reference Data，1990，19（3）：763-808.endprint

摘 要：溫度、壓力對(duì)二氧化碳物性參數(shù)影響較大，本文采用計(jì)算精度和使用范圍較廣的Span-Wagner模型和Vesovic模型對(duì)二氧化碳物性參數(shù)進(jìn)行。結(jié)算結(jié)果表明：二氧化碳臨界點(diǎn)附近各物性參數(shù)均表現(xiàn)出較明顯的奇異性。從氣態(tài)轉(zhuǎn)變到液態(tài)（液態(tài)轉(zhuǎn)變到其他），物性參數(shù)會(huì)出現(xiàn)突變；但從氣態(tài)或液態(tài)轉(zhuǎn)變到超臨界態(tài)，物性參數(shù)呈連續(xù)變化。因此，在二氧化碳?jí)毫严嚓P(guān)模型中不能將物性參數(shù)視為定值進(jìn)行計(jì)算。

關(guān)鍵詞：二氧化碳 壓裂 物性參數(shù) 溫度 壓力

中圖分類號(hào)：TE357 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1672-3791（2014）03（c）-0088-01

隨著我國(guó)非常規(guī)儲(chǔ)層的大力開(kāi)發(fā)，傳統(tǒng)壓裂技術(shù)弊端日益凸顯，無(wú)水壓裂技術(shù)優(yōu)勢(shì)明顯，其中二氧化碳?jí)毫研Ч燃裑1～2]。與常規(guī)壓裂液流體不同，二氧化碳物性參數(shù)受溫度、壓力影響較大[3]，因此，為更好優(yōu)化二氧化碳?jí)毫咽┕?shù)，需準(zhǔn)確計(jì)算不同溫度、壓力下二氧化碳各物性參數(shù)值。

目前，使用較廣的立方型狀態(tài)方程適用范圍及計(jì)算精度都不是很盡人意[4～5]，因此，Span-Wagner基于亥姆赫茲自由能和實(shí)驗(yàn)數(shù)據(jù)修正得到了針對(duì)二氧化碳物性參數(shù)的計(jì)算模型[6]；此外，Vesovic等針對(duì)二氧化碳粘度和導(dǎo)熱系數(shù)的計(jì)算模型具有較高的精度[7]。

1 Span-Wagner模型

亥姆赫茲自由能可表示為兩個(gè)獨(dú)立的變量密度ρ和溫度T的函數(shù)，無(wú)因次亥姆赫茲自由能可由兩部分構(gòu)成：

（1）

其中：Φo為理想部分；Φr為殘余部分；δ=ρ/ρc為對(duì)比密度；τ=Tc/T為逆對(duì)比溫度；ρc為臨界密度；Tc為臨界溫度。

理想部分無(wú)因次亥姆赫茲自由能Φo：

（2）

殘余部分的無(wú)因次亥姆赫茲自由能φr表示為：

（3）

其中：

。

2 Vesovic模型

采用Vesovic模型計(jì)算二氧化碳流體粘度和導(dǎo)熱系數(shù)具有較好的精度，其通式可寫(xiě)為：

（4）

其中：Xo（T）表示僅有兩個(gè)分子相互作用的零密度極限條件下性質(zhì)；ΔcX（ρ，T）用以修正臨界點(diǎn)附近的波動(dòng)；ΔX（ρ，T）表示所以其他所有作用對(duì)性質(zhì)的貢獻(xiàn)。

3 計(jì)算結(jié)果及分析

-20 ℃、0 ℃和20 ℃三條曲線在壓力從低到高的過(guò)程中出現(xiàn)密度值的躍升，而40 ℃、60 ℃和80 ℃三條曲線所對(duì)應(yīng)的密度值則相對(duì)平緩上升，未出現(xiàn)較大的跳躍點(diǎn)。在低于臨界溫度時(shí)，對(duì)氣體二氧化碳加壓會(huì)使其從氣態(tài)轉(zhuǎn)換為液態(tài)，因此密度值發(fā)生跳躍；而在高于臨界溫度時(shí)，對(duì)氣體二氧化碳加壓，會(huì)使其從液態(tài)向超臨界態(tài)過(guò)渡，密度呈連續(xù)變化，因而無(wú)密度跳躍。

二氧化碳在不同溫度下的熱容值隨著壓力增加都是呈先增加后減小的趨勢(shì)，表現(xiàn)出明顯的峰值性，且溫度越高的曲線峰值所對(duì)應(yīng)的壓力值也越高。在低壓段，已經(jīng)進(jìn)入液態(tài)的-20 ℃、0 ℃和20 ℃三條曲線較仍處于氣態(tài)的三條曲線所對(duì)應(yīng)熱容明顯更大；但高壓段之后，六條曲線幾乎重合；而介于二者之間的區(qū)域，由于存在臨界點(diǎn)的熱容奇異點(diǎn)，此時(shí)的熱容值為溫度、壓力協(xié)同控制。

二氧化碳在不同溫度下粘度隨壓力增大而逐漸增大。-20 ℃、0 ℃和20 ℃三條曲線從氣態(tài)變?yōu)橐簯B(tài)時(shí)，二氧化碳粘度表現(xiàn)出明顯的陡升，而剩下3條曲線從氣態(tài)變?yōu)槌R界態(tài)時(shí)，其粘度升高趨勢(shì)相對(duì)平緩，且大小介于氣態(tài)和液態(tài)之間。整體看來(lái)，二氧化碳粘度主要受溫度、壓力協(xié)同控制作用。

二氧化碳導(dǎo)熱系數(shù)在不同溫度下隨壓力變化趨勢(shì)與粘度變化趨勢(shì)類似，皆是受溫度、壓力協(xié)同控制，但溫度對(duì)導(dǎo)熱系數(shù)的影響程度較其對(duì)粘度的影響程度更弱。

4 結(jié)語(yǔ)

使用S-W模型和Vesovic模型能較精確的計(jì)算二氧化碳在不同溫度、壓力下各物性參數(shù)。在整個(gè)壓裂施工過(guò)程中，熱容表現(xiàn)為壓力的函數(shù)，焓值表現(xiàn)為溫度的函數(shù)，而其他各物性參數(shù)表現(xiàn)為溫度和壓力的函數(shù)。因此，在二氧化碳?jí)毫严嚓P(guān)模型中不能將其物性參數(shù)考慮為定值計(jì)算。

參考文獻(xiàn)

[1] 唐穎，唐玄，王廣源，等.頁(yè)巖氣開(kāi)發(fā)水力壓裂技術(shù)綜述[J].地質(zhì)通報(bào)，2011，30（2）：393-399.

[2] 侯向前，盧擁軍，方波，等.非常規(guī)儲(chǔ)集層低碳烴無(wú)水壓裂液[J].石油勘探與開(kāi)發(fā)，2013，40（5）：601-605.

[3] 王海柱，沈忠厚，李根生.超臨界CO2鉆井井筒壓力溫度耦合計(jì)算[J].石油勘探與開(kāi)發(fā)，2011，38（1）：97-102.

[4] 王海柱，沈忠厚，李根生，等.CO2氣體物性參數(shù)精確計(jì)算方法研究[J].石油鉆采工藝，2011，33（5）：65-67.

[5] 里德，芝芬，怡生，等.氣體和液體性質(zhì)[M].石油工業(yè)出版社，1994.

[6] Span R，Wagner W.A new equation of state for carbon dioxide covering the fluid region from the triple‐point temperature to 1100 K at pressures up to 800 MPa[J].Journal of physical and chemical reference data，1996，25（6）：1509-1596.

[7] Vesovic V，Wakeham W A，Olchowy G A，et al.The transport properties of carbon dioxide[J].Journal of Physical and Chemical Reference Data，1990，19（3）：763-808.endprint

摘 要：溫度、壓力對(duì)二氧化碳物性參數(shù)影響較大，本文采用計(jì)算精度和使用范圍較廣的Span-Wagner模型和Vesovic模型對(duì)二氧化碳物性參數(shù)進(jìn)行。結(jié)算結(jié)果表明：二氧化碳臨界點(diǎn)附近各物性參數(shù)均表現(xiàn)出較明顯的奇異性。從氣態(tài)轉(zhuǎn)變到液態(tài)（液態(tài)轉(zhuǎn)變到其他），物性參數(shù)會(huì)出現(xiàn)突變；但從氣態(tài)或液態(tài)轉(zhuǎn)變到超臨界態(tài)，物性參數(shù)呈連續(xù)變化。因此，在二氧化碳?jí)毫严嚓P(guān)模型中不能將物性參數(shù)視為定值進(jìn)行計(jì)算。

關(guān)鍵詞：二氧化碳 壓裂 物性參數(shù) 溫度 壓力

中圖分類號(hào)：TE357 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1672-3791（2014）03（c）-0088-01

隨著我國(guó)非常規(guī)儲(chǔ)層的大力開(kāi)發(fā)，傳統(tǒng)壓裂技術(shù)弊端日益凸顯，無(wú)水壓裂技術(shù)優(yōu)勢(shì)明顯，其中二氧化碳?jí)毫研Ч燃裑1～2]。與常規(guī)壓裂液流體不同，二氧化碳物性參數(shù)受溫度、壓力影響較大[3]，因此，為更好優(yōu)化二氧化碳?jí)毫咽┕?shù)，需準(zhǔn)確計(jì)算不同溫度、壓力下二氧化碳各物性參數(shù)值。

目前，使用較廣的立方型狀態(tài)方程適用范圍及計(jì)算精度都不是很盡人意[4～5]，因此，Span-Wagner基于亥姆赫茲自由能和實(shí)驗(yàn)數(shù)據(jù)修正得到了針對(duì)二氧化碳物性參數(shù)的計(jì)算模型[6]；此外，Vesovic等針對(duì)二氧化碳粘度和導(dǎo)熱系數(shù)的計(jì)算模型具有較高的精度[7]。

1 Span-Wagner模型

亥姆赫茲自由能可表示為兩個(gè)獨(dú)立的變量密度ρ和溫度T的函數(shù)，無(wú)因次亥姆赫茲自由能可由兩部分構(gòu)成：

（1）

其中：Φo為理想部分；Φr為殘余部分；δ=ρ/ρc為對(duì)比密度；τ=Tc/T為逆對(duì)比溫度；ρc為臨界密度；Tc為臨界溫度。

理想部分無(wú)因次亥姆赫茲自由能Φo：

（2）

殘余部分的無(wú)因次亥姆赫茲自由能φr表示為：

（3）

其中：

。

2 Vesovic模型

采用Vesovic模型計(jì)算二氧化碳流體粘度和導(dǎo)熱系數(shù)具有較好的精度，其通式可寫(xiě)為：

（4）

其中：Xo（T）表示僅有兩個(gè)分子相互作用的零密度極限條件下性質(zhì)；ΔcX（ρ，T）用以修正臨界點(diǎn)附近的波動(dòng)；ΔX（ρ，T）表示所以其他所有作用對(duì)性質(zhì)的貢獻(xiàn)。

3 計(jì)算結(jié)果及分析

-20 ℃、0 ℃和20 ℃三條曲線在壓力從低到高的過(guò)程中出現(xiàn)密度值的躍升，而40 ℃、60 ℃和80 ℃三條曲線所對(duì)應(yīng)的密度值則相對(duì)平緩上升，未出現(xiàn)較大的跳躍點(diǎn)。在低于臨界溫度時(shí)，對(duì)氣體二氧化碳加壓會(huì)使其從氣態(tài)轉(zhuǎn)換為液態(tài)，因此密度值發(fā)生跳躍；而在高于臨界溫度時(shí)，對(duì)氣體二氧化碳加壓，會(huì)使其從液態(tài)向超臨界態(tài)過(guò)渡，密度呈連續(xù)變化，因而無(wú)密度跳躍。

二氧化碳在不同溫度下的熱容值隨著壓力增加都是呈先增加后減小的趨勢(shì)，表現(xiàn)出明顯的峰值性，且溫度越高的曲線峰值所對(duì)應(yīng)的壓力值也越高。在低壓段，已經(jīng)進(jìn)入液態(tài)的-20 ℃、0 ℃和20 ℃三條曲線較仍處于氣態(tài)的三條曲線所對(duì)應(yīng)熱容明顯更大；但高壓段之后，六條曲線幾乎重合；而介于二者之間的區(qū)域，由于存在臨界點(diǎn)的熱容奇異點(diǎn)，此時(shí)的熱容值為溫度、壓力協(xié)同控制。

二氧化碳在不同溫度下粘度隨壓力增大而逐漸增大。-20 ℃、0 ℃和20 ℃三條曲線從氣態(tài)變?yōu)橐簯B(tài)時(shí)，二氧化碳粘度表現(xiàn)出明顯的陡升，而剩下3條曲線從氣態(tài)變?yōu)槌R界態(tài)時(shí)，其粘度升高趨勢(shì)相對(duì)平緩，且大小介于氣態(tài)和液態(tài)之間。整體看來(lái)，二氧化碳粘度主要受溫度、壓力協(xié)同控制作用。

二氧化碳導(dǎo)熱系數(shù)在不同溫度下隨壓力變化趨勢(shì)與粘度變化趨勢(shì)類似，皆是受溫度、壓力協(xié)同控制，但溫度對(duì)導(dǎo)熱系數(shù)的影響程度較其對(duì)粘度的影響程度更弱。

4 結(jié)語(yǔ)

使用S-W模型和Vesovic模型能較精確的計(jì)算二氧化碳在不同溫度、壓力下各物性參數(shù)。在整個(gè)壓裂施工過(guò)程中，熱容表現(xiàn)為壓力的函數(shù)，焓值表現(xiàn)為溫度的函數(shù)，而其他各物性參數(shù)表現(xiàn)為溫度和壓力的函數(shù)。因此，在二氧化碳?jí)毫严嚓P(guān)模型中不能將其物性參數(shù)考慮為定值計(jì)算。

參考文獻(xiàn)

[1] 唐穎，唐玄，王廣源，等.頁(yè)巖氣開(kāi)發(fā)水力壓裂技術(shù)綜述[J].地質(zhì)通報(bào)，2011，30（2）：393-399.

[2] 侯向前，盧擁軍，方波，等.非常規(guī)儲(chǔ)集層低碳烴無(wú)水壓裂液[J].石油勘探與開(kāi)發(fā)，2013，40（5）：601-605.

[3] 王海柱，沈忠厚，李根生.超臨界CO2鉆井井筒壓力溫度耦合計(jì)算[J].石油勘探與開(kāi)發(fā)，2011，38（1）：97-102.

[4] 王海柱，沈忠厚，李根生，等.CO2氣體物性參數(shù)精確計(jì)算方法研究[J].石油鉆采工藝，2011，33（5）：65-67.

[5] 里德，芝芬，怡生，等.氣體和液體性質(zhì)[M].石油工業(yè)出版社，1994.

[6] Span R，Wagner W.A new equation of state for carbon dioxide covering the fluid region from the triple‐point temperature to 1100 K at pressures up to 800 MPa[J].Journal of physical and chemical reference data，1996，25（6）：1509-1596.

[7] Vesovic V，Wakeham W A，Olchowy G A，et al.The transport properties of carbon dioxide[J].Journal of Physical and Chemical Reference Data，1990，19（3）：763-808.endprint