LNG接收終端工藝虛擬仿真實驗設計

朱建魯， 宋存永， 李玉星， 李其銳， 宋 釗

(1.中國石油大學(華東) 油氣儲運工程國家級虛擬仿真實驗教學中心，山東 青島 266580；2.中石化石油工程設計有限公司數字化工程設計所，山東 東營 257026)

0 引 言

將專業前沿的工程實踐情景融入教學環節，使學生有機會接觸工程實踐、應用所學專業知識[1-3]，有助于營造良好的科研氛圍，增強學生的研究能力、創新意識和專業知識掌握能力，激發學生投身科研的熱情[4-6]。

隨著人們環保意識的增強和經濟社會的發展，作為一種新型環保能源，液化天然氣(Liquefied Natural Gas, LNG)越來越受歡迎[7]。近年來，我國沿海地區 LNG 接收站建設和開發迅速發展，培養一批具有扎實的理論功底、熟練的實踐操作能力，能夠適應LNG接收站設計、生產、建設、管理等崗位需要的高能力應用型人才[8-10]，成為了石油高等院校的首要任務之一。由于LNG接收終端儲存有大量易燃易爆的LNG，一旦有違反安全操作規程的行為，將會引起爆炸等巨大安全事故，危險性極高，因此存在“實習無實操”的現象[11]。利用虛擬仿真實驗對LNG接收終端工藝流程進行模擬，不僅可以使學生形象、生動地理解和熟悉工藝流程，還能夠有效地提高學生的實操技能，為走向工作崗位、保證LNG接收終端安全穩定運行打下堅實的基礎[9,12]。

以油氣儲運工程國家級虛擬仿真實驗教學中心(以下簡稱“實驗教學中心”)為依托，充分利用山東省油氣儲運安全重點實驗室研究成果，緊密結合LNG接收終端現場實際，改革專業實驗教學方法，完善專業實驗教學模式，采用業內通用的化工過程軟件Aspen HYSYS，進行適用于本科教學中“軟件實訓”課程的實驗設計，提高教學效果。

1 LNG接收終端工藝流程

LNG接收終端主要是接受外來LNG、儲存和外輸天然氣，包括LNG接收裝置、LNG儲罐、BOG壓縮機、BOG再冷凝器、氣化器等關鍵設備。圖1所示為典型的LNG接收終端工藝流程[13-14]。

圖1 LNG接收終端工藝流程

LNG運輸船停靠LNG接收終端碼頭后，通過船上的卸料泵將LNG從船艙排出，經過卸料臂、卸船管線和再循環管線將LNG輸送到岸上的LNG儲罐存儲。同時，為了維持LNG船艙的壓力，岸上的BOG通過天然氣回氣臂返回艙內。LNG一部分由罐內潛液泵直接進入高壓外輸泵，另一部分與經BOG壓縮機壓縮后的BOG混合進入再冷凝器，將BOG冷凝為LNG，之后經外輸高壓泵升壓后分別進入氣化器氣化，氣化后的天然氣經計量后通過輸氣管道送至用戶[15-16]。

2 LNG接收終端動態仿真的理論基礎

在Aspen HYSYS中，儲罐、分離器和氣化器采用集中參數法進行模型簡化，并忽略動量的變化，假設流體物性參數在空間是相等的，在分析中只考慮時間梯度，利用常微分方程進行描述[14]。壓縮機、泵和閥門的響應時間與其他設備相比要小很多，可認為是準穩態部件，按照穩態模型處理[17]。

2.1 儲罐和分離器動態模型

質量守恒計算公式：

(1)

能量守恒計算公式：

(2)

式中：下標1為進口，下標2為出口；V為儲罐或分離器的容積，m3；h為流體比焓，J·kg-1，其中下標V為氣相，下標L為液相；ρ為流體密度，kg·m-3；W為流量，kg·s-1。

2.2 壓縮機動態模型

能頭動態計算公式：

(3)

溫度計算公式：

(4)

式中：p1為進口壓力，Pa；p2為出口壓力，Pa；ρ1為設備入口流體密度，kg·m-3；n為多變系數；ηP為多變效率。

2.3 泵動態模型

能頭動態計算公式：

(5)

式中,η為泵效率。

2.4 氣化器動態模型

流量守恒計算公式：

(6)

能量守恒計算公式：

(7)

式中：CV為流通系數，kg·s-1[(Pa·kg·m-3)-1]1/2；Q為氣化器換熱量，J/s。

2.5 節流閥動態模型

(8)

式中：k為閥門開度，無量綱。

3 虛擬仿真實驗設計的教學目標

各個物理量參數的變化對LNG接收終端流程的性能都具有重要的影響。本實驗選取3個參數——外輸天然氣負荷、儲罐內壓力、壓縮機功率，對其進行擾動，以此得到流程中幾個重要部位的響應曲線，并觀察幾個控制的實現。其中：

(1) 根據天然氣用戶的用氣量波動頻繁、不同時段變化量大的特點，實驗中通過模擬用戶一天用氣量的波動，研究外輸天然氣負荷變化對LNG接收終端流程的性能影響；

(2) 環境熱量不可避免地傳遞到儲罐內，被LNG液體吸收后蒸發成BOG，從而導致罐內壓力增高，實驗中通過對儲罐增加熱損失，研究儲罐內壓力變化對LNG接收終端流程的性能影響；

(3) 由于在啟動流程時壓縮機的功率會忽然增加，實驗中通過增加壓縮機功率擾動，研究其對LNG接收終端流程的性能影響。

4 LNG接收終端動態仿真模型的建立

根據工程實際，實驗設計時中通過設定設備的尺寸、選擇邊界條件等過程，建立了LNG接收終端動態仿真模型。模型中設立了兩套汽化系統，浸沒燃燒式氣化器作為備用。為了使流程更合理運行，需添加必要的控制系統。LNG接收終端動態仿真模型如圖2所示。模擬邊界條件為：CH497.0 mol%、C2H62.2 mol%、C3H80.4 mol%、N20.2 mol%,外輸天然氣壓力為8.5 MPa,流量為960 t/h;選用PR狀態方程。

圖2 LNG接收終端動態流程圖

4.1 安全閥的控制

當壓力達到安全閥的設定值時，安全閥打開，排出氣體，降低儲罐壓力，防止儲罐內壓力過大。安全閥的設定參數如下，開啟設定壓力：119 kPa；閥門全打開時的壓力：119.5 kPa；孔口截面積：9.00×10-2m2；閥門關閉壓力：118.5 kPa，回座壓力：118 kPa。

4.2 BOG再冷凝器的控制

4.2.1 BOG再冷凝器入口LNG流量的控制

為冷凝BOG所需的LNG流量由控制器FC1控制，控制點設定在進入再冷凝器的LNG管路上。FC1控制選用遠程、反作用方式，調整比例參數Kp為0.1，積分參數Ti為0.5。

4.2.2 BOG再冷凝器的壓力控制

BOG再冷凝器的高壓控制由控制器PC來完成。當壓力過高時，打開閥門，釋放天然氣到BOG總管。PC控制選用遠程、正作用方式，壓力設定值為620 kPa，壓力設定范圍是0～700 kPa，調整參數Kp為1.0；Ti為0.5。

4.2.3 BOG再冷凝器的液位控制

當再冷凝器的液位高出設定值60%時，通過控制器LC打開閥門，使天然氣進入再冷凝器，再冷凝器的壓力增高，從而其底部出口流速增加、再冷凝器液位降低。LC采用直接作用方式，液位設定范圍是(40～90)%，調整參數Kp為1.0，Ti為0.5。

4.3 外輸泵最小流量控制

為了控制外輸泵最小流量，模型中提供一條回流總管，由控制器FC2調節，回流返回再冷凝器。FC2控制選用遠程、反作用方式，流量設定值為16.74×106mol/h，流量設定范圍是(1～20)×106mol/h，調整參數Kp為0.005，Ti為0.1。

4.4 海水泵流量控制

海水泵流量控制由控制器FC3完成，在達到氣化的效果上，盡量減少能量的浪費。控制點設在海水泵出口管線上。FC3控制選用遠程、反作用方式，流量設定值為5.323×108mol/h；流量設定范圍是(3～6)×108mol/h，調節參數Kp為0.5，Ti為0.1。

4.5 外輸流量控制

輸送至每臺氣化器的LNG流量由控制器FC4完成，其設定值根據下游用氣需求量確定。控制點設在氣化器LNG進口管線上。FC4控制選用遠程、反作用方式，流量設定值為17×106mol/h，流量設定范圍是(15～18)×106mol/h，調節參數Kp為1.0，Ti為0.5。

5 LNG接收終端虛擬仿真實驗結果與分析

5.1 外輸天然氣負荷擾動對系統動態的影響

在模擬時，根據用戶的日用氣不均勻性來進行設定研究，得到了外輸天然氣負荷及外輸天然氣壓力的變化曲線，如圖3所示。由圖中可以看出，外輸天然氣輸出量變化曲線為類正弦曲線，其外輸壓力變化趨勢相反。

圖4所示為氣化后的LNG溫度和經過海水泵的海水流量的響應曲線。從圖中可以看出，對氣化LNG所需海水的流量應根據天然氣流量變化而變化，但是由于控制器FC3的作用，海水流量是個恒定值。當天然氣流量增加時，海水流量不變，則海水得到更深度的冷凝，所以海水溫度降低，反之，海水溫度升高。另外，經過氣化器后的LNG溫度變化與流量變化呈相反的變化趨勢。這是由于當LNG流量增加時，氣化器換熱量基本不變，所以溫度相應減小。

(a) 海水響應

(b) LNG溫度響應

圖4 氣化流程的響應曲線

圖5所示為再冷凝器入口BOG和LNG的響應。從圖中可以看出，流入再冷凝器的BOG流量以平滑的正弦曲線趨勢變化，這是由于外輸LNG流量變化引起再冷凝器壓力變化造成的。LNG流量變化有相同的變化趨勢，由此可看出遠程控制FC1起到了作用。

圖5 再冷凝器入口BOG與LNG的流量響應曲線

圖6所示為冷凝器液位控制和天然氣回流的響應。從圖中可以看出，在控制器FC1的控制下，再冷凝器的液位一直處于60%左右，變化幅度很小。當液位高于60%時，回流閥門打開，由于閥不能迅速閉合，所以NG回路中一直有流量，其流量大小由閥門的開度決定。

圖6 再冷凝器液位控制與NG回流響應曲線

5.2 儲罐內壓力擾動對系統動態的影響

實驗中，在t=20 min時，給儲罐添加2.82×107kJ/h的熱損失。

圖7所示為添加擾動后儲存系統的響應。從圖中可以看出，添加了熱損失后儲罐溫度壓力不斷上升，由于BOG壓縮機功率恒定，BOG流量變化較小，壓縮機出口壓力增大，使得罐內泵出口背壓升高，罐內泵功率不變，排出的LNG流量不斷下降。在t=60 min時刻，儲罐壓力達到119 kPa，安全閥打開，儲罐壓力上升速度減小，在t=95 min時刻，儲罐壓力恒定，維持在119.3 kPa左右，整個流程也再次達到穩定。

圖8所示為添加熱損失后，再冷凝工藝的響應。從圖中可以看出，隨著擾動的加入，再冷凝器入口BOG溫度降低0.4 ℃，壓力升高1.2 kPa。這是由于，BOG流量增加，同時壓縮機功率不變，所以壓縮機出口壓力即進入再冷凝器的BOG壓力增加。在液位控制器LC的作用下，再冷凝器的液位維持恒定，其出口流量增加。

5.3 壓縮機功率擾動對系統動態的影響

實驗中，壓縮機功率在t=10 min時刻以W設=1.437×105kJ/h為基點，添加了200%的擾動。

(a) 儲罐的響應

(b) 儲罐出口BOG與LNG流量響應

(a) 再冷凝器入口BOG溫度和壓力響應

(b) 再冷凝器壓力和液位控制響應

圖9所示為壓縮機出口BOG在添加了擾動后的響應。從圖中可以看出，隨著擾動的加入，壓縮機出口BOG流量和壓力都迅速增加。其中，流量增加到1 178 kmol/h，增幅182%，壓力增加后又減小，最后穩定在620 kPa，稍大于原來穩定值。系統在t=18 min時刻，再次達到平衡狀態。

圖10所示為再冷凝系統在添加了擾動后的響應。從圖中可以看出，擾動加入后，在控制器FC1的作用下，再冷凝器入口LNG流量與壓縮機出口BOG流量的變化趨勢相同。由于入口BOG與LNG流量的增加，尤其是BOG流量的迅速增加，導致BOG再冷凝器壓力的忽然升高，且變化幅度較大，從而導致再冷凝器出口流量的增加，再冷凝器液位的降低。隨后，控制器PC打開閥門釋放BOG，t=12 min時刻再冷凝器壓力下降，在出口和入口LNG流量的共同作用下，再冷凝器的壓力經過波動之后再次達到平衡，平衡壓力較擾動前提高了1%。同時，液位控制器LC調節LNG出口閥門，使再冷凝器的液位恢復到設定值。

圖9 壓縮機功率增加200%擾動后壓縮機出口BOG響應

(a) 再冷凝器入口LNG流量響應

(b) 再冷凝器壓力與液位控制響應

圖11所示為儲存系統在添加了擾動后的響應。從圖中可以看出，隨時間的變化，儲罐內壓力小幅度下降，溫度小幅度升高。當添加了擾動后，儲罐出口的BOG流量隨壓縮機出口BOG流量的變化而變化，8 min后的穩定值較之前提高了182%。因罐內潛液泵功率不變，所以儲罐出口的LNG流量的變化與再冷凝器入口LNG的壓力有相反的變化趨勢，穩定后較擾動前的流量下降1.07%。

(a) 儲罐響應

(b) 儲罐出口BOG與LNG流量響應

6 結 語

通過建立LNG接收終端的動態模型，設計了LNG接收終端虛擬仿真實驗，分析了外輸天然氣負荷、儲罐壓力和壓縮機功率發生擾動時系統的動態特性：外輸天然氣負荷變化時，再冷凝器的液位能夠控制在60%左右，入口LNG流量隨入口BOG流量變化而變化；當對儲罐增加了熱損失擾動時，工藝系統在75 min后再次達到穩定；當對壓縮機功率增加了200%的擾動時，工藝系統能夠在8 min后迅速穩定，體現了實驗中所設置的控制器都起到了很好的控制作用，從而證實了LNG接收終端的動態特性良好。通過虛擬仿真實驗，學生可以學習Aspen HYSYS建模、動態工藝仿真計算、工況研究分析等方法，通過虛擬仿真軟件研究實際問題。這種將軟件實操與工程實際相結合的教學模式，有利于提高教學質量，進一步培養學生學習新的思維方式和技術手段。