CN115240881A - 一种安全壳泄漏率的测量方法 - Google Patents

一种安全壳泄漏率的测量方法 Download PDF

Info

Publication number
CN115240881A
CN115240881A CN202210508364.5A CN202210508364A CN115240881A CN 115240881 A CN115240881 A CN 115240881A CN 202210508364 A CN202210508364 A CN 202210508364A CN 115240881 A CN115240881 A CN 115240881A
Authority
CN
China
Prior art keywords
time
containment
leakage rate
denotes
temperature
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
CN202210508364.5A
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.)
China Nuclear Power Engineering Co Ltd
Harbin Engineering University
Original Assignee
China Nuclear Power Engineering Co Ltd
Harbin Engineering University
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 Nuclear Power Engineering Co Ltd, Harbin Engineering University filed Critical China Nuclear Power Engineering Co Ltd
Priority to CN202210508364.5A priority Critical patent/CN115240881A/zh
Publication of CN115240881A publication Critical patent/CN115240881A/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/002—Detection of leaks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • 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

一种安全壳泄漏率的测量方法
技术领域
本发明属于核工程技术领域，具体涉及一种安全壳泄漏率的 测量方法。
背景技术
安全壳是确保核电厂安全的最后一道物理屏障，在核电厂中， 安全壳整体试验，即安全壳打压试验，是对安全壳施工质量的检 验。安全壳整体试验过程的严谨性、试验结果的准确性、数据分 析的有效性是保证核电厂安全、稳定运行的重要指标。
安全壳密封性是衡量安全壳安装质量的重要指标，安全壳密 封性可以通过安全壳泄漏率测量试验确定。目前，国内外的安全 壳泄漏率测量试验均以绝对压力衰减法(又称为“压降法”)为 主，该方法至少存在以下不足：
(1)压降法测量安全壳泄漏率必须停止充气，之后监测安 全壳内的压力在试验时段内的变化值，由于在试验过程中的安全 壳内的气体会不断向外部泄漏，安全壳内的压力在理论上是逐步 下降的，为此，安全壳内的初始压力必须高于设计的压力平台， 以保证试验结束时安全壳内的压力不低于设计的试验压力。
(2)压降法的原理是计算两个时刻之间的安全壳内的干空 气的质量变化，如果试验过程中的湿度存在大幅、反复的变化， 采用干空气的泄漏率计算安全壳内真实气体的泄漏率会使试验 结果的偏差变大。
(3)压降法是采用多个时刻的干空气质量进行最小二乘法 求得质量变化斜率作为安全壳泄漏率，计算过程复杂，不易实现 安全壳泄漏率的实时显示。
发明内容
本发明要解决的技术问题是针对现有技术存在的以上不足， 提供一种安全壳泄漏率的测量方法，以实时的体积变化为依据， 可在恒压环境下完成，是一种全新的安全壳泄漏率测量技术，其 计算过程简单，可便于实现实时显示。
本发明要解决上述技术问题的技术方案是：
本发明提供一种安全壳泄漏率的测量方法，其包括：
向安全壳内充入气体，以使安全壳内的压力达到预设的试验 压力；
对安全壳内的实时压力进行监测，并根据安全壳内的实时压 力的监测结果向安全壳内补充气体，以维持安全壳内的实时压力 等于预设的试验压力，并检测向安全壳补充的气体的流量，以确 定测量体积泄漏率；
对安全壳内的温度和/或湿度进行监测，根据安全壳内的温 度和/或湿度监测结果，确定补偿泄漏率；
根据测量体积泄漏率和补偿泄漏率，确定安全壳泄漏率。
优选的是，所述对安全壳内的温度和/或湿度进行监测，根 据安全壳内的温度和/或湿度监测结果，确定补偿泄漏率，包括 以下步骤:
将安全壳的内部空间分割为k个虚拟的温度分区，对各个温 度分区的温度分别进行监测，得到各个温度分区在不同时刻的实 时温度数据，和/或，将安全壳的内部空间分割为m个虚拟的湿 度分区，对各个湿度分区的湿度分别进行监测，得到各个湿度分 区在不同时刻的实时湿度数据；
根据各个温度分区在不同时刻的实时温度数据，和/或，根 据各个湿度分区在不同时刻的实时湿度数据，得到所述补偿泄漏 率。
优选的是，所述补偿泄漏率具体是采用体积计算得到，则所 述安全壳泄漏率的计算步骤包括：
计算补偿体积泄漏率，以作为所述补偿泄漏率；
根据补偿体积泄漏率和测量体积泄漏率计算实际体积泄漏 率；
根据实际体积泄漏率计算质量泄漏率，以作为所述安全壳泄 漏率。
优选的是，所述根据补偿体积泄漏率和测量体积泄漏率计算 实际体积泄漏率，具体包括：
根据试验工况环境下/标准工况环境下的各个湿度分区在不 同时刻的实时湿度数据和各个温度分区在不同时刻的实时温度 数据，计算试验工况环境下/标准工况环境下的补偿体积泄漏率；
根据试验工况环境下或供气环境下的各个时刻向安全壳内 补充的气体的流量，计算试验工况环境下/标准工况环境下的测 量体积泄漏率；
将试验工况环境下/标准工况环境下的补偿体积泄漏率和试 验工况环境下/标准工况环境下的测量体积泄漏率相加，得到试 验工况环境下/标准工况环境下的实际体积泄漏率。
优选的是，在标准工况环境下的测量体积泄漏率L测,N∑∑的计 算式为：
Figure BDA0003636969330000031
其中，n表示时段数或周期数，i表示ti时刻或ti-1至ti时 间段，L测,P∑i表示试验工况环境下在ti时刻所有分区累计的测量 体积泄漏率,P0表示试验压力，PN表示标准工况环境下的压 力,TN表示标准工况环境下的温度，
Figure BDA0003636969330000032
表示ti-1至ti时段内的安 全壳内的有效比温度；
或者：
Figure BDA0003636969330000033
其中，n表示时段数或周期数，i表示ti时刻或ti-1至ti时 间段,L测,A∑i表示供气环境下在ti时刻所有分区累计的测量体积 泄漏率，PAi表示供气环境下向安全壳内充入气体的管路中的实 时压力，PN表示标准工况环境下的压力,TN表示标准工况环境下 的温度，
Figure BDA0003636969330000041
表示供气环境下在ti-1至ti时段内的安全壳内的有效 比温度；
在试验工况环境下的测量体积泄漏率L测,p∑∑的计算式为:
Figure BDA0003636969330000042
其中，L测,P∑i表示试验工况环境下在ti时刻所有分区累计的 测量体积泄漏率。
优选的是，所述补偿体积泄漏率是基于各个温度/湿度分区 分别进行补偿计算得到，其中，在试验工况环境下的补偿体积泄 漏率L补,P∑∑的计算式为：
Figure BDA0003636969330000043
其中，n表示时段数或周期数，m表示湿度分区的数量，i 表示ti时刻或ti-1至ti时间段，j表示第j温度分区或第j湿度 分区，k表示温度分区的数量，Hji表示第j湿度分区在ti时刻的 相对湿度，Hji-1表示第j湿度分区在ti-1时刻的相对湿度，PHji表示 第j湿度分区在ti时刻的饱和水蒸气分压，PHji-1表示第j湿度分 区在ti-1时刻的饱和水蒸气分压，VHj表示第j湿度分区占安全壳 的自由容积的百分比，V0表示安全壳的自由容积，P0表示试验压 力，Δt表示ti-1至ti时刻的时间长度，Tji表示第j温度分区在ti时刻的绝对温度，Tji-1表示第j温度分区在ti-1时刻的绝对温度， VTj表示第j温度分区占安全壳的自由容积的百分比；
在标准工况环境下的补偿体积泄漏率L补,N∑∑的计算式为：
Figure BDA0003636969330000051
其中，n表示时段数或周期数，m表示湿度分区的数量，i 表示ti时刻或ti-1至ti时间段，j表示第j温度分区或第j湿度 分区，k表示温度分区的数量，Hji表示第j湿度分区在ti时刻的 相对湿度，Hji-1表示第j湿度分区在ti-1时刻的相对湿度，PHji表示 第j湿度分区在ti时刻的饱和水蒸气分压，PHji-1表示第j湿度分 区在ti-1时刻的饱和水蒸气分压，VHj表示第j湿度分区占安全壳 的自由容积的百分比，V0表示安全壳的自由容积，P0表示试验压 力，PN表示标准工况环境下的压力，Δt表示ti-1至ti时刻的时间 长度,TN表示标准工况环境下的温度，THji-1表示第j湿度分区在 ti-1时刻的绝对温度，THji表示第j湿度分区在ti时刻的绝对温度， Tji表示第j温度分区在ti时刻的绝对温度，Tji-1表示第j温度分 区在ti-1时刻的绝对温度，VTj表示第j温度分区占安全壳的自由 容积的百分比。
优选的是，所述根据实际体积泄漏率计算质量泄漏率，具体 包括：
先计算一个Δt时间内的质量泄漏率N∑i，其计算式为：
Figure BDA0003636969330000052
其中，L实,N∑i表示标准工况环境下在ti时刻所有分区累计的 实际体积泄漏率，m气表示空气摩尔质量，m水表示水蒸气的摩 尔质量，
Figure BDA0003636969330000053
表示ti时刻的平均水蒸气分压，
Figure BDA0003636969330000054
表示ti-1时刻的 平均水蒸气分压，
Figure BDA0003636969330000055
表示ti时刻的平均水蒸气分压，R表示理想 气体常数，P0表示试验压力，PN表示标准工况环境下的压力，TN表 示标准工况环境下的温度；
或者：
Figure BDA0003636969330000061
其中，L实,P∑i表示试验工况环境下在ti时刻所有分区累计的 实际体积泄漏率，m气表示空气摩尔质量，m水表示水蒸气的摩 尔质量，
Figure BDA0003636969330000062
表示ti时刻的平均水蒸气分压，
Figure BDA0003636969330000063
表示ti-1时刻的 平均水蒸气分压，R表示理想气体常数，P0表示试验压力，PN表 示标准工况环境下的压力，
Figure BDA0003636969330000064
表示ti至ti-1时段内的安全壳内的 有效比温度；
再计算连续多个Δt时间内总的质量泄漏率M∑∑，总的 质量泄漏率M∑∑即为所述安全壳泄漏率，其计算式为：
Figure BDA0003636969330000065
其中，n表示时段数或周期数，i表示ti时刻或ti-1至ti时 间段。
优选的是，所述补偿体积泄漏率是基于所述安全壳内各个温 度分区的平均温度和各个湿度分区的平均湿度进行补偿计算得 到，则安全壳泄漏率的采用以下计算式计算得到，其中，在试验 工况环境下，所述安全壳泄漏率的计算式包括：
L实,P∑i＝L测,P∑i+L补,P∑i
Figure BDA0003636969330000066
Figure BDA0003636969330000067
其中，L实,P∑i表示试验工况环境下在ti时刻所有分区累计的 实际体积泄漏率，L测,P∑i表示试验工况环境下在ti时刻所有分区 累计的测量体积泄漏率，L补,P∑i表示试验工况环境下在ti时刻所 有分区累计的补偿体积泄漏率,下标p表示试验工况条件,下标i 表示ti时刻，
Figure BDA0003636969330000071
表示在ti时刻所有分区的平均相对湿度，
Figure BDA0003636969330000072
表 示在ti时刻所有分区的平均饱和水蒸气分压，
Figure BDA0003636969330000073
表示在ti-1时 刻所有分区的平均相对湿度，
Figure BDA0003636969330000074
表示在ti-1时刻所有分区的平 均饱和水蒸气分压，
Figure BDA0003636969330000075
表示ti时刻的平均温度，
Figure BDA0003636969330000076
表示ti-1时 刻的平均温度，V0表示安全壳的自由容积，Δt表示ti-1至ti时刻 的时间长度，MP∑i表示试验工况环境下在ti时刻所有分区累计 的质量泄漏率,m气表示空气摩尔质量，m水表示水蒸气的摩尔质 量,P0表示试验压力,
Figure BDA0003636969330000077
表示ti时刻的平均水蒸气分压,R表示理 想气体常数；
在标准工况环境下，所述安全壳泄漏率的计算式包括：
L实,N∑i＝L测,N∑i+L补,N∑i
Figure BDA0003636969330000078
Figure BDA0003636969330000079
其中，L实N∑i表示标准工况环境下在ti时刻所有分区累计的 实际体积泄漏率，L测,N∑i表示标准工况环境下在ti时刻所有分区 累计的测量体积泄漏率，L补,N∑i表示标准工况环境下在ti时刻所 有分区累计的补偿体积泄漏率，
Figure BDA00036369693300000710
表示在ti时刻所有分区的平均 相对湿度，
Figure BDA00036369693300000711
表示在ti时刻所有分区的平均水蒸气分压，
Figure BDA00036369693300000712
表 示在ti-1时刻所有分区的平均相对湿度，
Figure BDA00036369693300000713
表示在ti-1时刻所 有分区的平均水蒸气分压，
Figure BDA00036369693300000714
表示ti时刻的平均温度，
Figure BDA00036369693300000715
表示 ti-1时刻的平均温度，V0表示安全壳的自由容积，Δt表示ti-1至 ti时刻的时间长度,P0表示试验压力,PN表示标准工况环境下的 压力,TN表示标准工况环境下的温度,MN∑i表示标准工况环境 下在ti时刻所有分区累计的质量泄漏率，m气表示空气摩尔质量， m水表示水蒸气的摩尔质量,
Figure BDA00036369693300000716
表示ti时刻的平均水蒸气分 压,R表示理想气体常数。
优选的是，所述对安全壳内的温度和/或湿度进行监测，根 据安全壳内的温度和/或湿度监测结果，确定补偿泄漏率，包括 以下步骤:
将安全壳内划分为a个分块，分别监测每个分块的温度和湿 度；
根据对各个分块温度和湿度的监测结果，以及各个分块的体 积，计算得到所述补偿泄漏率。
优选的是，所述安全壳泄漏率的计算步骤包括：
计算a个分块在ti至ti+1时间段的体积泄漏率；
根据a个分块在ti至ti+1时间段的体积泄漏率，计算a个分 块在ti至ti+1时间段的质量泄漏率，以作为所述安全壳泄漏率。
优选的是，所述a个分块在ti至ti+1时间段的体积泄漏率
Figure BDA0003636969330000081
的计算式为：
Figure BDA0003636969330000082
其中，a表示分块的数量,Lin,i+1表示第ti+1时刻充入气体的 管路的出口处的充气体积流量，Tc,i+1,j表示第j分块在ti+1时刻 的安全壳内气体的绝对温度,Tc,i,j表示第j分块在ti时刻的安全 壳内气体的绝对温度，mc,i,j表示第j分块在ti时刻的安全壳内气 体的质量，Rg,eq,i,j表示第j分块在ti时刻的安全壳内气体的折合 气体常数，Pc表示充入气体的管路的出口处/安全壳内气体的压 力，Δt表示ti时刻至ti+1时刻的时间长度，Vc,i,j表示第j分块所 对应的体积,Hc,i+1,j表示第j分块在ti+1时刻的安全壳内相对湿 度，Hc,i,j表示第j分块在ti时刻的安全壳内相对湿度，f(Tc,i+1,j)表 示第j分块在ti+1时刻的饱和水蒸气分压,f(Tc,i,j)表示第j分块 在ti+1时刻的安全壳内饱和水蒸气分压；
所述a个分块在ti至ti+1时间段的质量泄漏率Gout,i+1j的计算 式为：
Figure BDA0003636969330000091
其中，j表示第j分块，i表示第ti时刻或ti至ti+1时间段,a 表示分块的数量,Lin,i+1表示第ti+1时刻充入气体的管路的出口 处的充气体积流量，Tc,i+1,j表示第j分块在ti+1时刻的安全壳内 气体的绝对温度,Tc,i,j表示第j分块在ti时刻的安全壳内气体的 绝对温度，mc,i,j表示第j分块在ti时刻的安全壳内气体的质量， Rg,eq,i,j表示第j分块在ti时刻的安全壳内气体的折合气体常数， Pc表示充入气体的管路的出口处/安全壳内气体的压力，Δt表示 ti至ti+1时刻的时间长度，Vc,i,j表示第j分块所对应的体积, Hc,i+1,j表示第j分块在ti+1时刻的安全壳内相对湿度，Hc,i,j表示第 j分块在ti时刻的安全壳内相对湿度，f(Tc,i+1,j)表示第j分块在 ti+1时刻的安全壳内饱和水蒸气分压,f(Tc,i,j)表示第j分块在ti+1时刻的安全壳内饱和水蒸气分压,Rg,eq,i+1,j表示第j分块在ti+1时刻的安全壳内气体的折合气体常数。
本发明提供的安全壳泄漏率的测量方法，可在维持安全壳内 部压力恒定的条件下完成，是一种恒压法安全壳泄漏率测量技术， 在测量过程中要保持连续充气，以保持安全壳内压力恒定，其原 理是通过分析实时的气体体积变化，建立恒压法计算模型，从而 确定安全壳泄漏率，与传统的压降法的原理(以分析两个时刻之 间安全壳内的干空气的质量变化)完全不同。并且，本方法充分 地关注到了测量过程中的温度变化和湿度变化对安全壳泄漏率 测量的影响，并结合温度、湿度连续变化的特征，可给定合适的 采集周期计算温度补偿和湿度补偿，也就是说，本方法中安全壳 泄漏率的计算模型中分析的是安全壳内真实泄漏的气体，本方法 的计算模型考虑了湿度变化等因素对压力的影响，进行了这部分影响的补偿，最终测量结果不会受到安全壳内的湿度变化等因素 的影响，相比于压降法，准确性更高。
附图说明
图1为本发明实施例中安全壳泄漏率的测量方法的示意图；
图2为本发明实施例中安全壳泄漏率的测量方法的原理图；
图3为本发明实施例中安全壳泄漏率的测量方法的计算流 程示意图；
图4为本实施例中24h内的安全壳内温度的实时变化图；
图5为本实施例中24h内的安全壳内湿度的实时变化图；
图6为本实施例中24h内的安全壳内压力的实时变化图；
图7为本实施例中24h内的向安全壳补充的气体的流量的实 时变化图；
图8为本实施例中24h内的实际体积泄漏率L实,N∑i的实时变 化图；
图9为本实施例中24h内的质量泄漏率M∑i的实时变化图。
具体实施方式
下面将结合本发明中的附图，对发明中的技术方案进行清楚、 完整的描述，显然，所描述的实施例是本发明的一部分实施例， 而不是全部的实施例。基于本发明中的实施例，本领域普通技术 人员在没有做出创造性劳动的前提下所获得的所有其他实施例， 都属于本发明的范围。
实施例1
如图1所示，本实施例公开一种安全壳泄漏率的测量方法， 包括：
向安全壳内充入气体，以使安全壳内的压力达到预设的试验 压力；
对安全壳内的实时压力进行监测，并根据安全壳内的实时压 力的监测结果向安全壳内补充气体，以维持安全壳内的实时压力 等于预设的试验压力，并检测向安全壳补充的气体的流量，以确 定测量体积泄漏率；
对安全壳内的温度和/或湿度进行监测，根据安全壳内的温 度和/或湿度监测结果，确定补偿泄漏率；
根据测量体积泄漏率和补偿泄漏率，确定安全壳泄漏率。
如图2所示，本方法的原理为：将安全壳分为多个分区，假 设各个分区稳定在一个初始状态：存在一个对空的泄漏口，该泄 漏口的泄漏率为L0，同时，存在一个恒压为P0的无限大空间， 不断的向各个分区内补充压缩气体，以保证各个分区内的压力始 终为P0，补充的压缩气体的流量为LD；设各个分区内的自由容 积为V0，试验压力为P0，各个分区的初始温度为T0，各个分区 的初始湿度为H0，各个分区内的气体分子数为N0，在上述各项 参数稳定的状态下，则测量泄漏率与实际泄漏率保持一致，即 LD＝L0；并且，由于在各个分区为恒压的试验环境下，虽然各个 分区的压力和容积保持不变，但湿度H和温度T可能随时间发 生变化，因此，本方法还考虑了湿度和温度对测量结果的影响， 将充压管路上测量的LD作为安全壳泄漏率的基准值或基础值， 并在此基础上引入了安全壳内(简称壳内)的温度变化体积补偿 和湿度变化体积补偿，得到补偿泄漏率。其中，计算流程如图3 所示。
具体来说，在安全壳上设置充压管路和恒压测试管路，通过 充压管路向安全壳内充入气体(比如，空气、压缩空气)，以对 安全壳进行正常充压，并对安全壳内压力或充压管路内的压力进 行检测并进行反馈，直至安全壳内压力接近预设的压力平台时， 切换至恒压测试管路进行充入气体，以直至安全壳内的压力达到 预设的试验压力，或者，直接由充压管路通入气体至安全壳内压 力达到预设的试验压力，之后，开始安全壳泄漏率测量。在测量 过程中，由于安全壳不断向外部泄漏，在理论上，安全内的压力 是不断下降的，通过恒压测试管路向安全壳内补充气体，同时， 通过对安全壳内的实时压力进行监测，并根据安全壳内实时压力 监测信号反馈结果控制补充气体的流量，使安全壳内的实时压力 始终维持等于预设的试验压力，通过读取补充气体的流量，得到 所述测量体积泄漏率。
需要说明的是，维持安全壳内的实时压力等于预设的试验压 力的控制方式，还可以是：采用一个压力大于试验压力P0的气 源，通过安全壳壳内压力反馈实时调节向安全壳补充的气体的流 量，不断向安全壳内补充气体，从而维持安全壳内压力始终为 P0。与前文所述的控制方式相比，该控制方式可有效缩短首次达 到预设试验压力的时间以及在试验平台上压力微幅波动的周期。
对安全壳内的温度和/或湿度进行监测，根据安全壳内的温 度和/或湿度监测结果，确定补偿泄漏率，具体可以包括以下步 骤:
将安全壳的内部空间分割为k个虚拟的温度分区，对各个温 度分区的温度分别进行监测，得到各个温度分区在不同时刻的实 时温度数据，和/或，将安全壳的内部空间分割为m个虚拟的湿 度分区，对各个湿度分区的湿度分别进行监测，得到各个湿度分 区在不同时刻的实时湿度数据；根据各个温度分区在不同时刻的 实时温度数据，和/或，根据各个湿度分区在不同时刻的实时湿 度数据，得到所述补偿泄漏率。
具体来说，通过布置在安全壳内的温度传感器，将安全壳内 部的自由空间分割为若干个(如k个)虚拟的温度分区，其中， 将典型区域、气体流动不通畅的房间作为一个分区，每个温度分 区的体积系数不可超过0.1，通过上述布置的温度传感器对各个 温度分区的温度进行监测，得到各个温度分区的温度变化数据， 根据各个温度分区的温度变化数据，在k个温度分区内进行温度 补偿(即温度修正)计算，即计算温度对体积变化的补偿，其中， 第j分区的初始温度设为Tj0，在△t＝ti-ti-1时间段内，温度一直在 缓慢的变化，ti时刻温度设为Tji；通过布置在安全壳内的湿度传 感器，将安全壳内部的自由空间分割为若干个(如m个)虚拟 的湿度分区，其中，将典型区域、气体流动不通畅的房间作为一 个分区，通过上述布置的湿度传感器对各个湿度分区的湿度进行 监测，得到各个湿度分区的湿度变化数据，根据各个湿度分区的 湿度变化数据，在m个分区内进行湿度补偿(即湿度修正)计算，即计算湿度对体积变化的补偿，其中，第j分区的初始湿度 设为Hj0，在△t＝ti-ti-1时间段内，湿度一直在缓慢的变化，ti时 刻湿度设为Hji。
本实施例中，补偿泄漏率具体是采用体积计算得到，安全壳 泄漏率的计算步骤包括：计算补偿体积泄漏率，以作为所述补偿 泄漏率；根据补偿体积泄漏率和测量体积泄漏率计算实际体积泄 漏率；根据实际体积泄漏率计算质量泄漏率，以作为所述安全壳 泄漏率。
具体来说，根据补偿体积泄漏率和测量体积泄漏率计算实际 体积泄漏率，具体包括：根据试验工况环境下/标准工况环境下 的各个湿度分区在不同时刻的实时湿度数据和各个温度分区在 不同时刻的实时温度数据，计算试验工况环境下/标准工况环境 下的补偿体积泄漏率；根据试验工况环境下或供气环境下的各个 时刻向安全壳内补充的气体的流量，计算试验工况环境下/标准 工况环境下的测量体积泄漏率；将试验工况环境下/标准工况环 境下的补偿体积泄漏率和试验工况环境下/标准工况环境下的测 量体积泄漏率相加，得到试验工况环境下/标准工况环境下的实 际体积泄漏率。
以在m个湿度分区和k个温度分区内同时进行的泄漏率补 偿情况为例，下面提供几组计算模型，具体如下：
(1)第一组计算模型
(1-1)计算补偿泄漏率
补偿体积泄漏率基于各个温度/湿度分区分别进行补偿计算 得到，其中：
在试验工况环境下的补偿体积泄漏率L补,P∑∑的计算式为：
Figure BDA0003636969330000141
其中，n表示时段数或周期数，m表示湿度分区的数量，i 表示ti时刻或ti-1至ti时间段，j表示第j温度分区或第j湿度 分区，k表示温度分区的数量，Hji表示第j湿度分区在ti时刻的 相对湿度，Hji-1表示第j湿度分区在ti-1时刻的相对湿度，PHji表示 第j湿度分区在ti时刻的饱和水蒸气分压，PHji-1表示第j湿度分 区在ti-1时刻的饱和水蒸气分压，VHj表示第j湿度分区占安全壳 的自由容积的百分比，V0表示安全壳的自由容积，P0表示试验压 力(0.42MPa.g，计算中取绝对压力0.5213MPa.a，下文中不再 一一赘述)，Δt表示ti-1至ti时刻的时间长度，Tji表示第j温度 分区在ti时刻的绝对温度，Tji-1表示第j温度分区在ti-1时刻的 绝对温度，VTj表示第j温度分区占安全壳的自由容积的百分比；
在标准工况环境下的补偿体积泄漏率L补,N∑∑的计算式为：
Figure BDA0003636969330000142
其中，n表示时段数或周期数，m表示湿度分区的数量，i 表示ti时刻或ti-1至ti时间段，j表示第j温度分区或第j湿度 分区，k表示温度分区的数量，Hji表示第j湿度分区在ti时刻的 相对湿度，Hji-1表示第j湿度分区在ti-1时刻的相对湿度，PHji表示 第j湿度分区在ti时刻的饱和水蒸气分压，PHji-1表示第j湿度分 区在ti-1时刻的饱和水蒸气分压，VHj表示第j湿度分区占安全壳 的自由容积的百分比，V0表示安全壳的自由容积，P0表示试验压 力，PN表示标准工况环境下的压力(1个大气压，0.1013MPa.a， 下文中不再一一赘述)，Δt表示ti-1至ti时刻的时间长度,TN表 示标准工况环境下的温度(0℃，绝对温度为273.15K，下文中 不再一一赘述)，THji-1表示第j湿度分区在ti-1时刻的绝对温度， THji表示第j湿度分区在ti时刻的绝对温度，Tji表示第j温度分区 在ti时刻的绝对温度，Tji-1表示第j温度分区在ti-1时刻的绝对 温度，VTj表示第j温度分区占安全壳的自由容积的百分比。
(1-2)计算测量体积泄漏率
测量体积泄漏率一般采用具有温度、压力测量功能的流量计 直接累计测量得到，但对于不具有上述累计功能的流量计，则按 以下计算式计算得到，其中：
在标准工况环境下的测量体积泄漏率L测,N∑∑的计算式为：
Figure BDA0003636969330000151
该计算式是由安全壳内为试验工况环境转换为标准工况环 境得到，其中，n表示时段数或周期数，i表示ti时刻或ti-1至 ti时间段，L测,P∑i表示试验工况环境下在ti时刻所有分区累计的 测量体积泄漏率(即试验工况环境下在ti时刻检测到的向安全 壳补充的气体的流量值),P0表示试验压力，PN表示标准工况环 境下的压力,TN表示标准工况环境下的温度，
Figure BDA0003636969330000152
表示ti-1至ti时 段内的安全壳内的有效比温度，
Figure BDA0003636969330000153
或者，
Figure BDA0003636969330000154
该计算式是由充压管路的供气环境转换为标准工况环境得 到，其中，n表示时段数或周期数，i表示ti时刻或ti-1至ti时 间段,L测,A∑i表示供气环境下在ti时刻所有分区累计的测量体积 泄漏率(即供气环境下在ti时刻检测到的向安全壳补充的气体 的流量值)，PAi表示供气环境下向安全壳内充入气体的管路处 的实时压力，PN表示标准工况环境下的压力,TN表示标准工况环 境下的温度，
Figure BDA0003636969330000161
表示供气环境下在ti-1至ti时段内的安全壳内的 有效比温度。
在试验工况环境下的测量体积泄漏率L测,p∑∑的计算式为:
Figure BDA0003636969330000162
其中，L测,P∑i表示试验工况环境下在ti时刻所有分区累计的 测量体积泄漏率(即试验工况环境下在ti时刻检测到的向安全 壳补充的气体的流量值)。
(1-3)计算实际体积泄漏率
实际体积泄漏率L实等于测量体积泄漏率与补偿体积泄漏率 之和，则在试验工况环境下的实际体积泄漏率L实,P∑∑的计算式为：
L实,P∑∑＝L测,P∑∑+L补,P∑∑
在标准工况环境下的实际体积泄漏率L实,N∑∑的计算式为：
L实,N∑∑＝L测,N∑∑+L补,N∑∑
其中，需要说明的是：在本实施例中，下标N表示标准工 况环境，下标P表示试验工况环境，下标∑∑表所有时刻所有分 区累计，下文中不再一一赘述。
(1-4)计算质量泄漏率
根据实际体积泄漏率计算质量泄漏率，质量泄漏率的计算步 骤包括：
先计算一个Δt时间内的质量泄漏率M∑i，其计算式为：
Figure BDA0003636969330000163
其中，下标∑i表示ti时刻所有分区累计，下标N∑i表示标准 工况环境下ti时刻所有分区累计，L实,N∑i表示标准工况环境下在 ti时刻所有分区累计的实际体积泄漏率，m气表示空气摩尔质量， m水表示水蒸气的摩尔质量，
Figure BDA0003636969330000164
表示ti时刻的平均水蒸气分压，
Figure BDA0003636969330000171
表示ti-1时刻的平均水蒸气分压，
Figure BDA0003636969330000172
表示ti时刻的平均水蒸 气分压，R表示理想气体常数(8.314J·mol-1·k-1)，P0表示试 验压力，PN表示标准工况环境下的压力，TN表示标准工况环境下 的温度；
或者，
Figure BDA0003636969330000173
其中，下标∑i表示ti时刻所有分区累计，下标P∑i表示试验 工况环境下ti时刻所有分区累计，L实,P∑i表示试验工况环境下在 ti时刻所有分区累计的实际体积泄漏率，m气表示空气摩尔质量， m水表示水蒸气的摩尔质量，
Figure BDA0003636969330000174
表示ti时刻的平均水蒸气分压，
Figure BDA0003636969330000175
表示ti-1时刻的平均水蒸气分压，R表示理想气体常数，P0表 示试验压力，PN表示标准工况环境下的压力，
Figure BDA0003636969330000176
表示ti至ti-1时段 内的安全壳内的有效比温度，
Figure BDA0003636969330000177
再计算连续多个Δt时间内总的质量泄漏率M∑∑，总的质量泄 漏率即本方法的安全壳泄漏率，其计算式为：
Figure BDA0003636969330000178
其中，n表示时段数或周期数，i表示ti时刻或ti-1至ti时间 段。
(2)第二组计算模型
与(1)第一组计算模型相比，对计算过程进行了适当的简 化，补偿体积泄漏率是基于安全壳内各个温度分区的平均温度和 各个湿度分区的平均湿度进行补偿计算得到，则安全壳泄漏率的 采用以下计算式计算得到：
在试验工况环境下，安全壳泄漏率的计算式包括：
L实,P∑i＝L测,P∑i+L补,P∑i
Figure BDA0003636969330000179
Figure BDA0003636969330000181
其中，L实,P∑i表示试验工况环境下在ti时刻所有分区累计的 实际体积泄漏率，L测,P∑i表示试验工况环境下在ti时刻所有分区 累计的测量体积泄漏率(即试验工况环境下在ti时刻检测到的 向安全壳补充的气体的流量值)，L补,P∑i表示试验工况环境下在 ti时刻所有分区累计的补偿体积泄漏率,下标p表示试验工况条 件,下标i表示ti时刻，
Figure BDA0003636969330000182
表示在ti时刻所有分区的平均相对湿 度，
Figure BDA0003636969330000183
表示在ti时刻所有分区的平均饱和水蒸气分压，
Figure BDA0003636969330000184
表示 在ti-1时刻所有分区的平均相对湿度，
Figure BDA0003636969330000185
表示在ti-1时刻所有 分区的平均饱和水蒸气分压，
Figure BDA0003636969330000186
表示ti时刻的平均温度，
Figure BDA0003636969330000187
表 示ti-1时刻的平均温度，V0表示安全壳的自由容积，Δt表示ti-1至 ti时刻的时间长度，MP∑i表示试验工况环境下在ti时刻所有分区累计的质量泄漏率,m气表示空气摩尔质量，m水表示水蒸气的摩 尔质量,P0表示试验压力,
Figure BDA0003636969330000188
表示ti时刻的平均水蒸气分压,R表 示理想气体常数；
在标准工况环境下，安全壳泄漏率的计算式包括：
L实,N∑i＝L测,N∑i+L补,N∑i
Figure BDA0003636969330000189
Figure BDA00036369693300001810
其中，L实N∑i表示标准工况环境下在ti时刻所有分区累计的 实际体积泄漏率，L测,N∑i表示标准工况环境下在ti时刻所有分区 累计的测量体积泄漏率，L补,N∑i表示标准工况环境下在ti时刻所 有分区累计的补偿体积泄漏率，
Figure BDA00036369693300001811
表示在ti时刻所有分区的平均 相对湿度，
Figure BDA00036369693300001812
表示在ti时刻所有分区的平均饱和水蒸气分压，
Figure BDA00036369693300001813
表示在ti-1时刻所有分区的平均相对湿度，
Figure BDA00036369693300001814
表示在ti-1时刻所有分区的平均饱和水蒸气分压，
Figure BDA00036369693300001815
表示ti时刻的平均温度，
Figure BDA00036369693300001816
表示ti-1时刻的平均温度，V0表示安全壳的自由容积，Δt表 示ti-1至ti时刻的时间长度,P0表示试验压力,PN表示标准工况 环境下的压力,TN表示标准工况环境下的温度,MN∑i表示标准 工况环境下在ti时刻所有分区累计的质量泄漏率，m气表示空气 摩尔质量，m水表示水蒸气的摩尔质量,
Figure BDA0003636969330000191
表示ti时刻的平均水 蒸气分压,R表示理想气体常数。
需要说明的是，安全壳内的实际的平均压力都是围绕P0以微 小的幅度波动。因此，本实实施例中各计算式中的P0还可以替换 成安全壳内的实时的测量压力Pi(替换前后相对变化仅万分之五 左右)，并且，用Pi计算会与安全壳内压力的实际的变化过程更 一致，结果更精确一些。
下面以1000m3安全壳为例，其中，试验压力P0(绝压)为 533.65kPa、温度分区为28个，各个温度传感器的布置情况如表 1所示，湿度分区为10个，各个湿度传感器的布置情况如表2 所示，对上述第一组计算模型计算过程进行详述，具体如下：
表1
Figure BDA0003636969330000192
Figure BDA0003636969330000201
表2
Figure BDA0003636969330000202
由于本系统测量过程中安全壳内的温度、湿度、压力是实时 变化的，进气管线上补充的气体的流量也会发生相应的变化，通 过安全壳内压力计、温度传感器、湿度传感器等检测器连续采集24h内的数据来完成安全壳泄漏率的测量与分析，其中：
24h内的安全壳内温度的实时变化如图4所示；
24h内的安全壳内湿度的实时变化如图5所示；
24h内的安全壳内压力的实时变化如图6所示；
24h内的向安全壳补充的气体的流量的实时变化如图7所示。
以24h内的第一个时间段Δt＝10s(即0到10s)为例，说 明单时刻体积泄漏率L实,N∑i和单时刻质量泄漏率M∑i的计算：
T1-T28温度传感器在ti-1时刻采集得到的温度值数据(单 位℃)依次为：
32.5680,32.5840,32.6770,32.7190,32.8540,33.0000,32. 9560,32.9010,32.8570,33.0530,32.8720,32.8330,32.7520,33 .1960,32.8320,33.1960,33.0260,33.3470,32.9680,33.2100,3 3.3170,33.1310,33.2670,33.3580,33.2160,33.0400,32.3840, 32.4840；
H1-H10湿度传感器在ti-1时刻采集得到的湿度值数据(％) 依次为：
48.1898,46.8971,47.1527,47.7582,46.3522,48.7407,46. 4725,46.7261,46.1951,48.9654；
在ti-1时刻采集得到的安全壳内的实时压力为：533.652kpa；
在ti-1时刻采集得到向安全壳内充入气体的管路中气体的流 量值，即测量泄漏率(标准工况)为：0.7945m3/h；
T1-T28温度传感器在ti时刻(ti-ti-1＝10s)采集得到的温度 值数据(单位℃)依次为：
32.5670,32.5810,32.6790,32.7170,32.8600,32.9950,32. 9550,32.9060,32.8540,33.0530,32.8750,32.8310,32.7530,33 .1950,32.8340,33.1980,33.0270,33.3480,32.9680,33.2100,3 3.3160,33.1320,33.2700,33.3560,33.2150,33.0390,32.3850, 32.4870；
H1-H10湿度传感器在ti时刻采集得到的湿度值数据(％)依 次为：
48.1968,46.8991,47.1556,47.7602,46.3552,48.7447,46. 4715,46.7281,46.1990,48.9595；
在ti时刻采集得到的安全壳内的实时压力为：533.652kpa；
在ti时刻采集得到向安全壳内充入气体的管路中的气体的 流量值，即测量泄漏率(标准工况)为：0.7975m3/h；
将上述采集得到的数据，代入到上文所述的计算式中，得到：
ti时刻安全壳内有效比温度为：
Figure BDA0003636969330000221
在ti-1至ti时间段内的测量体积泄漏率为：
L测,N∑i＝0.796m3/h。
按照实施例1中的计算模型，分别计算：
在ti-1至ti时间段内的补偿泄漏率为：
L补,N∑i＝3.16703m3/h。
在ti-1至ti时间段内的实际体积泄漏率为：
L实,N∑i＝L测,N∑i+L补,N∑i＝3.96268m3/h
通过ti-1时刻安全壳内温度、湿度，可计算得到ti-1时刻的 安全壳内水蒸气分压为：
Figure BDA0003636969330000222
通过ti时刻安全壳内温度、湿度，可计算得到ti时刻的安 全壳内水蒸气分压为：
Figure BDA0003636969330000223
在ti-1至ti时间段内的实际质量泄漏率为：
Figure BDA0003636969330000224
将采集的24h内的数据划分为24×360个10s的时间段，每 个时间段均按照上述单时刻计算过程，得到24×360组单时刻的 测量体积泄漏率L测,N∑i、补偿体积泄漏率L补,N∑i、实际体积泄漏 率L实,N∑i和质量泄漏率M∑i。这里不在一一赘述其他单时刻时段 内的计算，仅以图例说明24h内的实际体积泄漏率L实,N∑i和质量 泄漏率M∑i变化过程，其中：
24h内的实际体积泄漏率L实,N∑i的实时变化如图8所示。
24h内的质量泄漏率M∑i的实时变化如图9所示。
对24h内所有时段的单时刻的测量体积泄漏率L测,N∑i累加计 算平均值，计算得到24h内的测量体积泄漏率为：
L测,N∑∑＝1.05534m3/h
对24h内所有时段的单时刻补偿体积泄漏率L补,N∑i累加计算平 均值，计算得到24h内的补偿体积泄漏率为：
L测,N∑∑＝0.15381m3/h
安全壳体容积1000m3，相对体积泄漏率为0.617553％安全 壳内气体总体积/天。
24h内的平均的实际体积泄漏率为：
L实,N∑∑＝L测,N∑∑+L补,N∑∑＝1.20915m3/h
对24h内所有时段的单时刻的质量泄漏率M∑i，累加计算平 均值，计算得到24h的质量泄漏率为：
Figure BDA0003636969330000231
在0时刻，安全壳内气体总质量6063.65Kg，相对质量泄漏 率为0.618286％安全壳内气体总质量/天。
本实施例的安全壳泄漏率的测量方法，可在维持安全壳内部 压力恒定的条件下完成，是一种恒压法安全壳泄漏率测量技术， 在测量过程中要保持连续充气，以保持安全壳内压力恒定，其原 理是通过分析实时的体积变化，建立恒压法计算模型，从而确定 安全壳泄漏率，与传统的压降法的原理(以分析两个时刻之间安 全壳内的干空气的质量变化)完全不同。并且，本方法充分地关 注到了测量过程中的温度变化和湿度变化对安全壳泄漏率测量 的影响，并结合温度、湿度连续变化的特征，可给定合适的采集 周期计算温度补偿和湿度补偿，也就是说，本方法中安全壳泄漏 率的计算模型中分析的是安全壳内真实泄漏的气体，本方法的计 算模型考虑了湿度变化等因素对压力的影响，进行了这部分影响 的补偿，最终测量结果不会受到安全壳内的湿度变化等因素的影 响，相比于压降法，准确性更高。
实施例2
本实施例公开一种安全壳泄漏率的测量方法，相比于实施例 1所述的安全壳泄漏率的测量方法，区别在于：
对安全壳内的温度和/或湿度进行监测，根据安全壳内的温 度和/或湿度监测结果，确定补偿泄漏率，包括以下步骤:将安全 壳内划分为a个分块，分别监测每个分块的温度和湿度；根据对 各个分块温度和湿度的监测结果，以及各个分块的体积，计算得 到补偿泄漏率。
本方法采用是与实施例1中的第一组计算模型和第二组计 算模型都不同的(3)第三组计算模型，其具体计算步骤包括：
计算a个分块在ti至ti+1时间段的体积泄漏率；
根据a个分块在ti至ti+1时间段的体积泄漏率，计算a个分 块在ti至ti+1时间段的质量泄漏率，质量泄漏率即本方法测量得 到的安全壳泄漏率。
具体来说，a个分块在ti至ti+1时间段的体积泄漏率
Figure BDA0003636969330000241
的 计算式为：
Figure BDA0003636969330000242
其中，a表示分块的数量,Lin,i+1表示第ti+1时刻充入气体的 管路(即本文中恒压测试管路)的出口处的充气体积流量，Tc,i+1,j表示第j分块在ti+1时刻的安全壳内气体的绝对温度,Tc,i,j表示第 j分块在ti时刻的安全壳内气体的绝对温度，mc,i,j表示第j分块 在ti时刻的安全壳内气体(混合气体)的质量，Rg,eq,i,j表示第j 分块在ti时刻的安全壳内气体(混合气体)的折合气体常数，Pc表示充入气体的管路的出口处/安全壳内气体的压力，Δt表示ti时刻至ti+1时刻的时间长度，Vc,i,j表示第j分块所对应的体积, Hc,i+1,j表示第j分块在ti+1时刻的安全壳内相对湿度，Hc,i,j表示第 j分块在ti时刻的安全壳内相对湿度，f(Tc,i+1,j)表示第j分块在 ti+1时刻的安全壳内饱和水蒸气分压,f(Tc,i,j)表示第j分块在ti+1时刻的安全壳内饱和水蒸气分压；
a个分块在ti至ti+1时间段的质量泄漏率Gout,i+1j的计算式为：
Figure BDA0003636969330000251
其中，j表示第j分块，i表示第ti时刻或ti至ti+1时间段,a 表示分块的数量,Lin,i+1表示第ti+1时刻充入气体的管路的出口 处的充气体积流量，Tc,i+1,j表示第j分块在ti+1时刻的安全壳内 气体的绝对温度,Tc,i,j表示第j分块在ti时刻的安全壳内气体的 绝对温度，mc,i,j表示第j分块在ti时刻的安全壳内气体(混合气 体)的质量，Rg,eq,i,j表示第j分块在ti时刻的安全壳内气体(混 合气体)的折合气体常数，Pc表示充入气体的管路的出口处/安全壳内气体的压力，Δt表示ti时刻至ti+1时刻的时间长度，Vc,i,j表 示第j分块所对应的体积,Hc,i+1,j表示第j分块在ti+1时刻的安 全壳内相对湿度，Hc,i,j表示第j分块在ti时刻的安全壳内相对湿 度，f(Tc,i+1,j)表示第j分块在ti+1时刻的安全壳内饱和水蒸气分 压,f(Tc,i,j)表示第j分块在ti+1时刻的安全壳内饱和水蒸气分压, Rg,eq,i+1,j表示第j分块在ti+1时刻的安全壳内气体的折合气体常 数。
本实施例的安全壳泄漏率的测量方法，可在维持安全壳内部 压力恒定的条件下完成，是一种恒压法安全壳泄漏率测量技术， 在测量过程中要保持连续充气，以保持安全壳内压力恒定，其原 理是通过分析实时的体积变化，建立恒压法计算模型，从而确定 安全壳泄漏率，与传统的压降法的原理(以分析两个时刻之间安 全壳内的干空气的质量变化)完全不同。并且，本方法充分地关 注到了测量过程中的温度变化和湿度变化对安全壳泄漏率测量 的影响，并结合温度、湿度连续变化的特征，可给定合适的采集 周期计算温度补偿和湿度补偿，也就是说，本方法中安全壳泄漏 率的计算模型中分析的是安全壳内真实泄漏的气体，本方法的计 算模型考虑了湿度变化等因素对压力的影响，进行了这部分影响 的补偿，最终测量结果不会受到安全壳内的湿度变化等因素的影 响，相比于压降法，准确性更高。
可以理解的是，以上实施方式仅仅是为了说明本发明的原理 而采用的示例性实施方式，然而本发明并不局限于此。对于本领 域内的普通技术人员而言，在不脱离本发明的精神和实质的情况 下，可以做出各种变型和改进，这些变型和改进也视为本发明的 保护范围。

Claims (11)

1.一种安全壳泄漏率的测量方法，包括：
向安全壳内充入气体，以使安全壳内的压力达到预设的试验压力；
对安全壳内的实时压力进行监测，并根据安全壳内的实时压力的监测结果向安全壳内补充气体，以维持安全壳内的实时压力等于预设的试验压力，并检测向安全壳补充的气体的流量，以确定测量体积泄漏率；
对安全壳内的温度和/或湿度进行监测，根据安全壳内的温度和/或湿度监测结果，确定补偿泄漏率；
根据测量体积泄漏率和补偿泄漏率，确定安全壳泄漏率。
2.根据权利要求1所述的安全壳泄漏率的测量方法，其特征在于，所述对安全壳内的温度和/或湿度进行监测，根据安全壳内的温度和/或湿度监测结果，确定补偿泄漏率，包括以下步骤:
将安全壳的内部空间分割为k个虚拟的温度分区，对各个温度分区的温度分别进行监测，得到各个温度分区在不同时刻的实时温度数据，和/或，
将安全壳的内部空间分割为m个虚拟的湿度分区，对各个湿度分区的湿度分别进行监测，得到各个湿度分区在不同时刻的实时湿度数据；
根据各个温度分区在不同时刻的实时温度数据，和/或，根据各个湿度分区在不同时刻的实时湿度数据，得到所述补偿泄漏率。
3.根据权利要求2所述的安全壳泄漏率的测量方法，其特征在于，所述补偿泄漏率具体是采用体积计算得到，则所述安全壳泄漏率的计算步骤包括：
计算补偿体积泄漏率，以作为所述补偿泄漏率；
根据补偿体积泄漏率和测量体积泄漏率计算实际体积泄漏率；
根据实际体积泄漏率计算质量泄漏率，以作为所述安全壳泄漏率。
4.根据权利要求3所述的安全壳泄漏率的测量方法，其特征在于，所述根据补偿体积泄漏率和测量体积泄漏率计算实际体积泄漏率，具体包括：
根据试验工况环境下/标准工况环境下的各个湿度分区在不同时刻的实时湿度数据和各个温度分区在不同时刻的实时温度数据，计算试验工况环境下/标准工况环境下的补偿体积泄漏率；
根据试验工况环境下或供气环境下的各个时刻向安全壳内补充的气体的流量，计算试验工况环境下/标准工况环境下的测量体积泄漏率；
将试验工况环境下/标准工况环境下的补偿体积泄漏率和试验工况环境下/标准工况环境下的测量体积泄漏率相加，得到试验工况环境下/标准工况环境下的实际体积泄漏率。
5.根据权利要求4所述的安全壳泄漏率的测量方法，其特征在于，在标准工况环境下的测量体积泄漏率L测,N∑∑的计算式为：
Figure RE-FDA0003865889630000021
其中，n表示时段数或周期数，i表示ti时刻或ti-1至ti时间段，L测,P∑i表示试验工况环境下在ti时刻所有分区累计的测量体积泄漏率,P0表示试验压力，PN表示标准工况环境下的压力,TN表示标准工况环境下的温度，
Figure RE-FDA0003865889630000022
表示ti-1至ti时段内的安全壳内的有效比温度，
或者，
Figure RE-FDA0003865889630000031
其中，n表示时段数或周期数，i表示ti时刻或ti-1至ti时间段,L测,A∑i表示供气环境下在ti时刻所有分区累计的测量体积泄漏率，PAi表示供气环境下向安全壳内充入气体的管路中的实时压力，PN表示标准工况环境下的压力,TN表示标准工况环境下的温度，
Figure RE-FDA0003865889630000032
表示供气环境下在ti-1至ti时段内的安全壳内的有效比温度；
在试验工况环境下的测量体积泄漏率L测,p∑∑的计算式为:
Figure RE-FDA0003865889630000033
其中，L测,P∑i表示试验工况环境下在ti时刻所有分区累计的测量体积泄漏率。
6.根据权利要求4所述的安全壳泄漏率的测量方法，其特征在于，所述补偿体积泄漏率是基于各个温度/湿度分区分别进行补偿计算得到，其中，
在试验工况环境下的补偿体积泄漏率L补,P∑∑的计算式为：
Figure RE-FDA0003865889630000034
其中，n表示时段数或周期数，m表示湿度分区的数量，i表示ti时刻或ti-1至ti时间段，j表示第j温度分区或第j湿度分区，k表示温度分区的数量，Hji表示第j湿度分区在ti时刻的相对湿度，Hji-1表示第j湿度分区在ti-1时刻的相对湿度，PHji表示第j湿度分区在ti时刻的饱和水蒸气分压，PHji-1表示第j湿度分区在ti-1时刻的饱和水蒸气分压，VHj表示第j湿度分区占安全壳的自由容积的百分比，V0表示安全壳的自由容积，P0表示试验压力，Δt表示ti-1至ti时刻的时间长度，Tji表示第j温度分区在ti时刻的绝对温度，Tji-1表示第j温度分区在ti-1时刻的绝对温度，VTj表示第j温度分区占安全壳的自由容积的百分比；
在标准工况环境下的补偿体积泄漏率L补,N∑∑的计算式为：
Figure RE-FDA0003865889630000041
其中，n表示时段数或周期数，m表示湿度分区的数量，i表示ti时刻或ti-1至ti时间段，j表示第j温度分区或第j湿度分区，k表示温度分区的数量，Hji表示第j湿度分区在ti时刻的相对湿度，Hji-1表示第j湿度分区在ti-1时刻的相对湿度，PHji表示第j湿度分区在ti时刻的饱和水蒸气分压，PHji-1表示第j湿度分区在ti-1时刻的饱和水蒸气分压，VHj表示第j湿度分区占安全壳的自由容积的百分比，V0表示安全壳的自由容积，P0表示试验压力，PN表示标准工况环境下的压力，Δt表示ti-1至ti时刻的时间长度,TN表示标准工况环境下的温度，THji-1表示第j湿度分区在ti-1时刻的绝对温度，THji表示第j湿度分区在ti时刻的绝对温度，Tji表示第j温度分区在ti时刻的绝对温度，Tji-1表示第j温度分区在ti-1时刻的绝对温度，VTj表示第j温度分区占安全壳的自由容积的百分比。
7.根据权利要求4所述的安全壳泄漏率的测量方法，其特征在于，所述根据实际体积泄漏率计算质量泄漏率，具体包括：
先计算一个Δt时间内的质量泄漏率M∑i，其计算式为：
Figure RE-FDA0003865889630000051
其中，L实,N∑i表示标准工况环境下在ti时刻所有分区累计的实际体积泄漏率，m气表示空气摩尔质量，m水表示水蒸气的摩尔质量，
Figure RE-FDA0003865889630000052
表示ti时刻的平均水蒸气分压，
Figure RE-FDA0003865889630000053
表示ti-1时刻的平均水蒸气分压，
Figure RE-FDA0003865889630000054
表示ti时刻的平均水蒸气分压，R表示理想气体常数，P0表示试验压力，PN表示标准工况环境下的压力，TN表示标准工况环境下的温度，
或者，
Figure RE-FDA0003865889630000055
其中，L实,P∑i表示试验工况环境下在ti时刻所有分区累计的实际体积泄漏率，m气表示空气摩尔质量，m水表示水蒸气的摩尔质量，
Figure RE-FDA0003865889630000056
表示ti时刻的平均水蒸气分压，
Figure RE-FDA0003865889630000057
表示ti-1时刻的平均水蒸气分压，R表示理想气体常数，P0表示试验压力，PN表示标准工况环境下的压力，
Figure RE-FDA0003865889630000058
表示ti至ti-1时段内的安全壳内的有效比温度；
再计算连续多个Δt时间内总的质量泄漏率M∑∑，总的质量泄漏率M∑∑即为所述安全壳泄漏率，其计算式为：
Figure RE-FDA0003865889630000059
其中，n表示时段数或周期数，i表示ti时刻或ti-1至ti时间段。
8.根据权利要求4所述的安全壳泄漏率的测量方法，其特征在于，所述补偿体积泄漏率是基于所述安全壳内各个温度分区的平均温度和各个湿度分区的平均湿度进行补偿计算得到，则安全壳泄漏率的采用以下计算式计算得到，
在试验工况环境下，所述安全壳泄漏率的计算式包括：
L实,P∑i＝L测,P∑i+L补,P∑i
Figure RE-FDA0003865889630000061
Figure RE-FDA0003865889630000062
其中，L实,P∑i表示试验工况环境下在ti时刻所有分区累计的实际体积泄漏率，L测,P∑i表示试验工况环境下在ti时刻所有分区累计的测量体积泄漏率，L补,P∑i表示试验工况环境下在ti时刻所有分区累计的补偿体积泄漏率,下标p表示试验工况条件,下标i表示ti时刻，
Figure RE-FDA0003865889630000063
表示在ti时刻所有分区的平均相对湿度，
Figure RE-FDA0003865889630000064
表示在ti时刻所有分区的平均饱和水蒸气分压，
Figure RE-FDA0003865889630000065
表示在ti-1时刻所有分区的平均相对湿度，
Figure RE-FDA0003865889630000066
表示在ti-1时刻所有分区的平均饱和水蒸气分压，
Figure RE-FDA0003865889630000067
表示ti时刻的平均温度，
Figure RE-FDA0003865889630000068
表示ti-1时刻的平均温度，V0表示安全壳的自由容积，Δt表示ti-1至ti时刻的时间长度，MP∑i表示试验工况环境下在ti时刻所有分区累计的质量泄漏率,m气表示空气摩尔质量，m水表示水蒸气的摩尔质量,P0表示试验压力,
Figure RE-FDA0003865889630000069
表示ti时刻的平均水蒸气分压,R表示理想气体常数；
在标准工况环境下，所述安全壳泄漏率的计算式包括：
L实,N∑i＝L测,N∑i+L补,N∑i
Figure RE-FDA00038658896300000610
Figure RE-FDA00038658896300000611
其中，L实N∑i表示标准工况环境下在ti时刻所有分区累计的实际体积泄漏率，L测,N∑i表示标准工况环境下在ti时刻所有分区累计的测量体积泄漏率，L补,N∑i表示标准工况环境下在ti时刻所有分区累计的补偿体积泄漏率，
Figure RE-FDA00038658896300000612
表示在ti时刻所有分区的平均相对湿度，
Figure RE-FDA0003865889630000071
表示在ti时刻所有分区的平均水蒸气分压，
Figure RE-FDA0003865889630000072
表示在ti-1时刻所有分区的平均相对湿度，
Figure RE-FDA0003865889630000073
表示在ti-1时刻所有分区的平均水蒸气分压，
Figure RE-FDA0003865889630000074
表示ti时刻的平均温度，
Figure RE-FDA0003865889630000075
表示ti-1时刻的平均温度，V0表示安全壳的自由容积，Δt表示ti-1至ti时刻的时间长度,P0表示试验压力,PN表示标准工况环境下的压力,TN表示标准工况环境下的温度,MN∑i表示标准工况环境下在ti时刻所有分区累计的质量泄漏率，m气表示空气摩尔质量，m水表示水蒸气的摩尔质量,
Figure RE-FDA0003865889630000076
表示ti时刻的平均水蒸气分压,R表示理想气体常数。
9.根据权利要求1所述的安全壳泄漏率的测量方法，其特征在于，所述对安全壳内的温度和/或湿度进行监测，根据安全壳内的温度和/或湿度监测结果，确定补偿泄漏率，包括以下步骤:
将安全壳划分为a个分块，分别监测每个分块的温度和湿度；
根据对各个分块温度和湿度的监测结果，以及各个分块的体积，计算得到所述补偿泄漏率。
10.根据权利要求9所述的安全壳泄漏率的测量方法，其特征在于，所述安全壳泄漏率的计算步骤包括：
计算a个分块在ti至ti+1时间段的体积泄漏率；
根据a个分块在ti至ti+1时间段的体积泄漏率计算a个分块在ti至ti+1时间段的质量泄漏率，以作为所述安全壳泄漏率。
11.根据权利要求10所述的安全壳泄漏率的测量方法，其特征在于，所述a个分块在ti至ti+1时间段的体积泄漏率
Figure RE-FDA0003865889630000077
的计算式为：
Figure RE-FDA0003865889630000081
其中，a表示分块的数量,Lin,i+1表示第ti+1时刻充入气体的管路的出口处的充气体积流量，Tc,i+1,j表示第j分块在ti+1时刻的安全壳内气体的绝对温度,Tc,i,j表示第j分块在ti时刻的安全壳内气体的绝对温度，mc,i,j表示第j分块在ti时刻的安全壳内气体的质量，Rg,eq,i,j表示第j分块在ti时刻的安全壳内气体的折合气体常数，Pc表示充入气体的管路的出口处/安全壳内气体的压力，Δt表示ti时刻至ti+1时刻的时间长度，Vc,i,j表示第j分块所对应的体积,Hc,i+1,j表示第j分块在ti+1时刻的安全壳内相对湿度，Hc,i,j表示第j分块在ti时刻的安全壳内相对湿度，f(Tc,i+1,j)表示第j分块在ti+1时刻的饱和水蒸气分压,f(Tc,i,j)表示第j分块在ti+1时刻的安全壳内饱和水蒸气分压；
所述a个分块在ti至ti+1时间段的质量泄漏率Gout,i+1j的计算式为：
Figure RE-FDA0003865889630000082
其中，j表示第j分块，i表示第ti时刻或ti至ti+1时间段,a表示分块的数量,Lin,i+1表示第ti+1时刻充入气体的管路的出口处的充气体积流量，Tc,i+1,j表示第j分块在ti+1时刻的安全壳内气体的绝对温度,Tc,i,j表示第j分块在ti时刻的安全壳内气体的绝对温度，mc,i,j表示第j分块在ti时刻的安全壳内气体的质量，Rg,eq,i,j表示第j分块在ti时刻的安全壳内气体的折合气体常数，Pc表示充入气体的管路的出口处/安全壳内气体的压力，Δt表示ti至ti+1时刻的时间长度，Vc,i,j表示第j分块所对应的体积,Hc,i+1,j表示第j分块在ti+1时刻的安全壳内相对湿度，Hc,i,j表示第j分块在ti时刻的安全壳内相对湿度，f(Tc,i+1,j)表示第j分块在ti+1时刻的安全壳内饱和水蒸气分压,f(Tc,i,j)表示第j分块在ti+1时刻的安全壳内饱和水蒸气分压,Rg,eq,i+1,j表示第j分块在ti+1时刻的安全壳内气体的折合气体常数。
CN202210508364.5A 2022-05-10 2022-05-10 一种安全壳泄漏率的测量方法 Pending CN115240881A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210508364.5A CN115240881A (zh) 2022-05-10 2022-05-10 一种安全壳泄漏率的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210508364.5A CN115240881A (zh) 2022-05-10 2022-05-10 一种安全壳泄漏率的测量方法

Publications (1)

Publication Number Publication Date
CN115240881A true CN115240881A (zh) 2022-10-25

Family

ID=83668111

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210508364.5A Pending CN115240881A (zh) 2022-05-10 2022-05-10 一种安全壳泄漏率的测量方法

Country Status (1)

Country Link
CN (1) CN115240881A (zh)

Similar Documents

Publication Publication Date Title
US5485754A (en) Apparatus and method for measuring the air flow component and water vapor component of air/water vapor streams flowing under vacuum
ES2771798T3 (es) Procedimiento y dispositivo de medición para la determinación de variables específicas para la constitución del gas
ES2638793T3 (es) Métodos y sistemas para el ensayo de fugas
US4941345A (en) Method and apparatus for the measurement of gas properties
CN115240881A (zh) 一种安全壳泄漏率的测量方法
US20230168171A1 (en) High throughput high resolution gas sorption screening
CN105203267B (zh) 一种空间挤出常压累积检漏系统与方法
CN115938623A (zh) 安全壳泄漏率测量的试验方法、系统、电子设备及介质
CN203534928U (zh) 测量吸附工质对吸附/解吸性能的实验装置
CN106289372A (zh) 焊接绝热气瓶日蒸发率测量仪测试方法
KR102661248B1 (ko) 증기 펌프 성능을 평가하기 위한 방법
US4916630A (en) Bed tester for molecular sieve oxygen concentrator
CN113933213A (zh) 基于气体替代法的二元混合气体混气比测量方法及装置
CN114864115A (zh) 一种安全壳泄漏率测量试验系统
van der Beek et al. Gas oil piston prover, primary reference values for gas-volume
RU2298774C1 (ru) Способ контроля герметичности емкостей
Belzile et al. PVT study of trifluoromethane by the Burnett Method
WO2020136476A1 (en) Method for determining and monitoring the gas consumption in a gas network under pressure or under vacuum and gas network
RU79998U1 (ru) Установка для поверки промышленных счетчиков газа
RU138370U1 (ru) Устройство для испытания изделий при быстром изменении давления
EA029297B1 (ru) Устройство для проверки герметичности полимерных топливопроводов транспортных средств
RU2246102C1 (ru) Способ контроля герметичности заправленной теплоносителем гидравлической системы терморегулирования космического аппарата, снабженной гидропневматическим компенсатором
RU2527659C1 (ru) Способ контроля герметичности
RU2710006C2 (ru) Способ испытания ёмкости на герметичность
Inoue The thermal transpiration measurement of He, H2, O2, N2, CO, and CH4

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