CN115240881A - Method for measuring leakage rate of containment - Google Patents

Method for measuring leakage rate of containment 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
Chinese (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 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/en
Publication of CN115240881A publication Critical patent/CN115240881A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G21NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21CNUCLEAR REACTORS
    • G21C17/00Monitoring; Testing ; Maintaining
    • G21C17/002Detection of leaks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Abstract

The invention discloses a method for measuring the leakage rate of a containment vessel, which comprises the following steps: filling gas into the containment vessel to enable the pressure in the containment vessel to reach the preset test pressure; monitoring the real-time pressure in the containment, supplementing gas into the containment according to the monitoring result of the real-time pressure in the containment to maintain the real-time pressure in the containment equal to a preset test pressure, and detecting the flow of the gas supplemented to the containment to determine the measured volume leakage rate; monitoring the temperature and/or humidity in the containment, and determining a compensation leakage rate according to the monitoring result of the temperature and/or humidity in the containment; and determining the leakage rate of the containment according to the measured volume leakage rate and the compensation leakage rate. The invention is a brand new containment leakage rate measuring technology which can be completed in a constant pressure environment on the basis of real-time volume change, has simple calculation process and can conveniently realize real-time display.

Description

Method for measuring leakage rate of containment vessel
Technical Field
The invention belongs to the technical field of nuclear engineering, and particularly relates to a method for measuring a leakage rate of a containment.
Background
The containment is the last physical barrier for ensuring the safety of the nuclear power plant, and in the nuclear power plant, the containment integral test, namely the containment pressurization test, is the test of the construction quality of the containment. The rigidness of the whole test process of the containment, the accuracy of the test result and the effectiveness of data analysis are important indexes for ensuring the safe and stable operation of the nuclear power plant.
Containment tightness is an important index for measuring the installation quality of the containment, and can be determined through a containment leakage rate measurement test. At present, the leakage rate measurement tests of the safety shell at home and abroad mainly adopt an absolute pressure attenuation method (also called a pressure drop method), and the method at least has the following defects:
(1) The pressure drop method is used for measuring the leakage rate of the containment vessel, the inflation must be stopped, then the change value of the pressure in the safety shell in the test period is monitored, and the pressure in the containment vessel is reduced gradually in theory because the gas in the safety shell can leak to the outside continuously in the test process, therefore, the initial pressure in the containment vessel must be higher than a designed pressure platform, and the pressure in the containment vessel is guaranteed not to be lower than the designed test pressure at the end of the test.
(2) The principle of the pressure drop method is to calculate the mass change of dry air in a containment vessel between two moments, and if the humidity in the test process has large and repeated changes, the leakage rate of the real gas in the containment vessel is calculated by adopting the leakage rate of the dry air, so that the deviation of the test result is increased.
(3) The pressure drop method adopts dry air quality at multiple moments to obtain the quality change slope as the containment leakage rate by a least square method, the calculation process is complex, and the real-time display of the containment leakage rate is not easy to realize.
Disclosure of Invention
The invention aims to solve the technical problem of the prior art, and provides a containment leakage rate measuring method, which can be completed in a constant pressure environment based on real-time volume change, is a brand new containment leakage rate measuring technology, has a simple calculation process, and can be used for conveniently realizing real-time display.
The technical scheme for solving the technical problems is as follows:
the invention provides a method for measuring a leakage rate of a containment vessel, which comprises the following steps:
filling gas into the containment vessel to enable the pressure in the containment vessel to reach a preset test pressure;
monitoring the real-time pressure in the containment, supplementing gas into the containment according to the monitoring result of the real-time pressure in the containment to maintain the real-time pressure in the containment equal to a preset test pressure, and detecting the flow of the gas supplemented to the containment to determine the measured volume leakage rate;
monitoring the temperature and/or the humidity in the containment, and determining a compensation leakage rate according to the temperature and/or humidity monitoring result in the containment;
and determining the leakage rate of the containment according to the measured volume leakage rate and the compensation leakage rate.
Preferably, the monitoring the temperature and/or humidity in the containment vessel and determining the compensated leakage rate according to the monitoring result of the temperature and/or humidity in the containment vessel comprise the following steps:
dividing the internal space of the containment into k virtual temperature partitions, monitoring the temperature of each temperature partition respectively to obtain real-time temperature data of each temperature partition at different moments, and/or dividing the internal space of the containment into m virtual humidity partitions, monitoring the humidity of each humidity partition respectively to obtain real-time humidity data of each humidity partition at different moments;
and obtaining the compensation leakage rate according to the real-time temperature data of each temperature partition at different moments and/or according to the real-time humidity data of each humidity partition at different moments.
Preferably, the compensated leakage rate is obtained by volume calculation, and the calculation step of the containment leakage rate includes:
calculating a compensated volumetric leak rate as the compensated leak rate;
calculating an actual volume leakage rate according to the compensated volume leakage rate and the measured volume leakage rate;
and calculating the mass leakage rate according to the actual volume leakage rate to be used as the containment leakage rate.
Preferably, the calculating of the actual volumetric leakage rate according to the compensated volumetric leakage rate and the measured volumetric leakage rate specifically comprises:
calculating the compensation volume leakage rate under the test working condition environment/the standard working condition environment according to the real-time humidity data of each humidity partition under the test working condition environment/the standard working condition environment at different moments and the real-time temperature data of each temperature partition under the test working condition environment/the standard working condition environment at different moments;
calculating the measured volume leakage rate under the test working condition environment/the standard working condition environment according to the flow of the gas supplemented into the containment vessel at each moment under the test working condition environment or the gas supply environment;
and adding the compensation volume leakage rate under the test working condition environment/standard working condition environment and the measurement volume leakage rate under the test working condition environment/standard working condition environment to obtain the actual volume leakage rate under the test working condition environment/standard working condition environment.
Preferably, the volume leakage rate L is measured under a standard working condition environment Measuring, N sigma The calculation formula of (A) is as follows:
Figure BDA0003636969330000031
wherein n represents the number of time segments or cycles, and i represents t i Time of day or t i-1 To t i Time period, L Measure, P ∑ i Showing the test condition at t i Cumulative measured volumetric leakage rate, P, for all zones at any time 0 Denotes the test pressure, P N Indicating the pressure, T, in the environment of the standard operating conditions N The temperature in the environment of the standard operating condition is shown,
Figure BDA0003636969330000032
represents t i-1 To t i Effective specific temperature within the containment over a period of time;
or:
Figure BDA0003636969330000033
wherein n represents the number of time segments or cycles, and i represents t i Time of day or t i-1 To t i Time period, L Measuring, A sigma i Denotes the gas supply environment at t i Cumulative measured volume leakage rate, P, for all zones at that time Ai Representing the real-time pressure in the pipeline filling the containment with gas in the environment of gas supply, P N Indicating the pressure, T, in the environment of standard conditions N Which represents the temperature in the environment of the standard operating conditions,
Figure BDA0003636969330000041
indicating at t in the environment of the supply of air i-1 To t i Effective specific temperature in the containment over a period of time;
volume leakage rate L measured under test working condition environment Measuring, p sigma The calculation formula of (A) is as follows:
Figure BDA0003636969330000042
wherein L is Measuring, P sigma i Showing the test condition at t i The cumulative measured volumetric leak rate for all zones at that time.
Preferably, the compensation volume leakage rate is obtained by respectively performing compensation calculation based on each temperature/humidity partition, wherein the compensation volume leakage rate L under the test working condition environment Complement, P sigma The calculation formula of (c) is:
Figure BDA0003636969330000043
wherein n represents the number of time periods or cycles, m represents the number of humidity zones, i represents t i Time of day or t i-1 To t i Time period, j denotes the jth temperature or humidity compartment, k denotes the number of temperature compartments, H ji Denotes the jth humidity zone at t i Relative humidity at the moment H ji-1 Denotes the jth humidity zone at t i-1 Relative humidity at the moment P Hji Denotes the jth humidity division at t i Saturated partial pressure of water vapor at time, P Hji-1 Denotes the j moisture partition at t i-1 Saturated water vapor partial pressure at time, V Hj Represents the percentage of the jth humidity zone in the free volume of the containment vessel, V 0 Representing the free volume of the containment vessel, P 0 Denotes the test pressure,. DELTA.t denotes t i-1 To t i Time length of time, T ji Denotes the j temperature zone at t i Absolute temperature at time, T ji-1 Denotes the j temperature zone at t i-1 Absolute temperature of time, V Tj Representing the percentage of the jth temperature zone in the free volume of the containment vessel;
compensating volume leakage rate L under standard working condition environment Complement, N sigma The calculation formula of (c) is:
Figure BDA0003636969330000051
wherein n represents the number of time periods or cycles, m represents the number of humidity zones, i represents t i Time of day or t i-1 To t i Time period, j denotes the jth temperature zone or jth humidity zone, k denotes the number of temperature zones,H ji denotes the jth humidity division at t i Relative humidity at the moment H ji-1 Denotes the jth humidity zone at t i-1 Relative humidity at the moment of time, P Hji Denotes the jth humidity zone at t i Saturated partial pressure of water vapor at time, P Hji-1 Denotes the j moisture partition at t i-1 Saturated partial pressure of water vapour at time, V Hj Represents the percentage of the jth humidity zone in the free volume of the containment vessel, V 0 Representing the free volume of the containment vessel, P 0 Denotes the test pressure, P N Denotes the pressure in the environment of the standard working condition, and Δ t denotes t i-1 To t i Time length of time, T N Indicating the temperature, T, in the environment of the standard operating mode Hji-1 Denotes the jth humidity zone at t i-1 Absolute temperature at time, T Hji Denotes the jth humidity division at t i Absolute temperature at time, T ji Denotes the j temperature zone at t i Absolute temperature at time, T ji-1 Denotes the jth temperature partition at t i-1 Absolute temperature of moment, V Tj Represents the percentage of the j temperature zone to the free volume of the containment vessel.
Preferably, the calculating of the mass leakage rate according to the actual volume leakage rate specifically includes:
firstly, the mass leakage rate N in a delta t time is calculated ∑i The calculation formula is as follows:
Figure BDA0003636969330000052
wherein L is Real, N Σ i Indicating the standard working condition environment at t i Cumulative actual volumetric leakage rate, m, for all partitions at that time Qi (Qi) Represents the molar mass of air, m Water (W) Which represents the molar mass of the water vapor,
Figure BDA0003636969330000053
denotes t i The average water vapour partial pressure at the moment,
Figure BDA0003636969330000054
denotes t i-1 The average partial pressure of water vapour at the moment,
Figure BDA0003636969330000055
represents t i Mean partial pressure of water vapor at the time, R represents the ideal gas constant, P 0 Denotes the test pressure, P N Indicating the pressure, T, in the environment of standard conditions N Indicating the temperature under a standard working condition environment;
or:
Figure BDA0003636969330000061
wherein L is True, P Σ i Showing the test condition at t i Cumulative actual volumetric leakage rate, m, for all partitions at that time Qi (Qi) Represents the molar mass of air, m Water (W) Represents the molar mass of the water vapor,
Figure BDA0003636969330000062
represents t i The average partial pressure of water vapour at the moment,
Figure BDA0003636969330000063
denotes t i-1 Mean partial pressure of water vapor at the time, R represents the ideal gas constant, P 0 Denotes the test pressure, P N Indicating the pressure in the environment of the standard operating conditions,
Figure BDA0003636969330000064
represents t i To t i-1 Effective specific temperature in the containment over a period of time;
then calculating the total mass leakage rate M in a plurality of continuous delta t times ∑∑ Total mass leakage rate M ∑∑ The containment leakage rate is obtained by the calculation formula as follows:
Figure BDA0003636969330000065
wherein n is shown in the tableIndicating the number of time segments or cycles, i indicating t i Time of day or t i-1 To t i A time period.
Preferably, the compensation volume leakage rate is obtained by performing compensation calculation based on the average temperature of each temperature zone and the average humidity of each humidity zone in the containment, and the containment leakage rate is obtained by using the following calculation formula, where in the test working condition environment, the calculation formula of the containment leakage rate includes:
L true, P Σ i =L Measuring, P sigma i +L Complement, P ∑ i
Figure BDA0003636969330000066
Figure BDA0003636969330000067
Wherein L is True, P Σ i Showing the test condition at t i Cumulative actual volumetric leakage rate, L, for all partitions at that time Measure, P ∑ i Showing the test condition at t i Cumulative measured volume leakage rate, L, for all zones at that time Complement, P ∑ i Showing the test condition at t i The accumulated compensation volume leakage rate of all partitions at the moment, subscript p represents the test working condition, subscript i represents t i At the moment of time, the time of day,
Figure BDA0003636969330000071
is shown at t i The average relative humidity of all the zones at the time,
Figure BDA0003636969330000072
is shown at t i The average saturated water vapor partial pressure of all the partitions at that time,
Figure BDA0003636969330000073
is shown at t i-1 The average relative humidity of all the zones at the time,
Figure BDA0003636969330000074
is shown at t i-1 The average saturated water vapor partial pressure of all the zones at that time,
Figure BDA0003636969330000075
represents t i The average temperature at the time of day is,
Figure BDA0003636969330000076
represents t i-1 Mean temperature at time, V 0 Represents the free volume of the containment vessel, Δ t represents t i-1 To t i Time length of time, M P∑i Showing the test condition at t i Accumulated mass leakage rate, m, of all sub-areas at the moment Qi (Qi) Represents the molar mass of air, m Water (W) Represents the molar mass of water vapor, P 0 The pressure of the test is shown as,
Figure BDA0003636969330000077
represents t i Average water vapor partial pressure at the time, R represents an ideal gas constant;
under the standard working condition environment, the calculation formula of the containment leakage rate comprises the following steps:
L real, N Σ i =L Test, N Σ i +L Complement, N Σ i
Figure BDA0003636969330000078
Figure BDA0003636969330000079
Wherein L is Real N Σ i Indicating the standard working condition environment at t i Cumulative actual volumetric leakage rate, L, for all partitions at that time Test, N Σ i Indicating the standard working condition environment at t i Cumulative measured volume leakage rate, L, for all zones at that time Complement, N Σ i Indicating the standard working condition environment at t i The cumulative compensated volumetric leak rate for all segments at that time,
Figure BDA00036369693300000710
is shown at t i The average relative humidity of all the zones at the time,
Figure BDA00036369693300000711
is shown at t i The average water vapor partial pressure of all the zones at the time,
Figure BDA00036369693300000712
is shown at t i-1 The average relative humidity of all the zones at the time,
Figure BDA00036369693300000713
is shown at t i-1 The average water vapor partial pressure of all the zones at that time,
Figure BDA00036369693300000714
represents t i The average temperature at the time of day is,
Figure BDA00036369693300000715
represents t i-1 Mean temperature at time, V 0 Represents the free volume of the containment vessel, Δ t represents t i-1 To t i Time length of time, P 0 Denotes the test pressure, P N Indicating the pressure, T, in the environment of the standard operating conditions N Indicating the temperature, M, in the environment of the standard operating conditions N∑i Indicating the standard working condition environment at t i Accumulated mass leakage rate, m, of all sub-areas at the moment Qi (Qi) Represents the molar mass of air, m Water (W) Which represents the molar mass of the water vapor,
Figure BDA00036369693300000716
represents t i The average water vapor partial pressure at that time, R represents the ideal gas constant.
Preferably, the monitoring the temperature and/or humidity in the containment vessel and determining the compensated leakage rate according to the monitoring result of the temperature and/or humidity in the containment vessel comprise the following steps:
dividing the interior of the containment into a blocks, and respectively monitoring the temperature and the humidity of each block;
and calculating to obtain the compensation leakage rate according to the monitoring result of the temperature and the humidity of each block and the volume of each block.
Preferably, the calculating of the containment leakage rate includes:
calculating a blocks at t i To t i+1 Volumetric leak rate over time;
at t according to a blocks i To t i+1 Volume leakage rate of time period, calculating a blocks at t i To t i+1 And the mass leakage rate of the time period is used as the containment leakage rate.
Preferably, the a blocks are at t i To t i+1 Volume leakage rate of time period
Figure BDA0003636969330000081
The calculation formula of (c) is:
Figure BDA0003636969330000082
wherein a represents the number of partitions, L in,i+1 Denotes the t-th i+1 Inflation volume flow, T, at the outlet of the line constantly inflated with gas c,i+1,j Denotes the jth block at t i+1 Absolute temperature of gas in containment vessel, T, at time of day c,i,j Denotes the jth block at t i Absolute temperature of gas in safety enclosure at time m c,i,j Denotes the jth block at t i Mass of gas in containment vessel at time, R g,eq,i,j Denotes the jth block at t i Reduced gas constant, P, of the gas in the containment at time c Representing the pressure of the gas at the outlet/in-containment of the gas-filled line, at representing t i Time to t i+1 Time length of time, V c,i,j Represents the volume corresponding to the jth block, H c,i+1,j Denotes the jth block at t i+1 Relative humidity in containment at time H c,i,j Denotes the jth block at t i Relative humidity in Containment at time, f (T) c,i+1,j ) Denotes the jth block at t i+1 Saturated water vapor partial pressure at time f (T) c,i,j ) Denotes the jth block at t i+1 The saturated steam partial pressure in the containment vessel at the moment;
the a blocks are at t i To t i+1 Mass leakage rate of time period G out,i+1j Is calculated as:
Figure BDA0003636969330000091
wherein j represents the jth block, i represents the tth block i Time of day or t i To t i+1 Time period, a denotes the number of blocks, L in,i+1 Denotes the t-th i+1 The inflation volume flow, T, at the outlet of the line which is constantly inflated with gas c,i+1,j Denotes the jth block at t i+1 Absolute temperature of gas in containment at time, T c,i,j Denotes the jth block at t i Absolute temperature of gas in containment vessel at time m c,i,j Denotes the jth block at t i Mass of gas in containment vessel at time, R g,eq,i,j Denotes the jth block at t i Reduced gas constant, P, of the gas in the containment at time c Representing the pressure of the gas at the outlet/in-containment of the gas-filled line, at representing t i To t i+1 Time length of time, V c,i,j Represents the volume corresponding to the jth block, H c,i+1,j Denotes the jth block at t i+1 Instantaneous relative humidity in the containment, H c,i,j Denotes the jth block at t i Relative humidity in Containment at time, f (T) c,i+1,j ) Denotes the jth block at t i+1 Saturated partial pressure of water vapor in the containment at time f (T) c,i,j ) Denotes the jth block at t i+1 Saturated vapor partial pressure, R, in the containment at that time g,eq,i+1,j Denotes the jth block at t i+1 The reduced gas constant of the gas in the containment at the moment.
The containment leakage rate measuring method provided by the invention can be completed under the condition of keeping the internal pressure of the containment constant, and is a constant pressure method containment leakage rate measuring technology. In addition, the method fully pays attention to the influence of temperature change and humidity change on containment leakage rate measurement in the measurement process, and combines the characteristics of continuous change of temperature and humidity, and can give a proper acquisition period to calculate temperature compensation and humidity compensation, namely, the calculation model of the containment leakage rate in the method analyzes the gas actually leaked in the containment.
Drawings
FIG. 1 is a schematic diagram of a containment leak rate measurement method according to an embodiment of the invention;
FIG. 2 is a schematic diagram of a containment leak rate measurement method according to an embodiment of the invention;
FIG. 3 is a schematic diagram of a calculation flow of a containment leak rate measurement method according to an embodiment of the present disclosure;
FIG. 4 is a diagram illustrating the real-time variation of the temperature in the containment vessel within 24h in the present embodiment;
FIG. 5 is a diagram illustrating the real-time variation of the humidity in the containment vessel within 24h in the present embodiment;
FIG. 6 is a diagram showing the change in the pressure in the containment vessel in 24h in real time in the present embodiment;
FIG. 7 is a diagram showing the change in the flow rate of the gas supplied to the containment vessel in 24h in real time according to the present embodiment;
FIG. 8 shows the actual volume leakage rate L within 24h in this embodiment Real, N Σ i Real-time change maps of;
FIG. 9 shows the mass leakage rate M within 24h in this example ∑i Real-time variation diagram of (2).
Detailed Description
The technical solutions in the present invention will be described clearly and completely with reference to the accompanying drawings, and it is obvious that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the scope of the present invention.
Example 1
As shown in fig. 1, the present embodiment discloses a method for measuring a containment leakage rate, including:
filling gas into the containment vessel to enable the pressure in the containment vessel to reach a preset test pressure;
monitoring the real-time pressure in the containment, supplementing gas into the containment according to the monitoring result of the real-time pressure in the containment to maintain the real-time pressure in the containment equal to a preset test pressure, and detecting the flow of the gas supplemented to the containment to determine the measured volume leakage rate;
monitoring the temperature and/or humidity in the containment vessel, and determining the compensation leakage rate according to the monitoring result of the temperature and/or humidity in the containment vessel;
and determining the leakage rate of the containment according to the measured volume leakage rate and the compensation leakage rate.
As shown in fig. 2, the principle of the method is: dividing the containment into a plurality of partitions, and assuming that each partition is stable in an initial state: there is an empty leak port with a leak rate L 0 At the same time, there is a constant voltage of P 0 The infinite space of each subarea is continuously supplemented with compressed gas to ensure that the pressure in each subarea is always P 0 The flow rate of the supplemented compressed gas is L D (ii) a Let the free volume in each partition be V 0 Test pressure of P 0 The initial temperature of each zone is T 0 Initial humidity of each zone is H 0 The number of gas molecules in each partition is N 0 Under the condition that the parameters are stable,the measured leak rate is kept the same as the actual leak rate, i.e. L D =L 0 (ii) a In addition, in the test environment that each subarea is constant in pressure, although the pressure and the volume of each subarea are kept unchanged, the humidity H and the temperature T can change along with time, so the method also considers the influence of the humidity and the temperature on the measurement result, and measures the L measured on the charging pipeline D The volume compensation of temperature change and the volume compensation of humidity change in the containment (shell for short) are introduced on the basis of the reference value or the basic value of the containment leakage rate to obtain the compensated leakage rate. The calculation flow is shown in fig. 3.
Specifically, a pressurizing pipeline and a constant pressure test pipeline are arranged on the containment, gas (such as air and compressed air) is filled into the containment through the pressurizing pipeline to pressurize the containment normally, the pressure in the containment or the pressure in the pressurizing pipeline is detected and fed back, when the pressure in the containment is close to a preset pressure platform, the constant pressure test pipeline is switched to fill gas until the pressure in the containment reaches a preset test pressure, or gas is directly introduced into the containment through the pressurizing pipeline until the pressure in the containment reaches the preset test pressure, and then containment leakage rate measurement is started. In the measuring process, because the containment continuously leaks to the outside, theoretically, the pressure in the containment continuously drops, gas is supplemented into the containment through the constant-pressure testing pipeline, meanwhile, the real-time pressure in the containment is monitored, the flow of the supplementing gas is controlled according to the feedback result of the real-time pressure monitoring signal in the containment, the real-time pressure in the containment is always kept equal to the preset testing pressure, and the measured volume leakage rate is obtained by reading the flow of the supplementing gas.
It should be noted that, the control manner for maintaining the real-time pressure in the containment vessel equal to the preset test pressure may also be: using a pressure greater than the test pressure P 0 The gas source adjusts the flow of gas supplemented to the containment vessel in real time through the pressure feedback in the containment vessel shell, and continuously supplements the gas to the containment vessel, thereby maintaining the pressure in the containment vessel to be P all the time 0 . Compared with the control mode, the control mode can effectively shorten the time for reaching the preset test pressure for the first time and the period of pressure micro-amplitude fluctuation on the test platform.
Monitoring the temperature and/or humidity in the containment, and determining the compensation leakage rate according to the temperature and/or humidity monitoring result in the containment, wherein the method specifically comprises the following steps:
dividing the internal space of the containment into k virtual temperature partitions, respectively monitoring the temperature of each temperature partition to obtain real-time temperature data of each temperature partition at different moments, and/or dividing the internal space of the containment into m virtual humidity partitions, respectively monitoring the humidity of each humidity partition to obtain real-time humidity data of each humidity partition at different moments; and obtaining the compensation leakage rate according to the real-time temperature data of each temperature partition at different moments and/or according to the real-time humidity data of each humidity partition at different moments.
Specifically, the free space in the containment is divided into a plurality of (such as k) virtual temperature partitions by a temperature sensor arranged in the containment, wherein a typical area and a room with obstructed gas flow are taken as one partition, the volume coefficient of each temperature partition cannot exceed 0.1, the temperature of each temperature partition is monitored by the temperature sensor arranged in the containment to obtain the temperature change data of each temperature partition, temperature compensation (namely temperature correction) calculation is carried out in the k temperature partitions according to the temperature change data of each temperature partition, namely the temperature compensation of volume change is calculated, wherein the initial temperature of the j-th partition is set as T j0 At Δ t = t i -t i-1 The temperature is slowly changing over a period of time, t i Time temperature is set to T ji (ii) a The free space in the containment is divided into a plurality of (such as m) virtual humidity subareas by the humidity sensor arranged in the containment, wherein a typical area and a room with unsmooth gas flow are taken as one subarea, and the humidity of each humidity subarea is monitored by the humidity sensor arranged to obtain each humidity subareaHumidity change data of the humidity subareas are calculated in m subareas according to the humidity change data of each humidity subarea, namely, humidity compensation (namely humidity correction) calculation is carried out, namely, the compensation of the humidity on volume change is calculated, wherein the initial humidity of the jth subarea is set as H j0 At Δ t = t i -t i-1 The humidity changes slowly over a period of time, t i Time humidity is set to H ji
In this embodiment, the compensated leakage rate is specifically obtained by volume calculation, and the calculation step of the containment leakage rate includes: calculating a compensated volumetric leak rate as the compensated leak rate; calculating an actual volume leakage rate according to the compensated volume leakage rate and the measured volume leakage rate; and calculating the mass leakage rate according to the actual volume leakage rate to be used as the containment leakage rate.
Specifically, calculate the actual volume leakage rate according to compensation volume leakage rate and measurement volume leakage rate, specifically include: calculating the compensation volume leakage rate under the test working condition environment/the standard working condition environment according to the real-time humidity data of each humidity partition at different moments and the real-time temperature data of each temperature partition at different moments under the test working condition environment/the standard working condition environment; calculating the measured volume leakage rate under the test working condition environment/the standard working condition environment according to the flow of the gas supplemented into the containment vessel at each moment under the test working condition environment or the gas supply environment; and adding the compensation volume leakage rate under the test working condition environment/standard working condition environment and the measured volume leakage rate under the test working condition environment/standard working condition environment to obtain the actual volume leakage rate under the test working condition environment/standard working condition environment.
Taking the leakage rate compensation situation performed simultaneously in the m humidity zones and the k temperature zones as an example, several sets of calculation models are provided below, specifically as follows:
(1) First set of computational models
(1-1) calculating the compensated leak Rate
The compensation volume leakage rate is obtained by respectively carrying out compensation calculation on the basis of each temperature/humidity subarea, wherein:
in the test working condition environmentLower compensated volumetric leakage rate L Complement, P sigma The calculation formula of (A) is as follows:
Figure BDA0003636969330000141
wherein n represents the number of time periods or cycles, m represents the number of humidity zones, i represents t i Time of day or t i-1 To t i Time period, j denotes the jth temperature or humidity compartment, k denotes the number of temperature compartments, H ji Denotes the jth humidity zone at t i Relative humidity at the moment H ji-1 Denotes the jth humidity zone at t i-1 Relative humidity at the moment P Hji Denotes the jth humidity zone at t i Saturated water vapor partial pressure at time, P Hji-1 Denotes the j-th moisture partition at t i-1 Saturated water vapor partial pressure at time, V Hj Represents the percentage of the jth humidity zone in the free volume of the containment vessel, V 0 Representing the free volume of the containment vessel, P 0 Denotes the test pressure (0.42MPa. G, calculated absolute pressure 0.5213MPa. A, which is not described in detail below), and Δ t denotes t i-1 To t i Time length of time of day, T ji Denotes the jth temperature zone at t i Absolute temperature at time, T ji-1 Denotes the j temperature zone at t i-1 Absolute temperature of moment, V Tj Representing the percentage of the jth temperature zone in the free volume of the containment vessel;
compensated volume leakage rate L under standard working condition environment Complement, N sigma The calculation formula of (c) is:
Figure BDA0003636969330000142
wherein n represents the number of time periods or cycles, m represents the number of humidity zones, i represents t i Time of day or t i-1 To t i Time period, j denotes the jth temperature or humidity compartment, k denotes the number of temperature compartments, H ji Denotes the jth humidity division at t i Relative humidity at the moment H ji-1 Denotes the jth humidity division at t i-1 Relative humidity at the moment P Hji Denotes the jth humidity zone at t i Saturated water vapor partial pressure at time, P Hji-1 Denotes the j-th moisture partition at t i-1 Saturated water vapor partial pressure at time, V Hj Denotes the percentage of the jth humidity zone to the free volume of the containment vessel, V 0 Representing the free volume of the containment vessel, P 0 Denotes the test pressure, P N Denotes the pressure under the standard working condition (1 atmosphere, 0.1013MPa. A, which is not described in detail below), and Δ t denotes t i-1 To t i Time length of time, T N Indicates the temperature (0 ℃, 273.15K for absolute temperature, and is not described in detail below) under the standard working condition environment, T Hji-1 Denotes the jth humidity zone at t i-1 Absolute temperature at time, T Hji Denotes the jth humidity zone at t i Absolute temperature at time, T ji Denotes the jth temperature zone at t i Absolute temperature at time, T ji-1 Denotes the j temperature zone at t i-1 Absolute temperature of time, V Tj Represents the percentage of the j temperature zone to the free volume of the containment vessel.
(1-2) calculation of measured volume leakage Rate
The volume leakage rate is generally measured by directly and cumulatively measuring a flowmeter with a temperature and pressure measuring function, but for a flowmeter without the cumulative function, the volume leakage rate is calculated according to the following calculation formula:
measuring volume leakage rate L under standard working condition environment Measuring, N sigma The calculation formula of (c) is:
Figure BDA0003636969330000151
the calculation formula is obtained by converting the test working condition environment in the containment into the standard working condition environment, wherein n represents the number of time periods or the number of cycles, i represents t i Time of day or t i-1 To t i Period of time, L Measure, P ∑ i Showing the test condition at t i Measurement of the accumulation of all zones at a timeVolume leakage rate (i.e. at t under test condition environment) i The flow rate of gas to be supplied to the containment vessel detected at that moment), P) 0 Denotes the test pressure, P N Indicating the pressure, T, in the environment of the standard operating conditions N Which represents the temperature in the environment of the standard operating condition,
Figure BDA0003636969330000152
represents t i-1 To t i The effective specific temperature in the containment in the time zone,
Figure BDA0003636969330000153
alternatively, the first and second liquid crystal display panels may be,
Figure BDA0003636969330000154
the calculation formula is obtained by converting the gas supply environment of the pressurizing pipeline into the standard working condition environment, wherein n represents the number of time segments or the number of cycles, and i represents t i Time of day or t i-1 To t i Time period, L Measuring, A sigma i Indicating at t in the environment of the supply of air i Cumulative volumetric leak rate measured for all zones at time (i.e. at t in a gas supply environment i The flow value of the gas replenished to the containment vessel detected at the moment), P) Ai Representing the real-time pressure, P, at the line charging the containment in the gas supply environment N Indicating the pressure, T, in the environment of the standard operating conditions N The temperature in the environment of the standard working condition is shown,
Figure BDA0003636969330000161
denotes the gas supply environment at t i-1 To t i Effective specific temperature in the containment over a period of time.
Measuring volume leakage rate L under test working condition environment Measuring, p sigma The calculation formula of (c) is:
Figure BDA0003636969330000162
wherein L is Measure, P ∑ i Showing the test condition at t i The cumulative volume leakage rate of all the subareas at the moment (namely t under the test working condition environment) i The flow value of the gas replenished to the containment vessel detected at the moment).
(1-3) calculation of actual volume leakage Rate
Actual volumetric leakage rate L Fruit of Chinese wolfberry The actual volume leakage rate L under the test working condition environment is equal to the sum of the measured volume leakage rate and the compensation volume leakage rate Real, P sigma The calculation formula of (A) is as follows:
L real, P sigma =L Measuring, P sigma +L Complement, P sigma
Actual volume leakage rate L under standard working condition environment Real, N sigma The calculation formula of (c) is:
L real, N sigma =L Measuring, N sigma +L Complement, N sigma
Among them, it should be noted that: in this embodiment, the subscript N represents a standard working condition environment, the subscript P represents a test working condition environment, and all the partitions of the subscript Σ table at all times are accumulated, which is not described in detail below.
(1-4) calculation of Mass leakage Rate
Calculating mass leakage rate according to actual volume leakage rate, wherein the calculating step of mass leakage rate comprises the following steps:
firstly, the mass leakage rate M in a delta t time is calculated ∑i The calculation formula is as follows:
Figure BDA0003636969330000163
wherein the subscript ∑ i denotes t i All the subareas at the moment are accumulated, and the subscript N sigma i represents t under the standard working condition environment i All partitions are accumulated at all times, L Real, N Σ i Indicating the standard working condition environment at t i Cumulative actual volumetric leakage rate, m, for all zones at any time Qi (Qi) Represents the molar mass of air, m Water (W) Which represents the molar mass of the water vapor,
Figure BDA0003636969330000164
represents t i The average partial pressure of water vapour at the moment,
Figure BDA0003636969330000171
represents t i-1 The average water vapour partial pressure at the moment,
Figure BDA0003636969330000172
represents t i The average partial pressure of water vapor at that time, R represents an ideal gas constant (8.314J. Mol.) -1 ·k -1 ),P 0 Representing the test pressure, P N Indicating the pressure, T, in the environment of the standard operating conditions N Indicating the temperature under the environment of standard working conditions;
alternatively, the first and second electrodes may be,
Figure BDA0003636969330000173
wherein the subscript ∑ i denotes t i All partitions are accumulated at the moment, and subscript P sigma i represents t under the test working condition environment i All partitions are accumulated at all times, L True, P Σ i Showing the test condition at t i Cumulative actual volumetric leakage rate, m, for all zones at any time Qi (Qi) Represents the molar mass of air, m Water (W) Which represents the molar mass of the water vapor,
Figure BDA0003636969330000174
denotes t i The average water vapour partial pressure at the moment,
Figure BDA0003636969330000175
represents t i-1 Average water vapor partial pressure at the time, R represents an ideal gas constant, P 0 Denotes the test pressure, P N The pressure in the environment of the standard working condition is shown,
Figure BDA0003636969330000176
represents t i To t i-1 The effective specific temperature within the containment over a period of time,
Figure BDA0003636969330000177
then calculating the total mass leakage rate M in a plurality of continuous delta t times ∑∑ The total mass leakage rate is containment leakage rate of the method, and the calculation formula is as follows:
Figure BDA0003636969330000178
wherein n represents the number of time segments or cycles, and i represents t i Time of day or t i-1 To t i A period of time.
(2) Second set of computational models
Compared with the first group of calculation models, the calculation process is properly simplified, the compensation volume leakage rate is obtained by performing compensation calculation based on the average temperature of each temperature partition and the average humidity of each humidity partition in the containment, and the containment leakage rate is obtained by adopting the following calculation formula:
under the test working condition environment, the calculation formula of the containment leakage rate comprises the following steps:
L true, P ∑ i =L Measuring, P sigma i +L Complement, P ∑ i
Figure BDA0003636969330000179
Figure BDA0003636969330000181
Wherein L is True, P ∑ i Showing the test condition at t i Cumulative actual volumetric leakage rate, L, for all partitions at that time Measure, P ∑ i Showing the test condition at t i The volume leakage rate is measured by the accumulation of all the subareas at the moment (namely at t under the test working condition environment) i The flow value of the gas replenished to the containment vessel detected at the moment), L) Complement, P ∑ i Showing the test condition at t i The accumulated compensation volume leakage rate of all the subareas at the moment, subscript p represents the test working condition, subscript i represents t i At the moment of time, the time of day,
Figure BDA0003636969330000182
is shown at t i The average relative humidity of all the zones at a time,
Figure BDA0003636969330000183
is shown at t i The average saturated water vapor partial pressure of all the partitions at that time,
Figure BDA0003636969330000184
is shown at t i-1 The average relative humidity of all the zones at the moment,
Figure BDA0003636969330000185
is shown at t i-1 The average saturated water vapor partial pressure of all the zones at that time,
Figure BDA0003636969330000186
represents t i The average temperature at the time of day is,
Figure BDA0003636969330000187
represents t i-1 Mean temperature at time, V 0 Represents the free volume of the containment vessel, Δ t represents t i-1 To t i Time length of time, M P∑i Showing the test condition at t i Accumulated mass leakage rate, m, of all sub-areas at that moment Qi (Qi) Represents the molar mass of air, m Water (W) Represents the molar mass of water vapor, P 0 The pressure of the test is shown as,
Figure BDA0003636969330000188
denotes t i The average water vapor partial pressure at the moment, R representing the ideal gas constant;
under the standard working condition environment, the calculation formula of the containment leakage rate comprises the following steps:
L real, N Σ i =L Test, N Σ i +L Complement, N Σ i
Figure BDA0003636969330000189
Figure BDA00036369693300001810
Wherein L is Real N Σ i Indicating the standard working condition environment at t i Cumulative actual volumetric leakage rate, L, for all partitions at that time Measuring, N ∑ i Indicating the standard working condition environment at t i Cumulative measured volumetric leakage rate, L, for all zones at that time Complement, N Σ i Indicating the standard working condition environment at t i The cumulative compensated volumetric leak rate for all zones at that time,
Figure BDA00036369693300001811
is shown at t i The average relative humidity of all the zones at the moment,
Figure BDA00036369693300001812
is shown at t i The average saturated water vapor partial pressure of all the partitions at the time,
Figure BDA00036369693300001813
is shown at t i-1 The average relative humidity of all the zones at the time,
Figure BDA00036369693300001814
is shown at t i-1 The average saturated water vapor partial pressure of all the partitions at the time,
Figure BDA00036369693300001815
represents t i The average temperature at the time of day is,
Figure BDA00036369693300001816
represents t i-1 Mean temperature at time, V 0 Denotes the free volume of the containment vessel, Δ t denotes t i-1 To t i Time length of time, P 0 Indicating testPressure, P N Indicating the pressure, T, in the environment of standard conditions N Denotes the temperature, M, in the environment of the standard operating mode N∑i Indicating the standard working condition environment at t i Accumulated mass leakage rate, m, of all sub-areas at that moment Qi (Qi) Represents the molar mass of air, m Water (I) Which represents the molar mass of the water vapor,
Figure BDA0003636969330000191
represents t i The average water vapor partial pressure at that time, R, represents the ideal gas constant.
It should be noted that the actual average pressure in the containment vessel is around P 0 Fluctuating with a slight amplitude. Therefore, P in each calculation formula in the present embodiment 0 Can also be replaced by a real-time measured pressure P in the safety housing i (relative change before and after replacement is only about five parts per million), and, with P i The calculation is more consistent with the actual change process of the pressure in the containment vessel, and the result is more accurate.
Below at 1000m 3 Example of a containment vessel, wherein the test pressure P 0 The first set of calculation model calculation processes are detailed as follows, wherein (absolute pressure) is 533.65kPa, temperature partitions are 28, the arrangement condition of each temperature sensor is shown in table 1, humidity partitions are 10, and the arrangement condition of each humidity sensor is shown in table 2:
TABLE 1
Figure BDA0003636969330000192
Figure BDA0003636969330000201
TABLE 2
Figure BDA0003636969330000202
Because temperature, humidity, pressure in the containment are real-time change in this system measurement process, the flow of the supplementary gas on the admission pipeline also can take place corresponding change, accomplish the measurement and the analysis of containment leakage rate through detectors such as pressure gauge in the containment, temperature sensor, humidity transducer in 24h of continuous collection, wherein:
the real-time change of the temperature in the containment vessel within 24h is shown in fig. 4;
the real-time change of the humidity in the containment vessel within 24h is shown in fig. 5;
the real-time change of the pressure in the containment vessel within 24h is shown in fig. 6;
the real-time change in the flow of the gas to the containment vessel in 24h is shown in fig. 7.
Taking the first time period Δ t =10s (i.e., 0 to 10 s) within 24h as an example, the volumetric leakage rate L at a single moment is illustrated Real, N Σ i Mass leakage rate M at sum time ∑i The calculation of (c):
the temperature value data (single-bit ℃) acquired by the T1-T28 temperature sensor at the ti-1 moment are as follows in sequence:
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 humidity sensor at t i-1 The humidity value data (%) acquired at the moment are as follows:
48.1898,46.8971,47.1527,47.7582,46.3522,48.7407,46. 4725,46.7261,46.1951,48.9654;
at t i-1 The real-time pressure in the containment vessel acquired at any moment is as follows: 533.652kpa;
at t i-1 The flow value of the gas in the pipeline for filling the gas into the containment vessel is acquired at any time, namely the measured leakage rate (standard working condition) is as follows: 0.7945m 3 /h;
T1-T28 temperature sensor at T i Time (t) i -t i-1 =10 s) the temperature value data (unit ℃) obtained by collection were:
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 humidity sensor at t i The humidity value data (%) collected at the time are as follows:
48.1968,46.8991,47.1556,47.7602,46.3552,48.7447,46. 4715,46.7281,46.1990,48.9595;
at t i The real-time pressure in the containment vessel acquired at any moment is as follows: 533.652kpa;
at t i The flow value of the gas in the pipeline for filling the gas into the containment vessel is acquired at any time, namely the measured leakage rate (standard working condition) is as follows: 0.7975m 3 /h;
Substituting the acquired data into the above calculation formula to obtain:
t i the effective specific temperature in the containment at the moment is as follows:
Figure BDA0003636969330000221
at t i-1 To t i The measured volumetric leak rate over the time period was:
L test, N Σ i =0.796m 3 /h。
According to the calculation model in example 1, respectively:
at t i-1 To t i The compensated leak rate over the time period is:
L complement, N Σ i =3.16703m 3 /h。
At t i-1 To t i The actual volumetric leak rate over the time period is:
L real, N Σ i =L Test, N Σ i +L Complement, N Σ i =3.96268m 3 /h
Passing through t i-1 The temperature and the humidity in the containment vessel at the moment can be calculated to obtain t i-1 The partial pressure of vapor in the containment vessel at the moment is as follows:
Figure BDA0003636969330000222
passing through t i The temperature and the humidity in the containment vessel at the moment can be calculated to obtain t i The partial pressure of water vapor in the safety shell at the moment is as follows:
Figure BDA0003636969330000223
at t i-1 To t i The actual mass leakage rate over the time period is:
Figure BDA0003636969330000224
dividing the acquired data in 24h into 24 multiplied by 360 time periods of 10s, and obtaining the measured volume leakage rate L of 24 multiplied by 360 groups of single time according to the single time calculation process in each time period Measuring, N ∑ i Compensating for volume leakage rate L Complement, N Σ i Actual volume leakage rate L Real, N Σ i And mass leakage rate M ∑i . Here, the calculation in other single-time periods is not repeated, and only the actual volume leakage rate L in 24h is illustrated by the figure Real, N Σ i And mass leakage rate M ∑i A process of variation wherein:
actual volume leakage rate L within 24h Real, N Σ i The real-time variation of (a) is shown in fig. 8.
Mass leakage rate M within 24h ∑i The real-time variation of (a) is shown in fig. 9.
Volumetric leak rate L measured at a single time for all periods within 24h Measuring, N ∑ i Accumulating and calculating the average value, and calculating to obtain the measured volume leakage rate within 24h as follows:
L measuring, N sigma =1.05534m 3 /h
Compensating for volumetric leak rate L for a single instant over all periods of 24h Complement, N Σ i And accumulating and calculating an average value, and calculating to obtain a compensation volume leakage rate within 24h as follows:
L measuring, N sigma =0.15381m 3 /h
Safety casing volume 1000m 3 The relative volume leakage rate is 0.617553 percent of the total volume of gas in the safety shell per day.
The average actual volumetric leak rate over 24h is:
L real, N sigma =L Measuring, N sigma +L Complement, N sigma =1.20915m 3 /h
Mass leakage rate M for a single moment of all periods within 24h ∑i And accumulating and calculating the average value, wherein the calculated mass leakage rate of 24h is as follows:
Figure BDA0003636969330000231
at the time of 0, the total mass of the gas in the containment vessel is 6063.65Kg, and the relative mass leakage rate is 0.618286 percent of the total mass of the gas in the containment vessel per day.
The containment leakage rate measuring method of the embodiment can be completed under the condition of maintaining the pressure inside the containment constant, and is a constant pressure method containment leakage rate measuring technology. In addition, the method fully focuses on the influence of temperature change and humidity change on containment leakage rate measurement in the measurement process, and combines the characteristics of continuous change of temperature and humidity, and can give a proper acquisition period to calculate temperature compensation and humidity compensation, namely, the calculation model of the containment leakage rate in the method analyzes the gas actually leaked in the containment.
Example 2
The embodiment discloses a containment vessel leakage rate measuring method, which is different from the containment vessel leakage rate measuring method described in embodiment 1 in that:
the method comprises dividing the interior of the containment into a blocks, and respectively monitoring the temperature and humidity of each block; and calculating to obtain the compensated leakage rate according to the monitoring result of the temperature and the humidity of each block and the volume of each block.
The method adopts a third group of calculation models which are different from the first group of calculation models and the second group of calculation models in the embodiment 1 (3), and the specific calculation steps comprise:
calculating a blocks at t i To t i+1 Volumetric leak rate over time;
at t according to a blocks i To t i+1 Volume leakage rate of time period, calculating a blocks at t i To t i+1 And (3) mass leakage rate of the time period, wherein the mass leakage rate is the containment leakage rate measured by the method.
Specifically, a blocks are at t i To t i+1 Volume leakage rate of time period
Figure BDA0003636969330000241
The calculation formula of (A) is as follows:
Figure BDA0003636969330000242
wherein a represents the number of partitions, L in,i+1 Denotes the t-th i+1 The inflation volume flow, T, at the outlet of the line which is constantly inflated (i.e. the constant pressure test line in this context) c,i+1,j Denotes the jth block at t i+1 Absolute temperature of gas in containment at time, T c,i,j Denotes the jth block at t i Absolute temperature of gas in containment vessel at time m c,i,j Denotes the jth block at t i Mass of gas (mixture of gases) in containment at time, R g,eq,i,j Denotes the jth block at t i Reduced gas constant, P, of the gas (mixture of gases) in the containment at that moment c Representing the pressure of the gas at the outlet/in-containment of the gas-filled line, at representing t i Time to t i+1 Time length of time, V c,i,j Denotes the volume corresponding to the jth block, H c,i+1,j Denotes the jth block at t i+1 Relative humidity in containment vessel at time H c,i,j Denotes the jth block at t i Relative humidity in Containment at time, f (T) c,i+1,j ) Denotes the jth block at t i+1 Saturated partial pressure of water vapor in the containment at time f (T) c,i,j ) Denotes the jth block at t i+1 The saturated steam partial pressure in the containment vessel at the moment;
a blocks are at t i To t i+1 Mass leakage rate of time segment G out,i+1j The calculation formula of (A) is as follows:
Figure BDA0003636969330000251
wherein j represents the jth block, i represents the tth block i Time of day or t i To t i+1 Time period, a denotes the number of blocks, L in,i+1 Denotes the t-th i+1 The inflation volume flow, T, at the outlet of the line which is constantly inflated with gas c,i+1,j Denotes the jth block at t i+1 Absolute temperature of gas in containment vessel, T, at time of day c,i,j Denotes the jth block at t i Absolute temperature of gas in containment vessel, m, at time of day c,i,j Denotes the jth block at t i Mass of gas (mixture of gases) in containment at time, R g,eq,i,j Denotes the jth block at t i Reduced gas constant, P, of the gas (mixed gas) in the containment at that moment c Indicating the pressure of the gas at the outlet of the gas-filled circuit/inside the containment, at i Time to t i+1 Time length of time, V c,i,j Represents the volume, H, corresponding to the jth block c,i+1,j Denotes the jth block at t i+1 Relative humidity within the safety shell at the moment, H c,i,j Denotes the jth block at t i Relative humidity in containment at time, f (T) c,i+1,j ) Denotes the jth block at t i+1 Partial pressure of saturated vapor in containment at time f (T) c,i,j ) Denotes the jth block at t i+1 Saturated partial pressure of water vapor in containment vessel at time, R g,eq,i+1,j Denotes the jth block at t i+1 The gas in the containment vessel at the moment is reduced to a gas constant.
The containment leakage rate measuring method can be completed under the condition of maintaining the pressure inside the containment constant, and is a constant pressure method containment leakage rate measuring technology. In addition, the method fully focuses on the influence of temperature change and humidity change on containment leakage rate measurement in the measurement process, and combines the characteristics of continuous change of temperature and humidity, and can give a proper acquisition period to calculate temperature compensation and humidity compensation, namely, the calculation model of the containment leakage rate in the method analyzes the gas actually leaked in the containment.
It will be understood that the above embodiments are merely exemplary embodiments taken to illustrate the principles of the present invention, which is not limited thereto. It will be apparent to those skilled in the art that various modifications and improvements can be made without departing from the spirit and substance of the invention, and these modifications and improvements are also considered to be within the scope of the invention.

Claims (11)

1. A containment leakage rate measuring method comprises the following steps:
filling gas into the containment vessel to enable the pressure in the containment vessel to reach the preset test pressure;
monitoring the real-time pressure in the containment, supplementing gas into the containment according to the monitoring result of the real-time pressure in the containment to maintain the real-time pressure in the containment equal to a preset test pressure, and detecting the flow of the gas supplemented to the containment to determine the measured volume leakage rate;
monitoring the temperature and/or humidity in the containment vessel, and determining the compensation leakage rate according to the monitoring result of the temperature and/or humidity in the containment vessel;
and determining the leakage rate of the containment according to the measured volume leakage rate and the compensation leakage rate.
2. The containment leak rate measurement method according to claim 1, wherein the monitoring of the temperature and/or humidity in the containment vessel and the determining of the compensated leak rate from the monitoring of the temperature and/or humidity in the containment vessel comprises the following steps:
dividing the internal space of the containment into k virtual temperature zones, respectively monitoring the temperature of each temperature zone to obtain real-time temperature data of each temperature zone at different moments, and/or,
dividing the internal space of the containment into m virtual humidity zones, and monitoring the humidity of each humidity zone respectively to obtain real-time humidity data of each humidity zone at different moments;
and obtaining the compensated leakage rate according to the real-time temperature data of each temperature partition at different moments and/or according to the real-time humidity data of each humidity partition at different moments.
3. The method for measuring the leakage rate of the containment vessel according to claim 2, wherein the compensated leakage rate is obtained by volume calculation, and the calculation step of the leakage rate of the containment vessel comprises:
calculating a compensated volumetric leak rate as the compensated leak rate;
calculating an actual volume leakage rate according to the compensated volume leakage rate and the measured volume leakage rate;
and calculating the mass leakage rate according to the actual volume leakage rate to be used as the containment leakage rate.
4. The method for measuring the containment leakage rate according to claim 3, wherein the calculating the actual volume leakage rate according to the compensated volume leakage rate and the measured volume leakage rate specifically comprises:
calculating the compensation volume leakage rate under the test working condition environment/the standard working condition environment according to the real-time humidity data of each humidity partition at different moments and the real-time temperature data of each temperature partition at different moments under the test working condition environment/the standard working condition environment;
calculating the measured volume leakage rate under the test working condition environment/the standard working condition environment according to the flow of the gas supplemented into the containment vessel at each moment under the test working condition environment or the gas supply environment;
and adding the compensation volume leakage rate under the test working condition environment/standard working condition environment and the measurement volume leakage rate under the test working condition environment/standard working condition environment to obtain the actual volume leakage rate under the test working condition environment/standard working condition environment.
5. Method for measuring containment leak rate according to claim 4, wherein the volumetric leak rate L is measured in a standard operating environment Measuring, N sigma The calculation formula of (c) is:
Figure RE-FDA0003865889630000021
wherein n represents the number of time segments or cycles, and i represents t i Time of day or t i-1 To t i Period of time, L Measure, P ∑ i Showing the test condition at t i Cumulative measured volume leakage rate, P, for all zones at that time 0 Denotes the test pressure, P N Indicating the pressure, T, in the environment of standard conditions N The temperature in the environment of the standard operating condition is shown,
Figure RE-FDA0003865889630000022
represents t i-1 To t i The effective specific temperature within the containment over a period of time,
alternatively, the first and second electrodes may be,
Figure RE-FDA0003865889630000031
wherein n represents the number of time segments or cycles, and i represents t i Time of day or t i-1 To t i Period of time, L Measuring, A sigma i Indicating at t in the environment of the supply of air i Cumulative measured volumetric leakage rate, P, for all zones at any time Ai Represents the real-time pressure, P, in the line that fills the containment vessel with gas in the gas supply environment N Indicating the pressure, T, in the environment of standard conditions N Which represents the temperature in the environment of the standard operating condition,
Figure RE-FDA0003865889630000032
denotes the gas supply environment at t i-1 To t i Effective specific temperature in the containment over a period of time;
measuring volume leakage rate L under test working condition environment Measuring, p sigma The calculation formula of (c) is:
Figure RE-FDA0003865889630000033
wherein L is Measure, P ∑ i Showing the test condition at t i The cumulative measured volumetric leak rate for all zones at that time.
6. The containment leak rate measurement method of claim 4, wherein the compensated volumetric leak rate is calculated by performing compensation based on each temperature/humidity zone, wherein,
compensation volume leakage rate L under test working condition environment Complement, P sigma The calculation formula of (A) is as follows:
Figure RE-FDA0003865889630000034
wherein n represents the number of time periods or cycles, m represents the number of humidity zones, i represents t i Time of day or t i-1 To t i Time period, j denotes the jth temperature zone or jth humidity zone, k denotes the number of temperature zones, H ji Denotes the jth humidity zone at t i Relative humidity at the moment H ji-1 Denotes the jth humidity division at t i-1 Relative humidity at the moment P Hji Denotes the jth humidity zone at t i Saturated water vapor partial pressure at time, P Hji-1 Denotes the jth humidity zone at t i-1 Saturated partial pressure of water vapour at time, V Hj Represents the percentage of the jth humidity zone in the free volume of the containment vessel, V 0 Representing the free volume of the containment vessel, P 0 Denotes the test pressure,. DELTA.t denotes t i-1 To t i Time length of time of day, T ji Denotes the jth temperature zone at t i Absolute temperature at time, T ji-1 Denotes the jth temperature zone at t i-1 Absolute temperature of time, V Tj Representing the percentage of the jth temperature zone in the free volume of the containment vessel;
compensated volume leakage rate L under standard working condition environment Complement, N sigma The calculation formula of (A) is as follows:
Figure RE-FDA0003865889630000041
wherein n represents the number of time periods or cycles, m represents the number of humidity zones, i represents t i Time of day or t i-1 To t i Time period, j denotes the jth temperature zone or jth humidity zone, k denotes the number of temperature zones, H ji Denotes the jth humidity zone at t i Relative humidity of time of day, H ji-1 Denotes the jth humidity zone at t i-1 Relative humidity at the moment P Hji Denotes the jth humidity zone at t i Saturated partial pressure of water vapor at time, P Hji-1 Denotes the jth humidity zone at t i-1 Saturated partial pressure of water vapour at time, V Hj Represents the percentage of the jth humidity zone in the free volume of the containment vessel, V 0 Representing the free volume of the containment vessel, P 0 Denotes the test pressure, P N Denotes the pressure in the environment of the standard working condition, and Δ t denotes t i-1 To t i Time length of time, T N Indicating the temperature, T, in the environment of the standard operating mode Hji-1 Denotes the jth humidity zone at t i-1 Absolute temperature at time, T Hji Denotes the jth humidity zone at t i Absolute temperature at time, T ji Denotes the jth temperature zone at t i Absolute temperature at time, T ji-1 Denotes the jth temperature zone at t i-1 Absolute temperature of moment, V Tj Represents the percentage of the j temperature zone to the free volume of the containment vessel.
7. The containment leak rate measurement method according to claim 4, wherein the calculating of the mass leak rate from the actual volume leak rate specifically includes:
firstly, the mass leakage rate M in a delta t time is calculated ∑i The calculation formula is as follows:
Figure RE-FDA0003865889630000051
wherein L is Real, N Σ i Indicating the standard working condition environment at t i Cumulative actual volume leakage rate, m, for all zones at that time Qi (Qi) Represents the molar mass of air, m Water (W) Which represents the molar mass of the water vapor,
Figure RE-FDA0003865889630000052
represents t i The average partial pressure of water vapour at the moment,
Figure RE-FDA0003865889630000053
represents t i-1 Mean water vapor fraction at timeThe pressure is applied to the inner wall of the cylinder,
Figure RE-FDA0003865889630000054
represents t i Mean partial pressure of water vapor at the time, R represents the ideal gas constant, P 0 Denotes the test pressure, P N Indicating the pressure, T, in the environment of the standard operating conditions N The temperature in the environment of the standard operating condition is shown,
alternatively, the first and second electrodes may be,
Figure RE-FDA0003865889630000055
wherein L is True, P Σ i Showing the test condition at t i Cumulative actual volume leakage rate, m, for all zones at that time Qi (Qi) Represents the molar mass of air, m Water (W) Represents the molar mass of the water vapor,
Figure RE-FDA0003865889630000056
denotes t i The average partial pressure of water vapour at the moment,
Figure RE-FDA0003865889630000057
represents t i-1 Mean partial pressure of water vapor at the time, R represents the ideal gas constant, P 0 Denotes the test pressure, P N The pressure in the environment of the standard working condition is shown,
Figure RE-FDA0003865889630000058
represents t i To t i-1 Effective specific temperature within the containment over a period of time;
then calculating the total mass leakage rate M in a plurality of continuous delta t times ∑∑ Total mass leakage rate M ∑∑ The containment leakage rate is obtained by the calculation formula as follows:
Figure RE-FDA0003865889630000059
wherein n represents the number of time segments or cycles, and i represents t i Time of day or t i-1 To t i A time period.
8. The containment leakage rate measuring method according to claim 4, wherein the compensated volume leakage rate is obtained by performing compensation calculation based on the average temperature of each temperature zone and the average humidity of each humidity zone in the containment, and the containment leakage rate is obtained by using the following calculation formula,
under the test working condition environment, the calculation formula of the containment leakage rate comprises the following steps:
L true, P Σ i =L Measure, P ∑ i +L Complement, P ∑ i
Figure RE-FDA0003865889630000061
Figure RE-FDA0003865889630000062
Wherein L is True, P ∑ i Showing the test condition at t i Cumulative actual volumetric leakage rate, L, for all zones at that time Measure, P ∑ i Showing the test condition at t i Cumulative measured volumetric leakage rate, L, for all zones at that time Complement, P ∑ i Showing the test condition at t i The accumulated compensation volume leakage rate of all the subareas at the moment, the subscript p represents the test working condition, the subscript i represents t i At the moment of time, the time of day,
Figure RE-FDA0003865889630000063
is shown at t i The average relative humidity of all the zones at the moment,
Figure RE-FDA0003865889630000064
is shown at t i The average saturated water vapor partial pressure of all the partitions at that time,
Figure RE-FDA0003865889630000065
is shown at t i-1 The average relative humidity of all the zones at the moment,
Figure RE-FDA0003865889630000066
is shown at t i-1 The average saturated water vapor partial pressure of all the partitions at that time,
Figure RE-FDA0003865889630000067
represents t i The average temperature at the time of day is,
Figure RE-FDA0003865889630000068
represents t i-1 Mean temperature at time, V 0 Represents the free volume of the containment vessel, Δ t represents t i-1 To t i Time length of time, M P∑i Showing the test condition at t i Accumulated mass leakage rate of all sub-areas at any moment,m Qi (Qi) Represents the molar mass of air, m Water (W) Denotes the molar mass of water vapor, P 0 The pressure of the test is shown as,
Figure RE-FDA0003865889630000069
represents t i The average water vapor partial pressure at that time, R representing the ideal gas constant;
under the standard working condition environment, the calculation formula of the containment leakage rate comprises the following steps:
L real, N Σ i =L Test, N Σ i +L Complement, N Σ i
Figure RE-FDA00038658896300000610
Figure RE-FDA00038658896300000611
Wherein L is Real N Σ i Indicating the standard working condition environment at t i Cumulative actual volumetric leakage rate, L, for all zones at that time Measuring, N ∑ i Indicating the standard working condition environment at t i Cumulative measured volumetric leakage rate, L, for all zones at that time Complement, N Σ i Indicating the standard working condition environment at t i The cumulative compensated volumetric leak rate for all zones at that time,
Figure RE-FDA00038658896300000612
is shown at t i The average relative humidity of all the zones at the moment,
Figure RE-FDA0003865889630000071
is shown at t i The average water vapor partial pressure of all the zones at the moment,
Figure RE-FDA0003865889630000072
is shown at t i-1 The average relative humidity of all the zones at the moment,
Figure RE-FDA0003865889630000073
is shown at t i-1 The average water vapor partial pressure of all the zones at the moment,
Figure RE-FDA0003865889630000074
represents t i The average temperature at the time of day is,
Figure RE-FDA0003865889630000075
represents t i-1 Mean temperature at time, V 0 Represents the free volume of the containment vessel, Δ t represents t i-1 To t i Time length of time, P 0 Denotes the test pressure, P N Indicating the pressure, T, in the environment of standard conditions N Indicating the temperature, M, in the environment of the standard operating conditions N∑i Indicating the standard working condition environment at t i Accumulated mass leakage rate, m, of all sub-areas at that moment Qi (Qi) Represents the molar mass of air, m Water (I) Represents the molar mass of the water vapor,
Figure RE-FDA0003865889630000076
represents t i The average water vapor partial pressure at the time, R, represents the ideal gas constant.
9. The containment leak rate measurement method according to claim 1, wherein the monitoring of the temperature and/or humidity in the containment vessel and the determining of the compensated leak rate from the monitoring of the temperature and/or humidity in the containment vessel comprises the following steps:
dividing the containment into a blocks, and respectively monitoring the temperature and the humidity of each block;
and calculating to obtain the compensation leakage rate according to the monitoring result of the temperature and the humidity of each block and the volume of each block.
10. The containment leak rate measurement method of claim 9, wherein the containment leak rate calculation step comprises:
calculating a blocks at t i To t i+1 Volumetric leak rate over time;
at t according to a blocks i To t i+1 Volume leak rate of time segment a blocks are calculated at t i To t i+1 And taking the mass leakage rate of the time period as the containment leakage rate.
11. The method of measuring containment leak rate of claim 10, wherein the a blocks are at t i To t i+1 Volumetric leak rate of time period
Figure RE-FDA0003865889630000077
The calculation formula of (c) is:
Figure RE-FDA0003865889630000081
wherein a represents the number of partitions, L in,i+1 Denotes the t-th i+1 Inflation volume flow, T, at the outlet of the line which is constantly inflated c,i+1,j Denotes the jth block at t i+1 Absolute temperature of gas in containment vessel, T, at time of day c,i,j Denotes the jth block at t i Absolute temperature of gas in containment vessel, m, at time of day c,i,j Denotes the jth block at t i Mass of gas in containment vessel at time of day, R g,eq,i,j Denotes the jth block at t i Reduced gas constant, P, of the gas in the containment at time c Representing the pressure of the gas at the outlet/in-containment of the gas-filled line, at representing t i Time to t i+1 Time length of time, V c,i,j Represents the volume corresponding to the jth block, H c,i+1,j Denotes the jth block at t i+1 Instantaneous relative humidity in the containment, H c,i,j Denotes the jth block at t i Relative humidity in Containment at time, f (T) c,i+1,j ) Denotes the jth block at t i+1 Saturated water vapor partial pressure at time f (T) c,i,j ) Denotes the jth block at t i+1 The saturated steam partial pressure in the containment vessel at the moment;
the a blocks are at t i To t i+1 Mass leakage rate of time period G out,i+1j The calculation formula of (c) is:
Figure RE-FDA0003865889630000082
wherein j represents the jth block, i represents the tth block i Time of day or t i To t i+1 Time period, a denotes the number of blocks, L in,i+1 Denotes the t-th i+1 Inflation volume flow, T, at the outlet of the line constantly inflated with gas c,i+1,j Denotes the jth block at t i+1 Absolute temperature of gas in containment vessel, T, at time of day c,i,j Denotes the jth block at t i Absolute temperature of gas in containment vessel at time m c,i,j Denotes the jth block at t i Mass of gas in containment vessel at time of day, R g,eq,i,j Denotes the jth block at t i Reduced gas constant, P, of gas in containment vessel at time of day c Indicating chargerPressure of gas at outlet of gas-in line/in containment, Δ t representing t i To t i+1 Time length of time, V c,i,j Represents the volume corresponding to the jth block, H c,i+1,j Denotes the jth block at t i+1 Relative humidity in containment vessel at time H c,i,j Denotes the jth block at t i Relative humidity in Containment at time, f (T) c,i+1,j ) Denotes the jth block at t i+1 Saturated vapor partial pressure in containment at time f (T) c,i,j ) Denotes the jth block at t i+1 Saturated vapor partial pressure, R, in the containment at that time g,eq,i+1,j Denotes the jth block at t i+1 The reduced gas constant of the gas in the containment at the moment.
CN202210508364.5A 2022-05-10 2022-05-10 Method for measuring leakage rate of containment Pending CN115240881A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210508364.5A CN115240881A (en) 2022-05-10 2022-05-10 Method for measuring leakage rate of containment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210508364.5A CN115240881A (en) 2022-05-10 2022-05-10 Method for measuring leakage rate of containment

Publications (1)

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

Family

ID=83668111

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210508364.5A Pending CN115240881A (en) 2022-05-10 2022-05-10 Method for measuring leakage rate of containment

Country Status (1)

Country Link
CN (1) CN115240881A (en)

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 (en) Procedure and measurement device for the determination of specific variables for the gas constitution
ES2638793T3 (en) Methods and systems for leak testing
US4941345A (en) Method and apparatus for the measurement of gas properties
CN115240881A (en) Method for measuring leakage rate of containment
US20230168171A1 (en) High throughput high resolution gas sorption screening
CN105203267B (en) A kind of space squeezes out normal pressure Integrating dosemeters System and method for
CN115938623A (en) Test method and system for containment leakage rate measurement, electronic equipment and medium
CN203534928U (en) Experimental device for testing adsorption/desorption performance of adsorbing working substance pair
CN106289372A (en) Welded insulated gas cylinder Daily boil-off-rate measuring instrument method of testing
KR102661248B1 (en) Methods for Evaluating Steam Pump Performance
US4916630A (en) Bed tester for molecular sieve oxygen concentrator
CN113933213A (en) Binary mixed gas mixing ratio measuring method and device based on gas substitution method
CN114864115A (en) Containment leakage rate measurement test system
van der Beek et al. Gas oil piston prover, primary reference values for gas-volume
RU2298774C1 (en) Method for controlling reservoir sealing tightness
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 (en) INSTALLATION FOR VERIFICATION OF INDUSTRIAL GAS METERS
RU138370U1 (en) DEVICE FOR TESTING PRODUCTS DURING QUICK PRESSURE CHANGE
EA029297B1 (en) Device for testing tightness of polymer fuel pipelines in vehicles
RU2246102C1 (en) Method of inspecting leak-proofness of spacecraft temperature control hydraulic system provided with hydropneumatic compensator and supplied with heat-transfer agent
RU2527659C1 (en) Method of tightness control
RU2710006C2 (en) Method for testing leakproof container
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