LNG接收站冷排水的溫降及余氯對水環境影響的數值模擬
——以湄洲灣東吳港區為例

王璐，熊樂航，張遠，陳焰，馬淑芹，張曉嬌，夏瑞*

1.中國環境科學研究院 2.北京師范大學 3.中海石油環保服務(天津)有限公司

近年來，我國能源結構不斷調整，清潔能源天然氣在一次能源消費中的占比持續提升，以液化天然氣(liquefied natural gas，LNG)方式進口的天然氣數量逐年增加[1]。LNG接收站作為液化天然氣的接收終端，承擔了液化天然氣的接收、存儲和氣化等功能。為了滿足市場燃氣需求，沿海地區已陸續建成多座LNG接收站[2]。海水作為最廉價又最容易獲取的流體介質，被大量用于LNG 的氣化換熱和消防[3]。LNG接收站運行過程中，使用海水作為加熱介質對LNG進行氣化，同時為了防止海洋生物在設備中的附著，投入定量的次氯酸鈉。因此，換熱后的海水溫度有所降低，且水體中含有余氯，這會對附近海域水環境造成一定的影響，主要包括冷排水的排放對海域溫度場的影響和余氯排放對海域生物環境的影響。LNG接收站冷排水對附近海域水溫及海洋生物影響的研究相對較多。如張繼周[4]研究了LNG接收站冷排水的排放對浮游生物、魚類、蝦類、貝類等海洋生物的定性影響，但沒有運用模型計算進行定量研究；柏育材等[5]在實驗室條件下模擬某LNG接收站冷排水中余氯對大黃魚仔魚和幼魚的毒理效應，結果表明，LNG工程余氯排放短期內對海洋魚類影響不顯著；奚泉等[6]通過三維水流溫度場模擬廣東珠海LNG 接收站在運行時溫度場的分布以及不同等溫線的影響范圍，結果表明，取水口附近為低溫帶，取水口應盡量遠離岸線。另外，楊青云等[7]運用ECOMSED三維水環境數學模型模擬大鵬灣LNG接收站附近海域溫度場分布情況，得出底層溫降區域面積大于表層的結論，且溫降區域面積在-5～0 ℃范圍內逐漸增大；唐俊逸等[8]利用主成分分析法對大鵬灣LNG接收站附近海域的水質狀況進行了分析，表明水體環境質量幾乎不受影響；畢亞梅等[9-10]研究了余氯造成的漁業資源損失量，計算得出余氯排放對工程附近海域浮游植物損失量理論最大值為17.71×1010個、浮游動物損失量理論最大值為8.06×106個。綜上，已有的對LNG接收站冷排水的溫降和余氯對附近海域水環境影響范圍及排水口位置選擇的研究較少。

隨著科學技術的發展，水環境模擬軟件不斷被開發利用于水環境預測，如水動力學模型MIKE[11]、三維水動力-水質模型系統(Delft3D)[12]、非結構網格海洋環流與生態模型(FVCOM)[13]等。如崔丹等[14]利用二維水動力-水質數值耦合模型(MIKE21)對工業園區擬建排污口排放污染物質的輸移擴散進行數值模擬研究；雷曉玲等[15]利用Delft3D數學模型對三峽航道環保疏浚水質變化進行了數值模擬。湄洲灣東吳港區擬建設哈納斯莆田LNG接收站碼頭，這將對附近海域的水環境造成一定的影響，但是否會對養殖區、湄洲島生態系統重點保護區產生不利影響，尚需要進一步模擬預測。筆者采用大海域套小海域、工程區加密的方法，運用丹麥水利研究所(DHI)研發的MIKE21軟件[11]對工程附近的海域進行平面二維的數值模擬計算。首先計算該海域的流速場，明確其潮流的運動特性，進而預測LNG接收站冷排水的溫降及余氯對附近海域的影響范圍，判斷排水口位置的合理性，以期為LNG接收站工程作業及位置選取提供科學依據。

1 研究區概況

湄洲灣港(118°51′E～119°10′E，24°56′N～25°18′N)位于我國東南沿海中部，北毗福州港，南鄰廈門港，與臺灣隔海相望，與東吳港區、羅嶼港區、秀嶼港區、莆頭港區進行海鐵聯運，構成東南最大的物流中心樞紐。哈納斯莆田LNG接收站碼頭建設于湄洲灣東吳港區(119°4′32.41″E，25°6′21.05″N)，總用海面積為159.117 5萬m2，其中工作船碼頭面積為0.23萬m2，接卸LNG 280萬t/a。排水口冷排水的設計流量為15 000 m3/h(4.17 m3/s)。LNG排水口分別距養殖區、湄洲島生態系統重點保護區1.0、3.2 km(圖1)。

2 模型的建立與驗證

采用大海域嵌套小海域、工程區加密的方法，建立LNG接收站附近海域二維水動力場，在模型模擬可靠的基礎上，分別計算得到夏季、冬季大潮時湄洲灣東吳港區LNG冷排水溫降和余氯濃度包絡線，分析冷排水對水環境及周圍海域敏感目標的影響，判斷排水口位置設置的合理性。

2.1 計算域及網格、水深

計算域模擬范圍及網格、水深情況見圖2。建立的海域數學模型計算域為118.62°E～119.56°E，24.77°N～25.29°N，即為圖2中A、B、C點以及岸線圍成的海域。模擬采用非結構三角網格，整個區域由 10 510 個節點和 19 971 個三角單元組成，在LNG接收站附近進行加密，加密后網格空間步長約20～50 m。水深摘自中國人民解放軍海軍航海保證部制作的湄洲灣最新水深地形海圖(編號14171，2012年監測，比例為1∶40 000)、南日群島至深滬灣水深地形海圖(編號14170，2010年監測，比例為1∶120 000)，LNG接收站附近水深根據福建省港航管理局勘探中心于2015年7月勘測的1∶5 000的地形圖繪制。

圖2 計算域網格及水深Fig.2 Calculate field grid and water depth

2.2 參數設置

根據湄洲灣海區口外的崇武海洋站的多年潮位觀測資料及1984年4—6月東吳潮位站3個月的潮位觀測結果，潮汐性質屬正規半日潮，灣內外潮時幾乎一致，各地潮位基本上同漲同落，高、低潮出現時間同步。模擬時長為1個月，2015年8月26日—9月27日。模型計算時間步長根據CFL條件進行動態調整，確保模型計算穩定進行，最小時間步長為0.05 s。底床糙率通過曼寧系數(M)進行控制，M取20～60 m1/3/s。夏季、冬季的溫度場均選擇連續半月潮進行計算。

根據2002年通過原國家環境保護總局審批的《廣東LNG接收站和輸氣干線項目一期工程環境影響報告書》中采納的標準，建議排放口冷排水溫差標準為≤5 ℃。結合排放口的溫差監控設置，要求在5 ℃以內警報，因此排水溫度根據取水口溫度減少5 ℃計算。按照設計文件，余氯排放濃度不大于0.2 mg/L，計算中以此作為計算源強，排放位置與冷排口位置相同。水面綜合散熱系數根據GB/T 50102—2014《工業循環水冷卻設計規范》計算，可得夏季的水面綜合散熱系數為41 W/(m2·℃)，冬季的水面綜合散熱系數為36 W/(m2·℃)。

2.3 邊界條件與初始條件

采用的數值模式中，水動力模塊需給定閉邊界和開邊界2種邊界條件。其中，開邊界選擇俄勒岡州立大學建立的TPXO全球潮波預報模型，其在中國海域附近的精度為1/30°×1/30°。用TPXO全球潮波預報模型提取8個主要分潮(M2、S2、N2、K2、K1、O1、P1、Q1)的調和常數，利用MIKE內置的潮位預報工具箱調和出水位邊界數據作為開邊界。選擇以大海域和工程周邊的岸線作為閉邊界。開邊界條件即水域邊界條件，本研究開邊界給定潮位，即η=η(x,y,t)；閉邊界條件即水陸交界條件，在該邊界上水質點的法向流速(Vn)為0。

計算開始時HD模塊“冷態”起動，設定初始條件為0，即所有網格點上的初始水位、流速均為0。采用穩態啟動，初值的誤差會隨著計算的進行逐步消失，模型驗證的數據取模型計算穩定后的模擬值。溫度場根據美國國家環境預報中心(NCEP)提供的溫度逐日資料，結合現狀監測資料給出。夏季因各月的溫度逐日變化不大，結合現狀資料取平均值，設為27.3 ℃。冬季溫度逐日變化可在1.0 ℃以上，故本底溫度場為計算多次后的穩定溫度場。

2.4 預測模式

2.4.1冷排水對水環境影響預測

采用水流模型與水溫擴散模式相結合的方法預測冷排水的溫降對水環境的影響范圍，計算公式如下：

(1)

(2)

式中：H為靜止時河道的水深，m；T為二維垂向平均溫度，℃；h為常規水位，m；u、v分別為流速笛卡爾坐標系中的分量；FT為水平溫度擴散項，g/s；Ts為冷排水源項，g/s；S為源項，g/s；Dh為溫度水平擴散系數，Dh=A/σT；σT為普朗特常數，取0.9；A為水平紊動黏性系數，取1.0。

2.4.2余氯對水環境影響預測

采用水流模型與耗散模式相結合的方法預測冷排水中余氯對水環境的影響范圍，耗散模式計算公式如下：

(3)

S=qC0

(4)

Q=KC

(5)

式中：C為水環境中余氯濃度，mg/L；Q為耗散項，g/s；q為排放量，m3；C0為排放余氯濃度，mg/L；K為衰減系數，K=ln 2/T1/2;T1/2為半衰期，取1 h。

2.5 模型驗證

選取有實測水文資料的P1、P2點作為潮位驗證點，C1、C2、C3、C4點為潮流驗證點(圖3)。利用2015年8月30—31日大潮期25 h同步連續觀測數據進一步驗證模型的模擬結果。潮位的驗證結果見圖4，P1、P2點的納什效率系數分別為0.90、0.87；潮流的驗證結果見圖5，C1、C2、C3、C4點的流速納什效率系數分別為0.75、0.63、0.65、0.81，流向納什效率系數分別為0.72、0.55、0.53、0.52。可見，此模型能夠比較準確地模擬湄洲灣東吳港區的潮流場情況。

圖3 驗證點位置Fig.3 Location of validation point

圖4 8月30日大潮期潮位驗證曲線Fig.4 Validation curve of tidal level during high tide period on August 30

圖5 8月30—31日大潮期潮流驗證曲線Fig.5 Validation curve of tidal current during high tide period from August 30 to 31

3 結果與討論

3.1 冷排水溫降對水環境影響范圍的預測

LNG接收站冷排水在海域內逐漸向外擴散，本次預測中，冷排水擴散方向主要是沿岸擴散。圖6為夏季和冬季一個潮周期內冷排水溫降包絡線與周邊保護目標的疊置圖。從圖6可以看出，不同溫降條件下的最大影響范圍不同，溫降條件越低，影響范圍越大。不同溫降條件下的最大影響范圍如表1所示。由表1可知，溫降≥4 ℃的等溫線在夏季擴散面積達0.012 km2，冬季為0.014 km2；溫降≥2 ℃的等溫線在夏季擴散面積達0.062 km2，冬季為0.068 km2；溫降≥1 ℃的等溫線在夏季擴散面積達0.142 km2，冬季為0.128 km2；溫降≥0.5 ℃的等溫線在夏季擴散面積達0.200 km2，冬季為0.202 km2。

圖6 夏、冬季冷排水溫降包絡線Fig.6 Envelope diagram of temperature drop of cold drainage in summer and winter

表1 不同溫降條件下的最大影響范圍Table 1 Maximum influence range under different temperature drop conditions

根據張繼周[4]的研究，水溫低于適宜溫度范圍，將會抑制魚類的新陳代謝和生長發育，如果超過其忍受限度，還會導致死亡。溫降達3 ℃以上時，蝦類幼體的生長可能會受到抑制，其存活率可能會降低。從圖6可以看出，在整個潮周期間內，溫降≥0.5 ℃條件下的影響范圍最大，此時距湄洲島生態系統重點保護區和養殖區分別約為2.81和0.53 km，不會對保護目標的水溫產生直接影響。

3.2 冷排水中余氯對水環境影響范圍的預測

根據耗散方程及源強，計算LNG接收站排水口冷排水中余氯的擴散范圍，統計1個月內余氯的擴散范圍，結果見圖7。從圖7可以看出，不同余氯濃度下的最大影響范圍不同，余氯濃度越低，影響范圍越大。不同季節余氯最大影響范圍如表2所示。由表2可知，夏季余氯濃度為0.02 mg/L的排水口包絡線面積為0.257 km2，而余氯濃度為0.01 mg/L的排水口包絡線面積為0.434 km2，此時排水口的影響范圍為0.628 km。故LNG接收站冷排水中的余氯不會對周邊敏感目標產生直接影響。

圖7 夏、冬季排水口余氯擴散包絡線Fig.7 Diffusion envelope diagram of residual chlorine in outfall in summer and winter

表2 不同季節余氯最大影響范圍Table 2 Maximum influence range of residual chlorine in different seasons

4 結論

(1)利用MIKE21，采用模型嵌套的方式精細化模擬了湄洲灣東吳港區LNG冷排水對水環境的影響，經驗證，潮位納什效率系數均達到0.8以上，流速、流向納什效率系數均達到0.5以上，模型能夠比較準確地模擬湄洲灣東吳港區的潮流場情況。

(2)在溫降≥0.5 ℃條件下，LNG接收站的冷排水影響范圍為0.202 km2，距離湄洲島生態系統重點保護區約2.81 km，距離最近養殖區約0.53 km,LNG接收站冷排水的溫降未對周邊敏感目標產生直接影響。

(3)在夏季余氯濃度為0.01 mg/L條件下，LNG接收站的冷排水中余氯影響范圍為0.434 km2，最大影響距離為0.628 km,LNG接收站冷排水的余氯未對周邊敏感目標產生直接影響。