國內外溫排水熱影響監測與評價研究方法探討

石繼香，官云蘭，趙利民，朱 利

(1.東華理工大學，江西 南昌330013;2.環境保護部衛星環境應用中心，北京100094; 3.中國科學院遙感與數字地球研究所，北京100094)

0 前言

核能技術是一項能夠提供大量能源而且不會釋放溫室效應氣體的能源技術，各國對于核電發展越來越重視，但是核能發展仍面臨著一些問題，其中溫排水問題已引起國際社會的普遍關注［1］。核電廠冷卻方式有直流冷卻、通風冷卻塔冷卻或是兩者結合使用進行冷卻，冷卻方式的不同，其造成的熱污染也不一樣。采用通風冷卻塔再循環冷卻，只需要考慮冷卻塔排污水的受納水體，熱影響較為微弱。由于條件限制采用冷卻塔冷卻的核電廠較少，多采用一次循環冷卻方式，將冷卻水直接排回自然水體，排水溫度一般高于受納水體環境溫度6～11℃，造成局部水域水體溫度的急劇升高，改變自然水體的水質，因此，保障海域環境和生態系統正常運行，加強對核電站溫排水熱污染分布范圍的探測以及對熱污染程度的定量評價，是核電站運行期間環境監測與評價的重要工作之一［2～5］。國內外的海洋環境學家們從20世紀70年代開始，就核電廠溫排水對附近海域的影響做過很多的調查和研究工作。

1 研究方法

目前，國際上用于核電廠溫排水的環境影響評價的主要技術手段有實地測量、遙感觀測和數學物理模擬。隨著微尺度(1 km)溫排水熱污染研究的需要，X波段航海雷達技術也逐漸被用于核電站的熱污染監測。一般情況下，較少單一使用某一種方法來進行核電站溫排水熱污染影響研究，多使用幾種方法相結合的方法，相互輔助，相互驗證來進行研究。

1.1 實地測量與遙感相結合方法

實地測量精度較高，但成本也高、時效性較差，難以反映較大范圍的空間分布，很少單獨采用此方法進行核電站熱污染監測。衛星遙感的手段，大面積同步觀測可以快速、及時、直觀地得到監測結果。2種方法相結合，可以開展衛星調查與地面調查的比對，地面調查數據可以驗證遙感數據溫度反演的精度。

孫戀君［6］等對田灣核電站排放口附近海域溫排水環境影響范圍和程度進行調查，就采用了以航天遙感測量為主、航空遙感測量為輔、近海面實地測量為補充的方式，獲得目前田灣核電站排放口附近海域海面溫度場特征分布。調查結果顯示，一期工程投入運行后，排放口附近海域環境有明顯影響，且隨著核電機組的擴建有繼續擴大的趨勢。周穎［7，8］等利用減災環境衛星B星(HJ-1B)和風云三號衛星(FY-3)的熱紅外數據，使用單通道溫度反演算法，反演田灣核電廠附近海域海面溫度，反演結果與MODIS海表溫度產品具有較好的一致性，說明HJ-1B和FY-3的熱紅外數據能夠進行溫排水的分布研究。進而分析了溫排水的時空分布變化規律，發現夏季溫升超過3℃的區域遠小于冬季，夏季溫排水沿海岸單向狹長擴散，冬季的溫升區域呈扇形分布;低潮時超過3℃的溫升區域面積是高潮時的4倍。X Y Dai［9］等基于Landsat TM和ETM+6波段熱紅外圖像量化田灣核電站周圍海水表面溫度(SST)變化的一般規律，及2003－2009年田灣核電站熱污染的影響程度，并進行了熱污染影響的強度和范圍的調查，以及不同的季節之間的熱污染面積、強度對比，調查得出這期間的溫升區主要分布在核電站附近的淺海水域，超出基準溫度最高達6.2℃，也得出冬季溫排水熱污染強度及范圍明顯大于夏季的結論。

1.2 物理模型模擬方法

數學物理模擬主要是用于核電廠發電前評價的主要技術手段，模型的選擇、水面綜合散熱系數、擴散系數等典型參數的不同選擇均影響了評價的準確度［10］，物理模型能較好地模擬排水口近區三維水流、熱水摻混和浮力效應。國內大多學者都是采用物理模型實驗來進行熱擴散近區的研究，利用物理模型實驗來模擬溫排水近區的溫升分布。辛殿文［11］等基于2005年枯水期施測的大、小潮水文資料及數值計算的成果，建立了澳門路環發電廠溫排水物理模型，模型潮流驗證試驗表明，各驗證點位的潮位、流速、流向驗證曲線均吻合較好，模型潮流模擬達到了較高的精度。模型除了需要對潮流驗證外，還對電廠現狀溫排水的影響進行驗證。詳細研究了澳門機場擴建及機場以西水域規劃實施后電廠溫排水造成的熱污染問題，通過改變不同排水口位置，對排水口附近水域表面溫升以及取水口的取水溫升分布進行對比分析，明確了工程布置的最優方案。

1.3 值模擬方法

1.3.1 二維數值模擬方法 數值模擬又分為二維和三維模型，淺水二維模型能較好地模擬大范圍水域潮流場，較好地反映溫排水的對流、擴散和累積等效果，并模擬計算水面熱量散失的影響。

J Yu［12］等采用分步有限元方法與三角形網格結合二維流體力學模型很好的模擬了日照海域的特點，在此基礎上進行日照電廠周圍海表的溫度熱傳輸和水平分布的數值模擬，并考慮排水口的位置、溫排水排放量等因素，設計不同方案預測溫排水對水環境的影響，得出從一期和二期項目P1排水口進行海水冷卻方案最佳。張舒羽［13］等采用MIKE21軟件對核電廠溫排水進行數值模擬，分析不同潮型、不同流量及不同季節下的最大溫升包絡面積和取水口的溫升變化。發現溫升包絡面積與電廠冷卻水的排放量不為線性關系，此研究成果可為核電廠的建設提供依據。孫艷濤［14］等采用Delft 3D軟件平面二維水流溫度場數學模型對常熟電廠一期、二期擴建工程溫排水進行數值模擬計算，模型計算的結果與實測值吻合良好，在此基礎上分析溫度場的分布情況。分析結果顯示，冬季小潮比大潮條件下的總的溫升影響面積要大，且向下游和離岸方向擴散的距離也更遠，取水口處的溫升面積均較小。

1.3.2 三維數值模擬相方法 二維海洋流體動力學模型無法知曉水平流速沿深度的變化、無法模擬垂向分布等，無法模擬溫排水近區廢熱水與周圍環境水體的卷吸稀釋混合現象，因此二維模型存在某些的局限性。為了精確地模擬熱擴散近區的溫升分布，必須進行三維熱擴散預測方法的研究［15］。三維模型研究多采用物理模型試驗或試驗和數值模擬相結合的方法。但由于受到實驗室條件限制，模型范圍相對較小，模型模擬的水域范圍較小，不易準確模擬環境水體蓄熱，熱水回歸影響，也很難精確調節模型中水面綜合散熱系數使其滿足散熱系數相似要求，因此會對結果產生一定影響，并且物理模型實驗實際工作中費用昂貴。隨著工程要求的提高，三維水流模型逐漸成為研究熱點，三維熱擴散的數值計算在近些年得到了眾多學者的重視。

我國目前采用三維數值模型進行研究已有一些進展。汪一航［16］等采用POM模式，建立了莊河電廠附近海域的包括4個主要分潮的潮汐潮流三維數值模擬，然后用所得流場對莊河電廠的溫排水問題進行了數值模擬。彭溢［17］等以武漢陽邏電廠為例，采用三維水動力與溫度輸運數學模型預測溫升場分布。采用電廠溫排水平面與垂向兩個方向上的實測溫升資料對三維溫升場數學模型進行驗證，模擬結果較好的反映了水域中熱量的輸移擴散規律。模擬結果顯示，溫排水排入受納水域后，由于密度較小，水體表面的擴散范圍要大于底層，從垂向分布來看，溫升分層現象明顯，溫水的擴散主要在中上層進行，距離排放口越遠，溫升影響越小。

國外大多數學者則是采用水利學實驗和數值實驗相結合的方法進行熱擴散的研究。D J Lowe［18］等專門對較大溫排量的情況進行了研究，功率為877 MW的CharlesPoletti核電廠和功率為1 253 MW的Astoria燃油電廠相距僅500 m遠，這導致了在研究區域內水流較大流速較高。在低海潮的情形下，由于較低的水流速度情況下將會有更少的水可以使用來稀釋溫排水，所以電廠溫排水處于水流最壞的情形下，D J Lowe基于非結構化三角形元素采用靈活的網格配方與三軸溫度場實地調查數據成功校準后，采用三維水動力模型MIKE3，對排放到紐約市皇后區東河冷卻水熱羽流的影響進行評估，表明在溫排量最大情況下產生的熱羽流也不超過熱污染監管標準。

三維熱擴散數值模型方法的優劣很大程度上在于湍流的參數化方法上，為了盡可能地消除系數經驗選取的人為性，常應用湍流模式理論。常用的有雷諾應力模型和κ－ε雙方程模型。雷諾應力方程模式是比較精細的模型，能正確的反映浮力的作用，模擬效果較好，具有較高的精度，是目前比較精確的數學模型。但其計算繁瑣，因此在復雜工程中的應用受到一定限制。H Ramirez-Leon［19］等從海逸豪園、韋拉克魯斯、墨西哥核電廠，考慮鑒于上游的自由表面與大氣相互作用的熱通量，使用CAL-數值模型解決用于計算淺水區的溫度變化的納維-斯托克斯雷諾應力方程和能量方程，來進行沿海核電廠溫排水的研究，并能量方程計算水深測量，海洋，氣象，水文和運行機組數目用于數值模擬，分析熱羽流擴散規律;并對旱季和雨季的溫排水溫升結果進行定性統計，通過模擬結果與觀測數據比對，顯示出較好的一致性。

近岸海域海水溫度變化除了流速、濃度、紊動特性、潮汐因素外，還受淺水分潮和近岸地形等因素的影響，因而近岸海水溫度的預測也是熱擴散預測方法研究中的一個重要課題。國內外學者有的采用修正的 k-ε模型和擴展的 k-ε模型。Z Li［20］等和Y L Cheng［21］等分別通過二維數值和三維數值模型，選擇k-ε湍流模型，模擬相關區域的水流速度，獲得不同水流速度和不同溫排水排水口角度條件下的溫度分布規律，并研究地形對電廠溫排水擴散的影響結果表明，兩人的研究表明，合理的選擇熱排水口角度和溫排水的流速可以減少溫排水對電廠排放口的影響。由于溫排水受到逆流的影響，一些溫排水又回到排放口附近的區域。由于再循環的影響，減小了溫排水的影響。斜率較大的地形，溫排水在水平方向上擴散較大，而在垂直方向上產生的影響更大。

為了評估溫排水污染的影響，通常核電站周圍的水域被分為兩類，鄰近區域和遠區域。在Jiang［22］等人的研究中，新開發出一個嵌入式三維網格模型，通過該模型可以直接對近區域和遠區域整個區域進行模擬，該模型已被用于韓國KIM和SEO核電站對平均深度數學模型的模擬結果進行評估。S Seddigh-Marvasti［23］等通過 NASIR軟件，對于來自來源點的底部溫排水建立模型并應用松散的三角形網格劃分。應用平均深度有限元體積流量求解器，模擬溫排水的排放，通過對比數值結果和現有的實驗室測量的表面溫度，求得平均深度和表面溫度之間的換算因子關系(深度平均溫度轉換到表面溫度)。為了解底部的熱污染射流對環境的影響和可能的自然環流的影響提供了更好的方法。

1.4 航海雷達方法

X波段航海雷達，這通常是用于檢測海岸線和海面上的障礙，是目前海洋遙感最常用的工具之一。基于海面風浪成像的原理，證明海面觀測的可行性。X波段航海雷達圖像通過不同的算法可以獲得海面風浪、潮流和水深等不同物理特征的重要信息，例如內波、碎波和河流羽狀河口區。除上述信息外同樣能夠觀測到溫排水分布情況。與機載和衛星的技術相比，基于地面的遙感監測主要是低掠角模式。基于地面的遙感監測也具有持續性和長期監測的優點。然而，使用低掠角海上雷達方式微掃描檢測溫排水熱污染特性一直以來雖很少受到關注，卻也有一些研究成果。

S V Nghiemet［24］等人研究墨西哥灣暖流環指出，通過KU波段散射，使用大的入射角(入射角大于20°和小于60°)的反向散射數據，可以檢測出各種風向下的冷風和暖風海洋表面溫升界線。并指出海洋表面溫度的變化引起海表空間風壓的變化，導致大氣層變化。因此，變化無常的壓力導致海浪的變化是中尺度峰面邊界影像中的主導因素之一。暖流環的面積比來自核電廠的溫排水熱污染要大得多。H Y Bau［25］研究指出雷達反向散射受海面風浪所影響，海表溫度的變化將會影響海面風浪，采用低掠角的航海雷達方法，微尺度檢測(不到1 km)溫排水，獲取平均連續的圖像序列，檢測核電站的排水口得到到一個明顯的回波帶圖案。這個區域的海表溫度數據回波帶與溫排水的空間分布形態相匹配。基于X-波頻段船用雷達獲得的海表圖像分析溫排水分布規律，確認溫排水熱污染的分布受海潮現狀所決定。

2 結論與建議

核電站溫排水對環境的影響較大，對海洋生物的生長和繁殖造成影響，對它進行動態監測是非常必要的，遙感技術具有成本低、速度快、資料同步性好、大范圍面狀觀測等優點，采用衛星遙感與大面測量相結合的方式對溫排水進行監測能夠反映核電廠運行期間溫排水的影響范圍。采用以航天遙感測量為主、航空遙感測量為輔、近海面實地測量為補充的方式進行核電站熱污染環境監測，并與數模模擬結果進行比對，將是未來溫排水物理影響監測與評價的方向。且進行足夠的不同時空和不同光譜分辨率的核電站溫排水熱污染信息采集，對海表溫升的時空格局變化進行量化分析，也將是未來溫排水熱污染監測的一個側重點。

核電廠廠址的合理選擇、取、排水口位置的合理布局以及核電站運行后溫排水熱影響的預測和評價工作，主要依賴于三維數值模型的模擬，而在國內將三維流體動力學數值模型用于熱擴散的研究中還較少，仍屬于探索的階段，因此三維熱擴散的數值計算在國內急需解決，仍需加強研究。

核電廠的熱污染是很常見的，通過余熱利用是減輕溫排水熱影響的根本性措施，在核電站溫排水物理影響監測與評價中，也應加強對余熱利用方式的研究和探索，在溫排水排入受納水體前，盡可能取出和利用部分余熱，使溫排水降低溫度后再排入收納水體，減小核電站溫排水對受納水體的熱影響。

［1］ 朱曉翔，劉建琳，王鳳英.核電站溫排水環境影響研究方法調查評價［J］.電力科技與環保，2010，26 (1):8－10.

［2］ 李靜晶，張永興，姜子英，等.濱海核電廠熱污染評價中幾個問題的探討［J］.核安全，2012，(1):30－34.

［3］ 朱 利，吳傳慶，等.基于環境一號紅外相機的田灣核電站溫排水遙感監測研究［J］.中國環境科學，2012，1(1):1－6.

［4］ 曹宇峰，曾松福，林春梅.淺談我國濱海電廠溫排水對海洋環境的影響狀況［J］.海洋開發與管理，2013:(2):72－75.

［5］ 陳子斌.國外內陸核電廠情況整理與分析［J］.能源研究與管理，2013，(1):15－21.

［6］ 孫戀君，王鳳英，朱曉翔.田灣核電站溫排水環境影響遙感調查［J］.中國輻射衛生，2011:20(3):330－332.

［7］ Zhou Y，Gong C L，Hu Y，et al.Distribution of thermal discharge from a power station based on HJ-1B and FY-3 thermal infrared data［C］.International Symposium on Photoelectronic Detection and Imaging，2011: Advances in Infrared Imaging and Applications，2011 (8193).

［8］ 周 穎，鞏彩蘭，匡定波，等.基于環境減災衛星熱紅外波段數據研究核電廠溫排水分布［J］.紅外與毫米波學報，2012，31(6):544－549.

［9］ Dai X Y，Guo Z Y，Lin Y，et al.An of Satellite Remote Sensing Data for Monitoring Thermal Discharge Pollution from Tian wan Nuclear Power Plant in Eastern China［C］.5th International Congress on Image and Signal Processing(Cisp):2012，1019－1023.

［10］張愛玲，汪 萍，陳瑩瑩.濱海核電廠溫排水物理影響模擬中幾個關注問題及審評見解①［J］.輻射防護通訊，2011，31(16):21－27.

［11］辛殿文，董志慧.澳門路環發電廠溫排水物理模型試驗研究［J］.人民長江，2012，43(20):62－65.

［12］Yu J，Liu Z L，Li Y L.Optimization of thermal discharge scheme for the phase II project of rizhao power plant［C］.Asia-Pacific Power and Energy Engineering Conference(Appeec)，2019，7(1):2177－2180.

［13］張舒羽，黃世昌，韓海騫.浙江蒼南電廠冷卻水溫排放的數值模擬［J］.海洋學研究，2009，27(3):61－66.

［14］孫艷濤，吳修鋒.長江水域電廠溫排水數值模擬［J］.南北水調與水利科技，2013，11(6):15－18.

［15］蔣 爽，端木琳，王樹剛.海水熱擴散研究進展與新問題分析［J］.能源與環境保護，2006，20(5):5－9.

［16］汪一航，魏澤勛，王永剛，等.潮汐潮流三維數值模擬在莊河電廠溫排水問題中的應用［J］.海洋通報，2006，25(1):8－15.

［17］彭 溢，張 舒.電廠溫排水三維模擬及特征分析［J］.電力科技與環保，2013，29(1):7－9.

［18］ Lowe S A.Schuepfer F，Dunning D J.Case study: three-dimensional hydrodynamic model of a power plant thermal discharge［J］.Journal of Hydraulic Engineering-Asce，2009，135(4):247－256.

［19］Ramirez-Leon H，Couder-Castaneda C，Herrera-Diaz I E，et al.Numerical modeling of the thermal discharge of the Laguna Verde power station［J］.Revista International De MetodosNumericos Para Calculo Y Diseno En Ingenieria，2013，29(2):114－121.

［20］Li Z，Lin D M，Shu H W，et al.Numerical simulation of thermal discharge for large-scale seawater source heat pump system［J］.Proceedings of the 5 th International Symposium on Heating，Ventilating and Air Conditioning，Vols I and Ii，2007:757－762.

［21］Cheng Y L，Hao Q Z，Hu J.Effects of Topography on Diffusion of Thermal Discharge in Power Plant［C］.2nd International Conference on Challenges in Environmental Science and Computer Engineering(Cesce 2011)2011，(11):618－623.

［22］Jiang J，Fissel D B，Topham D.3D numerical modeling of circulations associated with a submerged buoyant jet in a shallow coastal environment［J］.Estuarine，Coastal and Shelf Science，2003，58:475－486.

［23］Seddigh-Marvasti S，Sabbagh-Yazdi S R，Mastorakis N E.Extracting surface temperatures from depth average numerical solution near thermal water discharge from bottom diffusers［J］.Continuum Mechanics，Fluids，Heat，2010，225－230.

［24］Nghiem S V，Li F K，Walsh E J，et al.Radar backscat－ter across the Gulf Stream sea surface temperature front［J］.IEEE Trans.Geosci.Remote Sens，2000，38 (2):926－941.

［25］Bau H Y，Wu L C，Chuang L Z H.Detection of thermal effluent discharge using sea surface nautical radar images［J］.Ieee Geoscience and Remote Sensing Letters，2013，10(3):612－616.