A method and system for online monitoring and automatic compensation of He-Xe closed Brayton cycle component deviation
Patent Information
- Application Number
- CN202610696516.7
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2026-05-20
- Publication Date
- 2026-09-01
AI Technical Summary
[0005]有鉴于此,本发明旨在提出一种He-Xe闭式布雷顿循环组分偏移的在线监测与自动补偿方法及系统,以解决在He-Xe二元工质闭式布雷顿循环长期密闭运行过程中,He-Xe闭式布雷顿循环因差异性泄漏导致的工质组分单向漂移及现有库存系统无法维持组分平衡的问题
1.首次提出一种专门用于He-Xe二元惰性气体闭式布雷顿循环(CBC)工质组分长期维持的在线补偿系统,该系统独立于仅调节工质总摩尔数的库存控制系统,由“混合分子量在线辨识模块(模块A)—漂移模式与故障辨识模块(模块B)—单组分精量注入模块(模块C)—补偿效果验证与系统协同模块(模块D)”构成四模块闭环架构。区别于通用压缩机分子量测量或防喘振保护,本发明以补偿差异性泄漏导致的组分漂移为直接控制目标,填补了现有CBC控制方法在工质“质量”维度上的空白。
Smart Images

Figure CN122676969A_ABST
Abstract
Description
Technical Field
[0001] This invention belongs to the technical field of inert gas closed Brayton cycle power generation systems, and in particular relates to an online monitoring and automatic compensation method and system for component shift in a He-Xe closed Brayton cycle. Background Technology
[0002] In a He-Xe binary working fluid closed Brayton cycle (CBC), the working fluid operates continuously within a high-pressure sealed circuit. Because He (… ) and Xe ( The molecular weights differ by approximately 33 times. According to Graham's law of flow efficiency, the ratio of their differential leakage rates through sealing surfaces (metal seals, elastomer O-rings, shaft seal gaps, etc.) is:
[0003] The leakage molar flow rate of He is approximately 5.7 times that of Xe (under ideal orifice flow conditions). In engineering scenarios dominated by elastomer sealing permeation, the measured leakage difference may be significantly higher than the predicted value of the ideal orifice model, even exceeding 10 times. During long-term operation, the average molecular weight of the mixed working fluid in the loop... It continues to drift monotonically towards the Xe-rich side, deviating from the design value (typical value). This leads to the following series of coupling performance degradations: Corrected shaft speed of the compressor The systemic performance is low, and the operating point has shifted to the left, moving out of the high-efficiency zone. heat exchanger Nu The heat transfer coefficient corresponding to the Nusselt number varies with Rise and fall ( (Based on Chapman-Enskog transport theory) When the velocity triangle between turbine stages deviates from the similar design point, polytropic efficiency degrades.
[0004] The overall thermal efficiency and specific mass (kg / kWe) of the system continued to deteriorate. Existing inventory control systems only regulate the total inventory of the working fluid and cannot address the issue of maintaining the composition of the working fluid at the "quality" level. This is particularly problematic in extreme applications where the working fluid cannot be exchanged with the outside environment (space nuclear power, underwater propulsion). Summary of the Invention
[0005] In view of this, the present invention aims to propose an online monitoring and automatic compensation method and system for component shift in a He-Xe closed Brayton cycle, in order to solve the problems of unidirectional drift of working fluid components caused by differential leakage in the He-Xe closed Brayton cycle during long-term closed operation of the He-Xe binary working fluid closed Brayton cycle and the inability of the existing inventory system to maintain component balance.
[0006] To achieve the above objectives, as one aspect of the present invention, an online monitoring and automatic compensation method and control system for the shift in working fluid composition caused by differential leakage in a He-Xe binary closed Brayton cycle are provided, comprising the following steps: S1. Collect current compressor parameters and calculate the measured equivalent speed. and measured equivalent flow rate; S2. In the pre-stored design characteristic diagram, find the corresponding design equivalent speed under the same measured equivalent flow rate. ; S3, using measured equivalent rotational speed Design equivalent speed and the design of the average molecular weight of the mixed working fluid The average molecular weight of the first mixed working fluid was calculated. ; S4. By monitoring the cold end outlet temperature of the regenerator Deviation from design value This yields a second independent estimate of the average molecular weight of the mixed working fluid, i.e., the second average molecular weight of the mixed working fluid. ; S5, to and Perform a consistency check and select those that meet the threshold range. and By performing a weighted average, the average molecular weight estimate of the mixed working fluid is obtained. Output to S6; S6, Utilization Calculate the molecular weight drift rate With the total voltage decay rate of the system ratio , if Within the tolerance range, after calculating the amount of moles that need to be added, proceed to step S7; S7. Inject gas into the main circuit. After injection, verify the compensation effect and coordinate the collaborative operation of this subsystem and the inventory control system.
[0007] Furthermore, in S1, the current compressor parameters include the inlet temperature. Inlet pressure Export pressure and physical rotation speed .
[0008] Furthermore, in S1, the measured equivalent rotational speed...
[0009] Among them, for He-Xe single-atom ideal gas , Is it for the design of the average molecular weight of the mixed working fluid? R is the universal gas constant, 8.314 J / (mol·K).
[0010] Furthermore, in S3, the average molecular weight of the current mixed working fluid The calculation formula is: .
[0011] Furthermore, in S4, the average molecular weight of the second mixed working fluid The function is:
[0012] Where Nu is the Nusselt number, Re is the Reynolds number, Pr is the Prandtl number, cp is the specific heat capacity at constant pressure, μ is the dynamic viscosity, λ is the thermal conductivity, and R is the universal gas constant.
[0013] Furthermore, in S5, the convergence determination method is as follows:
[0014] like (e.g., 5%): Two-way fusion, Take the weighted average and output it to S6; like : Triggers sensor fault alarm, does not output gas replenishment command, waits for manual verification.
[0015] Furthermore, in S6, the aforementioned , , According to Graham's law of flow, the theoretical ratio of the two under normal differential leakage is a constant. :
[0016] like That is, within ±20% tolerance: normal differential leakage, output gas replenishment command to S7; like ,or Abrupt change: Component failure and leakage, output fault alarm, prohibit gas replenishment, trigger safety interlock; like Sensor drift triggers calibration request; in, This refers to the molecular weight drift rate. Let MXe = 131.3 g / mol, MHe = 4 g / mol, xHe,0 be the mole fraction of helium at the initial moment, and Ktheory be the theoretical ratio constant under normal differential leakage.
[0017] Furthermore, in S6, the amount of moles that need to be replenished...
[0018]
[0019] in, , The total volume of the loop is denoted as the loop average temperature, ntotal as the total number of moles of working fluid in the loop, psys as the total system pressure, Vloop as the total loop volume, R as the universal gas constant, and Tavg as the loop average temperature.
[0020] Furthermore, in S7, the method for verifying the compensation effect is as follows: After injection, wait for the mixing and homogenization time. :
[0021]
[0022] in For the working fluid mass flow rate in the loop, The characteristic time for the working fluid to complete one full circuit cycle. The average density of the working fluid is used; after the waiting period, the output will be restarted. ; like Compensation successful. Entering normal monitoring mode. Resetting the timer. If the deviation does not converge: calculate the residual deviation, return to S6 to perform the second gas replenishment; set the maximum number of gas replenishment attempts, and trigger a system alarm if the limit is exceeded.
[0023] As another aspect of the present invention, an online monitoring and automatic compensation system for component shift in a He-Xe closed Brayton cycle is provided, comprising: The online molecular weight identification module for the mixture uses measurable operating parameters from both the compressor and heat exchanger sides to determine the average molecular weight of the current working fluid mixture via at least two independent paths. After passing the consistency check, the estimated average molecular weight of the mixed working fluid is obtained and output to the drift mode identification and compensation decision module. The drift pattern identification and compensation decision module utilizes Calculate the molecular weight drift rate With the total voltage decay rate of the system ratio ,if Within the tolerance range, after calculating the required molar amount to be replenished, the gas replenishment command is output to the single-component gas precision injection execution module; The single-component gas precision injection execution module receives the gas replenishment command and injects the stored single-component pure gas into the main circuit in a controlled manner to achieve precise quality compensation. The compensation effect verification and system coordination module verifies the compensation effect after injection and coordinates the collaborative operation of this subsystem and the inventory control system to prevent mutual interference between their control actions.
[0024] Beneficial effects: 1. This invention proposes for the first time an online compensation system specifically designed for the long-term maintenance of working fluid components in a He-Xe binary inert gas closed Brayton cycle (CBC) system. This system is independent of inventory control systems that only adjust the total molar number of the working fluid. It comprises a four-module closed-loop architecture: "Online Mixed Molecular Weight Identification Module (Module A) – Drift Mode and Fault Identification Module (Module B) – Single-Component Precision Injection Module (Module C) – Compensation Effect Verification and System Coordination Module (Module D)". Unlike general compressor molecular weight measurement or surge protection, this invention directly controls component drift caused by differential leakage, filling a gap in existing CBC control methods regarding the "quality" dimension of the working fluid.
[0025] 2. A dual-path molecular weight fusion estimation method is proposed, combining the compressor similarity criterion method and the regenerator terminal difference verification method. The two paths are based on different thermodynamic transport mechanisms (compressor equivalent speed and heat exchanger Nu number deviation). Online decoupling diagnosis of "true component drift" and "sensor failure" is achieved by calculating the consistency index ε. When ε exceeds a threshold, the system prohibits gas injection and issues a sensor failure alarm, thereby avoiding erroneous compensation actions caused by misreading of a single sensor. This fault-tolerant mechanism has substantial advantages in long-life unattended operation.
[0026] 3. A method utilizing molecular weight drift rate is proposed. With the total voltage decay rate of the system A method for quantitatively identifying leakage modes using the ratio of [theoretical ratio] to [other parameters]. Based on Graham's effective flow law, under normal differential leakage conditions, the theoretical ratio [is...]. This is a definite constant determined solely by the initial composition. The control system calculates the measured ratio in real time. / and with Comparison: If (Within ±20% tolerance) is considered a normal differential leak, and gas replenishment is permitted; if or A sudden change indicates a component failure and leak; therefore, gas replenishment is prohibited, and the safety interlock is triggered. The system was identified as sensor drift and calibration was requested. This innovation expands simple component control to include system health diagnostics and safety protection functions, which is particularly critical for long-life, sealed systems such as space nuclear power plants.
[0027] 4. A method for accurately calculating feedforward gas injection volume based on loop state equations is proposed. During injection, the actual injection volume is obtained in real time using the tank pressure drop metering method. After injection, the mixture undergoes a waiting period for homogenization. The molecular weight convergence is then re-verified by module A. This forms a three-level control closed loop of "feedforward calculation - pressure drop measurement closed loop - post-compensation verification", which, unlike traditional timed pulse or simple feedback injection methods, achieves precise control of the injection volume and traceability of the effect.
[0028] 5. A collaborative decoupling control strategy is proposed between the component compensation system and the inventory control system. The orthogonal control objectives of the two systems are clearly defined (component compensation adjustment...). M Inventory control adjustment However, physical coupling exists (He supplementation slightly increases the total pressure). This invention designs a "lock-and-unlock" communication protocol: before performing component compensation injection, a "component compensation lock" signal is sent to the inventory controller to pause its automatic adjustment; pending verification by module D. After convergence, the lockout is released. If the total system pressure exceeds the inventory control limit due to gas replenishment, the pressure correction will be completed first by the pressure relief valve of the inventory control system, and this system will not intervene in the total pressure regulation. This strategy achieves orderly coordination between the "quantity" and "quality" control of the working fluid, avoiding action conflicts between multiple controllers. Attached Figure Description
[0029] Figure 1 For flowcharts; Figure 2 For pre-stored design feature diagrams. Detailed Implementation Specific implementation method one: Referring to the accompanying drawings, this embodiment provides an online monitoring and automatic compensation method for component shift in a He-Xe closed Brayton cycle, comprising the following steps: S1. Collect current compressor parameters and calculate the measured equivalent speed. and measured equivalent flow rate; S2. In the pre-stored design characteristic diagram, find the corresponding design equivalent speed under the same measured equivalent flow rate. ; S3, using measured equivalent rotational speed Design equivalent speed and the design of the average molecular weight of the mixed working fluid The average molecular weight of the first mixed working fluid was calculated. ; Current compressor parameters include inlet temperature Inlet pressure Export pressure Physical rotation speed ; For He-Xe monatomic ideal gas ( The compressor's equivalent speed is defined as:
[0031] physical rotation speed Unchanged, inlet temperature Under measurable conditions, the following relationship exists between the two states:
[0032] Calculate the measured equivalent rotational speed using the above formula. ; like Figure 2 As shown, in the pre-stored design characteristic map (compressor map), the inlet pressure is used. Export pressure Find the corresponding design equivalent rotational speed for the same corrected mass flow. ; The four dashed lines correspond to the 100% equivalent speed characteristic lines of the compressor with MW=15 / 28 / 40 / 55g / mol. Based on the similarity criterion... The larger the molecular weight, the more the characteristic curve shifts to the lower left (reducing the equivalent flow rate and pressure ratio).
[0033] Visualization of diagnostic logic: Orange dot: Design working point (MW=40, everything is normal) Purple square: Measured operating point after He leakage (equivalent physical rotation speed) n , n The speed remains unchanged, but the equivalent speed decreases, and the operating point falls along the "MW=55" characteristic line. Red dashed line: The actual drift trajectory of the working point. The difference in equivalent flow rate between two operating points is a measurable diagnostic indicator.
[0034] Based on the measured equivalent speed and the design equivalent speed Calculate .
[0035] S4. By monitoring the cold end outlet temperature of the regenerator Deviation from design value This yields a second independent estimate of the average molecular weight of the mixed working fluid, i.e., the second average molecular weight of the mixed working fluid. ; Based on Chapman-Enskog transport theory, the thermal conductivity of He-Xe mixed gases and dynamic viscosity All The function, which in turn affects the heat exchanger tube side Nu Numbers and Pr number:
[0036] By monitoring the cold end outlet temperature of the regenerator Deviation from design value ,get Second independent estimate .
[0037] S5, to and Perform a consistency check and select those that meet the threshold range. and By performing a weighted average, the average molecular weight estimate of the mixed working fluid is obtained. Output to S6; Two path pairs The sensitivity indices differ (compressor path index 1 / 2, heat exchanger path index approximately 0.4), and the two estimation results should converge under normal operating conditions. Define a diagnostic consistency index:
[0038] like Two routes merged. Take the weighted average and output it to S6; It is 5%.
[0039] like : Triggers sensor fault alarm, does not output gas replenishment command, waits for manual verification.
[0040] This mechanism decouples the identification of "real component drift" and "sensor misreading", avoiding false gas replenishment caused by a single sensor failure.
[0041] S6, Utilization Calculate the molecular weight drift rate With the total voltage decay rate of the system ratio ,if Within the tolerance range, after calculating the amount of moles that need to be added, proceed to step S7; Assuming that component drift is confirmed to exist, distinguish whether the drift is "normal differential flow leakage" or "abnormal component failure leakage", and output corresponding control commands accordingly.
[0042] Drift Pattern Identification – Dual-Rate Decoupling Judgment Define two time derivatives that can be computed online: Molecular weight drift rate: (Obtained by sliding fit of the time series output from module A) System total voltage decay rate: (The absolute pressure on the low-pressure side of the system is measured directly by the sensor) According to Graham's law of flow, the theoretical ratio of the two under normal differential leakage is a constant. (with initial components only) Related):
[0043] Judgment logic: If the measured ratio (Within ±20% tolerance): For normal differential leakage, output a gas replenishment command to module C. If the measured ratio or Abrupt change: Component failure and leakage, output fault alarm, gas supply prohibited, safety interlock triggered. like (Molecular weight shift but total pressure remains essentially unchanged): Sensor drift triggers calibration request. Feedforward calculation of gas replenishment volume After confirming normal drift, based on the system state equation (He-Xe compressibility factor under operating pressure) Calculate the total number of moles of working fluid in the current loop:
[0044] in This is the total volume of the loop (design constant). This is the average temperature of the circuit (weighted by multiple temperature measurements). The required molar amount of He to be added is:
[0045] This value is output to S7 as a feedforward control setpoint.
[0046] S7. Inject gas into the main circuit. After injection, verify the compensation effect and coordinate the collaborative operation of this subsystem and the inventory control system.
[0047] The gas replenishment command injects stored single-component pure gas into the main circuit in a controlled manner to achieve precise quality compensation. The injection metering method utilizes the pressure drop in the storage tank before and after injection. Calculate the actual number of moles injected without an additional flow meter:
[0048] This measurement value is fed back in real time to verify the compensation effect, forming a closed-loop verification of the injection volume.
[0049] After injection, wait for the mixing and homogenization time. :
[0050] in The characteristic time for the working fluid to complete one full circuit cycle. This represents the average density of the working fluid. After the waiting period, module A will output again. : like (like Compensation successful. Entering normal monitoring mode. Resetting the timer. If the deviation does not converge: calculate the residual deviation, return to module B to perform a second gas replenishment; set the maximum number of gas replenishments (e.g., 3 times), and trigger a system alarm if the limit is exceeded.
[0051] Decoupling from the inventory control system; The control quantity of the inventory control system is the total loop pressure. (Adjust the total number of moles of working fluid) The controlled quantity of this invention is the mixed molecular weight. (Adjusting the molar ratio of the working fluid components). The two control objectives are orthogonal, but they are coupled: the gas injection operation increases the molar number of He while also slightly improving... .
[0052] Cooperative logic: During the gas replenishment operation of this invention, a "component compensation lock" signal is sent to the inventory control system to suspend the automatic adjustment action of the inventory controller, pending... The lock will be released after convergence. (If this is due to Qi replenishment...) If the inventory exceeds the control limit, the pressure will be corrected by the pressure relief valve of the inventory control system first. The controller of this invention does not interfere with the total pressure regulation. Specific implementation method two; An online monitoring and automatic compensation system for component shift in a He-Xe closed Brayton cycle includes: The online molecular weight identification module for the mixture uses measurable operating parameters from both the compressor and heat exchanger sides to determine the average molecular weight of the current working fluid mixture via at least two independent paths. After passing the consistency check, the estimated average molecular weight of the mixed working fluid is obtained and output to the drift mode identification and compensation decision module. The drift pattern identification and compensation decision module utilizes Calculate the molecular weight drift rate With the total voltage decay rate of the system ratio ,if Within the tolerance range, after calculating the required molar amount to be replenished, the gas replenishment command is output to the single-component gas precision injection execution module; The single-component gas precision injection execution module receives the gas replenishment command and injects the stored single-component pure gas into the main circuit in a controlled manner to achieve precise quality compensation. The hardware controlled by the single-component gas precision injection execution module includes: He high-pressure gas storage assembly: high-pressure He reservoir + primary pressure regulator + shut-off valve. The reservoir's rated pressure is higher than the maximum operating pressure of the main circuit, ensuring self-pressure differential-driven injection under any operating condition.
[0054] Xe high-pressure gas storage component (optional): configured only in Xe priority leakage scenarios or when bidirectional compensation is required, with the same structure as above.
[0055] Precision mass flow control valve (MFC): Installed on the injection branch from each storage tank to the loop, executing the commands issued by module B. Commands; use pulse width modulation (PWM) or proportional-integral (PI) control modes.
[0056] One-way check valve: Installed on the injection branch to prevent the high-pressure working medium in the main circuit from flowing back into the storage tank when not being injected.
[0057] Injection interface location: Low-pressure section of the main circuit (upstream of the compressor inlet, downstream of the cooler / radiant radiator outlet). This location features low working fluid temperature, low pressure, and relatively uniform flow velocity, which facilitates rapid mixing of the injected gas with the mainstream and minimizes disturbance to the turbine's mechanical thermodynamic state.
[0058] The compensation effect verification and system coordination module verifies the compensation effect after injection and coordinates the collaborative operation of this subsystem and the inventory control system to prevent mutual interference between their control actions.
[0059] The embodiments of the present invention disclosed above are merely illustrative of the invention. These embodiments do not exhaustively describe all details, nor do they limit the invention to the specific implementations described. Many modifications and variations can be made based on the content of this specification. This specification selects and specifically describes these embodiments to better explain the principles and practical applications of the invention, thereby enabling those skilled in the art to better understand and utilize the invention.
Claims
1. An online monitoring and automatic compensation method for He-Xe closed Brayton cycle component deviation, characterized in that, Includes the following steps: S1, collect current compressor parameters, calculate measured equivalent speed and measured equivalent flow rate; S2. In the pre-stored design characteristic diagram, find the corresponding design equivalent speed under the same measured equivalent flow rate. ; S3, using measured equivalent rotational speed Design equivalent speed and the design of the average molecular weight of the mixed working fluid The average molecular weight of the first mixed working fluid was calculated. ; S4. By monitoring the cold end outlet temperature of the regenerator Deviation from design value This yields a second independent estimate of the average molecular weight of the mixed working fluid, i.e., the second average molecular weight of the mixed working fluid. ; S5, to and Perform a consistency check and select those that meet the threshold range. and By performing a weighted average, the average molecular weight estimate of the mixed working fluid is obtained. Output to S6; S6, Utilization Calculate the molecular weight drift rate With the total voltage decay rate of the system ratio , if Within the tolerance range, after calculating the amount of moles that need to be added, proceed to step S7; S7. Inject gas into the main circuit. After injection, verify the compensation effect and coordinate the collaborative operation of this subsystem and the inventory control system.
2. The method for online monitoring and automatic compensation of component shift in a He-Xe closed Brayton cycle according to claim 1, characterized in that: In S1, the current compressor parameters include the inlet temperature. Inlet pressure Export pressure and physical rotation speed .
3. The method for online monitoring and automatic compensation of component shift in a He-Xe closed Brayton cycle according to claim 2, characterized in that: In S1, the measured equivalent rotational speed Among them, for He-Xe single-atom ideal gas , Is it for the design of the average molecular weight of the mixed working fluid? , R The universal gas constant is 8.314 J / (mol·K).
4. The method for online monitoring and automatic compensation of component shift in a He-Xe closed Brayton cycle according to claim 1, characterized in that: In S3, the average molecular weight of the current mixed working fluid The calculation formula is: 。 5. The method for online monitoring and automatic compensation of component shift in a He-Xe closed Brayton cycle according to claim 1, characterized in that: In S4, the average molecular weight of the second mixed working fluid The function is: in, Nu For Nusselt numbers, Re The Reynolds number is... Pr For Prandtl numbers, c p For isobaric specific heat capacity, μ For dynamic viscosity, λ Thermal conductivity, R This is the universal gas constant.
6. The method for online monitoring and automatic compensation of component shift in a He-Xe closed Brayton cycle according to claim 1, characterized in that: In S5, the convergence determination method is as follows: like (e.g., 5%): Two-way fusion, Take the weighted average and output it to S6; like : Triggers sensor fault alarm, does not output gas replenishment command, waits for manual verification.
7. The method for online monitoring and automatic compensation of component shift in a He-Xe closed Brayton cycle according to claim 1, characterized in that: In S6, the , , According to Graham's law of flow, the theoretical ratio of the two under normal differential leakage is a constant. : like That is, within ±20% tolerance: normal differential leakage, output gas replenishment command to S7; like ,or Abrupt change: Component failure and leakage, output fault alarm, prohibit gas replenishment, trigger safety interlock; like Sensor drift triggers calibration request; in, This refers to the molecular weight drift rate. The total voltage decay rate of the system. M Xe =131.3 g / mol M He =4g / mol x He,0 The initial mole fraction of helium. K theory This is the theoretical ratio constant under normal differential leakage conditions.
8. The method for online monitoring and automatic compensation of component shift in a He-Xe closed Brayton cycle according to claim 1, characterized in that: In S6, the amount of moles that need to be replenished in, , The total volume of the loop is The average temperature of the circuit. n total The total number of moles of working fluid in the loop. p sys For the total system pressure, V loop The total volume of the loop is R This is the universal gas constant. T avg This represents the average temperature of the circuit.
9. The method for online monitoring and automatic compensation of component shift in a He-Xe closed Brayton cycle according to claim 1, characterized in that: In S7, the method for verifying the compensation effect is as follows: After injection, wait for the mixing and homogenization time. : in For the working fluid mass flow rate in the loop, The characteristic time for the working fluid to complete one full circuit cycle. The average density of the working fluid is used; after the waiting period, the output will be restarted. ; like Compensation successful. Entering normal monitoring mode. Resetting the timer. If the deviation does not converge: calculate the residual deviation, return to S6 to perform the second gas replenishment; set the maximum number of gas replenishment attempts, and trigger a system alarm if the limit is exceeded.
10. An online monitoring and automatic compensation system for component shift in a He-Xe closed Brayton cycle, characterized in that, include: The online molecular weight identification module for the mixture uses measurable operating parameters from both the compressor and heat exchanger sides to determine the average molecular weight of the current working fluid mixture via at least two independent paths. After passing the consistency check, the estimated average molecular weight of the mixed working fluid is obtained and output to the drift mode identification and compensation decision module. The drift pattern identification and compensation decision module utilizes Calculate the molecular weight drift rate With the total voltage decay rate of the system ratio ,if Within the tolerance range, after calculating the required molar amount to be replenished, the gas replenishment command is output to the single-component gas precision injection execution module; The single-component gas precision injection execution module receives the gas replenishment command and injects the stored single-component pure gas into the main circuit in a controlled manner to achieve precise quality compensation. The compensation effect verification and system coordination module verifies the compensation effect after injection and coordinates the collaborative operation of this subsystem and the inventory control system to prevent mutual interference between their control actions.