CN108876175B - 一种核电厂安全壳泄漏率计算评估方法及系统 - Google Patents

一种核电厂安全壳泄漏率计算评估方法及系统 Download PDF

Info

Publication number
CN108876175B
CN108876175B CN201810685983.5A CN201810685983A CN108876175B CN 108876175 B CN108876175 B CN 108876175B CN 201810685983 A CN201810685983 A CN 201810685983A CN 108876175 B CN108876175 B CN 108876175B
Authority
CN
China
Prior art keywords
containment
uncertainty
air
moment
leakage rate
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
CN201810685983.5A
Other languages
English (en)
Other versions
CN108876175A (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 General Nuclear Power Corp
China Nuclear Power Engineering Co Ltd
CGN Power Co Ltd
Shenzhen China Guangdong Nuclear Engineering Design Co Ltd
Original Assignee
China General Nuclear Power Corp
China Nuclear Power Engineering Co Ltd
CGN Power Co Ltd
Shenzhen China Guangdong Nuclear Engineering Design 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 China General Nuclear Power Corp, China Nuclear Power Engineering Co Ltd, CGN Power Co Ltd, Shenzhen China Guangdong Nuclear Engineering Design Co Ltd filed Critical China General Nuclear Power Corp
Priority to CN201810685983.5A priority Critical patent/CN108876175B/zh
Publication of CN108876175A publication Critical patent/CN108876175A/zh
Application granted granted Critical
Publication of CN108876175B publication Critical patent/CN108876175B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0635—Risk analysis of enterprise or organisation activities
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Physics & Mathematics (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Measuring Volume Flow (AREA)
  • Examining Or Testing Airtightness (AREA)

Abstract

本发明公开了一种核电厂安全壳泄漏率计算评估方法，包括：获取计算安全壳泄漏率的A类评定不确定度；根据安全壳内传感器实时测量的数据，获取计算安全壳泄漏率的B类评定不确定度；将所述A类评定不确定度和所述B类评定不确定度合成为安全壳泄漏率不确定度，以根据所述安全壳泄漏率不确定度对安全壳泄漏率的计算进行评估。本发明还公开了一种核电厂安全壳泄漏率计算评估系统。本发明能够提高评估的可靠性和准确度。

Description

一种核电厂安全壳泄漏率计算评估方法及系统
技术领域
本发明涉及核电站技术领域，尤其涉及一种核电厂安全壳泄漏率计算评估方法及系统。
背景技术
国际上普遍采用质量守恒法计算安全壳泄漏率，根据理想气体状态方程，空气标准体积(0、℃1.01325×105Pa下的体积)等效于空气质量，采用一天内若干标准体积变化量DVH与时间进行线性拟合，直线斜率即为当天的泄漏率Qld。安全壳泄漏率Qld用其最佳估计值和不确定度来表征，安全壳泄漏率不确定度是衡量泄漏率计算结果可靠性的重要指标。
统计学中安全壳泄漏率不确定度由A类评定不确定度和B类评定不确定度合成，根据被测量的一系列测得值得到的实验标准偏差称为A类评定标准不确定度，根据有关信息估计的先验概率分布得到的标准偏差估计值称为B类评定标准不确定度。
现有技术普遍采用一种已知的泄漏率监测软件对泄漏率计算进行评估，软件中自行定义了A类评定不确定度计算方法，但其与统计学中规定的斜率不确定度计算方法并不一致，不能判断其合理性，导致这种评估方法的可靠性低。对于B类评定标准不确定度，软件直接给出一个固定值，导致这种评估方法准确度低。
发明内容
本发明针对现有技术中存在的问题，提供了一种核电厂安全壳泄漏率计算评估方法及系统，能够提高评估的可靠性和准确度。
本发明就上述技术问题而提出的技术方案如下：
一方面，本发明提供一种核电厂安全壳泄漏率计算评估方法，包括：
获取安全壳泄漏率计算的A类评定不确定度；
根据安全壳内传感器实时测量的数据，获取安全壳泄漏率计算的B类评定不确定度；
将所述A类评定不确定度和所述B类评定不确定度合成为安全壳泄漏率不确定度，以根据所述安全壳泄漏率不确定度对安全壳泄漏率的计算进行评估。
进一步地，所述安全壳泄漏率是对每间隔预设时长的安全壳空气标准体积变化量进行线性拟合后的直线斜率；
所述安全壳泄漏率的计算公式如下：
Figure BDA0001711647920000021
Figure BDA0001711647920000022
Xi＝ti；
其中，Qld为安全壳泄漏率，ti为间隔i个预设时长后的时刻，t0为参考时刻，Yi为ti时刻的安全壳空气标准体积变化量，VH(ti)为ti时刻的安全壳空气标准体积，VH(t0)为t0时刻的安全壳空气标准体积，Qsarj为tj时刻测量到的注入安全壳的压缩空气流量，N为测量次数。
进一步地，所述A类评定不确定度的计算公式为：
Figure BDA0001711647920000023
Figure BDA0001711647920000024
其中，uA(Qld)为A类评定不确定度，
Figure BDA0001711647920000033
为所有Xi的平均值，σ2为Yi围绕拟合直线的方差，a1为拟合直线的斜率，a0为拟合直线的截距。
进一步地，所述传感器包括压力传感器、流量传感器和至少一个温度传感器；
所述根据安全壳内传感器实时测量的数据，获取安全壳泄漏率计算的B类评定不确定度，具体包括：
根据所有温度传感器在ti时刻测量的安全壳空气温度，计算获得ti时刻的安全壳空气平均温度不确定度；
根据所述压力传感器在ti时刻测量的安全壳空气压力和ti时刻的安全壳空气平均温度不确定度，计算获得ti时刻的安全壳空气标准体积不确定度；
根据所述流量传感器在ti时刻测量的压缩空气流量，计算获得ti时刻的流量传感器测量不确定度；
根据ti时刻的安全壳空气标准体积不确定度和流量传感器测量不确定度，计算获得安全壳泄漏率的B类评定不确定度。
进一步地，所述根据所有温度传感器在ti时刻测量的安全壳空气温度，计算获得ti时刻的安全壳空气平均温度不确定度，具体包括：
根据所有温度传感器在ti时刻所测量的安全壳空气温度，计算获得ti时刻的安全壳空气平均温度和每个温度传感器的测量不确定度；
根据ti时刻的安全壳空气平均温度和每个温度传感器的测量不确定度，计算获得ti时刻的安全壳空气平均温度不确定度；
其中，当前时刻的安全壳空气平均温度的计算公式为：
Figure BDA0001711647920000031
当前时刻的安全壳平均温度不确定度的计算公式为：
Figure BDA0001711647920000032
其中，Tavg为安全壳空气平均温度，VL为安全壳内自由空间体积，Vm为第m个温度传感器代表的自由空间体积，tmavg为第m个温度传感器测量的安全壳空气温度，n为温度传感器的个数，uB(Tavg)为安全壳平均温度不确定度，uB(tmavg)为第m个温度传感器的测量不确定度。
进一步地，所述根据所述压力传感器在ti时刻测量的安全壳空气压力和ti时刻的安全壳空气平均温度不确定度，计算获得ti时刻的安全壳空气标准体积不确定度，具体包括：
根据所述压力传感器在ti时刻测量的安全壳空气压力，计算ti时刻的压力传感器测量不确定度；
根据所述压力传感器在ti时刻测量的安全壳空气压力和ti时刻的安全壳空气平均温度，计算获得ti时刻的安全壳空气标准体积；
根据ti时刻的压力传感器测量不确定度、安全壳空气标准体积和安全壳平均温度不确定度，计算获得ti时刻的安全壳空气标准体积不确定度；
其中，所述ti时刻的安全壳空气标准体积的计算公式如下：
Figure BDA0001711647920000041
所述ti时刻的安全壳空气标准体积不确定度的计算公式如下：
Figure BDA0001711647920000042
其中，TN为标准状态温度，PN为标准状态绝对压力，piavg为压力传感器在ti时刻测量的安全壳空气压力，Tiavg为ti时刻的安全壳空气平均温度，uB(VH(ti))为ti时刻的安全壳空气标准体积不确定度，uB(piavg)为ti时刻的压力传感器测量不确定度，uB(Tiavg)为ti时刻的安全壳空气平均温度不确定度。
进一步地，所述根据ti时刻的安全壳空气标准体积不确定度和流量传感器测量不确定度，计算获得安全壳泄漏率的B类评定不确定度，具体包括：
设置所述安全壳泄漏率Qld＝Qlda-Qldb；其中，
Figure BDA0001711647920000043
根据ti时刻的安全壳空气标准体积不确定度，计算获得Qlda不确定度；
根据ti时刻的流量传感器测量不确定度，计算获得Qldb不确定度；
根据所述Qlda不确定度和所述Qldb不确定度，计算获得安全壳泄漏率的B类评定不确定度；
其中，所述Qlda不确定度的计算公式如下：
Figure BDA0001711647920000051
所述Qldb不确定度的计算公式如下：
Figure BDA0001711647920000052
所述B类评定不确定度的计算公式如下：
Figure BDA0001711647920000053
其中，uB(Qlda)为Qlda不确定度，uB(Qldb)为Qldb不确定度，uB(VH(t0))为t0时刻的安全壳空气标准体积不确定度，uB(Qsarj)为tj时刻的流量传感器测量不确定度，uB(Qld)为B类评定不确定度。
进一步地，所述安全壳泄漏率不确定度的合成公式如下：
Figure BDA0001711647920000054
其中，u(Qld)为安全壳泄漏率不确定度，uA(Qld)为A类评定不确定度，uB(Qld)为B类评定不确定度。
进一步地，所述根据所述安全壳泄漏率不确定度对安全壳泄漏率的计算进行评估，具体包括：
若所述安全壳泄漏率不确定度低于预设阈值，则评估安全壳泄漏率的计算结果可靠；
若所述安全壳泄漏率不确定度高于预设阈值，则评估安全壳泄漏率的计算结果不可靠，对所述安全壳泄漏率的计算数据进行处理并重新计算。
另一方面，本发明提供一种实现上述核电厂安全壳泄漏率计算评估方法的系统，所述系统包括：
A类评定不确定度获取模块，用于获取安全壳泄漏率计算的A类评定不确定度；
B类评定不确定度获取模块，用于根据安全壳内传感器实时测量的数据，获取安全壳泄漏率计算的B类评定不确定度；
合成评估模块，用于将所述A类评定不确定度和所述B类评定不确定度合成为安全壳泄漏率不确定度，以根据所述安全壳泄漏率不确定度对安全壳泄漏率的计算进行评估。
本发明实施例提供的技术方案带来的有益效果是：
基于统计学理论获取计算安全壳泄漏率的A类评定不确定度，使A类评定不确定度的获取有据可循，提高安全壳泄漏率计算评估的可靠性，根据安全壳内传感器实时测量的数据获取计算安全壳泄漏率的B类评定不确定度，实现动态评估，提高安全壳泄漏率计算评估的可靠性。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例一提供的核电厂安全壳泄漏率计算评估方法的流程示意图；
图2是本发明实施例二提供的核电厂安全壳泄漏率计算评估系统的结构示意图。
具体实施方式
为了解决现有技术在核电厂安全壳泄漏率计算评估的可靠性和准确度低等技术问题，本发明旨在提供一种核电厂安全壳泄漏率计算评估方法，其核心思想是：基于统计学理论，获取计算安全壳泄漏率的A类评定不确定度；根据安全壳内传感器实时测量的数据，获取计算安全壳泄漏率的B类评定不确定度；将所述A类评定不确定度和所述B类评定不确定度合成为安全壳泄漏率不确定度，以根据所述安全壳泄漏率不确定度对安全壳泄漏率的计算进行评估。本发明所提供的核电厂安全壳泄漏率计算评估方法能够有效提高评估的可靠性和准确度。
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
实施例一
本发明实施例提供了一种核电厂安全壳泄漏率计算评估方法，如图1所示，所述方法包括：
S1、获取安全壳泄漏率计算的A类评定不确定度；
S2、根据安全壳内传感器实时测量的数据，获取安全壳泄漏率计算的B类评定不确定度；
S3、将所述A类评定不确定度和所述B类评定不确定度合成为安全壳泄漏率不确定度，以根据所述安全壳泄漏率不确定度对安全壳泄漏率的计算进行评估。
需要说明的是，在机组运行期间，需对核电厂安全壳泄漏率进行实时监测，以监视安全壳密封性变化，而安全壳泄漏率的监测是通过对安全壳泄漏率进行计算来实现的，但安全壳泄漏率的计算结果是否可靠还需通过安全壳泄漏率不确定度进行评估，其中安全壳泄漏率的计算方法为现有技术中常用的方法。而安全壳泄漏率不确定度是通过A类评定不确定度和B类评定不确定度合成所得，其中，A类评定不确定度基于统计学原理来获取，B类评定不确定度通过安全壳内传感器实时测量的数据来获取，有效提高安全壳泄漏率不确定度计算的可靠性和准确度，进而提高安全壳泄漏率计算评估的可靠性和准确度。
具体的，步骤S1中的A类评定不确定度的计算方法如下：
假设安全壳泄漏率Qld根据一天内48个空气标准体积变化量DVH线性拟合而得，直线斜率即为泄漏率Qld。Qld的A类评定不确定度来自DVH点的分散度，由安全壳内热力条件的随机波动所致，表示为uA(Qld)。
令Xi＝ti，i＝1…N,N＝48；ti＝0h，0.5h，1h，1.5h...23.5h；
Yi＝DVH(ti)；
Figure BDA0001711647920000081
DVH(ti)：安全壳空气标准体积变化量(Nm3)；
VH(ti)：当前时刻ti安全壳空气标准体积(Nm3)；
VH(t0)：参考时刻t0安全壳空气标准体积(Nm3)；
Qsar：注入安全壳的压缩空气流量(Nm3/h)，由压缩空气注入所致安全壳空气体积变化对泄漏监测产生干扰，需扣除从参考时刻至当前时刻注入安全壳的压缩空气总体积。
其中，VH(ti)、VH(t0)、Qsar根据设置在安全壳内的温度传感器、压力传感器、流量传感器的测量数据而获得。
根据统计学理论，安全壳泄漏率Qld的A类评定不确定度计算如下：
Figure BDA0001711647920000082
其中，σ2为Yi围绕拟合直线的方差(又称残余方差)，表示为：
Figure BDA0001711647920000083
其中，
Figure BDA0001711647920000084
为所有Xi的平均值，a1，a0分别为拟合直线的斜率和截距。
本实施例中，A类评定不确定度的计算方法与统计学中规定的斜率不确定度计算方法一致，使A类评定不确定度的计算有据可循，提高A类评定不确定度计算的可靠性。
具体的，在步骤S2中，所述根据安全壳内传感器实时测量的数据，获取计算安全壳泄漏率的B类评定不确定度，具体包括：
根据所有温度传感器在ti时刻测量的安全壳空气温度，计算获得ti时刻的安全壳空气平均温度不确定度；
根据所述压力传感器在ti时刻测量的安全壳空气压力和ti时刻的安全壳空气平均温度不确定度，计算获得ti时刻的安全壳空气标准体积不确定度；
根据所述流量传感器在ti时刻测量的压缩空气流量，计算获得ti时刻的流量传感器测量不确定度；
根据ti时刻的安全壳空气标准体积不确定度和流量传感器测量不确定度，计算获得安全壳泄漏率的B类评定不确定度。
进一步地，所述根据所有温度传感器在ti时刻测量的安全壳空气温度，计算获得ti时刻的安全壳空气平均温度不确定度，具体包括：
根据所有温度传感器在ti时刻所测量的安全壳空气温度，计算获得ti时刻的安全壳空气平均温度和每个温度传感器的测量不确定度；
根据ti时刻的安全壳空气平均温度和每个温度传感器的测量不确定度，计算获得ti时刻的安全壳空气平均温度不确定度。
在实际应用中，假设系统使用n个温度传感器测量安全壳空气温度，且每半小时测量一次，即预设时长一般设置为半小时，每个温度传感器的测量值代表一部分自由空间空气平均温度，则安全壳空气平均温度Tavg由每个传感器所代表的自由空间体积加权计算而得，计算模型如下：
Figure BDA0001711647920000091
VL：安全壳内自由空间体积，单位m3；
Vm：第m个温度传感器代表的自由空间体积，单位m3；
tmavg：第m个温度传感器测量的安全壳空气温度，单位K。
由于安全壳内各温度传感器的测量值互不相关，任何两个温度传感器的测量结果不确定度的相关系数为零，根据不确定度合成理论，每半小时安全壳空气平均温度不确定度合成公式如下：
Figure BDA0001711647920000101
其中，uB(tmavg)：第m个温度传感器的测量不确定度
将式(6)代入式(7)可得：
Figure BDA0001711647920000102
令：
Figure BDA0001711647920000103
则式(8)可表示为：
Figure BDA0001711647920000104
根据式(10)每半小时计算一次安全壳空气平均温度不确定度，则一天内可计算出48个安全壳空气平均温度不确定度。
进一步地，所述根据所述压力传感器在ti时刻测量的安全壳空气压力和ti时刻的安全壳空气平均温度不确定度，计算获得ti时刻的安全壳空气标准体积不确定度，具体包括：
根据所述压力传感器在ti时刻测量的安全壳空气压力，计算ti时刻的压力传感器测量不确定度；
根据所述压力传感器在ti时刻测量的安全壳空气压力和ti时刻的安全壳空气平均温度，计算获得ti时刻的安全壳空气标准体积；
根据ti时刻的压力传感器测量不确定度、安全壳空气标准体积和安全壳平均温度不确定度，计算获得ti时刻的安全壳空气标准体积不确定度。
需要说明的是，根据压力传感器在ti时刻测量的安全壳空气压力piavg和ti时刻的安全壳空气平均温度Tiavg，计算获得ti时刻的安全壳空气标准体积VH(ti)，其计算模型如下：
Figure BDA0001711647920000105
其中，TN：标准状态温度，273.15K；
PN：标准状态绝对压力，1.01325×105Pa；
VL：安全壳自由空间体积，单位m3；
piavg：压力传感器在ti时刻测量的安全壳空气压力，单位Pa；
Tiavg：ti时刻的安全壳空气平均温度，单位K；
i：测量次数，一般取值为1至48。
令
Figure BDA0001711647920000111
则式(11)可表示为：
Figure BDA0001711647920000112
由于piavg和Tiavg不相关，根据不确定度合成理论，ti时刻的安全壳空气标准体积不确定度计算如下：
Figure BDA0001711647920000113
其中，uB(piavg)：ti时刻的压力传感器测量不确定度；
uB(Tiavg)：ti时刻的安全壳空气平均温度不确定度。
根据式(14)每半小时计算一次安全壳空气标准体积不确定度，则一天内可计算出48个安全壳空气标准体积不确定度。
进一步地，所述根据ti时刻的安全壳空气标准体积不确定度和流量传感器测量不确定度，计算获得安全壳泄漏率的B类评定不确定度，具体包括：
设置所述安全壳泄漏率Qld＝Qlda-Qldb；
根据ti时刻的安全壳空气标准体积不确定度，计算获得Qlda不确定度；
根据ti时刻的流量传感器测量不确定度，计算获得Qldb不确定度；
根据所述Qlda不确定度和所述Qldb不确定度，计算获得安全壳泄漏率的B类评定不确定度。
需要说明的是，根据最小二乘法计算直线斜率原理，安全壳泄漏率Qld的计算如下：
Figure BDA0001711647920000121
其中，Xi：间隔i个预设时长后的时刻，一般表示一天内48个半小时时刻，取值为0h,0.5h,1h,1.5h…23.5h；
Yi：ti时刻的安全壳空气标准体积相对t0时刻的变化量DVH(ti)；
N：测量次数，一般为常数48，表示一天内半小时数据点数量。
由于N、Xi为常数，则令：
Figure BDA0001711647920000122
Figure BDA0001711647920000123
则式(15)可表示为：
Figure BDA0001711647920000124
安全壳空气标准体积变化量Yi为ti时刻相对t0时刻的安全壳空气体积之差，并扣除t0至ti时刻注入安全壳的压缩空气体积，计算如下：
Figure BDA0001711647920000125
将式(19)代入式(18)，安全壳泄漏率Qld可表示为：
Figure BDA0001711647920000126
令：
Figure BDA0001711647920000127
Figure BDA0001711647920000128
由于Qlda根据安全壳内压力、温度计算而得，Qldb根据压缩空气流量计算而得，三种数据互不相关，因此Qld的B类评定标准不确定度计算如下：
Figure BDA0001711647920000129
其中，uB(Qlda)为Qlda不确定度，uB(Qlda)为Qlda不确定度。
在计算Qlda不确定度时，由于半小时空气标准体积采用同一组温度、压力传感器测量数据计算而得，因此49个空气标准体积数据中任何两个数据均成强相关性，保守估计，令相关系数为1，则uB(Qlda)计算如下：
Figure BDA0001711647920000131
在计算Qlda不确定度时，由于从参考时刻t0至当前时刻ti的压缩空气流量数据采用同一流量传感器测量而得，任何两个流量数据成强相关性，保守估计，令相关系数为1，则uB(Qlda)计算如下：
Figure BDA0001711647920000132
其中，uB(Qsarj)为tj时刻的流量传感器测量不确定度。
需要说明的是，安全壳泄漏率计算需使用安全壳内温度传感器、压力传感器、压缩空气流量传感器测量数据，而传感器测量不确定度的计算方法可根据样本资料获取或向厂家咨询。获取传感器在ti时刻测量的数据，结合传感器测量不确定度的计算方法，即可计算出传感器测量不确定度。
本实施例采用不同于A类评定的方法对测量不确定度分量进行的评定称为B类评定，本实施例对安全壳泄漏率计算过程进行分解，从仪表不确定度出发，分析不确定度在每个环节的计算方法，确定安全壳泄漏率的B类评定不确定度评估方法，提高评估的准确度。
具体的，在步骤S3中，根据不确定度合成原理，安全壳泄漏率不确定度由A类评定不确定度和B类评定不确定度合成而得，由于两类不确定度分量互不相关，安全壳泄漏率不确定度的合成公式如下：
Figure BDA0001711647920000133
进一步地，在步骤S3中，所述根据所述安全壳泄漏率不确定度对安全壳泄漏率的计算进行评估，具体包括：
若所述安全壳泄漏率不确定度低于预设阈值，则评估安全壳泄漏率的计算结果可靠；
若所述安全壳泄漏率不确定度高于预设阈值，则评估安全壳泄漏率的计算结果不可靠，对所述安全壳泄漏率的计算数据进行处理并重新计算。
本发明实施例基于统计学理论获取计算安全壳泄漏率的A类评定不确定度，使A类评定不确定度的获取有据可循，提高安全壳泄漏率计算评估的可靠性，根据安全壳内传感器实时测量的数据获取计算安全壳泄漏率的B类评定不确定度，实现动态评估，提高安全壳泄漏率计算评估的可靠性。
实施例二
本发明实施例提供了一种核电厂安全壳泄漏率计算评估系统，能够实现上述核电厂安全壳泄漏率计算评估方法的所有流程，如图2所示，所述系统包括：
A类评定不确定度获取模块1，用于获取计算安全壳泄漏率的A类评定不确定度；
B类评定不确定度获取模块2，用于根据安全壳内传感器实时测量的数据，获取计算安全壳泄漏率的B类评定不确定度；
合成评估模块3，用于将所述A类评定不确定度和所述B类评定不确定度合成为安全壳泄漏率不确定度，以根据所述安全壳泄漏率不确定度对安全壳泄漏率的计算进行评估。
本发明实施例基于统计学理论获取计算安全壳泄漏率的A类评定不确定度，使A类评定不确定度的获取有据可循，提高安全壳泄漏率计算评估的可靠性，根据安全壳内传感器实时测量的数据获取计算安全壳泄漏率的B类评定不确定度，实现动态评估，提高安全壳泄漏率计算评估的可靠性。
综上所述，本发明提出了一种核电厂安全壳泄漏率计算评估方法及系统，其具有较好的实用效果：根据统计学理论给出安全壳泄漏率的A类评定不确定度计算方法，使评价方法有据可循；根据统计学理论系统分析B类评定不确定度在安全壳泄漏率各个计算环节的评价方法，填补国内的研究空白；安全壳泄漏率B类评定不确定度采用实时动态评估方法，比采用静态值更为准确。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种核电厂安全壳泄漏率计算评估方法，其特征在于，包括：
获取安全壳泄漏率计算的A类评定不确定度；
所述A类评定不确定度的计算公式为：
Figure FDA0003232282130000011
Figure FDA0003232282130000012
其中，Xi＝ti，ti为间隔i个预设时长后的时刻，t0为参考时刻，Yi为ti时刻的安全壳空气标准体积变化量，N为测量次数，
uA(Qld)为A类评定不确定度，
Figure FDA0003232282130000013
为所有Xi的平均值，σ2为Yi围绕拟合直线的方差，a1为拟合直线的斜率，a0为拟合直线的截距；
根据安全壳内传感器实时测量的数据，获取安全壳泄漏率计算的B类评定不确定度，具体包括：
根据所有温度传感器在ti时刻测量的安全壳空气温度，计算获得ti时刻的安全壳空气平均温度不确定度；
根据压力传感器在ti时刻测量的安全壳空气压力和ti时刻的安全壳空气平均温度不确定度，计算获得ti时刻的安全壳空气标准体积不确定度；
根据流量传感器在ti时刻测量的压缩空气流量，计算获得ti时刻的流量传感器测量不确定度；
根据ti时刻的安全壳空气标准体积不确定度和流量传感器测量不确定度，计算获得安全壳泄漏率的B类评定不确定度；
将所述A类评定不确定度和所述B类评定不确定度合成为安全壳泄漏率不确定度，以根据所述安全壳泄漏率不确定度对安全壳泄漏率的计算进行评估；所述根据所述安全壳泄漏率不确定度对安全壳泄漏率的计算进行评估，具体包括：若所述安全壳泄漏率不确定度低于预设阈值，则评估安全壳泄漏率的计算结果可靠；若所述安全壳泄漏率不确定度高于预设阈值，则评估安全壳泄漏率的计算结果不可靠，对所述安全壳泄漏率的计算数据进行处理并重新计算。
2.如权利要求1所述的核电厂安全壳泄漏率计算评估方法，其特征在于，所述安全壳泄漏率是对每间隔预设时长的安全壳空气标准体积变化量进行线性拟合后的直线斜率；
所述安全壳泄漏率的计算公式如下：
Figure FDA0003232282130000021
Figure FDA0003232282130000022
Xi＝ti；
其中，Qld为安全壳泄漏率，ti为间隔i个预设时长后的时刻，t0为参考时刻，Yi为ti时刻的安全壳空气标准体积变化量，VH(ti)为ti时刻的安全壳空气标准体积，VH(t0)为t0时刻的安全壳空气标准体积，Qsarj为tj时刻测量到的注入安全壳的压缩空气流量，N为测量次数。
3.如权利要求1所述的核电厂安全壳泄漏率计算评估方法，其特征在于，所述根据所有温度传感器在ti时刻测量的安全壳空气温度，计算获得ti时刻的安全壳空气平均温度不确定度，具体包括：
根据所有温度传感器在ti时刻所测量的安全壳空气温度，计算获得ti时刻的安全壳空气平均温度和每个温度传感器的测量不确定度；
根据ti时刻的安全壳空气平均温度和每个温度传感器的测量不确定度，计算获得ti时刻的安全壳空气平均温度不确定度；
其中，当前时刻的安全壳空气平均温度的计算公式为：
Figure FDA0003232282130000031
当前时刻的安全壳平均温度不确定度的计算公式为：
Figure FDA0003232282130000032
其中，Tavg为安全壳空气平均温度，VL为安全壳内自由空间体积，Vm为第m个温度传感器代表的自由空间体积，tmavg为第m个温度传感器测量的安全壳空气温度，n为温度传感器的个数，uB(Tavg)为安全壳平均温度不确定度，uB(tmavg)为第m个温度传感器的测量不确定度。
4.如权利要求3所述的核电厂安全壳泄漏率计算评估方法，其特征在于，所述根据所述压力传感器在ti时刻测量的安全壳空气压力和ti时刻的安全壳空气平均温度不确定度，计算获得ti时刻的安全壳空气标准体积不确定度，具体包括：
根据所述压力传感器在ti时刻测量的安全壳空气压力，计算ti时刻的压力传感器测量不确定度；
根据所述压力传感器在ti时刻测量的安全壳空气压力和ti时刻的安全壳空气平均温度，计算获得ti时刻的安全壳空气标准体积；
根据ti时刻的压力传感器测量不确定度、安全壳空气标准体积和安全壳平均温度不确定度，计算获得ti时刻的安全壳空气标准体积不确定度；
其中，所述ti时刻的安全壳空气标准体积的计算公式如下：
Figure FDA0003232282130000033
所述ti时刻的安全壳空气标准体积不确定度的计算公式如下：
Figure FDA0003232282130000034
其中，TN为标准状态温度，PN为标准状态绝对压力，piavg为压力传感器在ti时刻测量的安全壳空气压力，Tiavg为ti时刻的安全壳空气平均温度，uB(VH(ti))为ti时刻的安全壳空气标准体积不确定度，uB(piavg)为ti时刻的压力传感器测量不确定度，uB(Tiavg)为ti刻的安全壳空气平均温度不确定度。
5.如权利要求4所述的核电厂安全壳泄漏率计算评估方法，其特征在于，所述根据ti时刻的安全壳空气标准体积不确定度和流量传感器测量不确定度，计算获得安全壳泄漏率的B类评定不确定度，具体包括：
设置所述安全壳泄漏率Qld＝Qlda-Qldb；其中，
Figure FDA0003232282130000041
根据ti时刻的安全壳空气标准体积不确定度，计算获得Qlda不确定度；
根据ti时刻的流量传感器测量不确定度，计算获得Qldb不确定度；
根据所述Qlda不确定度和所述Qldb不确定度，计算获得安全壳泄漏率的B类评定不确定度；
其中，所述Qlda不确定度的计算公式如下：
Figure FDA0003232282130000042
所述Qlab不确定度的计算公式如下：
Figure FDA0003232282130000043
所述B类评定不确定度的计算公式如下：
Figure FDA0003232282130000044
其中，uB(Qlda)为Qlda不确定度，uB(Qldb)为Qldb不确定度，uB(VH(t0))为t0时刻的安全壳空气标准体积不确定度，uB(Qsarj)为tj时刻的流量传感器测量不确定度，uB(Qld)为B类评定不确定度。
6.如权利要求1所述的核电厂安全壳泄漏率计算评估方法，其特征在于，所述安全壳泄漏率不确定度的合成公式如下：
Figure FDA0003232282130000051
其中，u(Qld)为安全壳泄漏率不确定度，uA(Qld)为A类评定不确定度，uB(Qld)为B类评定不确定度。
7.一种实现如权利要求1至6任一项所述的核电厂安全壳泄漏率计算评估方法的系统，其特征在于，所述系统包括：
A类评定不确定度获取模块，用于获取安全壳泄漏率计算的A类评定不确定度；
B类评定不确定度获取模块，用于根据安全壳内传感器实时测量的数据，获取安全壳泄漏率计算的B类评定不确定度；
合成评估模块，用于将所述A类评定不确定度和所述B类评定不确定度合成为安全壳泄漏率不确定度，以根据所述安全壳泄漏率不确定度对安全壳泄漏率的计算进行评估。
CN201810685983.5A 2018-06-28 2018-06-28 一种核电厂安全壳泄漏率计算评估方法及系统 Active CN108876175B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810685983.5A CN108876175B (zh) 2018-06-28 2018-06-28 一种核电厂安全壳泄漏率计算评估方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810685983.5A CN108876175B (zh) 2018-06-28 2018-06-28 一种核电厂安全壳泄漏率计算评估方法及系统

Publications (2)

Publication Number Publication Date
CN108876175A CN108876175A (zh) 2018-11-23
CN108876175B true CN108876175B (zh) 2021-12-03

Family

ID=64295490

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810685983.5A Active CN108876175B (zh) 2018-06-28 2018-06-28 一种核电厂安全壳泄漏率计算评估方法及系统

Country Status (1)

Country Link
CN (1) CN108876175B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110085335B (zh) * 2019-05-14 2022-03-08 广西防城港核电有限公司 安全壳泄漏率在线监测信号处理方法
CN110189840B (zh) * 2019-06-19 2021-11-16 广西防城港核电有限公司 核电站临界启机阶段安全壳泄漏率状态判断的方法
CN111189594B (zh) * 2020-01-15 2021-01-19 大连理工大学 一种基于流量守恒定理的多层承压壳泄漏率的快速评估方法
CN113639939B (zh) * 2021-07-09 2023-11-14 中国辐射防护研究院 一种建筑物内空气泄漏率异常变化的实时监测方法及系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102420021A (zh) * 2011-11-15 2012-04-18 中广核工程有限公司 一种核电站安全壳泄漏率监测方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6741955B2 (en) * 2002-03-05 2004-05-25 Uson, L.P. System and method for leak rate testing during adiabatic cooling

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102420021A (zh) * 2011-11-15 2012-04-18 中广核工程有限公司 一种核电站安全壳泄漏率监测方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
核电站安全壳泄漏率在线监测系统;陈亮 等;《中国仪器仪表》;20160229;第67-68页 *
核电站安全壳泄漏率计算模型与分析;李晓庚;《科技创新导》;20151231;第153-154页 *

Also Published As

Publication number Publication date
CN108876175A (zh) 2018-11-23

Similar Documents

Publication Publication Date Title
CN108876175B (zh) 一种核电厂安全壳泄漏率计算评估方法及系统
CN108388694A (zh) 一种塑封光电耦合器贮存寿命预测方法
US20130306382A1 (en) Weighing apparatus
CN101319940B (zh) 一种长杆铂电阻温度计的检测装置及其使用方法
CN105606501A (zh) 粉尘浓度检测方法及装置
CN117503204A (zh) 一种超声ai数字诊疗系统
CN106248167A (zh) 液位测量装置及利用该测量装置的测试方法
CN106568535B (zh) 一种机动车尾气排放温度传感器的精度检测系统
CN114061691A (zh) 一种航天供气系统用气量测量方法及系统
CN111352058A (zh) 用于对测量设备的电位传感器在过程中进行校准的方法
CN110865250A (zh) 融合电流监测的配电设备状态监测装置及发热检测方法
KR20190125872A (ko) 관측자료 품질검사 장치 및 이를 이용한 관측자료 품질검사 방법
Chen et al. The implementation of the data validation process in a gas turbine performance monitoring system
Lampasi et al. Generalized lambda distribution for the expression of measurement uncertainty
CN112016045B (zh) 纳秒级脉冲功率计数据处理方法
CN112345139A (zh) 一种基于电池箱体上的压力检测方法
Hashemian Applying online monitoring for nuclear power plant instrumentation and control
JP4786670B2 (ja) 計測器ドリフト検知装置
CN107796492B (zh) 一种超声波燃气表在线式校准方法
Menabde et al. Use of the theory of measurement uncertainty in procedures for data processing and results obtained by checking-calibration gas flow meters
CN110161996A (zh) 用于发电厂机组能耗分析的方法及系统
CN109754582A (zh) 一种燃气表的无线抄表方法
CN113237575B (zh) 一种基于双温度激励的温度传感器动态校准方法
Sparasci et al. Submillimetre thermistors for balloon‐borne applications up to lower stratosphere: preliminary characterization with 0.02 K uncertainty
CN108955743B (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