CN101399091B - 一种用于在线监测核反应堆堆芯中子通量分布的方法 - Google Patents
一种用于在线监测核反应堆堆芯中子通量分布的方法 Download PDFInfo
- Publication number
- CN101399091B CN101399091B CN2008102321942A CN200810232194A CN101399091B CN 101399091 B CN101399091 B CN 101399091B CN 2008102321942 A CN2008102321942 A CN 2008102321942A CN 200810232194 A CN200810232194 A CN 200810232194A CN 101399091 B CN101399091 B CN 101399091B
- Authority
- CN
- China
- Prior art keywords
- centerdot
- reactor
- sigma
- neutron
- reactor core
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Abstract
Description
技术领域
本发明涉及核反应堆堆芯功率监测领域,特别涉及一种用于在线监测核反应堆堆芯中子通量分布的方法。
背景技术
为了保证核反应堆的安全运行,需要监测反应堆在运行过程中堆芯的功率分布。以压水堆为例,一般通过监测核反应堆堆内中子通量分布来监测堆芯功率分布。
现有的压水堆中子测量系统分为移动式堆内中子测量系统和堆外中子测量系统。
移动式堆内中子测量系统是间歇式工作的,其工作频率一般为每月一次。移动式堆内中子测量系统的探测器为小型裂变室,系统通过指令将5组探测器快速送入堆芯底部,然后慢速升到堆芯顶部,利用反应堆堆芯中子与探测器中的UO2发生裂变使探测器内的气体电离,从而得到反应堆堆芯中子信息。该系统能够提供准确的反应堆堆内中子通量分布,但其工作间隔时间太长,不能连续监测堆内中子通量分布。
堆外中子测量系统利用布置在堆芯外部的堆外中子探测器监测堆芯功率。该方法采用点堆模型,认为堆内中子通量分布形状不随时间变化,总功率与堆内任一点的中子通量成正比,由某一点的中子通量探测值即可得到反应堆的功率。堆外监测系统中,利用堆外中子探测器的中子通量探测值监测反应堆功率。该方法简单易行,但由于假定堆内中子通量分布形状不随时间变化,因此不能提供堆内中子通量的动态分布信息。
西屋公司开发的BEACON(Best Estimate Analyzer for Core Operation-NUclear)系统可在线监测堆芯功率分布,它使用二维模型计算堆内中子通量分布,使用堆外探测器提供轴向通量分布。
发明内容
本发明的目的在于提供一种在线监测核反应堆堆芯通量分布的方法,它能够根据反应堆堆内中子探测器及布置在堆芯外部的堆外探测器测量得到的中子通量信息,快速计算出反应堆堆芯的三维中子通量分布,从而给出反应堆堆芯通量分布的实时信息,为反应堆堆芯的控制和保护系统提供参考。
为了实现上述目的,本发明采用以下技术方案予以实现:一种用于在线监
测核反应堆堆芯中子通量分布的方法,基于在反应堆布置M1个堆内中子探测器和M2个堆外中子探测器,M1个堆内中子探测器的读数分别记为 其中是以反应堆几何中心为原点的极坐标值,分别代表M1个堆内中子探测器的位置,M2个堆外中子探测器的读数分别记为 该方法包括以下步骤:
首先,根据反应堆的运行状态和控制历史,更新反应堆的参考堆芯模型,确定参考堆芯模型的多群中子扩散方程矩阵形式中的参量M、F,其中k为反应堆有效增值系数,Φ为堆芯中子通量,不同的特征值kn对应不同的特征向量Φn,有:
求解多群中子扩散方程,并根据kn由大到小的顺序进行排序Φn(n=1,…,∞);其中,
r为以反应堆几何中心为原点的极坐标值;
g=1,2,…,G;
G为总的能群数;
Dg(r)为第g能群的扩散系数(1/cm);
∑tg(r)为第g能群的宏观总截面(1/cm);
∑g′→g(r)为从g′能群散射到g能群宏观散射截面(1/cm);
xg为g能群中子裂变份额;
v∑fg(r)为g能群宏观ν-裂变截面(1/cm);
k为反应堆有效增值系数;
Φg(r)为g能群堆芯中子通量;
其次,采集堆内、堆外中子探测器的读数,并建立方程组:
其中,分别为M1个堆内中子探测器的响应函数,为M2个堆外中子探测器空间响应函数,V为反应堆堆芯空间的集合,Φi(r)为n=i特征向量Φn,即参考堆芯的第i阶谐波,ai(i=1~N)为展开系数,N为展开阶数,(M1+M2)≥N,r为以反应堆几何中心为原点的极坐标值,然后根据上述方程组,求解展开系数ai(i=1~N);
本发明的进一步特点在于:
所述求解多群中子扩散方程采用Krylov子空间方法。
所述求解展开系数ai(i=1~N)采用最小二乘法。
所述M1个堆内中子探测器在反应堆中对称布置,最好是中心对称布置。
所述M2个堆外中子探测器环绕反应堆对称布置,最好是中心对称布置。
由于本发明可利用现有的反应堆监测设备,不必增加冗余的设备,即可对核反应堆堆内三维中子通量分布的在线监测,实现方便,成本低。本发明只需计算有限的几个展开系数ai(i=1~N),并在求解多群中子扩散方程时采用krylov子空间方法,求解展开系数时采用最小二乘法,计算量小;堆内、堆外中子探测器相对反应堆对称布置,尤其是中心对称布置,计算简化;因此,采用本发明计算速度快,在线监测实时性好。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明。
图1为核反应堆堆芯中子通量分布在线监测系统;
图2为某核电站堆内中子探测器及堆外中子探测器径向布置图;
图3为某核电站堆芯中子通量分布重构轴向FRMS分布图。
具体实施方式
本发明的思想是将真实的核反应堆堆芯内中子通量分布用参考堆芯模型的高阶谐波线性展开,根据布置在堆芯内、外的中子探测器的读数得到展开系数,继而在线重构出真实的反应堆堆芯内中子通量分布,通过可视化系统显示反应堆堆芯内中子通量分布,为反应堆堆芯控制和保护提供参考。
下面以压水堆为例对本发明进行进一步说明。
本发明基于参考堆芯模型,为了便于描述,首先需要说明表征参考堆芯模型的反应堆多群中子扩散方程。
反应堆中子扩散方程是关于中子通量分布的方程,中子通量分布是空间变量r和中子能量E的函数。在反应堆物理计算中,通常对中子能量E采用“分群”的近似方法,即把中子能量区域按能量的大小划分成G个能量区间,每一个能量区间称为一个能群。G取大于1时的中子扩散方程称为多群中子扩散方程。
反应堆多群中子扩散方程为:
式中:
r为以反应堆几何中心为原点的极坐标值;
g=1,2,…,G;
G为总的能群数;
Dg(r)为第g能群的扩散系数(1/cm);
∑tg(r)为第g能群的宏观总截面(1/cm);
∑g′→g(r)为从g′能群散射到g能群宏观散射截面(1/cm);
xg为g能群中子裂变份额;
v∑fg(r)为g能群宏观ν-裂变截面(1/cm);
k为反应堆有效增值系数;
Φg(r)为g能群堆芯中子通量。
反应堆多群中子扩散方程(1)可写成矩阵形式:
其中,
方程(2)是一个特征值问题,不同的特征值对应不同的特征向量。即:
方程(3)的解Φn(n=1,…,∞)即为中子扩散方程(1)的高阶谐波。根据kn由大到小的顺序对谐波进行排序。
因此,本发明的的第一步,需要根据反应堆的运行状态和控制历史,更新反应堆的参考堆芯模型,确定参考堆芯模型的多群中子扩散方程矩阵形式中的参量M、F,其中k为反应堆有效增值系数,Φ为堆芯中子通量,上式中不同的特征值kn对应不同的特征向量Φn,有:
求解多群中子扩散方程,并根据kn由大到小的顺序进行排序Φn(n=1,…,∞)。参考堆芯模型的建立,需要根据反应堆的运行状态和控制历史,对参考堆芯的截面数据进行更新,使参考堆芯最大程度的接近实际堆芯,以提高对真实堆芯中子通量分布的重构精度。在反应堆运行和控制过程中,能够对堆芯截面数据产生比较大的影响的因素主要有三种:第一,控制棒的移动;第二,硼浓度的调节;第三,燃耗积累。本发明的参考堆芯模型的建立,可以采用传统的方式。
本发明采用Krylov子空间方法对方程(2)进行求解,这种方法与传统的方法相比,计算量小,实时性好,便于在线监测。krylov子空间方法是一种基于投影的迭代方法,它的主要思想是:将大规模矩阵方程的解投影到一个较小规模的方程上,通过求解小规模矩阵方程,得到大规模矩阵方程的解。其基本原理如下:
求解大型矩阵的特征值问题Ax=λx,其中A∈Rnxn,取m(m《n)维子空间Vm=Span{v1,Av1,…,Am-1v1},限制在子空间Vm中计算A的“近似”特征向量xm,则xm表示为:xm=Vmym。将xm的表达式代入Ax=λx可得:
式(4)为m×m的特征值问题。一旦成功计算出规模较小的特征值问题的特征值λm和特征向量ym,就可得到规模较大的矩阵A“近似”的特征值λm和特征向量xm=Vmym。方程(4)的特征值λm和特征向量ym一般通过QR分解得到。
由中子扩散方程的矩阵形式(2)经变换后得到:M-1FΦ=kΦ,给定初始向量w1,依次求解Mw2=Fw1,Mw3=Fw2…Mwm=Fwm-1,则Span{w1,w2,…,wm}构成矩阵M-1F的Krylov子空间,然后根据Krylov子空间方法的原理求解中子扩散方程的特征值kn和特征向量Φn(n=1,2,…)。
本发明的核心是在线重构出真实的反应堆堆芯内中子通量分布,其原理基于多群中子扩散方程高阶谐波具有完备性,即:反应堆中,凡是满足反应堆中子扩散方程边界条件的任意的通量分布均可由中子扩散方程的高阶谐波线性展开。因此真实的反应堆堆芯中子通量分布可以用参考堆芯的高阶谐波线性展开成如下形式:
(其中,Φreal(r)为真实的反应堆堆芯中子通量分布,Φi(r)为参考堆芯的第i阶谐波,ai为展开系数)。
具体的实施步骤如下:
首先,在反应堆布置M1个堆内中子探测器和M2个堆外中子探测器,M1个堆内中子探测器的读数分别记为其中是以反应堆几何中心为原点的极坐标值,分别代表M1个堆内中子探测器的位置,,M2个堆外中子探测器的读数分别记为
其次,采集堆内、堆外中子探测器的读数,并建立方程组:
其中,分别为M1个堆内中子探测器的响应函数,为M2个堆外中子探测器空间响应函数,V为反应堆堆芯空间的集合,Φi(r)为n=i特征向量Φn,即参考堆芯的第i阶谐波,ai(i=1~N)为展开系数,N为展开阶数,(M1+M2)≥N,r为以反应堆几何中心为原点的极坐标值,然后根据上述方程组,求解展开系数ai(i=1~N)。
参照图1,核反应堆堆芯中子通量分布在线监测系统,由堆芯中子通量重构子系统6和谐波更新子系统7。
堆芯中子通量重构子系统6的功能是:利用堆芯信息采集子系统1,根据堆内中子探测器及堆外中子探测器的探测值在线重构堆芯中子通量分布,并实时显示重构得到的堆芯中子通量分布。堆芯中子通量重构子系统6由两部分构成,一为通量重构单元2,一为堆芯通量分布实时显示单元3。
堆芯通量分布实时显示单元3用于实时显示通量重构单元2重构得到的真实堆芯中子通量分布,为反应堆操控人员提供堆芯通量分布信息并为反应堆控制提供参考。
对于局部扰动(如控制棒移动)引起的堆芯中子通量分布畸变,由于谐波本身的特点,畸变后的堆芯中子通量分布用参考堆芯谐波展开,N需要取的很大。由于探测器的个数(M1+M2)有限,通常无法满足(M1+M2)≥N的条件,因此无法对局部扰动引起畸变后的堆芯中子通量分布进行精确重构。同时,当探测器个数一定时,参考堆芯与实际堆芯越接近,堆芯中子通量分布重构的精度越高。为了解决无法精确重构局部扰动引起的通量畸变问题以及为了提高堆芯中子通量分布的重构精度,本发明采用更新谐波的方法,即更新反应堆的参考堆芯模型,计算更新谐波,由谐波更新子系统7完成。
谐波更新子系统7主要用于计算并向堆芯中子通量重构系统6提供参考堆芯高阶谐波。谐波更新子系统7由两部分组成,一为参考堆芯更新单元5,一为参考堆芯谐波计算单元4。
参考堆芯更新单元5主要分三种情况根据堆芯截面参数来更新参考堆芯模型:第一,控制棒的移动;第二,硼浓度的调节;第三,燃耗积累。
对于控制棒的移动,参考堆芯更新单元5根据控制棒驱动机构提供的棒位信息,获得控制棒移动前后的棒位差。由于控制棒移动单位距离对堆芯截面参数的影响是已知的,因此根据控制棒移动前的堆芯截面数据及控制棒移动前后的棒位差,即可得到控制棒移动后的堆芯截面数据,完成参考堆芯模型的更新。
对于硼浓度的调节,同样由于硼浓度的变化对堆芯截面参数的影响是已知的。根据硼浓度调节前的堆芯截面数据及调节前后的硼浓度差值,得到硼浓度调节后的堆芯截面参数,完成参考堆芯模型的更新。
对于燃耗积累引起的堆芯截面参数变化,可由现有的组件计算程序(如CASMO、APOLLO、TPFAP等)对燃耗后的堆芯截面参数进行更新。由于燃耗积累对堆芯截面参数变化的影响是一个缓慢的过程,因此间隔一段时间(例如10天)进行一次燃耗计算,燃耗计算过程中需要用到堆芯的中子通量分布信息,该信息由通量重构单元2提供。
参考堆芯谐波计算单元4的功能是根据参考堆芯更新单元5提供的参考堆芯截面参数信息,完成对参考堆芯高阶谐波的计算,并将计算得到的高阶谐波提供给堆芯中子通量重构子系统6用于重构实际堆芯的中子通量分布。
参照图2,为某核电站堆内中子探测器及堆外中子探测器径向布置图。反应堆堆外布置有八个中子测量通道,分别为两个源量程测量通道9、两个中间量程测量通道10、和四个功率量程测量通道8。源量程测量通道9提供反应堆初始启动时的中子通量测量,中间量程测量通道10提供冗余的中子通量测量。在径向上,14个组件11相对反应堆中心对称分布,每个组件在轴向上均匀布置4个堆内中子探测器,所有堆内中子探测器以中心对称分布在反应堆中;堆外探测器布置在功率测量通道8中,共计16个,16个堆外探测器环绕反应堆中心对称分布。
本发明以上述核电站为例,进行了仿真计算。计算采用2群模型(即:令方程(1)中的G=2),定义第一群为快群,第二群为热群。实际反应堆功率是由热群中子与235U发生裂变反应产生的,因此人们更加关注堆芯内热群中子通量分布。以下计算只对堆芯热群中子通量分布进行重构。
为了验证利用谐波重构堆芯通量分布的精度,进行如下试验:以上述核电站装料后满功率运行12.71天后的堆芯为参考堆芯,对经过满功率运行16.22天后的堆芯中子通量分布进行重构。满功率运行12.71天后的堆芯中子通量分布用Φ12.71表示,满功率运行16.22天后的堆芯中子通量分布用Φ16.22表示,经过重构得到的堆芯中子通量分布用Φrec表示。堆芯轴向划分为23区,对于第Z区,定义 (Z=1~23)。其中,RMS{Φ16.22(Z)-Φrec(Z)}表示{Φ16.22(Z)-Φrec(Z)}的均方根,表示该层的平均通量;同理,定义
选取轴向的四个区(分别为Z=23、Z=19、Z=16、Z=12),各区1/4堆芯的径向误差分布分别如表1、表2、表3、表4所示。
表1
表2
表3
表4
参照图3,为FRMSrec(Z)与FRMS0(Z)沿轴向的分布,其中横轴代表轴向分区,纵轴代表FRMSrec(Z)与FRMS0(Z)的值。该图说明,在反应堆真实堆芯中子通量分布相对于参考堆芯中子通量分布变化很大的情况下,依据本发明的方法重构得到的堆芯中子通量分布相对于真实堆芯中子通量分布具有很小的误差。
试验计算结果表明,依据本发明的方法重构得到的堆芯中子通量分布相对于真实堆芯中子通量分布的误差均在3.5%以下,具有很好的精度。
本发明可以应用到任意类型的反应堆上,将该型反应堆诱发裂变的中子能群的通量分布用参考堆芯的高阶谐波线性展开,则只要能够提供该型反应堆诱发裂变的中子能群在探测点的值,就可依据探测值求解展开系数,从而得到堆芯中子通量分布。
Claims (6)
1.一种用于在线监测核反应堆堆芯中子通量分布的方法,基于在反应堆布置M1个堆内中子探测器和M2个堆外中子探测器,M1个堆内中子探测器的读数分别记为其中是以反应堆几何中心为原点的极坐标值,分别代表M1个堆内中子探测器的位置,M2个堆外中子探测器的读数分别记为该方法包括以下步骤:
首先,根据反应堆的运行状态和控制历史,更新反应堆的参考堆芯模型,确定参考堆芯模型的多群中子扩散方程矩阵形式中的参量M、F,其中k为反应堆有效增值系数,Φ为堆芯中子通量,不同的特征值kn对应不同的特征向量Φn,有:
求解多群中子扩散方程,并根据kn由大到小的顺序进行排序Φn(n=1,…,∞);
其中,
r为以反应堆几何中心为原点的极坐标值;
g=1,2,…,G;
G为总的能群数;
Dg(r)为第g能群的扩散系数(1/cm);
∑tg(r)为第g能群的宏观总截面(1/cm);
∑g′→g(r)为从g′能群散射到g能群宏观散射截面(1/cm);
χg为g能群中子裂变份额;
v∑fg(r)为g能群宏观v-裂变截面(1/cm);
k为反应堆有效增值系数;
其次,采集堆内、堆外中子探测器的读数,并建立方程组:
其中,分别为M1个堆内中子探测器的响应函数,为M2个堆外中子探测器空间响应函数,V为反应堆堆芯空间的集合,Φi(r)为n=i特征向量Φn,即参考堆芯的第i阶谐波,ai(i=1~N)为展开系数,N为展开阶数,(M1+M2)≥N,r为以反应堆几何中心为原点的极坐标值,然后根据上述方程组,求解展开系数ai(i=1~N);
2.根据权利要求1所述的一种用于在线监测核反应堆堆芯中子通量分布的方法,其特征在于,所述求解多群中子扩散方程采用Krylov子空间方法。
3.根据权利要求1所述的一种用于在线监测核反应堆堆芯中子通量分布的方法,其特征在于,所述求解展开系数ai(i=1~N)采用最小二乘法。
4.根据权利要求1所述的一种用于在线监测核反应堆堆芯中子通量分布的方法,其特征在于,所述M1个堆内中子探测器在反应堆中对称布置。
5.根据权利要求1所述的一种用于在线监测核反应堆堆芯中子通量分布的方法,其特征在于,所述M2个堆外中子探测器环绕反应堆对称布置。
6.根据权利要求4、或5所述的一种用于在线监测核反应堆堆芯中子通量分布的方法,其特征在于,所述对称布置是中心对称布置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102321942A CN101399091B (zh) | 2008-11-07 | 2008-11-07 | 一种用于在线监测核反应堆堆芯中子通量分布的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102321942A CN101399091B (zh) | 2008-11-07 | 2008-11-07 | 一种用于在线监测核反应堆堆芯中子通量分布的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101399091A CN101399091A (zh) | 2009-04-01 |
CN101399091B true CN101399091B (zh) | 2012-02-01 |
Family
ID=40517563
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008102321942A Expired - Fee Related CN101399091B (zh) | 2008-11-07 | 2008-11-07 | 一种用于在线监测核反应堆堆芯中子通量分布的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101399091B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10438708B2 (en) * | 2011-10-04 | 2019-10-08 | Westinghouse Electric Company Llc | In-core instrument thimble assembly |
CN103294899A (zh) * | 2013-05-10 | 2013-09-11 | 西安交通大学 | 一种计算小型实验反应堆堆芯中子通量分布的方法 |
CN103294898A (zh) * | 2013-05-10 | 2013-09-11 | 西安交通大学 | 一种计算反应堆全堆芯单棒功率的方法 |
CN103794256B (zh) * | 2013-11-19 | 2015-12-09 | 国核(北京)科学技术研究院有限公司 | 反应堆异常监测方法及系统 |
CN105280253B (zh) * | 2013-11-19 | 2018-12-18 | 国核(北京)科学技术研究院有限公司 | 操作反应堆堆芯功率的方法及系统 |
CN104678424B (zh) * | 2013-12-03 | 2017-08-25 | 中广核(北京)仿真技术有限公司 | 核电站模拟机的中子计数系统及方法 |
CN103943158B (zh) * | 2013-12-31 | 2016-06-29 | 西安交通大学 | 一种消除自给能中子探测器延迟效应的方法 |
CN104931798B (zh) * | 2015-04-29 | 2017-12-15 | 中广核核电运营有限公司 | 堆外核仪表系统的源量程脉冲测量装置及其检测定位方法 |
CN105445780B (zh) * | 2015-11-10 | 2018-11-30 | 中国原子能科学研究院 | 一种中子单色器高次谐波的测量装置及方法 |
CN106782687B (zh) * | 2015-11-20 | 2023-11-28 | 国家电投集团科学技术研究院有限公司 | 装载193盒燃料组件的堆芯的探测器布置方法 |
CN106202611B (zh) * | 2016-06-24 | 2019-04-09 | 西安交通大学 | 一种适用于核反应堆组件物理热工耦合的在线计算方法 |
CN106024078B (zh) * | 2016-06-24 | 2017-08-25 | 西安交通大学 | 一种诊断反应堆中子探测器失效的方法 |
CN106024080B (zh) * | 2016-06-24 | 2017-07-28 | 西安交通大学 | 一种获取反应堆堆芯中子通量密度精细分布的方法 |
CN107644133B (zh) * | 2017-09-19 | 2021-05-14 | 中国核动力研究设计院 | 一种核反应堆堆芯吊篮振动刻度因子计算方法 |
CN108447574B (zh) * | 2018-03-28 | 2020-08-11 | 中广核工程有限公司 | 一种核电站堆外探测器中子计数率计算方法及其系统 |
CN109767854B (zh) * | 2019-01-21 | 2022-07-08 | 中国科学院合肥物质科学研究院 | 一种基于堆外测量数据的反应堆内中子三维分布测量系统 |
CN109859867B (zh) * | 2019-01-29 | 2023-03-24 | 中山大学 | 一种基于本征正交分解的堆芯三维中子通量实时重构方法 |
CN109830317B (zh) * | 2019-01-29 | 2023-03-24 | 中山大学 | 一种基于跟踪计算的堆芯功率本征正交分解在线重构方法 |
CN111584104B (zh) * | 2020-05-08 | 2022-03-22 | 西安交通大学 | 基于多次碰撞源-蒙特卡罗耦合获取反应堆外探测器响应的方法 |
CN111584019B (zh) * | 2020-05-08 | 2022-12-09 | 西安交通大学 | 基于首次碰撞源-蒙特卡罗耦合获取反应堆外探测器响应的方法 |
EP4174875A1 (en) | 2021-10-27 | 2023-05-03 | Framatome | Computer implemented method for simulating an operation of a reactor core |
WO2023141741A1 (zh) * | 2022-01-25 | 2023-08-03 | 中广核研究院有限公司 | 一种堆芯功率在线监测方法、系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1013321B (zh) * | 1987-04-24 | 1991-07-24 | 法国原子能委员会 | 用于核反应堆的控制方法 |
EP1571675A2 (en) * | 2004-01-09 | 2005-09-07 | Korea Electric Power Corporation | Digital control system and method for neutron flux mapping systems |
-
2008
- 2008-11-07 CN CN2008102321942A patent/CN101399091B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1013321B (zh) * | 1987-04-24 | 1991-07-24 | 法国原子能委员会 | 用于核反应堆的控制方法 |
EP1571675A2 (en) * | 2004-01-09 | 2005-09-07 | Korea Electric Power Corporation | Digital control system and method for neutron flux mapping systems |
Non-Patent Citations (4)
Title |
---|
JP特开平10-10275A 1998.01.16 |
叶青,吴宏春.R-Z几何中子输运方程的球谐函数方法.《核动力工程》.2008,第29卷(第4期),19-23. * |
曹良志,吴宏春,姚栋.用于复杂环境下中子时空动力学计算的简化球谐函数方法.《核动力工程》.2007,第28卷(第1期),18-21. * |
曹良志,吴宏春.非结构网格中子输运方程的球谐函数解法研究.《核动力工程》.2004,第25卷(第5期),395-398. * |
Also Published As
Publication number | Publication date |
---|---|
CN101399091A (zh) | 2009-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101399091B (zh) | 一种用于在线监测核反应堆堆芯中子通量分布的方法 | |
Li et al. | Short-term wind power prediction based on extreme learning machine with error correction | |
CN102663251B (zh) | 基于计算流体力学模型的风电场功率物理预测方法 | |
Li et al. | Hydraulic turbine governing system identification using T–S fuzzy model optimized by chaotic gravitational search algorithm | |
CN102792542A (zh) | 发电量预测装置和其方法及程序 | |
CN103116711B (zh) | 基于辐照度特征参数的天气类型辨识方法 | |
CN103700036B (zh) | 一种适于电力系统多时间尺度的暂态稳定性投影积分方法 | |
CN103617817A (zh) | 监测反应堆堆芯功率的方法及系统 | |
Nigussie et al. | Design, modeling, and CFD analysis of a micro hydro Pelton turbine runner: for the case of selected site in Ethiopia | |
CN107229824A (zh) | 光伏电站发电单元功率曲线建模方法及装置 | |
CN102790389B (zh) | 一种暂态稳定裕度概率分布的计算方法 | |
CN104156542A (zh) | 一种基于隐式投影的有源配电系统稳定性仿真方法 | |
CN113409975A (zh) | 基于模型降阶和数据同化的堆芯功率分布监测方法及系统 | |
CN105846425A (zh) | 一种基于通用风电预测误差模型的经济调度方法 | |
Han et al. | A novel wind farm equivalent model for high voltage ride through analysis based on multi-view incremental transfer clustering | |
Wang et al. | On-line reconstruction of in-core power distribution by harmonics expansion method | |
CN109102102A (zh) | 基于多变量相空间重构和svr的光伏出力短期预测方法 | |
CN102545218A (zh) | 基于电能质量监测系统的在线负荷建模并行计算方法 | |
Xia et al. | Performance evaluation of a 3-D kinetic model for CANDU reactors in a closed-loop environment | |
CN107276093A (zh) | 基于场景削减的电力系统概率潮流计算方法 | |
Li et al. | Proper orthogonal decomposition based online power-distribution reconstruction method | |
Lukhyswara et al. | Solar irradiation forecasting uses time series analysis | |
Zhou et al. | PV power characteristic modeling based on multi-scale clustering and its application in generation prediction | |
Collins et al. | Coupled Simulations of Fast-Spectrum Molten Salt Reactors | |
Nissan | Intelligent technologies for nuclear power systems: heuristic and neural tools |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120201 Termination date: 20211107 |