Safety state quantity screening and dynamic optimization control method for hydrogen-electricity coupling system

Abstract

The invention discloses a safety state quantity screening and dynamic optimization control method for a hydrogen-electric coupling system. The method comprises the following steps: establishing a joint constraint condition according to a hydrogen-electric coupling system model and an accident set; screening out a state quantity which does not meet a multi-energy flow balance constraint condition, a state quantity which does not meet a feasible domain constraint condition and a state quantity which does not meet a safety domain constraint condition of the hydrogen-electricity coupling system under an 'N-1' condition, and taking the state quantities as to-be-adjusted state quantities; by utilizing the state quantities except the to-be-adjusted state quantity, forming an optimization control model of the hydrogen-electricity coupling system by taking the optimal efficiency degree and safety degree of the feasible region as an optimization target model and the safety region constraint condition of the hydrogen-electricity coupling system under the N-1 condition; and based on the optimization control model, obtaining a state quantity optimization value which simultaneously meets the optimization target model and the safety domain constraint condition of the hydrogen-electric coupling system under the N-1 condition, and controlling the hydrogen-electric coupling system according to the state quantity optimization value, thereby realizing optimization control of the hydrogen-electric coupling system, and predicting the safety and stability of the system in advance.

Description

Technical Field

[0001] The present invention relates to the technical field of hydrogen-electricity coupling systems, and in particular to a safety state quantity screening and dynamic optimization control method for hydrogen-electricity coupling systems. Background Art

[0002] Hydrogen energy is an ideal clean secondary energy source. In recent years, as hydrogen energy application technology has gradually matured, electric-hydrogen coupling technology has continued to advance, enabling hydrogen to be quickly integrated into the entire micro-energy system.

[0003] In the prior art, the optimization control method of the hydrogen-electric coupling system includes: by conducting a specific analysis of the three aspects of hydrogen production, storage, and transmission in the target area, taking the hydrogen-electric coupling state operation rationality index as the optimization index, and improving the accuracy of data analysis; for the hydrogen energy conversion system, establishing the objective function and constraint conditions that consider the operating cost and efficiency penalty coefficient, taking the power in the hydrogen conversion system as the optimization variable, and finally operating the comprehensive energy conversion system with the optimized variable; optimizing the capacity configuration of hydrogen energy storage, establishing a double-layer model, and passing the frequency fluctuation data as output to the outer optimization model when the inner model is given, and obtaining the optimal hydrogen energy storage capacity configuration after calculating the objective function; obtaining the internal temperature distribution data of the solid oxide electrolyzer (SOEC) stack, analyzing the relationship between temperature and electrical performance under different parameter combinations, and using the neural network predictive control algorithm to control the temperature of the SOEC system to achieve the purpose of improving the safety and efficiency of the system. Modeling the hydrogen energy subsystem, most micro-energy networks currently do not consider hydrogen energy, but only consider simple wind and light models. At present, most safety domains are displayed as plane graphics or hyperplanes or hyperplane sections, which are not as intuitive as three-dimensional images. No effective means are used to optimize and adjust the nodes beyond the safety domain to ensure the normal operation of the entire system. Summary of the invention

[0004] In order to solve the deficiencies in the prior art, the present invention provides a safety state quantity screening and dynamic optimization control method for a hydrogen-electricity coupling system, so as to achieve optimized control of the hydrogen-electricity coupling system and predict the safety and stability of the system in advance.

[0005] The technical solution of the present invention is as follows.

[0006] The present invention proposes a safety state quantity screening and dynamic optimization control method for a hydrogen-electric coupling system, wherein the hydrogen-electric coupling system comprises: an electric energy subsystem, a thermal subsystem and a hydrogen energy subsystem, including: Collect the state quantity of the hydrogen-electric coupling system, establish a hydrogen-electric coupling system model, and obtain fault data under various fault scenarios based on the hydrogen-electric coupling system model to construct an accident set; Establish joint constraints based on the hydrogen-electricity coupling system model and the accident set, including: establishing multi-energy flow balance constraints based on the hydrogen-electricity coupling system model; establishing feasible domain constraints of the hydrogen-electricity coupling system under the "N-1" condition based on the hydrogen-electricity coupling system model and the multi-energy flow balance constraints, with the goal of maximizing the energy supply range of the hydrogen-electricity coupling system under the "N-1" condition; establishing safety domain constraints of the hydrogen-electricity coupling system under the "N-1" condition based on the feasible domain constraints and accident set of the hydrogen-electricity coupling system under the "N-1" condition; Screen out the state quantities that do not satisfy the multi-energy flow balance constraint conditions, screen out the state quantities that do not satisfy the feasible domain constraint conditions of the hydrogen-electric coupling system under the "N-1" condition, and screen out the state quantities that do not satisfy the safety domain constraint conditions of the hydrogen-electric coupling system under the "N-1" condition as the state quantities to be adjusted; By utilizing state quantities other than the state quantities to be adjusted, the optimization target model is taken as the optimal efficiency and safety degree of the feasible domain; the optimization control model of the hydrogen-electric coupling system is constructed by utilizing the optimization target model and the safety domain constraints of the hydrogen-electric coupling system under the "N-1" condition; based on the optimization control model, the state quantities to be adjusted are optimized to obtain the optimized values ​​of the state quantities that simultaneously satisfy the optimization target model and the safety domain constraints of the hydrogen-electric coupling system under the "N-1" condition, and the hydrogen-electric coupling system is controlled according to the optimized values ​​of the state quantities.

[0007] Preferably, the state quantities of the hydrogen-electricity coupling system include: grid voltage, grid frequency and heating network temperature.

[0008] Preferably, the electric energy subsystem includes conventional thermal power units and new energy units, and the electric energy subsystem model includes: 1) Output constraint model of conventional thermal power units: , In the formula, for t Conventional thermal power units of efforts, and Conventional thermal power units The lower and upper limits of output; 2) Output constraint model of new energy units: , In the formula, for t New energy units are always producing power. and They are respectively the lower and upper limits of the output of new energy units.

[0009] Preferably, the thermal subsystem includes a thermal network and a hydraulic network, and the thermal subsystem model includes: 1) Thermal network model: , In the formula, is the association matrix between nodes and branches in the thermal network; is the association matrix between loops and branches in the thermal network; is the mass flow column vector of the thermal pipeline; is the injected mass flow column vector of the source node and load node in the thermal network; is the pressure loss column vector of the thermal pipeline; is the resistance coefficient column vector of the thermal pipe, which depends on the pipe diameter; 2) Hydraulic network model: , In the formula, is the column vector of thermal power consumed or supplied to the node; is the specific heat capacity of water; is the column vector of water supply temperature; is the column vector of outlet water temperature; , are the initial end node temperature and the terminal node temperature of the hot water pipe respectively; is the ambient temperature; is the total heat transfer coefficient per unit length of the hot water pipe; is the length of the hot water pipe; are the mass flow column vectors of each hot water pipe flowing into and out of the node respectively; is the heat medium temperature of the node; is the heat medium temperature at the end of each pipe section flowing into the node.

[0010] Preferably, the hydrogen energy subsystem includes an electrolyzer, a fuel cell and a hydrogen storage tank, and the hydrogen energy subsystem model includes: 1) Electrolyzer model: , In the formula, is the actual input power of the electrolyzer; input power to the electrolyzer; is the efficiency of the electrolyzer; is the hydrogen production power of the electrolyzer; is the hydrogen production efficiency of the electrolyzer; The electrolyzer model includes power and ramp constraints in the operating state: , In the formula, To indicate that the electrolytic cell is t A 0-1 variable indicating the start and stop status at all times; for tThe electrolyzer input power at the moment; , are the lower and upper limits of the electrolyzer input power, respectively; 2) Fuel cell model: , In the formula, is the output power of the fuel cell; is the input power of the fuel cell; for the efficiency of the fuel cell; The fuel cell model also includes the operating constraints of the fuel cell: , In the formula, To indicate that the fuel cell t A 0-1 variable indicating the start and stop status at all times; , are the lower and upper limits of the fuel cell input power, respectively; 3) Hydrogen storage tank model: , In the formula, The energy stored in the hydrogen tank at the initial moment; is the capacity of the hydrogen storage tank; for t The energy stored in the hydrogen tank at all times; , They are t The charging and discharging power of hydrogen storage tank at all times; , They are the charging and discharging efficiency of the hydrogen storage tank respectively; For the operating period; is the maximum capacity factor of the hydrogen storage tank; , They represent the hydrogen storage tank t A 0-1 variable indicating the state of charge and discharge at any given moment; , They are the maximum values ​​of charging and discharging power of the hydrogen storage tank respectively.

[0011] Preferably, the failure scenarios include: failure of key pipeline outlets at energy hubs, failure of key equipment at energy hubs, and failure of imbalance between supply and demand of renewable energy.

[0012] Preferably, the fault data under each fault scenario is obtained based on the hydrogen-electricity coupling system model to construct an accident set, including: Probabilistic modeling of the output of new energy units is carried out, and a joint model of the output of new energy units is formed with the output constraint model of new energy units, thereby updating the electric energy subsystem model and obtaining an updated hydrogen-electric coupling system model; Among them, the probability modeling of the output of new energy units includes: establishing a probability distribution model of wind turbine output based on a two-parameter Weibull random probability model, and establishing a probability distribution model of photovoltaic generator output based on a Beta probability distribution model; Based on the updated hydrogen-electric coupling system model, the fault data under each fault scenario is obtained to construct an accident set.

[0013] Preferably, the probability distribution model of wind turbine output is established based on the two-parameter Weibull random probability model to satisfy the following relationship: , , In the formula, is the probability density function of the actual wind speed, k w is the shape factor of the fan unit, c w is the scale factor of the fan unit; , V h , V out , V in They are rated wind speed, actual wind speed, cut-out wind speed and cut-in wind speed respectively; P wmax is the maximum output power of the wind turbine, The output of wind turbines based on probability distribution; The probability distribution model of photovoltaic generator output based on the Beta probability distribution model satisfies the following relationship: , , In the formula, is the probability density function of actual illumination, is the shape factor of the photovoltaic module, β p is the scale factor of the PV module; , They are rated light intensity and actual light intensity respectively; is the Gamma function; Indicates the maximum output power of the photovoltaic generator set. Photovoltaic generator output based on probability distribution.

[0014] Preferably, the multi-energy flow balance constraints include: 1) Power balance constraint, satisfying the following relationship: , In the formula, , They are t Time Node i Active injection power and reactive injection power; N is the number of system nodes; , They are the node admittance matrices i OK j The real and imaginary parts of the columns; for t Time Node i The voltage amplitude of for t Time Node j The voltage amplitude of for t Time Branch ij The phase angle difference; , Node i The lower and upper limits of the voltage amplitude; 2) Hydrogen power balance constraint, satisfying the following relationship: , In the formula, for t The hydrogen production power of the electrolyzer at any moment; , They are t The charging and discharging power of hydrogen storage tank at all times; for t Hydrogen load at the time; for t The input power of the fuel cell at any moment; 3) Power balance constraints of the electric-hydrogen coupling system, including electric-hydrogen coupling operation constraints and electric hydrogen production power constraints. The electric-hydrogen coupling operation constraints are as follows: , In the formula, for t Conventional thermal power units contribution; is the number of conventional thermal power units; for t New energy units are always producing; for t The output power of the fuel cell at any moment; for t Time Node The electrical load; is the hydrogen load of the system; for t The electrolyzer input power at the moment; is the number of nodes; The power constraints for hydrogen production are as follows: , In the formula, for Shike New Energy Station Output power; for Shike New Energy Station The load power provided; for Shike New Energy Station The abandoned power; for Shike New Energy Station Input power of the inner electrolyzer; 4) Thermal energy balance constraint, satisfying the following relationship: , In the formula, C CHP is the thermoelectric ratio; P CHP,h is the thermal power; P CHP,e is the electrical power; 5) Flow constraints of hydraulic network: , In the formula, B is the correlation matrix of the hydraulic grid; w is the branch flow vector; K is the pipe coefficient.

[0015] Preferably, the feasible domain constraints of the hydrogen-electric coupling system under the "N-1" condition satisfy the following relationship: , In the formula, For the The state quantity of each pipeline segment, The sum of the state quantities of all pipeline segments is the largest; The equation constraints for the key pipeline segments and key equipment in the feasible domain of the hydrogen-electric coupling system under the “N-1” condition are: The key pipeline sections and key equipment constraints of the feasible domain of the hydrogen-electric coupling system under the "N-1" condition; , are the lower and upper bound vectors of the inequality constraints respectively; , M is the number of pipeline segments.

[0016] Preferably, the safety domain constraint of the hydrogen-electric coupling system under the "N-1" condition satisfies the following relationship: , In the formula, It is the security domain constraint; It is the equation constraint condition for the key pipeline segments and key equipment in the safety domain of the hydrogen-electric coupling system under the “N-1” condition; These are the inequality constraints for key pipeline segments and key equipment in the safety domain of the hydrogen-electric coupling system under the “N-1” condition.

[0017] Preferably, according to the multi-energy flow balance constraints, the feasible domain constraints of the hydrogen-electric coupling system under the "N-1" condition, and the safety domain constraints, the isolation forest algorithm is used to screen the state quantities, not only to remove abnormal data in the state quantities, but also to extract the state quantities to be adjusted.

[0018] Preferably, the efficiency function of the feasible region satisfies the following relationship: , In the formula, is the efficiency function of the feasible domain, is the total number of state quantities other than the state quantity to be adjusted. In the embodiment, after the state quantity is screened, the hydrogen-electric coupling system The state quantity of each pipeline segment is the object, The state quantity to be adjusted is The state quantity of each pipeline segment, For the The maximum value of the state quantity of each pipeline segment under the constraints of the feasible region; The safety function of the feasible region satisfies the following relationship: , In the formula, is the safety function of the feasible domain; is the safety distance from the safe working point c in the safety domain to the upper safety boundary; is the safety distance from the minimum energy supply capacity working point z in the safety domain to the upper safety boundary, where the minimum energy supply capacity refers to the minimum load carrying capacity of the system. is the coupling coefficient of the hydrogen-electric coupling system, .

[0019] Preferably, the proportion of the number of coupled devices in the hydrogen-electricity coupling system to the total number of devices is used as the coupling degree coefficient. When the proportion of the number of coupled devices in the hydrogen-electricity coupling system to the total number of devices is less than or equal to 20%, , when the number of coupled devices in the hydrogen-electricity coupling system accounts for more than 20% and less than or equal to 50% of the total number of devices .

[0020] Preferably, the optimization control model of the hydrogen-electric coupling system is constructed by using the optimization target model and the safety domain constraint of the hydrogen-electric coupling system under the "N-1" condition, satisfying the following relationship: , In the formula, The optimization target model that represents the optimal efficiency function and safety function of the feasible domain is represents the safety domain constraint of the hydrogen-electric coupling system under the “N-1” condition, where: O is the safety domain constraint of the hydrogen-electric coupling system under the “N-1” condition, L is the state quantity set of each pipeline segment.

[0021] Preferably, based on the optimization control model, the state quantity to be adjusted is optimized, including: Based on the three-dimensional image method, the three-dimensional safety domain corresponding to the safety domain constraints of the hydrogen-electric coupling system under the "N-1" condition is drawn; The Manhattan distance method is used to calculate the safety distance between the state quantity to be adjusted and the boundary of the three-dimensional safety zone as the adjustment distance of the state quantity to be adjusted.

[0022] The beneficial effects of the present invention are that, compared with the prior art, at least the following are included: the hydrogen-electricity coupling micro-energy network optimization control method proposed in the present invention can optimize the control of the state of the hydrogen-electricity coupling micro-energy network, build a model based on the micro-energy network, simulate the actual micro-energy network, and provide a basic framework support for subsequent steps; obtain multi-energy flow balance constraints, feasible domain constraints and safety domain constraints through the model and the accident set, and provide constraints for the subsequent safety domain and feasible domain; obtain the safety domain and feasible domain through the constraints, and can screen the state quantity data of the micro-energy network according to the safety domain and feasible domain, and can predict the safety and stability of the micro-energy network system in advance; according to the analysis and optimization results of the state quantity data, the optimization control of the micro-energy network can be realized. The present invention is helpful to study the optimization control of the safe and efficient operation of the hydrogen-electricity coupling system, can improve the safety of power operation, and has good practical significance. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] Figure 1 It is a flow chart of the safety state quantity screening and dynamic optimization control method for the hydrogen-electricity coupling system proposed by the present invention; Figure 2 In the embodiment of the present invention, the power pipeline segment is used , thermal pipeline segment , Hydrogen pipeline segment Schematic diagram of the three-dimensional safety domain composed of energy supply loads; Figure 3 It is a schematic diagram of a process of optimizing the state quantity to be adjusted in an embodiment of the present invention. DETAILED DESCRIPTION

[0024] In order to make the purpose, technical scheme and advantages of the present invention clearer, the technical scheme of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. The embodiments described in this application are only embodiments of a part of the present invention, rather than all embodiments. Based on the spirit of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work belong to the protection scope of the present invention.

[0025] The present invention proposes a safety state quantity screening and dynamic optimization control method for a hydrogen-electric coupling system, wherein the hydrogen-electric coupling system comprises: an electric energy subsystem, a thermal subsystem and a hydrogen energy subsystem; Figure 1 As shown, including: Step 1: Collect the state quantity of the hydrogen-electric coupling system.

[0026] In a non-limiting preferred embodiment, the state quantities of the hydrogen-electricity coupling system include but are not limited to: grid voltage, grid frequency and heating network temperature; those skilled in the art determine the state quantities of the hydrogen-electricity coupling system to be collected based on the content of the state data.

[0027] Specifically, the corresponding state quantity data is collected through the vector measurement unit, pressure transmitter, temperature sensor and voltage transformer in the hydrogen-electricity coupling system.

[0028] In the embodiment, the hydrogen-electricity coupling system is a hydrogen-electricity coupling micro energy network.

[0029] Step 2: Establish a hydrogen-electricity coupling system model; obtain fault data under various fault scenarios based on the hydrogen-electricity coupling system model to construct an accident set.

[0030] Specifically, step 2 includes: Step 2.1, build a hydrogen-electricity coupling system model; Specifically, the hydrogen-electricity coupling system model includes, but is not limited to: an electric energy subsystem model, a thermal subsystem model and a hydrogen energy subsystem model.

[0031] The power subsystem includes conventional thermal power units and new energy units. The power subsystem model includes: 1) Output constraint model of conventional thermal power units: , In the formula, for t Conventional thermal power units of efforts, and Conventional thermal power units The lower and upper limits of output.

[0032] 2) Output constraint model of new energy units: , In the formula, for t New energy units are always producing power. and They are respectively the lower and upper limits of the output of new energy units.

[0033] The thermal subsystem includes thermal network and hydraulic network. The thermal subsystem model includes: 1) Thermal network model: , In the formula, is the association matrix between nodes and branches in the thermal network; is the association matrix between loops and branches in the thermal network; is the mass flow column vector of the thermal pipeline; is the injected mass flow column vector of the source node and load node in the thermal network; is the pressure loss column vector of the thermal pipeline; is the resistance coefficient column vector of the thermal pipe, which generally depends on the pipe diameter.

[0034] 2) Hydraulic network model: , In the formula, is the column vector of thermal power consumed or supplied to the node; is the specific heat capacity of water; is the column vector of water supply temperature; is the column vector of outlet water temperature; , are the initial end node temperature and the terminal node temperature of the hot water pipe respectively; is the ambient temperature; is the total heat transfer coefficient per unit length of the hot water pipe; is the length of the hot water pipe; are the mass flow column vectors of each hot water pipe flowing into and out of the node respectively; is the heat medium temperature of the node; is the heat medium temperature at the end of each pipe section flowing into the node.

[0035] It is the product of the sum of the mass flow rates of all hot water pipes flowing out of the node and the temperature of the heat medium at the node, and is used to represent the heat provided by the heat medium; It is the sum of the products of the mass flow rate of each hot water pipe flowing into the node and the temperature of the heat medium at the end of each pipe section flowing into the corresponding node. It is used to represent the heat required by each node. When the two are balanced, the energy consumption is the lowest.

[0036] The hydrogen energy subsystem includes an electrolyzer, a fuel cell, and a hydrogen storage tank. The hydrogen energy subsystem model includes: 1) Electrolyzer model: , In the formula, is the actual input power of the electrolyzer; input power to the electrolyzer; is the efficiency of the electrolyzer; is the hydrogen production power of the electrolyzer; is the hydrogen production efficiency of the electrolyzer.

[0037] The electrolyzer model includes power and ramp constraints in the operating state: , In the formula, To indicate that the electrolytic cell is t A 0-1 variable indicating the start and stop status at all times; for t The electrolyzer input power at the moment; , are the lower and upper limits of the electrolytic cell input power respectively.

[0038] 2) Fuel cell model: , In the formula, is the output power of the fuel cell; is the input power of the fuel cell; The efficiency of the fuel cell.

[0039] The fuel cell model includes the operating constraints of the fuel cell: , In the formula, To indicate that the fuel cell t A 0-1 variable indicating the start and stop status at all times; , are the lower and upper limits of the fuel cell input power respectively.

[0040] 3) Hydrogen storage tank model: , In the formula, The energy stored in the hydrogen tank at the initial moment; is the capacity of the hydrogen storage tank; for tThe energy stored in the hydrogen tank at all times; , They are t The charging and discharging power of hydrogen storage tank at all times; , They are the charging and discharging efficiency of the hydrogen storage tank respectively; For the operating period; is the maximum capacity factor of the hydrogen storage tank; , They represent the hydrogen storage tank t A 0-1 variable indicating the state of charge and discharge at any given moment; , They are the maximum values ​​of charging and discharging power of the hydrogen storage tank respectively.

[0041] Step 2.2, based on the hydrogen-electricity coupling system model, obtain the fault data under each fault scenario to construct an accident set.

[0042] In a non-limiting preferred embodiment, the fault scenarios include but are not limited to: failure of key pipeline outlets at energy hubs, failure of key equipment at energy hubs, and failure of imbalance between supply and demand of renewable energy; those skilled in the art can determine the specific fault scenarios based on the probability of occurrence of each fault in the hydrogen-electricity coupling system.

[0043] When obtaining fault data under various fault scenarios based on the hydrogen-electric coupling system model, the randomness of the output of new energy units is comprehensively considered, and the system energy supply capacity is increased and the upper limit of capacity constraints is expanded when analyzing the accident set. Therefore, considering the feasible boundary and safety boundary under the condition of new energy unit access, the output of new energy units is probabilistically modeled, and a joint model of the output of new energy units is formed with the output constraint model of new energy units, thereby updating the electric energy subsystem model.

[0044] In a non-limiting preferred embodiment, the probability distribution model of wind turbine output is established based on the two-parameter Weibull random probability model to satisfy the following relationship: , , In the formula, is the probability density function of the actual wind speed, k w is the shape factor of the fan unit, c w is the scale factor of the fan unit; , V h , V out , V in They are rated wind speed, actual wind speed, cut-out wind speed and cut-in wind speed respectively; Pwmax is the maximum output power of the wind turbine, is the wind turbine output based on probability distribution.

[0045] In a non-limiting preferred embodiment, the probability distribution model of the photovoltaic power generation group output is established based on the Beta probability distribution model to satisfy the following relationship: , , In the formula, is the probability density function of actual illumination, is the shape factor of the photovoltaic module, β p is the scale factor of the PV module; , They are rated light intensity and actual light intensity respectively; is the Gamma function; Indicates the maximum output power of the photovoltaic generator set. Photovoltaic generator output based on probability distribution.

[0046] The present invention also considers the access of new energy to the micro-energy network. Most of the current micro-energy networks do not consider the part of new energy. Based on the randomness of new energy power generation, new energy probability modeling is added to enhance the precision of the model. In addition, the new energy probability modeling associates the randomness of the output of new energy units with the system energy supply capacity and capacity constraints, and realizes that the feasible domain and safety domain based on the hydrogen-electric coupling system model fully consider the randomness of the output of new energy units and have adaptive adjustment capabilities.

[0047] Step 3: Establish joint constraints based on the hydrogen-electricity coupling system model and the accident set, including multi-energy flow balance constraints, feasible domain constraints and safety domain constraints.

[0048] Specifically, step 3 includes: Step 3.1, establishing multi-energy flow balance constraints based on the hydrogen-electricity coupling system model; In a non-limiting preferred embodiment, the multi-energy flow balance constraints include but are not limited to: 1) Power balance constraint, satisfying the following relationship: , In the formula, , They are t Time Node i Active injection power and reactive injection power; N is the number of system nodes; , They are the node admittance matricesi OK j The real and imaginary parts of the columns; for t Time Node i The voltage amplitude of for t Time Node j The voltage amplitude of for t Time Branch ij The phase angle difference; , Node i The lower and upper limits of the voltage amplitude; 2) Hydrogen power balance constraint, satisfying the following relationship: , In the formula, for t The hydrogen production power of the electrolyzer at any moment; , They are t The charging and discharging power of hydrogen storage tank at all times; for t Hydrogen load at the time; for t The input power of the fuel cell at the moment.

[0049] 3) Power balance constraints of the electric-hydrogen coupling system: Since the electric energy subsystem produces hydrogen through the excess electricity of the new energy station to meet the energy demand of the hydrogen energy subsystem, and the hydrogen energy subsystem provides electricity to the electric energy subsystem through hydrogen fuel cells, the electric-hydrogen coupling operation constraints are as follows: The power balance constraints of the electric-hydrogen coupling system include the electric-hydrogen coupling operation constraints and the electric-hydrogen production power constraints, where the electric-hydrogen coupling operation constraints are as follows: , In the formula, for t Conventional thermal power units contribution; is the number of conventional thermal power units; for t New energy units are always producing; for t The output power of the fuel cell at any moment; for t Time Node The electrical load; is the hydrogen load of the system; for t The electrolyzer input power at the moment; is the number of nodes.

[0050] Using the excess electricity from new energy stations to produce hydrogen, the power constraints for hydrogen production are as follows: , In the formula, for Shike New Energy Station Output power; for Shike New Energy Station The load power provided; for Shike New Energy Station The abandoned power; for Shike New Energy Station Input power to the electrolyzer.

[0051] 4) Thermal energy balance constraints: , In the formula, C CHP is the thermoelectric ratio; P CHP,h is the thermal power; P CHP,e For electrical power.

[0052] 5) Flow constraints of hydraulic network: , In the formula, B It is the association matrix of the hydraulic grid, describing the association between the hydraulic loop and the pipe section; w is the branch flow vector; K is the pipe coefficient.

[0053] Step 3.2, based on the hydrogen-electric coupling system model and the multi-energy flow balance constraints, with the goal of maximizing the energy supply range of the hydrogen-electric coupling system under the "N-1" condition, establish the feasible domain constraints of the hydrogen-electric coupling system under the "N-1" condition; Specifically, there are a total of pipeline segments, numbered 1, 2, ..., M ; Using vectors in Euclidean space Indicates the working status of the hydrogen-electric coupling system.

[0054] With the goal of maximizing the energy supply range of the hydrogen-electric coupling system under the "N-1" condition, the feasible domain constraints are established to satisfy the following relationship: , In the formula, For the The state quantity of each pipeline segment, The sum of the state quantities of all pipeline segments is the largest, so that the energy supply range of the hydrogen-electric coupling system under the "N-1" condition is maximized; The equation constraints for the key pipeline segments and key equipment in the feasible domain of the hydrogen-electric coupling system under the “N-1” condition are: Inequality constraints for key pipeline segments and key equipment in the feasible domain of the hydrogen-electric coupling system under the “N-1” condition; , are the lower bound vector and upper bound vector in the inequality constraints respectively.

[0055] Step 3.3, determine the safety domain boundary according to the feasible domain constraints and accident set of the hydrogen-electric coupling system under the "N-1" condition, and establish the safety domain constraints of the hydrogen-electric coupling system under the "N-1" condition based on the safety domain boundary.

[0056] Security domain constraints satisfy the following relationship: , In the formula, It is the security domain constraint; It is the equation constraint condition for the key pipeline segments and key equipment in the safety domain of the hydrogen-electric coupling system under the “N-1” condition; Inequality constraints for key pipeline segments and key equipment in the safety domain of the hydrogen-electric coupling system under the “N-1” condition Under the union of the feasible region and the safety region, the inequality constraint is .

[0057] In the method proposed in the present invention, the feasible domain and safety domain of the hydrogen-electricity coupling system can be obtained according to the joint constraints. The following principles should be observed when determining the feasible domain: the flow capacity of the key pipeline segment is within the allowable operating boundary; the flow capacity of the key pipeline segment of the connecting equipment is not greater than the minimum capacity of the key pipeline segment; the minimum capacity of the key pipeline segment should ensure the energy supply of the connected equipment. The following principles should be observed when determining the safety domain: the accident set includes accident data of key pipeline outlet failures of energy hubs, accident data of key equipment failures of energy hubs, and accident data of renewable energy supply and demand imbalance failures; the randomness of renewable energy is comprehensively considered during calculation, and the system energy supply capacity is increased and the upper limit of the capacity constraint is expanded when analyzing the first accident set.

[0058] Step 4, screen out the state quantities that do not meet the multi-energy flow balance constraint conditions, screen out the state quantities that do not meet the feasible domain constraint conditions of the hydrogen-electric coupling system under the "N-1" condition, and screen out the state quantities that do not meet the safety domain constraint conditions of the hydrogen-electric coupling system under the "N-1" condition as the state quantities to be adjusted.

[0059] Specifically, according to the multi-energy flow balance constraints, the feasible domain constraints of the hydrogen-electric coupling system under the "N-1" condition, and the safety domain constraints, the isolation forest algorithm is used to screen the state quantities, which not only removes the abnormal data in the state quantities, but also extracts the state quantities to be adjusted.

[0060] In a non-limiting preferred embodiment, the method of screening state quantities using the feasible domain constraints and safety domain constraints of the hydrogen-electricity coupling system includes but is not limited to: the isolation forest method.

[0061] In a non-limiting preferred embodiment, abnormal data includes, but is not limited to, distorted data caused by external factors.

[0062] Specifically, the state quantity data set is screened by the isolation forest algorithm, including: Step S1, is a state quantity data set, from X Extract Samples as X Subset of Put it into the root node, where is the number of samples.

[0063] Step S2, specify the dimension and in the subset A random cutting point is generated in Step S3, based on the cutting point p Subset Divide into two subspaces, including: specify the dimension q Less than p The samples are placed in the left subspace, specifying the dimension q Not less than p into the right subspace.

[0064] Step S4, repeating steps S2 and S3 until all trees reach a specified height.

[0065] Step S5, repeating steps S1 to S4, and constructing an isolation forest using the generated r isolated trees.

[0066] For each sample , in its traversal After the number of isolated trees, we will calculate Average height in the forest , normalize the average height of all samples to satisfy the following relationship: , In the formula, is the normalized value of the average height of all samples; To include The path length of the isolated tree in the subset of samples (the path of the isolated tree, Monotonically increasing), used to standardize the recording of external nodes x The path length satisfies the following relationship: , in, , For variables The harmonic function of ξ is Euler's constant; From the root node to the external node x Average height The average value of satisfies the following relationship: , In the formula, is an outlier.

[0067] Furthermore, after the collected state quantities are screened using the feasible domain constraints and safety domain constraints of the hydrogen-electricity coupling system, the state quantities before and after the screening are obtained.

[0068] Step 5, using state quantities other than the state quantities to be adjusted, with the efficiency and safety levels of the feasible domain being the best as the optimization target model; using the optimization target model and the safety domain constraints of the hydrogen-electric coupling system under the "N-1" condition to form an optimization control model of the hydrogen-electric coupling system; based on the optimization control model, optimizing the state quantities to be adjusted, obtaining the optimized values ​​of the state quantities that simultaneously satisfy the optimization target model and the safety domain constraints of the hydrogen-electric coupling system under the "N-1" condition, and controlling the hydrogen-electric coupling system according to the optimized values ​​of the state quantities.

[0069] Specifically, step 5 includes: Step 5.1, using the state quantities other than the state quantities to be adjusted, taking the efficiency and safety of the feasible domain as the optimal target model; In the embodiment, the hydrogen-electric coupling system has pipeline segments, numbered 1, 2, ..., M ; Using vectors in Euclidean space Indicates the working status of the hydrogen-electric coupling system. For the The state quantity of each pipeline segment, .

[0070] Among them, the efficiency degree function of the feasible domain satisfies the following relationship: , In the formula, is the efficiency function of the feasible domain, is the total number of state quantities other than the state quantity to be adjusted. In the embodiment, after the state quantity is screened, the hydrogen-electric coupling system The state quantity of each pipeline segment is the object, The state quantity to be adjusted is The state quantity of each pipeline segment, For the The maximum value of the state quantity of a pipeline segment under the constraints of the feasible domain.

[0071] The safety function of the feasible region satisfies the following relationship: , In the formula, is the safety function of the feasible domain; is the safety distance from the safe working point c in the safety domain to the upper safety boundary; is the safety distance from the minimum energy supply capacity working point z in the safety domain to the upper safety boundary, where the minimum energy supply capacity refers to the minimum load carrying capacity of the system. is the coupling coefficient of the hydrogen-electric coupling system, .

[0072] In the embodiment, the proportion of the number of coupled devices in the total number of devices in the hydrogen-electricity coupling system is used as the coupling degree coefficient. When the proportion of the number of coupled devices in the total number of devices in the hydrogen-electricity coupling system is less than or equal to 20%, , when the number of coupled devices in the hydrogen-electricity coupling system accounts for more than 20% and less than or equal to 50% of the total number of devices .

[0073] Step 5.2, using the optimization target model and the safety domain constraints of the hydrogen-electric coupling system under the "N-1" condition, the optimization control model of the hydrogen-electric coupling system is constructed to satisfy the following relationship: , In the formula, The optimization target model that represents the optimal efficiency function and safety function of the feasible domain is represents the safety domain constraint of the hydrogen-electric coupling system under the “N-1” condition, where: O is the safety domain constraint of the hydrogen-electric coupling system under the “N-1” condition, L is the state quantity set of each pipeline segment.

[0074] Step 5.3, based on the optimization control model, optimize the state quantity to be adjusted, obtain the optimized value of the state quantity that satisfies both the optimization target model and the safety domain constraint of the hydrogen-electric coupling system under the "N-1" condition, and control the hydrogen-electric coupling system according to the optimized value of the state quantity.

[0075] In a non-limiting preferred embodiment, the state quantity to be adjusted includes but is not limited to: the load power on the pipeline segment.

[0076] Specifically, step 5.3 includes: Step 5.3.1, based on the three-dimensional image method, draw the three-dimensional safety domain corresponding to the safety domain constraint conditions of the hydrogen-electric coupling system under the "N-1" condition; Specifically, the power pipeline segment , thermal pipeline segment , Hydrogen pipeline segment The three-dimensional safety domain composed of energy supply loads is as follows Figure 2 In the three-dimensional safety domain, when the operating state of the micro energy grid is within the region, it means that the system is in a safe and reliable operating state. The surface BCEF is the common safety upper bound, the surface ABCD is the common safety lower bound, and CDF is the power pipeline segment The upper limit of the energy supply load, AEB is the power pipeline segment The safety lower limit of the upper energy load, ADEF is the thermal pipeline segment The safe lower limit of the upper energy load. Figure 2 For the thermal pipeline segment Energy supply load and hydrogen pipeline segment The energy load is converted into numerical value and compared with the power pipeline segment The upper energy load uses the same unit, MVA.

[0077] Step 5.3.2, using the Manhattan distance method, calculate the safety distance between the state quantity to be adjusted and the boundary of the three-dimensional safety zone as the adjustment distance of the state quantity to be adjusted, satisfying the following relationship: , In the formula, is the adjustment amount of the state quantity, Indicates The upper bound of the state quantity of a pipeline segment.

[0078] exist If the safety domain boundary is exceeded, it is necessary to first make safety adjustments, calculate the Mandaton distance between the operating point and the safety boundary, and then optimize the performance function. The operation adjustment process is as follows: Figure 3 As shown in the figure, the calculated data itself is in an unsafe range. Now, through the Mandarin distance, we can get the distance that can be adjusted back to the safe range, and then make it return to the safe range. In the process of adjustment, the efficiency function should be considered first to ensure that the adjusted operating point can achieve higher operating efficiency and safety.

[0079] In a non-limiting preferred embodiment, the adjustment data is returned through the adjustment strategy; based on whether the output status of the pipeline exceeds the safety boundary after the adjustment data is adjusted, the pipeline is returned to the safety domain and continues to participate in the electricity-hydrogen coupling to operate as an energy grid.

[0080] Furthermore, the screened state quantities are used to optimize the control of the hydrogen-electric coupling system based on the optimization control model of the hydrogen-electric coupling system, and then the state quantities before and after are adjusted.

[0081] In the method proposed in the present invention, the state quantities screened out are parameters that do not meet some conditions, and these parameters are often the predictive data sources of system operation safety problems. Simply deleting these parameters without using these parameters to predict the safety and stability of the system will lead to the risk of system control being unable to avoid safety hazards. In the method proposed in the present invention, various constraints are gradually constructed in depth, from multi-energy flow balance constraints to feasible domain constraints of hydrogen-electric coupling systems under "N-1" conditions to safety domain constraints of hydrogen-electric coupling systems under "N-1" conditions, and the state quantities that do not meet any of the above three conditions are screened out as the state quantities to be adjusted, which is the basis for realizing the prediction of system safety and stability risks.

[0082] After screening out the state quantities to be adjusted, the remaining state quantities need to achieve both efficient and safe operation of the system in the feasible domain. Therefore, the state quantities other than the state quantities to be adjusted are used to optimize the efficiency degree function and the safety degree function of the feasible domain as the optimization target model; then the optimization target model and the safety domain constraints of the hydrogen-electric coupling system under the "N-1" condition are used to form an optimization control model of the hydrogen-electric coupling system. The optimization control model is a preventive control model that can adjust the system from an unsafe state to a safe state, and then from a safe state to a safe and efficient state. Therefore, based on the optimization control model, the state quantities to be adjusted are optimized, so that the unsafe working point can be adjusted to a safe and efficient working point along the shortest path.

[0083] Moreover, the method for establishing the safety domain constraint conditions of the hydrogen-electric coupling system under the "N-1" condition in the present invention enables the safety domain to be intuitively displayed through a three-dimensional image, which helps to adjust the unsafe working point to a safe and efficient working point along the shortest path.

Claims

1. A method for screening and dynamically optimizing the control of safety state quantities for a hydrogen-electricity coupling system, the hydrogen-electricity coupling system comprising: The electric energy subsystem, the thermal subsystem and the hydrogen energy subsystem are characterized by comprising: Collect the state quantity of the hydrogen-electric coupling system, establish a hydrogen-electric coupling system model, and obtain fault data under various fault scenarios based on the hydrogen-electric coupling system model to construct an accident set; Establish joint constraints based on the hydrogen-electricity coupling system model and the accident set, including: establishing multi-energy flow balance constraints based on the hydrogen-electricity coupling system model; establishing feasible domain constraints of the hydrogen-electricity coupling system under "N-1" conditions based on the hydrogen-electricity coupling system model and the multi-energy flow balance constraints, with the goal of maximizing the energy supply range of the hydrogen-electricity coupling system under "N-1" conditions; establishing safety domain constraints of the hydrogen-electricity coupling system under "N-1" conditions based on the feasible domain constraints and accident set of the hydrogen-electricity coupling system under "N-1" conditions; Screen out the state quantities that do not satisfy the multi-energy flow balance constraint conditions, screen out the state quantities that do not satisfy the feasible domain constraint conditions of the hydrogen-electric coupling system under the "N-1" condition, and screen out the state quantities that do not satisfy the safety domain constraint conditions of the hydrogen-electric coupling system under the "N-1" condition as the state quantities to be adjusted; Using the state quantities other than the state quantities to be adjusted, the efficiency and safety of the feasible domain are both optimal as the optimization target model; the efficiency function of the feasible domain satisfies the following relationship: ， In the formula, is the efficiency function of the feasible domain, is the total number of state quantities other than the state quantity to be adjusted. In the embodiment, after the state quantity is screened, the hydrogen-electric coupling system The state quantity of each pipeline segment is the object, The state quantity to be adjusted is The state quantity of each pipeline segment, For the The maximum value of the state quantity of each pipeline segment under the constraints of the feasible region; The safety function of the feasible region satisfies the following relationship: ， In the formula, is the safety function of the feasible domain; is the safety distance from the safe working point c in the safety domain to the upper safety boundary; is the safety distance from the minimum energy supply capacity working point z in the safety domain to the upper safety boundary, where the minimum energy supply capacity refers to the minimum load carrying capacity of the system. is the coupling coefficient of the hydrogen-electric coupling system, ; An optimization control model of the hydrogen-electricity coupling system is constructed by using the optimization target model and the safety domain constraints of the hydrogen-electricity coupling system under the "N-1" condition. Based on the optimization control model, the state quantity to be adjusted is optimized to obtain the optimized value of the state quantity that satisfies both the optimization target model and the safety domain constraints of the hydrogen-electricity coupling system under the "N-1" condition. The hydrogen-electricity coupling system is controlled according to the optimized value of the state quantity.

2. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 1 is characterized in that: The state quantities of the hydrogen-electricity coupling system include: grid voltage, grid frequency and heating network temperature.

3. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 1 is characterized in that: The power subsystem includes conventional thermal power units and new energy units. The power subsystem model includes: 1) Output constraint model of conventional thermal power units: ， In the formula, for t Conventional thermal power units of efforts, and Conventional thermal power units The lower and upper limits of output; 2) Output constraint model of new energy units: ， In the formula, for t New energy units are always producing power. and They are respectively the lower and upper limits of the output of new energy units.

4. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 1 is characterized in that: The thermal subsystem includes thermal network and hydraulic network. The thermal subsystem model includes: 1) Thermal network model: ， In the formula, is the association matrix between nodes and branches in the thermal network; is the association matrix between loops and branches in the thermal network; is the mass flow column vector of the thermal pipeline; is the injected mass flow column vector of the source node and load node in the thermal network; is the pressure loss column vector of the thermal pipeline; is the resistance coefficient column vector of the thermal pipe, which depends on the pipe diameter; 2) Hydraulic network model: ， In the formula, is the column vector of thermal power consumed or supplied to the node; is the specific heat capacity of water; is the column vector of water supply temperature; is the column vector of outlet water temperature; , are the initial end node temperature and the terminal node temperature of the hot water pipe respectively; is the ambient temperature; is the total heat transfer coefficient per unit length of the hot water pipe; is the length of the hot water pipe; are the mass flow column vectors of each hot water pipe flowing into and out of the node respectively; is the heat medium temperature of the node; is the heat medium temperature at the end of each pipe section flowing into the node.

5. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 1 is characterized in that: The hydrogen energy subsystem includes an electrolyzer, a fuel cell, and a hydrogen storage tank. The hydrogen energy subsystem model includes: 1) Electrolyzer model: ， In the formula, is the actual input power of the electrolyzer; input power to the electrolyzer; is the efficiency of the electrolyzer; is the hydrogen production power of the electrolyzer; is the hydrogen production efficiency of the electrolyzer; The electrolyzer model includes power and ramp constraints in the operating state: ， In the formula, To indicate that the electrolytic cell is t A 0-1 variable indicating the start and stop status at all times; for t The electrolyzer input power at the moment; , are the lower and upper limits of the electrolyzer input power, respectively; 2) Fuel cell model: ， In the formula, is the output power of the fuel cell; is the input power of the fuel cell; for the efficiency of the fuel cell; The fuel cell model also includes the operating constraints of the fuel cell: ， In the formula, To indicate that the fuel cell t A 0-1 variable indicating the start and stop status at all times; , are the lower and upper limits of the fuel cell input power, respectively; 3) Hydrogen storage tank model: ， In the formula, The energy stored in the hydrogen tank at the initial moment; is the capacity of the hydrogen storage tank; for t The energy stored in the hydrogen tank at all times; , They are t The charging and discharging power of hydrogen storage tank at all times; , They are the charging and discharging efficiency of the hydrogen storage tank respectively; For the operating period; is the maximum capacity factor of the hydrogen storage tank; , They represent the hydrogen storage tank t A 0-1 variable indicating the state of charge and discharge at any given moment; , They are the maximum values ​​of charging and discharging power of the hydrogen storage tank respectively.

6. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 1 is characterized in that: Failure scenarios include: failure of key pipeline outlets at energy hubs, failure of key equipment at energy hubs, and imbalance between supply and demand of renewable energy.

7. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 6 is characterized in that: Based on the hydrogen-electric coupling system model, the fault data under each fault scenario is obtained to construct an accident set, including: Probabilistic modeling of the output of new energy units is carried out, and a joint model of the output of new energy units is formed with the output constraint model of new energy units, thereby updating the electric energy subsystem model and obtaining an updated hydrogen-electric coupling system model; Among them, the probability modeling of the output of new energy units includes: establishing a probability distribution model of wind turbine output based on a two-parameter Weibull random probability model, and establishing a probability distribution model of photovoltaic generator output based on a Beta probability distribution model; Based on the updated hydrogen-electric coupling system model, the fault data under each fault scenario is obtained to construct an accident set.

8. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 7 is characterized in that: The probability distribution model of wind turbine output is established based on the two-parameter Weibull random probability model to satisfy the following relationship: ， In the formula, is the probability density function of the actual wind speed, k w is the shape factor of the fan unit, c w is the scale factor of the fan unit; , V h , V out , V in They are rated wind speed, actual wind speed, cut-out wind speed and cut-in wind speed respectively; P wmax is the maximum output power of the wind turbine, The output of wind turbines based on probability distribution; The probability distribution model of photovoltaic generator output based on the Beta probability distribution model satisfies the following relationship: ， In the formula, is the probability density function of actual illumination, is the shape factor of the photovoltaic module, β p is the scale factor of the PV module; , They are rated light intensity and actual light intensity respectively; is the Gamma function; Indicates the maximum output power of the photovoltaic generator set. Photovoltaic generator output based on probability distribution.

9. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 1, characterized in that: The multi-energy flow balance constraints include: 1) Power balance constraint, satisfying the following relationship: ， In the formula, , They are t Time Node i Active injection power and reactive injection power; N is the number of system nodes; , They are the node admittance matrices i OK j The real and imaginary parts of the columns; for t Time Node i The voltage amplitude of for t Time Node j The voltage amplitude of for t Time Branch ij The phase angle difference; , Node i The lower and upper limits of the voltage amplitude; 2) Hydrogen power balance constraint, satisfying the following relationship: ， In the formula, for t The hydrogen production power of the electrolyzer at any moment; , They are t The charging and discharging power of hydrogen storage tank at all times; for t Hydrogen load at the time; for t The input power of the fuel cell at any moment; 3) Power balance constraints of the electric-hydrogen coupling system, including electric-hydrogen coupling operation constraints and electric hydrogen production power constraints. The electric-hydrogen coupling operation constraints are as follows: ， In the formula, for t Conventional thermal power units contribution; is the number of conventional thermal power units; for t New energy units are always producing; for t The output power of the fuel cell at any moment; for t Time Node The electrical load; is the hydrogen load of the system; for t The electrolyzer input power at the moment; is the number of nodes; The power constraints for hydrogen production are as follows: ， In the formula, for Shike New Energy Station Output power; for Shike New Energy Station The load power provided; for Shike New Energy Station The abandoned power; for Shike New Energy Station Input power of the inner electrolyzer; 4) Thermal energy balance constraint, satisfying the following relationship: ， In the formula, C CHP is the thermoelectric ratio; P CHP,h is the thermal power; P CHP,e is the electrical power; 5) Flow constraints of hydraulic network: ， In the formula, B is the correlation matrix of the hydraulic grid; w is the branch flow vector; K is the pipe coefficient.

10. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 1, characterized in that: The feasible domain constraints of the hydrogen-electric coupling system under the "N-1" condition satisfy the following relationship: ， In the formula, For the The state quantity of each pipeline segment, The sum of the state quantities of all pipeline segments is the largest; are the equation constraints for the key pipeline segments and key equipment in the feasible domain of the hydrogen-electric coupling system under the "N-1" condition. The key pipeline sections and key equipment constraints of the feasible domain of the hydrogen-electric coupling system under "N-1" conditions; , are the lower and upper bound vectors of the inequality constraints respectively; , M is the number of pipeline segments.

11. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 10, characterized in that: The safety domain constraints of the hydrogen-electric coupling system under the "N-1" condition satisfy the following relationship: ， In the formula, It is the security domain constraint; The equation constraints for the key pipeline segments and key equipment in the safety domain of the hydrogen-electric coupling system under the "N-1" condition; These are the inequality constraints for key pipeline segments and key equipment in the safety domain of the hydrogen-electric coupling system under the "N-1" condition.

12. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 11, characterized in that: According to the multi-energy flow balance constraints, the feasible domain constraints of the hydrogen-electric coupling system under the "N-1" condition, and the safety domain constraints, the isolation forest algorithm is used to screen the state quantities, which not only removes the abnormal data in the state quantities, but also extracts the state quantities to be adjusted.

13. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 1, characterized in that: The proportion of the number of coupled devices in the total number of devices in the hydrogen-electricity coupling system is taken as the coupling degree coefficient. When the proportion of the number of coupled devices in the total number of devices in the hydrogen-electricity coupling system is less than or equal to 20%, , when the number of coupled devices in the hydrogen-electricity coupling system accounts for more than 20% and less than or equal to 50% of the total number of devices .

14. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 1, characterized in that: The optimization target model and the safety domain constraints of the hydrogen-electric coupling system under the "N-1" condition are used to form the optimization control model of the hydrogen-electric coupling system, which satisfies the following relationship: ， In the formula, The optimization target model that represents the optimal efficiency function and safety function of the feasible domain is represents the safety domain constraint of the hydrogen-electric coupling system under the "N-1" condition, where: O is the safety domain constraint of the hydrogen-electric coupling system under the "N-1" condition, L is the state quantity set of each pipeline segment.

15. The method for screening and dynamically optimizing the safety state quantity of a hydrogen-electricity coupling system according to claim 14, characterized in that: Based on the optimization control model, the state quantity to be adjusted is optimized, including: Based on the three-dimensional image method, the three-dimensional safety domain corresponding to the safety domain constraints of the hydrogen-electric coupling system under the "N-1" condition is drawn; The Manhattan distance method is used to calculate the safety distance between the state quantity to be adjusted and the boundary of the three-dimensional safety zone as the adjustment distance of the state quantity to be adjusted.

Images