湖南桃花江核電廠取排水工程對通航影響分析

周美林 ，夏 波 ，羅 忠 ，陳 杰

（1.湖南女子學院，長沙410072；2.長沙理工大學水利工程學院，長沙410004；3.湖南省水沙科學與水災害防治重點實驗室，長沙410004；4.海軍南海工程設計院，湛江524003）

湖南省桃花江核電工程位于益陽市桃江縣沾溪鎮肖車塘村謝家塘沖，毗鄰湖南省洞庭湖水系資水下游干流右岸，所在航段為Ⅴ級航道，雙線通航300 t級船型。核電廠擬裝4臺1 000 MW機組，冷卻水采用二次循環，取水流量為6 m3/s，擬從資水干流取水。取排水工程所在河段河谷呈開闊的U形，河道呈彎曲的S形延伸，河道中間為江心洲和趙林洲，趙林洲洲灘高程較高，其尾部左側還有一個小洲，高程略低（圖1）。取水工程位于右汊右岸，排水口位于取水建筑下游約1 km的右岸，排水渠道寬20 m，排水斷面平均流速為0.15 m/s。取排水期間，局部水域將會形成較大流速，并在一定范圍內形成水面坡降，工程引起的流速流態變化將對船舶的航行造成一定的影響。取排水工程上游約117.5 km處修建有柘溪水庫，下游約9 km處修建有修山水庫，工程位于修山水庫庫區回水末端，取排水會受到柘溪以及修山兩大水庫調度運行的影響。因此有必要對取排水工程進行研究和分析，為工程設計和航運安全提供科學依據。

在研究取排水工程對河道通航影響時，一般可采用物理模型或數學模型進行模擬研究，也可以同時使用2種模型配合研究。對此許多學者進行了大量的研究工作，盧祥興［1］等利用平面二維數學模型結合物理模型，對秦山核電三期工程取水口修建后附近海床的沖淤變化進行了預測；趙洪波［2］等通過對惠安核電廠取排水口工程海區水文、泥沙及歷年地形資料的分析，對該海區的泥沙來源、運移規律以及海床的演變特征和發展趨勢進行了研究；周美林［3］等從平面二維水動力學方程組出發，建立丁壩繞流平面二維水流數學模型，對河龍高速清溪河段水毀防護工程的丁壩群護岸進行了模擬分析；曹繼文［4］等采用雷諾平均的N-S方程，以各向異性的雷諾應力模型（RSM）封閉方程，對橫向取水口水流的復雜三維流動特性進行了研究；彭曉慶［5］等在廣州大學城雜用水廠取水工程興建后，采用數學模型計算了取水對取水口附近河道航深、流速和流態等造成的影響；高峰［6］等利用二維水流數學模型，采用嵌套網格的差分格式對華能營口熱電廠取水工程進行了計算分析；周美林［7］等基于有限單元法，建立了平面二維水流數學模型，針對取水口邊界條件進行了研究。類似還有謝翠松［8］和闕志夏［9］等人的工作。可見二維平面數學模型已經廣泛的應用于取排水工程對河道通航影響的研究中。本文基于有限單元法，建立平面二維水流數學模型，對湖南省桃花江核電廠取排水工程對通航的影響進行了研究。

1 數學模型的建立和驗證

1.1 基本控制方程

考慮到取排水工程所在河道段為寬淺水域，水力參數在垂直方向的變化明顯小于水平方向的變化，工程需要了解和計算的主要是工程前后物理量沿平面變化情況，因此本研究采用沿水深平均的封閉淺水方程組描述二維水流運動，二維淺水基本控制方程如下

式中：x、y為水平坐標軸；u、v為x、y軸向流速；t為時間變量；g為重力加速度；h為水深；α0為河床底面高程；ρ 為水流密度；f為柯氏力參數（f=2sinφ為地球旋轉角速度，φ 為緯度）；εxx、εxy、εyx、εyy為紊動粘滯系數；c為謝才系數。

1.2 數值離散

方程的離散包括時間離散和空間離散，時間的離散采用差分法，空間的離散采用有限單元法。運用Galerkin加權余量法把淺水方程離散成非線性代數方程，然后采用Newton-Raphson方法求解。離散區域內采用三角形六節點等參單元和四邊形八節點等參單元相耦合。根據對二次與線性或三次與二次函數混合插值，與二次與二次函數插值相比較，使用混合插值法雖然會輕微影響流速場計算的精度，但可以提高壓力場（水位）計算精度，并且可以減少節點數量，提高計算速度。因此在求解時采用混合插值法。

1.3 邊界條件

（1）上、下游邊界條件。上游邊界條件一般取研究區域上游來水流量的大小和方向，下游邊界一般給出水位過程線即可。

（2）閉邊界條件。模型中閉邊界采用滑移邊界條件，不考慮滲透，取法向流速為0。考慮可通過加大邊壁附近單元粗糙系數的方法，來模擬河道邊壁對附近水體流動的影響。

（3）動邊界處理技術。當研究河段有洲或灘存在時，隨著水位的變化，需要采用動邊界處理技術。本模型采用干濕法對動邊界進行處理，當單元水深不低于給定最小水深，整個單元為濕單元而參與計算；否則，整個單元為干單元而不參與計算，按線性插值調整參與計算的單元個數。

（4）取水建筑邊界條件。取水建筑附近的水流呈強三維運動，在只需要研究取水工程對航道通航的影響，而不需要了解取水建筑附近的三維流結構時，取水工程邊界可取為出流邊界，一般流量大小取為取水流量，也可設置為匯項，本模型中將取水工程邊界設置為出流邊界，可以更好地反映取水時的局部流態和實際情況。

1.4 模型驗證

根據周美林［3，7］等運用此模型對丁壩群和取水工程進行數值模擬的結果，可知該數學模型具有很高的精度，能滿足本文的研究需要。

2 水文分析計算

2.1 設計通航水位計算

根據《內河航道與港口水文規范》的規定，對于Ⅴ級航道，樞紐上游河段設計最高通航水位采用10 a一遇的洪水敞泄時的水面線，以及壩前正常蓄水位與相應的入庫流量組合的水面曲線，取其上包線作為沿程各站點的設計最高通航水位；樞紐上游河段設計最低通航水位采用保證率為95%的入庫流量與壩前運行水位相組合，以及壩前最低運行水位與相應的流量組合，得出多組回水曲線，取其下包線作為沿程各站點的設計最低通航水位。根據修山電站調度方式及工程所在河段的水文關系，通過分析計算，確定取水泵房設計最高通航水位為46.23 m，設計最低通航水位為42.4 m。

2.2 水位流量關系計算

取排水工程均位于已建成的修山水庫庫區內，其中取水工程位于分汊河道的右汊右岸，距修山壩址約9 km。因核電專用水文站設立時間不長，暫無汊道分流比觀測資料，只能按全斷面確定水位流量關系。水位流量關系按修山建庫前、后2種情況考慮，修山水庫建成前為天然狀態水位流量關系，建成后為受回水頂托的水位流量關系。修山及荷葉山廠址天然水位流量關系見表2。由于修山水庫蓄水，取水口斷面天然水位流量關系受庫區回水影響，不同的壩前水位對應的取水口斷面水位流量關系不同，采用推回水的方法計算，以修山水庫壩前水位作控制，計算出各流量對應的工程處水位，推回水所采用的斷面為2005年3月實測大斷面，河段糙率采用修山水庫設計時所采用的糙率（圖 2）。

2.3 研究工況確定

根據《內河航道與港口水文規范》和《內河通航標準》中關于通航水位及通航條件的規定，結合修山電站影響下工程位置的水位流量關系曲線結果，確定本研究的水位-流量組合（表1）。

表1 工程河段通航條件計算水位-流量組合Tab.1 Stage-discharge relation of navigation condition

3 通航計算分析

3.1 計算區域及網格設置

根據計算需要，選擇取排水工程上下游約6 km的河段作為計算區域。計算網格劃分采用三角形和四邊形混合的非結構化網格，分別為6節點單元和8節點單元，保證對流速計算具有二階精度。模型中布置了6 550個網格，共計17 835個格點，最小網格約5 m。

為分析取排水建筑附近流場，在取水口及排水口位置各設置5個斷面（圖3），布置測流點測量工程前縱向及橫向流速。

3.2 工程前通航條件計算分析

運用數學模型，按照表1水位-流量組合關系對工程前通航條件進行數值計算，分析結果得出工程前取水口位置存在橫向流速，在各級水位下的橫向流速均不大，最大值為0.24 m/s。排水口斷面亦存在橫向流速，最大值0.29 m/s。取、排水口縱向流速因水位高低而異，水位越高，縱向流速越大，當流量為8 950 m3/s時（10 a一遇洪水流量），取水口斷面河中縱向流速為2.65 m/s（圖4），在這種流量下，上水船舶可沿近岸緩流區行駛，因此當河段出現10 a一遇洪水時，船舶仍可正常通航，說明桃花江核電廠取水口河段在天然狀態下具有較好的通航條件。

3.3 修山電站建成后取排水河段的通航條件分析

修山電站建成后，壩位以上河段變為庫區航道，在電站死水位42.4 m時，其回水與上游在建的白竹洲電站下游最低通航水位銜接。庫區死水位時，取水口位置滿足1.6 m水深的寬度約為170 m，泵房所在位置的水深為1.4～6.5 m，泵房前沿航道內水深可達8.0 m，航道尺度較工程前大幅度提高，河段具有較好的通航條件。

由于工程河段為庫區航道，右汊通航條件得到徹底改善，左汊航道也同時得以改善，工程前左汊在枯水期幾乎斷流，工程后左汊最小水深可達2.4 m，航寬約130 m，能滿足Ⅴ級航道尺度要求，但該汊出口條件較差，彎曲半徑較小，而右汊是歷史上的通航汊道，其通航條件優于左汊，因此仍應以右汊作為通航汊道。

3.4 工程后通航條件計算分析

按照表1水位-流量組合關系，運用數學模型對工程后通航條件進行數值計算，分析結果表明，在修山電站正常擋水位43.0 m時，取、排水位置縱向流速很小，均小于0.1 m/s，橫向流速在0.02 m/s以下。當取水口水位為44.73 m，流量為6 000 m3/s時，最大橫向流速為0.42 m/s，該流速發生在2-2斷面12號測點上，距取水泵房上游角點前沿線約22 m，該測點前沿11號測點的橫向流速為0.16 m/s，距取水泵房外邊線約40 m。在流量為7 700～8 950 m3/s時（10 a一遇的洪水通航流量），12號測點的橫向流速均在0.41～0.42 m/s。與取水泵房鄰近的下游18號測點（橫向流速為0.33～0.34 m/s），形成泵房前沿橫向流速超標帶，過泵房后的24號測點橫向流速較小，為0.23 m/s。這種橫向超標流速的產生，除了由于水流在上游彎道轉折處產生橫流外，其主要原因是泵房阻水產生挑流作用引起橫向流速增加。洪水期泵房前沿縱向流速有所增加，主流上的22號測點工程前的縱向流速為2.65 m/s，工程后縱向流速增至3.06 m/s，均由泵房束窄斷面過水面積所引起，但這種較大縱向流速出現在河道中央，因此，建議上水船舶應提前避開主流急流區（22號測點），靠近趙林洲岸邊行駛，下水船舶進入右汊后，需要根據助航標志或警示牌提示，及早調整航向，靠趙林洲一側行駛，避免進入橫向流速超標區或撞擊泵房。當該河段在泄洪期間流量達到7 000 m3/s以上時，建議對右汊實行封航，船舶改為左汊航行。

排水口附近橫向流速較小，10 a一遇的洪水位時，最大橫向流速為0.28 m/s，基本不影響航行，最大縱向流速為2.78 m/s，出現在河道中部位置。建議上水船舶應避開主流航行。

從數值模型的計算結果看，取水口處最大橫向流速的位置距右岸岸線約50 m，距取水泵房順水流方向的外測邊線約為22 m，距左岸航道邊線約140 m，取水泵房為近岸式取水，進水閘與取水泵房合建，使用的通航水域范圍不大，因此取水泵房修建后，工程河段仍具有較好的通航條件。

工程對流速流態有一定的影響，一方面取排水改變局部河段的流速流向分布；另一方面建筑物的存在形成繞流及回流結構，將改變附近河段的流場流態。在航道可通航水文情況下，本研究對工程前后取排水工程附近樣點的最大縱向流速和最大橫向流速進行了統計分析，統計結果表明，在靠近取水建筑的河段，縱向流速和橫向流速的量值和方向都發生改變，從最大縱向流速和最大橫向流速的統計結果來看，在橫向距離取水建筑50 m以外，工程后的最大橫向流速滿足通航要求。主流中縱向流速較大，船舶可避開最大縱向流速區，沿趙林洲一側行駛。

4 結論

本文基于有限元方法，建立了平面二維水流數學模型，對湖南省桃花江核電廠取排水工程對通航的影響進行了數值研究。利用模型計算了工程前后取排水工程附近水流特性，分析了16種水位流量組合下取水口取水和排水口排水對航道通航條件的影響，得到如下結論：

（1）取水建筑物的存在將改變工程所在河段取水口附近水域的流場流態，洪水期會引起工程附近河段縱向流速和橫向流速的加大，但通過規劃航道路線、設置助航標志、加強水上交通安全管理后，船舶沿趙林洲一側航行，對船舶的正常通航影響較小。

（2）當修山電站在洪水期開閘泄洪，工程所在河段流量大于等于7 000 m3/s時，為保證通航安全，建議封閉右汊航道，船舶改由左汊航行。

（3）排水口附近橫向流速較小，10 a一遇的洪水位時，最大縱向流速2.78 m/s，出現在河道中部位置，上水船舶可避開主流靠右側航行。

［1］盧祥興，史英標，王敏.秦山核電三期工程取水口沖淤變化研究和預測［J］.泥沙研究，2000（6）：10-15.LU X X，SHI Y B，WANG M.Prediction about Seabed Erosion-Siltation near the Water Intake of Third-phase Project for Qinshan Nuclear Power Plant［J］.Journal of Sediment Research，2000（6）：10-15.

［2］趙洪波，吳以喜，楊華.惠安核電廠取排水口工程海區沖淤演變分析［J］.水道港口，2007，25（3）：159-162.ZHAO H B，WU Y X，YANG H.Analysis on Scour and Silting of the Sea Area near Intake and Drainage Project of Hui-an Nuclear Power Plant［J］.Journal of Waterway and Harbor，2007，25（3）：159-162.

［3］周美林，肖政，蔣昌波.沿河公路丁壩群水毀防護平面二維水流數值模擬研究［J］.水運工程，2007（8）：17-20.ZHOU M L，XIAO Z，JIANG C B.Numerical Simulation of Planar 2-D Flow with Spur-dikes Protecting Riverside Subgrade［J］.Port&Waterway Engineering，2007（8）：17-20.

［4］曹繼文，陳惠泉，賀益英.明渠岸邊橫向取水口的三維數值計算［J］.水利學報，2004（2）：119-125.CAO J W，CHEN H Q，HE Y Y.3-D numerical modeling of lateral intake flow in open channel［J］.Journal of Hydraulic Engineering，2004（2）：119-125.

［5］彭曉慶.廣州大學城雜用水廠取水工程對航道影響數學模型［J］.水運工程，2005，374（3）：70-73.PENG X Q.Mathematical Model Analysis of Influences on Waterway Due to Implementation of the Miscellaneous-Purpose Waterworks Intake Project of Guangzhou Higher Education Mega Centre［J］.Port&Waterway Engineering，2005，374（3）：70-73.

［6］高峰，陳漢寶.華能營口電廠取水工程的流場數值模擬［J］.水運工程，2007（7）：403-408.GAO F，CHEN H B.Numerical simulation for water intake project of Huaneng Yingkou Power Plant［J］.Port&Waterway Engineering，2007（7）：403-408.

［7］周美林，夏波，李世峰.南海市第二水廠擴建取水工程對通航影響平面二維數學模型研究［J］.長沙交通學院學報，2007，23（3）：20-25.ZHOU M L，XIA B，LI S F.Numerical model analysis of influences on navigation conditions due to the water intake project of No.2［J］.Journal of Changsha Communications University，2007，23（3）：20-25.

［8］謝翠松，余蔚卿，余明輝.碼頭和取水口建設對防洪影響的計算與分析［J］.南水北調與水利科技，2005，3（4）：54-58.XIE C S，YU W Q，YU M H.Calculation and Analysis of Influence on Flood Control by Port and Intake Construction［J］.Southto-North Water Transfers and Water Science&Technology，2005，3（4）：54-58.

［9］闕志夏.平面二維水流數學模型在碼頭和取水口建筑物防洪影響計算分析中的應用［J］.水利科技與經濟，2007，13（9）：650-652.