海南核電一回路泄漏率的優化算法

李 靜，劉 逸

（海南核電有限公司，海南昌江 570105）

0 引言

海南核電一回路泄漏率通過容控箱液位、穩壓器液位和一回路平均溫度計算獲得，現有計算方法能夠簡單快速得到一回路泄漏率。這種方法也有缺點，計算結果波動大、可靠性不高，選點不同會導致計算結果偏差很大。產生這一現象的原因是計算過程中沒有充分利用該時間段內數據。為了解決這個問題，將時間段內的所有數據都用于計算一回路泄漏率。

1 一回路泄漏率改進前的算法

海南核電現有的一回路泄漏率計算方法是記錄T0時刻的穩壓器液位、容控箱液位和一回路平均溫度。然后記錄2 h 后的穩壓器液位、容控箱液位和一回路平均溫度，再通過式（1）計算得到一回路泄漏率Fp。

式中 N1——計算開始時刻的容控箱液位，m

N2——計算結束時刻的容控箱液位，m

N3——計算開始時刻的穩壓器液位，m

N4——計算結束時刻的穩壓器液位，m

T1——計算開始時刻的一回路平均溫度，℃

T2——計算結束時刻的一回路平均溫度，℃

Δt——計算開始到計算結束的時間間隔，固定值2 h

這種計算方法的前提條件是容控箱液位、穩壓器液位和一回路平均溫度為線性變化，或者近似為線性變化。但實際這些數據都不是按照線性變化，而且很多時刻波動較大，導致式（1）計算結果誤差較大。

以4 月14 日8：00—10：00 數據為例，說明現有計算方式的問題。容控箱液位2 h 內的變化曲線如圖1 所示，通過最小二乘曲線擬合得到的容控箱液位隨時間變化的解析表達式：y=0.010 5x-461.99，縱坐標單位m，橫坐標單位min。

穩壓器液位2 h 內的變化曲線如圖2 所示，通過最小二乘曲線擬合得到的穩壓器液位隨時間變化的解析表達式：y=-0.012 2x+536.82，縱坐標單位m，橫坐標單位min。

一回路平均溫度2 h 內的變化曲線如圖3 所示，通過最小二乘曲線擬合得到的一回路平均溫度隨時間變化的解析表達式：y=0.480 6x-20 808，縱坐標單位℃，橫坐標單位min。

圖1 容控箱液位變化曲線

圖2 穩壓器液位變化曲線

表1 通過PT1RCP010 規程計算得到的一回路泄漏率

圖3 一回路平均溫度變化曲線

通過式（1），對7：30—9：3 的泄漏率進行計算，為了說明現有計算方法的問題，分別計算每隔1 min 的一回路泄漏率。表1 的計算結果表明，不同取樣點計算結果偏差很大，在10 min 的間隔時間內，最小計算結果6.6 L/h，最大計算結果81 L/h，計算結果方差達到22.68。

2 一回路泄漏率改進后的算法

為了減小計算結果的方差，提高結算結果的穩定性，必須將2 h 內的所有數據都用于計算，采用曲線擬合方法可以解決這個問題。通過一階曲線擬合得到的容控箱液位解析表達式：y=0.010 5x-461.99，穩壓器液位解析表達式：y=-0.012 2x+536.82；一回路平均溫度解析表達式：y=0.480 6x-20 808。

通過曲線擬合方法獲得的相同初始時刻的一回路泄漏率，T0和T1時刻與表1 一致，泄漏率計算結果為6.182 4 L/h、6.410 6 L/h、7.307 6 L/h、7.902 3 L/h、8.428 9 L/h、8.885 9 L/h、9.097 5 L/h、9.316 0 L/h、9.811 8 L/h、10.751 L/h。計算結果的方差僅為1.468。通過對比可知，使用曲線擬合方法計算結果穩定性大幅提高。

3 結論

海南核電現有的一回路泄漏率計算方法具有較大的不穩定性，為提高計算結果的穩定性，本文將原有的兩點式計算方式改為曲線擬合計算方式。通過對比計算，可靠性大幅提高，計算結果的方差由22.68 下降到1.468。本方法不僅適用于海南核電，還適用于其他M310 機組。