梁賢金，田娟娟

(1.武漢凱迪電力工程有限公司，湖北武漢 430023；2.武漢天華嘉易建筑設(shè)計(jì)有限公司，湖北武漢 430070)

某濱海電站溫排水?dāng)?shù)值模擬研究

梁賢金1，田娟娟2

(1.武漢凱迪電力工程有限公司，湖北武漢 430023；2.武漢天華嘉易建筑設(shè)計(jì)有限公司，湖北武漢 430070)

濱海直流循環(huán)電站溫排水水量大，對(duì)鄰近海域的生態(tài)環(huán)境有重要影響。正確預(yù)測(cè)評(píng)估溫排水對(duì)鄰近海域的影響范圍，對(duì)于直流循環(huán)電站工程設(shè)計(jì)有重要意義。數(shù)值模擬是開展溫排水研究工作的重要手段。針對(duì)菲律賓某濱海電站循環(huán)水取排水工程的實(shí)際特點(diǎn)，建立水流數(shù)學(xué)模型和溫升場(chǎng)數(shù)學(xué)模型，對(duì)其溫排水?dāng)U散過(guò)程進(jìn)行數(shù)值模擬研究。對(duì)其河口海岸的復(fù)雜不規(guī)則邊界，離散方法采用非結(jié)構(gòu)網(wǎng)格。模擬了某濱海電站循環(huán)水取排水工程初步設(shè)計(jì)方案的溫排水?dāng)U散過(guò)程，模擬結(jié)果表明：設(shè)計(jì)方案下，向海中排放的溫排水不會(huì)對(duì)所在海域環(huán)境造成明顯影響，符合當(dāng)?shù)丨h(huán)保部門3℃溫升范圍要求。

濱海電站；溫排水；一次循環(huán)；數(shù)值模擬

隨著“一帶一路”戰(zhàn)略的深入實(shí)施，中資企業(yè)采用資本輸出的方式在沿線國(guó)家建設(shè)了許多大型濱海火電站。濱海電站運(yùn)行需要以冷卻水為載體將廢熱釋放到海洋中，大量的溫排水一方面影響了排水口區(qū)域的流場(chǎng)，另一方面會(huì)造成排水口區(qū)域水溫上升，對(duì)海洋生態(tài)環(huán)境產(chǎn)生不同程度的影響。目前我國(guó)對(duì)溫排水的認(rèn)識(shí)還存在分歧，對(duì)某些法律法規(guī)和標(biāo)準(zhǔn)中的有關(guān)條款的認(rèn)識(shí)也存在差異，還沒有專門的冷卻水排放標(biāo)準(zhǔn)，僅在一些水環(huán)境質(zhì)量標(biāo)準(zhǔn)中對(duì)水體的溫升提出了明確的規(guī)定，而一些海洋國(guó)家的標(biāo)準(zhǔn)中，已對(duì)溫排水有明確的要求。在此領(lǐng)域，國(guó)外一些學(xué)者采用CORMIX3模型[1]、神經(jīng)網(wǎng)絡(luò)算法[2]等對(duì)此有深入研究；國(guó)內(nèi)在這方面也有較多研究，應(yīng)用的模型有深度平均的二維淺水潮波模型[3]、平面二維溫排水?dāng)?shù)學(xué)模型等[4]，而大多研究都面向電站規(guī)劃選址階段，對(duì)電站運(yùn)行后溫排水溫升影響的數(shù)值研究相對(duì)較少[5]。

對(duì)于沿岸淺海水域，特別是半封閉海灣，其基本運(yùn)動(dòng)是由外來(lái)潮波引起的潮汐運(yùn)動(dòng)，即協(xié)振潮。表述潮波運(yùn)動(dòng)的參考坐標(biāo)系不考慮地球曲率的影響，顯然適用于水平范圍遠(yuǎn)小于地球半徑的海域，對(duì)于沿岸海域和海灣無(wú)疑是適用的[6]。本文采用不規(guī)則三角形網(wǎng)格的分布雜交方法建立水動(dòng)力學(xué)模型及二維熱輸運(yùn)數(shù)學(xué)模型，對(duì)菲律賓某海灣燃煤電站工程溫排水的擴(kuò)散過(guò)程進(jìn)行模擬，為其取排水工程初步設(shè)計(jì)方案的評(píng)估與驗(yàn)證提供依據(jù)，也為項(xiàng)目建成后的環(huán)境監(jiān)測(cè)提供必要的參照數(shù)據(jù)。

1 溫排水?dāng)?shù)學(xué)模型

1.1基本方程

針對(duì)某濱海電站所在海灣的地形、陸域邊界及潮流特征，選用沿水深平均的平面二維水流數(shù)學(xué)模型，其基本方程[7]為：

連續(xù)方程：

(1)

動(dòng)量方程：

(2)

(3)

式中，H、Z分別為水深和水位，m；u、v分別為x、y向的流速，m/s；u*、v*分別為源(匯)輸入(出)河道時(shí)x、y向的流速，m/s；q為源(匯)單位面積流量，源時(shí)q取正，匯時(shí)q取負(fù)，m3/(s·m2)；ρ為水體密度，kg/m3；γt為紊動(dòng)黏性系數(shù)，m2/s；c為謝才系數(shù)，c=1/n×R1/6；R為水力半徑，m；n為海床糙率；f=2ωsinφ，為柯氏力系數(shù)，ω為地球自轉(zhuǎn)角速度，φ為計(jì)算水域所在地理緯度。

所采用的溫升場(chǎng)數(shù)學(xué)模型基本方程為：

(4)

式中，ΔT為水體增溫量；Ex、Ey分別為x、y向擴(kuò)散系數(shù)；Ks為水面綜合散熱系數(shù)；ρCp為水體定壓比熱通量；ΔT*為沿程源匯水體超溫值。

1.2數(shù)值方法

水流模型基本方程中含有非線性混合算子，可采用剖開算子法進(jìn)行離散求解。這一數(shù)值方法可以根據(jù)方程所含算子的不同特性，將該方程剖分為幾個(gè)不同的子算子方程，各子算子方程可采用與之適應(yīng)的數(shù)值方法求解[8]。這種方法能有效地解決方程的非線性和自由表面確定問(wèn)題，具有良好的計(jì)算穩(wěn)定性和較高的計(jì)算精度[9-10]。

2 數(shù)值模擬計(jì)算

2.1取排水工程方案

某3×135 MW流化床燃煤電站機(jī)組位于菲律賓某海灣，主機(jī)冷卻水系統(tǒng)采用海水一次循環(huán)冷卻水系統(tǒng)，3臺(tái)機(jī)組冷卻水量約20 m3/s。當(dāng)?shù)丨h(huán)保部門要求，電站排出的溫排水產(chǎn)生的3 ℃溫升影響范圍不能超過(guò)半徑150 m圓形水域。

循環(huán)冷卻水取排水工程采取深取淺排布置方式，取水口布置在-7.0 m等深線附近，取水管線距岸邊約300 m，由3條引水鋼管組成，管徑DN2200；排水管線位于取水管線南側(cè)約250 m處，排水口布置在-4.0 m等深線附近，與取水口直線距離約285 m，排水管線為3根DN2200鋼管。取排水工程總布置如圖1所示，計(jì)算海域范圍及其地形變化如圖2所示。

圖1 取排水工程總平面布置Fig.1 Layout of water intake and drainage project

圖2 某電站計(jì)算海域范圍及其地形變化Fig.2 Computation of sea area range and topography change

2.2計(jì)算條件

水流溫升場(chǎng)數(shù)學(xué)模型中，主要包含河床糙率n、紊動(dòng)黏性系數(shù)γt、擴(kuò)散系數(shù)Ex和Ey、水面綜合散熱系數(shù)Ks共5個(gè)待定參數(shù)。根據(jù)電站所在海灣的地形條件，糙率采用0.016～0.020。紊動(dòng)黏性系數(shù)γt可采用零紊流模型計(jì)算：γt=αu*H，其中α為常數(shù)，u*為摩阻流速。擴(kuò)散系數(shù)可選用Ex=5.0 m2/s，Ey=0.5 m2/s。Ks根據(jù)《工業(yè)循環(huán)水冷卻設(shè)計(jì)規(guī)范》推薦的公式計(jì)算得出，約為34.0 W/(m2·℃)。

電站溫排水?dāng)?shù)值模擬計(jì)算選取計(jì)算范圍為廠址所在的海灣，整個(gè)計(jì)算區(qū)域南北長(zhǎng)約53 km，東西寬約35 km。計(jì)算區(qū)域采用三角形單元進(jìn)行網(wǎng)格剖分，共布置了9618個(gè)節(jié)點(diǎn)，18 687個(gè)單元，外海海域的網(wǎng)格間距約為1.0 km，由外海向廠址海域計(jì)算網(wǎng)格逐漸加密，在廠址近岸處網(wǎng)格間距約為100 m，在電站取排水工程附近對(duì)網(wǎng)格進(jìn)行進(jìn)一步加密，最小網(wǎng)格間距為10～20 m。電站溫排水所在海域計(jì)算網(wǎng)格剖分圖如圖3所示。

圖3 某電站附近海域計(jì)算網(wǎng)格剖分Fig.3 Meshes for the numerical model near the power plant

2.3結(jié)算結(jié)果及分析

2.3.1潮流場(chǎng)計(jì)算結(jié)果分析

電站附近海域潮流場(chǎng)分布直接影響到排水口溫升的輸運(yùn)擴(kuò)散，計(jì)算海域潮流場(chǎng)主要受外海潮汐的影響，計(jì)算水文條件下的潮流場(chǎng)是相似的。以典型大潮水文條件進(jìn)行研究。圖4為典型大潮水文條件下，計(jì)算海域典型瞬時(shí)潮流場(chǎng)分布。由圖4可知，計(jì)算海域海床較深、潮差較小，潮流動(dòng)力極弱，絕大部分海域流速小于10 cm/s，電站近岸海域水深較淺，流速相對(duì)稍大，但最大流速僅在20 cm/s左右。

電站廠址前沿流速較小，水流較為平順。在自然狀態(tài)下，全潮過(guò)程潮流均為由海灣南側(cè)流向海灣北側(cè)的近岸流，受河水流和電站排水的影響，在全潮的不同時(shí)段在排水管線的左右兩側(cè)會(huì)出現(xiàn)短時(shí)間的回流區(qū)。電站附近海域的低流速特性和回流區(qū)的存在必然會(huì)降低潮流對(duì)溫排水的輸移、擴(kuò)散效果。

電站附近海域深海與近岸海域的潮流變化不同步，深海的潮流變化滯后于近岸海域，同時(shí)潮位與潮流的變化也是不同步的，潮流的變化要滯后于潮位的變化。計(jì)算海域水下地形變幅較大，流速變化較明顯，外海海域水深超過(guò)1500 m，潮流流速最小，電站附近海域開闊，水深條件較好，潮流流速稍大，廠址南北兩側(cè)的淺灘受陸域邊界的影響，入海河流受上游來(lái)流和外海潮汐的影響，這些近岸水域的潮流流速最大。

圖4 電站近岸海域漲(落)急潮流場(chǎng)分布Fig.4 Instant flow velocity field near the power plant

2.3.2溫升場(chǎng)計(jì)算結(jié)果分析

采用沿水深平均的平面二維數(shù)學(xué)模型對(duì)電站3臺(tái)機(jī)組運(yùn)行后在3種典型水文條件下的溫排水影響進(jìn)行了預(yù)測(cè)計(jì)算，計(jì)算水文條件下的全潮最大溫升包絡(luò)線，計(jì)算結(jié)果如圖5所示，相應(yīng)的溫升分布統(tǒng)計(jì)如表1所示。

由圖5和表1可知，電站近岸海域潮流主流方向?yàn)橛赡舷虮钡某绷?，潮流?dòng)力微弱，全潮漲落流速不超過(guò)20 cm/s，且受河流入?？谒骱碗娬九潘绊懞?，在排水口附近水域有回流區(qū)存在。由于采取近岸直排方式，排水方向基本與潮流方向垂直，出口流速遠(yuǎn)大于海域流速，在漲潮時(shí)，溫排水被河口水流帶向南側(cè)海域，在落潮時(shí)，溫排水被輸運(yùn)至排水口北側(cè)海域。電站溫排水主要影響排水口南北側(cè)海域，且主要集中在近岸海域附近，呈扁長(zhǎng)狀分布。

計(jì)算水文條件下的溫排水影響面積差別不大，其中0.5 ℃和1.0 ℃等溫升分布范圍稍大，2.0 ℃及其以上高溫升分布范圍較小，僅局限于排水口附近水域。相對(duì)而言，典型大潮水文條件下溫排水影響范圍較大：0.5 ℃等溫升值全潮最大包絡(luò)面積1.246 km2，1.0 ℃等溫升值全潮最大包絡(luò)面積0.410 km2，2.0 ℃等溫升值全潮最大包絡(luò)面積0.140 km2，3.0 ℃等溫升值全潮最大包絡(luò)面積0.045 km2，4.0 ℃等溫升值全潮最大包絡(luò)面積0.016 km2，且3.0 ℃等溫升值包絡(luò)線最大長(zhǎng)度280 m，最大寬度230 m，均小于300 m，沒有超出半徑150 m的圓形水域。

圖5 全潮最大溫升包絡(luò)線圖Fig.5 Isothermal enveloping curves of maximum temperature rise

表1 3臺(tái)機(jī)組運(yùn)行溫升影響范圍特征值

2.3.3取水口溫升特征

排水口與取水口的相對(duì)位置關(guān)系對(duì)取水口溫升有直接影響，取水溫度又直接影響機(jī)組出力。模擬計(jì)算結(jié)果顯示，典型大潮水文條件下，電廠全潮最高取水溫升0.40 ℃、全潮平均取水溫升0.38 ℃；典型中潮水文條件下，電廠全潮最高取水溫升0.32 ℃、全潮平均取水溫升0.30 ℃；典型小潮水文條件下，電廠全潮最高取水溫升0.22 ℃、全潮平均取水溫升0.20 ℃。

綜合前文分析可知，電站附近海域漲落潮動(dòng)力微弱，機(jī)組取水主要受排水水流和河流入海水流的影響。取水溫升在全潮過(guò)程中有所變化，但取水水域潮流在全潮過(guò)程中均為由南向北流動(dòng)，因而溫升波動(dòng)幅度較小。取水溫升最大值為0.40 ℃，對(duì)機(jī)組運(yùn)行出力影響較小，滿足工程設(shè)計(jì)要求，即取、排水口布置位置關(guān)系合理。

3 結(jié)語(yǔ)

采用沿水深平均的平面二維水流溫升場(chǎng)數(shù)學(xué)模型對(duì)某濱海電站典型水文條件下溫排水分布規(guī)律進(jìn)行了計(jì)算分析研究，主要結(jié)論如下：

(1)從溫排水的角度考慮，電站近岸排水、近區(qū)深層取水的方案是可行的，滿足當(dāng)?shù)丨h(huán)保標(biāo)準(zhǔn)要求。

(2)電站所在海域潮流動(dòng)力較弱，計(jì)算水文條件下的溫排水影響面積不大，其中0.5 ℃和1.0 ℃溫升分布范圍稍大，2 ℃及以上高溫升分布范圍較小，僅局限于排水口附近水域。

(3)電站取水口最高溫升0.4 ℃，對(duì)機(jī)組運(yùn)行出力影響較小，符合工程設(shè)計(jì)要求。

從研究目的、效率等因素考慮，數(shù)值模擬采用了簡(jiǎn)化模型。但濱海電站所在海域水體密度差異顯著，水下地形及水動(dòng)力條件復(fù)雜多變，溫排水在其中的輸移擴(kuò)散具有明顯的三維特性。綜合考慮計(jì)入環(huán)境水體本底溫、鹽、流相互作用以及水氣間實(shí)時(shí)熱交換過(guò)程，才能更真實(shí)地反映溫排水在環(huán)境水體的擴(kuò)散過(guò)程。隨著模擬手段的進(jìn)一步發(fā)展，溫排水的數(shù)值模擬將會(huì)更多地考慮這類因素。

受基礎(chǔ)資料限制，數(shù)值模擬是在現(xiàn)有海床地形條件下進(jìn)行的，工程建成后可能會(huì)影響附近海床的沖淤變化，從而對(duì)電站取水產(chǎn)生一定影響，建議關(guān)注電站運(yùn)行后附近海域的地形變化趨勢(shì)，做到防患于未然。同時(shí)，開發(fā)應(yīng)用可真實(shí)反映溫排水對(duì)自然水體實(shí)時(shí)影響的溫排水預(yù)報(bào)模型，對(duì)于電站溫排水環(huán)境影響后評(píng)價(jià)和實(shí)時(shí)監(jiān)測(cè)具有重要意義，這也是今后提升溫排水?dāng)?shù)值預(yù)測(cè)技術(shù)水平的一個(gè)方向。

[1] Schreiner S P, Krebs T A, Strebel D E, et al. Testing the CORMIX model using thermal plume data from four Maryland power plants[J]. Environmental Modelling & Software, 2002, 17(3): 321- 331.

[2] Romero C E, Shan J. Development of an artificial neural network-based software for prediction of power plant canal water discharge temperature[J]. Expert Systems with Applications, 2005, 29(4):831- 838.

[3] 孫秀敏, 張令戈. 熱電廠溫排水排海環(huán)境影響預(yù)測(cè)方法及應(yīng)用[J]. 遼寧城鄉(xiāng)科技, 2001, 21(1): 30- 31.

[4] 張繼民. 電廠溫排水區(qū)流動(dòng)特性分析及模型參數(shù)的研究[J]. 東北水利水電, 2005, 23(253): 51- 56.

[5] 張貝貝, 周靜, 紀(jì)平. 濱海電廠溫排水?dāng)?shù)值模擬研究現(xiàn)狀[J]. 中國(guó)水利水電科學(xué)研究院學(xué)報(bào), 2014(4): 402- 409.

[6] 畢聞彬. 濱海電廠溫排水管理研究[D]. 青島: 中國(guó)海洋大學(xué), 2005: 32- 38.

[7] 張繼民. 電廠溫排水對(duì)水生生物的熱影響及水質(zhì)影響研究[D]. 南京: 河海大學(xué), 2006.

[8] 孫艷濤, 吳修鋒. 長(zhǎng)江水域電廠溫排水?dāng)?shù)值模擬[J]. 南水北調(diào)與水利科技, 2013, 11(6): 163- 166.

[9] 劉弘, 劉本芹. 南水北調(diào)東線分流區(qū)通航水流條件及改善措施[J]. 中國(guó)港灣建設(shè), 2008(5): 24- 27.

[10] 宣以飛. 廣州龍穴造船基地圍堤吹填工程設(shè)計(jì)方案研究[D]. 南京: 河海大學(xué), 2005.

Research on Numerical Simulation of Cooling Water Diffusion in Coastal Thermal Power Plant

LIANG Xian-jin1, TIAN Juan-juan2

(1.Wuhan Kaidi Electric Power Engineering Co., Ltd., Wuhan 430023, China; 2.Wuhan Tianhua Jiayi Architectural design institute Co., Ltd., Wuhan 430070, China)

It is important to predict and assess cooling water’s impact on the adjacent ecological environment in a one-through coastal thermal power plant with a large amount of cooling water. Numerical simulation is a significant means to carry out cooling water diffusion research. Considering the characteristics of a power plant’s circulating water system in the Philippines, models of flow and temperature rise are used to simulate the temperature and water drainage, and the unstructured grid is adopted for the discrete irregular boundary of the estuary and coast. The cooling water diffusion process of the circulating water intake and drainage preliminary scheme is simulated, and the simulation results show that the cooling water does not have obvious impact on the ecological environment of the sea area, which meets the requirements of the local environmental protection department for temperature rise range within 3 ℃.

coastal thermal power plant; cooling water; one-through circulating; numerical simulation

10.14068/j.ceia.2017.05.013

X57

： A

： 2095-6444(2017)05-0058-05

2017-01-20

梁賢金(1983—)，男，湖北宜昌人，高級(jí)工程師，碩士，主要從事火力發(fā)電廠水工工藝方面的設(shè)計(jì)管理及研究工作，E-mail：liangxianjin@kdpe.com.cn