本發(fā)明屬于核電接入電網(wǎng)的安全穩(wěn)定分析技術(shù)領(lǐng)域，涉及一種核電機(jī)組冷卻劑屏蔽主泵動(dòng)態(tài)建模及分析方法。

背景技術(shù)：

反應(yīng)堆冷卻劑主循環(huán)泵是核電站最重要的設(shè)備之一，也是整個(gè)冷卻劑系統(tǒng)中唯一的旋轉(zhuǎn)設(shè)備，具有功率大、工作條件苛刻的特點(diǎn)。

正常運(yùn)行情況下，核電廠的廠用電系統(tǒng)由電網(wǎng)供電，若電網(wǎng)側(cè)出現(xiàn)事故，導(dǎo)致核電廠不能從電網(wǎng)側(cè)獲取廠用電時(shí)，反應(yīng)堆會(huì)自動(dòng)安全停堆，但在反應(yīng)堆安全停堆過程中，需要反應(yīng)堆冷卻系統(tǒng)繼續(xù)工作，以將堆芯中殘余的熱量傳遞出去。以ap1000為代表的第三代先進(jìn)壓水堆核電機(jī)組冷卻劑主泵采用屏蔽電動(dòng)機(jī)，其飛輪和電動(dòng)機(jī)轉(zhuǎn)子一起包容在rcs的承壓部件內(nèi)，使其徑向尺寸受到極大限制，因此，相比于其他堆型，具有更小的轉(zhuǎn)動(dòng)慣量，冷卻劑流量對電網(wǎng)側(cè)波動(dòng)反映更為明顯，因此研究電網(wǎng)側(cè)波動(dòng)對冷卻劑流量造成的影響，針對根據(jù)分析結(jié)果，進(jìn)行繼電保護(hù)整定，停堆等措施提供系統(tǒng)安全可靠性。

目前，常用的冷卻劑主泵流量建模方法主要有，流體網(wǎng)格計(jì)算方法相關(guān)模型、基于四象限特性曲線的相關(guān)計(jì)算模型等。基于流體網(wǎng)格計(jì)算方法的冷卻劑主泵流量計(jì)算模型主要應(yīng)用于電站仿真系統(tǒng)中，將冷卻劑主泵的工作過程簡化為由管路、閥門等連接而成的流體網(wǎng)格系統(tǒng)，更多的是從泵的壓力特性入手分析建模，忽略了泵本身的電學(xué)特性，并不適用于電力系統(tǒng)分析，也無法反映電網(wǎng)側(cè)波動(dòng)對冷卻劑流量的影響。而基于四象限特性曲線的冷卻劑主泵流量計(jì)算模型雖然考慮了冷卻劑主泵的電學(xué)特性，但由于計(jì)算繁瑣，參數(shù)不易獲得，一般只用于獨(dú)立分析主泵特性采用，難以與其他設(shè)備模型及電力系統(tǒng)仿真進(jìn)行結(jié)合。

基于此，為了通過仿真研究電網(wǎng)波動(dòng)對核電機(jī)組安全運(yùn)行的影響，一方面要求冷卻劑屏蔽主泵仿真模型能夠精確快速地表征冷卻劑流量的波動(dòng)，以便于制定相應(yīng)控制措施；另一方面，仿真模型應(yīng)能夠接入電力系統(tǒng)仿真軟件和其他核電設(shè)備模型中進(jìn)行仿真運(yùn)算。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明的目的在于提供一種核電機(jī)組冷卻劑屏蔽主泵動(dòng)態(tài)建模及分析方法，該方法考慮了冷卻劑主泵的摩擦轉(zhuǎn)矩，能夠準(zhǔn)確實(shí)時(shí)反映電網(wǎng)側(cè)電壓、頻率波動(dòng)對冷卻劑流量影響，進(jìn)而準(zhǔn)確分析電網(wǎng)波動(dòng)對核電安全影響，彌補(bǔ)目前已公開方法對主泵惰轉(zhuǎn)和主泵卡死分析的不足，對接入電網(wǎng)的核電分析具有更高的實(shí)時(shí)性；同時(shí)，模型結(jié)構(gòu)清晰、參數(shù)意義明確、易獲取，具有很高的準(zhǔn)確性和實(shí)用性。

為實(shí)現(xiàn)上述目的，本發(fā)明的技術(shù)方案是：一種核電機(jī)組冷卻劑屏蔽主泵動(dòng)態(tài)建模及分析方法，包括如下步驟，

s1、基于核電機(jī)組冷卻劑屏蔽主泵特性，建立冷卻劑主泵流量計(jì)算模型，包括轉(zhuǎn)矩平衡方程、質(zhì)量流量方程；

s2、基于冷卻劑主泵設(shè)計(jì)參數(shù)及測試數(shù)據(jù)，計(jì)算或辨識得到主泵模型參數(shù)，并計(jì)算主泵電磁轉(zhuǎn)矩、水力轉(zhuǎn)矩和摩擦轉(zhuǎn)矩；

s3、基于步驟s1所建立的冷卻劑主泵流量計(jì)算模型，仿真計(jì)算主泵惰轉(zhuǎn)工況，電壓、頻率波動(dòng)工況時(shí)冷卻劑流量實(shí)時(shí)變化，與測試曲線作對比，校驗(yàn)?zāi)Ｐ图皡?shù)；

s4、將步驟s1所建立的冷卻劑主泵流量計(jì)算模型接入核電機(jī)組系統(tǒng)模型，并接入電力系統(tǒng)仿真軟件，實(shí)現(xiàn)核電機(jī)組整體動(dòng)態(tài)仿真分析。

在本發(fā)明一實(shí)施例中，在步驟s1中，所述冷卻劑主泵流量計(jì)算模型的轉(zhuǎn)矩平衡方程、質(zhì)量流量方程，是采用以下方法建立的：

冷卻劑屏蔽主泵轉(zhuǎn)矩平衡方程為：

式中，ip為冷卻劑主泵的轉(zhuǎn)動(dòng)慣量，ωp為冷卻劑主泵轉(zhuǎn)子角速度，mpe為冷卻劑主泵電磁轉(zhuǎn)矩，mpm為冷卻劑主泵受到的總阻力矩；

主泵電磁轉(zhuǎn)矩mpe計(jì)算表達(dá)式為：

式中，p為主泵電動(dòng)機(jī)磁極對數(shù)，r2為主泵電動(dòng)機(jī)轉(zhuǎn)子每相繞組電阻，s為主泵電動(dòng)機(jī)轉(zhuǎn)差率，m1為主泵電動(dòng)機(jī)磁極對數(shù)，u1n為電網(wǎng)側(cè)系統(tǒng)基準(zhǔn)電壓，f1n為電網(wǎng)側(cè)系統(tǒng)基準(zhǔn)頻率，上述各參數(shù)均為設(shè)計(jì)常數(shù)；xσ為主泵電動(dòng)機(jī)轉(zhuǎn)子漏電抗，u1*、f1*分別為實(shí)際電網(wǎng)電壓標(biāo)幺值和頻率標(biāo)幺值；k1、k2為與冷卻劑主泵相關(guān)的常系數(shù)；

屏蔽主泵因其飛輪和轉(zhuǎn)子固封于承壓部件內(nèi)，使得其轉(zhuǎn)子轉(zhuǎn)動(dòng)慣量小于常規(guī)泵；轉(zhuǎn)速對轉(zhuǎn)矩變化靈敏度高，摩擦轉(zhuǎn)矩不可忽略，因此屏蔽主泵的阻力矩應(yīng)同時(shí)計(jì)及冷卻劑流量產(chǎn)生的阻力矩mp1和由于機(jī)械摩擦造成的摩擦阻力矩mp2，即有

mpm＝mp1+mp2(3)

冷卻劑流量產(chǎn)生的阻力矩mp1表達(dá)式為：

式中，d為冷卻劑質(zhì)量流量，ωp為為冷卻劑主泵轉(zhuǎn)子角速度，η為冷卻劑主泵效率，h為冷卻劑主泵揚(yáng)程；

mp2表示主泵機(jī)械摩擦引起的阻力矩，采用經(jīng)驗(yàn)公式表示為

式中，k3為常數(shù)，由泵的實(shí)驗(yàn)值確定；n取1.7-2.4；

冷卻劑質(zhì)量流量正比于主泵轉(zhuǎn)子角速度，但由于系統(tǒng)慣性，冷卻劑流量變化滯后于轉(zhuǎn)子角速度變化，將冷卻劑質(zhì)量流量方程表示為

式中，為轉(zhuǎn)子角速度標(biāo)幺值，d^*為主泵質(zhì)量流量標(biāo)幺值，τ為常數(shù)，取0.2s。

在本發(fā)明一實(shí)施例中，在步驟s2中，基于冷卻劑主泵設(shè)計(jì)參數(shù)及測試數(shù)據(jù)，計(jì)算或辨識得到主泵模型參數(shù)的具體實(shí)現(xiàn)如下：

s21、基于冷卻劑主泵惰轉(zhuǎn)工況測試數(shù)據(jù)，采用matlab優(yōu)化設(shè)計(jì)工具箱進(jìn)行參數(shù)辨識計(jì)算，辨識出主泵摩擦阻力矩相關(guān)系數(shù)k3，n；

s22、基于電網(wǎng)側(cè)頻率或電壓波動(dòng)情況下，冷卻劑流量測試數(shù)據(jù)，并將主泵摩擦阻力矩相關(guān)系數(shù)k3，n作為已知參數(shù)，采用matlab優(yōu)化設(shè)計(jì)工具箱進(jìn)行參數(shù)辨識，得到冷卻劑主泵電磁轉(zhuǎn)矩相關(guān)系數(shù)k1、k2；

s23、基于辨識的參數(shù)，建立相應(yīng)的電磁轉(zhuǎn)矩，摩擦轉(zhuǎn)矩模型方程。

在本發(fā)明一實(shí)施例中，所述步驟s3的具體實(shí)現(xiàn)過程如下：

s31、仿真冷卻劑主泵惰轉(zhuǎn)工況下，根據(jù)冷卻劑流量變化曲線，求解惰轉(zhuǎn)時(shí)間常數(shù)，與相關(guān)模型及參數(shù)進(jìn)行比較分析，驗(yàn)證所述冷卻劑主泵流量計(jì)算模型的有效性；

s32、仿真電網(wǎng)側(cè)極端事故情況下，根據(jù)冷卻劑流量與轉(zhuǎn)速變化趨勢，驗(yàn)證主泵低轉(zhuǎn)速保護(hù)優(yōu)先級高于冷卻劑低流量保護(hù)，驗(yàn)證所述冷卻劑主泵流量計(jì)算模型的合理性；

s33、仿真電網(wǎng)電壓、頻率階躍波動(dòng)工況下，根據(jù)冷卻劑流量變化，驗(yàn)證冷卻劑主泵流量計(jì)算模型對電網(wǎng)側(cè)電壓、頻率擾動(dòng)響應(yīng)仿真的準(zhǔn)確性。

在本發(fā)明一實(shí)施例中，在步驟s4中，實(shí)現(xiàn)核電機(jī)組整體動(dòng)態(tài)仿真分析具體為分析電網(wǎng)側(cè)波動(dòng)情況下，核電側(cè)主泵流量，燃料溫度，蒸汽壓力相關(guān)參數(shù)動(dòng)態(tài)響應(yīng)過程。

相較于現(xiàn)有技術(shù)，本發(fā)明具有以下有益效果：

本發(fā)明方法考慮了冷卻劑主泵的的摩擦轉(zhuǎn)矩，能夠準(zhǔn)確實(shí)時(shí)反映電網(wǎng)側(cè)電壓、頻率波動(dòng)對冷卻劑流量影響，進(jìn)而準(zhǔn)確分析電網(wǎng)波動(dòng)對核電安全影響，彌補(bǔ)目前已公開方法對主泵惰轉(zhuǎn)和主泵卡死分析的不足，對接入電網(wǎng)的核電分析具有更高的實(shí)時(shí)性。同時(shí)，模型結(jié)構(gòu)清晰、參數(shù)意義明確、易獲取，具有很高的準(zhǔn)確性和實(shí)用性。

附圖說明

圖1是本發(fā)明的冷卻劑屏蔽主泵動(dòng)態(tài)建模分析流程圖。

圖2是本發(fā)明的核電機(jī)組冷卻劑屏蔽主泵電動(dòng)機(jī)等效示意圖。

圖3是摩擦阻力矩相關(guān)參數(shù)k3、n辨識結(jié)果圖。

圖4是摩擦轉(zhuǎn)矩參數(shù)擬合曲線圖。

圖5是電磁轉(zhuǎn)矩相關(guān)參數(shù)k1、k2擬合結(jié)果圖。

圖6是冷卻劑主泵流量計(jì)算模型傳輸框圖。

圖7是t＝10s時(shí)主泵失去外電源冷卻劑流量變化曲線圖。

圖8是冷卻劑主泵電動(dòng)機(jī)供電頻率以4hz/s下降事故時(shí)堆芯流量與轉(zhuǎn)速歸一化曲線圖。

圖9是電壓/頻率階躍擾動(dòng)下冷卻劑流量變化曲線圖。

圖10是電網(wǎng)頻率下降時(shí)，核電機(jī)組動(dòng)態(tài)響應(yīng)圖。

具體實(shí)施方式

下面結(jié)合附圖，對本發(fā)明的技術(shù)方案進(jìn)行具體說明。

本發(fā)明的一種核電機(jī)組冷卻劑屏蔽主泵動(dòng)態(tài)建模及分析方法，包括如下步驟，

s1、基于核電機(jī)組冷卻劑屏蔽主泵特性，建立冷卻劑主泵流量計(jì)算模型，包括轉(zhuǎn)矩平衡方程、質(zhì)量流量方程；

s2、基于冷卻劑主泵設(shè)計(jì)參數(shù)及測試數(shù)據(jù)，計(jì)算或辨識得到主泵模型參數(shù)，并計(jì)算主泵電磁轉(zhuǎn)矩、水力轉(zhuǎn)矩和摩擦轉(zhuǎn)矩；

s3、基于步驟s1所建立的冷卻劑主泵流量計(jì)算模型，仿真計(jì)算主泵惰轉(zhuǎn)工況，電壓、頻率波動(dòng)工況時(shí)冷卻劑流量實(shí)時(shí)變化，與測試曲線作對比，校驗(yàn)?zāi)Ｐ图皡?shù)；

s4、將步驟s1所建立的冷卻劑主泵流量計(jì)算模型接入核電機(jī)組系統(tǒng)模型，并接入電力系統(tǒng)仿真軟件，實(shí)現(xiàn)核電機(jī)組整體動(dòng)態(tài)仿真分析。

在步驟s1中，所述冷卻劑主泵流量計(jì)算模型的轉(zhuǎn)矩平衡方程、質(zhì)量流量方程，是采用以下方法建立的：

冷卻劑屏蔽主泵轉(zhuǎn)矩平衡方程為：

式中，ip為冷卻劑主泵的轉(zhuǎn)動(dòng)慣量，ωp為冷卻劑主泵轉(zhuǎn)子角速度，mpe為冷卻劑主泵電磁轉(zhuǎn)矩，mpm為冷卻劑主泵受到的總阻力矩；

主泵電磁轉(zhuǎn)矩mpe計(jì)算表達(dá)式為：

式中，p為主泵電動(dòng)機(jī)磁極對數(shù)，r2為主泵電動(dòng)機(jī)轉(zhuǎn)子每相繞組電阻，s為主泵電動(dòng)機(jī)轉(zhuǎn)差率，m1為主泵電動(dòng)機(jī)磁極對數(shù)，u1n為電網(wǎng)側(cè)系統(tǒng)基準(zhǔn)電壓，f1n為電網(wǎng)側(cè)系統(tǒng)基準(zhǔn)頻率，上述各參數(shù)均為設(shè)計(jì)常數(shù)；xσ為主泵電動(dòng)機(jī)轉(zhuǎn)子漏電抗，u1*、f1*分別為實(shí)際電網(wǎng)電壓標(biāo)幺值和頻率標(biāo)幺值；k1、k2為與冷卻劑主泵相關(guān)的常系數(shù)；

屏蔽主泵因其飛輪和轉(zhuǎn)子固封于承壓部件內(nèi)，使得其轉(zhuǎn)子轉(zhuǎn)動(dòng)慣量小于常規(guī)泵；轉(zhuǎn)速對轉(zhuǎn)矩變化靈敏度高，摩擦轉(zhuǎn)矩不可忽略，因此屏蔽主泵的阻力矩應(yīng)同時(shí)計(jì)及冷卻劑流量產(chǎn)生的阻力矩mp1和由于機(jī)械摩擦造成的摩擦阻力矩mp2，即有

mpm＝mp1+mp2(3)

冷卻劑流量產(chǎn)生的阻力矩mp1表達(dá)式為：

式中，d為冷卻劑質(zhì)量流量，ωp為冷卻劑主泵轉(zhuǎn)子角速度，η為冷卻劑主泵效率，h為冷卻劑主泵揚(yáng)程；

mp2表示主泵機(jī)械摩擦引起的阻力矩，采用經(jīng)驗(yàn)公式表示為

式中，k3為常數(shù)，由泵的實(shí)驗(yàn)值確定；n取1.7-2.4；

冷卻劑質(zhì)量流量正比于主泵轉(zhuǎn)子角速度，但由于系統(tǒng)慣性，冷卻劑流量變化滯后于轉(zhuǎn)子角速度變化，將冷卻劑質(zhì)量流量方程表示為

式中，為轉(zhuǎn)子角速度標(biāo)幺值，d^*為主泵質(zhì)量流量標(biāo)幺值，τ為常數(shù)，取0.2s。

在步驟s2中，基于冷卻劑主泵設(shè)計(jì)參數(shù)及測試數(shù)據(jù)，計(jì)算或辨識得到主泵模型參數(shù)的具體實(shí)現(xiàn)如下：

s21、基于冷卻劑主泵惰轉(zhuǎn)工況測試數(shù)據(jù)，采用matlab優(yōu)化設(shè)計(jì)工具箱進(jìn)行參數(shù)辨識計(jì)算，辨識出主泵摩擦阻力矩相關(guān)系數(shù)k3，n；

s22、基于電網(wǎng)側(cè)頻率或電壓波動(dòng)情況下，冷卻劑流量測試數(shù)據(jù)，并將主泵摩擦阻力矩相關(guān)系數(shù)k3，n作為已知參數(shù)，采用matlab優(yōu)化設(shè)計(jì)工具箱進(jìn)行參數(shù)辨識，得到冷卻劑主泵電磁轉(zhuǎn)矩相關(guān)系數(shù)k1、k2；

s23、基于辨識的參數(shù)，建立相應(yīng)的電磁轉(zhuǎn)矩，摩擦轉(zhuǎn)矩模型方程。

所述步驟s3的具體實(shí)現(xiàn)過程如下：

s31、仿真冷卻劑主泵惰轉(zhuǎn)工況下，根據(jù)冷卻劑流量變化曲線，求解惰轉(zhuǎn)時(shí)間常數(shù)，與相關(guān)模型及參數(shù)進(jìn)行比較分析，驗(yàn)證所述冷卻劑主泵流量計(jì)算模型的有效性；

s32、仿真電網(wǎng)側(cè)極端事故情況下，根據(jù)冷卻劑流量與轉(zhuǎn)速變化趨勢，驗(yàn)證主泵低轉(zhuǎn)速保護(hù)優(yōu)先級高于冷卻劑低流量保護(hù)，驗(yàn)證所述冷卻劑主泵流量計(jì)算模型的合理性；

s33、仿真電網(wǎng)電壓、頻率階躍波動(dòng)工況下，根據(jù)冷卻劑流量變化，驗(yàn)證冷卻劑主泵流量計(jì)算模型對電網(wǎng)側(cè)電壓、頻率擾動(dòng)響應(yīng)仿真的準(zhǔn)確性。

在步驟s4中，實(shí)現(xiàn)核電機(jī)組整體動(dòng)態(tài)仿真分析具體為分析電網(wǎng)側(cè)波動(dòng)情況下，核電側(cè)主泵流量，燃料溫度，蒸汽壓力相關(guān)參數(shù)動(dòng)態(tài)響應(yīng)過程。

以下為本發(fā)明的具體實(shí)現(xiàn)過程。

本實(shí)施例提供了一種核電機(jī)組冷卻劑屏蔽主泵動(dòng)態(tài)建模及分析方法，如圖1所示，包括以下步驟：

步驟s1：基于核電機(jī)組冷卻劑屏蔽主泵特性，建立冷卻劑主泵流量計(jì)算模型，包括轉(zhuǎn)矩平衡方程、質(zhì)量流量方程。

冷卻劑屏蔽主泵轉(zhuǎn)矩平衡方程為：

在(1)式中，ip為冷卻劑主泵的轉(zhuǎn)動(dòng)慣量，ωp為冷卻劑主泵轉(zhuǎn)子角速度，mpe為冷卻劑主泵電磁轉(zhuǎn)矩，mpm為冷卻劑主泵受到的總阻力矩，包含兩部分，一部分為冷卻劑流量產(chǎn)生的阻力矩mp1，另一部分為主泵電動(dòng)機(jī)摩擦阻力矩mp2。

本實(shí)施例中，冷卻劑屏蔽主泵等效電路如圖2所示，主泵電磁轉(zhuǎn)矩mpe計(jì)算表達(dá)式為：

式中，p為主泵電動(dòng)機(jī)磁極對數(shù)，r2為主泵電動(dòng)機(jī)轉(zhuǎn)子每相繞組電阻，s為主泵電動(dòng)機(jī)轉(zhuǎn)差率，m1為主泵電動(dòng)機(jī)磁極對數(shù)，u1n為電網(wǎng)側(cè)系統(tǒng)基準(zhǔn)電壓，f1n為電網(wǎng)側(cè)系統(tǒng)基準(zhǔn)頻率，均為設(shè)計(jì)常數(shù)。xσ為主泵電動(dòng)機(jī)轉(zhuǎn)子漏電抗，u1*、f1*分別為實(shí)際電網(wǎng)電壓標(biāo)幺值和頻率標(biāo)幺值。k1、k2為與冷卻劑主泵相關(guān)的常系數(shù)，在本發(fā)明中采用已有實(shí)驗(yàn)數(shù)據(jù)進(jìn)行辨識。

冷卻劑流量產(chǎn)生的阻力矩mp1表達(dá)式為：

式中d為冷卻劑質(zhì)量流量，ωp為冷卻劑主泵轉(zhuǎn)子角速度，η為主泵效率，h為主泵揚(yáng)程。

主泵摩擦阻力矩mp2表達(dá)式為：

式中，k3為常數(shù)，由泵的實(shí)驗(yàn)值確定；n為經(jīng)驗(yàn)值，一般取1.7-2.4，由實(shí)驗(yàn)值確定。

冷卻劑質(zhì)量流量正比于主泵轉(zhuǎn)子角速度，但由于系統(tǒng)慣性，冷卻劑流量變化滯后于轉(zhuǎn)子角速度變化，將冷卻劑質(zhì)量流量方程表示為

式中，為轉(zhuǎn)子角速度標(biāo)幺值，d^*為主泵質(zhì)量流量標(biāo)幺值，τ為常數(shù)，可取值0.2s。

因此冷卻劑主泵流量計(jì)算模型為：

步驟s2：基于冷卻劑主泵設(shè)計(jì)參數(shù)及測試數(shù)據(jù)，辨識得到主泵模型參數(shù)，并計(jì)算主泵電磁轉(zhuǎn)矩、水力轉(zhuǎn)矩和摩擦轉(zhuǎn)矩。

冷卻劑屏蔽主泵模型需辨識的參數(shù)為：主泵電磁轉(zhuǎn)矩相關(guān)參數(shù)k1，k2，主泵摩擦阻力矩相關(guān)參數(shù)k3、n；并由此獲得主泵電磁轉(zhuǎn)矩，主泵水力阻矩和主泵摩擦轉(zhuǎn)矩方程。本實(shí)施例中，對ap1000核電機(jī)組冷卻劑屏蔽主泵建立的動(dòng)態(tài)模型進(jìn)行參數(shù)識別，并運(yùn)用辨識出的參數(shù)進(jìn)行仿真分析運(yùn)算。

s2.1基于冷卻劑屏蔽主泵惰轉(zhuǎn)工況，對(4)式進(jìn)行參數(shù)識別。采用采用matlab優(yōu)化設(shè)計(jì)工具箱，進(jìn)行參數(shù)辨識，辨識方程為：

式(7)中ip＝931kg·m²,ωp初始值為ωp0＝183.26rad/s，本實(shí)施例對摩擦阻力矩參數(shù)k3、n擬合結(jié)果如圖3、4所示，可以看出k3在第2次擬合即開始收斂，收斂值為0.0617，n在第12次擬合開始收斂，收斂值為1.86。

s2.2將步驟s2.1獲得的摩擦阻力矩參數(shù)k3、n代入公式(1)，得到

本實(shí)施例基于已有的冷卻劑屏蔽主泵測試數(shù)據(jù)進(jìn)行參數(shù)辨識，辨識結(jié)果如圖5所示，在第30次擬合后，k1，k2即開始收斂，收斂值分別為k1＝3059.4，k2＝245.3。

步驟s3：模型仿真與校驗(yàn)

本實(shí)施例基于所建立的數(shù)學(xué)模型，本實(shí)施例在simulink工具箱中搭建仿真模型，如圖6所示。通過不同工況下冷卻劑流量特性曲線分析驗(yàn)證本模型的有效性和準(zhǔn)確性。

s3.1惰轉(zhuǎn)曲線分析。本實(shí)施例給出了t＝10s時(shí)刻冷卻劑屏蔽主泵失去外電源后惰轉(zhuǎn)工況下本發(fā)明模型仿真曲線與主泵設(shè)計(jì)曲線的對比驗(yàn)證，如圖7所示。對比驗(yàn)證結(jié)果表明本發(fā)明提出的模型仿真得到的主泵惰轉(zhuǎn)曲線與設(shè)計(jì)惰轉(zhuǎn)曲線基本一致，驗(yàn)證了本發(fā)明提出模型的有效性。

s3.2電網(wǎng)側(cè)極端事故下的仿真驗(yàn)證。本實(shí)施例給出了冷卻劑屏蔽主泵供電頻率以4hz/s下降時(shí)，根據(jù)本發(fā)明提出模型仿真獲得的標(biāo)幺化冷卻劑流量和主泵轉(zhuǎn)速變化曲線，如圖8所示，本實(shí)施例仿真結(jié)果表明，主泵轉(zhuǎn)速下降和環(huán)路流量下降是同步的，同時(shí)主泵下降轉(zhuǎn)速超前于冷卻劑流量下降，可以根據(jù)仿真結(jié)果制定主泵低轉(zhuǎn)速停堆整定值的動(dòng)作優(yōu)先級高于冷卻劑低流量停堆整定。

s3.3頻率、電壓波動(dòng)下的仿真驗(yàn)證。本實(shí)施例給出了冷卻劑屏蔽主泵供電電壓階躍下降到0.95pu時(shí)冷卻劑流量變化，如圖9所示。其中圖9(a)為供電電壓階躍下降到0.95pu時(shí)冷卻劑流量變化曲線，圖9(b)為電網(wǎng)頻率階躍下降到0.95pu時(shí)冷卻劑流量變化曲線。由本實(shí)施例仿真結(jié)果可以看出頻率波動(dòng)對冷卻劑流量影響更為明顯。

步驟s4：接入整體模型的仿真。

將搭建的模型接入已有的核電機(jī)組，并接入電力系統(tǒng)仿真軟件進(jìn)行全時(shí)域仿真，可以實(shí)時(shí)得到在電網(wǎng)側(cè)故障時(shí)，核電側(cè)的冷卻劑流量，堆芯燃料溫度等輸出，進(jìn)而根據(jù)結(jié)果作用于堆芯停堆等保護(hù)措施制定。

本實(shí)施例給出了在t＝2s時(shí)電網(wǎng)側(cè)頻率階躍下降到0.95pu后，核電機(jī)組相關(guān)參數(shù)變化曲線，如圖10所示。圖10(a)為冷卻劑流量變化，圖10(b)為中子通量變化，圖10(c)為蒸汽發(fā)生器蒸汽壓力變化，圖10(d)為冷卻劑堆芯燃料溫度變化。

從圖10中可以看出，當(dāng)t＝2s時(shí)電網(wǎng)側(cè)頻率階躍下降到0.95pu時(shí)，冷卻劑流量迅速降低，使得冷卻劑帶著堆芯內(nèi)部溫度減小，堆芯燃料溫度上升，同時(shí)蒸汽壓力降低，反應(yīng)堆功率增加，中子通量上升。

以上是本發(fā)明的較佳實(shí)施例，凡依本發(fā)明技術(shù)方案所作的改變，所產(chǎn)生的功能作用未超出本發(fā)明技術(shù)方案的范圍時(shí)，均屬于本發(fā)明的保護(hù)范圍。