CN115295179B - 反应堆功率测量的补偿方法 - Google Patents

反应堆功率测量的补偿方法 Download PDF

Info

Publication number
CN115295179B
CN115295179B CN202211004907.6A CN202211004907A CN115295179B CN 115295179 B CN115295179 B CN 115295179B CN 202211004907 A CN202211004907 A CN 202211004907A CN 115295179 B CN115295179 B CN 115295179B
Authority
CN
China
Prior art keywords
detector
compensation
power
reactor
function
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.)
Active
Application number
CN202211004907.6A
Other languages
English (en)
Other versions
CN115295179A (zh
Inventor
刘勇
王振忠
贾晓淳
贺长兴
胡赟
王晓坤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute of Atomic of Energy
Original Assignee
China Institute of Atomic of Energy
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Institute of Atomic of Energy filed Critical China Institute of Atomic of Energy
Priority to CN202211004907.6A priority Critical patent/CN115295179B/zh
Publication of CN115295179A publication Critical patent/CN115295179A/zh
Application granted granted Critical
Publication of CN115295179B publication Critical patent/CN115295179B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21C—NUCLEAR REACTORS
    • G21C7/00—Control of nuclear reaction
    • G21C7/28—Control of nuclear reaction by displacement of the reflector or parts thereof
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Plasma & Fusion (AREA)
  • General Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明的实施例公开了一种反应堆功率测量的补偿方法。所述反应堆的堆芯外围设置有可移动的反射层，通过所述反射层的位移来控制所述反应堆的反应性，所述反射层外设置有多个核探测器。所述补偿方法用于补偿反射层位移对反应堆功率测量值造成的偏差，其包括：获取所述多个核探测器中工作组探测器和补偿组探测器的功率测量值；根据所述工作组探测器和补偿组探测器的功率测量值之间的比值，确定反射层位置；根据反射层位置确定所述工作组探测器的功率补偿系数；利用所述功率补偿系数来补偿所述工作组探测器的功率测量值，得到所述反应堆的当前功率值。

Description

反应堆功率测量的补偿方法
技术领域
本发明的实施例涉及技术领域，具体涉及一种反应堆功率测量的补偿方法。
背景技术
在反应堆中，反射层可用于控制反应堆的反应性，其需要通过位移改变自身遮挡堆芯的程度来实现对反应性的控制。然而，当反射层上下移动时，空间上的中子注量率水平和分布均会剧烈变化，并且不同位置处中子注量率的变化规律具有显著差异。
反应堆的功率一般是采用核探测器测量的，核探测器通过测量反应堆的中子注量率来间接测量反应堆功率。对于通过反射层控制反应性的特定反应堆，核探测器布置在反射层之外，反射层的位移会明显地影响核探测器对反应堆功率的测量结果，造成核探测器测量得到的功率值偏离真实的反应堆功率。
发明内容
根据本发明的一个方面，提供了一种反应堆功率测量的补偿方法。所述反应堆的堆芯外围设置有可移动的反射层，所述反射层外设置有多个核探测器，通过所述反射层的位移来控制所述反应堆的反应性。所述补偿方法用于补偿反射层位移对反应堆功率测量值造成的偏差，其包括：获取所述多个核探测器中工作组探测器和补偿组探测器的功率测量值；根据所述工作组探测器和补偿组探测器的功率测量值之间的比值，确定反射层位置；根据所述反射层位置确定所述工作组探测器的功率补偿系数；利用所述功率补偿系数来补偿所述工作组探测器的功率测量值，得到所述反应堆的当前功率值。
附图说明
通过下文中参照附图对本发明的实施例所作的描述，本发明的其它目的和优点将显而易见，并可帮助对本发明有全面的理解。
图1是根据本发明一个实施例的反应堆的结构示意图。
图2是图1的反应堆的径向结构示意图。
图3是图1的反应堆的轴向结构示意图。
图4是根据本发明一个实施例反应堆中反射层运动至不同位置的示意图。
图5是根据本发明一个实施例的补偿方法的流程示意图。
图6是根据本发明另一个实施例的补偿方法的流程示意图。
图7a和图7b是根据本发明一个实施例的多个探测器位置处中子注量率的变化关系曲线的示意图。
图8是根据图7b中变化关系曲线G(x)确定的函数f1示意图。
图9是根据图7b中变化关系曲线G(x)确定的函数f2示意图。
需要说明的是，附图并不一定按比例来绘制，而是仅以不影响读者理解的示意性方式示出。
附图标记说明：
10、堆芯；20、反射层；1～6、核探测器。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚，下面将结合本申请实施例的附图，对本申请的技术方案进行清楚、完整地描述。显然，所描述的实施例是本申请的一个实施例，而不是全部的实施例。基于所描述的本申请的实施例，本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例，都属于本申请保护的范围。
需要说明的是，除非另外定义，本申请使用的技术术语或者科学术语应当为本申请所属领域内具有一般技能的人士所理解的通常意义。若全文中涉及“第一”、“第二”等描述，则该“第一”、“第二”等描述仅用于区别类似的对象，而不能理解为指示或暗示其相对重要性、先后次序或者隐含指明所指示的技术特征的数量，应该理解为“第一”、“第二”等描述的数据在适当情况下可以互换。若全文中出现“和/或”，其含义为包括三个并列方案，以“A和/或B”为例，包括A方案，或B方案，或A和B同时满足的方案。此外，为了便于描述，在这里可以使用空间相对术语，如“上方”、“下方”、“顶部”、“底部”等，仅用来描述如图中所示的一个器件或特征与其他器件或特征的空间位置关系，应当理解为也包含除了图中所示的方位之外的在使用或操作中的不同方位。
图1示出了根据本发明一个实施例的反应堆的结构示意图。图2和图3示出了图1中反应堆的不同视角的结构示意图。
参见图1至图3，反应堆包括堆芯10、反射层20和多个核探测器。可移动的反射层20设置于堆芯10的外围且贴近堆芯10，多个核探测器设置于反射层20外，且分布于反射层20周围的不同位置。其中，核探测器1至核探测器4沿反射层20的周向方向均匀布置，核探测器5布置于堆芯10的上方，核探测器6布置于堆芯10的下方。
图4示出了根据本发明一个实施例的反射层位移至不同位置的示意图。反射层20可以在如图4所示的范围内上下位移，即，从反射层20完全脱离堆芯10到完全遮挡堆芯10的范围内位移，从而通过反射层20的位移来控制反应堆的反应性。
当反射层20运动至不同位置时，其周围空间上的中子注量率大小和分布会发生明显改变，并且反射层20的位置对不同位置处的中子注量率有着不同的影响关系。因此，由于反射层20的位置的位移造成的影响，即使在反应堆的实际功率水平不变的情况下，各个不同位置处的核探测器(1～6)测量得到的反应堆功率值均会发生变化，导致核探测器(1～6)直接测得的反应堆功率由于反射层20的控制性位移而造成明显的偏离。
本发明的实施例提供了一种反应堆功率测量的补偿方法，可以用于补偿反射层位移对反应堆功率测量值造成的偏差，得到较为真实的反应堆功率值。图5示出了根据本发明一个实施例的补偿方法的流程示意图。如图5所述，本实施例中的补偿方法具体包括以下步骤。
步骤S11，获取反应堆中多个核探测器中工作组探测器和补偿组探测器的功率测量值。
步骤S12，根据工作组探测器和补偿组探测器的功率测量值之间的比值，确定反射层位置；再根据所述反射层位置确定所述工作组探测器的功率补偿系数。
步骤S13，利用功率补偿系数来补偿工作组探测器的功率测量值，得到反应堆的当前功率值。
采用本发明实施例的补偿方法，相比于直接通过核探测器测得的反应堆功率值由于受到反射层控制性位移影响而造成明显偏差，本实施例利用工作组探测器和补偿组探测器的功率测量值，经过计算得到功率补偿系数，可以用于补偿反射层控制性位移对反应堆功率测量结果的影响，提高了核探测器测量反应堆功率的测量精度。
在步骤S11中，可以接收工作组探测器和补偿组探测器的输出信号，并对两者的输出信号进行处理得到工作组探测器的功率测量值和补偿组探测器的功率测量值。
在步骤S12中，首先需要对工作组探测器的功率测量值和补偿组探测器的功率测量值进行比例计算。具体地，两者的比值在数值上可以为工作组探测器功率测量值(Pa)除以补偿组探测器功率测量值(Pb)，即Pa/Pb；或者，可以是补偿组探测器功率测量值(Pb)除以工作组探测器功率测量值(Pa)，即Pb/Pa；或者，可以是工作组探测器功率测量值(Pa)与补偿组探测器功率测量值(Pb)进行对数运算后的差值，即|logz Pa-logz Pb|，其中，z为任意底数。
在计算得到工作组探测器的功率测量值和补偿组探测器的功率测量值之间的比值后，即可根据该比值通过特定算法来确定反射层位置，根据反射层位置通过特定算法即可确定工作组探测器的功率补偿系数。
在步骤S13中，利用计算得到的功率补偿系数来补偿工作组探测器的功率测量值(Pa)，具体地，将工作组探测器的功率测量值乘以功率补偿系数，即可得到反应堆的当前功率值。本实施例中利用功率补偿系数补偿后的功率值接近反应堆的真实功率，与反应堆的真实功率之间的偏差较小，从而提高了核探测器测量反应堆功率时的测量精度。
在一些实施中，首先需要在反应堆中选定工作组探测器以及补偿组探测器的位置。其中，工作组探测器的输出信号可以用于计算最终的反应堆当前功率值，工作组探测器和补偿组探测器的输出信号共同用于确定功率补偿系数。
可选的，工作组探测器的数量可以是一个或者多个，补偿组探测器的数量也可以是一个或多个。其中，多个工作组探测器可以共用一个补偿组探测器，也可以每个工作组探测器对应设置一个补偿组探测器。
在一些实施例中，可以根据反应堆功率恒定时反应堆中各核探测器位置处的中子注量率随反射层位置的变化关系，来选定工作组探测器以及补偿组探测器的位置。其中，在反应堆功率恒定时，例如，反应堆在某个恒定的堆功率下，各核探测器位置处的中子注量率随反射层位置的变化关系可以记为函数G(x)，G(x)的自变量为反射层位置，函数值为核探测器位置处的中子注量率。本实施例中，可以根据各核探测器的函数G(x)来选定工作组探测器以及补偿组探测器的位置。
需要说明的是，核探测器通过测量中子注量率来间接测量反应堆功率，核探测器位置处中子注量率水平与其测量得到反应堆功率测量值大小相对应。因此，获取的反应堆功率恒定时各核探测器位置处的中子注量率随反射层位置变化关系的函数G(x)，也可以反映各核探测器所测量得到的反应堆功率测量值随反射层位置的变化关系。
具体地，在选定工作组探测器和补偿组探测器的位置时需要满足以下条件：当反应堆功率恒定时，工作组探测器位置处的G(x)的函数值变化范围小于所述补偿组探测器位置处的G(x)的函数值变化范围，即，工作组探测器位置处的G(x)的函数值变化范围应尽量较小，而补偿组探测器位置处的G(x)的函数值变化范围应尽量较大。同时，补偿组探测器位置处的中子注量率随反射层的变化应满足单值函数关系，即，补偿组探测器位置处的G(x)为单值函数。此外，补偿组探测器位置处的G(x)与工作组探测器位置处的G(x)之比也应为单值函数。
本实施例中选择G(x)的函数值变化范围较小的核探测器作为工作组探测器，其受到反射层控制性位移影响较小，因此工作组探测器直接测量得到的功率测量值的偏差也较小，在此基础上进行功率补偿，有利于进一步提高精度，得到更加精确的反应堆功率。
此外，本实施例中选择G(x)的函数变化范围较大的核探测器作为补偿组探测器，可以在后续计算补偿组探测器与工作组探测器的功率测量值的比值时，提高计算得到的比值的信噪比，防止补偿组探测器与工作组探测器的G(x)的函数值变化范围均较小而造成比值的信噪比较差。同时，本实施例中选定的补偿组探测器的G(x)为单值函数，并且补偿组探测器位置处的G(x)与工作组探测器位置处的G(x)之比也应为单值函数，可以避免采用探测器位置处的G(x)为多值函数或者工作组探测器与补偿组探测器位置处的G(x)比值为多值函数而导致后续无法计算得到函数f1和f2。
需要说明的是，本发明实施例中所述的“G(x)函数值变化范围”是指G(x)函数的最大函数值与最小函数值的差值。
在一些实施例中，当反应堆的堆芯外围设置有多个可移动的反射层时，选定工作组探测器的位置时还需要考虑核探测器在方向上均匀性。具体地，可以根据各反射层对各核探测器位置处中子注量率的影响关系，来选定工作组探测器的位置，其中，各反射层对工作组探测器位置处的中子注量率的影响关系应尽可能一致，从而在不同的反射层分别位移时，均能够通过相同的算法来计算补偿系数，可以在减小计算量的同时得到接近实际的反应堆功率值，并且在多个反射层一起位移时也能够稳定计算得到准确的补偿系数，进而得到较为精确的反应堆功率值。
同样地，在一些实施例中，当反应堆的堆芯外围设置有多个可移动的反射层时，选定补偿组探测器的位置时也需要考虑核探测器在方向上均匀性。具体地，可以根据各反射层对各核探测器位置处中子注量率的影响关系，来选定补偿组探测器的位置，其中，各反射层对补偿组探测器位置处的中子注量率的影响关系应尽可能一致，从而在不同的反射层分别位移以及多个反射层一起位移时，均能够稳定计算得到准确的补偿系数，进而得到较为精确的反应堆功率值。
以图1-4所示的反应堆为例，如图2所示，在堆芯10的外围，沿堆芯10的周向方向上均匀布置有3个可移动的反射层(反射层20-1、20-2、20-3)，可以沿堆芯10的径向方向布置4个核探测器(核探测器1～4)，沿堆芯10的轴向方向布置两个核探测器(核探测器5、6)。
其中，不同的反射层位移，对于核探测器1至核探测器4中任意一个核探测器位置处的中子注量率的影响关系均不同的，而对于核探测器5或者核探测器6位置处的中子注量率的影响关系是接近甚至一致的。例如，在核探测器6位置处，中子注量率随反射层20-1、20-2的位置变化关系的函数分别为G(x)1和G(x)2，而G(x)1和G(x)2非常接甚至是一致的。因而，在选定工作组探测器或补偿组探测器的位置时，可以选择核探测器5或者核探测器6的位置。
进一步地，如图4所示，当反应堆功率设定为额定功率不变时，反射层20从完全脱离堆芯10的位置(如图4中(a)所示)位移至反射层20完全遮挡堆芯10的位置(如图4中(d)所示)的过程中，核探测器1至核探测器6位置处的中子注量率变化曲线(即G(x)函数曲线)如图7a和7b所示，其中，横坐标表示反射层位置，即，堆芯被反射层遮挡的部分的占比，例如，反射层20完全脱离堆芯10时对应的横坐标为0％，堆芯10的一半被反射层20遮挡时对应的横坐标为50％，反射层20完全遮挡堆芯10时对应的横坐标为100％。
参见图7a和图7b，核探测器6位置处的G(x)函数值变化范围小于核探测器5位置处的G(x)函数值变化范围，并且核探测器5位置处的G(x)函数为单值函数。因此，可以选定核探测器6作为工作组探测器，核探测器5作为补偿组探测器。
此外，在一些实施例中，当选择工作组探测器和补偿组探测器时，还可以根据反应堆中的空间位置来选定补偿组探测器和/或工作组探测器。具体地，在设置工作组探测器和/或补偿组探测器时，其所在的空间位置与其他设备的位置之间不存在干涉，不会影响其他设备的布置。
在本实施例中，可以利用特定算法来计算得到补偿系数。当选定工作组探测器和补偿组探测器的位置后，可以获取反应堆功率恒定时工作组和补偿组探测器位置处的中子注量率随反射层位置变化关系的函数G(x)，并根据工作组探测器和补偿组探测器位置处的函数G(x)，经过相应的物理计算，来确定反射层位置关系函数f1以及补偿系数关系函数f2。
其中，函数f1为工作组探测器与补偿组探测器位置处的中子注量率的比值与反射层位置之间的数值关系式，可以描述补偿组和工作组探测器位置处的中子注量率的差异性变化与反射层位置之间的关系，其自变量为工作组探测器与补偿组探测器的功率测量值的比值，函数值为反射层位置。
具体地，函数f1可以为工作组探测器位置处G(x)函数与补偿组探测器位置处的G(x)函数的比值的反函数。当反应堆功率恒定在某个功率值不变时，可以通过计算反射层在不同位置时补偿组与工作组探测器的功率测量值的比值，或者，补偿组和工作组探测器位置处中子注量率的比值，再经过反函数计算，来得到函数f1。
以如图1-4所示的反应堆为例，当选定核探测器6作为工作组探测器并选定核探测器5作为补偿组探测器时，工作组探测器和补偿组探测器位置处中子注量率的随反射层位置变化关系曲线如图7b所示。将工作组探测器位置处的G(x)函数记为G(x)a，将补偿组探测器位置处的G(x)函数记为G(x)b，则函数G(x)a/G(x)b的反函数，即为函数f1，函数f1的曲线图如图8所示。
函数f2为反射层位置与工作组探测器的功率补偿系数之间的数值关系式，可以描述反射层位置与工作组探测器位置处的中子注量率的关系，其自变量为所述反射层位置，函数值为工作组探测器的功率补偿系数。
具体地，函数f2为反射层位于基准点处时位置处的中子注量率与工作组探测器位置处的G(x)函数的比值。在选定补偿的基准点后，可以将反射层位于基准点位置处时工作组探测器位置处的中子注量率除以工作组探测器位置处的G(x)函数，即可得到函数f2。在基准点位置处，函数f2输出的补偿系数为1，即不需要对工作组探测器的功率测量值进行补偿。
以如图1-4所示的反应堆为例，当选定核探测器6作为工作组探测器并选定核探测器5作为补偿组探测器时，工作组探测器位置处中子注量率的随反射层位置变化关系曲线如图7b所示。例如，当选定反射层20完全遮挡堆芯10的位置(即，图7b中横坐标为100％的位置)为基准点时，将工作组探测器位置处的中子注量率G0除以工作组探测器位置处的函数G(x)a，即G0/G(x)a，即为函数f2，函数f2的曲线图如图9所示。
本实施例中，在确定函数f2之前，需要选定一个补偿的基准点。核探测器通过测量中子注量率来间接测量反应堆功率，则工作组探测器位置处不同的中子注量率均对应一个不同的反应堆功率值。在工作组探测器位置处，反应堆功率为100％时所对应的中子注量率为n，本实施例中可以根据工作组探测器位置处的函数G(x)，选定工作组探测器位置处中子注量率为n时所对应的反射层位置，作为基准点。反射层位于基准点处时，工作组探测器的功率测量值为100％，当反射层位于其他位置时，即可相对于基准点所对应的100％的功率值，对工作组探测器所测量得到的功率测量值进行补偿。
当确定好函数f1、f2后，可以根据f1和f2来计算补偿系数。
参见图6，本实施例的补偿方法具体包括步骤S21～步骤S23。其中，步骤S21、S23分别与上述实施例中的步骤S11、S13相同，此处不再赘述。
在步骤S22中，首先根据工作组探测器和补偿组探测器的功率测量值之间的比值以及函数f1，确定当前的反射层位置；再根据确定的反射层位置以及函数f2，确定工作组探测器的功率补偿系数。当计算得到补偿系数后，即可对工作组探测器测量得到的功率测量值进行补偿，即，将工作组探测器测量得到的功率测量值乘以确定好的功率补偿系数，得到最终的反应堆当前的功率值。
其中，当反应堆中只设置一个反射层时，或者，当反应堆中所有反射层(反射层组)均同步位移时，本发明实施例中经过函数f1计算得到的反射层位置可以直接视为该反射层或反射层组的位置。当反应堆设置有多个反射层时，并且功率运行工况下仅有一个反射层(称为调节反射层)位移而其他反射层固定在正对堆芯的位置不动时，经过函数f1计算得到的反射层位置可以视为该调节反射层的位置。当反应堆中设置有多个反射层，并且功率运行工况下多个反射层均会位移时，经过函数f1计算得到的反射层位置可以视为所有反射层的平均有效位置。
在一些实施例中，可以设置数据处理系统来进行功率补偿。具体地，当设置好工作组和补偿组探测器后，可以将数据处理系统与工作组探测器和补偿组探测器连接，从而获取工作组和补偿组探测器的输出信号。数据处理系统中存储有计算机程序，该计算机程序被执行时可以使数据处理系统执行图6所示的实施例中的补偿方法。当确定好函数f1和f2中，可以将函数f1和f2输入至数据处理系统中，从而能够对功率测量值进行补偿。
具体地，当测量反应堆功率时，数据处理系统首先获得工作组和补偿组探测器的输出信号并转化为功率测量值，然后计算得到工作组和补偿组探测器的功率测量值的比值，并输入至函数f1中，使得函数f1输出反射层位置，再将反射层位置输入至函数f2，使得函数f2输功率补偿系数，最后将工作组探测器的功率测量值乘以功率补偿系数，即得到反应堆的当前功率值，并显示在数据处理系统的显示器上。此外，在实时测量反应堆功率时，还可以得到反应堆功率随时间变化的曲线。
需要说明的是，各核探测器位置处的中子注量率随反射层位置变化关系的函数G(x)是经过物理计算得到的。当反应堆的堆芯设计和反射层设计完成后，就可以计算反射层控制性位移对空间内中子注量率的影响关系，即，计算探测器位置处的中子注量率随反射层位置变化关系的函数G(x)。其中，可以利用反应堆设计的相关技术中的物理计算方法来计算得到函数G(x)，本实施例中并不限定具体的计算方法。示例地，可以通过模拟计算反应堆中反射层在不同位置时空间内各个位置处的中子注量率，得到函数G(x)。
为补偿反射层位移对移动功率测量的影响，传统的方法中，一般是直接测量反射层的实时位置，再根据反射层的位移对反应堆功率测量的数值进行针对性的补偿。该方法最大的弊端在于反射层或者控制棒的位置测量在一般的核电站和核设施中都是非安全级的测量参数，而反应堆功率测量一般都要求为安全级测量参数。如果采用直接测得反射层位置的方案需要将位置测量从非安全级提升至安全级，并且位置的测量精度也需要相应提升，而现有的技术条件难以满足该要求。
相比于传统方法，采用本发明实施例中的补偿方法，不会增加测量参数的类型，不用依赖其他测量系统工作。本发明实施例的补偿方法所使用的测量参数均为核探测器的测量参数，相比于不进行功率补偿时的测量方法，仅增加参数的数量，但不会增加参数的类型。相较于需要利用除核探测器测量参数以外的其他参数的一些潜在方法，本发明实施例的方法具有不增加测量参数类型的优势，对于反应堆功率的测量仍然可由核测量系统独立完成。
并且，本发明实施例中的方法，不改变其反应堆功率的测量原理。无论是工作组探测器，还是补偿组探测器，其测量原理均是将核探测器输出的电参数转换为中子注量率后再转换为反应堆功率，不需要借助其他的测量原理。
相较于直接将反射层位置进行测量继而补偿反射层位移对反应堆功率测量影响的方法而言，本发明实施例中的补偿方法能够满足对反应堆功率测量设备的安全级要求。
此外，本发明实施例中的补偿方法还可以附带获得反射层位置的参数。本发明实施例中的反射层位置参数是通过不同位置的核探测器参数间接计算得到的，其可以应用于多个方面。例如，将本发明实施例中间接计算得到的反射层位置与直接测得的反射层位置参数相互校核，或者，根据该反射层位置的变化情况判断反射层的运动方向，或者，根据该反射层位置和其他参数一起估算反应堆的燃耗深度等。
本发明实施例中的补偿方法，利用多个核探测器的功率测量值的差异性变化经过特定算法计算得到反射层位置，再利用反射层位置经过特定算法来补偿反应堆功率测量值的偏差，可以极大地减小反应堆功率测量时由于反射层控制性位移造成的偏差，提高了核探测器测量反应堆功率的测量精度。
对于本发明的实施例，还需要说明的是，在不冲突的情况下，本发明的实施例及实施例中的特征可以相互组合以得到新的实施例。
以上，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，本发明的保护范围应以权利要求的保护范围为准。

Claims (9)

1.一种反应堆功率测量的补偿方法，其特征在于，所述反应堆的堆芯外围设置有可移动反射层，通过所述反射层的位移来控制所述反应堆的反应性，所述反射层外设置有多个核探测器；所述补偿方法用于补偿反射层位移对反应堆功率测量值造成的偏差，其包括：
选定工作组探测器和补偿组探测器的位置；
获取所述多个核探测器中工作组探测器和补偿组探测器的功率测量值；
根据所述工作组探测器和补偿组探测器的功率测量值之间的比值，确定反射层位置；
根据所述反射层位置确定所述工作组探测器的功率补偿系数；
利用所述功率补偿系数来补偿所述工作组探测器的功率测量值，得到所述反应堆的当前功率值；
所述选定工作组探测器和补偿组探测器的位置，包括：
根据反应堆功率恒定时各核探测器位置处的中子注量率随反射层位置变化关系的函数G(x)，选定所述工作组探测器和补偿组探测器的位置；
其中，
当反应堆功率恒定时，所述工作组探测器位置处的G(x)的函数值变化范围小于所述补偿组探测器位置处的G(x)的函数值变化范围，所述函数值变化范围为最大函数值与最小函数值的差值；
所述补偿组探测器位置处的G(x)为单值函数，且所述补偿组探测器位置处的G(x)与所述工作组探测器位置处的G(x)之比也为单值函数。
2.根据权利要求1所述的补偿方法，其特征在于，还包括：
获取反应堆功率恒定时各核探测器位置处的中子注量率随反射层位置变化关系的函数G(x)；
根据工作组探测器和补偿组探测器位置处的函数G(x)，确定反射层位置关系函数f1以及补偿系数关系函数f2；
其中，
所述函数f1为所述工作组探测器与补偿组探测器位置处的中子注量率的比值与反射层位置之间的数值关系式；
所述函数f2为反射层位置与所述工作组探测器的功率补偿系数之间的数值关系式。
3.根据权利要求2所述的补偿方法，其特征在于，所述函数f1为工作组探测器的G(x)函数与补偿组探测器的G(x)函数的比值的反函数。
4.根据权利要求2所述的补偿方法，其特征在于，所述函数f2为反射层位于基准点处时工作组探测器位置处的中子注量率与工作组探测器的G(x)函数的比值。
5.根据权利要求2-4任一项所述的补偿方法，其特征在于，所述根据工作组探测器和补偿组探测器的功率测量值之间的比值，确定工作组探测器的功率补偿系数，包括：
根据工作组探测器和补偿组探测器的功率测量值之间的比值以及所述函数f1，确定当前的反射层位置；
根据所述反射层位置以及所述函数f2，确定所述工作组探测器的功率补偿系数。
6.根据权利要求1-4任一项所述的补偿方法，其特征在于，所述工作组探测器和补偿组探测器的功率测量值之间的比值为：
工作组探测器功率测量值除以补偿组探测器功率测量值；或者
补偿组探测器功率测量值除以工作组探测器功率测量值；或者
工作组探测器功率测量值与补偿组探测器功率测量值进行对数运算后的差值。
7.根据权利要求1所述的补偿方法，其特征在于，选定工作组探测器的位置，还包括：
当反应堆的堆芯外围设置有多个反射层时，根据各反射层对各核探测器位置处中子注量率的影响关系，选定工作组探测器的位置；
其中，各反射层对工作组探测器位置处的中子注量率的影响关系一致。
8.根据权利要求1所述的补偿方法，其特征在于，选定补偿组探测器的位置，还包括：
当反应堆的堆芯外围设置有多个反射层时，根据各反射层对各核探测器位置处中子注量率的影响关系，选定补偿组探测器的位置；
其中，各反射层对补偿组探测器位置处的中子注量率的影响关系一致。
9.根据权利要求1、7-8任一项所述的补偿方法，其特征在于，还包括：根据反应堆中的空间位置，选定补偿组探测器和/或工作组探测器的位置。
CN202211004907.6A 2022-08-22 2022-08-22 反应堆功率测量的补偿方法 Active CN115295179B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211004907.6A CN115295179B (zh) 2022-08-22 2022-08-22 反应堆功率测量的补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211004907.6A CN115295179B (zh) 2022-08-22 2022-08-22 反应堆功率测量的补偿方法

Publications (2)

Publication Number Publication Date
CN115295179A CN115295179A (zh) 2022-11-04
CN115295179B true CN115295179B (zh) 2023-12-12

Family

ID=83831070

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211004907.6A Active CN115295179B (zh) 2022-08-22 2022-08-22 反应堆功率测量的补偿方法

Country Status (1)

Country Link
CN (1) CN115295179B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62229098A (ja) * 1986-03-31 1987-10-07 株式会社東芝 制御棒引抜監視装置
US5091139A (en) * 1989-06-26 1992-02-25 General Electric Company Automated thermal limit monitor
JPH06194482A (ja) * 1992-12-24 1994-07-15 Hitachi Ltd 原子炉
KR20010038820A (ko) * 1999-10-27 2001-05-15 이종훈 개선된 노외계측기 교정방법
JP2005233751A (ja) * 2004-02-19 2005-09-02 Toshiba Corp 反射体制御方式の高速炉およびその中性子反射体
JP2008191069A (ja) * 2007-02-07 2008-08-21 Toshiba Corp 反射体制御方式高速炉の中性子検出器設置方法及び運転方法、並びに反射体制御方式高速炉
CN106297921A (zh) * 2016-08-31 2017-01-04 苏州热工研究院有限公司 一种利用自给能探测器瞬发γ信号进行信号补偿的方法
CN107230505A (zh) * 2017-06-21 2017-10-03 中国核动力研究设计院 一种反应堆核功率监测方法及系统
CN113421669A (zh) * 2021-06-17 2021-09-21 中国核动力研究设计院 基于局部非线性修正的堆芯功率分布在线重构方法及系统
WO2022113592A1 (ja) * 2020-11-30 2022-06-02 日立Ｇｅニュークリア・エナジー株式会社 原子炉制御装置および原子炉制御方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6453262B2 (ja) * 2016-03-09 2019-01-16 日立Ｇｅニュークリア・エナジー株式会社 制御棒操作監視方法及び制御棒操作監視システム

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62229098A (ja) * 1986-03-31 1987-10-07 株式会社東芝 制御棒引抜監視装置
US5091139A (en) * 1989-06-26 1992-02-25 General Electric Company Automated thermal limit monitor
JPH06194482A (ja) * 1992-12-24 1994-07-15 Hitachi Ltd 原子炉
KR20010038820A (ko) * 1999-10-27 2001-05-15 이종훈 개선된 노외계측기 교정방법
JP2005233751A (ja) * 2004-02-19 2005-09-02 Toshiba Corp 反射体制御方式の高速炉およびその中性子反射体
JP2008191069A (ja) * 2007-02-07 2008-08-21 Toshiba Corp 反射体制御方式高速炉の中性子検出器設置方法及び運転方法、並びに反射体制御方式高速炉
CN106297921A (zh) * 2016-08-31 2017-01-04 苏州热工研究院有限公司 一种利用自给能探测器瞬发γ信号进行信号补偿的方法
CN107230505A (zh) * 2017-06-21 2017-10-03 中国核动力研究设计院 一种反应堆核功率监测方法及系统
WO2022113592A1 (ja) * 2020-11-30 2022-06-02 日立Ｇｅニュークリア・エナジー株式会社 原子炉制御装置および原子炉制御方法
CN113421669A (zh) * 2021-06-17 2021-09-21 中国核动力研究设计院 基于局部非线性修正的堆芯功率分布在线重构方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Design Study of Sodium Cooled Small Fast Reactor;UEDA. et al.;GENES4/ANP2003;Paper114(1-10) *
Temporal Variation of the Neutron Flux in the Carousel Facility of a TRIGA Reactor;Jacimovie. et al.;International Conference Nuclear Energy for New Europe 2003;116.(1-8) *

Also Published As

Publication number Publication date
CN115295179A (zh) 2022-11-04

Similar Documents

Publication Publication Date Title
KR101651893B1 (ko) 신경망 회로를 이용한 원자로 노심의 축방향 출력분포 합성방법 및 그 방법이 적용된 노심감시계통
CN103411684B (zh) 测量反应腔室内薄膜温度的方法
CN107833568A (zh) 一种屏幕亮度调节方法、装置及显示终端
EP0560396A2 (en) Nuclear magnetic resonance imaging with improved image quality and operation efficiency
CN115295179B (zh) 反应堆功率测量的补偿方法
CN107966928A (zh) 背隙误差补偿方法、装置、系统、计算机设备及放射治疗设备
JP6469388B2 (ja) 電界強度算出プログラム、電界強度算出装置及び電界強度算出方法
CN102927925A (zh) 基于自适应加权融合算法的倾斜光纤光栅应变测量方法
KR101628404B1 (ko) 혼합된 인-코어 맵 생성 방법 및 고정형 계측 기구의 교정에 대한 적용 방법
JPH03214099A (ja) 原子炉炉心における出力分布の決定方法及び原子炉炉心周囲にある中性子検出器の較正方法
KR100856379B1 (ko) 유전알고리즘을 이용한 노외계측기의 설정치 결정방법
JPH0282196A (ja) 棒位置決定方法及び装置
CN109799391B (zh) 导体交流电阻测量方法、系统及计算机存储介质
Bettini et al. Iterative Axisymmetric Identification Algorithm (IAIA) for real-time reconstruction of the plasma boundary of ITER
KR100694436B1 (ko) 원자로 노심내 바나듐 및 백금계측기의 지연신호 보상장치
CN113936823B (zh) 获取堆外探测器刻度系数的方法、系统、装置及存储介质
CN106767952B (zh) 一种电感式位移传感器的干扰消除方法
CN112179550B (zh) 一种四支点压电测力仪的优化布置方法
CN110737016B (zh) 一种用于回旋加速器的主径向探针信号标定系统与方法
CN115656619B (zh) 一种发电机功率因数在dcs系统中的显示方法及系统
CN113345613A (zh) 高温气冷堆主蒸汽系统的监测控制方法及装置
CN103353277B (zh) 一种用于飞机装配工装数字化测量的激光跟踪仪设站方法
CN113570147A (zh) 一种薄壁件多工况加工误差快速预测方法及设备
JPH1020072A (ja) 炉心性能計算装置
CN110767322A (zh) 一种基于响应面模型的海洋浮式平台热点应力推算方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant