王明輝，金春林，關(guān)耀華

（1.國核工程有限公司，上海 200233；2.國網(wǎng)吉林省電力有限公司電力科學研究院，長春 130021；3.吉林電力技術(shù)開發(fā)公司，長春 130021）

非能動壓水堆核電機組的保護和安全監(jiān)視系統(tǒng)（PMS）作為執(zhí)行核反應(yīng)堆停堆保護和安全專設(shè)的重要系統(tǒng)，雖然經(jīng)過細致的設(shè)計、試驗、驗證和嚴格的核安全審批，但作為尚無同類型投入運行非能動壓水堆核電機組的專用系統(tǒng)，存在技術(shù)成熟但尚無實際運行經(jīng)驗的問題。基于上述原因，非能動壓水堆核電機組的PMS 軟件處于不斷完善升級過程中，因此，對PMS軟件邏輯的算法分析和試驗驗證尤為重要。壓力、溫度、液位等監(jiān)測與設(shè)備機械極限相關(guān)，直接影響反應(yīng)堆熱傳遞能力的停堆保護變量，而超溫ΔT保護和超功率ΔT保護更為復雜，ΔT不能直接測量，需要從其他參數(shù)計算得到，以下對超溫ΔT保護和超功率ΔT保護的邏輯、算法合理性、安全性和可靠性進行了分析。

1 超溫ΔT、超功率ΔT 停堆保護的意義

反應(yīng)堆熱工水力設(shè)計在堆正常運行期間基本任務(wù)是要保證在正常運行，包括瞬態(tài)條件下，有效帶走堆內(nèi)熱量，防止燃料破損。其主要考慮因素包括偏離泡核沸騰（DNB）、燃料和包殼溫度、熱管因子［1］，而超溫ΔT和超功率ΔT與這三個因素有直接關(guān)系，超溫ΔT停堆防止偏離泡核沸騰的發(fā)生，超功率ΔT停堆則防止燃料和包殼溫度達到燒毀值。超溫ΔT保護停堆的意義在于保證正常運行或瞬態(tài)條件下燃料棒上的熱流密度不會大于極限最大熱流密度，以防止傳熱能力下降導致加熱壁燒毀。超功率ΔT保護停堆的意義在于保證正常運行或瞬態(tài)條件下核燃料裂變產(chǎn)生的熱量不會大于極限最大熱量，以防止燃料本身和包殼溫度超限導致燃料和包殼燒毀。

從目的上看，超溫ΔT保護和超功率ΔT保護都屬于堆芯熱導出導致的停堆保護，以堆芯熱功率為過程參數(shù)。由于保護原理上的不同，超溫ΔT保護和超功率ΔT保護采用不同的設(shè)定值。

2 邏輯算法分析

2.1 超溫ΔT 停堆保護

超溫ΔT停堆保護防止由于壓力、功率、冷卻劑溫度、軸向功率分布瞬態(tài)組合造成的DNB。ΔT堆芯功率計算基于冷卻劑的熱段和冷段溫度及系統(tǒng)壓力的熱力學性能，包括從堆芯到熱段溫度探測器和冷段到堆芯的管道流體瞬態(tài)延遲動態(tài)補償。

變量設(shè)定值根據(jù)冷段溫度、穩(wěn)壓器壓力和軸向功率形狀計算，在正常功率分布情況下，超溫ΔT設(shè)定值不超過DNB堆芯熱設(shè)計極限。如果上下功率量程中子通量偏差代表的軸向功率分布不均勻使局部DNB峰值超過堆芯熱設(shè)計極限，則超溫ΔT設(shè)定值自動減小。如果堆芯熱功率有一個變量超過設(shè)定值，則發(fā)生超溫ΔT停堆動作。功能邏輯見圖1?；舅惴ǎ?/p>

a.超溫ΔT停堆裕量小于0觸發(fā)停堆；

b.超溫ΔT停堆裕量等于超溫ΔT設(shè)定值（ΔTOTSP）減手動可調(diào)偏置值（COT）減根據(jù)熱段溫度TH、冷段溫度TC、穩(wěn)壓器壓力PPRZ計算出的堆芯熱功率過程值（qΔT）；

c.超溫ΔT設(shè)定值計算式：

2.1.1 超溫ΔT設(shè)定值ΔTOTSP

DNB堆芯熱設(shè)計極限：

式中：PPRZ為穩(wěn)壓器壓力相關(guān)的變量設(shè)定值；TC為冷段溫度相關(guān)的變量設(shè)定值。

使用插入法計算DNB 熱設(shè)計極限的壓力和冷段溫度折線函數(shù)。f（PPRZ，TC）折線函數(shù)定義為在6組壓力下，每組壓力點和7點溫度對應(yīng)的堆芯設(shè)計極限功率構(gòu)成的設(shè)計數(shù)據(jù)表。在6組壓力和每組7點溫度內(nèi)，采用線性插入法計算當前熱設(shè)計極限。當壓力和溫度超過設(shè)計數(shù)據(jù)表上下限后，采用上限或下限數(shù)據(jù)，不采用外推法［2］。

圖1 超溫ΔT 超功率ΔT 保護停堆功能邏輯

圖2 DNB堆芯熱設(shè)計極限

f1（ΔI）用來校正正常軸向功率分布情況下由軸向功率分布偏差引起的DNB堆芯熱設(shè)計極限偏移。因為軸向功率分布的正向和負向變化都會造成軸向局部溫度明顯高于平均溫度，所以，需要對正負軸向功率分布偏差進行懲罰性補償校正。軸向核功率（中子通量）偏差懲罰補償f1（ΔI）是關(guān)于ΔI的函數(shù)，用線性插入法實現(xiàn)。

ΔI為經(jīng)過密度修正后的上部和下部功率量程中子通量探測器測量值之間的偏差，也就是實際的核功率偏差。中子通量由功率量程核測儀表測量，并由反應(yīng)堆入口介質(zhì)密度變化影響造成的中子衰減進行補償。軸向中子通量偏差補償見圖3。

圖3 軸向中子通量偏差補償

由于壓力容器內(nèi)流經(jīng)堆芯和堆外核測裝置間降水段水密度的改變會造成的吸收中子能力的變化，從而導致堆外中子流量探測器測量值偏離實際值，所以進行中子通量的密度修正是必須的［3］。

式中：FNA為衰減因子；A、B為系數(shù)；ρ（PPRZ，TC）為流體密度。

系數(shù)A、B由兩點確定，其中一點為正常滿功率時衰減因子FNA＝1，ρ（PPRZ，TC）是與穩(wěn)壓器壓力（代表系統(tǒng)壓力）及冷段溫度（代表反應(yīng)堆內(nèi)降水段流體溫度）相對應(yīng)的流體密度。

2.1.2 堆芯熱功率過程值qΔT

堆芯熱功率根據(jù)熱段、冷段溫度和穩(wěn)壓器壓力計算，其代表意義為堆芯出入口流體焓變所吸收的熱量。

式中：qΔT為ΔT熱功率/額定熱功率×100%；ρ（TC，PPRZ）為特定溫度和絕對壓力下冷段冷卻劑密度；h（TH，PPRZ）為特定溫度和絕對壓力下熱段冷卻劑焓值；h（TC，PPRZ）為特定溫度和絕對壓力下冷段冷卻劑焓值；C為偏置系數(shù)，用于零功率時調(diào)零；ΔT為轉(zhuǎn)換系數(shù)，等于正常額定熱功率。

焓值和密度算法：根據(jù)國際水和水蒸氣熱力性質(zhì)學會通過的工業(yè)用水和水蒸氣熱力性質(zhì)計算公式IAPWS－IF97，并結(jié)合本保護計算所適用的區(qū)間計算。

窄量程冷段溫度信號TC的處理：使用使用控制系統(tǒng)中的冷段“溫度生成/選擇模塊”用來對傳感器輸入溫度信號進行處理，以生成真實的溫度信號。模塊主要功能有模數(shù)轉(zhuǎn)換、量程和硬件品質(zhì)檢查、濾波、信號質(zhì)量判斷、兩點輸入溫度選擇處理。

溫度信號選擇處理過程中對兩點溫度質(zhì)量良好且彼此偏差在允許范圍內(nèi)信號做加權(quán)平均運算：

式中：WC為TC加權(quán)系數(shù)；為經(jīng)過預處理的溫度信號。

對彼此偏差較大信號或單點質(zhì)量良好需經(jīng)過蒸汽管線壓力來驗證（檢查）合格后代表冷段溫度TC。用蒸汽管線壓力信號驗證冷管段溫度信號含義如下：穩(wěn)態(tài)工況下蒸汽發(fā)生器熱量交換時，一回路放熱后冷段冷卻劑溫度和二回路吸收熱后飽和蒸汽溫度成正相關(guān)，而二回路蒸汽飽和溫度和飽和壓力具有單一對應(yīng)關(guān)系，因此，用二回路蒸汽管線計算出的冷段溫度可簡單表述如下：

式中：TSG（PSG）為蒸汽管線飽和壓力對應(yīng)的飽和溫度；KSG為正常滿功率穩(wěn)定狀態(tài)運行時冷段溫度和TSG（PSG）溫度差。

窄量程熱段溫度信號TH的處理：熱段溫度信號TH的處理邏輯和冷段溫度處理邏輯類似。

當前熱功率百分比qΔT的作用，是對熱段溫度傳感器輸入信號預處理過程中根據(jù)不同功率修正溫度的測量偏差。

熱段溫度選擇計算根據(jù)信號質(zhì)量進行加權(quán)平均計算：

式中：WH為TH加權(quán)系數(shù)；為熱段溫度信號。

2.1.3 手動可調(diào)偏置值COT

為防止反應(yīng)性當量最大的一束控制棒卡在堆芯上部不能下落，保證安全停堆，超溫ΔT停堆保護必須設(shè)置最小停堆裕量。COT可以減小超溫ΔT設(shè)定值，其目的是人為調(diào)整正常額定熱功率停堆設(shè)定值和DNB設(shè)計極限的差值，以保證足夠停堆裕量。

2.2 超功率ΔT 停堆保護

超功率ΔT停堆保護防止由于燃料過度熱量產(chǎn)生導致的燃料棒及燃料包殼破損。

如果堆芯熱功率接近燃料燒毀臨界值，則發(fā)生超功率ΔT停堆動作。算法表達如下：

a.超功率ΔT停堆裕量小于0觸發(fā)停堆；

b.超功率ΔT停堆裕量等于超功率ΔT設(shè)定值（ΔTOPSP）減 根 據(jù)TH、TC、PPRZ計 算 出 的 堆 芯 功 率（qΔT）；

c.超功率ΔT設(shè)定值（ΔTOPSP）由一個常量和軸向功率分布補償組成：

3 非能動壓水堆超溫ΔT、超功率ΔT 保護的其他技術(shù)特點

以上保護邏輯運算由PMS 中的一個序列實現(xiàn)，完整的保護具有4個序列，從輸入到邏輯運算再到輸出在軟件和硬件上完全獨立，最終輸出采用四取二邏輯處理控制停堆設(shè)備；每個序列的邏輯運算由互為冗余的兩個雙穩(wěn)態(tài)處理邏輯系統(tǒng)完成。雖然執(zhí)行相同的保護目的，但是不同于第二代核電反應(yīng)堆普遍采用過程變量比例關(guān)系對比算法，非能動壓水堆ΔT熱功率計算采用了不同的數(shù)學模型，主要使用焓和焓變來計算ΔT的熱量變化［4－5］。

非能動壓水堆ΔT保護研究和應(yīng)用的前提，是針對反應(yīng)堆正常運行、功率變化平緩情況下的超功率保護，不適合事故情況下瞬態(tài)突變?？焖俚墓β首兓善渌斯β屎秃斯β首兓实韧６压δ軐崿F(xiàn)。

壓水堆功率運行模式時，反應(yīng)堆冷卻劑泵運行在額定轉(zhuǎn)速，冷段體積流量基本恒定，基于此種原因，ΔT功率計算沒有將冷卻劑流量和泵轉(zhuǎn)速納入計算，極大地簡化了數(shù)學模型。由于事故導致冷卻劑流量改變的瞬態(tài)必然會影響實際熱導出，造成qΔT計算值偏離真實值，非能動壓水堆采用的解決方案是通過另外的流量低和泵轉(zhuǎn)速低停堆保護功能實現(xiàn)停堆，良好的范圍重疊可以很好地解決流量瞬態(tài)問題。

4 結(jié)論

通過對超溫ΔT和超功率ΔT功能和數(shù)學模型的研究，分析了邏輯算法在反應(yīng)堆偏離正常運行工況下的保護計算原理和實現(xiàn)方式，有獨特優(yōu)越性。

a.采用獨特的焓變和熱功率算法，比國內(nèi)眾多核電廠使用溫升、轉(zhuǎn)速比等直接比例計算方式有更好的綜合計算精度。

b.數(shù)學模型和邏輯算法合理，能滿足精度和響應(yīng)時間要求，方式可靠，能保證核反應(yīng)堆安全運行。

c.有效的參數(shù)極限值處理，除溫度外直接或間接參數(shù)取用保守值，對溫度參數(shù)進行多種方式質(zhì)量判斷和選擇計算，有效降低了極端情況下運行工況惡化的可能性。

d.序列內(nèi)邏輯算法相對獨立，包括邏輯算法功能的獨立性，序列間四取二邏輯的最終觸發(fā)，保證了任何一個序列發(fā)生問題都不會錯誤觸發(fā)停堆。

［1］ 林誠格.非能動安全先進核電廠AP1000［M］.北京：原子能出版社，2008.

［2］ 于平安，朱瑞安.核反應(yīng)堆熱工分析［M］.上海：上海交通大學出版社，1979.

［3］ 謝仲生.核反應(yīng)堆物理分析［M］.西安：西安交通大學出版社，2004.

［4］ 劉宏春，王濤濤.嶺澳二期核電站數(shù)字化反應(yīng)堆保護系統(tǒng)［J］.核動力工程，2008，（1），1－2.

［5］ 周海翔，王衛(wèi)國.田灣核電站數(shù)字化反應(yīng)堆保護系統(tǒng)可靠性分析［J］.核電子學與探測技術(shù)，2009，（11）：1272.