Method and system for calculating a pressurized water reactor primary coolant leakage rate

By determining the initial moment data and real-time parameter data in the calculation of the primary coolant leakage rate of pressurized water reactors, and combining the moving average method and the leakage rate calculation formula, a leakage rate calculation system was constructed. This system solved the problem of inaccurate calculation of the primary coolant leakage rate of pressurized water reactors, and achieved real-time monitoring and accurate leakage rate calculation.

CN116344081BActive Publication Date: 2026-07-21YANGJIANG NUCLEAR POWER +2
View PDF 1 Cites 0 Cited by

Patent Information

Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
YANGJIANG NUCLEAR POWER
Filing Date
2023-03-07
Publication Date
2026-07-21

Smart Images

  • Figure CN116344081B_ABST
    Figure CN116344081B_ABST
Patent Text Reader

Abstract

The present application relates to a pressurized water reactor primary coolant leakage rate calculation method and system, comprising the following steps: determining initial time data; obtaining real-time parameter data of the pressurized water reactor primary circuit; calculating the leakage rate according to the initial time data and the real-time parameter data and combining the leakage rate calculation formula to obtain the pressurized water reactor primary coolant leakage rate. Through automatic acquisition of real-time parameter data of the primary circuit and corresponding leakage rate calculation, compared with artificial selection of time points for calculation, the present application can effectively improve the accuracy of primary circuit leakage rate calculation, greatly reduce the daily workload of operators, and through real-time calculation of the primary circuit leakage rate, real-time monitoring of the primary circuit leakage rate can be realized, and changes in the primary circuit leakage rate can be found in time.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This invention relates to the technical field of pressurized water reactor primary loop leakage rate calculation, and more specifically, to a method and system for calculating pressurized water reactor primary loop coolant leakage rate. Background Technology

[0002] Periodic leak rate tests of the pressurized water reactor primary coolant circuit involve monitoring the amount of coolant leakage through precise measurement methods. This test is performed daily during normal unit operation, with each test lasting more than two hours. During the test, the main control operator needs to record the start and end times of the test and transcribe various parameters from the digital instrumentation and control system (DCS) screen for leak rate calculation. However, manually selecting time points may lead to deviations in the calculation results due to fluctuations in unit parameters, and real-time monitoring of the leak rate is not possible during the test. Summary of the Invention

[0003] The technical problem to be solved by the present invention is to provide a method and system for calculating the leakage rate of the primary coolant in a pressurized water reactor.

[0004] The technical solution adopted by this invention to solve its technical problem is: to construct a method for calculating the leakage rate of the primary coolant in a pressurized water reactor, including the following steps:

[0005] Determine the initial time data;

[0006] Obtain real-time parameter data for the primary loop of the pressurized water reactor;

[0007] Based on the initial time data and the real-time parameter data, and combined with the leakage rate calculation formula, the leakage rate is calculated to obtain the primary coolant leakage rate of the pressurized water reactor.

[0008] In the method for calculating the leakage rate of the primary coolant in a pressurized water reactor as described in this invention, the initial time data includes: the average temperature of the primary circuit at the initial time, the liquid level in the control tank at the initial time, the liquid level in the pressurizer at the initial time, the time value at the initial time, the water volume in the TEP head tank at the initial time, the dilution amount at the initial time, and the borylation amount at the initial time.

[0009] The real-time parameter data includes: the real-time temperature of the primary circuit, the real-time liquid level of the capacity control box, the real-time liquid level of the voltage regulator, the real-time time value, the real-time water volume of the TEP head box, the real-time dilution amount, and the real-time boration amount.

[0010] In the pressurized water reactor primary loop coolant leakage rate calculation method of the present invention, the pressurized water reactor primary loop coolant leakage rate includes: the first state leakage rate of coolant when the pressurized water reactor primary loop system is in a first state and the second state leakage rate of coolant when the pressurized water reactor primary loop system is in a second state.

[0011] The first state leakage rate includes: the leakage rate under full-power steady-state conditions in the first state; the leakage rate calculation formula includes: the first leakage rate calculation formula;

[0012] The step of calculating the leakage rate based on the initial time data and the real-time parameter data, combined with the leakage rate calculation formula, to obtain the primary coolant leakage rate of the pressurized water reactor includes:

[0013] The leakage rate under full-power stable conditions is calculated based on the average temperature of the primary circuit at the initial moment, the liquid level of the control tank at the initial moment, the liquid level of the voltage regulator at the initial moment, the time value at the initial moment, the real-time temperature of the primary circuit, the real-time liquid level of the control tank, the real-time liquid level of the voltage regulator, and the time value at the real-time moment, combined with the first leakage rate calculation formula; wherein, the real-time temperature of the primary circuit, the real-time liquid level of the control tank, and the real-time liquid level of the voltage regulator are values ​​under full-power stable conditions in the first state.

[0014] In the pressurized water reactor primary loop coolant leakage rate calculation method of the present invention, the first state leakage rate further includes: leakage rate under any state; the leakage rate calculation formula includes: second leakage rate calculation formula; the arbitrary state is all states of the pressurized water reactor primary loop except the full power stable state in the first state;

[0015] The step of calculating the leakage rate based on the initial time data and the real-time parameter data, combined with the leakage rate calculation formula, to obtain the primary coolant leakage rate of the pressurized water reactor includes:

[0016] The density of the primary loop subcooled water under the initial state, the density of the primary loop subcooled water under the real-time state, the density of the pressurizer saturated water under the initial state, and the density of the pressurizer saturated steam under the real-time state are calculated using the water vapor graph function.

[0017] The leakage rate under any given state is calculated based on the initial liquid level in the control tank, the initial liquid level in the pressure regulator, the initial time value, the real-time liquid level in the control tank, the real-time liquid level in the pressure regulator, the real-time time value, the density of the primary loop subcooled water in the initial state, the density of the primary loop subcooled water in the real-time state, the saturated water density of the pressure regulator in the initial state, and the saturated steam density of the pressure regulator in the real-time state, combined with the second leakage rate calculation formula.

[0018] In the pressurized water reactor primary coolant leakage rate calculation method of the present invention, the second state leakage rate includes: the leakage rate under full power stable state in the second state; the leakage rate calculation formula includes: the third leakage rate calculation formula;

[0019] The step of calculating the leakage rate based on the initial time data and the real-time parameter data, combined with the leakage rate calculation formula, to obtain the primary coolant leakage rate of the pressurized water reactor includes:

[0020] The leakage rate under full-power stable conditions is calculated based on the initial average temperature of the primary circuit, the initial liquid level of the control tank, the initial liquid level of the voltage regulator, the initial time value, the initial water volume of the TEP head tank, the initial dilution amount, the initial boronization amount, the real-time temperature of the primary circuit, the real-time liquid level of the control tank, the real-time liquid level of the voltage regulator, the real-time time value, the real-time water volume of the TEP head tank, the real-time dilution amount, and the real-time boronization amount, and in conjunction with the third leakage rate calculation formula. The real-time temperature of the primary circuit, the real-time liquid level of the control tank, the real-time liquid level of the voltage regulator, the real-time water volume of the TEP head tank, the real-time dilution amount, and the real-time boronization amount are values ​​under the second state of full-power stable conditions.

[0021] In the pressurized water reactor primary coolant leakage rate calculation method of the present invention, the second state leakage rate includes: leakage rate in any state under the second state; the leakage rate calculation formula includes: fourth leakage rate calculation formula;

[0022] The step of calculating the leakage rate based on the initial time data and the real-time parameter data, combined with the leakage rate calculation formula, to obtain the primary coolant leakage rate of the pressurized water reactor includes:

[0023] The density of the primary loop subcooled water under the initial state, the density of the primary loop subcooled water under the real-time state, the density of the pressurizer saturated water under the initial state, and the density of the pressurizer saturated steam under the real-time state are calculated using the water vapor graph function.

[0024] The leakage rate under any given state is calculated based on the initial liquid level in the control tank, the initial liquid level in the pressure regulator, the initial time value, the initial water volume in the TEP head tank, the initial dilution amount, the initial boronization amount, the real-time liquid level in the control tank, the real-time liquid level in the pressure regulator, the real-time time value, the real-time water volume in the TEP head tank, the real-time dilution amount, the real-time boronization amount, the density of the primary loop subcooled water in the initial state, the density of the primary loop subcooled water in the real-time state, the saturated water density of the pressure regulator in the initial state, and the saturated steam density of the pressure regulator in the real-time state, combined with the fourth leakage rate calculation formula.

[0025] In the pressurized water reactor primary coolant leakage rate calculation method of the present invention, the step of calculating the leakage rate based on the initial time data and the real-time parameter data and in combination with the leakage rate calculation formula to obtain the pressurized water reactor primary coolant leakage rate includes:

[0026] The real-time parameter data are averaged using the moving average method.

[0027] In the method for calculating the primary coolant leakage rate of a pressurized water reactor according to the present invention, the method includes:

[0028] Based on the leakage rate of the primary coolant in the pressurized water reactor, a curve fitting was performed to obtain a trend diagram of the leakage rate of the primary coolant in the pressurized water reactor.

[0029] The trend graph of the primary coolant leakage rate of the pressurized water reactor and the real-time data of the primary coolant leakage rate of the pressurized water reactor are displayed.

[0030] The present invention also provides a system for calculating the primary coolant leakage rate of a pressurized water reactor, comprising:

[0031] The initial value determination unit is used to determine the initial time data;

[0032] The real-time data acquisition unit is used to acquire real-time parameter data of the primary loop of the pressurized water reactor.

[0033] The leakage rate calculation unit is used to calculate the leakage rate based on the initial time data and the real-time parameter data, combined with the leakage rate calculation formula, to obtain the leakage rate of the pressurized water reactor primary coolant.

[0034] The pressurized water reactor primary coolant leakage rate calculation system of the present invention further includes:

[0035] The human-computer interaction unit is used to display the initial time data, the real-time parameter data of the pressurized water reactor primary loop, the real-time data of the pressurized water reactor primary loop coolant leakage rate, and the trend graph of the pressurized water reactor primary loop coolant leakage rate.

[0036] The method and system for calculating the primary coolant leakage rate of a pressurized water reactor (PWR) according to the present invention have the following beneficial effects: The method includes the following steps: determining initial time data; acquiring real-time parameter data of the PWR primary circuit; calculating the leakage rate based on the initial time data and real-time parameter data, combined with the leakage rate calculation formula, to obtain the PWR primary coolant leakage rate. By automatically acquiring real-time parameter data of the primary circuit and performing corresponding leakage rate calculations, compared to manually selecting time points for calculation, the present invention can effectively improve the accuracy of primary circuit leakage rate calculation, significantly reduce the daily workload of operators, and, through real-time calculation of the primary circuit leakage rate, achieve real-time monitoring of the primary circuit leakage rate, enabling timely detection of changes in the primary circuit leakage rate. Attached Figure Description

[0037] The present invention will be further described below with reference to the accompanying drawings and embodiments. In the accompanying drawings:

[0038] Figure 1 This is a flowchart illustrating the method for calculating the primary coolant leakage rate of a pressurized water reactor provided in an embodiment of the present invention.

[0039] Figure 2 This is a schematic diagram of the logic configuration for calculating the primary coolant leakage rate of a pressurized water reactor provided in an embodiment of the present invention;

[0040] Figure 3 This is a schematic diagram of the pressurized water reactor primary coolant leakage rate calculation system provided in the embodiments of the invention;

[0041] Figure 4 This is a diagram of the human-computer interaction unit of the pressurized water reactor primary loop coolant leakage rate calculation system provided in this embodiment of the invention. Detailed Implementation

[0042] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0043] To address the problems of inaccurate calculation of primary coolant leakage rate and high human error rate in current pressurized water reactor (PWR) systems, this invention provides a method for calculating the primary coolant leakage rate of PWRs. This method is based on a DCS (Digital Control System) platform for leakage rate calculation.

[0044] In a preferred embodiment, such as Figure 1 As shown, the method for calculating the primary coolant leakage rate of a pressurized water reactor provided by this invention includes the following steps:

[0045] Step S10: Determine the initial time data.

[0046] In this embodiment of the invention, the initial time data includes: the initial average temperature of the primary loop, the initial liquid level in the control tank, the initial liquid level in the pressurizer, the initial time value, the initial water volume in the TEP head tank, the initial dilution amount, and the initial boronization amount. The initial time is when the test begins, and is directly acquired through the DCS platform based on the user's trigger operation, monitoring the primary loop of the pressurized water reactor.

[0047] Step S20: Obtain real-time parameter data of the primary loop of the pressurized water reactor.

[0048] In this embodiment of the invention, the real-time parameter data includes: the real-time temperature of the primary circuit, the real-time liquid level of the capacity control box, the real-time liquid level of the voltage regulator, the real-time time value, the real-time water volume of the TEP head box, the real-time dilution amount, and the real-time boration amount.

[0049] In this embodiment of the invention, the real-time parameter data of the pressurized water reactor primary loop can be obtained through a DCS platform. Specifically, relevant parameters can be acquired using temperature sensors, level sensors, and other sensors within the DCS platform.

[0050] Furthermore, in this embodiment of the invention, before calculating the leakage rate based on the initial time data and real-time parameter data and in conjunction with the leakage rate calculation formula, the following steps are taken: the real-time parameter data is averaged using the moving average method.

[0051] Specifically, to eliminate the impact of fluctuations in real-time parameter data, this invention employs a moving average method for smoothing the input real-time parameter data. This involves averaging the input real-time parameter data over a rolling period. The resulting value represents the median parameter value over a given period. The moving average method is characterized by its speed and variability; therefore, using it to smooth the data effectively eliminates the impact of parameter fluctuations. For example, in this embodiment, the real-time parameter data required for calculation is collected every second. After 60 collections (60 seconds), the average is calculated. This collection and averaging process is continuous; each time data is collected, the newly collected data is summed to calculate the average, thus achieving the smoothing effect of the moving average.

[0052] The specific code for the moving average processing is as follows:

[0053]

[0054]

[0055] The above algorithms will be encapsulated in DCS and can be called as needed for parameter fitting.

[0056] Step S30: Calculate the leakage rate based on the initial time data and real-time parameter data, combined with the leakage rate calculation formula, to obtain the primary coolant leakage rate of the pressurized water reactor.

[0057] In this embodiment of the invention, the pressurized water reactor primary loop system can generally be divided into two states: a first state and a second state. Specifically, the first state is a state where the boron concentration in the pressurized water reactor primary loop is stable. Within this stable boron concentration state, there are three operating conditions: RP full-power mode (i.e., operating at full power and stable state), RP non-full-power mode (i.e., not operating at full power), and non-RP mode (i.e., non-operating state).

[0058] The second state is the state where the boron concentration in the primary loop of the pressurized water reactor is unstable. Similarly, when the boron concentration in the primary loop of the pressurized water reactor is unstable, there are three operating conditions: RP full power mode (i.e., stable operation at full power), RP non-full power mode (i.e., non-operational full power state), and non-RP mode (i.e., non-operational state).

[0059] Therefore, in this embodiment of the invention, the coolant leakage rate of the pressurized water reactor primary loop includes: the first-state leakage rate of the coolant when the pressurized water reactor primary loop system is in a first state and the second-state leakage rate of the coolant when the pressurized water reactor primary loop system is in a second state.

[0060] Optionally, in this embodiment of the invention, the leakage rate in the first state includes: the leakage rate in the full-power stable state of the first state (RP full-power leakage rate) and the leakage rate in any state (RP non-full-power leakage rate and non-RP leakage rate). Wherein, any state refers to all states of the pressurized water reactor primary loop other than the full-power stable state in the first state. The leakage rate in the second state includes: the leakage rate in the full-power stable state of the second state (RP full-power leakage rate) and the leakage rate in any state of the second state (RP non-full-power leakage rate and non-RP leakage rate).

[0061] Optionally, in this embodiment of the invention, the leakage rate calculation formula includes: a first leakage rate calculation formula, a second leakage rate calculation formula, a third leakage rate calculation formula, and a fourth leakage rate calculation formula.

[0062] Specifically, when the pressurized water reactor primary loop system is in the first state, the first-state leakage rate can be calculated. The leakage rate under full-power steady-state conditions within the first-state leakage rate can be calculated using the following method:

[0063] The leakage rate is calculated based on initial time data and real-time parameter data, combined with the leakage rate calculation formula. The leakage rate of the pressurized water reactor primary coolant is obtained by: using the initial time average temperature of the primary loop, the initial time level of the control tank, the initial time level of the pressurizer, the initial time value, the real-time temperature of the primary loop, the real-time level of the control tank, the real-time level of the pressurizer, and the real-time time value, combined with the first leakage rate calculation formula, to obtain the leakage rate under full-power stable conditions; where the real-time temperature of the primary loop, the real-time level of the control tank, and the real-time level of the pressurizer are values ​​under the first state of full-power stable conditions.

[0064] Specifically, the leakage rate under the full-power stable state in the first state can be calculated based on the initial average temperature of the primary loop, the initial liquid level of the control tank, the initial liquid level of the voltage regulator, the initial time value, the real-time temperature of the primary loop, the real-time liquid level of the control tank, the real-time liquid level of the voltage regulator, and the real-time time value, combined with the first leakage rate calculation formula.

[0065] Optionally, in this embodiment of the invention, the formula for calculating the first leakage rate satisfies:

[0066]

[0067] Among them, F P1 The leakage rate is given by T1, the average temperature of the primary loop at the initial moment, H1, the liquid level in the control tank at the initial moment, L1, the liquid level in the voltage regulator at the initial moment, and t1, the time value at the initial moment. T2 represents the real-time temperature of the primary loop at full power stability, H2 represents the real-time liquid level in the control tank at full power stability, L2 represents the real-time liquid level in the voltage regulator at full power stability, and t2 represents the time value at the termination moment. Δt = t2 - t1. It should be noted that t1 and t2 are user-defined times. When the test begins, the user triggers the "t1 SET" button, and the system directly obtains the relevant parameters at time t1 (i.e., H1, T1, L1, P1, etc.). When the user decides to end the test, they trigger the "t2 SET" button, and the system directly obtains the relevant parameters at time t2 (i.e., H2, T2, L2, P2, etc.).

[0068] The leakage rate of RP (non-full power) and the leakage rate of non-RP in the first state leakage rate can be calculated using the following method:

[0069] The leakage rate is calculated based on initial and real-time data and the leakage rate calculation formula. This includes: calculating the primary coolant leakage rate of the pressurized water reactor using a water vapor graph function; calculating the primary coolant subcooling water density under initial conditions, the primary coolant subcooling water density under real-time conditions, the pressurizer saturated water density under initial conditions, and the pressurizer saturated steam density under real-time conditions; and calculating the leakage rate under any condition using the initial tank level, the pressurizer level, the initial time value, the real-time tank level, the pressurizer level, the real-time time value, the primary coolant subcooling water density under initial conditions, the primary coolant subcooling water density under real-time conditions, the pressurizer saturated water density under initial conditions, and the pressurizer saturated steam density under real-time conditions, combined with the second leakage rate calculation formula.

[0070] Optionally, in this embodiment of the invention, the formula for calculating the second leakage rate satisfies:

[0071]

[0072] Among them, F P2 Let H1 be the leakage rate under any given state, H1 be the initial liquid level in the control tank, L1 be the initial liquid level in the pressure regulator, and t1 be the initial time value; H2 be the real-time liquid level in the control tank under any given state, L2 be the real-time liquid level in the pressure regulator under any given state, t2 be the time value at the termination time, and Δt = t2 - t1; ρ L The density of saturated water in the voltage regulator (kg / m3) under state P1 / T1 (i.e., initial state) is ρ. G The saturated steam density of the pressurizer (kg / m3) under the P2 / T2 state (i.e., real-time state), ρ1 is the density of the primary loop subcooled water (kg / m3) under the P1 / T1 state, and ρ2 is the density of the primary loop subcooled water (kg / m3) under the P2 / T2 state.

[0073] Specifically, when the pressurized water reactor primary loop system is in the second state, the second-state leakage rate can be calculated. The leakage rate of RP under full-power steady-state conditions within the second-state leakage rate can be calculated using the following method:

[0074] The leakage rate is calculated based on initial time data and real-time parameter data, combined with the leakage rate calculation formula. The leakage rate of the pressurized water reactor primary loop coolant is obtained by: using the initial time data of the primary loop average temperature, initial time data of the control tank level, initial time data of the pressurizer level, initial time value, initial time data of the TEP head tank water volume, initial time data of the dilution amount, initial time data of the boron amount, the real-time temperature of the primary loop, the real-time level of the control tank, the real-time level of the pressurizer, the real-time time value, the real-time water volume of the TEP head tank, the real-time dilution amount, and the real-time boron amount, combined with the third leakage rate calculation formula, to obtain the leakage rate under full power stable conditions. Among them, the real-time temperature of the primary loop, the real-time level of the control tank, the real-time level of the pressurizer, the real-time water volume of the TEP head tank, the real-time dilution amount, and the real-time boron amount are the values ​​under the full power stable conditions in the second state.

[0075] Optionally, in this embodiment of the invention, the third leakage rate calculation formula satisfies:

[0076]

[0077]

[0078] (3) In the formula, F P3 The leakage rate is given by H1, the initial liquid level of the control tank, L1, the initial liquid level of the regulator, and t1, the initial time value. H2 is the real-time liquid level of the control tank under full power stability, L2 is the real-time liquid level of the regulator under full power stability, t2 is the time value at the termination time, M1 is the initial water volume of the TEP head tank, M2 is the real-time water volume of the TEP head tank, W1 is the initial dilution amount, W2 is the real-time dilution amount, B1 is the initial boronization amount, and B2 is the real-time boronization amount.

[0079] The leakage rate of RP (non-full power) and the leakage rate of non-RP in the second state leakage rate can be calculated using the following method:

[0080] The leakage rate is calculated based on initial and real-time data and the leakage rate calculation formula. This includes: calculating the primary coolant leakage rate of the pressurized water reactor using a water vapor graph function; calculating the primary coolant subcooling water density under initial conditions, the primary coolant subcooling water density under real-time conditions, the pressurizer saturated water density under initial conditions, and the pressurizer saturated steam density under real-time conditions; and calculating the leakage rate under any condition using the initial values ​​of the control tank level, pressurizer level, time, TEP head tank water volume, dilution amount, boronization amount, real-time control tank level, pressurizer level, time, TEP head tank water volume, dilution amount, boronization amount, primary coolant subcooling water density under initial conditions, primary coolant subcooling water density under real-time conditions, pressurizer saturated water density under initial conditions, and pressurizer saturated steam density under real-time conditions, combined with the fourth leakage rate calculation formula.

[0081] Optionally, in this embodiment of the invention, the fourth leakage rate calculation formula satisfies:

[0082]

[0083]

[0084] (4) In the formula, F P4 Let H1 be the leakage rate under any given state, H1 be the initial liquid level in the control tank, L1 be the initial liquid level in the pressure regulator, and t1 be the initial time value; H2 be the real-time liquid level in the control tank under any given state, L2 be the real-time liquid level in the pressure regulator under any given state, and t2 be the time value at the termination time; ρ L The density of saturated water in the voltage regulator (kg / m3) under state P1 / T1 (i.e., initial state) is ρ. G The saturated steam density of the pressurizer (kg / m3) under P2 / T2 state (i.e., real-time state), ρ1 is the subcooled water density of the first loop under P1 / T1 state (kg / m3), ρ2 is the subcooled water density of the first loop under P2 / T2 state (kg / m3), M1 is the water volume of the TEP head tank at the initial moment, M2 is the real-time water volume of the TEP head tank, W1 is the dilution amount at the initial moment, W2 is the real-time dilution amount, B1 is the boronization amount at the initial moment, and B2 is the real-time boronization amount.

[0085] In this embodiment of the invention, the software codes for the three water densities—saturated water density, saturated vapor density, and subcooled water density—are independent of each other, enabling accurate calculation of water vapor density at pressures ranging from 1 to 165 bar. This allows for primary loop leakage rate calculations under full-power RP, non-full-power RP, and non-RP modes (states). The specific code is as follows:

[0086] Specifically, the algorithm for calculating the density of subcooled water is as follows:

[0087] FUNCTION_BLOCK V1_PTGL1_CAL

[0088] VAR_INPUT

[0089] P:REAL;

[0090] T:REAL;

[0091] END_VAR

[0092] VAR_OUTPUT

[0093] V1_PTGL1:REAL;

[0094] END_VAR

[0095] VAR

[0096] I1:ARRAY[1..34]OF REAL:=0,0,0,0,0,0,0,0,1,1,1,1,1,1,2,2,2,2,2,3,3,3,4,4,4,5,8,8,21,23,29,30,31,32;

[0097] J1:ARRAY[1..34]OF REAL:=-2,-1,0,1,2,3,4,5,-9,-7,-1,0,1,3,-3,0,1,3,17,-4,0,6,-5,-2,10,-8,-11,-6,-29,-31,-38,-39,-40,-41;

[0098] n1:ARRAY[1..34]OF REAL:=0.14632971213167,-0.84548187169114,-3.756360367204,3.3855169168385,-0.9579196 3387872,0.15772038513228,-0.016616417199501,8.1214629983568E-04,2.8319080123804E-04,

[0099] -6.0706301565874E-04,-0.018990068218419,-0.032529748770505,

[0100] -0.021841717175414,-5.283835796993E-05,-4.7184321073267E-04,

[0101] -3.0001780793026E-04,4.7661393906987E-05,-4.4141845330846E-06,

[0102] -7.2694996297594E-16,-3.1679644845054E-05,-2.8270797985312E-06,

[0103] -8.5205128120103E-10,-2.2425281908E-06,-6.5171222895601E-07,

[0104] -1.4341729937924E-13,-4.0516996860117E-07,-1.2734301741641E-09,

[0105] -1.7424871230634E-10,-6.8762131295531E-19,1.4478307828521E-20,

[0106] 2.6335781662795E-23,

[0107] -1.1947622640071E-23,1.8228094581404E-24,-9.3537087292458E-26;

[0108] TAU:REAL;

[0109] R1:REAL;

[0110] G_P:REAL;

[0111] G_P1:REAL;

[0112] PS:REAL;

[0113] i:INT;

[0114] G_P2:REAL;

[0115] V1_PTG:REAL;

[0116] U:INT;

[0117] END_VAR

[0118] IF MIDUCAL_START=TRUE THEN IF I>=20THEN

[0119] I:=0;

[0120] R1:=0.461526;

[0121] PS: = P / 16.53;

[0122] TAU: = 1386 / T;

[0123] G_P:=0;

[0124] FOR U:=1TO 34BY+1DO G_P1:=EXPT((7.1-PS),(I1[U]-1)); G_P2:=EXPT((TAU-1.222),J1[U]); G_P:=G_P-n1[U]*I1[U]*G_P1*G_P2; END_FOR;

[0125] V1_PTG:=R1*T / P*PS*G_P / 1000;

[0126] V1_PTGL1:=1 / V1_PTG;

[0127] END_IFI:=I+1;

[0128] The algorithm for calculating the saturated water density using the END_IFMIDUCAL_START instruction is as follows: FUNCTION_BLOCKV1_PTBS_CAL VAR_INPUT

[0129] P:REAL;

[0130] END_VAR

[0131] VAR_OUTPUT

[0132] V1_PTBS:REAL;

[0133] END_VAR

[0134] VAR

[0135] BETA:REAL;

[0136] F:REAL;

[0137] H:REAL;

[0138] E:REAL;

[0139] G:REAL;

[0140] K:REAL;

[0141] D:REAL;

[0142] L:REAL;

[0143] J:REAL;

[0144] T4_P:REAL;

[0145] T:REAL;

[0146] R1:REAL;

[0147] PS:REAL;

[0148] G_P1:REAL;

[0149] U:INT;

[0150] G_P:REAL;

[0151] G_P2:REAL;

[0152] TAU:REAL;

[0153] I1:ARRAY[1..34]OF REAL:=0,0,0,0,0,0,0,0,1,1,1,1,1,1,2,2,2,2,2,3,3,3,4,4,4,5,8,8,21,23,29,30,31,32;

[0154] J1:ARRAY[1..34]OF REAL:=-2,-1,0,1,2,3,4,5,-9,-7,-1,0,1,3,-3,0,1,3,17,-4,0,6,-5,-2,10,-8,-11,-6,-29,-31,-38,-39,-40,-41;

[0155] n1:ARRAY[1..34]OF REAL:=0.14632971213167,-0.84548187169114,-3.756360367204,3.3855169168385,-0.95791963387872,0.15772038513228,-0.016616417199501,8.1214629983568E-04,2.8319080123804E-04,

[0156] -6.0706301565874E-04,-0.018990068218419,-0.032529748770505,

[0157] -0.021841717175414,-5.283835796993E-05,-4.7184321073267E-04,

[0158] -3.0001780793026E-04,4.7661393906987E-05,-4.4141845330846E-06,

[0159] -7.2694996297594E-16,-3.1679644845054E-05,-2.8270797985312E-06,

[0160] -8.5205128120103E-10,-2.2425281908E-06,-6.5171222895601E-07,

[0161] -1.4341729937924E-13,-4.0516996860117E-07,-1.2734301741641E-09,

[0162] -1.7424871230634E-10,-6.8762131295531E-19,1.4478307828521E-20,

[0163] 2.6335781662795E-23,

[0164] -1.1947622640071E-23,1.8228094581404E-24,-9.3537087292458E-26;

[0165] V1_pt:REAL;

[0166] I:INT;

[0167] END_VAR

[0168] IF MIDUCAL_START=TRUE THEN

[0169] IF I>=20 THEN

[0170] I:=0;

[0171] BETA:=EXPT(P,0.25);

[0172] E:=EXPT(BETA,2)-17.073846940092*BETA+14.91510861353;

[0173] F:=1167.0521452767

[0174] *EXPT(BETA,2)+12020.82470247*BETA-4823.2657361591;

[0175] G:=-724213.16703206*EXPT(BETA,2)-3232555.0322333*BETA+405113.40542057;

[0176] H:=EXPT(F,2)-4*E*G;

[0177] J:=-F-EXPT(H,0.5);

[0178] D:=2*G / J;

[0179] K:=EXPT((650.17534844798+D),2);

[0180] L:=EXPT((K-4*(-0.23855557567849+650.17534844798*D)),0.5);

[0181] T4_P:=(650.17534844798+D-L) / 2;

[0182] T:=T4_P;

[0183] R1:=0.461526;

[0184] PS:=P / 16.53;

[0185] TAU:=1386 / T;

[0186] G_P:=0;

[0187] FOR U:=1 TO 34 BY+1 DO

[0188] G_P1:=EXPT((7.1-PS),(I1[U]-1));

[0189] G_P2:=EXPT((TAU-1.222),J1[U]);

[0190] G_P:=G_P-n1[U]*I1[U]*G_P1*G_P2;

[0191] END_FOR;

[0192] V1_pt:=R1*T / P*PS*G_P / 1000;

[0193] V1_PTBS:=1 / V1_pt;

[0194] END_IF

[0195] I:=I+1;

[0196] END_IF

[0197] MIDUCAL_START is the instruction to start the density algorithm.

[0198] The algorithm for saturated vapor density is as follows:

[0199] FUNCTION_BLOCK V2_PTBZ_CAL

[0200] VAR_INPUT

[0201] P:REAL;

[0202] END_VAR

[0203] VAR_OUTPUT

[0204] V2_PTBZ:REAL;

[0205] END_VAR

[0206] VAR

[0207] Ir:ARRAY[1..43]OF REAL:=1,1,1,1,1,2,2,2,2,2,3,3,3,3,3,4,4,4,5,6,6,6,7,7,7,8,8,9,10,10,10,16,16,18,20,20,20,21,22,23,24,24,24;

[0208] Jr:ARRAY[1..43]OF REAL:=0,1,2,3,6,1,2,4,7,36,0,1,3,6,35,1,2,3,7,3,16,35,0,11,25,8,36,13,4,10,14,29,50,57,20,35,48,21,53,39,26,40,58;

[0209] nr:ARRAY[1..43]OFREAL:=-1.7731742473213E-03,

[0210] -0.017834862292358,-0.045996013696365,-0.057581259083432,

[0211] -0.05032527872793,-3.3032641670203E-05,-1.8948987516315E-04,-3.9392777243355E-03,-0.043797295650573,-2.6674547914087E-05,2.0481737692309E-08,4.3870667284435E-07,-3.227767723857E-05,-1.5033924542148E-03,-0.040668253562649,-7.8847309559367E-10,1.2790717852285E-08,4.8225372718507E-07,2.2922076337661E-06,-1.6714766451061E-11,-2.1171472321355E-03,-23.895741934104,

[0212] -5.905956432427E-18,-1.2621808899101E-06,-0.038946842435739,1.1256211360459E-11,-8.2311340897998,1.9809712802088E-08,1.0406965210174E-19,-1.0234747095929E-13,-1.0018179379511E-09,-8.0882908646985E-11,0.10693031879409,-0.33662250574171,8.9185845355421E-25,3.0629316876232E-13,-4.2002467698208E-06,-5.9056029685639E-26,3.7826947613457E-06,-1.2768608934681E-15,7.3087610595061E-29,5.5414715350778E-17,-9.436970724121E-07;

[0213] G:REAL;

[0214] BETA:REAL;

[0215] H:REAL;

[0216] F:REAL;

[0217] E:REAL;

[0218] P1:REAL;

[0219] T1:REAL;

[0220] T4_P:REAL;

[0221] i:INT;

[0222] GR_PI1:REAL;

[0223] GR_PI2:REAL;

[0224] TAU1:REAL;

[0225] GR_PI:REAL;

[0226] V2_PT:REAL;

[0227] R2:REAL;

[0228] G0_PI:REAL;

[0229] D:REAL;

[0230] J:REAL;

[0231] K:REAL;

[0232] L:REAL;

[0233] U:INT;

[0234] END_VAR

[0235] IF MIDUCAL_START=TRUE THEN

[0236] IF I>=20 THEN

[0237] I:=0;

[0238] BETA:=EXPT(P,0.25);

[0239] E:=EXPT(BETA,2)-17.073846940092*BETA+14.91510861353;

[0240] F:=1167.0521452767

[0241] *EXPT(BETA,2)+12020.82470247*BETA-4823.2657361591;

[0242] G:=-724213.16703206*EXPT(BETA,2)-3232555.0322333*BETA+405113.40542057;

[0243] H:=EXPT(F,2)-4*E*G;

[0244] J:=-F-EXPT(H,0.5);

[0245] D:=2*G / J;

[0246] K:=EXPT((650.17534844798+D),2);

[0247] L:=EXPT((K-4*(-0.23855557567849+650.17534844798*D)),0.5);T4_P:=(650.17534844798+D-L) / 2;

[0248] P1:=P;

[0249] T1:=T4_P;

[0250] R2:=0.461526;

[0251] TAU1:=540 / T1;

[0252] G0_PI:=1 / P1;

[0253] GR_PI:=0;

[0254] FOR U:=1TO 43BY+1DOGR_PI1:=EXPT(P1,(Ir[U]-1));

[0255] GR_PI2:=EXPT((TAU1-0.5),Jr[U]);

[0256] GR_PI:=GR_PI+nr[U]*Ir[U]*GR_PI1*GR_PI2;

[0257] END_FORV2_PT:=R2*T1 / P1*P1*(G0_PI+GR_PI) / 1000;

[0258] V2_PTBZ:=1 / V2_PT;

[0259] END_IF

[0260] I:=I+1;

[0261] END_IFMIDUCAL_START is the density input instruction. The above algorithm will be encapsulated in DCS and can be called as needed.

[0262] Furthermore, in this embodiment of the invention, the method for calculating the primary coolant leakage rate of a pressurized water reactor further includes: performing curve fitting based on the primary coolant leakage rate of the pressurized water reactor to obtain a trend graph of the primary coolant leakage rate; and displaying the trend graph of the primary coolant leakage rate and the real-time data of the primary coolant leakage rate.

[0263] Specifically, in this embodiment of the invention, when the real-time data of the pressurized water reactor primary coolant leakage rate is calculated, the calculated real-time data is also displayed in real time to achieve real-time monitoring of the pressurized water reactor primary coolant leakage rate.

[0264] Furthermore, in this embodiment of the invention, curve fitting can be performed based on the real-time data of the calculated pressurized water reactor primary coolant leakage rate to obtain a trend graph of the pressurized water reactor primary coolant leakage rate, and the trend graph of the pressurized water reactor primary coolant leakage rate can be displayed in real time so that the operator can intuitively and in real time grasp the changing trend of the pressurized water reactor primary coolant leakage rate and promptly detect abnormal changes.

[0265] Furthermore, in this embodiment of the invention, the method for calculating the primary coolant leakage rate of a pressurized water reactor further includes: receiving a user-input instruction to output the final calculation result; outputting the final value of the primary coolant leakage rate of the pressurized water reactor according to the instruction to output the final calculation result and printing it out.

[0266] refer to Figure 2 , Figure 2 This is a logic configuration diagram of the full-power RP leakage rate provided by the present invention.

[0267] in, Figure 2 In the file, PY34CV011MN_TEST1.AV, PY3RCP007MN_TEST1.AV, and PY3RCP940KM_TEST1.AV represent the initial time data;

[0268] PY3RCV011MN_NET01, LEAKAGE_003.SLIDAVE_201.OUT.CV, and LEAKAGE_002.SLIDAVE_201.OUT.CV represent real-time parameter data.

[0269] Furthermore, Figure 2In this logic configuration, SUB is the subtraction module, used to perform subtraction operations; MUL is the multiplication module, used to perform multiplication operations; DIV is the division module, used to perform division operations; ADD is the addition module, used to perform addition operations; and SEL is the selection module, used to perform corresponding selection operations according to user instructions. This logic configuration diagram calculates the RP full-power leakage rate according to formula (1). Specifically, when the pressurized water reactor primary loop is in the RP full-power state, the user triggers the SEL module. At this time, the SEL module outputs the leakage rate at the current moment, which is the RP full-power leakage rate.

[0270] In a preferred embodiment, such as Figure 3 As shown, the present invention provides a system for calculating the primary coolant leakage rate of a pressurized water reactor.

[0271] Specifically, such as Figure 3 As shown, the pressurized water reactor primary coolant leakage rate calculation system includes:

[0272] The initial value determination unit 301 is used to determine the initial time data.

[0273] The real-time data acquisition unit 302 is used to acquire real-time parameter data of the pressurized water reactor primary loop.

[0274] The leakage rate calculation unit 303 is used to calculate the leakage rate based on the initial time data and real-time parameter data, combined with the leakage rate calculation formula, to obtain the leakage rate of the pressurized water reactor primary coolant.

[0275] Furthermore, such as Figure 3 As shown, the pressurized water reactor primary loop coolant leakage rate calculation system also includes a human-machine interface unit 304. This human-machine interface unit 304 is used to display initial time data, real-time parameter data of the pressurized water reactor primary loop, real-time data of the pressurized water reactor primary loop coolant leakage rate, and a trend graph of the pressurized water reactor primary loop coolant leakage rate.

[0276] In this embodiment of the invention, the human-machine interface unit 304 of the pressurized water reactor primary coolant leakage rate calculation system realizes specific functions around the three dimensions of "low human-caused risk, few operation steps, and wide information coverage". It can display real-time data of primary coolant leakage rate. In specific application, it is only necessary to perform the corresponding operation on the interactive interface of the human-machine interface unit 304 (such as clicking start, end, and reset). At the same time, the operation on the interactive interface will not bring any risk to the unit.

[0277] like Figure 4 As shown, this invention provides a schematic diagram of a preferred embodiment of the interactive interface of the human-computer interaction unit 304. Specifically, as... Figure 4As shown, "t1 SET" is the start button (for acquiring initial data) and "t2 SET" is the end button. Specifically, when the user presses the "t1 SET" button, the experiment begins, and the system directly acquires the relevant parameters at time t1 through the DCS platform. When the user determines the end time, they press the "t2 SET" button, and the system directly acquires the relevant parameters at time t2 through the DCS platform. The currently acquired real-time data is the data at the end time. It should be noted that after "t1 SET" is pressed, the real-time parameter data related to the primary coolant leakage rate calculation is processed using the moving average method, and the leakage rate calculation is performed while maintaining the initial data. This achieves the real-time calculation function of the primary coolant leakage rate, which is displayed in real-time through the interface of the human-machine interface unit 304, and also displayed in the form of a curve (trend graph), thus achieving visualized monitoring of the primary coolant leakage rate.

[0278] The human-machine interface unit 304 displays the parameters required for calculating the primary loop leakage rate of the pressurized water reactor. It also records the test start time, duration, and stop time, and displays real-time leakage rate data and trend graphs. Figure 4 In the diagram, B.1 displays the relevant parameters and leakage rate in the first state; B.2 displays the relevant parameters and leakage rate in the second state.

[0279] Specifically, the specific coordination and operation process between the units in the pressurized water reactor primary coolant leakage rate calculation system can be referred to the above-mentioned pressurized water reactor primary coolant leakage rate calculation method, and will not be repeated here.

[0280] This invention uses the moving average method to smooth the calculation of relevant instrument parameters (real-time parameter data), which can avoid the influence of random sampling on the test results and can be accurate to the second, resulting in more accurate results. At the same time, by obtaining real unit data and relying on the leakage rate calculation formula, the leakage rate configuration of the pressurized water reactor primary loop is accurately calculated, which further proves that the leakage rate calculation accuracy of this invention is far higher than that of traditional manual sampling calculation.

[0281] Furthermore, this invention designs a single screen (the interactive interface of the human-computer interaction unit 304) that centralizes the parameters required for calculating the pressurized water reactor primary coolant leakage rate. Once the initial and final times (the final time being triggered by the user when the system is in a stable full-power state) are determined, the leakage rate under full-power operation can be calculated quickly and accurately. Simultaneously, this interactive interface also displays the real-time leakage rate of the pressurized water reactor primary coolant, presenting it as a curve. This allows users to monitor the leakage rate in real-time using the curve, making the monitoring of the pressurized water reactor primary coolant leakage rate visible.

[0282] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the apparatus disclosed in the embodiments, since it corresponds to the method disclosed in the embodiments, the description is relatively simple; relevant parts can be referred to the method section.

[0283] Those skilled in the art will further recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

[0284] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein can be implemented directly by hardware, a software module executed by a processor, or a combination of both. The software module can be located in random access memory (RAM), main memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art.

[0285] The above embodiments are only for illustrating the technical concept and features of the present invention, and are intended to enable those skilled in the art to understand the content of the present invention and implement it accordingly. They do not limit the scope of protection of the present invention. All equivalent changes and modifications made within the scope of the claims of the present invention should fall within the scope of the claims of the present invention.

Claims

1. A method for calculating the primary coolant leakage rate of a pressurized water reactor, characterized in that, Includes the following steps: Determine the initial time data; the initial time is when the test begins, and it is obtained directly from the primary loop of the pressurized water reactor through the DCS platform based on the user's trigger operation; Real-time parameter data of the primary loop of the pressurized water reactor is acquired through the DCS platform. Based on the initial time data and the real-time parameter data, and combined with the leakage rate calculation formula, the leakage rate is calculated to obtain the primary coolant leakage rate of the pressurized water reactor. The initial time data includes: the initial average temperature of the primary loop, the initial liquid level of the control tank, the initial liquid level of the voltage regulator, the initial time value, the initial water volume of the TEP head tank, the initial dilution amount, and the initial borylation amount. The real-time parameter data includes: the real-time temperature of the primary circuit, the real-time liquid level of the capacity control box, the real-time liquid level of the voltage regulator, the real-time time value, the real-time water volume of the TEP head box, the real-time dilution amount, and the real-time boronization amount. The pressurized water reactor primary loop coolant leakage rate includes: the first-state leakage rate of the coolant when the pressurized water reactor primary loop system is in a first state and the second-state leakage rate of the coolant when the pressurized water reactor primary loop system is in a second state; the first state is a state in which the boron concentration in the pressurized water reactor primary loop is stable, and the second state is a state in which the boron concentration in the pressurized water reactor primary loop is unstable. The leakage rate in the first state includes: the leakage rate in the full-power stable state under the first state, and the leakage rate in any state; the arbitrary state refers to all states of the pressurized water reactor primary loop other than the full-power stable state under the first state. The leakage rate in the second state includes: the leakage rate in the full-power steady state in the second state, and the leakage rate in any state in the second state; The leakage rate calculation formulas include: a first leakage rate calculation formula, a second leakage rate calculation formula, a third leakage rate calculation formula, and a fourth leakage rate calculation formula.

2. The method for calculating the primary coolant leakage rate of a pressurized water reactor according to claim 1, characterized in that, The step of calculating the leakage rate based on the initial time data and the real-time parameter data, combined with the leakage rate calculation formula, to obtain the primary coolant leakage rate of the pressurized water reactor includes: The leakage rate under full-power stable conditions is calculated based on the average temperature of the primary circuit at the initial moment, the liquid level of the control tank at the initial moment, the liquid level of the voltage regulator at the initial moment, the time value at the initial moment, the real-time temperature of the primary circuit, the real-time liquid level of the control tank, the real-time liquid level of the voltage regulator, and the time value at the real-time moment, combined with the first leakage rate calculation formula; wherein, the real-time temperature of the primary circuit, the real-time liquid level of the control tank, and the real-time liquid level of the voltage regulator are values ​​under full-power stable conditions in the first state.

3. The method for calculating the primary coolant leakage rate of a pressurized water reactor according to claim 1, characterized in that, The step of calculating the leakage rate based on the initial time data and the real-time parameter data, combined with the leakage rate calculation formula, to obtain the primary coolant leakage rate of the pressurized water reactor includes: The density of the primary loop subcooled water under the initial state, the density of the primary loop subcooled water under the real-time state, the density of the pressurizer saturated water under the initial state, and the density of the pressurizer saturated steam under the real-time state are calculated using the water vapor graph function. The leakage rate under any given state is calculated based on the initial liquid level in the control tank, the initial liquid level in the pressure regulator, the initial time value, the real-time liquid level in the control tank, the real-time liquid level in the pressure regulator, the real-time time value, the density of the primary loop subcooled water in the initial state, the density of the primary loop subcooled water in the real-time state, the saturated water density of the pressure regulator in the initial state, and the saturated steam density of the pressure regulator in the real-time state, combined with the second leakage rate calculation formula.

4. The method for calculating the primary coolant leakage rate of a pressurized water reactor according to claim 2, characterized in that, The step of calculating the leakage rate based on the initial time data and the real-time parameter data, combined with the leakage rate calculation formula, to obtain the primary coolant leakage rate of the pressurized water reactor includes: The leakage rate under full-power stable conditions is calculated based on the initial average temperature of the primary circuit, the initial liquid level of the control tank, the initial liquid level of the voltage regulator, the initial time value, the initial water volume of the TEP head tank, the initial dilution amount, the initial boronization amount, the real-time temperature of the primary circuit, the real-time liquid level of the control tank, the real-time liquid level of the voltage regulator, the real-time time value, the real-time water volume of the TEP head tank, the real-time dilution amount, and the real-time boronization amount, and in conjunction with the third leakage rate calculation formula. The real-time temperature of the primary circuit, the real-time liquid level of the control tank, the real-time liquid level of the voltage regulator, the real-time water volume of the TEP head tank, the real-time dilution amount, and the real-time boronization amount are values ​​under the second state of full-power stable conditions.

5. The method for calculating the primary coolant leakage rate of a pressurized water reactor according to claim 4, characterized in that, The step of calculating the leakage rate based on the initial time data and the real-time parameter data, combined with the leakage rate calculation formula, to obtain the primary coolant leakage rate of the pressurized water reactor includes: The density of the primary loop subcooled water under the initial state, the density of the primary loop subcooled water under the real-time state, the density of the pressurizer saturated water under the initial state, and the density of the pressurizer saturated steam under the real-time state are calculated using the water vapor graph function. The leakage rate under any given state is calculated based on the initial liquid level in the control tank, the initial liquid level in the pressure regulator, the initial time value, the initial water volume in the TEP head tank, the initial dilution amount, the initial boronization amount, the real-time liquid level in the control tank, the real-time liquid level in the pressure regulator, the real-time time value, the real-time water volume in the TEP head tank, the real-time dilution amount, the real-time boronization amount, the density of the primary loop subcooled water in the initial state, the density of the primary loop subcooled water in the real-time state, the saturated water density of the pressure regulator in the initial state, and the saturated steam density of the pressure regulator in the real-time state, combined with the fourth leakage rate calculation formula.

6. The method for calculating the primary coolant leakage rate of a pressurized water reactor according to claim 1, characterized in that, Before calculating the leakage rate of the pressurized water reactor primary coolant based on the initial time data and the real-time parameter data, combined with the leakage rate calculation formula, the following steps are included: The real-time parameter data are averaged using the moving average method.

7. The method for calculating the primary coolant leakage rate of a pressurized water reactor according to claim 1, characterized in that, The method includes: Based on the leakage rate of the primary coolant in the pressurized water reactor, a curve fitting was performed to obtain a trend diagram of the leakage rate of the primary coolant in the pressurized water reactor. The trend graph of the primary coolant leakage rate of the pressurized water reactor and the real-time data of the primary coolant leakage rate of the pressurized water reactor are displayed.

8. A pressurized water reactor primary coolant leakage rate calculation system for implementing the pressurized water reactor primary coolant leakage rate calculation method according to any one of claims 1-7, characterized in that, include: The initial value determination unit is used to determine the initial time data; The real-time data acquisition unit is used to acquire real-time parameter data of the primary loop of the pressurized water reactor. The leakage rate calculation unit is used to calculate the leakage rate based on the initial time data and the real-time parameter data, combined with the leakage rate calculation formula, to obtain the leakage rate of the pressurized water reactor primary coolant.

9. The pressurized water reactor primary coolant leakage rate calculation system according to claim 8, characterized in that, Also includes: The human-computer interaction unit is used to display the initial time data, the real-time parameter data of the pressurized water reactor primary loop, the real-time data of the pressurized water reactor primary loop coolant leakage rate, and the trend graph of the pressurized water reactor primary loop coolant leakage rate.