CN112528445A - Electricity-gas interconnection comprehensive energy system collaborative planning method and system considering reliability constraint - Google Patents

Electricity-gas interconnection comprehensive energy system collaborative planning method and system considering reliability constraint Download PDF

Info

Publication number
CN112528445A
CN112528445A CN202011464888.6A CN202011464888A CN112528445A CN 112528445 A CN112528445 A CN 112528445A CN 202011464888 A CN202011464888 A CN 202011464888A CN 112528445 A CN112528445 A CN 112528445A
Authority
CN
China
Prior art keywords
planning
comprehensive energy
gas
natural gas
energy system
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
CN202011464888.6A
Other languages
Chinese (zh)
Other versions
CN112528445B (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.)
NARI Group Corp
Nari Technology Co Ltd
State Grid Electric Power Research Institute
Original Assignee
NARI Group Corp
Nari Technology Co Ltd
State Grid Electric Power Research Institute
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 NARI Group Corp, Nari Technology Co Ltd, State Grid Electric Power Research Institute filed Critical NARI Group Corp
Priority to CN202011464888.6A priority Critical patent/CN112528445B/en
Publication of CN112528445A publication Critical patent/CN112528445A/en
Application granted granted Critical
Publication of CN112528445B publication Critical patent/CN112528445B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0637Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/06Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/14Pipes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

The invention provides an electric-gas interconnection comprehensive energy system collaborative planning method considering reliability constraint, which comprises the steps of inputting parameter values required by electric-gas interconnection comprehensive energy system collaborative planning, constructing a collaborative planning main problem model, constructing a power system operation feasibility sub-problem model of a verification collaborative planning scheme, constructing a natural gas system operation feasibility sub-problem model of the verification collaborative planning scheme, constructing an electric-gas interconnection comprehensive energy system operation reliability sub-problem model of the verification collaborative planning scheme, outputting an electric-gas interconnection comprehensive energy system collaborative planning scheme considering reliability constraint and the like, the problems that the model structure is incomplete, the solving method is unreasonable and the like in the traditional collaborative planning process are solved, and the basic guarantee effect of collaborative planning on safe, reliable and economic operation of the electricity-gas interconnection comprehensive energy system is improved. The invention also provides an electric-gas interconnection comprehensive energy system collaborative planning system considering the reliability constraint.

Description

Electricity-gas interconnection comprehensive energy system collaborative planning method and system considering reliability constraint
Technical Field
The invention relates to a comprehensive energy system planning technology, in particular to an electric-gas interconnection comprehensive energy system collaborative planning method considering reliability constraint.
Background
At present, the contradiction between energy demand increase and ecological environment protection in the world is increasingly prominent, natural gas becomes one of main power generation resources with the advantages of economy and environmental protection, and a power system and a natural gas system gradually show a deep coupling trend. Therefore, the operational reliability of the power system is also affected by the risks of natural gas in the supply and transmission processes, and new challenges are provided for the short-term operation and the long-term planning of the electricity-gas interconnected comprehensive energy system. When the system reliability is considered in the conventional collaborative planning method for the electricity-gas interconnection comprehensive energy system, a penalty function reflecting the shortage of system energy supply is generally added in a planning objective function to reflect the influence of the energy supply reliability, so that a compromise scheme comprehensively considering the system economy and the energy supply reliability is obtained.
In the prior art, for example, chinese patent with publication number CN111416349A discloses a collaborative planning method for an electrical interconnection comprehensive energy system. The method takes newly added equipment constraint, power distribution system constraint (externally connected power grid output constraint, distributed gas unit output constraint, distribution network radial constraint, rotating reserve capacity constraint, distribution network constraint and reliability constraint) and natural gas system constraint (urban gate station output constraint, natural gas network constraint and natural gas abundance constraint) as constraint conditions, takes the minimum total cost (sum of investment cost and operation cost) of the comprehensive energy system as an optimization target, constructs a collaborative planning model of the comprehensive energy system, and can determine a collaborative planning scheme with lower cost and stable operation.
However, the prior art has disadvantages. In the aspect of collaborative planning model construction, the consideration of reliability is embodied in that an economic cost penalty function corresponding to load loss is added into a target function, then reliability analysis is carried out on a proposed scheme and adjustment is carried out as required, the probability that the initially selected scheme obtained by the collaborative planning model constructed by the method does not meet the system reliability requirement is extremely high, and scientific and reasonable guide basis is lacked in scheme adjustment; in the prior art, in the aspect of collaborative planning model solution, continuous variable and integer variable decoupling is carried out on a model, a large-scale problem is decomposed, and sub-module alternate iterative solution is carried out.
Therefore, a new technical solution is needed to solve the above problems.
Disclosure of Invention
The purpose of the invention is as follows: the invention provides an electric-gas interconnection comprehensive energy system collaborative planning method considering reliability constraint aiming at the defects of collaborative planning model construction and solution in the prior art, and solves the problems of low reliability of a primary selection scheme, lack of basis for scheme adjustment and the like caused by unreasonable model construction and the problems of inaccurate solution result, unreasonable scheme and the like caused by simple emphasis on mathematical simplification skills due to neglecting the modularization characteristic of a planning problem in the model solving process.
The technical scheme is as follows: in order to achieve the purpose, the invention discloses an electric-gas interconnection comprehensive energy system collaborative planning method considering reliability constraint, which is characterized by comprising the following steps of:
(1) constructing a collaborative planning main problem model taking the lowest total cost of the electricity-gas interconnected comprehensive energy system as a planning target, and obtaining a primary planning scheme according to the collaborative planning main problem model;
(2) constructing a sub-problem model for checking the operation feasibility of the collaborative planning scheme power system to judge whether the primarily selected planning scheme obtained in the step (1) meets the operation feasibility requirement of the power system;
if the requirements are met, the main problem model constructed in the step (1) is considered to be a collaborative planning model considering the operation feasibility of the power system, the obtained primary selection planning scheme is an optimized planning scheme meeting the operation feasibility requirements of the power system, and the next step is carried out;
if the requirement is not met, constructing a power system operation feasibility correction model, circularly iterating the step (1) to the step (2), iteratively adding the power system operation feasibility correction model in a power system operation feasibility test criterion in the circulating process to form a power system operation feasibility test correction criterion, adding the power system operation feasibility test correction criterion into a collaborative planning main problem model to serve as a newly added constraint condition until the obtained optimized planning scheme meets the power system operation feasibility requirement, when the circulation is terminated, forming the main problem model which is the collaborative planning model considering the power system operation feasibility, obtaining the planning scheme which is the optimized planning scheme meeting the power system operation feasibility requirement, and entering the next step;
(3) constructing a sub-problem model for testing the operation feasibility of the natural gas system of the collaborative planning scheme, and judging whether the obtained optimized planning scheme meeting the operation feasibility requirement of the power system meets the operation feasibility requirement of the natural gas system according to the operation feasibility test criterion of the natural gas system;
if the requirement is met, the collaborative planning model considering the operation feasibility of the power system is considered to be the collaborative planning model considering the operation feasibility of the electricity-gas interconnected comprehensive energy system, the obtained optimized planning scheme meeting the operation feasibility requirement of the power system is the optimized planning scheme meeting the operation feasibility requirement of the electricity-gas interconnected comprehensive energy system, and the next step is carried out;
if the requirement is not met, constructing a natural gas system operation feasibility correction model, circularly iterating the step (1) to the step (3), iteratively adding the natural gas system operation feasibility correction model in a natural gas system operation feasibility test criterion in the circulating process to form a natural gas system operation feasibility test correction criterion, adding the natural gas system operation feasibility test correction criterion into a collaborative planning main problem model to serve as a newly added constraint condition until the obtained optimization planning scheme meets the operation feasibility requirement of the electricity-gas interconnection comprehensive energy system, when the circulation is terminated, forming the main problem model which is the collaborative planning model considering the operation feasibility of the electricity-gas interconnection comprehensive energy system, and turning to the next step, wherein the obtained planning scheme is the optimization planning scheme meeting the operation feasibility requirement of the electricity-gas interconnection comprehensive energy system;
(4) constructing a sub-problem model for testing the operation reliability of the electric-gas interconnection comprehensive energy system of the collaborative planning scheme, and judging whether the obtained optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnection comprehensive energy system meets the operation reliability requirement of the electric-gas interconnection comprehensive energy system or not according to an operation reliability test criterion of the electric-gas interconnection comprehensive energy system;
if the requirement is met, the collaborative planning model considering the operation feasibility of the electric-gas interconnection comprehensive energy system is considered to be the collaborative planning model considering the operation reliability of the electric-gas interconnection comprehensive energy system, and the obtained optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnection comprehensive energy system is the optimized planning scheme meeting the operation reliability requirement of the electric-gas interconnection comprehensive energy system;
if the requirement is not met, constructing an electric-gas interconnection comprehensive energy system operation reliability correction model, circularly iterating the step (1) to the step (4), iteratively adding an electric-gas interconnection comprehensive energy system operation reliability correction model in the electric-gas interconnection comprehensive energy system operation reliability test criterion to form an electric-gas interconnection comprehensive energy system operation reliability test correction criterion, and is added into the collaborative planning main problem model to be used as a newly added constraint condition until the obtained optimized planning scheme meets the requirement of the operation reliability of the electricity-gas interconnection comprehensive energy system, when the circulation is terminated, the formed main problem model is a collaborative planning model considering the operation reliability of the electric-gas interconnection comprehensive energy system, and the obtained planning scheme is an optimized planning scheme meeting the operation reliability requirement of the electric-gas interconnection comprehensive energy system.
Corresponding to the collaborative planning method, the invention also provides an electric-gas interconnection comprehensive energy system collaborative planning system considering reliability constraint, which comprises the following steps:
the system comprises a first module, a second module and a third module, wherein the first module is internally provided with a collaborative planning main problem model with the lowest total cost of an electric-gas interconnected comprehensive energy system as a planning target, and obtains a primary planning scheme according to the collaborative planning main problem model;
the second module is internally provided with a sub-problem model for checking the operation feasibility of the power system and a power system operation feasibility correction model;
judging whether the primarily selected planning scheme obtained in the first module meets the requirement of the operation feasibility of the power system or not through a sub-problem model for checking the operation feasibility of the power system of the collaborative planning scheme;
if the requirements are met, the main problem model constructed by the first module is considered to be a collaborative planning model considering the operation feasibility of the power system, and the obtained primary selection planning scheme is an optimized planning scheme meeting the operation feasibility requirements of the power system;
if the requirement is not met, constructing a power system operation feasibility correction model, performing loop iteration, in the loop process, iteratively adding the power system operation feasibility correction model in the power system operation feasibility test criterion to form a power system operation feasibility test correction criterion, adding the power system operation feasibility test correction criterion into the collaborative planning main problem model to serve as a newly added constraint condition until the obtained optimized planning scheme meets the power system operation feasibility requirement, when the loop is terminated, the formed main problem model is the collaborative planning model considering the power system operation feasibility, and the obtained planning scheme is the optimized planning scheme meeting the power system operation feasibility requirement;
the third module is internally provided with a sub-problem model for testing the operation feasibility of the natural gas system and a natural gas system operation feasibility correction model;
judging whether the obtained optimized planning scheme meeting the operation feasibility requirement of the power system meets the operation feasibility requirement of the natural gas system or not according to the operation feasibility test criterion of the natural gas system by a sub-problem model for testing the operation feasibility of the natural gas system of the collaborative planning scheme;
if the requirement is met, the collaborative planning model considering the operation feasibility of the power system is considered to be the collaborative planning model considering the operation feasibility of the electricity-gas interconnected comprehensive energy system, and the obtained optimized planning scheme meeting the operation feasibility requirement of the power system is the optimized planning scheme meeting the operation feasibility requirement of the electricity-gas interconnected comprehensive energy system;
if the requirement is not met, iteratively adding a natural gas system operation feasibility correction model in a natural gas system operation feasibility test criterion through a natural gas system operation feasibility correction model and a loop iteration, wherein the natural gas system operation feasibility test correction criterion is formed and added to a collaborative planning main problem model to serve as a newly added constraint condition until the obtained optimization planning scheme meets the operation feasibility requirement of the electricity-gas interconnected comprehensive energy system, when the loop is terminated, the formed main problem model is the collaborative planning model considering the operation feasibility of the electricity-gas interconnected comprehensive energy system, and the obtained planning scheme is the optimization planning scheme meeting the operation feasibility requirement of the electricity-gas interconnected comprehensive energy system and is transferred to the next step;
the fourth module is internally provided with a sub-problem model for testing the operation reliability of the electric-gas interconnection comprehensive energy system and an electric-gas interconnection comprehensive energy system operation reliability correction model of the collaborative planning scheme;
judging whether the obtained optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnected comprehensive energy system meets the operation reliability requirement of the electric-gas interconnected comprehensive energy system or not according to a sub-problem model for testing the operation reliability of the electric-gas interconnected comprehensive energy system of the collaborative planning scheme and an operation reliability test criterion of the electric-gas interconnected comprehensive energy system;
if the requirement is met, the collaborative planning model considering the operation feasibility of the electric-gas interconnection comprehensive energy system is considered to be the collaborative planning model considering the operation reliability of the electric-gas interconnection comprehensive energy system, and the obtained optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnection comprehensive energy system is the optimized planning scheme meeting the operation reliability requirement of the electric-gas interconnection comprehensive energy system;
if the requirement is not met, the model is corrected through the operation reliability of the electric-gas interconnection comprehensive energy system, and in the circulation iteration and circulation process, iteratively adding an electric-gas interconnection comprehensive energy system operation reliability correction model in the electric-gas interconnection comprehensive energy system operation reliability test criterion to form an electric-gas interconnection comprehensive energy system operation reliability test correction criterion, and is added into the collaborative planning main problem model to be used as a newly added constraint condition until the obtained optimized planning scheme meets the requirement of the operation reliability of the electricity-gas interconnection comprehensive energy system, when the circulation is terminated, the formed main problem model is a collaborative planning model considering the operation reliability of the electric-gas interconnection comprehensive energy system, and the obtained planning scheme is an optimized planning scheme meeting the operation reliability requirement of the electric-gas interconnection comprehensive energy system.
The invention has the beneficial effects that: the electric-gas interconnection comprehensive energy system collaborative planning method considering the reliability constraint is provided, the electric-gas interconnection comprehensive energy system collaborative planning problem is decomposed into a main problem taking the lowest total cost of the electric-gas interconnection comprehensive energy system as a planning target, a sub-problem for checking the operation feasibility of an electric power system of a collaborative planning scheme, a sub-problem for checking the operation feasibility of a natural gas system of the collaborative planning scheme and a sub-problem for checking the operation reliability of the electric power system of the collaborative planning scheme, a main problem model and sub-correction models corresponding to the sub-problems are established, and iterative solution is repeated until a collaborative planning investment scheme meeting the reliability, feasibility and economic constraint of the electric-gas interconnection comprehensive energy system is generated. The method has the advantages that the medium-long term planning process considering the reliability, feasibility and economy of the electric-gas interconnected comprehensive energy system is established, the problems of incomplete model structure, unreasonable solving method and the like in the traditional collaborative planning process are solved, and the basic guarantee level of the collaborative planning optimization scheme on the safe, reliable and economic operation of the electric-gas interconnected comprehensive energy system is improved.
The invention also provides a technical scheme of the electronic equipment, which comprises the following steps:
one or more processors; and a storage device for storing one or more programs which, when executed by the one or more processors, cause the one or more processors to implement the co-planning method described above.
The invention also provides a technical scheme of a computer readable medium, wherein a computer program is stored on the computer readable medium, and the computer program is executed by a processor to realize the collaborative planning method.
Drawings
FIG. 1 is a flow chart of a collaborative planning method for an electric-gas interconnected integrated energy system in consideration of reliability constraints according to the present invention;
fig. 2 is a schematic diagram of an IEEE 118 node system according to embodiment 1 of the present invention;
fig. 3 is a schematic view of a 14-node natural gas system in embodiment 1 of the present invention.
Detailed Description
The present invention will be further described with reference to the accompanying drawings.
Example one
Referring to fig. 1, the present embodiment provides a collaborative planning method for an electrical-gas interconnection energy system considering reliability constraints, which includes the following steps.
Step 1: inputting parameter values required by the electric-gas interconnection comprehensive energy system collaborative planning;
step 2: constructing a collaborative planning main problem model taking the lowest total cost of the electricity-gas interconnected comprehensive energy system as a planning target, and obtaining a primary planning scheme according to the parameter values required by collaborative planning input in the step 1;
and step 3: constructing a sub-problem model for checking the operation feasibility of the power system of the collaborative planning scheme, judging whether the primary planning scheme obtained in the step 2 meets the operation feasibility requirement of the power system according to the operation feasibility check criterion of the power system, if so, considering that the main problem model constructed in the step 2 is the collaborative planning model for considering the operation feasibility of the power system, obtaining the primary planning scheme which is the optimized planning scheme meeting the operation feasibility requirement of the power system, turning to the step 4, if not, constructing a power system operation feasibility correction model, circularly iterating the steps 2 to 3, iteratively adding the power system operation feasibility correction model in the operation feasibility check criterion of the power system in the circulating process to form the operation feasibility check correction criterion of the power system, and adding the criterion into the collaborative planning main problem model as a newly added constraint condition, until the obtained optimized planning scheme meets the operation feasibility requirement of the power system, when the cycle is terminated, the formed main problem model is a collaborative planning model considering the operation feasibility of the power system, the obtained planning scheme is the optimized planning scheme meeting the operation feasibility requirement of the power system, and the step 4 is carried out;
and 4, step 4: constructing a sub-problem model for checking the operation feasibility of the natural gas system of the collaborative planning scheme, judging whether the obtained optimized planning scheme meeting the operation feasibility requirement of the power system meets the operation feasibility requirement of the natural gas system according to an operation feasibility check criterion of the natural gas system, if so, considering that the collaborative planning model considering the operation feasibility of the power system is a collaborative planning model considering the operation feasibility of the electric-gas interconnected comprehensive energy system, obtaining the optimized planning scheme meeting the operation feasibility requirement of the power system is an optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnected comprehensive energy system, turning to step 5, if not, constructing a natural gas system operation feasibility correction model, circularly iterating steps 2 to 4, and iteratively adding the natural gas system operation feasibility correction model in the operation feasibility check criterion of the natural gas system in the circulating process, forming a natural gas system operation feasibility test and correction criterion, adding the criterion into the collaborative planning main problem model as a newly added constraint condition until the obtained optimization planning scheme meets the operation feasibility requirement of the electric-gas interconnection comprehensive energy system, when the circulation is terminated, forming the main problem model which is the collaborative planning model considering the operation feasibility of the electric-gas interconnection comprehensive energy system, obtaining the planning scheme which is the optimization planning scheme meeting the operation feasibility requirement of the electric-gas interconnection comprehensive energy system, and turning to the step 5;
and 5: constructing a sub-problem model for checking the operation reliability of the electric-gas interconnection comprehensive energy system of the collaborative planning scheme, judging whether the obtained optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnection comprehensive energy system meets the operation reliability requirement of the electric-gas interconnection comprehensive energy system according to an operation reliability check criterion of the electric-gas interconnection comprehensive energy system, if so, considering that the collaborative planning model considering the operation feasibility of the electric-gas interconnection comprehensive energy system is the collaborative planning model considering the operation reliability of the electric-gas interconnection comprehensive energy system, and the obtained optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnection comprehensive energy system is the optimized planning scheme meeting the operation reliability requirement of the electric-gas interconnection comprehensive energy system, turning to step 6, if not, constructing an electric-gas interconnection comprehensive energy system operation reliability correction model, circularly iterating the step 2 to the step 5, iteratively adding the electric-gas interconnection comprehensive energy system operation reliability correction model in the operation reliability check criteria of the electric-gas interconnection comprehensive energy system in the circulating process to form an electric-gas interconnection comprehensive energy system operation reliability check correction criteria, adding the electric-gas interconnection comprehensive energy system operation reliability check correction criteria into a collaborative planning main problem model to serve as a newly added constraint condition until the obtained optimization planning scheme meets the operation reliability requirement of the electric-gas interconnection comprehensive energy system, forming the main problem model when the circulation is terminated, namely the collaborative planning model considering the operation reliability of the electric-gas interconnection comprehensive energy system, and obtaining the planning scheme, namely the optimization planning scheme meeting the operation reliability requirement of the electric-gas interconnection comprehensive energy system, turning to step 6;
step 6: and outputting a power-gas interconnection comprehensive energy system collaborative planning scheme considering reliability constraint.
According to the electric-gas interconnection comprehensive energy system collaborative planning method considering the reliability constraint, the parameter values required by the electric-gas interconnection comprehensive energy system collaborative planning in the step 1 comprise: whether the existing generator set is represented on a node or not, whether the existing power transmission line is represented on the node or not, whether the existing natural gas pipeline set is represented on the node or not, whether the newly added generator set is represented on the node or not, whether the newly added power transmission line set is represented between the nodes or not, whether the newly added natural gas pipeline set is represented between the nodes or not, a natural gas source set, a natural gas load set, a node association matrix, operation time, a carbon emission limit value, a carbon emission rate, a load reducible limit value, a load reducible requirement, a system rotation standby requirement, a planning period, new equipment installation and debugging duration, line reactance, a discount rate, a residual value coefficient, a generator set investment cost, a power transmission line investment cost, a natural gas pipeline investment cost, a bus-branch association matrix, a bus-generator association matrix, the system comprises a bus-load incidence matrix, a node natural gas source incidence matrix, a gas unit fuel performance coefficient and the like.
Further, the collaborative planning main problem model with the lowest total cost of the electricity-gas interconnected comprehensive energy system as the planning target in the step 2 takes the lowest comprehensive of the newly-increased investment cost and the operation cost of the collaborative planning as the target, and considers the constraints of the new equipment installation and debugging time, the equipment operation state time sequence change, the total installed capacity of the system meeting the peak demand and the standby demand of load prediction, and the like, and the model can be expressed as follows:
Min C
Figure BDA0002833732330000081
Figure BDA0002833732330000082
Figure BDA0002833732330000083
Figure BDA0002833732330000084
Figure BDA0002833732330000085
Figure BDA0002833732330000086
Figure BDA0002833732330000087
Figure BDA0002833732330000088
Figure BDA0002833732330000089
in the formula, C is newly increased investment cost of collaborative planning, and the unit is ten thousand yuan; t is a planning cycle year index; i is a generator set index; l is a transmission line set index; p isNatural gas pipeline set index; b is a load partition set index; h is a running time set index in a planning period; CG is a set for representing whether a newly added generator set exists on a node or not; CL represents whether a newly added transmission line set exists between nodes; CP is the representation node whether there is a new natural gas pipeline set; EG represents whether the existing generator set is collected on the characterization node; d is the discount rate; GIC is the investment cost of the generator set, and the unit is ten thousand yuan; FIC is the investment cost of the transmission line, and the unit is ten thousand yuan; PIC is the investment cost of the natural gas pipeline, and the unit is ten thousand yuan; z is a binary variable representing the newly added state of the generator set, 1 is newly added, and 0 is not newly added; y is a binary variable representing the newly added state of the power transmission line, 1 is newly added, and 0 is not newly added; x is a binary variable representing the newly added state of the natural gas pipeline, 1 is newly added, and 0 is not newly added; DT is the running time, and the unit is h; OC is the operation cost, and the unit is ten thousand yuan/kW.h; p is the rated power of the generator and has the unit of kW; pmaxThe unit is kW which is the maximum output of the generator; SV is the residual value; gamma is a residual coefficient; t is a planning period and the unit is year; t iscomThe installation and debugging time for new equipment is long, and the unit is year; PD is the load demand, and the unit is kW; r is the system rotation standby requirement, and the unit is kW.
Further, the sub-problem model for checking the operation feasibility of the collaborative planning scheme electric power system in the step 3 takes the minimum loss of the system electric power load as a target, and considers the constraints of the generator set capacity constraint, the system emission limit, the branch power flow constraint of the power transmission line and the like, and the model can be expressed as follows:
Figure BDA0002833732330000091
S.t.
K*PLlbht=A*Pibht-B*(PDbht-DLbht) (11)
Figure BDA0002833732330000092
Figure BDA0002833732330000093
Figure BDA0002833732330000094
Figure BDA0002833732330000095
Figure BDA0002833732330000096
Figure BDA0002833732330000097
Figure BDA0002833732330000098
θref=0 (19)
in the formula, t is a planning cycle year index; i is a generator set index; l is a transmission line set index; b is a load partition set index; h is a running time set index in a planning period; m is a bus head end node set index (bus set index); n is a bus terminal node set index; EG represents whether the existing generator set is collected on the characterization node; EL represents whether the existing transmission line set exists on the node; CG is a set for representing whether a newly added generator set exists on a node or not; CL represents whether a newly added transmission line set exists between nodes; DL is reducible load, and the unit is kW; k is a bus-branch incidence matrix; a is a bus-generator incidence matrix; b is a bus-load incidence matrix; PL is branch power flow, and the unit is kW; p is the rated power of the generator and has the unit of kW; PD is the load demand, and the unit is kW; pmaxThe unit is kW which is the maximum output of the generator;
Figure BDA0002833732330000101
a binary variable for representing the newly added state of the generator set in the last iteration is added, wherein 1 is added, and 0 is not added;
Figure BDA0002833732330000102
the binary variable representing the newly added state of the power transmission line in the last iteration is 1, and 0 is not newly added; EM is carbon emission rate, and the unit is kg/(h.kW); DT is the running time, and the unit is h; EML is the limit of carbon emission in kg; theta is a voltage phase angle; x is the line reactance; m is a large number.
The collaborative planning method for the electricity-gas interconnection comprehensive energy system considering the reliability constraint is characterized in that the electric power system operation feasibility test criterion in the step 3 can be expressed as:
Figure BDA0002833732330000103
in the formula, t is a planning cycle year index; b is a load partition set index; h is a running time set index in a planning period; DL is reducible load, and the unit is kW; DLmaxThe unit is kW to limit the load to be reduced.
The collaborative planning method for the electricity-gas interconnection comprehensive energy system considering the reliability constraint is characterized in that the electric power system operation feasibility correction model in the step 3 can be expressed as:
Figure BDA0002833732330000104
in the formula, V is the operation feasibility correction of the power system, and the unit is kW; t is a planning cycle year index; i is a generator set index; b is a load partition set index; h is a running time set index in a planning period; pmaxThe unit is kW which is the maximum output of the generator; z is a binary variable representing the newly added state of the generator set, 1 is newly added, and 0 is not newly added;
Figure BDA0002833732330000105
and (4) adding 1 and not adding 0 to the binary variable for representing the newly added state of the generator set in the last iteration.
Further, the criterion for checking and correcting the operation feasibility of the power system in step 3 may be represented as:
Figure BDA0002833732330000106
in the formula, t is a planning cycle year index; i is a generator set index; b is a load partition set index; h is a running time set index in a planning period; DL is reducible load, and the unit is kW; DLmaxThe unit is kW which is the limit value capable of reducing the load; v is the operation feasibility correction of the power system, and the unit is kW.
According to the electricity-gas interconnection comprehensive energy system collaborative planning method considering the reliability constraint, the subproblem model for checking the operation feasibility of the collaborative planning scheme natural gas system in the step 4 aims at minimizing the imbalance of supply and demand of the natural gas nodes, considers the constraints of zero net injection of the nodes, load capacity constraint of the natural gas nodes, injection capacity constraint of the natural gas nodes, flow constraint of the natural gas pipeline and the like, and can be expressed as
Min Wbht=1T1+1T2 (23)
Figure BDA0002833732330000111
Fibht=ai+biPibht+ci(Pibht)2 (25)
Lgbht=Fibht(Pibht) (26)
Figure BDA0002833732330000112
Figure BDA0002833732330000113
Figure BDA0002833732330000114
Figure BDA0002833732330000115
Wherein W is the supply and demand unbalance of the natural gas node and the unit is m3H; t is a planning cycle year index; i is a generator set index; l is a transmission line set index; sp is natural gas source set index; g is natural gas load set index; p is a natural gas pipeline set index; j is a natural gas pipe network node index; b is a load partition set index; h is a running time set index in a planning period; CP is the representation node whether there is a new natural gas pipeline set; EP is a set for representing whether existing natural gas pipelines exist on the node; NGS is a natural gas source set; NGL is a natural gas load set; g (j) is a node set communicated with the node j; t is a planning period and the unit is year; omega is a relaxation variable vector; e is a node natural gas source incidence matrix; d is a node natural gas load incidence matrix; v is the gas supply quantity of a natural gas source and is in m3H; l is the natural gas load in m3H; f is the natural gas pipeline flow rate and the unit is m3H; f is the natural gas consumption of the gas unit, and the unit is m3H; p is the rated power of the generator and has the unit of kW; a, b and c are fuel performance coefficients of the gas turbine set;
Figure BDA0002833732330000121
and (3) representing a binary variable of the newly added state of the natural gas pipeline in the last iteration, wherein 1 is added, and 0 is not added.
Further, the natural gas system operation feasibility test criterion in step 4 can be expressed as:
Wbht≤εbht (31)
in the formula, t is a planning cycle year index; b is a load partition set index; h is a running time set index in a planning period; w is the supply and demand unbalance of the natural gas node, and the unit is m3H; epsilon is a preset difference of supply and demand unbalance of the natural gas node and is m3/h。
Further, the natural gas system operation feasibility correction model in step 4 may be expressed as:
Figure BDA0002833732330000122
wherein Y is the correction quantity of the natural gas system operation feasibility and is expressed in m3H; t is a planning cycle year index; i is a generator set index; p is a natural gas pipeline set index; b is a load partition set index; h is a running time set index in a planning period; CG is a set for representing whether a newly added generator set exists on a node or not; CP is the representation node whether there is a new natural gas pipeline set; f is the natural gas consumption of the gas unit, and the unit is m3H; p is the rated power of the generator and has the unit of kW; f is the natural gas pipeline flow rate and the unit is m3H; z is a binary variable representing the newly added state of the generator set, 1 is newly added, and 0 is not newly added;
Figure BDA0002833732330000123
a binary variable for representing the newly added state of the generator set in the last iteration is added, wherein 1 is added, and 0 is not added; x is a binary variable representing the newly added state of the natural gas pipeline, 1 is newly added, and 0 is not newly added;
Figure BDA0002833732330000124
and (3) representing a binary variable of the newly added state of the natural gas pipeline in the last iteration, wherein 1 is added, and 0 is not added.
Further, the natural gas system operation feasibility test correction criterion in step 4 may be expressed as:
Wbht+Ybht≤εbht (33)
in the formula, t is a planning cycle year index; i is a generator set index; b is a load partition set index; h is a running time set index in a planning period; y is the correction quantity of the natural gas system operation feasibility, and the unit is m3H; w is the supply and demand unbalance of the natural gas node, and the unit is m3H; y is the correction quantity of the natural gas system operation feasibility, and the unit is m3H; epsilon is a preset difference of supply and demand unbalance of the natural gas node and is m3/h。
Further, the sub-problem model for checking the operation reliability of the electrical-gas interconnection comprehensive energy system of the collaborative planning scheme in the step 5 aims at minimizing the imbalance of the supply and demand of the power grid node, considers the constraints of the capacity of the generator set, the branch flow constraint of the power transmission line and the like, and can be expressed as follows:
Min Sbht=1T*s1+1T*s2 (34)
S.t.
K*PLlbht+s1-s2=A*Pibht-B*PDbht (35)
Figure BDA0002833732330000131
Figure BDA0002833732330000132
Figure BDA0002833732330000133
Figure BDA0002833732330000134
Figure BDA0002833732330000135
Figure BDA0002833732330000136
θref=0 (42)
in the formula, S is the unbalance amount of the power supply and demand of the power grid node, and the unit is kW; t is a planning cycle year index; i is a generator set index; m is a bus head end node set index (bus set index); n is a bus terminal node set index; l is a transmission line set index; b is a load partition set index; h is a running time set index in a planning period; CG is a set for representing whether a newly added generator set exists on a node or not; CL represents whether a newly added transmission line set exists between nodes; EG represents whether the existing generator set is collected on the characterization node; EL represents whether the existing transmission line set exists on the node; z is a binary variable representing the newly added state of the generator set, 1 is newly added, and 0 is not newly added; y is a binary variable representing the newly added state of the power transmission line, 1 is newly added, and 0 is not newly added; t is a planning period and the unit is year; s is a relaxation variable vector; k is a bus-branch incidence matrix; a is a bus-generator incidence matrix; b is a bus-load incidence matrix; PL is branch power flow, and the unit is kW; p is the rated power of the generator and has the unit of kW; PD is the load demand, and the unit is kW; pmaxThe unit is kW which is the maximum output of the generator; theta is a voltage phase angle; x is the line reactance; m is a large number.
Further, the criterion for testing the operational reliability of the electrical-gas interconnection comprehensive energy system in the step 5 can be expressed as:
Figure BDA0002833732330000141
Figure BDA0002833732330000142
in the formula, t is a planning cycle year index; b is a load partition set index; h is a running time set index in a planning period; EENS is the system energy shortage in kW.h; DT is the running time, and the unit is h; and S is the unbalance amount of the power grid node in kW.
Further, the electric-gas interconnection comprehensive energy system operation reliability correction model in step 5 can be expressed as:
Figure BDA0002833732330000143
in the formula, Z is the operation reliability correction quantity of the electric-gas interconnection comprehensive energy system, and the unit is kW; t is a planning cycle year index; i is a generator set index; l is a transmission line set index; b is a load partition set index; h is a running time set index in a planning period; CG is a set for representing whether a newly added generator set exists on a node or not; CL represents whether a newly added transmission line set exists between nodes; PL is branch power flow, and the unit is kW; p is the rated power of the generator and has the unit of kW; z is a binary variable representing the newly added state of the generator set, 1 is newly added, and 0 is not newly added;
Figure BDA0002833732330000144
a binary variable for representing the newly added state of the generator set in the last iteration is added, wherein 1 is added, and 0 is not added; y is a binary variable representing the newly added state of the power transmission line, 1 is newly added, and 0 is not newly added;
Figure BDA0002833732330000145
and (4) representing a binary variable of the newly added state of the power transmission line in the last iteration, wherein 1 is added, and 0 is not added.
Further, the criterion for checking and correcting the operation reliability of the integrated energy system considering the electrical-electrical interconnection in step 5 can be expressed as:
Figure BDA0002833732330000146
in the formula, t is a planning cycle year index; b is a load partition set index; h is a running time set index in a planning period; EENS is the system energy shortage in kW.h; and Z is the operation reliability correction quantity of the electric-gas interconnection comprehensive energy system, and the unit is kW.
Further, the electric-gas interconnection energy system collaborative planning scheme considering the reliability constraint in step 6 includes: the node where the newly increased generator set is located, the rated capacity of the newly increased generator set, the numbers of the head and tail nodes of the newly increased power transmission line, the rated current-carrying capacity of the newly increased power transmission line, the numbers of the head and tail nodes of the newly increased natural gas pipeline, the rated flow of the newly increased natural gas pipeline, the investment cost, the operation cost and the like.
Based on the collaborative planning method provided by the embodiment, the method is applied to a specific case for verification:
an electric-gas interconnection comprehensive energy system formed by coupling an IEEE 118 node power system (shown in figure 2) and a 14 node natural gas system (shown in figure 3) is selected, and the research and planning period is 20 years. The power system comprises 118 buses, 186 branches, 54 power nodes and 91 load nodes, 16 candidate natural gas generator sets and 8 candidate power transmission lines. The natural gas system comprises 14 nodes, 11 pipelines and 3 gas source nodes, wherein the candidate natural gas pipelines are 3, and the system supplies power for the natural gas generator set and the non-power generation gas load. The load in each month in the planning year is supposed to be divided into three typical load values (valley value, waist value and peak value), the time ratio of the three typical load values to each medium load value in each month is different, the initial power load value is 5400MW, the average annual increase rate is 3%, the rotating standby consideration is 5% of the load value, and the average annual increase rate of the non-power generation natural gas load is 1.5%. In the first planned year, the system energy shortage limit is considered to be 160MWh, and the annual increase is 3%. The total power generation amount in the planning period is limited to 4300 MWh. The average electricity price was 0.05rmb/kWh and was considered to remain unchanged for the planning period. The residual coefficient was 5%. The total investment cost and the newly increased equipment quantity are not limited. According to the electric-gas interconnection comprehensive energy system collaborative planning method considering the reliability constraint, a collaborative planning optimization scheme can be obtained through calculation.
The input parameter values required for partial co-planning are shown in tables 1 to 6.
The values of the parameters covered by the partial co-planning scheme are shown in tables 7 to 10.
Table 1 newly added transmission line parameters
Figure BDA0002833732330000151
TABLE 2 newly added natural gas pipeline parameters
Figure BDA0002833732330000152
Figure BDA0002833732330000161
TABLE 3 Current Natural gas pipeline parameters
Figure BDA0002833732330000162
TABLE 4 Natural gas load parameters
Year t Natural gas load L (m)3/h) Year t Natural gas load L (m)3/h)
1 96986 11 111498
2 98430 12 113636
3 100044 13 115186
4 101460 14 117030
5 102677 15 118647
6 103810 16 120321
7 105141 17 121683
8 106840 18 123868
9 107973 19 125910
10 110278 20 128322
TABLE 5 Natural gas Source parameters
Figure BDA0002833732330000171
TABLE 6 carbon emissions parameters
Figure BDA0002833732330000172
TABLE 7 newly-added parameters of generator set and year of delivery
Figure BDA0002833732330000173
Figure BDA0002833732330000181
Table 8 newly added transmission line parameters and year of commissioning
Figure BDA0002833732330000182
TABLE 9 newly added natural gas pipeline parameters and year of delivery
Figure BDA0002833732330000183
TABLE 10 collaborative planning optimization solution costs
Figure BDA0002833732330000184
Example two
Corresponding to the collaborative planning method provided by the first embodiment, the invention further provides an electric-gas interconnection comprehensive energy system collaborative planning system considering reliability constraint, which includes:
the system comprises a first module, a second module and a third module, wherein the first module is internally provided with a collaborative planning main problem model with the lowest total cost of an electric-gas interconnected comprehensive energy system as a planning target, and obtains a primary planning scheme according to the collaborative planning main problem model;
the second module is internally provided with a sub-problem model for checking the operation feasibility of the power system and a power system operation feasibility correction model;
judging whether the primarily selected planning scheme obtained in the first module meets the requirement of the operation feasibility of the power system or not through a sub-problem model for checking the operation feasibility of the power system of the collaborative planning scheme;
if the requirements are met, the main problem model constructed by the first module is considered to be a collaborative planning model considering the operation feasibility of the power system, and the obtained primary selection planning scheme is an optimized planning scheme meeting the operation feasibility requirements of the power system;
if the requirement is not met, constructing a power system operation feasibility correction model, performing loop iteration, in the loop process, iteratively adding the power system operation feasibility correction model in the power system operation feasibility test criterion to form a power system operation feasibility test correction criterion, adding the power system operation feasibility test correction criterion into the collaborative planning main problem model to serve as a newly added constraint condition until the obtained optimized planning scheme meets the power system operation feasibility requirement, when the loop is terminated, the formed main problem model is the collaborative planning model considering the power system operation feasibility, and the obtained planning scheme is the optimized planning scheme meeting the power system operation feasibility requirement;
the third module is internally provided with a sub-problem model for testing the operation feasibility of the natural gas system and a natural gas system operation feasibility correction model;
judging whether the obtained optimized planning scheme meeting the operation feasibility requirement of the power system meets the operation feasibility requirement of the natural gas system or not according to the operation feasibility test criterion of the natural gas system by a sub-problem model for testing the operation feasibility of the natural gas system of the collaborative planning scheme;
if the requirement is met, the collaborative planning model considering the operation feasibility of the power system is considered to be the collaborative planning model considering the operation feasibility of the electricity-gas interconnected comprehensive energy system, and the obtained optimized planning scheme meeting the operation feasibility requirement of the power system is the optimized planning scheme meeting the operation feasibility requirement of the electricity-gas interconnected comprehensive energy system;
if the requirement is not met, iteratively adding a natural gas system operation feasibility correction model in a natural gas system operation feasibility test criterion through a natural gas system operation feasibility correction model and a loop iteration, wherein the natural gas system operation feasibility test correction criterion is formed and added to a collaborative planning main problem model to serve as a newly added constraint condition until the obtained optimization planning scheme meets the operation feasibility requirement of the electricity-gas interconnected comprehensive energy system, when the loop is terminated, the formed main problem model is the collaborative planning model considering the operation feasibility of the electricity-gas interconnected comprehensive energy system, and the obtained planning scheme is the optimization planning scheme meeting the operation feasibility requirement of the electricity-gas interconnected comprehensive energy system and is transferred to the next step;
the fourth module is internally provided with a sub-problem model for testing the operation reliability of the electric-gas interconnection comprehensive energy system and an electric-gas interconnection comprehensive energy system operation reliability correction model of the collaborative planning scheme;
judging whether the obtained optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnected comprehensive energy system meets the operation reliability requirement of the electric-gas interconnected comprehensive energy system or not according to a sub-problem model for testing the operation reliability of the electric-gas interconnected comprehensive energy system of the collaborative planning scheme and an operation reliability test criterion of the electric-gas interconnected comprehensive energy system;
if the requirement is met, the collaborative planning model considering the operation feasibility of the electric-gas interconnection comprehensive energy system is considered to be the collaborative planning model considering the operation reliability of the electric-gas interconnection comprehensive energy system, and the obtained optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnection comprehensive energy system is the optimized planning scheme meeting the operation reliability requirement of the electric-gas interconnection comprehensive energy system;
if the requirement is not met, the model is corrected through the operation reliability of the electric-gas interconnection comprehensive energy system, and in the circulation iteration and circulation process, iteratively adding an electric-gas interconnection comprehensive energy system operation reliability correction model in the electric-gas interconnection comprehensive energy system operation reliability test criterion to form an electric-gas interconnection comprehensive energy system operation reliability test correction criterion, and is added into the collaborative planning main problem model to be used as a newly added constraint condition until the obtained optimized planning scheme meets the requirement of the operation reliability of the electricity-gas interconnection comprehensive energy system, when the circulation is terminated, the formed main problem model is a collaborative planning model considering the operation reliability of the electric-gas interconnection comprehensive energy system, and the obtained planning scheme is an optimized planning scheme meeting the operation reliability requirement of the electric-gas interconnection comprehensive energy system.
The models and test criteria in the second embodiment are the same as those in the first embodiment, and are not described herein again.
The above embodiments may be implemented by software in combination with necessary hardware. Therefore, the technical solution according to the embodiment of the present disclosure may be embodied in the form of a software product, where the software product may be stored in a nonvolatile storage medium (which may be an optical disc, a usb disk, a mobile hard disk, or the like) or on a network server or a cloud server, and includes a plurality of instructions to enable a computing device (which may be a personal computer, a server, or a network device, or the like) to execute the method according to the first embodiment of the present disclosure, which is not described herein again.

Claims (10)

1. An electric-gas interconnection comprehensive energy system collaborative planning method considering reliability constraint is characterized by comprising the following steps:
(1) constructing a collaborative planning main problem model taking the lowest total cost of the electricity-gas interconnected comprehensive energy system as a planning target, and obtaining a primary planning scheme according to the collaborative planning main problem model;
(2) constructing a sub-problem model for checking the operation feasibility of the collaborative planning scheme power system to judge whether the primarily selected planning scheme obtained in the step (1) meets the operation feasibility requirement of the power system;
if the requirements are met, the main problem model constructed in the step (1) is considered to be a collaborative planning model considering the operation feasibility of the power system, the obtained primary selection planning scheme is an optimized planning scheme meeting the operation feasibility requirements of the power system, and the next step is carried out;
if the requirement is not met, constructing a power system operation feasibility correction model, circularly iterating the step (1) to the step (2), iteratively adding the power system operation feasibility correction model in a power system operation feasibility test criterion in the circulating process to form a power system operation feasibility test correction criterion, adding the power system operation feasibility test correction criterion into a collaborative planning main problem model to serve as a newly added constraint condition until the obtained optimized planning scheme meets the power system operation feasibility requirement, when the circulation is terminated, forming the main problem model which is the collaborative planning model considering the power system operation feasibility, obtaining the planning scheme which is the optimized planning scheme meeting the power system operation feasibility requirement, and entering the next step;
(3) constructing a sub-problem model for testing the operation feasibility of the natural gas system of the collaborative planning scheme, and judging whether the obtained optimized planning scheme meeting the operation feasibility requirement of the power system meets the operation feasibility requirement of the natural gas system according to the operation feasibility test criterion of the natural gas system;
if the requirement is met, the collaborative planning model considering the operation feasibility of the power system is considered to be the collaborative planning model considering the operation feasibility of the electricity-gas interconnected comprehensive energy system, the obtained optimized planning scheme meeting the operation feasibility requirement of the power system is the optimized planning scheme meeting the operation feasibility requirement of the electricity-gas interconnected comprehensive energy system, and the next step is carried out;
if the requirement is not met, constructing a natural gas system operation feasibility correction model, circularly iterating the step (1) to the step (3), iteratively adding the natural gas system operation feasibility correction model in a natural gas system operation feasibility test criterion in the circulating process to form a natural gas system operation feasibility test correction criterion, adding the natural gas system operation feasibility test correction criterion into a collaborative planning main problem model to serve as a newly added constraint condition until the obtained optimization planning scheme meets the operation feasibility requirement of the electricity-gas interconnection comprehensive energy system, when the circulation is terminated, forming the main problem model which is the collaborative planning model considering the operation feasibility of the electricity-gas interconnection comprehensive energy system, and turning to the next step, wherein the obtained planning scheme is the optimization planning scheme meeting the operation feasibility requirement of the electricity-gas interconnection comprehensive energy system;
(4) constructing a sub-problem model for testing the operation reliability of the electric-gas interconnection comprehensive energy system of the collaborative planning scheme, and judging whether the obtained optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnection comprehensive energy system meets the operation reliability requirement of the electric-gas interconnection comprehensive energy system or not according to an operation reliability test criterion of the electric-gas interconnection comprehensive energy system;
if the requirement is met, the collaborative planning model considering the operation feasibility of the electric-gas interconnection comprehensive energy system is considered to be the collaborative planning model considering the operation reliability of the electric-gas interconnection comprehensive energy system, and the obtained optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnection comprehensive energy system is the optimized planning scheme meeting the operation reliability requirement of the electric-gas interconnection comprehensive energy system;
if the requirement is not met, constructing an electric-gas interconnection comprehensive energy system operation reliability correction model, circularly iterating the step (1) to the step (4), iteratively adding an electric-gas interconnection comprehensive energy system operation reliability correction model in the electric-gas interconnection comprehensive energy system operation reliability test criterion to form an electric-gas interconnection comprehensive energy system operation reliability test correction criterion, and is added into the collaborative planning main problem model to be used as a newly added constraint condition until the obtained optimized planning scheme meets the requirement of the operation reliability of the electricity-gas interconnection comprehensive energy system, when the circulation is terminated, the formed main problem model is a collaborative planning model considering the operation reliability of the electric-gas interconnection comprehensive energy system, and the obtained planning scheme is an optimized planning scheme meeting the operation reliability requirement of the electric-gas interconnection comprehensive energy system.
2. The collaborative planning method according to claim 1, wherein: the collaborative planning main problem model taking the lowest total cost of the electricity-gas interconnection comprehensive energy system as a planning target in the step (1) is expressed as follows:
Min C
Figure FDA0002833732320000021
Figure FDA0002833732320000022
Figure FDA0002833732320000023
Figure FDA0002833732320000024
Figure FDA0002833732320000025
Figure FDA0002833732320000031
Figure FDA0002833732320000032
Figure FDA0002833732320000033
Figure FDA0002833732320000034
in the formula, C is newly increased investment cost of collaborative planning; t is a planning cycle year index; i is a generator set index; l is a transmission line set index; p is a natural gas pipeline set index; b is a load partition set index; h is a running time set index in a planning period; CG is a set for representing whether a newly added generator set exists on a node or not; CL represents whether a newly added transmission line set exists between nodes; CP is the representation node whether there is a new natural gas pipeline set; EG represents whether the existing generator set is collected on the characterization node; d is the discount rate; GIC is the generator set investment cost; FIC is the investment cost of the transmission line; PIC is the investment cost of the natural gas pipeline; z is a binary variable representing the newly added state of the generator set, 1 is newly added, and 0 is not newly added; y is a binary variable representing the newly added state of the power transmission line, 1 is newly added, and 0 is not newly added; x is a binary variable representing the newly added state of the natural gas pipeline, 1 is newly added, and 0 is not newly added; DT is the running time; OC is the operating cost; p is the rated power of the generator; pmaxThe maximum output of the generator is obtained; SV is the residual value; gamma is a residual coefficient; t is a planning period; t iscomInstalling a debugging time for the new equipment; PD is the load demand; and R is the system rotation standby requirement.
3. The collaborative planning method according to claim 1 or 2, wherein: the electric power system operation feasibility test criterion in the step 3 in the step (2) is shown as follows:
Figure FDA0002833732320000035
in the formula, t is a planning cycle year index; b is a load partition set index; h is a running time aggregation cable in a planning periodLeading; DL is to reduce the load; DLmaxA limit value for reducing the load;
the power system operation feasibility correction model is expressed as:
Figure FDA0002833732320000036
in the formula, V is the correction quantity of the operation feasibility of the power system; t is a planning cycle year index; i is a generator set index; b is a load partition set index; h is a running time set index in a planning period; pmaxThe maximum output of the generator is obtained; z is a binary variable representing the newly added state of the generator set, 1 is newly added, and 0 is not newly added;
Figure FDA0002833732320000037
a binary variable for representing the newly added state of the generator set in the last iteration is added, wherein 1 is added, and 0 is not added;
the power system operation feasibility test correction criterion is expressed as:
Figure FDA0002833732320000041
in the formula, t is a planning cycle year index; i is a generator set index; b is a load partition set index; h is a running time set index in a planning period; DL is to reduce the load; DLmaxA limit value for reducing the load; and V is the operation feasibility correction quantity of the power system.
4. The collaborative planning method according to claim 1, wherein: the subproblem model for checking the operation feasibility of the collaborative planning scheme natural gas system in the step (3) is expressed as
Min Wbht=1T1+1T2
Figure FDA0002833732320000042
Fibht=ai+biPibht+ci(Pibht)2
Lgbht=Fibht(Pibht)
Figure FDA0002833732320000043
Figure FDA0002833732320000044
Figure FDA0002833732320000045
Figure FDA0002833732320000046
In the formula, W is the unbalance amount of the natural gas node supply and demand; t is a planning cycle year index; i is a generator set index; l is a transmission line set index; sp is natural gas source set index; g is natural gas load set index; p is a natural gas pipeline set index; j is a natural gas pipe network node index; b is a load partition set index; h is a running time set index in a planning period; CP is the representation node whether there is a new natural gas pipeline set; EP is a set for representing whether existing natural gas pipelines exist on the node; NGS is a natural gas source set; NGL is a natural gas load set; g (j) is a node set communicated with the node j; t is a planning period and the unit is year; omega is a relaxation variable vector; e is a node natural gas source incidence matrix; d is a node natural gas load incidence matrix; v is the gas supply quantity of the natural gas source; l is the natural gas load; f is the natural gas pipeline flow; f is the natural gas consumption of the gas turbine set; p is the rated power of the generator; a, b and c are fuel of gas turbine setA coefficient of performance;
Figure FDA0002833732320000047
and (3) representing a binary variable of the newly added state of the natural gas pipeline in the last iteration, wherein 1 is added, and 0 is not added.
5. The collaborative planning method according to claim 1 or 4, wherein: the natural gas system operation feasibility test criterion in the step (3) is represented as:
Wbht≤εbht
in the formula, t is a planning cycle year index; b is a load partition set index; h is a running time set index in a planning period; w is the supply and demand unbalance of the natural gas node; epsilon is a preset difference of the supply and demand unbalance of the natural gas node;
the natural gas system operation feasibility correction model is expressed as:
Figure FDA0002833732320000051
in the formula, Y is the correction quantity of the operation feasibility of the natural gas system; t is a planning cycle year index; i is a generator set index; p is a natural gas pipeline set index; b is a load partition set index; h is a running time set index in a planning period; CG is a set for representing whether a newly added generator set exists on a node or not; CP is the representation node whether there is a new natural gas pipeline set; f is the natural gas consumption of the gas turbine set; p is the rated power of the generator; f is the natural gas pipeline flow; z is a binary variable representing the newly added state of the generator set, 1 is newly added, and 0 is not newly added;
Figure FDA0002833732320000058
a binary variable for representing the newly added state of the generator set in the last iteration is added, wherein 1 is added, and 0 is not added; x is a binary variable representing the newly added state of the natural gas pipeline, 1 is newly added, and 0 is not newly added;
Figure FDA0002833732320000059
the binary variable representing the newly added state of the natural gas pipeline in the last iteration is 1, and 0 is not newly added;
the natural gas system operation feasibility test correction criterion is expressed as:
Wbht+Ybht≤εbht
in the formula, t is a planning cycle year index; i is a generator set index; b is a load partition set index; h is a running time set index in a planning period; y is the operation feasibility correction of the natural gas system; w is the supply and demand unbalance of the natural gas node; y is the operation feasibility correction of the natural gas system; epsilon is a preset difference of the supply and demand unbalance amount of the natural gas node.
6. The collaborative planning method according to claim 1, wherein: the subproblem model for checking the operation reliability of the electric-gas interconnection comprehensive energy system of the collaborative planning scheme in the step (4) is expressed as follows:
Min Sbht=1T*s1+1T*s2
S.t.
K*PLlbht+s1-s2=A*Pibht-B*PDbht
Figure FDA0002833732320000052
Figure FDA0002833732320000053
Figure FDA0002833732320000054
Figure FDA0002833732320000055
Figure FDA0002833732320000056
Figure FDA0002833732320000057
θref=0
in the formula, S is the amount of unbalance of the power supply and demand of the power grid node; t is a planning cycle year index; i is a generator set index; m is a bus head end node set index; n is a bus terminal node set index; l is a transmission line set index; b is a load partition set index; h is a running time set index in a planning period; CG is a set for representing whether a newly added generator set exists on a node or not; CL represents whether a newly added transmission line set exists between nodes; EG represents whether the existing generator set is collected on the characterization node; EL represents whether the existing transmission line set exists on the node; z is a binary variable representing the newly added state of the generator set, 1 is newly added, and 0 is not newly added; y is a binary variable representing the newly added state of the power transmission line, 1 is newly added, and 0 is not newly added; t is a planning period; s is a relaxation variable vector; k is a bus-branch incidence matrix; a is a bus-generator incidence matrix; b is a bus-load incidence matrix; PL is branch power flow; p is the rated power of the generator; PD is the load demand; pmaxThe maximum output of the generator is obtained; theta is a voltage phase angle; x is the line reactance; m is a large number.
7. The collaborative planning method according to claim 1 or 6, wherein: the operation reliability test criterion of the electricity-gas interconnection comprehensive energy system in the step (4) is represented as follows:
Figure FDA0002833732320000061
Figure FDA0002833732320000062
in the formula, t is a planning cycle year index; b is a load partition set index; h is a running time set index in a planning period; EENS is the system energy deficit; DT is the running time, and the unit is h; s is the amount of unbalance of the power supply and demand of the power grid node;
the operation reliability correction model of the electricity-gas interconnection comprehensive energy system is expressed as follows:
Figure FDA0002833732320000063
in the formula, Z is the operation reliability correction quantity of the electric-gas interconnection comprehensive energy system; t is a planning cycle year index; i is a generator set index; l is a transmission line set index; b is a load partition set index; h is a running time set index in a planning period; CG is a set for representing whether a newly added generator set exists on a node or not; CL represents whether a newly added transmission line set exists between nodes; PL is branch power flow; p is the rated power of the generator; z is a binary variable representing the newly added state of the generator set, 1 is newly added, and 0 is not newly added;
Figure FDA0002833732320000064
a binary variable for representing the newly added state of the generator set in the last iteration is added, wherein 1 is added, and 0 is not added; y is a binary variable representing the newly added state of the power transmission line, 1 is newly added, and 0 is not newly added;
Figure FDA0002833732320000065
the binary variable representing the newly added state of the power transmission line in the last iteration is 1, and 0 is not newly added;
the operation reliability test correction criterion considering the electricity-gas interconnection comprehensive energy system is expressed as follows:
Figure FDA0002833732320000066
in the formula, t is a planning cycle year index; b is a load partition set index; h is a running time set index in a planning period; EENS is the system energy deficit; and Z is the operation reliability correction quantity of the electric-gas interconnection comprehensive energy system.
8. An electric-gas interconnected integrated energy system collaborative planning system considering reliability constraint is characterized by comprising the following steps:
the system comprises a first module, a second module and a third module, wherein the first module is internally provided with a collaborative planning main problem model with the lowest total cost of an electric-gas interconnected comprehensive energy system as a planning target, and obtains a primary planning scheme according to the collaborative planning main problem model;
the second module is internally provided with a sub-problem model for checking the operation feasibility of the power system and a power system operation feasibility correction model;
judging whether the primarily selected planning scheme obtained in the first module meets the requirement of the operation feasibility of the power system or not through a sub-problem model for checking the operation feasibility of the power system of the collaborative planning scheme;
if the requirements are met, the main problem model constructed by the first module is considered to be a collaborative planning model considering the operation feasibility of the power system, and the obtained primary selection planning scheme is an optimized planning scheme meeting the operation feasibility requirements of the power system;
if the requirement is not met, constructing a power system operation feasibility correction model, performing loop iteration, in the loop process, iteratively adding the power system operation feasibility correction model in the power system operation feasibility test criterion to form a power system operation feasibility test correction criterion, adding the power system operation feasibility test correction criterion into the collaborative planning main problem model to serve as a newly added constraint condition until the obtained optimized planning scheme meets the power system operation feasibility requirement, when the loop is terminated, the formed main problem model is the collaborative planning model considering the power system operation feasibility, and the obtained planning scheme is the optimized planning scheme meeting the power system operation feasibility requirement;
the third module is internally provided with a sub-problem model for testing the operation feasibility of the natural gas system and a natural gas system operation feasibility correction model;
judging whether the obtained optimized planning scheme meeting the operation feasibility requirement of the power system meets the operation feasibility requirement of the natural gas system or not according to the operation feasibility test criterion of the natural gas system by a sub-problem model for testing the operation feasibility of the natural gas system of the collaborative planning scheme;
if the requirement is met, the collaborative planning model considering the operation feasibility of the power system is considered to be the collaborative planning model considering the operation feasibility of the electricity-gas interconnected comprehensive energy system, and the obtained optimized planning scheme meeting the operation feasibility requirement of the power system is the optimized planning scheme meeting the operation feasibility requirement of the electricity-gas interconnected comprehensive energy system;
if the requirement is not met, iteratively adding a natural gas system operation feasibility correction model in a natural gas system operation feasibility test criterion through a natural gas system operation feasibility correction model and a loop iteration, wherein the natural gas system operation feasibility test correction criterion is formed and added to a collaborative planning main problem model to serve as a newly added constraint condition until the obtained optimization planning scheme meets the operation feasibility requirement of the electricity-gas interconnected comprehensive energy system, when the loop is terminated, the formed main problem model is the collaborative planning model considering the operation feasibility of the electricity-gas interconnected comprehensive energy system, and the obtained planning scheme is the optimization planning scheme meeting the operation feasibility requirement of the electricity-gas interconnected comprehensive energy system and is transferred to the next step;
the fourth module is internally provided with a sub-problem model for testing the operation reliability of the electric-gas interconnection comprehensive energy system and an electric-gas interconnection comprehensive energy system operation reliability correction model of the collaborative planning scheme;
judging whether the obtained optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnected comprehensive energy system meets the operation reliability requirement of the electric-gas interconnected comprehensive energy system or not according to a sub-problem model for testing the operation reliability of the electric-gas interconnected comprehensive energy system of the collaborative planning scheme and an operation reliability test criterion of the electric-gas interconnected comprehensive energy system;
if the requirement is met, the collaborative planning model considering the operation feasibility of the electric-gas interconnection comprehensive energy system is considered to be the collaborative planning model considering the operation reliability of the electric-gas interconnection comprehensive energy system, and the obtained optimized planning scheme meeting the operation feasibility requirement of the electric-gas interconnection comprehensive energy system is the optimized planning scheme meeting the operation reliability requirement of the electric-gas interconnection comprehensive energy system;
if the requirement is not met, the model is corrected through the operation reliability of the electric-gas interconnection comprehensive energy system, and in the circulation iteration and circulation process, iteratively adding an electric-gas interconnection comprehensive energy system operation reliability correction model in the electric-gas interconnection comprehensive energy system operation reliability test criterion to form an electric-gas interconnection comprehensive energy system operation reliability test correction criterion, and is added into the collaborative planning main problem model to be used as a newly added constraint condition until the obtained optimized planning scheme meets the requirement of the operation reliability of the electricity-gas interconnection comprehensive energy system, when the circulation is terminated, the formed main problem model is a collaborative planning model considering the operation reliability of the electric-gas interconnection comprehensive energy system, and the obtained planning scheme is an optimized planning scheme meeting the operation reliability requirement of the electric-gas interconnection comprehensive energy system.
9. An electronic device, comprising:
one or more processors; and storage means for storing one or more programs which, when executed by the one or more processors, cause the one or more processors to implement the co-planning method according to any one of claims 1 to 7.
10. A computer-readable medium, on which a computer program is stored which, when being executed by a processor, carries out the co-planning method according to any one of claims 1 to 7.
CN202011464888.6A 2020-12-14 2020-12-14 Collaborative planning method and system for electricity-gas interconnection comprehensive energy system Active CN112528445B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011464888.6A CN112528445B (en) 2020-12-14 2020-12-14 Collaborative planning method and system for electricity-gas interconnection comprehensive energy system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011464888.6A CN112528445B (en) 2020-12-14 2020-12-14 Collaborative planning method and system for electricity-gas interconnection comprehensive energy system

Publications (2)

Publication Number Publication Date
CN112528445A true CN112528445A (en) 2021-03-19
CN112528445B CN112528445B (en) 2022-11-11

Family

ID=74999366

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011464888.6A Active CN112528445B (en) 2020-12-14 2020-12-14 Collaborative planning method and system for electricity-gas interconnection comprehensive energy system

Country Status (1)

Country Link
CN (1) CN112528445B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111416349A (en) * 2020-04-22 2020-07-14 南方电网科学研究院有限责任公司 Collaborative planning method for electrical interconnection comprehensive energy system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107609312A (en) * 2017-10-17 2018-01-19 上海燃气市北销售有限公司 City gas transmission network gas supply capacity analogy method and simulation system
JP2018159997A (en) * 2017-03-22 2018-10-11 株式会社シンテックホズミ Circuit diagram presentation device

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018159997A (en) * 2017-03-22 2018-10-11 株式会社シンテックホズミ Circuit diagram presentation device
CN107609312A (en) * 2017-10-17 2018-01-19 上海燃气市北销售有限公司 City gas transmission network gas supply capacity analogy method and simulation system

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111416349A (en) * 2020-04-22 2020-07-14 南方电网科学研究院有限责任公司 Collaborative planning method for electrical interconnection comprehensive energy system
CN111416349B (en) * 2020-04-22 2021-12-03 南方电网科学研究院有限责任公司 Collaborative planning method for electrical interconnection comprehensive energy system

Also Published As

Publication number Publication date
CN112528445B (en) 2022-11-11

Similar Documents

Publication Publication Date Title
Li et al. Coordinated scheduling for improving uncertain wind power adsorption in electric vehicles—Wind integrated power systems by multiobjective optimization approach
Liu et al. A Lagrange-multiplier-based reliability assessment for power systems considering topology and injection uncertainties
Yang et al. Integrated grid, coal-fired power generation retirement and GESS planning towards a low-carbon economy
CN111416349A (en) Collaborative planning method for electrical interconnection comprehensive energy system
Lin et al. Optimal ESS allocation in distribution network using accelerated generalised Benders decomposition
Cong et al. Optimal allocation of soft open points in active distribution network with high penetration of renewable energy generations
CN111523247A (en) Discrete solution-based comprehensive energy system multi-energy flow joint calculation method and device
CN113256152A (en) Reliability evaluation method and system for gas-electricity coupling system based on Bayesian network
CN112528445B (en) Collaborative planning method and system for electricity-gas interconnection comprehensive energy system
Li et al. Dynamic environmental/economic scheduling for microgrid using improved MOEA/D-M2M
Ashari et al. Optimal location, size and type of DGs to reduce power losses and voltage deviation considering THD in radial unbalanced distribution systems
CN112564160B (en) Wind power uncertainty-based random configuration method for energy storage system, terminal and storage medium
Jia et al. Robust expansion planning model of integrated energy system with energy hubs integrated
Zhu et al. Distributionally robust optimization model considering deep peak shaving and uncertainty of renewable energy
CN112418732B (en) Planning method and system for comprehensive energy system with multi-energy hub
CN111079982A (en) Planning method, system, medium and electronic device for cable path of wind power plant
Yang et al. Steady-state security region-based chance-constrained optimization for integrated energy systems
CN115189378A (en) Distributed power supply and energy storage grid-connected locating and sizing method and device and electronic equipment
Murali et al. Estimation of locational marginal price in a restructured electricity market with different loss cases using seed genetic algorithm
Li et al. A fast linearized AC power flow-constrained robust unit commitment approach with customized redundant constraint identification method
Nikolaidis et al. Indices to assess the integration of renewable energy resources on transmission systems
Tee et al. Battery energy storage system sizing using PSO algorithm in DIgSILENT PowerFactory
Matcha et al. Regular paper LMP Calculation with Distributed Loss using GA based DCOPF
Hong et al. Multi-objective siting and sizing of distributed generation based on the mixed integer second-order cone programming
Souheyla et al. Golden Jackal Algorithm for Optimal Size and Location of Distributed Generation in Unbalanced Distribution Networks

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