核動力裝置一回路自然循環能力最優化設計

郝承明，曲自信，賴建永，喻 巧，張 皓，羅涵禹，趙 京，孫 燕，梁鐵波，汪 宇，嚴思偉

(中國核動力研究設計院 反應堆系統設計技術國家級重點實驗室，成都 610213)

世界上幾次重大核事故已經證明，當反應堆發生事故時，依靠能動的安全系統、復雜的控制邏輯、以及操縱人員的現場判斷往往存在一定的風險[1].對于壓水反應堆核動力裝置而言，自然循環能力是其安全性能的重要指標[2-3].一般情況下反應堆堆芯與蒸汽發生器之間存在顯著的標高差，這樣一回路冷卻劑具備了形成自然循環的必要條件[4-5].此外，堆芯流體溫度、蒸汽發生器壓力、堆芯功率等參數也是影響反應堆自然循環能力的重要參數[6-7].探討如何合理確定設計參數，提升壓水反應堆一回路冷卻劑的自然循環能力，特別是提升自然循環工況下排走的堆芯功率水平，無疑對提升反應堆的固有安全性具有重要意義.

為獲得核動力裝置最佳設計參數，可將最優化理論運用到參數設計中[8-9].粒子群優化算法因其較好的非線性優化能力適用于多參數耦合的核動力系統[10].因此，本文以某壓水反應堆核動力裝置為研究對象，在采用堆芯冷卻劑平均溫度不變運行方案[11]條件下，建立了一回路冷卻劑自然循環能力評價模型，完成一回路運行壓力、堆芯冷卻劑進出口溫度以及蒸汽發生器壓力等參數對一回路冷卻劑自然循環能力的影響敏感性分析，并以上述變量構成設計參數向量，采用粒子群優化算法，對目標核動力系統的一回路冷卻劑自然循環排熱能力進行了優化設計.

1 自然循環計算模型

1.1 基本假設

在建立如圖1所示壓水堆核動力裝置一回路冷卻劑自然循環計算模型時，本文采用了如下主要假設：

圖1 壓水反應堆一回路系統示意圖

1)只考慮給定邊界條件下一回路冷卻劑達到的穩態自然循環狀態，即按照穩態工況來計算自然循環；

2)不考慮冷卻劑在流動截面上的參數分布，即認為冷卻劑在一回路中的流動是一維流動；

3)若反應堆內冷卻劑出現沸騰現象，則汽相處于飽和狀態；

4)無論冷卻劑是單相流動還是兩相流動，流體的可壓縮性可以忽略不計；

5)除了反應堆堆芯活性區以及蒸汽發生器U型傳熱管區域外，認為一回路冷卻劑在其他區域內的流動為絕熱流動.

1.2 物理模型

在本文所建立的自然循環流動模型中，采用一維兩相均相流模型描述兩相流動，基本守恒方程組由質量守恒方程、動量守恒方程及能量守恒方程組成.在所采用的基本假設條件下，基本守恒方程包括：

1)質量守恒方程

(1)

其中：W為質量流量，kg/s；x為幾何位置，m.單相流動中，下標“m”表示單相液體參數；兩相流動中，下標“m”表示氣液混合物參數.

2)動量守恒方程

(2)

其中：ρ為流體密度，kg/m3；A為流通橫截面積，m2；p為壓力，Pa；(?p/?z)f為摩擦壓降梯度，Pa/m；(?p/?z)g為重位壓降梯度，Pa/m；(?p/?z)l為局部壓降梯度，Pa/m.

3)能量守恒方程

(3)

其中：i為流體比焓，J/kg；q為管內線功率密度，W/m，絕熱段中q=0.

當流體處于單相液體狀態或過冷沸騰狀態時，汽相的物性參數采用飽和物性，而液相的主要物性參數可以由當地壓力和比焓查水和水蒸汽熱力性質表獲得，如液相密度可表示為：

ρL=ρ(p,i)

(4)

當流體處于飽和沸騰兩相流動狀態時，飽和液體和飽和蒸汽的物性參數均為當地壓力下的飽和參數，即：

ρg=ρ(p)ρf=ρ(p)

(5)

ig=i(p)if=i(p)

(6)

根據守恒原理，容易獲得兩相混合物的折算參數，主要折算參數可表示為：

Wm=ρmuA=GA

(7)

im=xig+(1-x)if

(8)

ρm=αρg+(1-α)ρf

(9)

其中：G為質量流速，kg/(m2·s)；u為流體流速，m/s；α為空泡份額；x為質量含氣率；ρg、ρf分別為飽和蒸汽和飽和水的密度，kg/m3；ig、if分別為飽和蒸汽和飽和水的比焓，J/kg.

在均相流模型中，空泡份額與質量含汽率之間的關系可表示為：

(10)

除上述守恒方程外，為了確定一回路冷卻劑在自然循環條件下的流動和換熱情況，還需要補充動量方程和能量方程中各項分壓降及換熱相關的本構關系，其中關于壓降的本構關系式主要包括摩擦壓降、局部壓降和重位壓降.以下簡述各項壓降的計算方法.

4)摩擦壓降關系式

當流體處于單相流動狀態時，管內摩擦壓降為：

(11)

其中：d為管道內徑，m；fL為單相摩擦阻力系數，可通過下式計算：

(12)

當一回路中出現兩相流動時，冷卻劑沿流動通道上的摩擦壓降梯度可采用下面兩種方法的任一種計算：

(13a)

(13b)

其中：fm為氣液兩相摩擦阻力系數，可通過式(14)計算：

(14)

5)局部壓降關系式

回路中局部阻力損失的統一表達式為：

(15)

其中：ξ為局部阻力系數.

6)重位壓降關系式

回路中所有管道重位壓降梯度可統一表示為：

(16)

7)管內單相對流換熱系數

對于蒸發器傳熱管內的冷卻劑，當其流動形態為單相對流換熱時，管內換熱系數按式(17)計算.當流體的雷諾數在2 320～10 000之間時，對換熱系數進行線性插值處理.

(17)

其中：λL為管內流體的導熱系數，W/(m·K).

8) 棒束通道內的單相對流換熱

對于冷卻劑流過堆芯棒束型組件時的單相對流換熱，采用Weisman關系式[11]進行計算，即：

Nu=CRe0.8Pr1/3

(18)

其中：C為取決于燃料組件柵格的排列方式，對于正方形柵格，當1.1≤Pn/dcs≤1.3時，C=0.042(Pn/dcs)-0.024，Pn燃料元件中心距，dcs為燃料元件包殼外徑；對于三角形柵格，當1.1≤Pn/dcs≤1.5時，C=0.026(Pn/dcs)-0.006.

9) 泡核沸騰傳熱

對于處于欠熱泡核沸騰區和飽和泡核沸騰區的堆芯冷卻劑，采用Jens-Lottes經驗公式[12]預測熱負荷與壁面過熱度之間的關系：

(19)

其中：Tcs為包殼外表面溫度，℃；Tfs為流體的飽和溫度，℃；q為表面熱流密度，W/m2；p為流體工作壓強，MPa.

10)大容積沸騰

對于蒸汽發生器二次側的蒸發段，本文選用如下模型[13]進行計算其沸騰換熱系數：

α2b=0.557p0.15q0.7

(20)

其中：α2b為二次側飽和沸騰換熱系數，W/(m2K)；p為二次側運行絕對壓力，Pa；q為二次側蒸發段表面熱負荷，W/m2.

1.3 計算流程圖

圖2為本文進行壓水反應堆核動力裝置一回路冷卻劑自然循環能力計算的流程圖.

圖2 自然循環能力評估流程圖

2 粒子群優化算法

2.1 算法流程

粒子群優化算法(PSO)最早是由Kennedy和Eberhart于1995年提出的，其主要思想是借鑒鳥群飛行覓食行為，依賴個體(單只鳥)自身信息和群體(鳥群)信息，調整飛行方向和飛行速度，最終實現對目標(食物)的最優化，本文采用的PSO算法[14]流程圖如圖3所示.

圖3 粒子群算法流程圖

2.2 算法測試

為測試算法的有效性，本文利用PSO算法，對已公布的7個標準例題進行50次求解[15]，計算結果及其統計情況列于表1中.結果表明，粒子群算法可以正確求解表1中所示的標準例題.

表1 粒子群算法求解標準例題時的性能

3 結果與討論

3.1 參數敏感性分析

利用本文所建立的計算模型，對參考核動力裝置一回路冷卻劑的自然循環特性進行了參數敏感性計算，即在計算中保持其他參數不變，僅在給定范圍內改變單一變量的取值，確定在一回路冷卻劑達到穩態自然循環運行工況時所能帶走的堆芯熱功率，亦即本文所指的自然循環能力.本文在進行參數敏感性分析和優化分析時，對參數均采用了無量綱化處理，即當各無量綱參數為1時，表示與參考核動力裝置相應參數相同.圖4為單獨改變反應堆一回路運行壓力時，一回路冷卻劑在穩態自然循環工況下能夠排走的堆芯功率水平的相對變化情況.結果表明，當保持設計要求不變時，隨著一回路運行壓力的升高，無量綱自然循環能力也升高.主要原因是隨著一回路運行壓力升高，堆芯出口冷卻劑的允許工作溫度范圍也加大，在控制堆芯平均溫度不變的情況下，堆芯進出口冷卻劑溫度差變大，自然循環驅動力有所增加，進而使自然循環工況下的循環流量增加，在同等安全準則條件下，可以帶走的堆芯功率水平也增加，自然循環能力提升.

圖4 一回路壓力對自然循環能力影響

圖5所示為一回路冷卻劑的自然循環能力隨堆芯進口溫度的變化規律.在本文計算參數范圍內，一回路冷卻劑的自然循環能力隨堆芯進口溫度升高而升高.隨著額定工況下堆芯進口溫度升高，保持一回路壓力不變和堆芯出口冷卻劑溫度不變情況下，將導致主管道內徑增加、蒸汽發生器傳熱管長度減小，從而降低一回路總阻力；另一方面，隨著堆芯進口溫度的升高，堆芯進出口溫升減小，這又有使自然循環驅動力減小的趨勢.綜合來看，在本文參數范圍內，增加堆芯進口溫度可以提高一回路冷卻劑的自然循環能力.

圖5 堆芯進口溫度對自然循環能力影響

一回路冷卻劑自然循環能力隨堆芯出口冷卻劑溫度的升高逐漸降低，計算結果示見圖6.這是因為，當額定工況下的堆芯出口溫度升高時，堆芯平均溫度也隨之提高，額定工況下一回路冷卻劑的流量減小，相應的一回路主管道尺寸減小、蒸汽發生器傳熱管增長，一回路總阻力有增加趨勢；盡管增加堆芯出口溫度有使一回路冷熱芯之間的驅動力有增加趨勢，但在計算參數范圍內，堆芯出口冷卻劑溫度升高的總體效果是使自然循環能力降低.

圖6 堆芯出口溫度對自然循環能力影響

圖7為蒸汽發生器壓力變化時，一回路冷卻劑自然循環能力的變化情況.由計算結果可以看出，蒸汽發生器內的工作壓力對于一回路冷卻劑自然循環特性的影響較小.在計算參數范圍內，隨著蒸汽發生器工作壓力的增加，一回路冷卻劑自然循環能力略有降低.其原因是蒸汽發生器內工作壓力的增加導致二次側工作介質的飽和溫度增加，在保持一回路冷卻劑參數和二次側給水溫度不變情況下，這將使得蒸汽發生器傳熱管的預熱段增長，傳熱管的總長度也略有增加，進而使得一回路冷卻劑在流過蒸汽發生器傳熱管束時的流動阻力也有所增加，這樣使一回路冷卻劑在自然循環工況下的循環流量有降低趨勢，相應地使一回路冷卻劑自然循環能力也有所降低.

圖7 蒸汽發生器壓力對自然循環能力影響

3.2 優化計算結果

利用粒子群算法和本文開發的針對壓水反應堆核動力裝置一回路冷卻劑自然循環能力評估程序，開展了以一回路運行壓力、堆芯進口冷卻劑溫度、堆芯出口冷卻劑溫度以及蒸汽發生器壓力為優化變量的一回路冷卻劑自然循環能力最優化設計，計算結果見表2.

表2 PSO算法對一回路冷卻劑自然循環能力的優化設計結果

計算結果表明，在設計空間內，適當增加一回路系統壓力、反應堆堆芯進口溫度和蒸汽發生器壓力，適當減小堆芯出口溫度，可顯著提高參考核動力裝置一回路冷卻劑的自然循環能力.相比于參考核動力裝置，優化后的自然循環能力可提高約26.46%，優化效果顯著.

4 結 語

本文以壓水反應堆核動力裝置為研究對象，建立了給定邊界條件下一回路冷卻劑自然循環運行工況的評價模型，開發了相應的計算程序；完成了典型熱工參數對一回路冷卻劑自然循環能力的敏感性分析；利用粒子群算法實現了對一回路冷卻劑自然循環能力的參數最優化設計.通過本文研究，得到如下結論：

1)粒子群算法可以用于解決一回路冷卻劑自然循環能力這種多參數、強耦合最優化問題；

2)在額定工況下，提高一回路壓力、增加堆芯進口冷卻劑溫度、降低堆芯出口冷卻劑溫度有使一回路冷卻劑自然循環能力增加的趨勢；

3)在額定工況下，在所研究參數范圍內增加蒸汽發生器工作壓力對于一回路冷卻劑自然循環能力的影響很??；

4)在設計參數空間內，適當調整一回路壓力、堆芯進口冷卻劑溫度、堆芯出口冷卻劑溫度和蒸汽發生器壓力，可使參考核動力裝置的一回路冷卻劑自然循環能力增加25.46%.