CN112434402A - Interval practical safety domain modeling method - Google Patents

Interval practical safety domain modeling method Download PDF

Info

Publication number
CN112434402A
CN112434402A CN202011139944.9A CN202011139944A CN112434402A CN 112434402 A CN112434402 A CN 112434402A CN 202011139944 A CN202011139944 A CN 202011139944A CN 112434402 A CN112434402 A CN 112434402A
Authority
CN
China
Prior art keywords
energy
supply capacity
representing
pipeline
energy supply
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011139944.9A
Other languages
Chinese (zh)
Other versions
CN112434402B (en
Inventor
王丹
李思源
贾宏杰
雷杨
果营
李家熙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN202011139944.9A priority Critical patent/CN112434402B/en
Publication of CN112434402A publication Critical patent/CN112434402A/en
Application granted granted Critical
Publication of CN112434402B publication Critical patent/CN112434402B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector

Abstract

The invention discloses a regional practical safety domain modeling method, which is based on a regional comprehensive energy system containing renewable energy, is used for processing by considering the uncertainty of the renewable energy and adopting a regional mathematical method, combines engineering practice, establishes a regional practical safety domain model of the regional comprehensive energy system containing the renewable energy based on the N-1 safety criterion of the regional comprehensive energy system and the pipeline segment load at the output side of an energy junction, further establishes a regional practical safety boundary model, and provides a regional practical safety distance model and a regional maximum energy supply capacity model. On the basis, a system dimension reduction observation method and a full-dimension observation method are established, and the safety of the regional comprehensive energy system containing renewable energy can be accurately and quickly analyzed.

Description

Interval practical safety domain modeling method
Technical Field
The invention relates to the field of comprehensive energy systems, in particular to a regional practical safety domain modeling method for a regional comprehensive energy system containing renewable energy sources.
Background
The energy is a power source for promoting national economic development and guaranteeing human survival and production, and how to realize efficient, green and sustainable supply of the energy is the key point of attention of all countries in the world. The energy Internet, the Internet plus, the intelligent energy and other concepts are put forward, and the importance of the cross-integration and the diversification common development of various energy sources in different subject fields is emphasized. Meanwhile, key energy conversion technologies represented by combined cooling heating and power supply, gas boilers and the like are developed rapidly, and a technical basis is laid for the provision of a novel energy supply mode. Under the background, the comprehensive energy system breaks through independent planning, design and operation of the original single energy system. The integrated energy system is an important physical carrier of an energy internet, and relates to a plurality of links such as energy production, transmission, distribution, conversion, consumption and the like, and a regional integrated energy system (rees) is a typical application scenario among the links.
The comprehensive energy system takes various heterogeneous energy subsystems such as an electric power system, a natural gas system, a regional heat supply (medium is steam or water) system, a water supply system and the like as backbones, and has the characteristics of wide source load distribution, diversified energy consumption characteristics, intensive energy interaction, strong physical information coupling and the like. On one hand, the multi-energy coupling can promote the cascade utilization of energy, green low-carbon energy utilization and the consumption of renewable energy; on the other hand, however, due to the different structures, characteristics and compositions of the energy subsystems, the structure of the comprehensive energy system is also more complicated due to the multi-energy coupling, and the local disturbance (load fluctuation or element fault) of a certain element is transferred to the whole multi-energy flow network through the coupling component. Therefore, the safety analysis of the integrated energy system is more complex, and due to the mutual influence of multiple energy flows, the normally operating system may be in an unsafe state after a fault occurs in a certain element, and research on the safety analysis of the integrated energy system is urgently needed.
The safety and reliability of the energy system are the most basic requirements of the operation of the energy system, are the basis of research on aspects such as planning, operation, transaction and the like, and are one of the important research directions of the comprehensive energy system. At present, the research on the safety of the comprehensive energy system is mostly based on a traditional point-by-point method, the safety state is judged by locally limiting operating points through point-by-point simulation check, the obtained safety information is one-sided, the calculated amount is large, the consumed time is long, and the method is not suitable for online safety analysis. The 'domain' method is an effective tool for solving the problems, the safety boundary of the system can be effectively and conveniently observed by utilizing the safety domain, various safety information such as the safety state, the safety distance, the safety margin, the adjustment direction and the like of the system can be obtained according to the relative position of the working point in the safety domain, the safety evaluation efficiency is greatly improved, and the solution of various optimization problems related to safety is simplified. The research on the security domain in the regional integrated energy system has a certain foundation, but the existing research fails to consider the uncertain influence of renewable energy, and meanwhile, the full-dimensional observation on the security domain is not researched, so that the modeling research on the security domain of the regional integrated energy system containing renewable energy is necessary, and the purpose of providing a convenient and efficient analysis tool for ensuring the safe and reliable energy supply of the regional integrated energy system containing renewable energy is achieved.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provides a regional practical safety domain modeling method of a regional integrated energy system containing renewable energy, which considers the uncertainty of the renewable energy and adopts a regional mathematical method for processing, combines engineering practice, and establishes a regional practical safety domain definition and model of the regional integrated energy system containing the renewable energy based on the N-1 safety criterion of the regional integrated energy system and the pipeline section load on the output side of an energy junction, and is described in detail as follows:
the purpose of the invention is realized by the following technical scheme:
an interval practical safety domain modeling method based on a regional comprehensive energy system containing renewable energy sources comprises the following steps:
considering the uncertainty of renewable energy sources, processing by adopting a region mathematical method, combining engineering practice, and establishing a region practical safety domain model of the regional integrated energy source system containing the renewable energy sources based on the N-1 safety criterion of the regional integrated energy source system and the pipeline segment load on the output side of the energy junction;
establishing a regional comprehensive energy system interval practical safety boundary model containing renewable energy sources based on the regional comprehensive energy system interval practical safety boundary model containing renewable energy sources;
establishing a regional comprehensive energy system interval maximum energy supply capacity model containing renewable energy sources and a simulation solving method based on the provided safety boundary model;
establishing a regional comprehensive energy system interval dimension reduction projection observation method containing renewable energy sources based on the maximum energy supply working point solved by the maximum energy supply capacity model and according to a safe boundary simulation fitting solving method;
establishing a practical safe distance model and a simulation solving method of a regional comprehensive energy system interval containing renewable energy sources based on the proposed safe boundary model;
establishing a regional comprehensive energy system interval full-dimensional observation method containing renewable energy sources based on the proposed safe distance model;
further, the regional integrated energy system interval practical safety domain containing renewable energy is defined as:
considering the safety inequality constraint and the energy balance constraint of the N-1 interval, and collecting all working points which can meet the safety constraint in the regional comprehensive energy system containing renewable energy;
the working point refers to the minimum set of state variables representing the system safety when the regional comprehensive energy system containing renewable energy is in normal operation; the state variable is selected as the load of the pipeline section on the output side of the energy hub, and the load comprises the load of the pipeline section on the output side of energy coupling equipment (a gas boiler, cogeneration, an electric boiler, an air conditioning unit and a gas turbine) and associated components (a transformer, a compressor, a heat exchanger and energy storage facilities) in the energy hub.
Further, the regional integrated energy system interval practical security domain model containing renewable energy sources is as follows:
Figure BDA0002737910820000031
wherein, L represents a working point vector; h (L) 0 represents a regional comprehensive energy system multipotency flow energy balance equation; λ represents a pipeline scaling coefficient vector;
Figure BDA0002737910820000032
respectively representing an upper limit energy supply capacity matrix and a lower limit energy supply capacity matrix;
Figure BDA0002737910820000033
Figure BDA0002737910820000034
respectively representing the upper and lower boundaries of the upper and lower limit energy supply capacities.
The safety inequality constraint of the N-1 interval comprises the following steps:
(1) safety constraint for transition between output side of key pipeline and upper limit N-1 of key equipment
Figure BDA0002737910820000035
Figure BDA0002737910820000036
Lm+Ln≤Cn
Figure BDA0002737910820000037
In the formula, Ln、LmRepresents the load carried by the pipelines n and m; hjRepresents the load carried by the device j; cn、CjRepresents the rated capacity of pipeline n and equipment j; cpro,zRepresenting the probabilistic capacity of the new energy device z; cRE,zRepresenting the rated capacity of the new energy device z; s represents the total number of new energy devices accessed;
Figure BDA0002737910820000038
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with distributed renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure BDA0002737910820000039
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of the key pipeline with the distributed renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure BDA00027379108200000310
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with centralized renewable energy sources corresponding to the y load transfer scheme under the x fault.
(2) Safety constraint for transition between output side of key pipeline and lower limit N-1 of key equipment
Figure BDA00027379108200000311
PC+PCP≤Ln+Lm
In the formula (I), the compound is shown in the specification,
Figure BDA00027379108200000312
respectively representing the upper limit and the lower limit of the minimum energy supply capacity of the renewable energy source corresponding to the y load transfer scheme; l isn、LmRepresents the load carried by the pipelines n and m; cRE,zRepresenting the rated capacity of the new energy device z; pC、PCPRepresenting the energy consumption of the compressor and the heat pump. s represents the total number of new energy devices accessed;
further, the regional integrated energy system interval practical safety boundary model containing renewable energy sources is as follows:
Figure BDA0002737910820000041
in the formula (I), the compound is shown in the specification,
Figure BDA0002737910820000042
representing the upper and lower limits of the practical safety boundary corresponding to the pipeline m;
Figure BDA0002737910820000043
representing the upper and lower limits of the practical safety boundary corresponding to the pipeline m; lambda [ alpha ]kIndicating that pipeline k corresponds to a pipeline scale factor. L iskRepresents the load carried by the pipeline k;
Figure BDA0002737910820000044
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure BDA0002737910820000045
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of the key pipeline with the distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure BDA0002737910820000046
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with centralized renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure BDA0002737910820000047
respectively representing the upper limit and the lower limit of the minimum energy supply capacity of the renewable energy source corresponding to the y load transfer scheme; pC、PCPRepresenting the energy consumption of the compressor and the heat pump.
Further, the model of the maximum energy supply capacity of the renewable energy-containing regional integrated energy system interval is as follows:
Figure BDA0002737910820000048
Figure BDA0002737910820000049
in the formula (I), the compound is shown in the specification,
Figure BDA0002737910820000051
respectively representing the upper and lower limits of the maximum energy supply capacity. L represents a working point vector; h (L) 0 represents a regional comprehensive energy system multipotency flow energy balance equation; λ represents a pipeline scaling coefficient vector;
Figure BDA0002737910820000052
Figure BDA0002737910820000053
upper and lower boundaries representing an upper limit energy supply capacity and a lower limit energy supply capacity, respectively;
the simulation solving method comprises the following steps:
(1) initializing a primal-dual interior point method based on an interval practical safety domain of a regional comprehensive energy system containing renewable energy, and respectively acquiring initial working points and constraints corresponding to upper and lower limits of the maximum energy supply capacity by using the above formula;
(2) decoupling and calculating the power flow based on the fact that the upper formula is converted into a Lagrange function by a source-dual interior point method, calculating KKT conditions corresponding to the upper limit and the lower limit of the maximum energy supply capacity by the upper formula, and solving the disturbed KKT conditions by a Newton method;
(3) and outputting the upper and lower limits of the maximum energy supply capacity and the corresponding working points.
Further, the interval dimension reduction projection observation method of the regional comprehensive energy system containing renewable energy comprises the following steps:
and solving an interval practical safety boundary in a 2-dimensional or 3-dimensional space by using a simulation fitting method, and realizing interval dimension reduction projection observation by drawing the boundary. Taking a practical security domain in a three-dimensional interval as an example, the method specifically comprises the following steps:
(1): solving the upper and lower limits of the maximum energy supply capacity and corresponding working points;
(2): based onITSCA critical point array is obtained. Selecting any one group of multi-energy pipeline outlet load combination Lb＝(Lm,Ln,Lo) As a free variable, the remaining M-3 variables are fixed toITSCPipeline section loading when distributed; order (L)m,Ln) Are limited to upper limits by step sizes DeltaL respectively
Figure BDA0002737910820000054
Approaching until reaching the upper limit of the working point, solving the working point to satisfy the constraint of (L)m+Ln,+Lo) Maximum time of LoAnd load distribution, and setting the operating point Lb＝(Lm,Ln,Lo) Record to the critical point array B1Performing the following steps;
(3): based on
Figure BDA0002737910820000055
A critical point array is obtained. Selecting any one group of multi-energy pipeline outlet load combination Lb＝(Lm,Ln,Lo) As a free variable, the remaining M-3 variables are fixed to
Figure BDA0002737910820000056
Pipeline section loading when distributed; order (L)m,Ln) Are limited to upper limits by step sizes DeltaL respectively
Figure BDA0002737910820000057
Approaching until reaching the upper limit of the working point, solving the working point to satisfy the constraint of (L)m+Ln,+Lo) Maximum time of LoAnd load distribution, and setting the operating point Lb＝(Lm,Ln,Lo) Record to the critical point array B2Performing the following steps;
(4) the method comprises the following steps: fitting critical point array B by using least square method1And B2And the interval dimension reduction projection observation of the interval practical security domain is realized.
Further, the regional integrated energy system interval practical safety distance model containing renewable energy sources is as follows:
Figure BDA0002737910820000061
in the formula (I), the compound is shown in the specification,
Figure BDA0002737910820000062
representing the upper and lower limits of practical safety distance corresponding to the pipeline m;
Figure BDA0002737910820000063
representing the upper and lower limits of practical safety distance corresponding to the pipeline m; lambda [ alpha ]kRepresenting pipeline k corresponding to the pipeline scale factor;
Figure BDA0002737910820000064
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure BDA0002737910820000065
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of the key pipeline with the distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure BDA0002737910820000066
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with centralized renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure BDA0002737910820000067
respectively representing the upper limit and the lower limit of the minimum energy supply capacity of the renewable energy source corresponding to the y load transfer scheme; pC、PCPRepresenting the energy consumption of the compressor and the heat pump.
The regional comprehensive energy system interval practical safe distance simulation solving model comprises the following steps:
Figure BDA0002737910820000068
Figure BDA0002737910820000069
Figure BDA00027379108200000610
Figure BDA0002737910820000071
Figure BDA0002737910820000072
Figure BDA0002737910820000073
Figure BDA0002737910820000074
Figure BDA0002737910820000075
in the formula, Lq,cur、Lq,optRepresenting the load of the q-th pipeline corresponding to the current working point and the optimal working point; phimRepresenting a pipeline segment set corresponding to the constraint of the lower limit safety inequality on the mth pipeline;
Figure BDA0002737910820000076
IPSB m,u(Lq,opt) Representing the upper and lower bounds of the upper safety margin, with the argument Lq,opt；
Figure BDA0002737910820000077
IPSB m,u(Lq,cur) Representing the upper and lower bounds of the upper safety margin, with the argument Lq,cur；
Figure BDA0002737910820000078
IPSB m,l(Lq,opt) Representing the upper and lower bounds of the lower safety margin, with the argument Lq,opt；
Figure BDA0002737910820000079
IPSB m,l(Lq,cur) Representing the upper and lower bounds of the lower safety margin, with the argument Lq,cur. L represents a working point vector; h (L) 0 represents a regional comprehensive energy system multipotency flow energy balance equation;
the practical safety distance between the zones is solved by utilizing a primal-dual interior point method, and the practical safety distance between the zones corresponding to the No. m pipeline is taken as an example, and the practical safety distance between the zones is specifically:
(1) initializing a primal-dual interior point method based on an interval practical safety domain of a regional comprehensive energy system containing renewable energy, and respectively acquiring initial working points and constraints corresponding to upper and lower limits of maximum energy supply capacity by utilizing the 1 st and 2 nd formulas and the 3 rd and 4 th formulas;
(2) the method comprises the steps of converting an upper formula into a Lagrange function based on a source-dual interior point method, decoupling and calculating a power flow, calculating KKT conditions corresponding to upper and lower limits of practical safety distance by using the upper formula, and solving the disturbed KKT conditions by using a Newton method;
(3) the upper and lower limits of the practical safe distance are output.
Further, the regional comprehensive energy system interval full-dimensional observation method containing renewable energy comprises the following steps of firstly modeling an interval full-dimensional radius:
Figure BDA0002737910820000081
in the formula (I), the compound is shown in the specification,
Figure BDA0002737910820000082
representing the upper and lower limits of the upper full-dimensional radius corresponding to the pipeline m;
Figure BDA0002737910820000083
representing the upper and lower limits of the lower full-dimensional radius for line m.
Figure BDA0002737910820000084
Respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure BDA0002737910820000085
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of the key pipeline with the distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure BDA0002737910820000086
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with centralized renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure BDA0002737910820000087
respectively representing the upper limit and the lower limit of the minimum energy supply capacity of the renewable energy source corresponding to the y load transfer scheme; pC、PCPRepresenting the energy consumption of the compressor and the heat pump.
And solving the full-dimensional radius of the interval by using an interior point method, and drawing a full-dimensional radius boundary in a radar map form to realize full-dimensional observation of the interval.
Compared with the prior art, the technical scheme of the invention has the following beneficial effects:
1. compared with the existing regional comprehensive energy system security domain modeling method, the regional comprehensive energy system security domain modeling method considers the influence of renewable energy uncertainty, adopts a regional mathematics method to process the uncertainty, establishes a regional practical security domain model, a regional practical security boundary, a safe distance model and a regional maximum energy supply capacity model, can quantify the influence of renewable energy on the security domain, and obtains the practical region boundary of the security domain.
2. The invention provides a visual observation method of a security domain under a full-dimensional visual angle and a dimensionality reduction visual angle, and a 2-dimensional or 3-dimensional visual interval practical security domain of a system can be obtained through an interval dimensionality reduction projection observation method; by the interval full-dimensional observation method, the radar graph type system full-dimensional visual security domain can be obtained, and the system full-dimensional characteristics are represented.
3. The invention has wide application prospect, the interval dimension reduction projection observation method can be directly applied to safety control, and the visual characteristics of the method are helpful for operating personnel to visually determine the control scheme by observing whether the working point is in the area and the relative distance. The interval full-dimensional observation method can be applied to system planning research, and a planning scheme with good safety, new energy influence degree and investment sum is provided by maximizing relevant full-dimensional characteristics.
Drawings
FIG. 1 is a flow chart of solving a maximum energy supply capacity of a regional integrated energy system containing renewable energy;
FIG. 2 is a flowchart illustrating a regional integrated energy system interval practicalization safety boundary solution process with renewable energy;
FIG. 3 is a flowchart illustrating a process of solving a practical safe distance between regional integrated energy systems with renewable energy;
FIG. 4 is a schematic diagram of a test example topology and key equipment of a regional integrated energy system with renewable energy;
FIG. 5 is based on (L)8,L1,L2)The interval practical safety domain reduced dimension projection observation;
FIG. 6 is a full-dimensional observation of an interval practicalized security domain;
Detailed Description
The invention is described in further detail below with reference to the figures and specific examples. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The embodiment of the invention provides a method for modeling an interval practical safety domain oriented to a regional integrated energy system containing renewable energy, and the method is described in detail in the following with reference to fig. 1 to 3:
101: considering the uncertainty of renewable energy sources, processing by adopting a region mathematical method, combining engineering practice, and establishing a region practical safety domain model of the regional integrated energy source system containing the renewable energy sources based on the N-1 safety criterion of the regional integrated energy source system and the pipeline segment load on the output side of the energy junction;
102: establishing a regional comprehensive energy system interval practical safety boundary model containing renewable energy sources based on the proposed model;
103: establishing a regional comprehensive energy system interval maximum energy supply capacity model containing renewable energy sources and a simulation solving method based on the provided safety boundary model;
104: establishing a regional comprehensive energy system interval dimension reduction projection observation method containing renewable energy sources based on the maximum energy supply working point solved by the maximum energy supply capacity model and according to a safe boundary simulation fitting solving method;
105: establishing a practical safe distance model and a simulation solving method of a regional comprehensive energy system interval containing renewable energy sources based on the proposed safe boundary model;
106: establishing a regional comprehensive energy system interval full-dimensional observation method containing renewable energy sources based on the proposed safe distance model;
the scheme of example 1 is further described below with reference to fig. 1-3, calculation formulas, and examples, and is described in detail below:
with respect to step 101: considering the uncertainty of renewable energy sources, processing by adopting a region mathematical method, combining engineering practice, and establishing a region practical safety domain model of the regional integrated energy source system containing the renewable energy sources based on the N-1 safety criterion of the regional integrated energy source system and the pipeline segment load on the output side of the energy junction;
considering the safety inequality constraint and the energy balance constraint of the N-1 interval, and collecting all working points which can meet the safety constraint in the regional comprehensive energy system containing renewable energy; the working point refers to the minimum set of state variables representing the system safety when the regional comprehensive energy system containing renewable energy is in normal operation; the state variable is selected as the load of the pipeline section on the output side of the energy hub, and comprises energy coupling equipment (a gas boiler, cogeneration, an electric boiler and an air conditioning unit) in the energy hubAnd gas turbine) and associated components (transformers, compressors, heat exchangers, and energy storage facilities). Assuming that there are M pipeline segments in the energy hub, which are numbered 1,2, …, M, …, M, the operating point can be expressed as a vector L ═ L in euclidean space1,L2,…,Lm,…,LMIn which L ismThe pipeline section load of the No. m pipeline is obtained.
The regional comprehensive energy system interval practical safety domain model containing renewable energy sources is as follows:
Figure BDA0002737910820000101
wherein, L represents a working point vector; h (L) 0 represents a regional comprehensive energy system multipotency flow energy balance equation; λ represents a pipeline scaling coefficient vector;
Figure BDA0002737910820000102
respectively representing an upper limit energy supply capacity matrix and a lower limit energy supply capacity matrix;
Figure BDA0002737910820000103
Figure BDA0002737910820000104
respectively representing the upper and lower boundaries of the upper and lower limit energy supply capacities.
The safety inequality constraint of the N-1 interval comprises the following steps:
(1) safety constraint for transition between output side of key pipeline and upper limit N-1 of key equipment
When the load is switched to the belt, the renewable energy equipment is accessed in a distributed mode, and the safety inequality constraint of the corresponding interval is as follows:
Figure BDA0002737910820000105
Figure BDA0002737910820000106
when the load is switched to the belt, the renewable energy equipment is accessed in a centralized mode, and the safety inequality constraint of the corresponding interval is as follows:
Lm+Ln≤Cn (4)
Figure BDA0002737910820000107
in the formula, Ln、LmRepresents the load carried by the pipelines n and m; hjRepresents the load carried by the device j; cn、CjRepresents the rated capacity of pipeline n and equipment j; cpro,zRepresenting the probabilistic capacity of the new energy device z; cRE,zRepresenting the rated capacity of the new energy device z; s represents the total number of new energy devices accessed;
Figure BDA0002737910820000111
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with distributed renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure BDA0002737910820000112
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of the key pipeline with the distributed renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure BDA0002737910820000113
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with centralized renewable energy sources corresponding to the y load transfer scheme under the x fault.
(2) Safety constraint for transition between output side of key pipeline and lower limit N-1 of key equipment
Figure BDA0002737910820000114
PC+PCP≤Ln+Lm (7)
In the formula (I), the compound is shown in the specification,
Figure BDA0002737910820000115
respectively representing the upper limit and the lower limit of the minimum energy supply capacity of the renewable energy source corresponding to the y load transfer scheme; l isn、LmRepresents the load carried by the pipelines n and m; cRE,zRepresenting the rated capacity of the new energy device z; pC、PCPRepresenting the energy consumption of the compressor and the heat pump. s represents the total number of new energy devices accessed;
with respect to step 102: establishing a regional comprehensive energy system interval practical safety boundary model containing renewable energy sources based on the proposed model;
considering that the operating point is composed of a set of EH pipeline section loads, all the N-1 interval transition safety constraints can be equivalent to M practical safety boundaries. The interval practical safety domain modeling method based on the regional comprehensive energy system containing renewable energy sources has the advantages that the interval practical safety boundary corresponding to the mth pipe section can be expressed as follows:
Figure BDA0002737910820000121
in the formula (I), the compound is shown in the specification,
Figure BDA0002737910820000122
representing the upper and lower limits of the practical safety boundary corresponding to the pipeline m;
Figure BDA0002737910820000123
representing the upper and lower limits of the practical safety boundary corresponding to the pipeline m; lambda [ alpha ]kIndicating that pipeline k corresponds to a pipeline scale factor. L iskRepresents the load carried by the pipeline k;
Figure BDA0002737910820000124
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure BDA0002737910820000125
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of the key pipeline with the distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure BDA0002737910820000126
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with centralized renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure BDA0002737910820000127
respectively representing the upper limit and the lower limit of the minimum energy supply capacity of the renewable energy source corresponding to the y load transfer scheme; pC、PCPRepresenting the energy consumption of the compressor and the heat pump.
With respect to step 103: establishing a regional comprehensive energy system interval maximum energy supply capacity model containing renewable energy sources and a simulation solving method based on the provided safety boundary model;
the sum of the maximum loads which can be supplied by the system is defined as the maximum energy supply capacity of the region, and the region comprehensive energy system region maximum energy supply capacity model containing renewable energy sources is as follows:
Figure BDA0002737910820000128
Figure BDA0002737910820000131
in the formula (I), the compound is shown in the specification,
Figure BDA0002737910820000132
respectively representing the upper and lower limits of the maximum energy supply capacity. L represents a working point vector; h (L) 0 represents a regional comprehensive energy system multipotency flow energy balance equation; λ represents a pipeline scaling coefficient vector;
Figure BDA0002737910820000133
Figure BDA0002737910820000134
upper and lower boundaries representing an upper limit energy supply capacity and a lower limit energy supply capacity, respectively;
the simulation solving method is shown in fig. 1, and the process is as follows:
(1) initializing a primal-dual interior point method based on an interval practical safety domain of a regional comprehensive energy system containing renewable energy, and respectively acquiring initial working points and constraints corresponding to upper and lower limits of the maximum energy supply capacity by using the above formula;
(2) decoupling and calculating the power flow based on the fact that the upper formula is converted into a Lagrange function by a source-dual interior point method, calculating KKT conditions corresponding to the upper limit and the lower limit of the maximum energy supply capacity by the upper formula, and solving the disturbed KKT conditions by a Newton method;
(3) and outputting the upper and lower limits of the maximum energy supply capacity and the corresponding working points.
With respect to step 104: establishing a regional comprehensive energy system interval dimension reduction projection observation method containing renewable energy sources based on the maximum energy supply working point solved by the maximum energy supply capacity model and according to a safe boundary simulation fitting solving method;
and solving an interval practical safety boundary in a 2-dimensional or 3-dimensional space by using a simulation fitting method, and realizing interval dimension reduction projection observation by drawing the boundary. Taking a practical security domain of a three-dimensional interval as an example, the solving is shown in the flowchart 2, and specifically includes:
(1): solving the upper and lower limits of the maximum energy supply capacity and corresponding working points;
(2): based onITSCA critical point array is obtained. Selecting any one group of multi-energy pipeline outlet load combination Lb＝(Lm,Ln,Lo) As a free variable, the remaining M-3 variables are fixed toITSCPipeline section loading when distributed; order (L)m,Ln) Are limited to upper limits by step sizes DeltaL respectively
Figure BDA0002737910820000135
Approaching until reaching the upper limit of the working point, solving the working point to satisfy the constraint of (L)m+Ln,+Lo) Maximum time of LoAnd loadDistribution of operating points Lb＝(Lm,Ln,Lo) Record to the critical point array B1Performing the following steps;
(3): based on
Figure BDA0002737910820000136
A critical point array is obtained. Selecting any one group of multi-energy pipeline outlet load combination Lb＝(Lm,Ln,Lo) As a free variable, the remaining M-3 variables are fixed to
Figure BDA0002737910820000137
Pipeline section loading when distributed; order (L)m,Ln) Are limited to upper limits by step sizes DeltaL respectively
Figure BDA0002737910820000138
Approaching until reaching the upper limit of the working point, solving the working point to satisfy the constraint of (L)m+Ln,+Lo) Maximum time of LoAnd load distribution, and setting the operating point Lb＝(Lm,Ln,Lo) Record to the critical point array B2Performing the following steps;
(4) the method comprises the following steps: fitting critical point array B by using least square method1And B2And the interval dimension reduction projection observation of the interval practical security domain is realized.
With respect to step 105: establishing a practical safe distance model and a simulation solving method of a regional comprehensive energy system interval containing renewable energy sources based on the proposed safe boundary model;
since the interval utility safety margin can be expressed by a hyperplane, the interval utility safety distance based on the interval utility safety margin of the regional integrated energy system including renewable energy is defined as a vertical distance from the operating point to the corresponding safety margin. By using an analytic method, the practical safety distance of the interval corresponding to the mth pipeline segment can be expressed as:
Figure BDA0002737910820000141
in the formula (I), the compound is shown in the specification,
Figure BDA0002737910820000142
representing the upper and lower limits of practical safety distance corresponding to the pipeline m;
Figure BDA0002737910820000143
representing the upper and lower limits of practical safety distance corresponding to the pipeline m; lambda [ alpha ]kRepresenting pipeline k corresponding to the pipeline scale factor;
Figure BDA0002737910820000144
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure BDA0002737910820000145
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of the key pipeline with the distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure BDA0002737910820000146
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with centralized renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure BDA0002737910820000147
respectively representing the upper limit and the lower limit of the minimum energy supply capacity of the renewable energy source corresponding to the y load transfer scheme; pC、PCPRepresenting the energy consumption of the compressor and the heat pump.
The regional comprehensive energy system interval practical safe distance simulation solving model comprises the following steps:
Figure BDA0002737910820000151
Figure BDA0002737910820000152
Figure BDA0002737910820000153
Figure BDA0002737910820000154
Figure BDA0002737910820000155
Figure BDA0002737910820000156
Figure BDA0002737910820000157
Figure BDA0002737910820000158
in the formula, Lq,cur、Lq,optRepresenting the load of the q-th pipeline corresponding to the current working point and the optimal working point; phimRepresenting a pipeline segment set corresponding to the constraint of the lower limit safety inequality on the mth pipeline;
Figure BDA0002737910820000159
IPSB m,u(Lq,opt) Representing the upper and lower bounds of the upper safety margin, with the argument Lq,opt；
Figure BDA00027379108200001510
IPSB m,u(Lq,cur) Representing the upper and lower bounds of the upper safety margin, with the argument Lq,cur；
Figure BDA00027379108200001511
IPSB m,l(Lq,opt) Representing the upper and lower bounds of the lower safety margin, with the argument Lq,opt；
Figure BDA00027379108200001512
IPSB m,l(Lq,cur) Representing the upper and lower bounds of the lower safety margin, with the argument Lq,cur. L represents a working point vector; h (L) 0 represents a regional comprehensive energy system multipotency flow energy balance equation;
the practical safety distance between the intervals is solved by utilizing the primal-dual interior point method, and taking the practical safety distance between the corresponding upper intervals of the No. m pipeline as an example, as shown in FIG. 3, the practical safety distance specifically comprises the following steps:
(1) initializing a primal-dual interior point method based on an interval practical safety domain of a regional comprehensive energy system containing renewable energy, and respectively acquiring initial working points and constraints corresponding to upper and lower limits of maximum energy supply capacity by utilizing the 1 st and 2 nd formulas and the 3 rd and 4 th formulas;
(2) the method comprises the steps of converting an upper formula into a Lagrange function based on a source-dual interior point method, decoupling and calculating a power flow, calculating KKT conditions corresponding to upper and lower limits of practical safety distance by using the upper formula, and solving the disturbed KKT conditions by using a Newton method;
(3) the upper and lower limits of the practical safe distance are output.
With respect to step 106: establishing a regional comprehensive energy system interval full-dimensional observation method containing renewable energy sources based on the proposed safe distance model;
firstly, modeling is carried out on the full-dimensional radius of the interval, and the practical safety distance from a zero-load working point to a practical safety boundary of each interval is defined as:
Figure BDA0002737910820000161
in the formula (I), the compound is shown in the specification,
Figure BDA0002737910820000162
representing the upper and lower limits of the upper full-dimensional radius corresponding to the pipeline m;
Figure BDA0002737910820000163
representing the upper and lower limits of the lower full-dimensional radius for line m.
Figure BDA0002737910820000164
Respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure BDA0002737910820000165
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of the key pipeline with the distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure BDA0002737910820000166
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with centralized renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure BDA0002737910820000167
respectively representing the upper limit and the lower limit of the minimum energy supply capacity of the renewable energy source corresponding to the y load transfer scheme; pC、PCPRepresenting the energy consumption of the compressor and the heat pump.
And solving the full-dimensional radius of the interval by using the method shown in FIG. 3, and drawing the full-dimensional radius boundary in a radar map form to realize full-dimensional observation of the interval.
The feasibility of the optimization method provided by the embodiment of the invention is verified by specific experiments as follows:
an example scene is set by referring to a regional comprehensive energy system which takes electric power, natural gas and heat as energy demand in a certain engineering case, so that the number of transformer substations and regional energy stations and the number of pipelines are reduced for convenience of analysis, and an example topological structure is shown in fig. 4.
The maximum energy supply capacity of the region comprehensive energy system containing renewable energy sources is calculated to be [16.3500,21.9055] MW, and the ITSC allowable load of each key pipeline segment is shown in tables 1 and 2.
TABLE 1 Key line segmentsIs/are as followsITSCAllowable load
Figure BDA0002737910820000171
TABLE 2. of the respective Key line segments
Figure BDA0002737910820000172
Allowable load
Figure BDA0002737910820000173
Selecting a combination of multiple energy pipelines (L)8,L1,L2) The interval utility security domain is observed for free variables in a dimensionality reduction manner, as shown in fig. 5. As can be seen from fig. 5, due to renewable energy uncertainty, unlike the traditional security domain determinacy boundary, the security domain exists in an upper boundary distribution area, a lower boundary distribution area, and an area that is not affected by renewable energy.
Full-dimensional observation is carried out on the safety domain of the example, as shown in fig. 6, it can be seen that the safety domain of the full-dimensional angle is not smooth enough, the area affected by the renewable energy is large in proportion, and a large promotion space is provided.
Those skilled in the art will appreciate that the drawings are only schematic illustrations of preferred embodiments, and the above-described embodiments of the present invention are merely provided for description and do not represent the merits of the embodiments.
The present invention is not limited to the above-described embodiments. The foregoing description of the specific embodiments is intended to describe and illustrate the technical solutions of the present invention, and the above specific embodiments are merely illustrative and not restrictive. Those skilled in the art can make many changes and modifications to the invention without departing from the spirit and scope of the invention as defined in the appended claims.

Claims (8)

1. An interval practical safety domain modeling method is based on a regional comprehensive energy system containing renewable energy sources, and is characterized by comprising the following steps:
considering the uncertainty of renewable energy sources, processing by adopting a region mathematical method, combining engineering practice, and establishing a region practical safety domain model of the regional integrated energy source system containing the renewable energy sources based on the N-1 safety criterion of the regional integrated energy source system and the pipeline segment load on the output side of the energy junction;
establishing a regional comprehensive energy system regional practical safety boundary model containing renewable energy sources based on the regional comprehensive energy system regional practical safety boundary model;
establishing a regional comprehensive energy system interval maximum energy supply capacity model containing renewable energy sources and a simulation solving method based on a safety boundary model;
solving a maximum energy supply working point based on a maximum energy supply capacity model, and establishing a regional comprehensive energy system interval dimension reduction projection observation method containing renewable energy according to a safe boundary simulation solving method;
establishing a practical safe distance model and a simulation solving method of a regional comprehensive energy system interval containing renewable energy sources based on a safe boundary model;
based on a safe distance model, a regional comprehensive energy system interval full-dimensional observation method containing renewable energy is established.
2. The method for modeling an interval utility security domain according to claim 1, wherein the interval utility security domain of regional integrated energy system with renewable energy is defined as:
the method comprises the following steps that (1) safety inequality constraint and energy balance constraint of an N-1 interval are involved, and all working point sets which can meet the safety constraint in a regional comprehensive energy system containing renewable energy sources are collected;
the working point refers to a minimum set of state variables representing the safety of the regional integrated energy system when the regional integrated energy system containing renewable energy is in normal operation; selecting the state variable as the pipeline section load on the output side of the energy hub, wherein the pipeline section load on the output side of the energy coupling equipment and the associated components in the energy hub is included; the energy coupling equipment in the energy hub comprises a gas boiler, cogeneration, an electric boiler, an air conditioning unit and a gas turbine; the associated components include a transformer, a compressor, a heat exchanger, and an energy storage facility.
3. The method according to claim 1, wherein the model of the interval practical safety domain (IPSR-RIES) of the regional integrated energy system containing renewable energy comprises:
Figure FDA0002737910810000011
wherein, L represents a working point vector; h (L) 0 represents a regional comprehensive energy system multipotency flow energy balance equation; λ represents a pipeline scaling coefficient vector;
Figure FDA0002737910810000012
respectively representing an upper limit energy supply capacity matrix and a lower limit energy supply capacity matrix;
Figure FDA0002737910810000013
Figure FDA0002737910810000021
upper and lower boundaries representing an upper limit energy supply capacity and a lower limit energy supply capacity, respectively;
the safety inequality constraint of the N-1 interval comprises the following steps:
(1) safety constraint for transition between output side of key pipeline and upper limit N-1 of key equipment
Figure FDA0002737910810000022
Figure FDA0002737910810000023
Lm+Ln≤Cn
Figure FDA0002737910810000024
In the formula, Ln、LmRepresents the load carried by the pipelines n and m; hjRepresents the load carried by the device j; cn、CjRepresents the rated capacity of pipeline n and equipment j; cpro,zRepresenting the probabilistic capacity of the new energy device z; cRE,zRepresenting the rated capacity of the new energy device z; s represents the total number of new energy devices accessed;
Figure FDA0002737910810000025
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with distributed renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure FDA0002737910810000026
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of the key pipeline with the distributed renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure FDA0002737910810000027
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with centralized renewable energy sources corresponding to the y load transfer scheme under the x fault;
(2) safety constraint for transition between output side of key pipeline and lower limit N-1 of key equipment
Figure FDA0002737910810000028
PC+PCP≤Ln+Lm
In the formula (I), the compound is shown in the specification,
Figure FDA0002737910810000029
respectively representing the upper limit and the lower limit of the minimum energy supply capacity of the renewable energy source corresponding to the y load transfer scheme; l isn、LmRepresents the load carried by the pipelines n and m; cRE,zRepresenting the rated capacity of the new energy device z; s represents the total number of new energy devices accessed; pC、PCPRepresenting the energy consumption of the compressor and the heat pump.
4. The method for modeling an interval utility security domain according to claim 1, wherein the renewable energy-containing regional integrated energy system interval utility security boundary model is:
Figure FDA0002737910810000031
in the formula (I), the compound is shown in the specification,
Figure FDA0002737910810000032
representing the upper and lower limits of the practical safety boundary corresponding to the pipeline m;
Figure FDA0002737910810000033
representing the upper and lower limits of the practical safety boundary corresponding to the pipeline m; lambda [ alpha ]kRepresenting pipeline k corresponding to the pipeline scale factor; l iskRepresents the load carried by the pipeline k;
Figure FDA0002737910810000034
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure FDA0002737910810000035
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of the key pipeline with the distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure FDA0002737910810000036
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with centralized renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure FDA0002737910810000037
respectively representing the upper limit and the lower limit of the minimum energy supply capacity of the renewable energy source corresponding to the y load transfer scheme; pC、PCPRepresenting the energy consumption of the compressor and the heat pump.
5. The interval utility security domain modeling method of claim 1, wherein the renewable energy-containing regional integrated energy system interval maximum energy supply capability model is:
Figure FDA0002737910810000038
Figure FDA0002737910810000039
in the formula (I), the compound is shown in the specification,
Figure FDA00027379108100000310
respectively representing the upper limit and the lower limit of the maximum energy supply capacity; l represents a working point vector; h (L) 0 represents a regional comprehensive energy system multipotency flow energy balance equation; λ represents a pipeline scaling coefficient vector;
Figure FDA0002737910810000041
Figure FDA0002737910810000042
upper and lower boundaries representing an upper limit energy supply capacity and a lower limit energy supply capacity, respectively;
the simulation solving method comprises the following steps:
(1) initializing a primal-dual interior point method based on an interval practical safety domain of a regional comprehensive energy system containing renewable energy, and respectively acquiring initial working points and constraints corresponding to upper and lower limits of the maximum energy supply capacity by using the above formula;
(2) decoupling and calculating the power flow based on the fact that the upper formula is converted into a Lagrange function by a source-dual interior point method, calculating KKT conditions corresponding to the upper limit and the lower limit of the maximum energy supply capacity by the upper formula, and solving the disturbed KKT conditions by a Newton method;
(3) and outputting the upper and lower limits of the maximum energy supply capacity and the corresponding working points.
6. The interval practical safety domain modeling method according to claim 1, wherein the interval dimension reduction projection observation method of the regional integrated energy system containing renewable energy sources is as follows:
solving an interval practical safety boundary in a 2-dimensional or 3-dimensional space by using a simulation fitting method, and realizing interval dimension reduction projection observation by drawing the boundary; taking a practical security domain in a three-dimensional interval as an example, the method specifically comprises the following steps:
(1): solving the upper and lower limits of the maximum energy supply capacity and corresponding working points;
(2): based onITSCObtaining the critical point array, selecting any one group of multi-energy pipeline outlet load combination Lb＝(Lm,Ln,Lo) As a free variable, the remaining M-3 variables are fixed toITSCPipeline section loading when distributed; order (L)m,Ln) Are limited to upper limits by step sizes DeltaL respectively
Figure FDA0002737910810000043
Approaching until reaching the upper limit of the working point, solving the working point to satisfy the constraint of (L)m+Ln,+Lo) Maximum time of LoAnd load distribution, and setting the operating point Lb＝(Lm,Ln,Lo) Record to the critical point array B1Performing the following steps;
(3): based on
Figure FDA0002737910810000044
Obtaining a critical point array; selecting any group of moreEnergy pipeline outlet load combination Lb＝(Lm,Ln,Lo) As a free variable, the remaining M-3 variables are fixed to
Figure FDA0002737910810000045
Pipeline section loading when distributed; order (L)m,Ln) Are limited to upper limits by step sizes DeltaL respectively
Figure FDA0002737910810000046
Approaching until reaching the upper limit of the working point, solving the working point to satisfy the constraint of (L)m+Ln,+Lo) Maximum time of LoAnd load distribution, and setting the operating point Lb＝(Lm,Ln,Lo) Record to the critical point array B2Performing the following steps;
(4): fitting critical point array B by using least square method1And B2And the interval dimension reduction projection observation of the interval practical security domain is realized.
7. The method for modeling an interval utility safety domain according to claim 1, wherein the interval utility safety distance model of the regional integrated energy system with renewable energy comprises:
Figure FDA0002737910810000051
Figure FDA0002737910810000052
in the formula (I), the compound is shown in the specification,
Figure FDA0002737910810000053
representing the upper and lower limits of practical safety distance corresponding to the pipeline m;
Figure FDA0002737910810000054
shows the practical application of the pipeline m under the corresponding conditionUpper and lower limits of safe distance; lambda [ alpha ]kRepresenting pipeline k corresponding to the pipeline scale factor;
Figure FDA0002737910810000055
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure FDA0002737910810000056
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of the key pipeline with the distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure FDA0002737910810000057
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with centralized renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure FDA0002737910810000058
respectively representing the upper limit and the lower limit of the minimum energy supply capacity of the renewable energy source corresponding to the y load transfer scheme; pC、PCPRepresenting the energy consumption of the compressor and the heat pump;
the regional comprehensive energy system interval practical safe distance simulation solving model comprises the following steps:
Figure FDA0002737910810000059
Figure FDA00027379108100000510
Figure FDA00027379108100000511
Figure FDA0002737910810000061
Figure FDA0002737910810000062
Figure FDA0002737910810000063
Figure FDA0002737910810000064
Figure FDA0002737910810000065
in the formula, Lq,cur、Lq,optRepresenting the load of the q-th pipeline corresponding to the current working point and the optimal working point; phimRepresenting a pipeline segment set corresponding to the constraint of the lower limit safety inequality on the mth pipeline;
Figure FDA0002737910810000066
IPSB m,u(Lq,opt) Representing the upper and lower bounds of the upper safety margin, with the argument Lq,opt；
Figure FDA0002737910810000067
IPSB m,u(Lq,cur) Representing the upper and lower bounds of the upper safety margin, with the argument Lq,cur；
Figure FDA0002737910810000068
IPSB m,l(Lq,opt) Representing the upper and lower bounds of the lower safety margin, with the argument Lq,opt；
Figure FDA0002737910810000069
IPSB m,l(Lq,cur) Representing the upper and lower bounds of the lower safety margin, with the argument Lq,cur(ii) a L represents a working point vector; h (L) 0 represents a regional comprehensive energy system multipotency flow energy balance equation;
the practical safety distance between the zones is solved by utilizing a primal-dual interior point method, and the practical safety distance between the zones corresponding to the No. m pipeline is taken as an example, and the practical safety distance between the zones is specifically:
(1) initializing a primal-dual interior point method based on an interval practical safety domain of a regional comprehensive energy system containing renewable energy, and respectively acquiring initial working points and constraints corresponding to upper and lower limits of maximum energy supply capacity by utilizing the 1 st and 2 nd formulas and the 3 rd and 4 th formulas;
(2) the method comprises the steps of converting an upper formula into a Lagrange function based on a source-dual interior point method, decoupling and calculating a power flow, calculating KKT conditions corresponding to upper and lower limits of practical safety distance by using the upper formula, and solving the disturbed KKT conditions by using a Newton method;
(3) the upper and lower limits of the practical safe distance are output.
8. The regional utility security domain modeling method of claim 1, wherein the regional integrated energy system regional full-dimensional observation method with renewable energy comprises the following steps of firstly modeling a regional full-dimensional radius:
Figure FDA0002737910810000071
in the formula (I), the compound is shown in the specification,
Figure FDA0002737910810000072
representing the upper and lower limits of the upper full-dimensional radius corresponding to the pipeline m;
Figure FDA0002737910810000073
representing the upper and lower limits of the lower full-dimensional radius corresponding to the pipeline m;
Figure FDA0002737910810000074
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure FDA0002737910810000075
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of the key pipeline with the distributed renewable energy sources corresponding to the y load transfer scheme under the 1 st fault;
Figure FDA0002737910810000076
respectively representing the upper limit energy supply capacity and the lower limit energy supply capacity of key equipment with centralized renewable energy sources corresponding to the y load transfer scheme under the x fault;
Figure FDA0002737910810000077
respectively representing the upper limit and the lower limit of the minimum energy supply capacity of the renewable energy source corresponding to the y load transfer scheme; pC、PCPRepresenting the energy consumption of the compressor and the heat pump;
and solving the full-dimensional radius of the interval by using an interior point method, and drawing a full-dimensional radius boundary in a radar map form to realize full-dimensional observation of the interval.
CN202011139944.9A 2020-10-22 2020-10-22 Interval practical safety domain modeling method Active CN112434402B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011139944.9A CN112434402B (en) 2020-10-22 2020-10-22 Interval practical safety domain modeling method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011139944.9A CN112434402B (en) 2020-10-22 2020-10-22 Interval practical safety domain modeling method

Publications (2)

Publication Number Publication Date
CN112434402A true CN112434402A (en) 2021-03-02
CN112434402B CN112434402B (en) 2022-08-12

Family

ID=74695889

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011139944.9A Active CN112434402B (en) 2020-10-22 2020-10-22 Interval practical safety domain modeling method

Country Status (1)

Country Link
CN (1) CN112434402B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113327044A (en) * 2021-06-10 2021-08-31 国家电网有限公司 Collaborative operation quality analysis system of energy hub
CN113469424A (en) * 2021-06-22 2021-10-01 天津大学 Multi-target planning method for comprehensive energy system
CN113486509A (en) * 2021-07-01 2021-10-08 天津大学 Multi-objective optimization control method for comprehensive energy system

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120266209A1 (en) * 2012-06-11 2012-10-18 David Jeffrey Gooding Method of Secure Electric Power Grid Operations Using Common Cyber Security Services
US20130238294A1 (en) * 2012-03-07 2013-09-12 Siemens Corporation Global solutions of smart building-grid energy management models
CN105356447A (en) * 2015-10-27 2016-02-24 河海大学 Analysis method for steady-state security region of electric-pneumatic interconnected integrated energy system
CN108494015A (en) * 2018-02-09 2018-09-04 中国科学院电工研究所 The integrated energy system design method of one introduces a collection-lotus-storage coordination and interaction
CN109447360A (en) * 2018-11-06 2019-03-08 国网电力科学研究院(武汉)能效测评有限公司 A kind of energy of cool and thermal power multiple-energy-source energy supplying system purchases scheme optimization method
CN110232217A (en) * 2019-05-16 2019-09-13 天津大学 A kind of comprehensive energy distribution system operation domain modeling method
CN110472830A (en) * 2019-07-12 2019-11-19 中国电力科学研究院有限公司 A kind of Transmission Expansion Planning in Electric method and system considering new energy access
US20200057417A1 (en) * 2018-08-14 2020-02-20 Shanghai Jiaotong University Blockchain-based optimization method for complex scenarios in energy system
WO2020062633A1 (en) * 2018-09-30 2020-04-02 东北电力大学 N-1 static safety analysis method for integrated energy system
CN111245012A (en) * 2020-02-14 2020-06-05 重庆大学 Link line power security domain characterization method considering new energy uncertainty
CN111416340A (en) * 2020-03-10 2020-07-14 天津大学 Regional comprehensive energy system optimization control method based on practical security domain
CN111415271A (en) * 2020-03-17 2020-07-14 清华大学 Comprehensive energy system planning method based on security domain identification

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130238294A1 (en) * 2012-03-07 2013-09-12 Siemens Corporation Global solutions of smart building-grid energy management models
US20120266209A1 (en) * 2012-06-11 2012-10-18 David Jeffrey Gooding Method of Secure Electric Power Grid Operations Using Common Cyber Security Services
CN105356447A (en) * 2015-10-27 2016-02-24 河海大学 Analysis method for steady-state security region of electric-pneumatic interconnected integrated energy system
CN108494015A (en) * 2018-02-09 2018-09-04 中国科学院电工研究所 The integrated energy system design method of one introduces a collection-lotus-storage coordination and interaction
US20200057417A1 (en) * 2018-08-14 2020-02-20 Shanghai Jiaotong University Blockchain-based optimization method for complex scenarios in energy system
US20210216055A1 (en) * 2018-09-30 2021-07-15 Northeast Electric Power University N-1 static security analysis method for integrated energy system
WO2020062633A1 (en) * 2018-09-30 2020-04-02 东北电力大学 N-1 static safety analysis method for integrated energy system
CN109447360A (en) * 2018-11-06 2019-03-08 国网电力科学研究院(武汉)能效测评有限公司 A kind of energy of cool and thermal power multiple-energy-source energy supplying system purchases scheme optimization method
CN110232217A (en) * 2019-05-16 2019-09-13 天津大学 A kind of comprehensive energy distribution system operation domain modeling method
CN110472830A (en) * 2019-07-12 2019-11-19 中国电力科学研究院有限公司 A kind of Transmission Expansion Planning in Electric method and system considering new energy access
CN111245012A (en) * 2020-02-14 2020-06-05 重庆大学 Link line power security domain characterization method considering new energy uncertainty
CN111416340A (en) * 2020-03-10 2020-07-14 天津大学 Regional comprehensive energy system optimization control method based on practical security domain
CN111415271A (en) * 2020-03-17 2020-07-14 清华大学 Comprehensive energy system planning method based on security domain identification

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LIU JUNYONG等: ""Clean energy consumption of power systems towards smart agriculture: roadmap, bottlenecks and technologies"", 《CSEE JOURNAL OF POWER AND ENERGY SYSTEM》 *
刘柳 等: ""面向区域综合能源系统的安全域模型"", 《电力自动化设备》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113327044A (en) * 2021-06-10 2021-08-31 国家电网有限公司 Collaborative operation quality analysis system of energy hub
CN113327044B (en) * 2021-06-10 2023-11-21 国家电网有限公司 Collaborative operation quality analysis system of energy hub
CN113469424A (en) * 2021-06-22 2021-10-01 天津大学 Multi-target planning method for comprehensive energy system
CN113486509A (en) * 2021-07-01 2021-10-08 天津大学 Multi-objective optimization control method for comprehensive energy system
CN113486509B (en) * 2021-07-01 2023-02-28 天津大学 Multi-objective optimization control method for comprehensive energy system

Also Published As

Publication number Publication date
CN112434402B (en) 2022-08-12

Similar Documents

Publication Publication Date Title
Xu et al. Smart energy systems: A critical review on design and operation optimization
CN112434402B (en) Interval practical safety domain modeling method
Zhou et al. An engineering approach to the optimal design of distributed energy systems in China
Bracco et al. A mathematical model for the optimal operation of the University of Genoa Smart Polygeneration Microgrid: Evaluation of technical, economic and environmental performance indicators
Rismanchi District energy network (DEN), current global status and future development
Yan et al. Flexibility improvement and stochastic multi-scenario hybrid optimization for an integrated energy system with high-proportion renewable energy
CN106447529A (en) Distributed energy system modeling and running optimization method considering hot water pipe network
CN109919480B (en) Three-layer target energy Internet planning method and equipment
Chen et al. Economic and environmental operation of power systems including combined cooling, heating, power and energy storage resources using developed multi-objective grey wolf algorithm
Rogers et al. The 20% house–An integrated assessment of options for reducing net carbon emissions from existing UK houses
Wang et al. Improving wind power integration by a novel short-term dispatch model based on free heat storage and exhaust heat recycling
Chen et al. Operation optimization of combined cooling, heating, and power superstructure system for satisfying demand fluctuation
CN109816201A (en) A kind of characterizing method of polyisocyanate mass-energy source homogeneity coupling
Jayasuriya et al. Energy hub modelling for multi-scale and multi-energy supply systems
CN110400090B (en) Multi-target random optimization-based intelligent park multi-energy micro-grid configuration method
Qin et al. Increasing flexibility of combined heat and power systems through optimal dispatch with variable mass flow
Gabbar Smart energy grid infrastructures and interconnected micro energy grids
Tan et al. Exploiting integrated flexibility from a local smart energy hub
CN105976063A (en) Multi-target area multi-microgrid structure and capacity optimization planning method
Li et al. An optimal planning method for CCHP systems based on operation simulation
Wang et al. An integrated energy distribution system planning method with multiple energy storage systems
CN111798044B (en) RIES operation planning simulation method based on improved minimum cross entropy
CN113469424B (en) Multi-target planning method for comprehensive energy system
Li et al. Research on Integrated Smart Energy Evaluation Index System Based on Optimal Algorithm
Pike et al. District heating: a practical solution for reducing fossil fuel dependency in Quebec’s remote communities

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant