工業(yè)規(guī)模CO2管道泄放過程中的壓力響應(yīng)及相態(tài)變化

喻健良，郭曉璐，閆興清，張永春，陳紹云

（1大連理工大學(xué)化工機(jī)械學(xué)院，遼寧 大連 116024；2大連理工大學(xué)化工學(xué)院，遼寧 大連 116024）

引言

CO2捕集和封存技術(shù)（carbon capture and storage，CCS）是實現(xiàn)全球二氧化碳減排的最有效方式[1]。管道運(yùn)輸是其重要環(huán)節(jié)，也是較安全、環(huán)保和經(jīng)濟(jì)的運(yùn)輸方式[2]。大量CO2管道的應(yīng)用使其安全運(yùn)行問題備受關(guān)注[3-4]。當(dāng)CO2長輸管道發(fā)生泄漏時，由于焦耳-湯姆遜效應(yīng)產(chǎn)生的低溫容易導(dǎo)致管道脆性斷裂擴(kuò)展，同時泄漏的CO2擴(kuò)散后會在附近形成高濃度區(qū)，危及區(qū)域內(nèi)人員健康和安全[5-7]。

CO2管道泄漏過程包括管內(nèi)壓力和溫度變化、射流膨脹、干冰生成和升華、重氣擴(kuò)散等[8]。而管內(nèi)壓力和溫度變化是研究低溫脆斷以及泄放區(qū)域內(nèi)CO2擴(kuò)散過程的關(guān)鍵[9-12]。國內(nèi)外對CO2管道泄漏時壓力響應(yīng)及相變的實驗研究數(shù)量不多。Ahmad等[13]介紹了一次大規(guī)模全口徑密相CO2埋地管道斷裂實驗，回路管道規(guī)格φ219.1 mm×12.7 mm，總長226.8 m，發(fā)現(xiàn)管內(nèi)壓力經(jīng)歷快速下降、準(zhǔn)靜態(tài)水平、后期泄放3個階段。Cosham等[14]介紹了3次全尺寸CO2管道斷裂實驗，管道規(guī)格φ914.0 mm×25.4 mm，長度分別為16.16、16.97、22.71 m，發(fā)現(xiàn)減壓過程都出現(xiàn)了壓力平臺，且?guī)щs質(zhì)CO2壓力平臺更高。Koeijera等[15-16]通過規(guī)格φ 12 mm×1 mm、長140 m管道裝置研究了過冷液相CO2減壓過程，發(fā)現(xiàn)液相CO2很快轉(zhuǎn)變?yōu)閮上嘈狗?，最后轉(zhuǎn)變?yōu)闅庀嘈狗?。Cosham等[17]采用規(guī)格φ168.3 mm×10.97 mm、長144 m實驗管道，研究了氣相、液相CO2及含雜質(zhì)CO2減壓行為，發(fā)現(xiàn)液相較氣相CO2的減壓曲線平臺更長，而含雜質(zhì)CO2減壓曲線平臺則又高又短。Han等[18-19]利用細(xì)長管的圓柱形容器裝置，研究發(fā)現(xiàn)細(xì)長管內(nèi)的液態(tài)CO2轉(zhuǎn)變?yōu)闅夤虄上嗔骺赡馨l(fā)生在CO2進(jìn)入大氣之前。Xie等[20]利用規(guī)格φ40 mm×5 mm、長23 m回路管道裝置進(jìn)行了CO2垂直泄放實驗，認(rèn)為超臨界CO2比氣相CO2壓降速率大得多?？傊?，目前對CO2管道泄漏時壓力響應(yīng)及相變研究的文獻(xiàn)報道并不一致，還需要進(jìn)一步研究。

本文采用的管道規(guī)格為φ273 mm×20 mm，長256 m，屬于工業(yè)規(guī)模級別，相對于文獻(xiàn)報道中的管道尺寸，更加適合和接近于模擬實際CO2運(yùn)輸管道?；诖?，對CO2管道泄放時的管內(nèi)壓力和溫度變化進(jìn)行了測量和分析，研究了管道內(nèi)壓力響應(yīng)情況以及相態(tài)變化規(guī)律，以期為CO2管道斷裂擴(kuò)展和泄漏擴(kuò)散研究提供實驗依據(jù)。

1 CO2管道實驗裝置及方案

1.1 實驗裝置

圖1為實驗裝置和現(xiàn)場圖，該裝置由主管道、雙膜爆破裝置、加熱裝置以及數(shù)據(jù)采集系統(tǒng)等組成。主管道規(guī)格φ273 mm×20 mm，長度256 m，材質(zhì)16Mn。

該裝置采用雙膜爆破裝置實現(xiàn)可控泄放，并采用不同口徑孔板實現(xiàn)不同口徑泄放；利用加熱裝置為管道升溫以達(dá)到不同的初始溫度和壓力；實驗過程中管道壓力和溫度采用分布式同步數(shù)據(jù)采集系統(tǒng)測量記錄。圖2是雙膜爆破裝置原理。該裝置由兩個爆破片及夾持器、法蘭及管道組成。通過在兩個爆破片（A和B）之間的空腔注入高壓氣體來控制爆破。該裝置長度約0.9 m，相對于整條管道來說很短，對壓力和溫度等測量數(shù)據(jù)的影響可忽略。該方法可達(dá)到快速泄壓的目的，泄壓速率范圍為1×104～4×104MPa·s-1[21]。

1.2 實驗方案

為測量管道壓力和溫度，在管道上分布了12個壓力傳感器、18個管頂熱電偶、6個管底熱電偶以及12個管外壁熱電偶。本文涉及的CO2減壓波速度為100～300 m·s-1，而管道長度256 m足夠長，因此壓力傳感器采集頻率1 kHz，精度0.06%FS～0.6%FS，可以支持本文研究內(nèi)容；熱電偶響應(yīng)時間200 ms，誤差±1℃，均滿足高壓低溫的環(huán)境條件。采用美國NI公司的CompactRIO實時控制器進(jìn)行數(shù)據(jù)采集，包括1個NI cRIO-9025嵌入式實時控制器和5個NI 9144以太網(wǎng)擴(kuò)展機(jī)箱，各機(jī)箱內(nèi)包括13個NI 9219通用模塊和3個NI 9213熱電偶輸入模塊。數(shù)據(jù)采集系統(tǒng)采用美國NI公司的LabVIEW軟件進(jìn)行編程。

1.3 實驗條件

本文選取了3組初始相態(tài)分別為氣相、氣液兩相、超臨界相態(tài)實驗進(jìn)行分析。圖3是CO2相態(tài)圖。圖中R是三相點(diǎn)，C是臨界點(diǎn)，臨界壓力Pc為7.38 MPa。RC、RB、RA分別為氣液、氣固、液固平衡線。圖示了3組泄放實驗Test 1、Test 2和Test 3的初始相態(tài)所在的位置。表1是實驗初始條件，表2為實驗所取測量點(diǎn)位置。

圖1 實驗裝置和現(xiàn)場圖 Fig.1 Schematic and scene graph of experimental apparatus

圖2 雙膜爆破裝置原理 Fig.2 Schematic diagram of dual-disc blasting device

圖3 CO2相態(tài)圖 Fig.3 CO2phase graph

2 實驗結(jié)果和討論

2.1 氣相CO2管道泄放中壓力響應(yīng)及相態(tài)變化

圖4為Test 1實驗壓力響應(yīng)曲線。由圖可知，1.48 t CO2且泄放口徑100 mm的泄放時間約55.4 s，平均壓降速率約95.7 kPa·s-1；管道減壓曲線變化過程整體上差別很小，而實際在泄放初始時有明顯不同[22]。圖4局部放大圖是管道上5個測點(diǎn)（P2、P3、P5、P6、P8）在泄放初始時的壓力響應(yīng)過程。本次實驗的減壓波傳播分為2個時間段，每個時間段為1.1 s（1.1 s是減壓波以238.2 m·s-1介質(zhì)聲速自P2點(diǎn)傳播至管道末端的時間）。在第1個1.1 s內(nèi)，當(dāng)減壓波傳播至P2時，該處壓力快速下降至4.63 MPa，降幅0.67 MPa；此時該處短暫穿過相包絡(luò)線進(jìn)入兩相區(qū)，壓力略有反彈后緩慢下降[23]。減壓波到達(dá)P3、P5、P6、P8時，壓力依次開始突降，降幅依次減小，分別為0.62、0.41、0.37、0.01 MPa。在第2個1.1 s內(nèi)，減壓波自管道末端傳播至泄放端，并依次導(dǎo)致P8、P6、P5、P3、P2再次大幅度降低，降幅依次減小，分別為0.72、0.49、0.37、0.23、0.09 MPa。經(jīng)過此次減壓波反射，管道泄放端和末端壓力差逐漸縮小，對數(shù)壓力與時間呈線性關(guān)系下降。

圖5為Test 1實驗壓力-溫度對應(yīng)曲線。泄放開始后，距離泄放端較近的T1、T2處壓力-溫度曲線向右偏離了CO2飽和線，始終處于氣相狀態(tài)；而處 于管道末端附近T3-top和T3-bot的壓力和溫度分別降至3.22 MPa和-2.9℃、3.62 MPa和1.5℃時，開始貼著氣液飽和線下降，說明該處底部和頂部CO2先后進(jìn)入了兩相區(qū)。T3-top最低溫度降至-43℃，T3-bot最低溫度降至-64℃，說明距離泄放端越遠(yuǎn)，溫降幅度越大，且底部較頂部降幅大，處于兩相區(qū)時間越長。

表1 實驗初始條件 Table 1 Experimental conditions

表2 實驗取測量點(diǎn)位置 Table 2 Experimental measurement point locations

圖4 壓力響應(yīng)曲線 Fig.4 Pressure response curve

圖5 壓力-溫度對應(yīng)曲線 Fig.5 Corresponding curve of pressure-temperature

通過以上分析，對于氣相CO2泄放，泄放端氣相CO2快速膨脹排出，使該處發(fā)生壓力突降，形成的減壓波向管道末端傳播，所經(jīng)之處壓力依次發(fā)生突降。壓力突降導(dǎo)致的氣體膨脹作用，使壓降停滯或減慢。隨著管內(nèi)壓力和溫度下降，泄放端附近CO2相態(tài)始終為氣態(tài)，而管道末端附近溫度的大幅下降使該處出現(xiàn)氣液均相CO2，且管道底部較頂部的氣液均相CO2維持時間更長。圖6為氣相CO2泄放相態(tài)變化過程。

2.2 氣液兩相CO2管道泄放中壓力響應(yīng)及相態(tài)變化

圖6 氣相CO2泄放相態(tài)變化 Fig.6 Phase transition in process of gas CO2discharge

圖7 壓力響應(yīng)曲線 Fig.7 Pressure response curve

Test 2實驗的CO2初始狀態(tài)為氣液分層，泄放口徑50 mm，氣液界面處于泄放口之下，其中液相CO2體積為2.25 m3，氣相CO2體積為8.56 m3。圖7是Test 2實驗壓力響應(yīng)曲線。由圖可知，3.32 t CO2充裝量的泄放時間約777.2 s，平均壓降速率約7.2 kPa·s-1。圖7局部放大圖是P2、P3、P5、P7、P8測點(diǎn)在泄放初始時的壓力響應(yīng)過程。圖中“→”代表減壓波自管道泄放端向末端傳播，“←”則代表減壓波自管道末端向泄放端傳播。本次實驗減壓波傳播過程分為8個時間段，每個時間段為1.2 s（1.2 s是減壓波以198.5 m·s-1介質(zhì)聲速自P2處傳播至管道末端的時間）。表3是P2、P3、P5、P7、P8測量點(diǎn)的壓力響應(yīng)參量。表中包含3種壓力響應(yīng)參量：壓力突降幅度ΔPf、壓力反彈幅度ΔPr和準(zhǔn)靜態(tài)壓力水平Pqs。在第1個1.2 s內(nèi)，當(dāng)減壓波依次傳播至P2、P3、P5、P7、P8時，各處發(fā)生第1次壓力突降和反彈，達(dá)到第1個準(zhǔn)靜態(tài)壓力水平Pqs1。ΔPf1和ΔPr1隨距離增大而逐次減小，Pqs1隨距離增大而逐次降低。在第2個1.2 s內(nèi)，減壓波在管道末端發(fā)生反射，向泄放端開始傳播，所經(jīng)各處導(dǎo)致第2次壓力突降和反彈后達(dá)到第2個準(zhǔn)靜態(tài)壓力水平Pqs2。ΔPf2和ΔPr1隨距離增大而逐次增大，Pqs2隨距離增大而逐次升高。第3次時間段1.2 s的減壓波傳播過程與第1次情況相同。從表中數(shù)據(jù)看出，ΔPf和ΔPr隨傳播次數(shù)增加整體不斷減小，各個測點(diǎn)準(zhǔn)靜態(tài)壓力水平逐漸接近。當(dāng)減壓波開始第4次傳播時，由于管道泄放速率已大于氣泡生成率，各處壓力不再發(fā)生反彈；并且隨著減壓波傳播次數(shù)增加，各處壓力越來越接近。

圖8為Test 2實驗壓力-溫度對應(yīng)曲線。泄放開始時，經(jīng)過減壓波多次反射，所有測點(diǎn)處CO2由氣液分層轉(zhuǎn)變?yōu)闅庖壕?；?dāng)T1-top、T2-top、T3-top、T1-bot、T2-bot、T3-bot依次降至對應(yīng)壓力和溫度4.7 MPa和14℃、4.7 MPa和14℃、4.2 MPa和9℃、2.9 MPa和-4℃、2.2 MPa和-15℃、1.1 MPa和-38℃時，依次脫離氣液飽和線，即由氣液均相變?yōu)闅庀?。氣液兩相與氣相CO2泄放管內(nèi)溫度變化具有相似規(guī)律。

通過分析得知此次氣液兩相CO2泄放，管道泄放端上層氣相CO2首先快速膨脹排出，使管道上層壓力突降，形成的減壓波向管道末端傳播，同時向管道底部液相CO2傳播。當(dāng)減壓波進(jìn)入液相CO2時，飽和液態(tài)CO2過熱汽化，使壓力突降停滯并發(fā)生反彈。減壓波沿管道的多次傳播使整條管道底部液相CO2變?yōu)檫^熱液體而沸騰汽化為氣液均相。隨著CO2流出和壓力下降，氣液均相向氣相CO2的轉(zhuǎn)變從管道泄放端頂部開始，向管道末端頂部傳播，而后向管道底部傳播，直到整條管道轉(zhuǎn)變?yōu)闅庀嘈狗拧D9為氣液兩相CO2泄放相態(tài)變化過程。

2.3 超臨界CO2管道泄放中壓力響應(yīng)及相態(tài)變化

圖10是Test 3實驗的壓力響應(yīng)曲線。由圖可知，5.87 t CO2且泄放口徑50 mm的泄放時間約312.4 s，平均壓降速率約25.3 kPa·s-1。圖10局部放大圖是泄放開始時管道上P1、P2、P3、P4、P5、P6測點(diǎn)的壓力響應(yīng)過程。本次實驗減壓波傳播過程分為2個時間段，每個時間段為1.3 s（1.3s是減壓波以198.5 m·s-1介質(zhì)聲速自P1傳播至管道末端的時間）。在第一個1.3 s內(nèi)，P1處壓力突降至臨界壓力Pc（7.38 MPa）時發(fā)生氣泡成核使壓降停滯或減慢，停滯時間約0.015 s；壓力繼續(xù)突降至7.035 MPa時快速反彈；反彈壓力超過Pc時，壓升速率減慢。減壓波傳播至P2時的壓降過程與P1處相似，但到達(dá)Pc時壓降停滯時間較長，約0.025 s。當(dāng)減壓波傳播至P3和P4時，該兩處壓力都是在Pc處開始反彈，而P5和P6處的壓力反彈點(diǎn)則要高于Pc，故不會發(fā)生氣泡成核；該處壓力突降和反彈是由前端傳播來的減壓波和壓力反彈波決定的。減壓波經(jīng)過P1、P2、P3、P4、P5、P6時導(dǎo)致的壓降幅度依次減小，分別為0.875、0.859、0.535、0.535、0.517、0.432 MPa；在壓力反彈之后，所有測點(diǎn)壓力幾乎達(dá)到同一個準(zhǔn)靜態(tài)壓力水平7.62 MPa，這個壓力水平是由泄放端的壓力反彈決定的。在第2個1.3 s內(nèi)，減壓波自管道末端向泄放端傳播，并依次導(dǎo)致P6、P5、P3、P3、P2、P1再次大幅度降低，且經(jīng)過Pc時，壓降出現(xiàn)停滯或減緩。經(jīng)過此次減壓波反射，管內(nèi)壓力逐漸低于Pc，CO2相態(tài)由超臨界相態(tài)轉(zhuǎn)變?yōu)闅庖壕?，不再發(fā)生壓力突降和反彈。

表3 壓力響應(yīng)參量 Table 3 Pressure response parameters/MPa

圖8 壓力-溫度對應(yīng)曲線 Fig.8 Corresponding curve of pressure-temperature

圖9 管道內(nèi)氣液兩相CO2泄放相態(tài)變化 Fig.9 Phase transition in process of gas-liquid CO2discharge

圖10 壓力響應(yīng)曲線 Fig.10 Pressure response curve

圖11 第1個時間段1.3 s內(nèi)的壓降速率 Fig.11 Pressure drop rate of first 1.3 s

圖12 壓力-溫度對應(yīng)曲線 Fig.12 Corresponding curve of pressure-temperature

圖11是第1個時間段1.3 s內(nèi)壓變速率曲線。從圖中得知，P1處壓變速率降至-51.358 MPa·s-1時，由于壓力穿過Pc，壓變速率發(fā)生突變，回復(fù)至零附近，這時壓降停滯；由于泄壓速率高于氣泡生成率，導(dǎo)致該處壓變速率再次降低至-68.536 MPa·s-1，這時氣泡生成率已高于泄壓速率，壓力以指數(shù)增長；壓變速率增長至26.40 MPa·s-1時，壓力再次穿過Pc，使該處壓變速率減小至3.340 MPa·s-1附近，而后增長至7.758 MPa·s-1。P2處壓變速率變化和P1處相似，但兩次壓變速率降幅更小，分別為44.568 MPa·s-1和44.071 MPa·s-1；兩次壓變速率增幅也更小，分別為16.187 MPa·s-1和9.397 MPa·s-1。P3、P4、P5、P6處壓力突降沒有穿過Pc，因此在壓變速率突降和增長中都只有一個降幅和一個增幅；壓變速率降幅隨距離增長依次變小，分別為9.245、5.586、5.346、3.741 MPa·s-1；壓變速率增幅隨距離增長依次變小，分別為8.946、7.955、7.720、6.5462 MPa·s-1。由分析可知，泄壓速率和氣泡生成率是影響減壓波傳播過程的重要因素。

圖12為Test 3實驗壓力-溫度對應(yīng)曲線。泄放開始的CO2仍處于超臨界狀態(tài)，當(dāng)壓力低于Pc時，所有測點(diǎn)壓力和溫度沿著飽和線下降，說明超臨界CO2很快轉(zhuǎn)變?yōu)闅庖壕唷1-top、T2-top、T3-top、T1-bot、T2-bot、T3-bot依次脫離氣液飽和線，即由氣液兩相轉(zhuǎn)變?yōu)闅庀嗨鶎?yīng)的壓力和溫度點(diǎn)分別為5 MPa和15.4℃、4.7 MPa和13℃、3.7 MPa和-4.4℃、3 MPa和-4.2℃、2.9 MPa和-6℃、0.5 MPa和-39℃。由此可知，超臨界CO2在壓力降至Pc時轉(zhuǎn)變?yōu)闅庖壕?，隨著管內(nèi)壓力和溫度下降，氣液均相向氣相CO2的轉(zhuǎn)變先從管道頂部開始，而后向管道底部傳播。超臨界CO2與氣液兩相以及氣相CO2管道泄放中的管內(nèi)溫度變化具有相似規(guī)律。

通過以上分析，將超臨界CO2管道泄放過程分為超臨界CO2泄放、氣液均相CO2泄放和氣相CO2泄放3個過程。在Pc之上時進(jìn)行的是超臨界CO2泄放，發(fā)生了壓力突降和反彈、準(zhǔn)靜態(tài)壓力水平以及再次壓降和反彈等過程。當(dāng)壓力普遍低于Pc，超臨界CO2全部轉(zhuǎn)變?yōu)闅庖壕郈O2。隨著管內(nèi)壓力和溫度下降，氣液均相CO2向氣相CO2的轉(zhuǎn)變，從管道頂部開始，而后從管道底部傳播。圖13為超臨界CO2泄放相態(tài)變化過程。

圖13 管道內(nèi)超臨界CO2泄放相態(tài)變化 Fig.13 Phase transition in process of supercritical CO2discharge

3 結(jié)論

（1）氣相CO2管道泄放中，減壓波前沿所經(jīng)之處，各點(diǎn)壓力立即開始下降，壓降之后，由于氣體膨脹作用，使壓降發(fā)生停滯或減慢。減壓波傳播使管道壓力和溫度不斷下降，泄放端附近CO2相態(tài)始終為氣態(tài)；而管道末端附近則出現(xiàn)氣液均相CO2，且底部較頂部氣液均相CO2維持時間更長。

（2）泄放口位于氣液界面之上的氣液兩相CO2管道泄放中，減壓波傳播發(fā)生多次反射并導(dǎo)致多次壓力突降和反彈。減壓波自管道上層向底部傳播使管底液相CO2變?yōu)檫^熱液體而汽化，使壓降停滯并反彈；減壓波自泄放端向管道末端傳播使整條管道的氣液分層CO2轉(zhuǎn)變?yōu)闅庖壕唷庖壕郈O2向氣相的轉(zhuǎn)變由管道頂部開始，而后傳播至底部。

（3）超臨界CO2管道泄放中，壓力突降和反彈以及準(zhǔn)靜態(tài)壓力水平發(fā)生在臨界區(qū)域附近，且壓力穿過Pc時，壓變速率會停滯和減慢。管道壓力高于Pc時進(jìn)行的是超臨界CO2泄放；當(dāng)管道壓力普遍低于Pc時，超臨界CO2全部轉(zhuǎn)變?yōu)闅庖壕郈O2進(jìn)行泄放。氣液均相CO2向氣相的轉(zhuǎn)變由管道頂部開始，而后傳播至底部。

[1] IEA.Energy Technology Perspectives 2012: Pathways to a Clean Energy System[R].2012

[2] Wareing C J, Fairweather M, Falle S A E G, Woolley R M.Validation of a model of gas and dense phase CO2jet releases for carbon capture and storage application [J].International Journal of Greenhouse Gas Control, 2014, 20: 254-271.

[3] Koornneef J, Spruijt M, Molag M, Ramirez A, Faaij A, Turkenburg W.Uncertainties in risk assessment of CO2pipelines [J].Energy Procedia, 2009, 1: 1587-1594.

[4] Duncan I J, Wang H.Estimating the likelihood of pipeline failure in CO2transmission pipelines: new insights on risks of carbon capture and storage [J].International Journal of Greenhouse Gas Control, 2014, 21: 49-60.

[5] Woolley R M, Fairweather M, Wareing C J, et al.An integrated, multi-scale modelling approach for the simulation of multiphase dispersion from accidental CO2pipeline releases in realistic terrain [J].International Journal of Greenhouse Gas Control, 2014, 27: 221-238.

[6] Woolley R M, Fairweather M, Wareing C J, et al.CO2PipeHaz: quantitative hazard assessment for next generation CO2pipelines [J].Energy Procedia, 2014, 63: 2510-2529.

[7] Rian K E, Grimsmo B, Laksa B, Vembe B E, Lilleheie N I, Brox E, Evanger T.Advanced CO2dispersion simulation technology for improved CCS safety [J].Energy Procedia, 2014, 63: 2596-2609.

[8] Molag M, Dam C.Modelling of accidental from a high pressure CO2pipelines [J].Energy Procedia, 2011, 4: 2301-2307.

[9] Lund H, Flatten T, Munkejord S T.Depressurization of carbon dioxide in pipelines — models and methods [J].Energy Procedia, 2011, 4: 2984-2991.

[10] Aursand E, Aursand P, Berstad T, D?rum C, Hammer M, Munkejord S T, Nordhagen H O.CO2pipeline integrity: a coupled fluid-structure model using a reference equation of state for CO2[J].Energy Procedia, 2013, 37: 3113-3122.

[11] Martynov S, Brown S, Mahgerefteh H, Sundara V.Modelling choked flow for CO2from the dense phase to below the triple point [J].International Journal of Greenhouse Gas Control, 2013, 19: 552-558.

[12] Brown S, Martynov S, Mahgerefteh H, Chen Shaoyun, Zhang Yongchun.Modelling the non-equilibrium two-phase flow during depressurisation of CO2pipelines [J].International Journal of Greenhouse Gas Control, 2014, 30: 9-18.

[13] Ahmad M, Lowesmith B, Koeijer G D, Nilsen S, Tonda H, Spinelli C, Cooper R, Clausen S, Mendes R, Florisson O.COSHER joint industry project: large scale pipeline rupture tests to study CO2release and dispersion [J].International Journal of Greenhouse Gas Control, 2015, 37: 340-353.

[14] Cosham A, Jones D G, Armstrong K, Allason D, Barnett J.Ruptures in gas pipelines, liquid pipelines and dense phase carbon dioxide pipelines//Proceedings of the 2012 9th International Pipeline Conference[C].2012.

[15] Koeijera G d, Borch J H, Jakobsenb J, Drescher M.Experiments and modeling of two-phase transient flow during CO2pipeline depressurization [J].Energy Procedia, 2009, (1): 1683-1689.

[16] Drescher M, Varholm K, Munkejord S T, Hammer M, Held R, Koeijer G d, Barnett J.Experiments and modelling of two-phase transient flow during pipeline depressurization of CO2with various N2compositions [J].Energy Procedia, 2014, 63: 2448-2457.

[17] Cosham A, Jones D G, Armstrong K, Allason D, Barnett J.The decompression behaviour of carbon dioxide in the dense phase// Proceedings of the 2012 9th International Pipeline Conference[C].2012.

[18] Han S H, Kim J, Chang D.An experimental investigation of liquid CO2release through a capillary tube [J].Energy Procedia, 2013, 37: 4724-4730.

[19] Han S H, Chang D, Kim J, Chang W.Experimental investigation of the flow characteristics of jettisoning in a CO2carrier [J].Process Safety and Environmental Protection, 2014, 92: 60-69.

[20] Xie Qiyuan, Tu Ran, Jiang Xi, Li Kang, Zhou Xuejin.The leakage behavior of supercritical CO2flow in an experimental pipeline system [J].Applied Energy, 2014, 130: 574-580.

[21] Bartak J.A study of the rapid depressurization of hot water and the dynamics of vapour bubble generation in superheated water [J].International Journal of Multiphase Flow, 1990, 16(5): 789-798.

[22] Martynov S, Brown S, Mahgerefteh H, Sundara V, Chen Shaoyun, Zhang Yongchun.Modelling three-phase releases of carbon dioxide from high-pressure pipelines [J].Process Safety and Environmental Protection, 2014, 92: 36-46.

[23] Botros K K, Studzinski W, Geerligs J, Glover A.Determination of decompression wave speed in rich gas mixtures [J].The Canadian Journal of Chemical Engineering, 2004, 82: 880-891.

[24] Brown S, Beck J, Mahgerefteh H, Fraga E S.Global sensitivity analysis of the impact of impurities on CO2pipeline failure [J].Reliability Engineering and System Safety, 2013, 115: 43-54.