Boundary determination method and device of steady-state security domain, equipment, medium and product

Abstract

The invention discloses a boundary determination method and device of a steady-state security domain, equipment, a medium and a product. The method comprises the following steps: establishing a comprehensive energy station model under an electric heating system, and obtaining related parameter information of the comprehensive energy station model and load information of each energy station; according to the load information and the related parameter information, determining operation constraint information satisfying the operation of the electric heating system; according to the operation constraint information and a preset quotient gradient system model, constructing a steady-state security domain model of the comprehensive energy station model; and solving the steady-state security domain model based on a numerical trajectory unification algorithm to obtain the boundary of the steady-state security domain of the comprehensive energy station model. A steady-state safety domain model is constructed through a preset quotient gradient system model, an accurate model of a multi-energy station steady-state safety region is completely represented, solving is carried out through a numerical trajectory unification algorithm, and a boundary is determined. The calculation of a large number of meaningless points in the safety area is avoided, the calculation efficiency is improved, and the accuracy of boundary determination is ensured.

Description

Technical Field

[0001] The present invention relates to the technical field of integrated energy systems, and in particular to a method, device, equipment, medium and product for determining the boundary of a steady-state safety domain. Background Art

[0002] There is a strong correlation between the energy sources in the integrated energy system, and the interaction between them will affect the stability and efficiency of the system operation. In order to ensure that the system operation can meet the energy needs of users, it is very necessary to explore the safety model and safety analysis method of regional integrated energy stations. The energy safety of the integrated energy system is based on the production energy demand and the joint planning and coordinated operation of the integrated energy station, and the key equipment of the energy hub and the distribution network operation safety constraints and multi-energy coupling component constraints are established under the interconnection of multiple energy stations. The interdependence between different energy supply systems may lead to more safety problems due to disturbance conduction.

[0003] There are many methods for estimating the steady-state safety region of power systems. The "steady-state safety region" was proposed in the 1970s and is defined as a set of realistic power injections that satisfy the safety constraints imposed by the power flow equations and equipment operating limitations. In large power systems, an approximate safety region represented by a hyperplane is usually developed.

[0004] However, when constructing the steady-state safety region, some approximations and simplifications are assumed, such as assuming that the steady-state safety region is convex and connected, or simplifications are made in the power flow model, which leads to incomplete representation of the approximate steady-state safety region and inaccurate boundary determination. Summary of the invention

[0005] The present invention provides a method, device, equipment, medium and product for determining the boundary of a steady-state security domain, so as to ensure the accuracy of boundary determination.

[0006] According to a first aspect of the present invention, a method for determining a boundary of a stable security domain is provided, characterized by comprising:

[0007] Establishing an integrated energy station model under the electric heating system, and obtaining relevant parameter information of the integrated energy station model and load information of each energy station;

[0008] Determining, based on the load information and the relevant parameter information, operation constraint information that satisfies the operation of the electric heating system;

[0009] Constructing a steady-state safety domain model of the integrated energy station model according to the operation constraint information and the preset quotient gradient system model;

[0010] The steady-state safety domain model is solved based on the unified numerical trajectory algorithm to obtain the boundary of the steady-state safety domain of the integrated energy station model.

[0011] According to a second aspect of the present invention, there is provided a device for determining a boundary of a stable security domain, comprising:

[0012] An information acquisition module is used to establish a comprehensive energy station model under the electric heating system, and obtain relevant parameter information of the comprehensive energy station model and load information of each energy station;

[0013] An information determination module, used to determine operation constraint information satisfying the operation of the electric heating system according to the load information and the relevant parameter information;

[0014] A model building module, used to build a steady-state safety domain model of the integrated energy station model according to the operation constraint information and a preset quotient gradient system model;

[0015] The boundary determination module is used to solve the steady-state safety domain model based on the numerical trajectory unified algorithm to obtain the boundary of the steady-state safety domain of the integrated energy station model.

[0016] According to a third aspect of the present invention, there is provided an electronic device, the electronic device comprising:

[0017] at least one processor; and

[0018] a memory communicatively connected to the at least one processor; wherein,

[0019] The memory stores a computer program that can be executed by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for determining the boundary of the steady-state security domain described in any embodiment of the present invention.

[0020] According to a fourth aspect of the present invention, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the method for determining the boundary of a steady-state security domain described in any embodiment of the present invention when executed.

[0021] According to another aspect of the present invention, an embodiment of the present invention further provides a computer program product, the computer program product includes a computer program, and when the computer program is executed by a processor, it implements the method for determining the boundary of a steady-state security domain of any embodiment of the present invention.

[0022] The technical solution of the embodiment of the present invention is to establish an integrated energy station model under the electric heating system to obtain relevant parameter information of the integrated energy station model and load information of each energy station; determine the operation constraint information that satisfies the operation of the electric heating system based on the load information and relevant parameter information; construct a steady-state safety domain model of the integrated energy station model based on the operation constraint information and the preset quotient gradient system model; solve the steady-state safety domain model based on the numerical trajectory unified algorithm to obtain the boundary of the steady-state safety domain of the integrated energy station model. The steady-state safety domain model is constructed by the preset quotient gradient system model, and the accurate model of the steady-state safety area of ​​multiple energy stations is completely characterized. The boundary is determined by solving the boundary through the numerical trajectory unified algorithm. The calculation of a large number of meaningless points in the safety area is avoided, the calculation efficiency is improved, and the accuracy of the boundary determination is guaranteed.

[0023] It should be understood that the contents described in this section are not intended to identify the key or important features of the embodiments of the present invention, nor are they intended to limit the scope of the present invention. Other features of the present invention will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0025] Figure 1 is a flow chart of a method for determining the boundary of a steady-state security domain provided according to Embodiment 1 of the present invention;

[0026] Figure 2 is a flow chart of a method for determining the boundary of a steady-state security domain provided according to Embodiment 2 of the present invention;

[0027] Figure 3 This is an example topology diagram of a comprehensive energy station in a method for determining the boundary of a steady-state security zone provided in Embodiment 2 of the present invention;

[0028] Figure 4 This is a schematic diagram of boundary calculation of a method for determining the boundary of a steady-state security domain provided according to Embodiment 2 of the present invention;

[0029] Figure 5 2 is a schematic diagram of a structure of a device for determining a boundary of a steady-state security domain provided according to a third embodiment of the present invention;

[0030] Figure 6 It is a schematic diagram of the structure of an electronic device implementing an embodiment of the present invention. DETAILED DESCRIPTION

[0031] In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present invention.

[0032] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0033] Embodiment 1

[0034] Figure 1 A flowchart of a method for determining the boundary of a steady-state security domain is provided for the first embodiment of the present invention. This embodiment is applicable to the determination of the boundary of a steady-state security domain of multiple integrated energy stations. The method can be executed by a device for determining the boundary of a steady-state security domain. The device for determining the boundary of a steady-state security domain can be implemented in the form of hardware and / or software. The device for determining the boundary of a steady-state security domain can be configured in an electronic device. Figure 1 As shown, the method includes:

[0035] S110, establishing a comprehensive energy station model under the electric heating system, and obtaining relevant parameter information of the comprehensive energy station model and load information of each energy station.

[0036] In this embodiment, the electric heating system can be understood as a system that includes the integration of electricity and heat. The integrated energy station model can be understood as a system model composed of systems and nodes that integrate electric and thermal energy forms. The energy station can be understood as the main body that provides energy management and supply. The load information can be understood as the load demand of the energy station. The relevant parameter information can be understood as the relevant technical parameters that constitute the integrated energy station model, such as node load (including active and reactive loads), active and reactive output of photovoltaics, network topology, resistance and reactance and other related technical parameters.

[0037] Specifically, relevant personnel can input the integrated energy station information under the actual electric heating system, and the processor can obtain the integrated energy station information and use it to establish an integrated energy station model under the electric heating system, and obtain the load information of each energy station in the integrated energy station model, as well as the relevant parameter information under each node in the integrated energy station model.

[0038] S120. Determine operation constraint information that satisfies the operation of the electric heating system according to load information and related parameter information.

[0039] In this embodiment, the operation constraint information can be understood as the operation constraint conditions for ensuring the safe operation of the electric heating system.

[0040] Specifically, the processor can determine the access nodes and node parameters of each access node in the integrated energy station model based on relevant parameter information, determine the operation constraints corresponding to the nodes through the node parameters, determine the operation constraints corresponding to the equipment based on the load information, and then obtain the operation constraint information that meets the different dimensions of the operation of the electric heating system.

[0041] S130. Construct a steady-state safety domain model of the integrated energy station model based on the operation constraint information and the preset commercial gradient system model.

[0042] In this embodiment, the preset quotient gradient system model (QGS) can be understood as a mathematical model for analyzing and solving the optimal power flow problem in the power system and calculating the unstable equilibrium point. The steady-state security domain model can be understood as a mathematical model for analyzing and ensuring that the energy system remains stable and safe under various operating conditions.

[0043] Specifically, the processor can determine the operational constraints corresponding to the power system and the thermal system respectively based on the operational constraints of the integrated energy station model, determine the representation equation of the steady-state safety domain of the integrated energy station model, and based on the operational constraints of the two systems, explore the complete representation of the steady-state safety domain through the quotient gradient system (QGS) in the nonlinear power system to construct a steady-state safety domain model of the integrated energy station model.

[0044] S140. Solve the steady-state safety domain model based on the unified numerical trajectory algorithm to obtain the boundary of the steady-state safety domain of the integrated energy station model.

[0045] In this embodiment, the Trajectory-Unified Method (TJU) can be understood as an algorithm for calculating the feasible region of a general large-scale optimal power flow problem. The boundary can be understood as a boundary value to ensure the stable operation of the integrated energy station model.

[0046] Specifically, the processor can construct a QGS based on the constraints in the steady-state safety domain model, and then determine a feasible solution by searching along the path. By changing the search direction, the complete boundary of the steady-state safety area is calculated to obtain the boundary of the steady-state safety domain of the integrated energy station model.

[0047] The technical solution of the embodiment of the present invention is to establish an integrated energy station model under the electric heating system to obtain relevant parameter information of the integrated energy station model and load information of each energy station; determine the operation constraint information that satisfies the operation of the electric heating system based on the load information and relevant parameter information; construct a steady-state safety domain model of the integrated energy station model based on the operation constraint information and the preset quotient gradient system model; solve the steady-state safety domain model based on the numerical trajectory unified algorithm to obtain the boundary of the steady-state safety domain of the integrated energy station model. The steady-state safety domain model is constructed by the preset quotient gradient system model, and the accurate model of the steady-state safety area of ​​multiple energy stations is completely characterized. The boundary is determined by solving the boundary through the numerical trajectory unified algorithm. The calculation of a large number of meaningless points in the safety area is avoided, the calculation efficiency is improved, and the accuracy of the boundary determination is guaranteed.

[0048] Embodiment 2

[0049] Figure 2 This is a flow chart of a method for determining the boundary of a stable security domain provided in the second embodiment of the present invention. This embodiment is a further refinement of the above embodiment. Figure 2 As shown, the method includes:

[0050] S201. Establish a comprehensive energy station model under the electric heating system, and obtain relevant parameter information of the comprehensive energy station model and load information of each energy station.

[0051] S202. Determine the number of access nodes, transmission lines, energy stations, heating units and power generation equipment in the integrated energy station model according to the network topology in the relevant parameter information.

[0052] In this embodiment, the network topology can be understood as a physical layout used to describe the connection between various devices. Access nodes can be understood as physical devices or logical entities used to connect energy stations, energy supply pipelines, and load centers. Transmission lines can be understood as media for energy transmission. Heating units can be understood as equipment that provides thermal energy. Information on the number of power generation equipment can be understood as the number of equipment used for power generation, including the number of different types of power generation equipment, such as the number of thermal load nodes, the number of photovoltaics, and the number of branches.

[0053] Specifically, the processor can classify the types and quantities of equipment included in the integrated energy station model based on the network topology in the relevant parameter information, and determine information such as the corresponding quantities of access nodes, transmission lines, energy stations, heating units and different power generation equipment in the integrated energy station model.

[0054] For example, a specific topological structure diagram is used to show: Figure 3 This is an example diagram of a topology of an integrated energy station in a method for determining the boundary of a steady-state security zone provided in Embodiment 2 of the present invention, such as Figure 3 As shown, the integrated energy station includes 33 nodes, 1 heating unit, and 3 energy stations. Energy station a provides electrical load at access node 6, and energy stations b and c provide electrical load and thermal load. Energy station b is at access node 24, and energy station c is at access node 29.

[0055] S203: Determine node constraints related to the access nodes according to the quantity information of the power generation equipment, the node attribute information of each access node in the related parameter information, and the equivalent parameter information of each transmission line.

[0056] In this embodiment, the node attribute information can be understood as reflecting the output attributes of the access node, for example, it can include photovoltaic active output, photovoltaic reactive output, active load, reactive load, node voltage angle and node voltage amplitude, etc. The node constraint condition can be understood as the condition constraining the relevant attributes of the access node.

[0057] Specifically, the processor may determine the node constraint conditions related to the access node according to the quantity information of the power generation equipment, the node attribute information of each access node in the relevant parameter information, and the equivalent parameter information of each transmission line through the pre-set maximum and minimum values.

[0058] S204. Determine output power constraints of the heating unit according to the quantity information of power generation equipment and the heat load demand of each energy station in the load information.

[0059] In this embodiment, the heat load requirement can be understood as: The output power constraint can be understood as:

[0060] Specifically, the processor may accumulate the heat load demand of the energy station in the load information based on the number of heat load nodes in the power generation equipment quantity information to determine the output power constraint condition of the heating unit.

[0061] Exemplarily, the node constraint condition and the output power constraint condition can be expressed by the following formula:

[0062]

[0063] Where P Gi , QGi are the photovoltaic active and reactive outputs of access node i respectively; P Li , Q Li are the active and reactive loads of access node i respectively; P m is the output power of the heating unit, P mLi is the heat load demand of each energy station, where G and B are the equivalent conductance and equivalent magnetic resistance of the transmission line. θ is the node voltage angle, V is the node voltage amplitude, and S ij is the branch flow through branch ij. e is the number of heat load nodes, N PG is the number of photovoltaic nodes, N L Number of branches. They represent the minimum and maximum values ​​of the photovoltaic active output of the access node i, They represent the minimum and maximum values ​​of the photovoltaic reactive output of the access node i, represents the maximum value of the branch flow, They represent the minimum and maximum values ​​of the access voltage amplitude of the access node i, respectively. The minimum and maximum values ​​can be set according to requirements.

[0064] S205. Determine the operation constraint information that satisfies the operation of the electric heating system by using the output power constraint condition and the constraint conditions of each node.

[0065] Specifically, the processor may use the output power constraint condition and all node constraint conditions together as operation constraint information that satisfies the operation of the electric heating system.

[0066] S206. Based on the operation constraint information, the voltage angle and voltage amplitude of the active power-reactive power PQ node in the integrated energy station model are used as state variable information.

[0067] In this embodiment, the active and reactive power node (PQ) can be understood as a node in the power system where the active power P and reactive power Q are given, and the node voltage and phase are to be determined. The state variable information can be understood as control variables including the power subsystem and the heating subsystem.

[0068] Specifically, based on the operation constraint information, the processor can divide the problem variables of the integrated energy station model into two groups: state variables that depend on the control variables, and control variables that can be adjusted in practice. The processor can set the state variables to the voltage angle and voltage amplitude of the PQ node.

[0069] S207. Based on the operation constraint information, the actual power generation and actual voltage value of the photovoltaic node in the integrated energy station model are used as control variable information.

[0070] In this embodiment, the photovoltaic node can be understood as an access node for photovoltaic power generation. The control variable information includes control variables of the power subsystem and the heating subsystem.

[0071] Specifically, based on the operation constraint information, the processor can set the control variables to the actual power generation and voltage of the photovoltaic nodes in the integrated energy station model.

[0072] S208. Determine the power constraint information of the power subsystem and the heating constraint information of the heating subsystem in the electric heating system according to the state variable information and the control variable information.

[0073] In this embodiment, the power subsystem can be understood as the partial system of the electric heating system that only includes power. The heating subsystem can be understood as the partial system of the electric heating system that only includes heating. The power constraint information can be understood as constraint information related to power, which can include two forms of equality constraints and inequality constraints. The heating constraint information can be understood as constraint information related to heating, which can include two forms of equality constraints and inequality constraints.

[0074] Specifically, the processor may determine the power constraint information of the power subsystem and the heating constraint information of the heating subsystem in the electric heating system through the state variable information and the control variable information.

[0075] Exemplarily, it can be expressed by the following formula:

[0076]

[0077] In the formula, and are the constraint equations for the power subsystem and the heating subsystem, and denote the constraint inequalities of the power system and the heating system respectively, u e represents the control variable of the power subsystem, u g represents the control variable of the heating subsystem, y e represents the state variable of the power subsystem, y g Represents the state variables of the heating subsystem.

[0078] S209. The power inequality constraints in the power constraint information and the heat supply inequality constraints in the heat supply constraint information are respectively converted by presetting slack variable information to obtain slack constraint conditions of the integrated energy station model.

[0079] In this embodiment, the preset slack variable information can be understood as preset variables used for equation conversion, which can include slack variables of the power subsystem and slack variables of the heating subsystem. The slack constraint condition can be understood as the converted slack constraint condition.

[0080] Specifically, the processor can convert the power inequality constraints in the power constraint information and the heating inequality constraints in the heating constraint information respectively by presetting the slack variable information, so as to obtain the slack constraint conditions of the integrated energy station model.

[0081] Exemplarily, the relaxed constraint condition H0(x) after conversion may be:

[0082]

[0083] x=(u e ,u g ,y e ,y g ,S e ,S g )

[0084] In the formula, H e (x) represents the power constraint information of the power subsystem, H g (x) represents the heating constraint information of the heating subsystem. and are the preset slack variables for the power subsystem and the heating subsystem respectively. represents the number of equations after the introduction of slack variables.

[0085] S210. Determine the steady-state safety domain model of the integrated energy station model based on the relaxed constraint conditions and the preset quotient gradient system model.

[0086] Specifically, the processor can determine the representation of the steady-state safety domain according to the relaxed constraints, and determine the steady-state safety domain model of the integrated energy station model by presetting the quotient gradient system model.

[0087] Further, based on the above embodiment, the steps of determining the steady-state safety domain model of the integrated energy station model according to the relaxed constraint conditions and the preset quotient gradient system model can be refined as follows:

[0088] Construct the power output variables of renewable energy power generation equipment in the integrated energy station model and the scenario power output variables under different scenarios; determine the constraint comprehensive set based on the relaxed constraints, power output variables, scenario power output variables, preset relaxation conditions and the ramp limit information of renewable energy power generation equipment; determine the steady-state safety domain model of the integrated energy station model based on the regular stable equilibrium manifold representation of the preset quotient gradient system model and the constraint comprehensive set.

[0089] In this embodiment, the renewable energy power generation equipment may include conventional generators and photovoltaic power generation. The power output variable may be understood as the power generation output of the renewable energy power generation equipment. The scenario power output variable may be understood as the power generation output corresponding to different scenario changes. The preset relaxation condition may be understood as a relaxation variable used to constrain inequalities related to the renewable energy power generation equipment as equality constraints. The ramp limit information may be understood as the downward / upward ramp limit of the generator. The constraint comprehensive set may be understood as a set including all constraint conditions. The regular stable equilibrium manifold representation may be understood as a mathematical formula for describing the stability near the equilibrium point.

[0090] Specifically, due to the highly intermittent and volatile characteristics of renewable energy, the penetration level of renewable energy continues to increase, which brings challenges to the integrated energy station model. The uncertainty and variability of renewable energy lead to the rescheduling of the system. Considering the impact of renewable energy power generation equipment, the processor can construct the power output variables of renewable energy power generation equipment in the integrated energy station model and the scenario power output variables under different scenarios. The processor can determine the constraints related to renewable energy based on the power output variables, scenario power output variables and the ramp limit information of renewable energy power generation equipment, and convert the inequalities in the constraints related to renewable energy into equations through preset relaxation conditions to obtain the constraints related to renewable energy in the form of equations, and combine them with the relaxation constraints to obtain a comprehensive set of constraints. The processor can determine the steady-state safety domain model of the integrated energy station model based on the regular stable equilibrium manifold representation of the preset quotient gradient system model and the comprehensive set of constraints.

[0091] For example, the steady-state safety region SSR of the integrated energy station can be expressed as:

[0092]

[0093] Each point in the steady-state safety region represents a feasible control that satisfies the energy flow equations and safety operation constraints of the Integrated Energy System (IES) model. Due to nonlinear constraints, the IES steady-state safety region is usually non-convex and can be composed of multiple disjoint components. Based on the above relaxed constraints and the representation of the steady-state safety region, the complete representation of the steady-state safety region is explored through QGS in nonlinear dynamic systems:

[0094]

[0095] Where DH(x) T is the Jacobian matrix from the vector H(x).

[0096] The stable equilibrium manifold (SEM) of QGS is 0 if the set of equations is satisfied, that is, H(Σ s )=0, it is called regular SEM (RSEM), otherwise it is called degenerate SEM (DSEM).

[0097] According to the properties of nonlinear dynamic systems, a path-connected set, such as SSR, is a steady-state security component of the integrated energy system model if and only if it is a canonical stable equilibrium manifold corresponding to the relevant quotient gradient system (QGS) of the integrated energy system, that is:

[0098]

[0099] In addition, the steady-state safety region of the coupled energy system is equal to the union of the canonical stable equilibrium manifolds of the relevant quotient gradient systems corresponding to the coupled energy system, that is,

[0100]

[0101] In the formula, is the i-th regular stable equilibrium manifold of the associated quotient gradient system.

[0102] The relaxed constraint formulation is extended by including the uncertainty of renewable energy sources. Then, the formulation can provide a safe operating point based on the IES steady-state safety region.

[0103] The vector of renewable energy output can be represented as a set of scenario power output variables

[0104]

[0105] Where S is the total number of representative scenarios, s is the scenario index, and each scenario is a set of power outputs of PV power generation.

[0106] Formally, the constraint synthesis set can be concisely expressed as follows:

[0107]

[0108] In the formula, and is the power output variable of conventional generators and photovoltaic power generation in renewable energy power generation equipment under basic conditions, and are the scenario power output variables of conventional generators and PV generation in each representative scenario. The coupling between the base case and the representative scenario is introduced by the ramp limit of the generator. RD / RU are the ramp-down / upward limits of the generator. and is the slack variable that converts the inequality ramp constraints into equality constraints.

[0109] Similarly, QGS can be constructed based on the constraint comprehensive set, and we get:

[0110]

[0111] The steady-state safety region with respect to the uncertainty of renewable energy can be characterized by the canonical stable equilibrium manifold of the corresponding QGS.

[0112] S211. Obtain relevant network data and variable boundary information of the integrated energy station model, and determine the initial boundary point of the steady-state safety domain.

[0113] In this embodiment, the relevant network data can be understood as data related to building a comprehensive energy station model. The variable boundary information can be understood as the maximum and minimum values ​​related to the control variables and state variables. The initial boundary point can be understood as the boundary point determined for the first time.

[0114] Specifically, the processor can obtain relevant network data and variable boundary information of the integrated energy station model input by relevant personnel, and obtain the initial boundary point by constructing the QGS corresponding to the steady-state safety domain and integrating it.

[0115] S212. According to each initial boundary point and the steady-state safety domain model, a boundary search direction and feasible solutions under each boundary search direction are determined by a unified numerical trajectory algorithm.

[0116] In this embodiment, the boundary search direction can be understood as a direction for searching and determining a feasible solution. A feasible solution can be understood as a solution that satisfies the constraint conditions.

[0117] Specifically, the processor can connect the initial boundary points, determine the boundary search direction, construct QGS through the steady-state safety domain model, and determine the feasible solution of QGS in each boundary search direction through the numerical trajectory unified algorithm.

[0118] Further, based on the above embodiment, the steps of determining the boundary search direction and the feasible solutions under each boundary search direction by the unified numerical trajectory algorithm according to each initial boundary point and the steady-state safety domain model can be refined as follows:

[0119] Connect the initial boundary points, determine the boundary search direction, and the initial search point in the boundary search direction; construct the target quotient gradient system model corresponding to the steady-state safety domain model; perform trajectory integration on the target quotient gradient system model based on the initial search point and the unified algorithm of the numerical trajectory to determine the energy function value; determine the trajectory endpoint value based on the energy function values ​​on the search trajectory and the preset threshold; solve based on the preset numerical algorithm and the trajectory endpoint value to obtain a feasible solution in the boundary search direction.

[0120] In this embodiment, the target quotient gradient system model can be understood as a quotient gradient system model constructed by a steady-state safety domain model. The initial search point can be understood as the first point used for searching. The energy function value can be understood as a function expression used to describe the state variables of the system energy domain. The preset threshold can be understood as a threshold set for determining the end point of the trajectory. The trajectory end point value can be understood as the value corresponding to the end point of the boundary search trajectory. The preset numerical algorithm can be understood as an algorithm for calculating a feasible solution, such as Newton's method or interior point method.

[0121] Specifically, the processor can connect each initial boundary point, determine the boundary search direction, and determine the initial point search point (the midpoint of the connecting line between the boundary points) in the search direction. The processor can construct the corresponding target quotient gradient system model through the steady-state safety domain model, and give the initial search point x0, where H'(x0)≠0. Starting from the initial search point x0, the corresponding QGS trajectory is integrated to calculate the energy function value H'(x(t)); the energy function value H'(x(t)) is monitored along the trajectory, and the integration is stopped when H'(x(t)) is less than a preset threshold to obtain the end point of the trajectory. The processor can apply a numerical method (such as Newton's method or interior point method), use the trajectory end point as an initial guess to solve the equation H'(x), and output the solution as a feasible solution.

[0122] S213. According to the number of feasible solutions and the convergence region of the regular stable equilibrium manifold in the target quotient gradient system model, the boundary points in different boundary search directions are determined as the boundaries of the steady-state safety domain of the integrated energy station model.

[0123] In this embodiment, the convergence region can be understood as a region where all trajectories tend to the manifold over time.

[0124] Specifically, the processor can compare the number of feasible solutions with a preset number threshold. If the number of feasible solutions reaches the preset number threshold, it determines whether all non-overlapping safety components are obtained. If so, it stops; if not, it tries another initial search point located in the convergence area of ​​an adjacent regular stable manifold and returns to execute step S211. If the number of feasible solutions does not reach the preset number threshold, it tries another boundary search direction and returns to continue executing the step of giving an initial search point in another boundary search direction (i.e., step S212).

[0125] The technical solution of the embodiment of the present invention, by establishing the equivalence between the steady-state safety domain and the regular stable equilibrium manifold of the non-hyperbolic dynamic system constructed by the multi-integrated energy station model, takes into account the influence of renewable energy uncertainty on the steady-state safety domain, and can fully characterize the accurate model of the steady-state safety region of multiple energy stations. The proposed method does not require matrix inversion, and therefore avoids the divergence problem caused by the Jacobian matrix along the trajectory along possible pathological conditions. In addition, the global convergence of the TJU method ensures that any initial point must converge to a feasible solution (if it exists). Combined with QGS and local methods, such as the interior point method, the TJU method is powerful and effective, fully describing the IES steady-state safety domain and having a certain robustness. This method is to calculate the steady-state solution of QGS, not the trajectory itself. In the process of exploring new boundary points, the calculation of a large number of meaningless points in the safety region is avoided. By calculating the boundary points of different boundary search directions in parallel, the computational efficiency can be improved.

[0126] For example, in order to facilitate understanding of the search process, a specific example can be used to illustrate: Figure 4 The following is a schematic diagram of boundary calculation of a method for determining the boundary of a steady-state security domain provided in the second embodiment of the present invention. Figure 4 As shown, the boundary points include X1, X2, X3, X4 and X5, V represents the boundary search direction, and the symbols in the above example are used to construct the QGS corresponding to the steady-state security domain model, where QGS represents the direction V. T Parameters of the projection onto :

[0127]

[0128] By continuously changing the moving step length λ to λ = λ + Δλ, we can T A series of feasible solutions are obtained in the direction until the boundary of the steady-state safety region in that direction is reached, such as Figure 4 As shown. Calculate some feasible points on the boundary of the steady-state safety region. V T is the perpendicular direction of the line connecting the two boundary points. Then, we can follow the direction V 12 Arrived at the new boundary point X5. A new search direction can be generated by connecting X5 with X1 or X2. By continuously changing the direction, the complete boundary of the steady-state safety area can be calculated.

[0129] Embodiment 3

[0130] Figure 5 This is a schematic diagram of the structure of a device for determining the boundary of a stable security domain provided in Embodiment 3 of the present invention. Figure 5 As shown, the device includes: an information acquisition module 51, an information determination module 52, a model construction module 53 and a boundary determination module 54.

[0131] The information acquisition module 51 is used to establish a comprehensive energy station model under the electric heating system, and obtain relevant parameter information of the comprehensive energy station model and load information of each energy station;

[0132] An information determination module 52, configured to determine operation constraint information satisfying the operation of the electric heating system according to the load information and the relevant parameter information;

[0133] A model building module 53, used to build a steady-state safety domain model of the integrated energy station model according to the operation constraint information and a preset quotient gradient system model;

[0134] The boundary determination module 54 is used to solve the steady-state safety domain model based on the numerical trajectory unified algorithm to obtain the boundary of the steady-state safety domain of the integrated energy station model.

[0135] The technical solution of the embodiment of the present invention is to establish an integrated energy station model under the electric heating system to obtain relevant parameter information of the integrated energy station model and load information of each energy station; determine the operation constraint information that satisfies the operation of the electric heating system based on the load information and relevant parameter information; construct a steady-state safety domain model of the integrated energy station model based on the operation constraint information and the preset quotient gradient system model; solve the steady-state safety domain model based on the numerical trajectory unified algorithm to obtain the boundary of the steady-state safety domain of the integrated energy station model. The steady-state safety domain model is constructed by the preset quotient gradient system model, and the accurate model of the steady-state safety area of ​​multiple energy stations is completely characterized. The boundary is determined by solving the boundary through the numerical trajectory unified algorithm. The calculation of a large number of meaningless points in the safety area is avoided, the calculation efficiency is improved, and the accuracy of the boundary determination is guaranteed.

[0136] Furthermore, the information determination module is specifically used for:

[0137] Determine the number of access nodes, transmission lines, energy stations, heating units and power generation equipment in the integrated energy station model according to the network topology in the relevant parameter information;

[0138] Determine a node constraint condition related to the access node according to the quantity information of the power generation equipment, the node attribute information of each access node in the relevant parameter information, and the equivalent parameter information of each transmission line;

[0139] Determine the output power constraint of the heating unit according to the quantity information of the power generation equipment and the heat load demand of each of the energy stations in the load information;

[0140] The output power constraint condition and each of the node constraint conditions are used to determine the operation constraint information that satisfies the operation of the electric heating system.

[0141] Furthermore, the model building module includes:

[0142] A first determining unit is used to use the voltage angle and voltage amplitude of the active power-reactive power PQ node in the integrated energy station model as state variable information based on the operation constraint information;

[0143] A second determination unit is used to use the actual power generation and the actual voltage value of the photovoltaic node in the integrated energy station model as control variable information based on the operation constraint information;

[0144] A third determining unit, configured to determine power constraint information of the power subsystem and heating constraint information of the heating subsystem in the electric heating system according to the state variable information and the control variable information;

[0145] a fourth determining unit, configured to convert the power inequality constraint in the power constraint information and the heat supply inequality constraint in the heat supply constraint information respectively through preset slack variable information, so as to obtain a slack constraint condition of the integrated energy station model;

[0146] The fifth determination unit is used to determine the steady-state safety domain model of the integrated energy station model according to the relaxed constraint conditions and the preset quotient gradient system model.

[0147] The fifth determining unit is specifically used for:

[0148] Constructing power output variables of renewable energy power generation equipment in the integrated energy station model and scenario power output variables under different scenarios;

[0149] Determining a comprehensive constraint set according to the relaxed constraint condition, the power output variable, the scenario power output variable, the preset relaxation condition and the ramp-up limit information of the renewable energy power generation equipment;

[0150] According to the canonical stable equilibrium manifold representation of the preset quotient gradient system model and the comprehensive set of constraints, the steady-state safety domain model of the integrated energy station model is determined.

[0151] Furthermore, the boundary determination module includes:

[0152] a sixth determination unit, configured to obtain relevant network data and variable boundary information of the integrated energy station model and determine an initial boundary point of a steady-state safety domain;

[0153] a seventh determination unit, configured to determine a boundary search direction and a feasible solution under each boundary search direction by a unified numerical trajectory algorithm according to each of the initial boundary points and the steady-state safety domain model;

[0154] The eighth determination unit is used to determine the boundary points of different boundary search directions as the boundaries of the steady-state safety domain of the integrated energy station model according to the number of feasible solutions and the convergence area of ​​the regular stable equilibrium manifold in the target quotient gradient system model.

[0155] Wherein, the seventh determining unit is specifically used for:

[0156] Connecting the initial boundary points, determining a boundary search direction, and an initial search point in the boundary search direction;

[0157] Constructing a target quotient gradient system model corresponding to the steady-state safety domain model;

[0158] Based on the initial search point and the numerical trajectory unified algorithm, the target quotient gradient system model is trajectory-integrated to determine the energy function value;

[0159] Determine the trajectory endpoint value according to each of the energy function values ​​on the search trajectory and a preset threshold;

[0160] A solution is obtained based on a preset numerical algorithm and the trajectory endpoint value pair to obtain a feasible solution in the boundary search direction.

[0161] The device for determining the boundary of a steady-state security domain provided in an embodiment of the present invention can execute the method for determining the boundary of a steady-state security domain provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0162] Embodiment 4

[0163] Figure 6 A schematic diagram of the structure of an electronic device 60 that can be used to implement an embodiment of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or required herein.

[0164] like Figure 6As shown, the electronic device 60 includes at least one processor 61, and a memory connected to the at least one processor 61, such as a read-only memory (ROM) 62, a random access memory (RAM) 63, etc., wherein the memory stores a computer program that can be executed by at least one processor, and the processor 61 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 62 or the computer program loaded from the storage unit 68 to the random access memory (RAM) 63. In the RAM 63, various programs and data required for the operation of the electronic device 60 can also be stored. The processor 61, the ROM 62, and the RAM 63 are connected to each other via a bus 64. An input / output (I / O) interface 65 is also connected to the bus 64.

[0165] A number of components in the electronic device 60 are connected to the I / O interface 65, including: an input unit 66, such as a keyboard, a mouse, etc.; an output unit 67, such as various types of displays, speakers, etc.; a storage unit 68, such as a disk, an optical disk, etc.; and a communication unit 69, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 69 allows the electronic device 60 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0166] The processor 61 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the processor 61 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any appropriate processor, controller, microcontroller, etc. The processor 61 executes the various methods and processes described above, such as a boundary determination method for a steady-state security domain.

[0167] In some embodiments, the boundary determination method of the steady-state security domain may be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 68. In some embodiments, part or all of the computer program may be loaded and / or installed on the electronic device 60 via the ROM 62 and / or the communication unit 69. When the computer program is loaded into the RAM 63 and executed by the processor 61, one or more steps of the boundary determination method of the steady-state security domain described above may be performed. Alternatively, in other embodiments, the processor 61 may be configured to execute the boundary determination method of the steady-state security domain in any other appropriate manner (e.g., by means of firmware).

[0168] Various implementations of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on chips (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0169] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, so that when the computer program is executed by the processor, the functions / operations specified in the flow chart and / or block diagram are implemented. The computer program may be executed entirely on the machine, partially on the machine, partially on the machine and partially on a remote machine as a stand-alone software package, or entirely on a remote machine or server.

[0170] In the context of the present invention, a computer-readable storage medium may be a tangible medium that may contain or store a computer program for use by or in conjunction with an instruction execution system, device, or equipment. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices, or equipment, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. A more specific example of a machine-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0171] To provide interaction with a user, the systems and techniques described herein may be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices may also be used to provide interaction with the user; for example, the feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form (including acoustic input, voice input, or tactile input).

[0172] The systems and techniques described herein may be implemented in a computing system that includes backend components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes frontend components (e.g., a user computer with a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such backend components, middleware components, or frontend components. The components of the system may be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0173] A computing system may include a client and a server. The client and the server are generally remote from each other and usually interact through a communication network. The client and server relationship is generated by computer programs running on the corresponding computers and having a client-server relationship with each other. The server may be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system to solve the defects of difficult management and weak business scalability in traditional physical hosts and VPS services.

[0174] In one embodiment, the embodiment of the present invention further includes a computer program product, the computer program product includes a computer program, and when the computer program is executed by a processor, it implements the boundary determination method of the steady-state security domain of any embodiment of the present invention.

[0175] In the process of implementation, the computer program product can be written in one or more programming languages ​​or a combination thereof to perform the computer program code of the present invention, including object-oriented programming languages, such as Java, Smalltalk, C++, and conventional procedural programming languages, such as "C" language or similar programming languages. The program code can be executed entirely on the user's computer, partially on the user's computer, as a separate software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In the case of a remote computer, the remote computer can be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or can be connected to an external computer (for example, using an Internet service provider to connect through the Internet).

[0176] It should be understood that the various forms of processes shown above can be used to reorder, add or delete steps. For example, the steps described in the present invention can be executed in parallel, sequentially or in different orders, as long as the desired results of the technical solution of the present invention can be achieved, and this document does not limit this.

[0177] The above specific implementations do not constitute a limitation on the protection scope of the present invention. It should be understood by those skilled in the art that various modifications, combinations, sub-combinations and substitutions can be made according to design requirements and other factors. Any modification, equivalent substitution and improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for determining the boundary of a steady-state security domain, characterized in that: include: Establishing an integrated energy station model under the electric heating system, and obtaining relevant parameter information of the integrated energy station model and load information of each energy station; Determining, based on the load information and the relevant parameter information, operation constraint information that satisfies the operation of the electric heating system; Constructing a steady-state safety domain model of the integrated energy station model according to the operation constraint information and the preset quotient gradient system model; The steady-state safety domain model is solved based on the unified numerical trajectory algorithm to obtain the boundary of the steady-state safety domain of the integrated energy station model.

2. The method according to claim 1, characterized in that The determining, based on the load information and the related parameter information, operation constraint information satisfying the operation of the electric heating system includes: Determine the number of access nodes, transmission lines, energy stations, heating units and power generation equipment in the integrated energy station model according to the network topology in the relevant parameter information; Determine a node constraint condition related to the access node according to the quantity information of the power generation equipment, the node attribute information of each access node in the relevant parameter information, and the equivalent parameter information of each transmission line; Determine the output power constraint of the heating unit according to the quantity information of the power generation equipment and the heat load demand of each of the energy stations in the load information; The output power constraint condition and each of the node constraint conditions are used to determine the operation constraint information that satisfies the operation of the electric heating system.

3. The method according to claim 1, characterized in that The step of constructing a steady-state safety domain model of the integrated energy station model according to the operation constraint information and the preset quotient gradient system model includes: Based on the operation constraint information, the voltage angle and voltage amplitude of the active power-reactive power PQ node in the integrated energy station model are used as state variable information; Based on the operation constraint information, the actual power generation and the actual voltage value of the photovoltaic node in the integrated energy station model are used as control variable information; Determining power constraint information of the power subsystem and heating constraint information of the heating subsystem in the electric heating system according to the state variable information and the control variable information; The power inequality constraints in the power constraint information and the heat supply inequality constraints in the heat supply constraint information are respectively converted by presetting slack variable information to obtain slack constraint conditions of the integrated energy station model; According to the relaxed constraint conditions and the preset quotient gradient system model, the steady-state safety domain model of the integrated energy station model is determined.

4. The method according to claim 3, characterized in that Determining the steady-state safety domain model of the integrated energy station model according to the relaxed constraint conditions and the preset quotient gradient system model includes: Constructing power output variables of renewable energy power generation equipment in the integrated energy station model and scenario power output variables under different scenarios; Determining a comprehensive constraint set according to the relaxed constraint condition, the power output variable, the scenario power output variable, the preset relaxation condition and the ramp-up limit information of the renewable energy power generation equipment; According to the canonical stable equilibrium manifold representation of the preset quotient gradient system model and the comprehensive set of constraints, the steady-state safety domain model of the integrated energy station model is determined.

5. The method according to claim 1, characterized in that The method of solving the steady-state safety domain model based on the unified numerical trajectory algorithm to obtain the boundary of the steady-state safety domain of the integrated energy station model includes: Obtain relevant network data and variable boundary information of the integrated energy station model, and determine the initial boundary point of the steady-state safety domain; According to each of the initial boundary points and the steady-state safety domain model, a boundary search direction and a feasible solution under each of the boundary search directions are determined by a unified numerical trajectory algorithm; According to the number of feasible solutions and the convergence region of the regular stable equilibrium manifold in the target quotient gradient system model, boundary points in different boundary search directions are determined as the boundaries of the steady-state safety domain of the integrated energy station model.

6. The method according to claim 5, characterized in that The step of determining the boundary search direction and the feasible solutions under each boundary search direction by using a unified numerical trajectory algorithm according to each of the initial boundary points and the steady-state safety domain model includes: Connecting the initial boundary points, determining a boundary search direction, and an initial search point in the boundary search direction; Constructing a target quotient gradient system model corresponding to the steady-state safety domain model; Based on the initial search point and the numerical trajectory unified algorithm, the target quotient gradient system model is trajectory-integrated to determine the energy function value; Determine the trajectory endpoint value according to each of the energy function values ​​on the search trajectory and a preset threshold; A solution is obtained based on a preset numerical algorithm and the trajectory endpoint value pair to obtain a feasible solution in the boundary search direction.

7. A device for determining the boundary of a steady-state security zone, characterized in that: include: An information acquisition module is used to establish a comprehensive energy station model under the electric heating system, and obtain relevant parameter information of the comprehensive energy station model and load information of each energy station; An information determination module, used to determine operation constraint information satisfying the operation of the electric heating system according to the load information and the relevant parameter information; A model building module, used to build a steady-state safety domain model of the integrated energy station model according to the operation constraint information and a preset quotient gradient system model; The boundary determination module is used to solve the steady-state safety domain model based on the numerical trajectory unified algorithm to obtain the boundary of the steady-state safety domain of the integrated energy station model.

8. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for determining the boundary of a steady-state security domain according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the method for determining the boundary of a steady-state security domain according to any one of claims 1 to 6 when executed.

10. A computer program product, characterized in that The computer program product comprises a computer program, which, when executed by a processor, implements the method for determining the boundary of a steady-state security domain according to any one of claims 1 to 6.

Images