喬加飛，李 卓，廖海燕，王兵兵，劉航滔，段遠(yuǎn)源

（1.國家能源集團(tuán)新能源技術(shù)研究院有限公司，北京 100024；2.清華大學(xué)熱科學(xué)與動(dòng)力工程教育部重點(diǎn)實(shí)驗(yàn)室，北京 100084）

2020年9月，國家主席習(xí)近平提出2030年前實(shí)現(xiàn)碳達(dá)峰、2060年前實(shí)現(xiàn)碳中和的宏偉目標(biāo)?！半p碳”目標(biāo)下，對(duì)于清潔、低碳、安全、高效的能源利用方式的需求愈發(fā)迫切[1]。超臨界二氧化碳（S-CO2）循環(huán)具有成本低、效率高、設(shè)備尺寸小、靈活性強(qiáng)等優(yōu)點(diǎn)[2]，在包括核能[3]、太陽能[4]、燃煤[5]、沼氣發(fā)電[6]等方面具有廣闊的應(yīng)用空間，主要分為閉式（間接加熱式）和半閉式（直接加熱式）循環(huán)。閉式S-CO2循環(huán)多采用純CO2為工質(zhì)，通過外部熱源與換熱器對(duì)工質(zhì)加熱；半閉式S-CO2循環(huán)使燃料在燃燒室直接燃燒，產(chǎn)生富CO2混合體系作為工質(zhì)循環(huán)做功，同時(shí)結(jié)合碳捕集、利用與封存（carbon capture, utilization and storage，CCUS）技術(shù)可以捕集高純度的CO2。

S-CO2循環(huán)的重要特點(diǎn)之一是循環(huán)工質(zhì)跨臨界或始終保持在超臨界狀態(tài)。與其他工質(zhì)相比，CO2及其混合體系具有較低的臨界參數(shù)，易達(dá)到超臨界狀態(tài)，具備黏度小、密度大、可壓縮性小等一系列優(yōu)點(diǎn)，因而適用于多場(chǎng)景的動(dòng)力循環(huán)。對(duì)于閉式循環(huán)系統(tǒng)，Liu等人[7]、Guo等人[8]提出引入CO2混合工質(zhì)，調(diào)整工質(zhì)臨界參數(shù)，實(shí)現(xiàn)不同組元的優(yōu)勢(shì)互補(bǔ)。此外，結(jié)合富氧燃燒技術(shù)的半閉式S-CO2循環(huán)，由于其高效和低碳的特性，也越來越受到研究者的重視，目前已提出多種循環(huán)方式，如半閉式富氧燃燒聯(lián)合循環(huán)（semi-closed oxy-fuel combustion combined cycle，SCOC-CC）循環(huán)[9-11]、E-MATIANT循環(huán)[12]、Allam循環(huán)[13]等。對(duì)于這些循環(huán)的深入研究和工程應(yīng)用，均需要深刻認(rèn)識(shí)CO2及其混合體系的熱力學(xué)性質(zhì)。

結(jié)合閉式S-CO2動(dòng)力系統(tǒng)的設(shè)計(jì)需求，楊富方等[14]綜述了CO2在循環(huán)范圍內(nèi)的實(shí)驗(yàn)與模型研究現(xiàn)狀，認(rèn)為已有多參數(shù)狀態(tài)方程[15]可在一定的溫度壓力范圍內(nèi)準(zhǔn)確描述CO2純工質(zhì)的熱力學(xué)性質(zhì)，但在近臨界區(qū)的定性和定量表現(xiàn)不夠理想，可采用跨接方法[16]改進(jìn)其表現(xiàn)。何欣欣等[17]以S-CO2循環(huán)為應(yīng)用場(chǎng)景介紹了CO2熱物性方程研究進(jìn)展及應(yīng)用，推薦使用GERG-2008模型[18]計(jì)算CO2熱力學(xué)性質(zhì)。針對(duì)CO2混合工質(zhì)，Li等人[19-20]圍繞CCUS技術(shù)介紹了CO2混合體系實(shí)驗(yàn)及理論研究情況，而在半閉式S-CO2循環(huán)設(shè)計(jì)關(guān)注的組分、溫度、壓力范圍內(nèi)，缺乏對(duì)混合工質(zhì)熱力學(xué)性質(zhì)的評(píng)估。混合體系的熱物性并非純工質(zhì)熱物性的線性加權(quán)，當(dāng)組元性質(zhì)差異較大時(shí)，非理想性更突出，其理論描述也需要結(jié)合特定混合規(guī)則，構(gòu)建合適的熱力學(xué)模型。

對(duì)此，本文分析了CO2混合工質(zhì)熱力學(xué)性質(zhì)實(shí)驗(yàn)數(shù)據(jù)和理論模型的研究現(xiàn)狀，并根據(jù)半閉式S-CO2循環(huán)研究需求，指出現(xiàn)有實(shí)驗(yàn)數(shù)據(jù)和模型研究的不足及發(fā)展方向。

1 典型系統(tǒng)和工況

不同于閉式循環(huán)，半閉式S-CO2循環(huán)工質(zhì)在燃燒室內(nèi)被直接加熱，通常透平進(jìn)口溫度更高，以實(shí)現(xiàn)更高的熱效率。圖1為簡(jiǎn)單回?zé)岚腴]式S-CO2循環(huán)系統(tǒng)示意。循環(huán)中，壓縮機(jī)出口處的工質(zhì)溫度和壓力略高于CO2純工質(zhì)的臨界點(diǎn)，以防止混合體系發(fā)生液化，影響系統(tǒng)的安全穩(wěn)定運(yùn)行[21]。壓縮后，工質(zhì)與由空氣分離器制備的O2混合后進(jìn)入回?zé)崞魑鼰?，然后在燃燒室中與燃料反應(yīng)吸收熱量；在部分循環(huán)設(shè)計(jì)中，工質(zhì)與O2也會(huì)直接在燃燒室中混合。高溫高壓的混合工質(zhì)膨脹做功后，在回?zé)崞鲀?nèi)放熱，隨后在冷凝器，混合工質(zhì)中的水被分離，從而獲得高純度的CO2，其中，燃燒產(chǎn)生的CO2被儲(chǔ)存，其余作為工質(zhì)進(jìn)入壓縮機(jī)繼續(xù)循環(huán)，保證循環(huán)過程中CO2工質(zhì)量恒定。

圖1 半閉式S-CO2循環(huán)示意Fig.1 Schematic diagram of semi-closed S-CO2 cycle

與閉式S-CO2循環(huán)系統(tǒng)相同，半閉式S-CO2循環(huán)系統(tǒng)的工質(zhì)會(huì)經(jīng)過CO2的近臨界區(qū)和高溫高壓區(qū)，臨界參數(shù)、氣液相平衡、密度、比熱容、聲速等熱力學(xué)性質(zhì)的變化對(duì)熱力過程與系統(tǒng)設(shè)計(jì)、優(yōu)化及運(yùn)行穩(wěn)定產(chǎn)生較大影響。此外，由于涉及混合體系，明確循環(huán)所關(guān)注的工質(zhì)組成和CO2含量也非常重要。表1列出了部分學(xué)者研究的半閉式S-CO2循環(huán)關(guān)鍵溫度壓力及組分參數(shù)，表2列出了對(duì)于循環(huán)關(guān)鍵節(jié)點(diǎn)計(jì)算的工質(zhì)組分，混合工質(zhì)主要組分為CO2、O2、H2O，還混有摩爾分?jǐn)?shù)均不超過2%的N2和Ar。其中，O2組分峰值從回?zé)徇^程到進(jìn)入燃燒室前呈現(xiàn)低溫高壓（壓縮后的CO2與O2混合）及高溫高壓（回?zé)岷螅┑奶攸c(diǎn)。H2O是燃燒產(chǎn)物，其組分峰值出現(xiàn)于燃燒室后，呈現(xiàn)高溫高壓（燃燒后）及高溫低壓（膨脹做功及回?zé)岷螅┑奶攸c(diǎn)。由于N2與Ar化學(xué)性質(zhì)穩(wěn)定，在循環(huán)過程中將作為雜質(zhì)始終存在。在壓縮機(jī)進(jìn)、出口處，混合體系與純工質(zhì)近似，僅含有極少雜質(zhì)，因此不在本文重點(diǎn)研究范圍內(nèi)。

表1 半閉式S-CO2循環(huán)最高參數(shù)Tab.1 The highest parameters of semi-closed S-CO2 cycle

表2 半閉式S-CO2循環(huán)關(guān)鍵組分Tab.2 Key components of semi-closed S-CO2 cycle

綜上，CO2摩爾分?jǐn)?shù)超80%的混合工質(zhì)在半閉式S-CO2循環(huán)中經(jīng)歷了從CO2臨界點(diǎn)到最高溫度1 400 ℃、壓力40 MPa的狀態(tài)變化，因此需要在廣闊的溫度壓力范圍內(nèi)掌握工質(zhì)準(zhǔn)確可靠的熱力學(xué)性質(zhì)。由于高溫高壓的工質(zhì)性質(zhì)偏離理想氣體規(guī)律較遠(yuǎn)，因此該多元混合體系不能視為理想氣體混合物，需要慎重考慮不同組元之間的相互作用。針對(duì)S-CO2循環(huán)工質(zhì)的特殊性，要建立可靠的理論模型，需要CO2與其他組元，特別是H2O和O2的二元交互系數(shù)，并使用包含多組元的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行驗(yàn)證。因此，本文將重點(diǎn)關(guān)注含CO2的二元與多元體系實(shí)驗(yàn)與理論研究現(xiàn)狀，并給出研究建議。

2 實(shí)驗(yàn)研究現(xiàn)狀

2.1 臨界性質(zhì)與飽和性質(zhì)

純凈物的臨界點(diǎn)為氣液共存的最高溫度和壓力狀態(tài)。混合體系不僅增加了組分變量，使系統(tǒng)自由度上升，其真實(shí)臨界軌跡線也更復(fù)雜。對(duì)于混合物的熱力學(xué)性質(zhì)，采用混合規(guī)則計(jì)算得到的臨界參數(shù)，例如最簡(jiǎn)單的是將不同組元臨界溫度以摩爾分?jǐn)?shù)為權(quán)重線性加和的Kay規(guī)則[31]，得到的并非真實(shí)的臨界軌跡線。

圖2展示了CO2混合體系臨界性質(zhì)的實(shí)驗(yàn)數(shù)據(jù)分布，根據(jù)Abdulagatov等人[32]的研究，在p-T圖上，CO2+O2[33-36]與CO2+N2[35-38]體系的臨界特性表現(xiàn)為上凸曲線，CO2+Ar[32]體系在低Ar摩爾分?jǐn)?shù)時(shí)表現(xiàn)為趨近線性，而CO2+H2O[39-42]體系的臨界性質(zhì)表現(xiàn)則較為復(fù)雜，其臨界曲線不連續(xù)且向超高壓方向延伸。對(duì)于S-CO2循環(huán)，N2與Ar的臨界溫度與壓力均低于CO2，但雜質(zhì)的引入會(huì)增加混合體系的臨界壓力，在系統(tǒng)循環(huán)設(shè)計(jì)時(shí)需格外注意。目前對(duì)CO2+H2O體系的臨界性質(zhì)研究還很不充分。

圖2 CO2混合體系臨界性質(zhì)實(shí)驗(yàn)數(shù)據(jù)分布Fig.2 Distribution of experimental data on critical properties of CO2 mixtures

氣液相平衡（vapour-liquid equilibrium，VLE）性質(zhì)描述的是混合工質(zhì)氣液兩相達(dá)到平衡時(shí)的溫度T、壓力p和氣液兩相組分yi、xi之間的關(guān)系，是混合體系最基本的平衡性質(zhì)之一?；旌象w系的壓力并非簡(jiǎn)單地隨組分摩爾分?jǐn)?shù)呈線性變化，多元體系的相平衡情況更加復(fù)雜[43]。

圖3展示了氣液相平衡實(shí)驗(yàn)數(shù)據(jù)在p-T圖上的分布。由于其余組元臨界溫度壓力遠(yuǎn)低于CO2，故在循環(huán)范圍內(nèi)只有CO2+H2O的氣液相平衡數(shù)據(jù)。

圖3 CO2混合體系氣液相平衡實(shí)驗(yàn)數(shù)據(jù)分布Fig.3 VLE experimental data distribution for CO2 mixtures

圖4展示了CO2氣相摩爾分?jǐn)?shù)大于75%的CO2+H2O混合體系氣液相平衡數(shù)據(jù)在p-T圖上的分布情況。目前尚無針對(duì)S-CO2循環(huán)工質(zhì)開展的近臨界區(qū)域相平衡研究，但基于CO2的超臨界萃取[44]、超臨界水氧化技術(shù)[45]、碳捕集封存（carbon capture storage，CCS）及石油工業(yè)[46-49]等背景的CO2+H2O體系相平衡研究實(shí)驗(yàn)數(shù)據(jù)較為豐富。未來應(yīng)開展更多CO2+H2O體系的相平衡實(shí)驗(yàn)，并獲得其他體系的相平衡數(shù)據(jù)，以提升混合體系狀態(tài)方程的精度。

圖4 CO2+H2O體系氣液相平衡實(shí)驗(yàn)數(shù)據(jù)分布(x(CO2)>75%)Fig.4 VLE experimental data distribution for CO2+H2O mixtures (x(CO2)>75%))

2.2 單相區(qū)密度、比熱容與聲速

混合物的單相區(qū)密度，又稱為pvTx性質(zhì)，是最基礎(chǔ)的熱力學(xué)性質(zhì)之一，對(duì)工質(zhì)的壓縮、流動(dòng)和膨脹過程均有重要影響。圖5展示了含CO2體系的pvTx實(shí)驗(yàn)數(shù)據(jù)分布?？梢钥吹?，含CO2的三元及以上體系的pvTx實(shí)驗(yàn)數(shù)據(jù)非常稀少。對(duì)于CO2+N2+O2三元體系，僅有Muirbrook等人[35]于1962年測(cè)定的34個(gè)數(shù)據(jù)點(diǎn)，溫度壓力范圍覆蓋273～298 K和0.1～10.6 Mpa，尚未達(dá)到S-CO2循環(huán)所關(guān)注的溫度范圍。對(duì)于CO2+N2+Ar三元體系，Yang等人[50]采用目前密度測(cè)量精度最高的磁懸浮密度計(jì)法獲得了323.15～423.15 K、3～31 MPa的富CO2體系的150個(gè)數(shù)據(jù)點(diǎn)，測(cè)量不確定度在0.2%以內(nèi)。而對(duì)于CO2+N2+O2+Ar四元體系，僅有Nazeri等人[51]使用振動(dòng)管密度計(jì)測(cè)定的43個(gè)數(shù)據(jù)點(diǎn)，溫度壓力范圍覆蓋273～423 K和1.5～10.6 MPa，x(CO2)=0.898 3，氣相測(cè)量數(shù)據(jù)不確定度為1%。對(duì)于溫度高于423 K的區(qū)域，僅有CO2+N2與CO2+H2O體系有少量實(shí)驗(yàn)數(shù)據(jù)，而其余二元或三元體系均缺乏實(shí)驗(yàn)研究。大量實(shí)驗(yàn)數(shù)據(jù)集中于低CO2摩爾分?jǐn)?shù)的區(qū)域，包括CO2在水中的溶解度實(shí)驗(yàn)[52-56]，以及CCS技術(shù)中低CO2摩爾分?jǐn)?shù)體系[57]的研究。而在S-CO2循環(huán)范圍的實(shí)驗(yàn)數(shù)據(jù)，特別是重點(diǎn)關(guān)注的CO2+O2以及CO2+H2O實(shí)驗(yàn)數(shù)據(jù)較為缺乏。

圖5 CO2混合體系pvTx性質(zhì)實(shí)驗(yàn)數(shù)據(jù)分布Fig.5 Experimental data distribution for pvTx property of CO2 mixtures

Kunz等人[58]對(duì)與天然氣相關(guān)的混合體系進(jìn)行評(píng)估時(shí)，指出使用磁懸浮密度計(jì)法測(cè)量的數(shù)據(jù)代表了pvTx實(shí)驗(yàn)數(shù)據(jù)的不確定度下限水平，約為0.1%。在S-CO2循環(huán)涉及的體系中，Yang等人[50,57]、Souissi等人[59]、Mondejar等人[60]采用此方法進(jìn)行了測(cè)量，其可靠性應(yīng)受到重視。由于實(shí)驗(yàn)管路和實(shí)驗(yàn)設(shè)備的限制，磁懸浮密度計(jì)法目前較難測(cè)量CO2+H2O的相關(guān)體系，導(dǎo)致該體系的pvTx(y)實(shí)驗(yàn)精度受到限制。近年來，純氧燃燒結(jié)合S-CO2循環(huán)技術(shù)越來越受到熱物性實(shí)驗(yàn)研究者的關(guān)注，數(shù)據(jù)空白最嚴(yán)重的CO2+O2體系也陸續(xù)有新的實(shí)驗(yàn)開展[61-62]，Park等人[63]還專門開展了面向S-CO2循環(huán)需求的物性實(shí)驗(yàn)研究。

比熱容和聲速實(shí)驗(yàn)數(shù)據(jù)會(huì)影響理論模型在S-CO2循環(huán)重點(diǎn)溫度壓力范圍的精度，也會(huì)作為參數(shù)參與循環(huán)的設(shè)計(jì)與優(yōu)化。對(duì)于CO2純工質(zhì)，比定壓熱容cp和比定容熱容cv在臨界點(diǎn)趨于無窮大；聲速w在臨界點(diǎn)則趨于0[14]。由于實(shí)驗(yàn)數(shù)據(jù)存在大量空白，目前對(duì)混合體系的比熱容與聲速缺乏清晰的認(rèn)知。目前尚無比定容熱容的實(shí)驗(yàn)數(shù)據(jù)，僅有Bishnoi等人[64]測(cè)定的CO2+N2體系在313～363 K、0.2～16.5 MPa下的比定壓熱容數(shù)據(jù)。此外，Barbero等人[65]和Hnědkovsky等人[66]分別在常溫和高溫高壓下測(cè)量了CO2水溶液的表觀比熱容，以研究非電解質(zhì)稀溶液的比熱容表現(xiàn)，但由于采用了不同的表觀體積與表觀比熱容定義，對(duì)S-CO2熱力學(xué)性質(zhì)的實(shí)驗(yàn)數(shù)據(jù)庫不具有參考價(jià)值。對(duì)于聲速，CO2與O2[67]、N2[67-69]、Ar[67,70]的二元混合體系均有實(shí)驗(yàn)數(shù)據(jù)，其中Wegge等人[70]使用球共鳴法測(cè)量得到了CO2+Ar體系的聲速數(shù)據(jù)，此方法也被Kunz等人[58]認(rèn)為是最精確的聲速測(cè)量方法之一，不確定度小于0.1%。在溫度高于500 K的區(qū)域，缺乏關(guān)于CO2混合體系的比熱容與聲速的實(shí)驗(yàn)研究。CO2混合體系聲速與比定壓熱容實(shí)驗(yàn)數(shù)據(jù)分布情況見圖6。

圖6 CO2混合體系聲速與比定壓熱容實(shí)驗(yàn)數(shù)據(jù)分布Fig.6 Experimental data distribution of w and cp for CO2 mixtures

表3匯總了現(xiàn)有密度、比熱容與聲速的實(shí)驗(yàn)研究現(xiàn)狀及研究建議。

表3 CO2混合物密度、比熱容與聲速研究現(xiàn)狀及研究建議Tab.3 Experimental data description and research recommendations of theomodynamic properties of CO2 mixtures

總之，對(duì)S-CO2循環(huán)，當(dāng)前含CO2混合體系的pvTx、比熱容和聲速實(shí)驗(yàn)數(shù)據(jù)較不充分，在高溫區(qū)存在大量數(shù)據(jù)空白，同時(shí)CO2+H2O與CO2+O2的實(shí)驗(yàn)數(shù)據(jù)數(shù)量與精度均未滿足需求。針對(duì)S-CO2循環(huán)需求，對(duì)于單相區(qū)密度應(yīng)重點(diǎn)關(guān)注所有二元及多元體系，特別是當(dāng)前尚無實(shí)驗(yàn)數(shù)據(jù)體系的400 K以上區(qū)域，以及CO2+H2O體系在10 MPa以下區(qū)域的密度；對(duì)于比熱容與聲速，則需要對(duì)全部二元及多元體系開展更多實(shí)驗(yàn)研究，以填補(bǔ)現(xiàn)存的大量數(shù)據(jù)空白區(qū)域。

3 狀態(tài)方程研究現(xiàn)狀

狀態(tài)方程是計(jì)算與導(dǎo)出工質(zhì)熱力學(xué)性質(zhì)的理論模型，是S-CO2循環(huán)進(jìn)行系統(tǒng)設(shè)計(jì)與優(yōu)化的基礎(chǔ)。最簡(jiǎn)單的狀態(tài)方程是理想氣體狀態(tài)方程，仍有學(xué)者采用此方程對(duì)半閉式S-CO2循環(huán)進(jìn)行計(jì)算[71]。適用于S-CO2循環(huán)相關(guān)計(jì)算的狀態(tài)方程主要分為半理論方程和多參數(shù)方程。表4展示了近年來部分學(xué)者對(duì)S-CO2循環(huán)開展研究時(shí)采用的熱力學(xué)模型。

表4 研究半閉式S-CO2循環(huán)使用的狀態(tài)方程Tab.4 Equations of state used in studies of the semi-closed S-CO2 cycle

半理論方程的構(gòu)建基于一定的理論基礎(chǔ)，利用數(shù)據(jù)對(duì)有理論背景的參數(shù)進(jìn)行優(yōu)化，方程一般適用于多種物質(zhì)，有較好的外推性能，且可以從純凈物推廣到混合物而不對(duì)方程形式有巨大改變。目前，較多應(yīng)用于S-CO2循環(huán)計(jì)算的是立方型方程，如PR、SRK方程[73]，PR方程形式如下：

SRK方程形式如下：

式中：p為壓力，T為溫度，v為比體積，R為氣體常數(shù)，a、b為與溫度相關(guān)的系數(shù)。

多參數(shù)狀態(tài)方程是一類基于大量實(shí)驗(yàn)數(shù)據(jù)構(gòu)建而成的熱力學(xué)模型，通常含有大量無明確物理意義的參數(shù)，利用實(shí)驗(yàn)數(shù)據(jù)優(yōu)化參數(shù)與函數(shù)形式，可以在實(shí)驗(yàn)數(shù)據(jù)覆蓋的范圍內(nèi)精確復(fù)現(xiàn)實(shí)驗(yàn)數(shù)據(jù)，但外推預(yù)測(cè)性能較差。對(duì)于CO2的混合體系，目前應(yīng)用最廣泛的是LKP方程[75-76]、GERG-2008模型[18]和EOS-CG模型[77]。

LKP方程的構(gòu)建基于對(duì)比態(tài)原理，表達(dá)形式為：

式中：下標(biāo)r代表對(duì)比參數(shù)，bi、ci、di、β、γ是可調(diào)參數(shù)。對(duì)于混合體系，應(yīng)用混合規(guī)則，調(diào)節(jié)相互作用系數(shù)，即可計(jì)算熱力學(xué)性質(zhì)。

GERG-2008模型是對(duì)GERG-2004模型的改進(jìn)，以Helmholtz函數(shù)為顯式進(jìn)行表達(dá)：

式中：α為Helmholtz函數(shù)，上標(biāo)o代表理想氣體項(xiàng)，上標(biāo)r代表剩余項(xiàng)；δ、τ、x為對(duì)比密度、對(duì)比溫度與摩爾分?jǐn)?shù)。利用實(shí)驗(yàn)數(shù)據(jù)對(duì)剩余項(xiàng)進(jìn)行擬合，即可高精度復(fù)現(xiàn)實(shí)驗(yàn)數(shù)據(jù)。EOS-CG方程在GERG-2008模型的基礎(chǔ)上作進(jìn)一步優(yōu)化，將適用氣體種類從18種縮小到6種（CO2、O2、N2、H2O、Ar、CO），幾乎與S-CO2循環(huán)所需氣體組分完全相同。

針對(duì)S-CO2循環(huán)的研究需求，各類方程精度有所不同。Chowdhury等人[78]對(duì)比了PR、LKP和理想氣體狀態(tài)方程對(duì)S-CO2循環(huán)燃燒過程的計(jì)算結(jié)果，認(rèn)為在高溫高壓條件下，LKP方程計(jì)算效果更好，適合超臨界系統(tǒng)的計(jì)算。Scaccabarozzi等人[24]使用Aspen Plus軟件內(nèi)置的二元交互系數(shù)對(duì)比了PR、SRK、BWR方程后，認(rèn)為PR方程計(jì)算半閉式S-CO2循環(huán)更準(zhǔn)確，特別是CO2+H2O體系，密度性質(zhì)相對(duì)偏差小于2%，而SRK方程則在計(jì)算CO2+O2體系時(shí)更準(zhǔn)確。White等人[79]則對(duì)CO2+H2O體系展開了重點(diǎn)研究，推薦LKP方程作為計(jì)算工具。Li等人[80]研究認(rèn)為，PR方程計(jì)算CO2+N2與CO2+Ar體系的整體精確度優(yōu)于SRK方程，特別是液相區(qū)域（PR：2.07%，SRK：5.42%）。Mazzoccoli等人[81]對(duì)比了常溫下多種狀態(tài)方程對(duì)富CO2體系的熱力學(xué)性質(zhì)計(jì)算，認(rèn)為L(zhǎng)KP方程和GERG模型計(jì)算表現(xiàn)相對(duì)較優(yōu)，而PR方程在立方型方程中表現(xiàn)最好，對(duì)密度性質(zhì)的預(yù)測(cè)偏差普遍在1%以上。對(duì)于當(dāng)前精度最高多參數(shù)狀態(tài)方程的GERG-2008和EOS-CG，其聲稱在有效范圍內(nèi)pvTx性質(zhì)的不確定度在0.1%～1%，然而其擴(kuò)展有效溫度范圍為700 K，遠(yuǎn)低于S-CO2循環(huán)的最高溫度。此外，Lozano-Martín[61]、Yang[50,57]、Nazeri[82]、Souissi[59]、Wegge[70]等人的研究表明，GERG-2008和EOS-CG與實(shí)驗(yàn)數(shù)據(jù)的偏差往往會(huì)超出其聲明的不確定度，對(duì)于富CO2體系還會(huì)出現(xiàn)偏差增大的現(xiàn)象。需要注意的是，GERG-2008和EOS-CG模型在擬合參數(shù)時(shí)，由于實(shí)驗(yàn)研究的空白，沒有使用CO2+O2體系的pvtx實(shí)驗(yàn)數(shù)據(jù)，使用CO2+H2O體系pvtx實(shí)驗(yàn)數(shù)據(jù)的最高溫度也沒有超過700 K。因此，使用GERG-2008和EOS-CG模型進(jìn)行S-CO2循環(huán)的計(jì)算時(shí)，必須考慮其對(duì)特定組分以及超出有效溫度壓力范圍計(jì)算的可靠性。根據(jù)當(dāng)前狀態(tài)方程的研究現(xiàn)狀，對(duì)于溫度低于700 K的循環(huán)部件研究，推薦使用當(dāng)前計(jì)算精度最高的GERG-2008和EOS-CG模型；而當(dāng)溫度超過當(dāng)前實(shí)驗(yàn)數(shù)據(jù)的覆蓋范圍時(shí)，應(yīng)當(dāng)優(yōu)先使用具有理論背景的PR、SRK等立方型方程或LKP等維里型方程。特別地，針對(duì)重點(diǎn)組元體系，推薦使用LKP方程計(jì)算CO2+H2O體系物性，PR方程計(jì)算CO2+O2體系物性。此外，如果有基于前述數(shù)據(jù)空白所開展的實(shí)驗(yàn)研究，得到了新的二元交互系數(shù)，或者對(duì)現(xiàn)有經(jīng)驗(yàn)?zāi)Ｐ偷耐馔颇Ｐ?，?yīng)當(dāng)及時(shí)對(duì)現(xiàn)有模型進(jìn)行更新。

整體而言，目前廣泛應(yīng)用于CCUS技術(shù)的高精度多參數(shù)狀態(tài)方程并不能很好地適應(yīng)S-CO2循環(huán)的高溫高壓計(jì)算需求，而已經(jīng)有充分研究成果的立方型狀態(tài)方程的精度并不夠高。因此，狀態(tài)方程研究除了依靠新的實(shí)驗(yàn)數(shù)據(jù)對(duì)已有方程進(jìn)行優(yōu)化外，還應(yīng)考慮充分利用現(xiàn)有數(shù)據(jù)與物質(zhì)共性規(guī)律，如構(gòu)建S-CO2循環(huán)專用的基于深刻理論研究的截?cái)嗑S里型狀態(tài)方程，以達(dá)到對(duì)現(xiàn)有實(shí)驗(yàn)數(shù)據(jù)的良好外推效果。

4 結(jié)論與展望

1）相比其他性質(zhì)，CO2混合體系氣液相平衡數(shù)據(jù)較為豐富，但針對(duì)S-CO2循環(huán)研究需求，CO2+H2O體系近臨界及高溫區(qū)域氣液相平衡的實(shí)驗(yàn)研究還較缺乏；CO2的部分二元體系及多元體系的pvTx實(shí)驗(yàn)數(shù)據(jù)相當(dāng)缺乏，全部體系均無超過800 K的pvTx實(shí)驗(yàn)研究。

2）缺乏對(duì)于CO2混合體系不同組分與溫度壓力條件下的比熱容與聲速實(shí)驗(yàn)研究。

3）目前廣泛應(yīng)用的半理論與多參數(shù)狀態(tài)方程在氣相區(qū)與超臨界區(qū)存在精度及有效溫度壓力范圍不夠的缺陷，這是由混合體系高溫實(shí)驗(yàn)數(shù)據(jù)缺乏造成的，應(yīng)考慮采用能夠充分利用現(xiàn)有數(shù)據(jù)并具有良好外推性能的狀態(tài)方程形式。