CN107941436B - 一种计算逆止阀泄漏率的方法及装置 - Google Patents

一种计算逆止阀泄漏率的方法及装置 Download PDF

Info

Publication number
CN107941436B
CN107941436B CN201610892667.6A CN201610892667A CN107941436B CN 107941436 B CN107941436 B CN 107941436B CN 201610892667 A CN201610892667 A CN 201610892667A CN 107941436 B CN107941436 B CN 107941436B
Authority
CN
China
Prior art keywords
closed space
check valve
pressure
leakage
tested
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
CN201610892667.6A
Other languages
English (en)
Other versions
CN107941436A (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
Fujian Ningde Nuclear Power Co Ltd
Original Assignee
China General Nuclear Power Corp
Fujian Ningde Nuclear Power 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, Fujian Ningde Nuclear Power Co Ltd filed Critical China General Nuclear Power Corp
Priority to CN201610892667.6A priority Critical patent/CN107941436B/zh
Publication of CN107941436A publication Critical patent/CN107941436A/zh
Application granted granted Critical
Publication of CN107941436B publication Critical patent/CN107941436B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00—Investigating fluid-tightness of structures
    • G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/26—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors
    • G01M3/28—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds
    • G01M3/2876—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds for valves

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Examining Or Testing Airtightness (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本申请提供了一种计算逆止阀的泄漏率的方法及装置，对于边界包括待测的逆止阀的封闭空间与外界未达到热平衡的情况，利用所述封闭空间的压力随时间的变化数据进行曲线拟合，并依据拟合出的曲线公式，排除所述封闭空间与外界的热量传递引起的封闭空间的水体积的变化因素，得到因所述待测的逆止阀的泄漏引起的所述封闭空间的压力变化速率，并由此得到所述待测的逆止阀的泄漏率。因为使用曲线拟合这一新的方式得到逆止阀的泄漏率可以排除所述封闭空间与外界的热量传递对测试结果的影响，所以使得测试结果比温度未稳定时获得的数据更接近真实值，避免测试结果出现负值导致结果无效或者测试结果偏大而导致结果不满意。

Description

一种计算逆止阀泄漏率的方法及装置
技术领域
本申请涉及自动控制领域，尤其涉及一种计算逆止阀泄漏率的方法及装置。
背景技术
对于核电站而言，保证核安全的三大功能(也称作三要素)包括：反应性控制、堆芯冷却和放射性产物的包容。其中，压水堆核电厂设置了三道屏障：燃料元件包壳、一回路压力边界和反应堆厂房(安全壳)。发生事故后，用来限制事故后果并将反应堆恢复到可靠的安全状态的系统称之为专设安全设施。安全注入系统(系统代码为RIS)就是其中之一。当发生一回路或二回路破口事件或事故时，RIS向一回路注入含硼水，除了保证堆芯的冷却外，同时保证堆芯处于次临界状态。安全注入系统包括下列三个子系统，其功能是在不同的压力条件下将含硼水注入反应堆冷却剂系统(RCP，是一回路的主系统)中。
如图1所示，安全注入系统的子系统中压安注系统主要由三个安注箱RIS001BA、RIS002 BA和RIS 003BA组成，分别接到RCP系统三个环路的冷管段上。安注箱内存CB≈2400ppm的含硼水，用压力约为4.2MPa.a的氮气覆盖。一回路大破口和中等破口事故时，当RCP压力降到安注箱压力以下时，由氮气压将含硼水注入RCP冷段，能在短时间内淹没堆芯，避免燃料棒熔化。
中压安注箱与RCP系统的每条连接管线上设有一个电动隔离阀(分别为RIS001VP/RIS002 VP/RIS003VP)和两道逆止阀(RIS004VP/RIS005 VP/RIS006VP和RCP121 VP/RCP221 VP/RCP 321VP)。在正常停堆期间，当一回路压力低于70bar.a时，关闭电动隔离阀RIS001VP、RIS002VP和RIS003VP，防止安注箱向RCP注入硼水。机组启动期间，一回路压力高于7.0bar.a后，电动隔离阀打开，由每条注入管线上的两个串联的逆止阀来保证一回路的高压 冷却剂不会进入安注箱。
为了验证逆止阀RCP121 VP、RCP221 VP和RCP 321 VP的密封性，将RRA系统与一回路隔离(RRA014/015VP关闭)并降温到60℃，一回路升压到70bar.a，一回路维持170℃左右。
打开阀门RS101VP、RIS119VP、RIS120VP、RIS 121 VP、RIS 122 VP和RIS 124VP，关闭阀门RIS113 VP、RIS114VP和RIS115VP、RIS125VP，与密封性完好的逆止阀RIS004 VP、RIS005 VP和RIS 006VP一起，在待测逆止阀RCP121 VP、RCP221 VP和RCP 321 VP上游形成一个封闭空间，即此封闭空间除待测逆止阀之外，其它边界的密封性均完好。验证开始时，封闭空间压力在42～48bar.g之间，假定水不可压缩的前提下，如果待测的逆止阀有泄漏，封闭空间的压力就会出现上涨，且其上涨速率与泄漏率大小有关，可通过封闭空间的压力表080LP反映压力上涨速率。
使用上述测量方法实际验证验证逆止阀的密封性时，常常会得到为负的泄漏率，如表1和表2所示。
表1
(L核电厂历史测试结果)
测试时间 L3RCP121VP L3RCP221VP L3RCP321VP
2010.06.07 -11.76 -11.76 -11.76
2010.10.27 -12 4.5 -18
2012.08.01 -10.2 -10.2 -10.2
2012.08.31 -12 1.2 -24
L4RCP121VP L4RCP221VP L4RCP321VP
2011.02.03 -24 -24 -24
2011.02.23 -43.2 NA -48
2011.06.23 -10.5 -10.5 -10.5
2011.07.19 -12 -63.5 -13.2
2012.07.11 -33 -33 -33
表2
(宁德核电厂历史测试结果)
测试时间 N1RCP121VP N1RCP221VP N1RCP321VP
N101 -11.76 -11.76 -11.76
N102 -12 -12 -12
从物理意义上说，逆止阀泄漏率是负值意味着被测的逆止阀存在反向泄漏，即从低压区泄漏到高压区，实际上这是不可能的。可见，得到的测试结果是无意义的。
发明内容
本申请提供了一种计算逆止阀泄漏率的方法及装置，目的在于解决现有的逆止阀的泄漏率的检测方法得到的结果无意义的问题。
为了实现上述目的，本申请提供了以下技术方案：
一种计算逆止阀泄漏率的方法，包括：
形成封闭空间，所述封闭空间的边界包括待测的逆止阀；
获取所述封闭空间的压力随时间的变化数据；
依据所述封闭空间的压力随时间的变化数据进行曲线拟合；
依据拟合出的曲线公式，得到因所述待测的逆止阀的泄漏引起的所述封闭空间的压力变化速率；
由所述封闭空间的压力变化速率得到所述待测的逆止阀的泄漏率。
可选地，所述依据拟合出的曲线公式，得到因所述待测的逆止阀的泄漏引起的所述封闭空间的压力变化速率包括：
依据拟合出的曲线公式P(t)＝-(P0-C)*e-λt+βt+C，得到因逆止阀的泄漏引起的封闭空间的压力变化速率β，其中，λ、β和C均为拟合所获取的常数，P(t)为所述封闭空间的压力值，t为时间。
可选地，除了所述待测的逆止阀之外，所述封闭空间的边界的其它部分无泄漏，但所述封闭空间与外界未达到热平衡，所述封闭空间与外界进行热传递。
可选地，所述获取封闭空间的压力随时间的变化数据包括：
通过压力表获取所述封闭空间每隔预设时间的压力值。
可选地，在所述依据所述封闭空间的压力随时间的变化数据进行曲线拟合之后，还包括：
计算曲线拟合优度，所述曲线拟合优度用于表征所述待测的逆止阀的泄漏率的准确性。
一种计算逆止阀泄漏率的装置，包括：
数据获取单元，用于获取封闭空间的压力随时间的变化数据，所述封闭空间的边界包括待测的逆止阀；
计算单元，用于依据所述封闭空间的压力随时间的变化数据进行曲线拟合，依据拟合出的曲线公式，得到因所述待测的逆止阀的泄漏引起的所述封闭空间的压力变化速率，并由所述封闭空间的压力变化速率得到所述待测的逆止阀的泄漏率。
可选地，所述计算单元用于依据拟合出的曲线公式，得到因所述待测的逆止阀的泄漏引起的所述封闭空间的压力变化速率包括：
所述计算单元具体用于，依据拟合出的曲线公式P(t)＝-(P0-C)*e-λt+βt+C，得到因逆止阀的泄漏引起的封闭空间的压力变化速率β，其中，λ、β和C均为拟合所获取的常数，P(t)为所述封闭空间的压力值，t为时间。
可选地，所述封闭空间的边界除所述待测的逆止阀之外，所述边界的其它部分无泄漏，但所述封闭空间与外界未达到热平衡，所述封闭空间与外界进行热传递。
可选地，所述数据获取单元用于获取封闭空间的压力随时间的变化数据包括：
所述数据获取单元具体用于，通过压力表获取所述封闭空间每隔预设时间的压力值。
可选地，所述计算单元还用于：
在所述依据所述封闭空间的压力随时间的变化数据进行曲线拟合之后，计算曲线拟合优度，所述曲线拟合优度用于表征所述待测的逆止阀的泄漏率的准确性。
本申请所述的方法及装置，利用边界包括待测的逆止阀的封闭空间的压力随时间的变化数据进行曲线拟合，并依据拟合出的曲线公式，得到因所述待测的逆止阀的泄漏引起的所述封闭空间的压力变化速率，由所述封闭空间的压力变化速率得到所述待测的逆止阀的泄漏率，因为使用曲线拟合这一新的方式得到逆止阀的泄漏率，所以，能够有效避免测试结果出现负值。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为测试系统简图；
图2为本申请实施例公开的一种计算逆止阀泄漏率的方法的流程图；
图3为本申请实施例公开的计算逆止阀泄漏率的GUI界面图；
图4为本申请实施例公开的一种计算逆止阀泄漏率的装置结构示意图。
具体实施方式
发明人在研究的过程中发现，测试出现负值，是由于测试条件不满足。测试要求在RCP及RRA温度稳定的状态下进行，但是机组上行阶段，测试前一回路和RRA系统的温度均为170℃左右，RRA系统隔离后冷却降温到60℃以下，而RCP121/221/321VP上游测试区域温度仍保持约170℃，6小时后进行 测试时RRA系统与一回路还没有达到温度平衡，不满足测试条件。此时进行测试，一、三环路会向RRA系统传递热量，导致封闭空间内温度降低，空间内液体收缩速率很大，超过了逆止阀的泄漏速率，从而导致封闭空间的压力不升反降。
分析历史测试数据发现，三个机组二环路RCP221VP泄漏率为负值的概率比一、三环路出现负值的概率要小，这也证明了这种现象是客观存在的。
实际上，封闭空间的压力变化由两方面因素引起，一部分是逆止阀的泄漏△QL，另一部分是热量转移导致的体积变化△QS。要想获得逆止阀准确的泄漏率△QL，就要排除热量传递引起的体积变化△QS对测试结果的干扰。这可以通过两种途径实现：
1)等△QS＝0后进行测试
等到封闭区域与周围达到热平衡(动态平衡)后，由热传递引起的收缩△QS就会消除。但是RRA隔离后，在缺少对流的情况下，测试区域向RRA的热传导很慢，达到热平衡需要很长时间，需要数小时甚至2至3天，对大修工期产生较大影响。
2)计算测试时的△QS
中压安注管线两道逆止阀之间的管线与余热排出系统(RRA)相连，
由于测试区域向余热排出系统(RRA)的传热过程十分复杂的动态过程，需要考虑温度、压力、传热系数、水密度变化、材料的弹性系数等很多因素，无法通过简单的计算获取某一时间由热量传递所引入的压力变化。
本申请的发明人创造性地提出了以下解决途径：
根据定期测试导则PTR的描述，逆止阀前后压差对泄漏的影响具有双重性，一方面压差增高，导致泄漏增大，另一方面，压差增高会使逆止阀阀瓣向密封面的推力更大，密封更严，泄漏率更小，所以测试允许在压差低于正常压差的情况下进行，测试结果对压差不太敏感，可以认为测试过程中逆止阀的泄漏引起的压力变化速率△QL是一个常数，这就简化了测试空间内的压力变化的影响因素，为通过寻求温度效应导致的容积收缩规律提供了便利。
根据傅立叶定律和牛顿冷却定律，发明人创造性地提出以下规律：
P(t)＝-(P0-C)*e-λt+βt+C
其中，P0为测试开始时封闭空间的初始压力，λ为系统综合传热系数，也是一个常数，β为逆止阀的泄漏引起的压力变化速率，C为积分常数。
推导过程如下：
傅立叶定律：在导热现象中，单位时间内通过给定截面的热量，正比例于垂直于该截面方向上的温度变化率和截面面积。
牛顿冷却定律：当物体表面与周围存在温度差时，单位时间从单位面积散失的热量与温度差成正比。
令此物之温度为T(t),在t＝t0时之温度为T0,又室温为H,则牛顿冷却定律为：
T′(t)＝-α(T(t)-H) 公式(1)
积分可得：
Figure BDA0001130158670000071
其中，α为热传递系数，是一个大于0的常数，它取决于物体的特性。
低温下，非饱和液态水的比体积与温度大致是呈线性关系，低温低压下，温度、压力变化不大的情况下，我们假定定压比热容是一个常数，测试区域管道的综合弹性系数也是一个常数，那么测试区域的压力变化就和测试区域的热量传递大致呈线性正比的关系。忽略低温水的可压缩性，根据以上两个定律，可以通过公式推导，得出封闭空间内由于热传导而产生的压力变化规律为：
P’(t)＝-δe-λt 公式(3)
其中，λ为系统综合传热系数，由系统本身的特性决定，是一个常数；δ是一个与测试开始时封闭空间的初始压力和λ决定的常数。
由于逆止阀的泄漏率对压差不敏感，可以假定其为常数β，那么封闭空间的压力变化可以表述为：
P(t)＝-(P0-C)*e-λt+βt+C 公式(4)
只要时间足够长，封闭区域内热传导进出达到动态平衡后，温度就会稳定下来，因为热传递而导致的收缩效应就会消除。体现到上面的公式(4)上就是，t足够大时，-(P0-C)*e-λt就会接近0，P(t)≈βt+C。
以上公式中的，λ、β、C都是一个常数，我们可以通过测试过程中采集到的多组压力变化数据，利用曲线拟合，拟合出一条有常数项的自然指数函数曲线，从而得到β。
依据上述原理，本申请中，依据采集到的压力随时间的变化数据，拟合出公式(4)，得到β，可以进一步得到待测的逆止阀的泄漏率。
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
本申请实施例公开的计算逆止阀泄漏率的方法如图2所示，包括以下步骤：
S201：形成封闭空间，所述封闭空间的边界包括待测的逆止阀，除了待测的逆止阀之外，封闭空间的边界的其它部分无泄漏，但所述封闭空间与外界未达到热平衡，所述封闭空间与外界进行热传递，即密闭性完好。
需要说明的是，封闭空间的一个示例如图1(加粗区域)所示。在实际中，管道的布置可能比图1所示要复杂，封闭空间也可以有出图1之外的其它形式。只要满足边界包括待测的逆止阀，除了待测的逆止阀之外，边界的其它部分的密闭性完好的封闭空间均在本申请的保护范围内。
本实施例中，封闭空间的边界的其它部分的密闭性完好的目的在于，避免边界中其它部分的泄漏对试验结果产生干扰。
S202：获取封闭空间的压力随时间的变化数据。
例如，N102大修时机组上行期间进行T-RIS041测试时进行了压力变化数据采集，每隔10min取一次数据，共采集10个试验数据，见表3。
表3
Figure BDA0001130158670000091
以图1为例，封闭空间(加粗区域)的压力值通过压力表080LP读取。
S203：依据封闭空间的压力随时间的变化数据进行曲线拟合，得到常数已知的公式(4)，即拟合出的曲线公式。
S204：依据拟合出的曲线公式，以排除所述封闭空间与外界的热量传递引起的封闭空间的水体积的变化因素，得到因逆止阀的泄漏引起的封闭空间的压力变化速率。
S205：由封闭空间的压力变化速率得到逆止阀的泄漏率。
本实施例中，S203～S205可以利用MATLAB实现：
利用MATLAB提供的GUI图形界面功能设计对话框，如图3所示。对话框左侧两列为数据输入框，第一个输入框输入试验开始时的初始压力，下 面依次输入每隔10min采集的试验数据，对话框中间的一列输出框输出测试结果(即逆止阀的泄漏率及拟合优度)，最下面为开始进行计算的“START”按钮；对话框右侧输出拟合曲线图，图中圆形点为数据点，曲线为拟合曲线，拟合曲线的公式显示在拟合曲线的右上角。
将S202的数据输入后，点击“START”按钮，进行自定义拟合，得到试验结果为3.626e-9bar/h，满足程序的验收准则。拟合曲线的函数公式为：
P(t)＝36.20*e-0.07t+0.00t+8.3 公式(5)
其中P(t)为试验区域的压力，得到的系统综合传热系数λ为0.07，代表泄漏率的常数β为0.00bar/10min(取小数点后两位)，逆止阀的泄漏率为0.00×6(换算为每小时)＝0bar/h，输出在对话框中间的结果输出框。
在得到逆止阀的泄漏率后，还可以执行以下步骤得到表征结果的准确性的参数。
S206：计算曲线拟合优度。
拟合优度(Goodness of Fit)是指拟合曲线对观测值的拟合程度，可由判定系数R2表征：
Figure BDA0001130158670000101
其中，TSS为Total Square Sum/总离差平方和。RSS为Regression Square Sum/回归平方和。ESS为Error Square Sum/误差平方和(残差平方和)
数值上，判定系数R2等于观测值与拟合值之间简单相关系数的平方，取值范围是[0，1]，R^2的值越接近1，说明回归直线对观测值的拟合程度越好。
结合图3，曲线拟合优度在对话框中间的一列的第二个输出框输出，为0.99794。
从以上步骤可以看出，本实施例所述方法，通过分析热量传递对压力变化的影响规律，建立数学模型，利用MATLAB软件对压力变化数据进行自定义拟合，通过拟合曲线排除热量传递对试验结果的影响，外推计算出逆止阀 的泄漏率。由于改进方案考虑了热量的传递过程对试验结果的影响因素并对其进行了工程量化，试验结果远比温度未稳定时获得的数据更接近真实值，能够有效避免测试结果出现负值。
需要说明的是，由于拟合曲线外推法是基于试验区域的水与外界存在热传递(总体为向外传递)，且未达到热平衡状态的分析，一旦测试现象与我们分析的不符，拟合曲线外推法就不再适用。现根据以上的分析，制定试验方法的选定及结果判断的准则如下：
1、如果所测试验区域的压力上涨速率为正值，说明试验管路已经达到热平衡，无需采用本方法，可采用现有方法直接判定试验结果。
也就是说，本申请实施例所述的方法，可以与现有方法结合使用：先使用现有方法测量逆止阀的泄漏率，如果结果为负值，则使用图2所示的方法测量逆止阀的泄漏率。
2、如果测量得到的封闭空间的压力数值(精确到小数点后一位)连续两次出现停止变化，说明系统达到热平衡状态，试验结果以最后稳定的数据为准。
3、如果拟合优度较差，则得到的结果误差较大，可信度下降，不建议采用。
本申请实施例还公开了一种计算逆止阀泄漏率的装置，如图4所示，包括：数据获取单元和计算单元。
数据获取单元用于获取封闭空间的压力随时间的变化数据，所述封闭空间的边界包括待测的逆止阀，除了待测的逆止阀之外，封闭空间的边界的其它部分的密闭性完好。
具体地，数据获取单元获取封闭空间的压力随时间的变化数据的具体实现方式为：通过压力表获取所述封闭空间每隔预设时间的压力值。
计算单元用于依据封闭空间的压力随时间的变化数据进行曲线拟合，依据拟合出的曲线公式，得到因逆止阀的泄漏引起的封闭空间的压力变化速率，再由封闭空间的压力变化速率得到逆止阀的泄漏率。
具体地，计算单元依据拟合出的曲线公式，得到因所述待测的逆止阀的泄漏引起的所述封闭空间的压力变化速率的具体实现方式为：依据拟合出的曲线公式P(t)＝-(P0-C)*e-λt+βt+C，得到因逆止阀的泄漏引起的封闭空间的压力变化速率β，其中，λ、β和C均为拟合所获取的常数。
进一步地，计算单元还可以用于计算曲线拟合优度，以反应结果的准确性。
本申请实施例方法所述的功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算设备可读取存储介质中。基于这样的理解，本申请实施例对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该软件产品存储在一个存储介质中，包括若干指令用以使得一台计算设备(可以是个人计算机，服务器，移动计算设备或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其它实施例的不同之处，各个实施例之间相同或相似部分互相参见即可。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下，在其它实施例中实现。因此，本申请将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (8)

1.一种计算逆止阀泄漏率的方法，其特征在于，包括：
形成封闭空间，所述封闭空间的边界包括待测的逆止阀；
获取所述封闭空间的压力随时间的变化数据；
依据所述封闭空间的压力随时间的变化数据进行曲线拟合；
依据拟合出的曲线公式，得到因所述待测的逆止阀的泄漏引起的所述封闭空间的压力变化速率；
由所述封闭空间的压力变化速率得到所述待测的逆止阀的泄漏率；
其中，所述依据拟合出的曲线公式，得到因所述待测的逆止阀的泄漏引起的所述封闭空间的压力变化速率包括：
依据拟合出的曲线公式P(t)＝-(P0-C)*e-λt+βt+C，得到因逆止阀的泄漏引起的封闭空间的压力变化速率β，其中，λ、β和C均为拟合所获取的常数，P(t)为t时刻所述封闭空间的压力，t为时间。
2.根据权利要求1所述的方法，其特征在于，除了所述待测的逆止阀之外，所述封闭空间的边界的其它部分无泄漏，但所述封闭空间与外界未达到热平衡，所述封闭空间与外界进行热传递。
3.根据权利要求1所述的方法，其特征在于，所述获取封闭空间的压力随时间的变化数据包括：
通过压力表获取所述封闭空间每隔预设时间的压力值。
4.根据权利要求1至3任一项所述的方法，其特征在于，在所述依据所述封闭空间的压力随时间的变化数据进行曲线拟合之后，还包括：
计算曲线拟合优度，所述曲线拟合优度用于表征所述待测的逆止阀的泄漏率的准确性。
5.一种计算逆止阀泄漏率的装置，其特征在于，包括：
数据获取单元，用于获取封闭空间的压力随时间的变化数据，所述封闭空间的边界包括待测的逆止阀；
计算单元，用于依据所述封闭空间的压力随时间的变化数据进行曲线拟合，依据拟合出的曲线公式，得到因所述待测的逆止阀的泄漏引起的所述封闭空间的压力变化速率，并由所述封闭空间的压力变化速率得到所述待测的逆止阀的泄漏率；
其中，所述计算单元用于依据拟合出的曲线公式，得到因所述待测的逆止阀的泄漏引起的所述封闭空间的压力变化速率包括：
所述计算单元具体用于，依据拟合出的曲线公式P(t)＝-(P0-C)*e-λt+βt+C，得到因逆止阀的泄漏引起的封闭空间的压力变化速率β，其中，λ、β和C均为拟合所获取的常数，P(t)为t时刻所述封闭空间的压力值，t为时间。
6.根据权利要求5所述的装置，其特征在于，所述封闭空间的边界除所述待测的逆止阀之外，所述边界的其它部分无泄漏，但所述封闭空间与外界未达到热平衡，所述封闭空间与外界进行热传递。
7.根据权利要求5所述的装置，其特征在于，所述数据获取单元用于获取封闭空间的压力随时间的变化数据包括：
所述数据获取单元具体用于，通过压力表获取所述封闭空间每隔预设时间的压力值。
8.根据权利要求5至7任一项所述的装置，其特征在于，所述计算单元还用于：
在所述依据所述封闭空间的压力随时间的变化数据进行曲线拟合之后，计算曲线拟合优度，所述曲线拟合优度用于表征所述待测的逆止阀的泄漏率的准确性。
CN201610892667.6A 2016-10-13 2016-10-13 一种计算逆止阀泄漏率的方法及装置 Active CN107941436B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610892667.6A CN107941436B (zh) 2016-10-13 2016-10-13 一种计算逆止阀泄漏率的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610892667.6A CN107941436B (zh) 2016-10-13 2016-10-13 一种计算逆止阀泄漏率的方法及装置

Publications (2)

Publication Number Publication Date
CN107941436A CN107941436A (zh) 2018-04-20
CN107941436B true CN107941436B (zh) 2020-01-31

Family

ID=61928918

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610892667.6A Active CN107941436B (zh) 2016-10-13 2016-10-13 一种计算逆止阀泄漏率的方法及装置

Country Status (1)

Country Link
CN (1) CN107941436B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109282953B (zh) * 2018-11-09 2023-10-27 成都珂睿科技有限公司 一种用于检测单向阀的内渗漏速率的装置及其测试方法
CN112098073B (zh) * 2020-08-28 2021-11-09 清华大学 逆止组合阀冷态逆止性能检测方法及其检测装置

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102735408A (zh) * 2012-02-28 2012-10-17 宝山钢铁股份有限公司 一种基于标准压力变化曲线的密闭容器密封性能检测方法
CN102759434A (zh) * 2011-04-29 2012-10-31 通用汽车环球科技运作有限责任公司 诊断车辆中的阀泄漏的系统和方法
CN102853261A (zh) * 2011-06-27 2013-01-02 国际商业机器公司 确定输送管道中的流体泄漏量的方法和装置
CN103900803A (zh) * 2012-12-26 2014-07-02 北京谊安医疗系统股份有限公司 一种机控安全阀检测方法和装置
DE102013220379A1 (de) * 2013-06-03 2014-12-04 Hyundai Motor Company System und Verfahren zum Steuern der Wasserstoffbefüllung eines Fahrzeug
CN104568344A (zh) * 2014-10-09 2015-04-29 芜湖扬宇机电技术开发有限公司 Egr管路故障检测诊断方法
CN104865026A (zh) * 2015-04-12 2015-08-26 中国计量学院 聚乙烯阀门密封性室温低压试验方法
JP2015206686A (ja) * 2014-04-21 2015-11-19 株式会社ケーヒン 検出装置及び検出方法
CN105675231A (zh) * 2016-02-22 2016-06-15 中航光电科技股份有限公司 一种负压密封检测装置
CN105715956A (zh) * 2014-12-22 2016-06-29 通用汽车环球科技运作有限责任公司 用于检查高压罐阀的方法和诊断装置、高压罐系统和具有该系统的汽车

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102759434A (zh) * 2011-04-29 2012-10-31 通用汽车环球科技运作有限责任公司 诊断车辆中的阀泄漏的系统和方法
CN102853261A (zh) * 2011-06-27 2013-01-02 国际商业机器公司 确定输送管道中的流体泄漏量的方法和装置
CN102735408A (zh) * 2012-02-28 2012-10-17 宝山钢铁股份有限公司 一种基于标准压力变化曲线的密闭容器密封性能检测方法
CN103900803A (zh) * 2012-12-26 2014-07-02 北京谊安医疗系统股份有限公司 一种机控安全阀检测方法和装置
DE102013220379A1 (de) * 2013-06-03 2014-12-04 Hyundai Motor Company System und Verfahren zum Steuern der Wasserstoffbefüllung eines Fahrzeug
JP2015206686A (ja) * 2014-04-21 2015-11-19 株式会社ケーヒン 検出装置及び検出方法
CN104568344A (zh) * 2014-10-09 2015-04-29 芜湖扬宇机电技术开发有限公司 Egr管路故障检测诊断方法
CN105715956A (zh) * 2014-12-22 2016-06-29 通用汽车环球科技运作有限责任公司 用于检查高压罐阀的方法和诊断装置、高压罐系统和具有该系统的汽车
CN104865026A (zh) * 2015-04-12 2015-08-26 中国计量学院 聚乙烯阀门密封性室温低压试验方法
CN105675231A (zh) * 2016-02-22 2016-06-15 中航光电科技股份有限公司 一种负压密封检测装置

Also Published As

Publication number Publication date
CN107941436A (zh) 2018-04-20

Similar Documents

Publication Publication Date Title
Di Marcello et al. The TRANSURANUS mechanical model for large strain analysis
Mao et al. Creep deformation and damage behavior of reactor pressure vessel under core meltdown scenario
CN107941436B (zh) 一种计算逆止阀泄漏率的方法及装置
Muhammad et al. LES study on the turbulent thermal stratification and thermo-mechanical fatigue analysis for NPP surge line
Zhang et al. A theoretical model for delayed hydride cracking velocity considering the temperature history and temperature gradients
Xu et al. ATHLET simulation of PKL I2. 2 IB-loca benchmark test and quantitative assessment
Garcia et al. Mono-dimensional mechanical modelling of fuel rods under normal and off-normal operating conditions
Manngård et al. Evaluation of loss-of-coolant accident simulation tests with the fuel rod analysis code FRAPTRAN-1.4
Gao et al. Development of PWR lower head failure model for severe accident analysis
Brunetto et al. Extension of the OFFBEAT fuel performance code to finite strains and validation against LOCA experiments
Jianfeng et al. Investigation on multilayer failure mechanism of RPV with a high temperature gradient from core meltdown scenario
Sagrado et al. Modeling RIA benchmark cases with FRAPTRAN and SCANAIR: A comparative exercise
US20210027901A1 (en) Methods and systems for facilitating the management of reactor transient conditions associated with reactors
Zheng et al. Damage performance analysis of prestressed concrete containments following a loss of coolant accident considering different external temperatures
Dragunov et al. Stress–strain kinetics in calculations of high-temperature strength and longevity of reactor structures
Gao et al. Study on the lower head failure in severe accident Part I: Model development and validation
Cancemi et al. Inverse Heat Conduction Problem in Estimating Nuclear Power Plant Piping Performance
Jeury et al. CABRI experimental reactor: experimental reassessment of CABRI core power and measurement uncertainties
Qibao et al. Fracture Assessment for Pressurized Water Reactor Vessel Nozzles Subjected to Pressure and Thermal Loading
Qian et al. Probabilistic Pressurized Thermal Shock Analysis for a Reactor Pressure Vessel Considering Plume Cooling Effect
Ruan et al. 3D CFD and FEM Evaluations of RPV Stress Intensity Factor During PTS Loading
Trewin et al. Best-Estimate Plus Uncertainty in Thermal-Hydraulic Mixing Calculations Using KWU-MIX for Pressurized Thermal Shock
Gilman et al. Environmentally-Assisted Fatigue Analysis Using a Strain-Based Approach
Syed et al. Identification of constitutive parameters for high temperature deformation of pressure tube of Indian PHWR considering multi-axial state of stress
Cadiou et al. Qualification of the MEXIICO loop dedicated to nuclear power transients: An experimental and modelling approach

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