Integrated energy system distributed scheduling method and system considering typhoon uncertainty

Abstract

The invention provides an integrated energy system distributed scheduling method and system considering typhoon uncertainty, and relates to the technical field of power-gas integrated energy cooperative operation and toughness enhancement, and the method comprises the steps: constructing an offshore wind power operation model, and employing a Monte Carlo simulation method to carry out the sampling and generation of a typhoon scene and an offshore wind power output scene; based on a typhoon scene and an offshore wind power output scene, establishing a scene-based transmission-distribution-gas integrated energy system three-layer collaborative scheduling model; the method comprises the following steps of: performing solvability conversion on a scene-based transmission-distribution-gas integrated energy system three-layer collaborative scheduling model based on mixed sample average approximation into a determined continuous linear model; and solving the determined continuous linear model by using an efficient distributed collaborative algorithm of a multi-parameter planning theory and a critical domain accelerated search algorithm based on global optimal solution estimation to obtain a global optimal solution. The solving efficiency can meet the timeliness requirement of emergency response in a typhoon scene.

Description

Technical Field

[0001] The present disclosure relates to the technical field of coordinated operation and resilience enhancement of electricity-gas integrated energy, and specifically to a distributed scheduling method and system for an integrated energy system taking into account typhoon uncertainty. Background Art

[0002] The statements in this section merely provide background information related to the present disclosure and do not necessarily constitute prior art.

[0003] Offshore Wind Power (OWP) has entered a rapid development stage as an important renewable energy source. However, the impact of typhoon extreme weather on OWP is significantly stronger than onshore wind power. Although most typhoons do not directly land, they will still force wind turbines to shut down at excessively high wind speeds, which will lead to wind power ramp events (WPRE). This short-term, large-scale power shortage poses a huge challenge to the power balance of the power system.

[0004] At present, there are many solutions to the problem of safe and economic dispatch of coastal power systems under typhoon weather, but there are still the following shortcomings: 1) Existing solutions usually rely on pre-set typhoon scenarios for modeling, and fail to fully consider the uncertainty of typhoon path, movement speed, central air pressure and other parameters; 2) The existing scheme fails to take into account the risk of damage to wind turbines during typhoons, which will actually increase the demand for system backup capacity; 3) The traditional dispatching mode adopted by the existing scheme is limited to the transmission network level, and it is not sufficient to rely solely on the limited regulation capacity of thermal power units to cope with the uncertainty of WPRE and wind power output.

[0005] 4) With the popularization of flexible resources such as distributed power sources, energy storage equipment, and gas turbines on the distribution network side, as well as the addition of electric-gas coupling equipment, the distribution network and natural gas network are more involved in the dispatching of the transmission network side, providing a promising solution for dealing with the impact of typhoons. However, although the three-layer collaborative dispatching mode of this transmission-distribution-gas integrated energy system is more flexible, it is also more complex. Because the transmission network, distribution network, and natural gas network belong to different operating entities, the use of centralized optimization methods cannot guarantee privacy security and the computational burden is too heavy; in addition, due to the nested interactions between different levels in the three-layer collaborative dispatching mode, the traditional distributed algorithm solution efficiency is low, and it is difficult to meet the timeliness requirements of emergency response in typhoon scenarios. Summary of the invention

[0006] In order to solve the above problems, the present disclosure proposes a distributed scheduling method and system for an integrated energy system taking into account typhoon uncertainty, establishes an OWP operation model taking into account typhoon parameter uncertainty and wind turbine damage risk, so as to fully characterize the active output characteristics of the OWP in a typhoon scenario; proposes a three-layer collaborative scheduling model for a transmission-distribution-gas integrated energy system, so as to make full use of the flexibility resources in different networks to cope with the uncertainty of WPRE and wind power output, and introduces an efficient distributed collaborative algorithm based on multi-parameter planning theory and a critical domain acceleration search algorithm based on global optimal solution estimation to further improve its convergence speed, so as to ensure that the solution efficiency of the above three-layer collaborative scheduling model can meet the timeliness requirements of emergency response in typhoon scenarios.

[0007] According to some embodiments, the present disclosure adopts the following technical solutions: The distributed dispatching method of integrated energy system considering typhoon uncertainty includes: Considering the impact of typhoon uncertainty on offshore wind power, an offshore wind power operation model is constructed. Based on the offshore wind power operation model, the Monte Carlo simulation method is used to sample and generate typhoon scenarios and offshore wind power output scenarios; Based on typhoon scenarios and offshore wind power output scenarios, the coordinated operation of the transmission network, distribution network and natural gas network is considered, and a scenario-based three-layer coordinated dispatching model of the transmission-distribution-gas integrated energy system is established; Based on the mixed sample average approximation, the three-layer collaborative scheduling model of the scenario-based transmission-distribution-gas integrated energy system is transformed into a deterministic continuous linear model. The deterministic continuous linear model is solved using an efficient distributed collaborative algorithm based on multi-parameter planning theory and a critical domain accelerated search algorithm based on global optimal solution estimation to obtain the global optimal solution for the distributed scheduling of the integrated energy system under typhoon scenarios.

[0008] According to some embodiments, the present disclosure adopts the following technical solutions: The distributed dispatching system of the integrated energy system taking into account the uncertainty of typhoons includes: The scenario generation module is used to consider the impact of typhoon uncertainty on offshore wind power, build an offshore wind power operation model, and use the Monte Carlo simulation method to sample and generate typhoon scenarios and offshore wind power output scenarios based on the offshore wind power operation model; The model building module is used to establish a scenario-based three-layer coordinated dispatching model of the transmission-distribution-gas integrated energy system based on typhoon scenarios and offshore wind power output scenarios, taking into account the coordinated operation of the transmission network, distribution network and natural gas network. The scheduling solution module is used to transform the scenario-based transmission-distribution-gas integrated energy system three-layer collaborative scheduling model into a deterministic continuous linear model based on the mixed sample average approximation, and use the efficient distributed collaborative algorithm of multi-parameter planning theory and the critical domain acceleration search algorithm based on the global optimal solution estimation to solve the deterministic continuous linear model, and obtain the global optimal solution for the distributed scheduling of the integrated energy system under the typhoon scenario.

[0009] According to some embodiments, the present disclosure adopts the following technical solutions: A computer program product includes a computer program, which, when executed by a processor, implements the distributed scheduling method for an integrated energy system taking into account typhoon uncertainty.

[0010] According to some embodiments, the present disclosure adopts the following technical solutions: A non-transitory computer-readable storage medium is used to store computer instructions. When the computer instructions are executed by a processor, the distributed scheduling method of the integrated energy system taking into account typhoon uncertainty is implemented.

[0011] According to some embodiments, the present disclosure adopts the following technical solutions: An electronic device comprises: a processor, a memory and a computer program; wherein the processor is connected to the memory, the computer program is stored in the memory, and when the electronic device is running, the processor executes the computer program stored in the memory so that the electronic device executes the distributed scheduling method of the integrated energy system taking into account the uncertainty of typhoons.

[0012] Compared with the prior art, the present invention has the following beneficial effects: The disclosed distributed scheduling method for an integrated energy system taking into account typhoon uncertainty establishes an OWP operation model taking into account typhoon parameter uncertainty and wind turbine damage risk, so as to fully characterize the active output characteristics of the OWP in a typhoon scenario; proposes a scenario-based three-layer collaborative scheduling model for a transmission-distribution-gas integrated energy system, so as to make full use of the flexibility resources in different networks to cope with the uncertainty of WPRE and wind power output; proposes an efficient distributed collaborative algorithm based on multi-parameter planning theory, and introduces a critical domain acceleration search algorithm based on global optimal solution estimation to further improve its convergence speed, so as to ensure that the solution efficiency of the above three-layer collaborative scheduling model can meet the timeliness requirements of emergency response in typhoon scenarios.

[0013] The disclosed distributed scheduling method for the integrated energy system taking into account typhoon uncertainty is a scenario-based three-layer coordinated scheduling model for the transmission-distribution-gas integrated energy system, including a transmission network-side coordinated scheduling model, a distribution network-side coordinated scheduling model, and a natural gas network-side coordinated scheduling model. It can realize distributed coordinated scheduling of the transmission network, distribution network, and natural gas network, and make full use of the flexibility resources in different networks to ensure the safe operation of the system in typhoon scenarios.

[0014] The distributed scheduling method for the integrated energy system taking into account typhoon uncertainty disclosed in the present invention is based on an efficient distributed collaborative algorithm based on multi-parameter planning theory, and introduces a critical domain acceleration search algorithm based on global optimal solution estimation to solve the scenario-based transmission-distribution-gas integrated energy system three-layer collaborative scheduling model. Compared with traditional distributed algorithms, such as the alternating direction multiplier method, the Benders decomposition method, etc., it has a significant advantage in solution speed, and the calculation time can be shortened by more than 2-3 times. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] The accompanying drawings constituting a part of the present disclosure are used to provide a further understanding of the present disclosure. The illustrative embodiments of the present disclosure and their descriptions are used to explain the present disclosure and do not constitute an improper limitation on the present disclosure.

[0016] Figure 1 This is a flowchart for implementing the application of the distributed scheduling method for an integrated energy system taking into account typhoon uncertainty according to an embodiment of the present disclosure. DETAILED DESCRIPTION

[0017] The present disclosure is further described below in conjunction with the accompanying drawings and embodiments.

[0018] It should be noted that the following detailed descriptions are all illustrative and are intended to provide further explanation of the present disclosure. Unless otherwise specified, all technical and scientific terms used herein have the same meanings as those commonly understood by those skilled in the art to which the present disclosure belongs.

[0019] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present disclosure. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprising" and / or "including" are used in this specification, it indicates the presence of features, steps, operations, devices, components and / or combinations thereof.

[0020] Example 1 In one embodiment of the present disclosure, a distributed scheduling method for an integrated energy system taking into account typhoon uncertainty is provided, which can realize distributed coordinated scheduling of a transmission network, a distribution network, and a natural gas network, and at the same time fully utilize the flexibility resources in different networks to ensure the safe operation of the system in a typhoon scenario, and the steps include: Step 1: Considering the impact of typhoon uncertainty on offshore wind power, an offshore wind power operation model is constructed. Based on the offshore wind power operation model, the Monte Carlo simulation method is used to sample and generate typhoon scenarios and offshore wind power output scenarios; Step 2: Based on typhoon scenarios and offshore wind power output scenarios, consider the coordinated operation of the transmission network, distribution network, and natural gas network, and establish a scenario-based three-layer coordinated dispatching model for the transmission-distribution-gas integrated energy system; Step 3: Based on the mixed sample average approximation, the scenario-based transmission-distribution-gas integrated energy system three-layer coordinated scheduling model is transformed into a deterministic continuous linear model; Step 4: Use the efficient distributed collaborative algorithm of multi-parameter planning theory and the critical domain acceleration search algorithm based on global optimal solution estimation to solve the determined continuous linear model and obtain the global optimal solution for distributed scheduling of integrated energy systems under typhoon scenarios.

[0021] As an embodiment, the implementation process of the distributed scheduling method of the integrated energy system taking into account typhoon uncertainty disclosed in the present invention is as follows: Figure 1 The specific implementation process is as follows: Step 1: Construct an offshore wind power operation model that considers the uncertainty of typhoon parameters and the risk of wind turbine damage, use the Monte Carlo simulation method to sample and generate typhoon scenarios and offshore wind power output scenarios, fully characterize the active output characteristics of OWP in typhoon scenarios, and lay the foundation for the subsequent construction of a three-layer coordinated dispatching model of the transmission-distribution-gas integrated energy system, including: considering the impact of typhoon uncertainty on offshore wind power, using the Holland vortex model to model the typhoon wind farm, and modeling the uncertainty of typhoon parameters. When the probability distribution satisfied by the uncertain parameters of the typhoon is determined, a typhoon scenario is generated using a Monte Carlo simulation-based method; for any offshore wind farm, in a specific typhoon scenario, a wind speed-power active output model is constructed. According to the relative distance between the typhoon and the offshore wind farm, the entire process of the offshore wind farm encountering a typhoon event is divided into three stages, namely, the typhoon approach stage, the typhoon transit stage, and the typhoon departure stage. The wind power output vectors of the three stages are integrated to obtain the complete output of the offshore wind farm in the typhoon scenario. The specific implementation is as follows: Step 1-1: Typhoon wind field modeling, as follows: Step 1-1-1: Use the Holland vortex model to model the typhoon wind field. Its mathematical form can be expressed as follows: (1) (2) (3) in, is the set of all scheduling periods; For a specific scheduling period; The straight-line distance from the typhoon center Gradient wind speed value at the measuring point; is the Coriolis force parameter; is the typhoon center pressure value, Characterizes the pressure difference between the typhoon center pressure and the ambient atmospheric pressure; is the air density parameter; is the latitude of the typhoon eye; is the correction factor, and its typical value range is between 0 and 0.4. is the Holland radial pressure curve parameter; It is the maximum wind speed radius of the typhoon.

[0022] Step 1-2: Modeling of typhoon parameter uncertainty, as follows: Step 1-2-1: The uncertainty of the typhoon path comes from the direction angle With translation speed The prediction errors are denoted as and . Prediction Error and It can be considered to satisfy the following normal distribution: (4) in, and Respectively represent the mean parameter and variance parameter of the normal distribution, both of which can be estimated through historical data; * is used to refer to and , is a commonly used mathematical representation.

[0023] Step 1-2-2: Typhoon center pressure The uncertainty comes from its own prediction error. Here, the probability distribution satisfied by the typhoon center pressure is directly modeled, which can be modeled as the following lognormal distribution: (5) in, Represents the predicted value of the typhoon's central pressure, and They respectively characterize the mean parameter and variance parameter of the lognormal distribution, both of which can be estimated through historical data.

[0024] Step 1-3: Typhoon scenario generation. When the probability distribution satisfied by the uncertain parameters of the typhoon is determined, the typhoon scenario can be generated by Monte Carlo simulation; the details are as follows: Step 1-3-1: Get the typhoon eye position for the current period , the predicted position for the next period And the predicted typhoon center pressure value for the next period .

[0025] Step 1-3-2: Based on the current typhoon eye position and the predicted position for the next period , calculate the predicted azimuth change value and predicted translational velocity .

[0026] Step 1-3-3: Based on the Monte Carlo simulation method, the typhoon path prediction error is calculated according to formula (4). and Sampling is performed, and the possible typhoon center pressure value in the next period is calculated based on formula (5). Conduct sampling.

[0027] Step 1-3-4: Calculate the azimuth change value of the simulated sampling and the translation speed of the simulated sampling , and then calculate the possible typhoon eye position in the next period .in, and The azimuth and translation speed of the current period Step 1-3-5: Return to step 1-3-1 and repeat steps 1-3-1 to 1-3-4 until the last period, and finally obtain a complete simulated typhoon scenario.

[0028] Step 1-3-6: Assume the number of scenes to be generated is , repeat steps 1-3-1 to 1-3-5 until you get simulated typhoon scenarios. The index of each typhoon scenario is defined as , the set of all scenes is defined as .

[0029] Step 1-4: Generate offshore wind farm output scenario, define the set of offshore wind farms as , as follows: Step 1-4-1: For any offshore wind farm , in a specific typhoon scenario Under this condition, its active output can be described by the following wind speed-power model: (6) in, , and Respectively represent the rated wind speed, cut-in wind speed and cut-out wind speed; is the rated active power of the wind turbine, For offshore wind farms The wind speed at the place.

[0030] Step 1-4-2: According to the relative distance between the typhoon and the offshore wind farm, the whole process of the offshore wind farm encountering a typhoon event can be divided into three stages, namely, the approach stage , transit stage and the departure stage .in, , , as well as They are the time when the typhoon encounters, the time when offshore wind power is shut down, the time when offshore wind power resumes power generation, and the time when the typhoon impact disappears.

[0031] Step 1-4-3: In the approach phase , Formula (6) can be used to directly calculate the offshore wind farm The active output of .

[0032] Step 1-4-4: In transit , offshore wind farm The wind speed is always greater than the cut-out wind speed of the wind turbines inside it, that is, At this time, according to formula (6), the offshore wind farm The active output is .

[0033] Step 1-4-5: In transit , if the wind turbine is within the critical damage radius of the typhoon The distance between the wind turbine and the typhoon center , then the wind turbine may be damaged. Its damage probability can be described by the following probability distribution function: (7) in, For a single The damage probability of wind turbines at the moment; represents the standard normal cumulative distribution function; Characterizes the standard deviation of the natural logarithm of the wind speed design corresponding to when the wind turbine reaches the damage threshold; is the median value of the design wind speed; and They are respectively the moments of entering and leaving the critical destructive radius of the typhoon.

[0034] Step 1-4-6: During the departure phase As the typhoon center gradually moves away, offshore wind farms The wind speed at the location gradually weakens to below the cut-out wind speed. At this time, the wind turbines that were not damaged during the transit phase will resume power generation, and their active output can be calculated by the following formula: (8) in, It can be calculated by formula (6).

[0035] Step 1-4-7: Vectorize the wind power output of the three phases , , Put together, you get an offshore wind farm In the typhoon scene The complete force vector for (9) Step 1-4-8: For offshore wind farms ,right All scenarios in the above example are executed by executing steps 1-4-1 to 1-4-7 to obtain the offshore wind farm. The number of complete output vectors in all scenarios should be equal to the total number of scenarios .

[0036] Step 1-4-9: Define the offshore wind farm The complete output vector in the prediction scenario is , then the offshore wind farm In any scene The complete output vector under the prediction scenario can be expressed as With the scene The deviation value of the complete output vector compared to the predicted scenario The summation form is: (10) Step 1-4-10: For all the For all offshore wind farms, execute steps 1-4-1 to 1-4-9 to obtain the complete output vectors of all offshore wind farms in all typhoon scenarios, as shown below: (11) Step 2: Based on the stochastic optimization theory and the typhoon scenario and offshore output scenario generated in step 1, a three-layer coordinated dispatching model of the transmission-distribution-gas integrated energy system is constructed, including: establishing a scenario-based three-layer coordinated dispatching model of the transmission-distribution-gas integrated energy system, wherein the scenario-based three-layer coordinated dispatching model of the transmission-distribution-gas integrated energy system includes a transmission network-side coordinated dispatching model, a distribution network-side coordinated dispatching model, and a natural gas network-side coordinated dispatching model; wherein the objective function of the transmission network-side coordinated dispatching model is to minimize the transmission network-side operating cost, the objective function of the distribution network-side coordinated dispatching model is to minimize the distribution network-side operating cost, and the objective function of the natural gas network-side coordinated dispatching model is to minimize the natural gas network-side operating cost, and the constraints of each coordinated dispatching model are set separately. The details are as follows: Step 2-1: Establish a transmission network-side coordinated dispatching model, which includes objective functions and constraints, as follows: Step 2-1-1: The objective function of the transmission network side coordinated dispatching model is to minimize the transmission network side operating cost, and its mathematical form is as follows: (12) in, is the set of all scheduling times; It is the set of all thermal power generating units on the transmission grid side; is the power generation cost coefficient of thermal power units; It is the benchmark output of thermal power units; is the reserve cost coefficient on the transmission network side; and They are the benchmark upward spinning reserve capacity and the benchmark downward spinning reserve capacity on the transmission network side respectively.

[0037] Step 2-1-2: Affine decision criterion constraint. In step 1, it has been derived that the impact of typhoon parameter uncertainty on offshore wind farms can be converted into a deviation value compared to the output of offshore wind farms in the forecast scenario. It is assumed that the transmission network operator will allocate the adjustment power that the transmission-side thermal power units and the distribution network interconnected with the transmission side need to bear according to the following affine decision criterion: (13) in, is the set of all thermal power generating units on the transmission grid side, It is the set of all distribution network tie lines interconnected with the transmission network; The actual output of the thermal power unit; It is the actual interactive active power for transmission and distribution; Predict interactive active power for transmission and distribution; The actual output of the thermal power unit; is the power adjustment factor of the thermal power unit; is the power adjustment factor of the distribution network; is a random variable that characterizes the output deviation of offshore wind farms.

[0038] Step 2-1-3: Transmission grid power balance constraints: (14) in, is the set of all nodes of the transmission network; It is the set of all distribution network tie lines interconnected with the transmission network; Contribute to forecasting of offshore wind farms; It is the load on the transmission grid side; is the benchmark interaction power of the distribution network tie line.

[0039] Step 2-1-4: Thermal power unit operation constraints: (15) (16) in, is the probability that the inequality in brackets holds true; is a risk factor; and They are the upward spin reserve capacity and downward spin reserve capacity on the transmission grid side respectively; and The upper limit of the upward and downward reserve capacity on the transmission network side; and They are the upper and lower limits of the output of thermal power units respectively.

[0040] Step 2-1-5: Thermal power unit climbing constraints: (17) In the above formula, It is the upper limit of the climbing capacity of thermal power units.

[0041] Step 2-1-6: Transmission and distribution tie line exchange power constraints: (18) In the above formula, is the upper bound of the interconnection line exchange power.

[0042] Step 2-1-7: Transmission grid side line flow constraints: (19) in, It is the set of all lines on the transmission network side; is the upper limit of the power flow of the transmission line on the transmission network side; , , , They are the power transmission distribution factors of offshore wind farms, transmission grid side loads, thermal power units, and transmission and distribution interconnection lines.

[0043] Step 2-2: Define the set of all distribution networks interconnected with the upper transmission network as , then for any distribution network , a coordinated dispatching model on the distribution network side can be established, which includes objective functions and constraints; the details are as follows: Step 2-2-1: The objective function of the coordinated dispatching model on the distribution network side is to minimize the operating cost of the distribution network side. Its mathematical form is as follows: (20) in, It is the collection of all distributed power sources in the distribution network; It is the set of all gas generating units in the distribution network; It is the collection of all power-to-gas devices in the distribution network; It is the collection of all energy storage devices in the distribution network; is the cost coefficient of distributed generation; is the reserve capacity cost coefficient on the distribution network side; is the cost coefficient of the power-to-gas device; and are the charging and discharging cost coefficients of the energy storage device respectively; is the power generation cost coefficient of the gas unit; Provides benchmark active power output for distributed power sources; It is the benchmark active power of the power-to-gas device; and are the benchmark charging and discharging active powers of the energy storage device respectively; It is the benchmark active power output of the gas generating unit; and They are the benchmark upward spin reserve capacity and the benchmark downward spin reserve capacity on the distribution network side respectively.

[0044] Step 2-2-2: Affine decision criterion constraint. The distribution network operator will redistribute the adjustment power allocated in the transmission network dispatch to the distributed power sources, gas units, energy storage equipment and power-to-gas devices in the distribution network according to the affine decision criterion. Its mathematical form is as follows: (twenty one) (twenty two) in, Represents the adjustment power allocated to the distribution network in the transmission network side dispatch; , , , , They are the adjustment factors for distributed power sources, gas-fired units, energy storage device discharge, energy storage device charging, and power-to-gas devices.

[0045] Step 2-2-3: Distribution network flow constraints: (twenty three) in, For distribution network nodes The set of all child nodes of ; It is the set of all lines in the distribution network; It is the actual active power output of distributed power sources; It is the actual active power output of the gas unit; is the actual active power of the power-to-gas device; and are the actual charging and discharging active powers of the energy storage device respectively; is the active load on the distribution network side; Flow to the distribution network node Active power of For the distribution network node Flow to its child nodes Active power of It is the active power injected into the root node of the distribution network by the transmission and distribution tie line; exchanging power between the transmission network and the tie line of the distribution network; Flow to the distribution network node Reactive power; For the distribution network node Flow to its child nodes Reactive power; Provide reactive power for distributed electricity; Provide reactive power for gas generator sets; It is the reactive load on the distribution network side; For distribution network nodes The voltage amplitude of For distribution network nodes The voltage amplitude of the parent node; is the voltage amplitude of the root node of the distribution network; and Distribution network lines resistance and reactance.

[0046] Step 2-2-4: Distribution network node voltage constraints: (twenty four) in, and are the upper and lower bounds of the voltage amplitude of the distribution network nodes, respectively.

[0047] Step 2-2-5: Power-to-gas device operation constraints: (25) in, It is the maximum operating active power of the power-to-gas device.

[0048] Step 2-2-6: Distribution network side reserve capacity constraints: (26) (27) In the above formula, and They are the upward spin reserve capacity and downward spin reserve capacity on the distribution network side respectively; and It is the upper limit of the upward spin reserve capacity and downward spin reserve capacity on the distribution network side; and are the upper and lower bounds of the distributed generation output respectively; and They are the upper and lower limits of the gas unit output respectively.

[0049] Step 2-2-7: Energy storage equipment operation constraints: (28) (29) In the above formula, and They are the upper limits of the charging and discharging power of the energy storage device respectively; The amount of electricity stored for the energy storage device; The upper limit of the capacity of the energy storage device; is the conversion efficiency of the energy storage device.

[0050] Step 2-3: Establish a natural gas grid-side coordinated dispatch model, which includes objective functions and constraints; the details are as follows: Step 2-3-1: The objective function of the coordinated dispatching model on the natural gas grid side is to minimize the operating cost of the natural gas grid side. Its mathematical form is as follows: (30) In the above formula, A collection of gas wells in a natural gas grid; It is the set of all nodes in the natural gas grid; is the gas production cost coefficient of the natural gas well; is the gas production of the natural gas well; It is the gas consumption of gas unit.

[0051] Step 2-3-2: Gas unit coupling constraints: (31) In the above formula, is the gas consumption of the gas unit, is the first-order term of the gas consumption coefficient of the gas generator set; is the constant term of gas consumption coefficient of gas generator set; It is the adjustment power allocated to the gas unit in the worst scenario.

[0052] Step 2-3-3: Coupling constraints of the power-to-gas device: (32) In the above formula, Gas consumption of the power-to-gas device; is the conversion efficiency of the power-to-gas device.

[0053] Step 2-3-4: Natural gas flow balance constraints: (33) In the above formula, is the gas production of the natural gas well; Gas consumption of the power-to-gas device; is the gas consumption of the gas unit; and Respectively represent pipelines and Gas flow rate; is the natural gas node load.

[0054] Step 2-3-5: Steady-state gas flow equation for natural gas grid: (34) In the above formula, Represents a pipeline Weymouth characteristic parameters; and is the node pressure; For pipeline The upper limit of natural gas flow.

[0055] Step 2-3-6: Natural gas well production constraints: (35) In the above formula, It is the upper limit of gas production of natural gas wells.

[0056] Step 2-3-7: Node pressure constraint: (36) In the above formula, and are the upper and lower bounds of the node pressure, respectively.

[0057] Step 3: Based on the mixed sample average approximation, the three-layer coordinated scheduling model of transmission-distribution-gas with random variables is transformed into a solvable model, and it is rewritten into a deterministic continuous linear model to facilitate subsequent distributed solution. Including: Based on the mixed sample average approximation, the three-layer coordinated scheduling model of transmission-distribution-gas integrated energy system based on the scenario is transformed into a solvable model, including: the solvability transformation of the opportunity constraints containing random variables, based on each generated scenario, the sample average approximation is used to describe the opportunity constraints containing random variables, and a binary indicator variable is set. If the internal constraint is established in a certain scenario, the binary indicator variable takes the value of 1, otherwise it takes the value of 0, and the binary indicator variable is relaxed to be in the interval Continuous variables, and introduce mixed inequalities to tighten the feasible domain expanded by relaxation, and the general form of chance constraints is transformed into a set of continuous linear constraints without random variables; the scenario-based transmission-distribution-gas integrated energy system three-layer coordinated scheduling model after solvability transformation is expressed as a linear objective function and a set of linear constraints. The details are as follows: Step 3-1: For the chance constraints containing random variables, such as formulas (16)-(19), (24), (25), (27)-(29), perform solvability transformation and rewrite them into deterministic continuous linear constraints; the details are as follows: Step 3-1-1: The general form of a chance constraint involving random variables is as follows: (37) In the above formula, and Respectively represent the active / reactive power decision vectors of the controllable units. and They represent the decision vectors for increasing / decreasing the spare capacity respectively. Characterizes linear terms involving random variables. is the decision vector , , , The linear function corresponding to the corresponding joint chance constraint is constructed.

[0058] Step 3-1-2: Based on the generated in step 1 In this scenario, we use the sample average approximation to rewrite formula (37) into the following form: (38) In the above formula, is the upper bound of the uncertainty term; is a non-negative parameter; For the scene The probability of occurrence, and the condition is met ; is a slack variable; formula (a) in formula (39) is used to determine the scenario Internal constraints Is it true, if it is in the scene If the following holds, then the binary indicator variable The value is 1; otherwise The value is 0. Formula (b) is used to ensure that the maximum total probability of scenarios violating internal constraints does not exceed the risk factor , that is, formula (b) can ensure At confidence level Established.

[0059] Step 3-1-3: Binary indicator variables Relaxation is in the interval Continuous variables, and introduce mixed inequalities to tighten the relaxation For the expanded feasible domain, the general form of the chance constraint, such as formula (37), can be transformed into a set of continuous linear constraints without random variables as shown below: (39) In the above formula, is the mixed inequality introduced; represents the feasible domain after tightening; , , To follow the sequence Reordered non-negative parameters , and satisfies and ; A continuous indicator variable representing a mixed inequality.

[0060] Step 3-2: Build a compact model; details are as follows: Step 3-2-1: After the solvability transformation, the three-layer coordinated scheduling model of the transmission-distribution-gas integrated energy system can be expressed as the following linear objective function (40) and a set of linear constraints (41)-(43): (40) (41) (42) (43) In the above formula, , , They are the compact objective functions of the transmission grid side, distribution grid side, and natural gas grid side respectively; , , They represent the internal decision variables on the transmission network side, distribution network side, and natural gas network side respectively; is the transmission network and the index is The boundary variables between the distribution networks, The index is The boundary variables between the distribution network and the local gas grid. , , , , , , , , , , , , , , , , , , , They represent the coefficient matrices corresponding to the constraints.

[0061] Step 4: Use an efficient distributed collaborative algorithm based on multi-parameter planning theory and a critical domain acceleration search algorithm based on global optimal solution estimation to solve the three-layer transmission-distribution-gas collaborative scheduling model after solvability transformation, including: using an efficient distributed collaborative algorithm based on multi-parameter planning theory and a critical domain acceleration search algorithm based on global optimal solution estimation to solve the scenario-based transmission-distribution-gas integrated energy system three-layer collaborative scheduling model after solvability transformation, divided into two solution levels: the transmission level and the distribution level. The transmission level is an iterative solution between the transmission network and the distribution network, and the distribution level is an iterative solution between the distribution network and the local natural gas network. The iterative solution at the transmission level must be performed after convergence at the distribution level. Then calculate the critical domain and local objective function. The calculation method of the critical domain and local objective function at the transmission level is consistent with that at the distribution level. If the iterative solution does not converge to the global optimal solution of the subproblem, the obtained local optimal solution and optimal target value are used as the data set and the optimal target value respectively. The data set is used to fit the global objective function, and the estimated value of the global optimal solution is calculated based on the fitted value. Set the step vector, search for the next critical domain based on the step vector, until the subproblem is infeasible or the local optimal value begins to increase, return to the last feasible critical domain, and switch to a one-by-one search method to find the global optimal solution. The details are as follows: Step 4-1: Use an efficient distributed collaborative algorithm based on multi-parameter planning theory to solve the three-layer collaborative scheduling model of the transmission-distribution-gas integrated energy system after solvability transformation, which can be divided into two solution levels: the transmission level and the distribution level. The transmission level is the iterative solution between the transmission network and the distribution network, and the distribution level is the iterative solution between the distribution network and the local natural gas network. The iterative solution at the transmission level must be performed after the distribution level converges. The details are as follows: Step 4-1-1: Distributed collaborative solution at the distribution level. Note that this distributed collaborative solution process is for all the sets The distribution network is carried out in parallel; the specific steps are as follows: Step 4-1-1-1: For any Initialize the number of iterations for the indexed distribution network .

[0062] Step 4-1-1-2: Input the power distribution network Initial values ​​of boundary variables with the local gas grid , note that in the subsequent calculations Treated as parameters rather than variables.

[0063] Step 4-1-1-3: Distribution network The operator will bound the variable initial value Sent to the local gas grid operator.

[0064] Step 4-1-1-4: Local gas grid operator Solve the gas grid side problem as shown in equation (45), where For the The step vector for the iteration is used to remove the local gas grid operator from the distribution network The bounds variable received by the operator Move to the next critical region. Sometimes ,when hour, The calculation method of is shown in step 4-3.

[0065] (44) in, is the cost coefficient of the decision variable on the gas grid side.

[0066] Step 4-1-1-5: If the grid-side subproblem (44) is feasible, the local gas grid operator calculates the critical region at the distribution level and the local objective function and sends it back to the distribution network operator. and the local objective function The specific calculation method of is shown in step 4-2.

[0067] Step 4-1-1-6: Distribution network The operator receives the distribution level critical domain and the local objective function Calculate the distribution side main problem (45) to obtain the updated boundary variables : (45) Step 4-1-1-7: Determine convergence. If ,in is the convergence criterion, which is a small positive constant. If the iterative solution at the distribution level converges, the iterative solution at the transmission level can be performed. Otherwise, update , distribution network The operator will Send it to the local natural gas grid operator and return to step 4-1-1-4 to continue iterating until convergence.

[0068] Step 4-1-2: Distributed collaborative solution at the transmission level; the specific steps are as follows: Step 4-1-2-1: Initialize the number of transmission level iterations .

[0069] Step 4-1-2-2: Input the initial values ​​of the boundary variables between the transmission network and the distribution network , note that in the subsequent calculations Treated as parameters rather than variables.

[0070] Step 4-1-2-3: The transmission network operator sets the initial values ​​of the boundary variables Sent to each distribution network operator.

[0071] Step 4-1-2-4: For the items belonging to the set All distribution networks, their respective operators The distribution network sub-problem is solved in parallel as shown in Equation (46), where For the The step vector of the iteration is used to convert the boundary variables received by each distribution network operator from the transmission network operator into Move to the next critical region. Sometimes ,when hour, The calculation method of is shown in step 4-3.

[0072] (46) in, and are the distribution-gas boundary variables and local objective functions that converge in the iterative solution at the distribution level, both of which are steady values.

[0073] Step 4-1-2-5: If the distribution network side subproblem 46 is feasible, each distribution network operator calculates the critical region at the transmission level and the local objective function and sends it back to the transmission grid operator. and the local objective function The specific calculation method of is shown in step 4-2.

[0074] Step 4-1-2-6: The transmission network operator receives the critical domain and the local objective function Calculate the distribution side main problem (47) to obtain the updated boundary variables .

[0075] (47) Step 4-1-2-7: Determine convergence. If , then the iterative solution at the transmission level converges, and the global optimal solution has been found; otherwise, update , the transmission grid operator will Send it to each distribution network operator and return to step 4-1-2-4 to continue iterating until convergence.

[0076] Step 4-2: Calculate the critical region and local objective function. Note that the calculation method of the critical region and local objective function at the transmission level is the same as that at the distribution level; the details are as follows: Step 4-2-1: Define the Lagrangian function of the gas grid side subproblem (44) at the distribution level as follows: (48) In the above formula, is the dual multiplier of the equality constraint, is the dual multiplier of the inequality constraint.

[0077] Step 4-2-2: After solving the gas network side sub-problem (44), according to the optimal solution of this sub-problem The effective and ineffective constraints can be expressed as follows: (49) (50) In the above formula, and are the sets of active constraints and inequality constraints respectively.

[0078] Step 4-2-3: To solve the Karush-Kuhn-Tucker equations (49), the coefficient matrix on the left side of equation (49) can be inverted block by block and then multiplied by the right side of the equation. The calculation result is as follows: (51) In the above formula, and The sub-block inverse matrix after block inversion of the coefficient matrix.

[0079] Step 4-2-4: From formula (51), we can see that the optimal solution of the sub-problem is It can be expressed as Function: (52) Step 4-2-5: Although It will change during the iterative solution process, that is, although It is considered as a parameter in solving the sub-problem on the gas grid side, but it is considered as a variable in the complete solution of the coordinated dispatch model at the distribution level. , at this time, as long as the combination of the effective constraints and ineffective constraints in (50) remains unchanged, equation (52) always holds. The feasible region of is called the critical region, which must satisfy the following conditions: the dual multiplier remains greater than zero, and the ineffective constraint remains less than zero. The mathematical form of the critical region can be expressed as: (53) (54) (55) Step 4-2-6: Under the limitation of , by substituting equation (52) into the objective function of the gas grid side subproblem (44), we can get the current critical region The local objective function is: (56) Step 4-3: Calculate the step vector. Note that the calculation method of the step vector at the transmission level is the same as that at the distribution level; the details are as follows: Step 4-3-1: Set the iteration threshold to , when the number of iterations hour, Can be set as the local objective function used in this iteration exist Take a small step forward in the negative gradient direction. The mathematical form is as follows: (57) In the above formula, is the step length, is the local objective function exist The negative gradient unit vector of .

[0080] Step 4-3-2: When the number of iterations If the iterative solution does not converge to the global optimal solution of the subproblem, the local optimal solution obtained is and the optimal target value As data sets and save.

[0081] Step 4-3-3: Utilize the dataset and Fitting the global objective function Here, the convex quadratic objective function shown below is used for fitting: (58) In the above formula, , , is a constant coefficient matrix.

[0082] Step 4-3-4: Calculate the estimated value of the global optimal solution according to formula (58) At this time, the step vector can be set as shown in the following formula: (59) In the above formula, is an appropriate step size parameter.

[0083] Step 4-3-5: Based on step vector Search for the next critical region. If the subproblem is feasible and the local optimal value decreases, and Update Add and , and return to step 4-3-2.

[0084] Step 4-3-6: Repeat steps 4-3-2 to 4-3-5 until the subproblem is infeasible or the local optimal value begins to increase, return to the last feasible critical region, and switch to the one-by-one search method of step 4-3-1 to find the global optimal solution.

[0085] Example 2 In one embodiment of the present disclosure, a distributed scheduling system for an integrated energy system taking into account typhoon uncertainty is provided, comprising: The scenario generation module is used to consider the impact of typhoon uncertainty on offshore wind power, build an offshore wind power operation model, and use the Monte Carlo simulation method to sample and generate typhoon scenarios and offshore wind power output scenarios based on the offshore wind power operation model; The model building module is used to establish a scenario-based three-layer coordinated dispatching model of the transmission-distribution-gas integrated energy system based on typhoon scenarios and offshore wind power output scenarios, taking into account the coordinated operation of the transmission network, distribution network and natural gas network. The scheduling solution module is used to transform the scenario-based transmission-distribution-gas integrated energy system three-layer collaborative scheduling model into a deterministic continuous linear model based on the mixed sample average approximation, and use the efficient distributed collaborative algorithm of multi-parameter planning theory and the critical domain acceleration search algorithm based on the global optimal solution estimation to solve the deterministic continuous linear model, and obtain the global optimal solution for the distributed scheduling of the integrated energy system under the typhoon scenario.

[0086] Example 3 In one embodiment of the present disclosure, a computer program product is provided, including a computer program, which, when executed by a processor, implements the distributed scheduling method of the integrated energy system taking into account typhoon uncertainty.

[0087] Example 4 In one embodiment of the present disclosure, a non-transitory computer-readable storage medium is provided, wherein the non-transitory computer-readable storage medium is used to store computer instructions. When the computer instructions are executed by a processor, the distributed scheduling method of the integrated energy system taking into account typhoon uncertainty is implemented.

[0088] Example 5 In one embodiment of the present disclosure, an electronic device is provided, including: a processor, a memory, and a computer program; wherein the processor is connected to the memory, the computer program is stored in the memory, and when the electronic device is running, the processor executes the computer program stored in the memory, so that the electronic device executes the distributed scheduling method of the integrated energy system taking into account typhoon uncertainty.

[0089] The present disclosure is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present disclosure. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of the processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0090] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1A process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0091] Although the above describes the specific implementation methods of the present disclosure in conjunction with the accompanying drawings, it is not intended to limit the scope of protection of the present disclosure. Technical personnel in the relevant field should understand that on the basis of the technical solution of the present disclosure, various modifications or variations that can be made by those skilled in the art without creative work are still within the scope of protection of the present disclosure.

Claims

1. A distributed scheduling method for an integrated energy system taking into account typhoon uncertainty, characterized in that: include: Considering the impact of typhoon uncertainty on offshore wind power, an offshore wind power operation model is constructed. Based on the offshore wind power operation model, the Monte Carlo simulation method is used to sample and generate typhoon scenarios and offshore wind power output scenarios; Based on typhoon scenarios and offshore wind power output scenarios, the coordinated operation of the transmission network, distribution network and natural gas network is considered, and a scenario-based three-layer coordinated dispatching model of the transmission-distribution-gas integrated energy system is established; Based on the mixed sample average approximation, the three-layer collaborative scheduling model of the scenario-based transmission-distribution-gas integrated energy system is transformed into a deterministic continuous linear model. The deterministic continuous linear model is solved using an efficient distributed collaborative algorithm based on multi-parameter planning theory and a critical domain accelerated search algorithm based on global optimal solution estimation to obtain the global optimal solution for the distributed scheduling of the integrated energy system under typhoon scenarios.

2. The distributed scheduling method for an integrated energy system taking into account typhoon uncertainty according to claim 1, characterized in that: Considering the impact of typhoon uncertainty on offshore wind power, the Holland vortex model is used to model the typhoon wind farm, and the typhoon parameter uncertainty is modeled. When the probability distribution satisfied by the typhoon's uncertain parameters is determined, the typhoon scenario is generated by Monte Carlo simulation. For any offshore wind farm, a wind speed-power active output model is constructed in a specific typhoon scenario. According to the relative distance between the typhoon and the offshore wind farm, the whole process of the offshore wind farm encountering a typhoon event is divided into three stages, namely, the typhoon approach stage, the typhoon transit stage and the typhoon departure stage. The wind power output vectors of the three stages are integrated to obtain the complete output of the offshore wind farm in the typhoon scenario.

3. The distributed scheduling method for an integrated energy system taking into account typhoon uncertainty according to claim 1, characterized in that: A scenario-based three-layer coordinated dispatching model of a transmission-distribution-gas integrated energy system is established, wherein the scenario-based three-layer coordinated dispatching model of a transmission-distribution-gas integrated energy system includes a transmission network side coordinated dispatching model, a distribution network side coordinated dispatching model and a natural gas network side coordinated dispatching model; wherein the objective function of the transmission network side coordinated dispatching model is to minimize the transmission network side operating cost, the objective function of the distribution network side coordinated dispatching model is to minimize the distribution network side operating cost, and the objective function of the natural gas network side coordinated dispatching model is to minimize the natural gas network side operating cost, and the constraints of each coordinated dispatching model are set respectively.

4. The distributed scheduling method for an integrated energy system taking into account typhoon uncertainty according to claim 1, characterized in that: Based on the hybrid sample average approximation, the scenario-based transmission-distribution-gas integrated energy system three-layer coordinated scheduling model is transformed into solvability, including: the solvability transformation of the opportunity constraints containing random variables is carried out, based on each generated scenario, the opportunity constraints containing random variables are described using the sample average approximation, and a binary indicator variable is set. If the internal constraint is established in a certain scenario, the binary indicator variable takes the value of 1, otherwise it takes the value of 0, and the binary indicator variable is relaxed to be in the interval , and introduces mixed inequalities to tighten the feasible domain expanded by relaxation. The general form of chance constraints is transformed into a set of continuous linear constraints without random variables. The scenario-based transmission-distribution-gas integrated energy system three-layer coordinated scheduling model after solvability transformation is expressed as a linear objective function and a set of linear constraints.

5. The distributed scheduling method for an integrated energy system taking into account typhoon uncertainty as claimed in claim 1, characterized in that: An efficient distributed collaborative algorithm based on multi-parameter planning theory and a critical domain acceleration search algorithm based on global optimal solution estimation are used to solve the scenario-based three-layer collaborative scheduling model of the transmission-distribution-gas integrated energy system after solvability transformation. The model is divided into two solution levels: the transmission level and the distribution level. The transmission level is an iterative solution between the transmission network and the distribution network, and the distribution level is an iterative solution between the distribution network and the local natural gas network. The iterative solution at the transmission level must be performed after convergence at the distribution level.

6. The distributed scheduling method for an integrated energy system taking into account typhoon uncertainty as claimed in claim 5, characterized in that: Calculate the critical domain and local objective function. The calculation method of the critical domain and local objective function at the transmission level is consistent with that at the distribution level. If the iterative solution does not converge to the global optimal solution of the subproblem, the obtained local optimal solution and optimal target value are used as the data set and the optimal target value respectively. The data set is used to fit the global objective function, and the estimated value of the global optimal solution is calculated based on the fitted value. Set the step vector and search the next critical region based on the step vector until the subproblem is infeasible or the local optimal value begins to increase, return to the last feasible critical region, and switch to a one-by-one search mode to find the global optimal solution.

7. A distributed dispatching system for an integrated energy system taking into account typhoon uncertainty, characterized in that: include: The scenario generation module is used to consider the impact of typhoon uncertainty on offshore wind power, build an offshore wind power operation model, and use the Monte Carlo simulation method to sample and generate typhoon scenarios and offshore wind power output scenarios based on the offshore wind power operation model; The model building module is used to establish a scenario-based three-layer coordinated dispatching model of the transmission-distribution-gas integrated energy system based on typhoon scenarios and offshore wind power output scenarios, taking into account the coordinated operation of the transmission network, distribution network and natural gas network. The scheduling solution module is used to transform the scenario-based transmission-distribution-gas integrated energy system three-layer collaborative scheduling model into a deterministic continuous linear model based on the mixed sample average approximation, and use the efficient distributed collaborative algorithm of multi-parameter planning theory and the critical domain acceleration search algorithm based on the global optimal solution estimation to solve the deterministic continuous linear model, and obtain the global optimal solution for the distributed scheduling of the integrated energy system under the typhoon scenario.

8. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the distributed scheduling method for an integrated energy system taking into account typhoon uncertainty as described in any one of claims 1 to 6 is implemented.

9. A non-transitory computer-readable storage medium, characterized in that: The non-transitory computer-readable storage medium is used to store computer instructions. When the computer instructions are executed by the processor, the distributed scheduling method of the integrated energy system taking into account typhoon uncertainty as described in any one of claims 1-6 is implemented.

10. An electronic device, characterized in that: include: A processor, a memory and a computer program; wherein the processor is connected to the memory, the computer program is stored in the memory, and when the electronic device is running, the processor executes the computer program stored in the memory so that the electronic device executes the distributed scheduling method for the integrated energy system taking into account the typhoon uncertainty as described in any one of claims 1-6.

Images