CN111799005A - 一种不同量程归一化方法、存储介质和实时在线反应性仪 - Google Patents

一种不同量程归一化方法、存储介质和实时在线反应性仪 Download PDF

Info

Publication number
CN111799005A
CN111799005A CN202010810975.6A CN202010810975A CN111799005A CN 111799005 A CN111799005 A CN 111799005A CN 202010810975 A CN202010810975 A CN 202010810975A CN 111799005 A CN111799005 A CN 111799005A
Authority
CN
China
Prior art keywords
range
data
reactivity
irc
control system
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.)
Pending
Application number
CN202010810975.6A
Other languages
English (en)
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.)
Cai Yueqing
Original Assignee
Fuzhou Qizhenggu Technology Co ltd
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 Fuzhou Qizhenggu Technology Co ltd filed Critical Fuzhou Qizhenggu Technology Co ltd
Priority to CN202010810975.6A priority Critical patent/CN111799005A/zh
Publication of CN111799005A publication Critical patent/CN111799005A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21C—NUCLEAR REACTORS
    • G21C17/00—Monitoring; Testing ; Maintaining
    • G21C17/10—Structural combination of fuel element, control rod, reactor core, or moderator structure with sensitive instruments, e.g. for measuring radioactivity, strain
    • G21C17/104—Measuring reactivity
    • 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)
  • Plasma & Fusion (AREA)
  • General Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明涉及反应堆测试领域，特别涉及一种不同量程归一化方法、存储介质和实时在线反应性仪。所述一种不同量程归一化方法，包括步骤：获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值。通过获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值，保证了全量程反应性连续计算。

Description

一种不同量程归一化方法、存储介质和实时在线反应性仪
技术领域
本发明涉及反应堆测试领域，特别涉及一种不同量程归一化方法、存储介质和实时在线反应性仪。
背景技术
反应性是反应堆有效增值系数Keff对临界值的相对偏离量，是反应堆物理启动时反映反应堆运行情况的重要参数。它的准确测量对于反应堆的安全运行和发挥其经济效益有其重要意义。反应性仪是直接测量反应堆反应性的一种设备，也是核电厂反应堆启动物理试验的关键设备。
核反应堆的功率水平变化十几个数量级，一般用3种探测器来覆盖，分别是功率量程(PRC)、中间量程(IRC)、源量程(SRC)。这几个量程给出的功率水平的单位都不一样，在量程切换时，反应性无法连续计算。故而如何在接收不同量程信号后，将数据转换为相同单位的数值，保证全量程反应性连续计算成了亟需解决的问题。
发明内容
为此，需要提供一种不同量程归一化方法，用以解决在接收到不同量程后，量程功率水平的单位不一样，在量程切换时，反应性无法连续计算的问题。具体技术方案如下：
一种不同量程归一化方法，包括步骤：获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值。
进一步的，所述不同量程通道数据包括：功率量程、中间量程和源量程；
所述“将所述不同量程通道数据转换至相同单位的数值”，还包括步骤：
中间量程功率量程转换系数为：IRC2PRC＝PRC/IRC；
源量程中间量程转换系数为：SRC2IRC＝IRC/SRC；
初始化所述IRC2PRC和所述SRC2IRC；
所述PRC为功率量程测量数据，所述IRC为中间量程测量数据，所述SRC为源量程测量数据。
为解决上述技术问题，还提供了一种存储介质，具体技术方案如下：
一种存储介质，其中存储有指令集，所述指令集用于执行：获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值。
进一步的，所述指令集还用于执行：所述不同量程通道数据包括：功率量程、中间量程和源量程；
所述“将所述不同量程通道数据转换至相同单位的数值”，还包括步骤：
中间量程功率量程转换系数为：IRC2PRC＝PRC/IRC；
源量程中间量程转换系数为：SRC2IRC＝IRC/SRC；
初始化所述IRC2PRC和所述SRC2IRC；
所述PRC为功率量程测量数据，所述IRC为中间量程测量数据，所述SRC为源量程测量数据。
为解决上述技术问题，还提供了一种实时在线反应性仪，具体技术方案如下：
一种实时在线反应性仪，所述反应性仪用于：获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值。
进一步的，所述反应性仪还用于：所述不同量程通道数据包括：功率量程、中间量程和源量程；
所述“将所述不同量程通道数据转换至相同单位的数值”，还包括步骤：
中间量程功率量程转换系数为：IRC2PRC＝PRC/IRC；
源量程中间量程转换系数为：SRC2IRC＝IRC/SRC；
初始化所述IRC2PRC和所述SRC2IRC；
所述PRC为功率量程测量数据，所述IRC为中间量程测量数据，所述SRC为源量程测量数据。
进一步的，所述反应性仪通过通信连接计算机控制系统，通信连接方式包括以下中的一种或多种：使用modbus RTU协议、使用modbus TCP协议、使用modbus ASCII协议；
所述反应性仪用于：向所述计算机控制系统发送数据读取请求，并接收所述计算机控制系统返回的数据，根据所述返回的数据计算反应性，并发送计算结果至计算机控制系统进行输出显示。
进一步的，所述反应性仪通过通信连接计算机控制系统，通信连接方式包括以下中的一种或多种：使用modbus RTU协议、使用modbus TCP协议、使用modbus ASCII协议；
所述反应性仪用于：接收所述计算机控制系统写入的数据，根据所述写入的数据计算反应性，并接收所述计算机控制系统的读取数据指令，返回计算结果至计算机控制系统进行输出显示。
进一步的，所述反应性仪还用于：计算中子源强度；
所述“计算中子源强度”，包括步骤：
测量得深次临界稳定状态下的功率水平n；
计算在此深次临界状态下的keff或反应性ρ，及瞬发中子平均代时间l或瞬发中子寿命l0；
根据下式计算出中子源强度S：
Figure BDA0002630950730000031
或
Figure BDA0002630950730000032
所述反应性ρ定义如下：
Figure BDA0002630950730000033
所述keff本征值用来描述核反应堆系统的次临界度，所述次临界方程如下：
Figure BDA0002630950730000041
其中l0＝keff*l；
在Keff<0.98或ρ<-2000pcm时，中子源强度的测量误差符合工程应用需求。
进一步的，所述反应性仪通过串口与所述计算机控制系统进行通信，或所述反应性仪通过TCP/IP网络与所述计算机控制系统进行通信。
本发明的有益效果是：通过获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值，保证了全量程反应性连续计算。
附图说明
图1为具体实施方式所述一种不同量程归一化方法的流程图；
图2为具体实施方式所述归一化后的功率的计算流程图；
图3为具体实施方式所述一种存储介质的模块示意图；
图4为具体实施方式所述一种实时在线反应性仪的模块示意图；
图5为具体实施方式所述相对误差随反应性的变化示意图；
图6为具体实施方式所述功率、反应性随时间变化示意图。
附图标记说明：
300、存储介质，
400、实时在线反应性仪。
具体实施方式
为详细说明技术方案的技术内容、构造特征、所实现目的及效果，以下结合具体实施例并配合附图详予说明。
请参阅图1至图2，在本实施方式中，一种不同量程归一化方法可应用在存储介质上，所述存储介质包括但不限于：反应性仪、个人计算机、服务器、通用计算机、专用计算机、网络设备、嵌入式设备、可编程设备等。具体实施方式如下：
步骤S101：获取不同量程通道数据。
步骤S102：将所述不同量程通道数据转换至相同单位的数值。
通过获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值，保证了全量程反应性连续计算。
其中因核反应堆的功率水平变化十几个数量级，一般用3种探测器来覆盖，即所述不同量程通道数据包括：功率量程、中间量程和源量程；
所述“将所述不同量程通道数据转换至相同单位的数值”，还包括步骤：
中间量程功率量程转换系数为：IRC2PRC＝PRC/IRC；
源量程中间量程转换系数为：SRC2IRC＝IRC/SRC；
初始化所述IRC2PRC和所述SRC2IRC；
所述PRC为功率量程测量数据，所述IRC为中间量程测量数据，所述SRC为源量程测量数据。具体可如下：
a)定义中间量程功率量程转换系数IRC2PRC＝PRC/IRC，源量程中间量程转换系数SRC2IRC＝IRC/SRC，将这两个参数初始化。
式中，IRC2PRC是指中间量程功率量程转换系数，SRC2IRC是指源量程中间量程转换系数，PRC指功率量程测量数据，IRC指中间量程测量数据，SRC指源量程测量数据。
b)定义中间量程功率量程转换系数重叠切换点P10，源量程中间量程重叠切换点P6，将这两个参数初始化。
量程重叠是指在某个区域，不同的测量通道都可以覆盖此区域，那在此区域定义一个点，大于此值用探测器1的数据，小于此值用探测器2的数据，这个点就是量程重叠切换点。比如，某反应堆源量程测量范围是1～1E6cps，覆盖4E-10％FP～4E-4％FP功率范围；中间量程测量范围是1E-11A～1E-3A，覆盖4E-7％FP～40％FP功率范围；则SRC与IRC有3个量级的重叠，重叠切换点P6设定为1E-5％FP，则P6在重叠区域。但功率大于P6时用IRC数据，而小于P6时用SRC数据；在P6～P6/2区间进行量程转换系数计算，保证不同量程归一化的正确性。
其中参数powerPRC用来表示归一化后的功率，其计算流程如图2所示：
步骤S201：IRC2PRC、SRC2IRC、P10、P6参数初始化，输入数据PRC、IRC、SRC。
步骤S202：PRC>P10？
是，执行步骤S203：powerPRC＝PRC。并跳至执行步骤S211：返回powerPRC。
否，执行步骤S204：PRC>P10/2？
是，执行步骤S205：IRC2PRC＝PRC/IRC，powerPRC＝IRC*IRC2PRC。并跳至执行步骤S211：返回powerPRC。
否，执行步骤S206：IRC>P6？
是，执行步骤S207：powerPRC＝IRC*IRC2PRC。并跳至执行步骤S211：返回powerPRC。
否，执行步骤S208：IRC>P6/2？
是，执行步骤S209：SRC2IRC＝IRC/SRC,powerPRC＝SRC*SRC2IRC*IRC2PRC。并跳至执行步骤S211：返回powerPRC。
否，执行步骤S210：powerPRC＝SRC*SRC2IRC*IRC2PRC。执行步骤S211：返回powerPRC。
请参阅图3，在本实施方式中，所述存储介质300包括但不限于：反应性仪、个人计算机、服务器、通用计算机、专用计算机、网络设备、嵌入式设备、可编程设备等。具体实施方式如下
一种存储介质300，其中存储有指令集，所述指令集用于执行：获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值。
通过获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值，保证了全量程反应性连续计算。
进一步的，所述指令集还用于执行：所述不同量程通道数据包括：功率量程、中间量程和源量程；
所述“将所述不同量程通道数据转换至相同单位的数值”，还包括步骤：
中间量程功率量程转换系数为：IRC2PRC＝PRC/IRC；
源量程中间量程转换系数为：SRC2IRC＝IRC/SRC；
初始化所述IRC2PRC和所述SRC2IRC；
所述PRC为功率量程测量数据，所述IRC为中间量程测量数据，所述SRC为源量程测量数据。
具体可如下：
a)定义中间量程功率量程转换系数IRC2PRC＝PRC/IRC，源量程中间量程转换系数SRC2IRC＝IRC/SRC，将这两个参数初始化。
式中，IRC2PRC是指中间量程功率量程转换系数，SRC2IRC是指源量程中间量程转换系数，PRC指功率量程测量数据，IRC指中间量程测量数据，SRC指源量程测量数据。
b)定义中间量程功率量程转换系数重叠切换点P10，源量程中间量程重叠切换点P6，将这两个参数初始化。
量程重叠是指在某个区域，不同的测量通道都可以覆盖此区域，那在此区域定义一个点，大于此值用探测器1的数据，小于此值用探测器2的数据，这个点就是量程重叠切换点。比如，某反应堆源量程测量范围是1～1E6cps，覆盖4E-10％FP～4E-4％FP功率范围；中间量程测量范围是1E-11A～1E-3A，覆盖4E-7％FP～40％FP功率范围；则SRC与IRC有3个量级的重叠，重叠切换点P6设定为1E-5％FP，则P6在重叠区域。但功率大于P6时用IRC数据，而小于P6时用SRC数据；在P6～P6/2区间进行量程转换系数计算，保证不同量程归一化的正确性。
其中参数powerPRC用来表示归一化后的功率，其计算流程如图2所示：
步骤S201：IRC2PRC、SRC2IRC、P10、P6参数初始化，输入数据PRC、IRC、SRC。
步骤S202：PRC>P10？
是，执行步骤S203：powerPRC＝PRC。并跳至执行步骤S211：返回powerPRC。
否，执行步骤S204：PRC>P10/2？
是，执行步骤S205：IRC2PRC＝PRC/IRC，powerPRC＝IRC*IRC2PRC。并跳至执行步骤S211：返回powerPRC。
否，执行步骤S206：IRC>P6？
是，执行步骤S207：powerPRC＝IRC*IRC2PRC。并跳至执行步骤S211：返回powerPRC。
否，执行步骤S208：IRC>P6/2？
是，执行步骤S209：SRC2IRC＝IRC/SRC,powerPRC＝SRC*SRC2IRC*IRC2PRC。并跳至执行步骤S211：返回powerPRC。
否，执行步骤S210：powerPRC＝SRC*SRC2IRC*IRC2PRC。执行步骤S211：返回powerPRC。
请参阅图4，在本实施方式中，一种实时在线反应性仪400的具体实施方式如下：
一种实时在线反应性仪400，所述反应性仪用于：获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值。
通过获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值，保证了全量程反应性连续计算。
进一步的，所述反应性仪还用于：所述不同量程通道数据包括：功率量程、中间量程和源量程；
所述“将所述不同量程通道数据转换至相同单位的数值”，还包括步骤：
中间量程功率量程转换系数为：IRC2PRC＝PRC/IRC；
源量程中间量程转换系数为：SRC2IRC＝IRC/SRC；
初始化所述IRC2PRC和所述SRC2IRC；
所述PRC为功率量程测量数据，所述IRC为中间量程测量数据，所述SRC为源量程测量数据。具体可如下：
a)定义中间量程功率量程转换系数IRC2PRC＝PRC/IRC，源量程中间量程转换系数SRC2IRC＝IRC/SRC，将这两个参数初始化。
式中，IRC2PRC是指中间量程功率量程转换系数，SRC2IRC是指源量程中间量程转换系数，PRC指功率量程测量数据，IRC指中间量程测量数据，SRC指源量程测量数据。
b)定义中间量程功率量程转换系数重叠切换点P10，源量程中间量程重叠切换点P6，将这两个参数初始化。
量程重叠是指在某个区域，不同的测量通道都可以覆盖此区域，那在此区域定义一个点，大于此值用探测器1的数据，小于此值用探测器2的数据，这个点就是量程重叠切换点。比如，某反应堆源量程测量范围是1～1E6cps，覆盖4E-10％FP～4E-4％FP功率范围；中间量程测量范围是1E-11A～1E-3A，覆盖4E-7％FP～40％FP功率范围；则SRC与IRC有3个量级的重叠，重叠切换点P6设定为1E-5％FP，则P6在重叠区域。但功率大于P6时用IRC数据，而小于P6时用SRC数据；在P6～P6/2区间进行量程转换系数计算，保证不同量程归一化的正确性。
其中参数powerPRC用来表示归一化后的功率，其计算流程如图2所示：
步骤S201：IRC2PRC、SRC2IRC、P10、P6参数初始化，输入数据PRC、IRC、SRC。
步骤S202：PRC>P10？
是，执行步骤S203：powerPRC＝PRC。并跳至执行步骤S211：返回powerPRC。
否，执行步骤S204：PRC>P10/2？
是，执行步骤S205：IRC2PRC＝PRC/IRC，powerPRC＝IRC*IRC2PRC。并跳至执行步骤S211：返回powerPRC。
否，执行步骤S206：IRC>P6？
是，执行步骤S207：powerPRC＝IRC*IRC2PRC。并跳至执行步骤S211：返回powerPRC。
否，执行步骤S208：IRC>P6/2？
是，执行步骤S209：SRC2IRC＝IRC/SRC,powerPRC＝SRC*SRC2IRC*IRC2PRC。并跳至执行步骤S211：返回powerPRC。
否，执行步骤S210：powerPRC＝SRC*SRC2IRC*IRC2PRC。执行步骤S211：返回powerPRC。
进一步的，所述反应性仪通过通信连接计算机控制系统，通信连接方式包括以下中的一种或多种：使用modbus RTU协议、使用modbus TCP协议、使用modbus ASCII协议；
所述反应性仪用于：向所述计算机控制系统发送数据读取请求，并接收所述计算机控制系统返回的数据，根据所述返回的数据计算反应性，并发送计算结果至计算机控制系统进行输出显示。需要说明的是modbus协议是一个master/slave或客户端/服务器架构的协议，在该方式中，反应性仪作为master或客户端节点，计算机控制系统作为slave或服务器节点，反应性仪用于向计算机控制系统发送数据请求，根据返回的数据计算反应性，并发送计算结果给计算机控制系统，计算机控制系统对计算结果进行显示。
或
进一步的，所述反应性仪通过通信连接计算机控制系统，通信连接方式包括以下中的一种或多种：使用modbus RTU协议、使用modbus TCP协议、使用modbus ASCII协议；
所述反应性仪用于：接收所述计算机控制系统写入的数据，根据所述写入的数据计算反应性，并接收所述计算机控制系统的读取数据指令，返回计算结果至计算机控制系统进行输出显示。在该方式中，反应性仪做为slave或或服务器节点，计算机控制系统作为master或客户端节点，计算机控制系统可直接向反应性仪写入数据，反应性仪无需自己发起数据请求来获取数据，反应性仪被写入数据后，会根据所述写入数据计算反应性，计算完后，计算机控制系统会主动读取计算结果，并对计算结果进行显示。
上述两种方式在modbus协议下均是可行的，可根据实际情况根据具体需要进行个性化设置。
在实际应用场景中，反应性仪通过通信连接的计算机控制系统是核电厂分布式计算机控制系统(DCS)，在本实施方式中，所述反应性仪设置有串口，所述反应性仪可通过串口与所述计算机控制系统进行通信，其中modbus RTU、modbus ASCII是基于串口通讯的，modbus TCP是基于TCP/IP网络通讯的。所述反应性仪通过上述任意一种通信协议向所述计算机控制系统发送数据读取请求，并接收所述计算机控制系统返回的数据；根据所述返回的数据计算反应性，将反应性等计算结果写入指定的寄存器，DCS可以实时显示反应性等重要参数。其中通信协议与串口的选择如图6所示，设置有UI界面。
所述反应性仪也通过上述任意一种通信协议接收所述计算机控制系统写入的数据，根据所述写入的数据计算反应性，并接收所述计算机控制系统的读取数据指令，DCS可以实时显示反应性等重要参数。
所述反应性仪上设置有串口，通过该串口所述反应性仪可通过modbus RTU协议和/或modbus ASCII协议与计算机控制系统进行通信，用于向所述计算机控制系统发送数据读取请求，并接收所述计算机控制系统返回的数据；根据所述返回的数据计算反应性，并发送计算结果至计算机控制系统进行输出显示；或反应性仪用于：接收所述计算机控制系统写入的数据，根据所述写入的数据计算反应性，并接收所述计算机控制系统的读取数据指令，返回数据给计算机控制系统。因所述反应性仪可使用多种modbus协议，并可根据需要设置为master/slave或客户端/服务器，故可根据需要适配不同的计算机控制系统。且所述反应性仪支持modbus RTU协议和modbus ASCII协议，可以避免modbus TCP模式接入带来的安全风险。
进一步的，所述反应性仪用于：计算中子源强度；
所述“计算中子源强度”，包括步骤：
测量得深次临界稳定状态下的功率水平n；
计算在此深次临界状态下的keff或反应性ρ，及瞬发中子平均代时间l或瞬发中子寿命l0；
根据下式计算出中子源强度S：
Figure BDA0002630950730000121
或
Figure BDA0002630950730000122
所述反应性ρ定义如下：
Figure BDA0002630950730000123
所述keff本征值用来描述核反应堆系统的次临界度，所述次临界方程如下：
Figure BDA0002630950730000124
其中l0＝keff*l；在Keff<0.98或ρ<-2000pcm时，中子源强度的测量误差符合工程应用需求。在本实施方式中，所谓符合工程应用需求即：在Keff<0.98或ρ<-2000pcm时，用该方法确定的中子源强度的最大相对误差小于25％。而后续由此反应性的测量(计算)最大相对误差也小于25％，符合工程应用需求。
通过上述反应性仪计算中子源强度，在Keff<0.98或ρ<-2000pcm时，中子源强度的测量误差符合工程应用需求，从而保证反应性在次临界状态下计算的正确性和全量程覆盖。
误差分析表明，在Keff<0.98(或ρ<-2000pcm)时，通过上述中子源强度的计算方法，可保证反应性在次临界状态下的测量误差满足工程应用要求，从而保证反应性在次临界状态下计算的正确性和全量程覆盖。其中反应堆次临界度是反应堆动力学行为的重要物理特征量，通常，用keff本征值来描述核反应堆系统的次临界度。在描述次临界度时，一般此时keff<1，ρ<0。
以下进行具体分析：一般中子源强度S是固定的，根据式(3)可知，反应性ρ与中子密度n(功率水平)成反比。可以分几种情况讨论：
a)接近临界，即反应性ρ接近0，中子密度n很大，此时n的测量相对误差比较小；但是反应性ρ计算有一个比较固定的计算误差即dρ，则此时其相对误差极大。即此时用此方法计算中子源强度S的误差极大，是不合适的。
b)深次临界，即反应性ρ绝对值很大(ρ小于0)，中子密度n比较小，此时n的测量相对误差比较大；反应性ρ计算有一个比较固定的计算误差即dρ，由于ρ绝对值大，此时其相对误差大幅减小。即此时用此方法计算中子源强度S的误差较小，是合适的。
c)反应性在以上情况a)和情况b)中间。相比较与情况b)，n的测量相对误差减少一半；反应性ρ计算的相对误差增大一倍。n的测量相对误差本身就不大，减少一半贡献有限，因此此时计算中子源强度S的误差比b)大。
例如，根据某反应堆的实测数据，源量程测量的相对误差大约2％，即
Figure BDA0002630950730000131
深次临界状态下，源量程测量值低，保守假设其测量的相对误差为4％，即
Figure BDA0002630950730000132
另外，核设计软件计算反应堆某状态下的反应性ρ或keff的行业允许误差为500pcm，即dρ＝500pcm。如果在接近临界状态，如ρ＝-100pcm，此时
Figure BDA0002630950730000141
如果深次临界状态，如ρ＝-10000pcm，此时
Figure BDA0002630950730000142
如果ρ＝-5000pcm，此时
Figure BDA0002630950730000143
如图5所示，反应性在小于-2000pcm时，反应性计算的相对误差变化比较缓慢；而在大于-2000pcm时，相对误差急剧增加。
以上分析可以知道，用次临界方程获得中子源强度的误差主要取决于反应性计算的相对误差，而这又取决于反应性(次临界度)。当反应性小于-2000pcm时，中子源强度的相对误差小于25％，满足工程应用的要求。一般可以在反应堆刚开始启堆阶段按本方法计算，即反应性在-10000～-5000pcm时，中子源强度的相对误差小于10％，满足工程应用的要求。
例如，某反应堆中子源强度S＝2.5e-4，有效增值系数keff＝0.9523(反应性ρ＝-5000pcm)，此时反应堆功率水平为1E-7％FP，瞬发中子平均代时间l＝2E-5s。假设计算误差-600pcm，此时计算得到的有效增值系数keff＝0.9470(反应性ρ＝-5600pcm)，测量得到的反应堆功率水平为1E-7％FP，则中子源强度
Figure BDA0002630950730000144
此中子源强度在反应性仪中参与反应性计算。如果此时反应堆达临界，反应性仪测量的反应性见图6中的“反应性1”，反应性的真值见图6中的“反应性2”，二者的差异见图6中的“反应性差”。此时反应性测量误差最大600pcm，随着反应堆的逐渐逼近临界，误差也逐渐减小。
故本反应性仪解决了中子源强度的计算的误差问题，也解决了次临界状态下的反应性测量问题，测量误差满足工程应用要求。
需要说明的是，尽管在本文中已经对上述各实施例进行了描述，但并非因此限制本发明的专利保护范围。因此，基于本发明的创新理念，对本文所述实施例进行的变更和修改，或利用本发明说明书及附图内容所作的等效结构或等效流程变换，直接或间接地将以上技术方案运用在其他相关的技术领域，均包括在本发明的专利保护范围之内。

Claims (10)

1.一种不同量程归一化方法，其特征在于，包括步骤：获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值。
2.根据权利要求1所述的一种不同量程归一化方法，其特征在于，
所述不同量程通道数据包括：功率量程、中间量程和源量程；
所述“将所述不同量程通道数据转换至相同单位的数值”，还包括步骤：
中间量程功率量程转换系数为：IRC2PRC＝PRC/IRC；
源量程中间量程转换系数为：SRC2IRC＝IRC/SRC；
初始化所述IRC2PRC和所述SRC2IRC；
所述PRC为功率量程测量数据，所述IRC为中间量程测量数据，所述SRC为源量程测量数据。
3.一种存储介质，其中存储有指令集，其特征在于，所述指令集用于执行：获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值。
4.根据权利要求3所述的一种存储介质，其特征在于，所述指令集还用于执行：所述不同量程通道数据包括：功率量程、中间量程和源量程；
所述“将所述不同量程通道数据转换至相同单位的数值”，还包括步骤：
中间量程功率量程转换系数为：IRC2PRC＝PRC/IRC；
源量程中间量程转换系数为：SRC2IRC＝IRC/SRC；
初始化所述IRC2PRC和所述SRC2IRC；
所述PRC为功率量程测量数据，所述IRC为中间量程测量数据，所述SRC为源量程测量数据。
5.一种实时在线反应性仪，其特征在于，所述反应性仪用于：获取不同量程通道数据，将所述不同量程通道数据转换至相同单位的数值。
6.根据权利要求5所述的一种实时在线反应性仪，其特征在于，所述反应性仪还用于：所述不同量程通道数据包括：功率量程、中间量程和源量程；
所述“将所述不同量程通道数据转换至相同单位的数值”，还包括步骤：
中间量程功率量程转换系数为：IRC2PRC＝PRC/IRC；
源量程中间量程转换系数为：SRC2IRC＝IRC/SRC；
初始化所述IRC2PRC和所述SRC2IRC；
所述PRC为功率量程测量数据，所述IRC为中间量程测量数据，所述SRC为源量程测量数据。
7.根据权利要求5所述的一种实时在线反应性仪，其特征在于，所述反应性仪通过通信连接计算机控制系统，通信连接方式包括以下中的一种或多种：使用modbus RTU协议、使用modbus TCP协议、使用modbus ASCII协议；
所述反应性仪用于：向所述计算机控制系统发送数据读取请求，并接收所述计算机控制系统返回的数据，根据所述返回的数据计算反应性，并发送计算结果至计算机控制系统进行输出显示。
8.根据权利要求5所述的一种实时在线反应性仪，其特征在于，所述反应性仪通过通信连接计算机控制系统，通信连接方式包括以下中的一种或多种：使用modbus RTU协议、使用modbus TCP协议、使用modbus ASCII协议；
所述反应性仪用于：接收所述计算机控制系统写入的数据，根据所述写入的数据计算反应性，并接收所述计算机控制系统的读取数据指令，返回计算结果至计算机控制系统进行输出显示。
9.根据权利要求5所述的一种实时在线反应性仪，其特征在于，所述反应性仪还用于：计算中子源强度；
所述“计算中子源强度”，包括步骤：
测量得深次临界稳定状态下的功率水平n；
计算在此深次临界状态下的keff或反应性ρ，及瞬发中子平均代时间l或瞬发中子寿命l0；
根据下式计算出中子源强度S：
Figure FDA0002630950720000031
或
Figure FDA0002630950720000032
所述反应性ρ定义如下：
Figure FDA0002630950720000033
所述keff本征值用来描述核反应堆系统的次临界度，所述次临界方程如下：
Figure FDA0002630950720000034
其中l0＝keff*l；
在Keff<0.98或ρ<-2000pcm时，中子源强度的测量误差符合工程应用需求。
10.根据权利要求5或6所述的一种实时在线反应性仪，其特征在于，
所述反应性仪通过串口与所述计算机控制系统进行通信，或所述反应性仪通过TCP/IP网络与所述计算机控制系统进行通信。
CN202010810975.6A 2020-08-13 2020-08-13 一种不同量程归一化方法、存储介质和实时在线反应性仪 Pending CN111799005A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010810975.6A CN111799005A (zh) 2020-08-13 2020-08-13 一种不同量程归一化方法、存储介质和实时在线反应性仪

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010810975.6A CN111799005A (zh) 2020-08-13 2020-08-13 一种不同量程归一化方法、存储介质和实时在线反应性仪

Publications (1)

Publication Number Publication Date
CN111799005A true CN111799005A (zh) 2020-10-20

Family

ID=72834593

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010810975.6A Pending CN111799005A (zh) 2020-08-13 2020-08-13 一种不同量程归一化方法、存储介质和实时在线反应性仪

Country Status (1)

Country Link
CN (1) CN111799005A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112687411A (zh) * 2020-12-24 2021-04-20 中国核动力研究设计院 一种基于多探测器接力信号的反应性测量方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106411870A (zh) * 2016-09-20 2017-02-15 中国核动力研究设计院 一种数字反应性仪

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106411870A (zh) * 2016-09-20 2017-02-15 中国核动力研究设计院 一种数字反应性仪

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘强: "核电厂中间量程探测器信号处理", 《仪器仪表用户》 *
潘泽飞: "《压水堆核动力电厂反应堆物理实验方法》", 31 January 2016 *
蒋国元: "《WWER-1000核电厂基础知识》", 30 September 2009 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112687411A (zh) * 2020-12-24 2021-04-20 中国核动力研究设计院 一种基于多探测器接力信号的反应性测量方法

Similar Documents

Publication Publication Date Title
CN113904375B (zh) 一种新能源并网系统电压支撑强度评估方法及系统
US11169189B2 (en) Systems and methods for operating generators based on generator steady state stability limits
CN111799005A (zh) 一种不同量程归一化方法、存储介质和实时在线反应性仪
CN108877969B (zh) 核功率理论模型建立及验证方法、系统及终端设备
CN113158589B (zh) 一种电池管理系统的仿真模型校准方法以及装置
WO2024125194A1 (zh) 反应堆一回路中流量确定方法、装置和计算机设备
CN111916238A (zh) 一种实时在线反应性仪
CN111797362B (zh) 一种中子源强度计算方法、存储介质和实时在线反应性仪
CN115267641B (zh) 同塔双回输电线路中电流互感器误差异常识别方法、系统
Arnold et al. Automated equilibrium titrator based on a personal computer
CN115144807B (zh) 差分滤噪和载流分级的电流互感器在线评估方法、装置
CN111309576A (zh) 一种压力测试方法和装置
Chen et al. Generalised‐fast decoupled state estimator
CN116070520A (zh) 水流阻力预测模型的构建方法、流量预测方法及装置
CN115296905A (zh) 一种基于移动终端的数据采集分析方法及系统
CN114708081A (zh) 信用风险预测方法及装置、电子设备和可读存储介质
Lin et al. Predictive capability and maturity assessment with bayesian network
CN110378041B (zh) 一种电网状态估计方法、系统、设备及计算机存储介质
Kim et al. Uncertainty estimation of the ADCP velocity measurements from the moving vessel method,(I) development of the framework
Weese et al. Deconvolution of neutron scattering data: a new computational approach
JP2006090944A (ja) 未臨界度評価装置及び未臨界度評価方法、並びに未臨界度評価用プログラム
CN117936134A (zh) 功率探测器的增益系数修正方法、装置、设备及介质
CN115469136A (zh) 一种基于电场传感器阵列的非侵入式三相电压测量方法
CN116430296A (zh) 电能表计量结果准确性的检验方法以及装置、介质、终端
CN116972781A (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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20230406

Address after: Room 504, Building 4, Jinhong Jiayuan, Youxi County, Sanming City, Fujian Province, 365100

Applicant after: Cai Yueqing

Address before: 301A, 3rd Floor, Building 1, No. 2, Chuangye Road, Shangjie Town, Minhou County, Fuzhou City, Fujian Province, 350000

Applicant before: Fuzhou qizhenggu Technology Co.,Ltd.