CN112288592B - Gas-thermal coupling system SCUC optimization scheduling method, device and storage medium - Google Patents
Gas-thermal coupling system SCUC optimization scheduling method, device and storage medium Download PDFInfo
- Publication number
- CN112288592B CN112288592B CN202011125210.5A CN202011125210A CN112288592B CN 112288592 B CN112288592 B CN 112288592B CN 202011125210 A CN202011125210 A CN 202011125210A CN 112288592 B CN112288592 B CN 112288592B
- Authority
- CN
- China
- Prior art keywords
- gas
- unit
- node
- power
- chp
- 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.)
- Active
Links
- 230000008878 coupling Effects 0.000 title claims abstract description 47
- 238000010168 coupling process Methods 0.000 title claims abstract description 47
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000003860 storage Methods 0.000 title claims abstract description 9
- 238000005457 optimization Methods 0.000 title claims description 13
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 71
- 239000007789 gas Substances 0.000 claims abstract description 68
- 239000003345 natural gas Substances 0.000 claims abstract description 33
- 230000009194 climbing Effects 0.000 claims abstract description 11
- 239000003245 coal Substances 0.000 claims abstract description 8
- 239000000446 fuel Substances 0.000 claims abstract description 8
- 238000012887 quadratic function Methods 0.000 claims abstract description 7
- 238000010977 unit operation Methods 0.000 claims abstract description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 43
- 230000006870 function Effects 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 230000005611 electricity Effects 0.000 claims description 5
- 239000002918 waste heat Substances 0.000 claims description 5
- 125000004432 carbon atom Chemical group C* 0.000 claims description 4
- 238000010438 heat treatment Methods 0.000 claims description 4
- 238000004134 energy conservation Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 238000010248 power generation Methods 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 claims description 2
- 229910052731 fluorine Inorganic materials 0.000 claims description 2
- 125000001153 fluoro group Chemical group F* 0.000 claims description 2
- 230000011218 segmentation Effects 0.000 claims description 2
- 238000005485 electric heating Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 230000002093 peripheral effect Effects 0.000 description 3
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000009440 infrastructure construction Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 238000009987 spinning Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Power Engineering (AREA)
- Human Resources & Organizations (AREA)
- Mathematical Optimization (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Primary Health Care (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- General Business, Economics & Management (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses an SCUC (substation configuration unit) optimal scheduling method and device for a gas-thermoelectric coupling system and a storage medium. The optimized scheduling method comprises the following steps: 1. reading in original data of the gas-thermoelectric coupling system; 2. defining a system state decision variable; 3. establishing a quadratic function piecewise linearization model; 4. setting unit operation state parameters, unit climbing, start-stop time constraint and thermal power unit output constraint; 5. setting natural gas consumption balance constraint, natural gas unit output constraint and each node flow air pressure constraint; 6. setting flow, heat and temperature balance constraints, and heat load and heat power balance constraints of each node; 7. setting power balance of a coupling link and equipment output constraint; 8. establishing an SCUC optimal scheduling model of the gas-thermal electric coupling system; the invention considers the safe operation constraint condition of the generator set, and considers two generator sets of coal and gas, and the difference of fuel price can show the using time period of the generator set in the result.
Description
Technical Field
The invention relates to the field of optimization scheduling of an integrated energy system, in particular to an SCUC optimization scheduling method and device of a gas-thermal coupling system and a storage medium.
Background
Along with the rapid development of social economy, the urbanization process is deepened continuously, the contradictions between energy consumption and environmental protection, energy infrastructure construction and space limitation and the like are more and more prominent, and meanwhile, the problems of energy supply bottleneck, low energy utilization rate and the like are increased due to independent planning and design and independent operation of different types of energy systems. In the comprehensive energy system, various energy devices are closely coupled to realize interaction and conversion of multiple energy sources, and a proper scheduling operation scheme is formulated through research and calculation, so that the energy utilization efficiency can be greatly improved, the environmental pollution is reduced, the load regulation rate and the equipment utilization rate are improved, and the high-efficiency and gradient utilization of the energy sources is realized. At present, comprehensive energy systems researched at home and abroad comprise an electricity-gas coupling system, an electricity-heat mixed energy network, an electricity-heat-gas coupling system or a combined cooling heating and power system and the like, and the purpose is to achieve the effects of cooperative optimization and advantage complementation.
The main contents of the research of the gas-thermal electric coupling system comprise multi-energy flow calculation, economic optimization scheduling, optimal power flow and the like, and currently, most of the researches are carried out by scholars, and the researches are less for the unit combination optimization problem with safety constraint conditions. Under the constraint conditions of generator set climbing, descending, equipment start-stop and on-off time and the like, the gas-thermal-electric coupling system not only meets the requirements of gas, heat and electricity loads, but also can ensure that the fuel cost is minimum and meets the requirements of economic and safe operation by adjusting the output of the coal-fired unit and the gas-fired unit at different time periods throughout the day.
Disclosure of Invention
In order to solve the defects mentioned in the background technology, the invention aims to provide a combined optimal scheduling method for a gas-thermoelectric coupling system generator set with safety constraint conditions, the invention establishes a constraint condition model in the aspect of unit safety, and comprises the aspects of climbing, starting and stopping, running time, output limitation and the like, and adopts a piecewise linearization method to carry out linearization treatment on nonlinear problems encountered in the calculation process, so that the solution can be accelerated, the probability of obtaining an optimal value is increased, the dynamic process of natural gas flow can be accurately simulated, and the calculation result is closer to an actual value; the output of the coal-fired unit and the gas-fired unit in different time periods throughout the day is adjusted, so that the gas-thermoelectric coupling system can meet the requirements of gas, heat and electricity loads and ensure the minimum fuel cost under the condition of ensuring safe operation.
The purpose of the invention can be realized by the following technical scheme:
a gas-thermoelectric coupling system SCUC optimal scheduling method comprises the following steps:
1) Reading in original data of the gas-thermal-electric coupling system and carrying out initialization processing; the power system adopts an IEEE39 node system, the natural gas system adopts a 20-node network, the thermal system adopts a 6-node network, and the coupling link adopts CHP and electric boiler equipment;
2) Defining system state variables and decision variables, calculating the power flow of the power system by using a direct current model, and setting power balance and power flow out-of-limit safety constraints;
3) Establishing a quadratic function piecewise linearization model;
4) Setting unit operation state parameters, and setting unit climbing and start-stop time constraints and thermal power unit output constraints;
5) Calculating the load flow of a gas network according to a natural gas linear model, calculating a correlation matrix of related node branches, and setting a natural gas consumption balance constraint, a natural gas unit output constraint and a flow and air pressure balance constraint of each node;
6) Calculating the flow of a heat supply network according to a hydraulic model and a thermodynamic model, calculating a related node branch incidence matrix, setting flow, heat and temperature balance constraints according to energy conservation and a thermodynamic theorem, and adding heat load and heat power balance constraints of each node;
7) The power balance and the equipment output constraint of a coupling link are set, the combined heat and power CHP comprises a gas turbine set and a waste heat boiler, power generation and heat supply are realized, the electric boiler is used as a heat source, and the purpose of system optimization scheduling is achieved by adjusting the output of the coupling equipment;
8) And establishing an SCUC optimization scheduling model of the gas-thermoelectric coupling system to calculate a target function value, adopting a piecewise linearization processing method for a nonlinear function, programming through Matlab software, and solving an optimal value by using a Cplex solver.
Further, the system state variables in the step 2) comprise voltage and phase angle in the power system, node pressure of the natural gas system, flow and temperature of the thermodynamic system, CHP efficiency, low heat value and electric boiler efficiency; the system decision variables comprise the output of a generator in the power system, the transformation ratio of a natural gas compressor, the output of a gas source point, the output of a heat source of a thermodynamic system and the temperature; calculating the power flow of the power system by using a direct current model, and setting power balance and power flow out-of-limit safety constraints; the branch and node active power flow balance and constraint are as follows:
in the formula, P ij Is the active power flow of the line i, j; theta.theta. i 、θ j The voltage phase angles of nodes at two ends of the line are respectively; x is the number of ij Is a line reactance;
for the line current to allow a maximum value, ij θ、
is a voltage phase angle limit value; m G 、M EB Coefficient matrixes of nodes where the generator set and the electric boiler are located are respectively set; the electricity consumption of the electric boiler is regarded as the load,
in order to load the node of the electric boiler,
is the node load;
the power is output by the generator set,
respectively the upper and lower limit values of the output of the generator; n is the node number, B ij A susceptance matrix for the power system.
Further, the quadratic function curve piecewise linearization model in the step 3) is as follows:
in the formula, m is the number of the sections divided by the variable value interval; 1 is a segment node number;
variable values segmented for t period 1;
is the total length of the first 1 segment of the variable interval;
for the segmentation value of the 1 st subinterval, l ∈ {1,2, \8230;, m };
1 segment of binary variable;
is the slope of segment 1; w is a ij,t Is the natural gas flow between nodes i, j.
Further, the unit operation state parameters are set in the step 4), and unit climbing and start-stop time constraints and unit output constraints are set as follows:
in the above-mentioned formula, the compound has the following structure,
the generated power of the unit j at the moment t is obtained;
the binary state variable of the unit is represented, the value of 1 represents that the unit j is in the running state at the moment of t, and the value of 0 represents that the unit j is in the shutdown state; j Pthe minimum value of the output of the unit j is shown,
representing the maximum value of the j output of the unit;
representing the maximum value that the unit j can obtain at the moment t;
the upper formula is climbing and start-up restriction, RU j For unit j ramp limit, SU j Starting a limit value for the unit j;
the upper formula is shutdown constraint, SD j Stopping the limiting value for the unit j;
the upper form being downhill restraint, RD j A downhill limit value for the unit j;
the above formula is the minimum start time constraint of the unit, UT j The minimum startup time of the unit j;
the starting time of the unit j before scheduling;
the running state of the unit j before dispatching; g j The time for which the unit j must operate;
the above formula is the minimum shutdown time constraint of the unit, DT j The minimum shutdown time of the unit j;
the shutdown time of the unit j before scheduling; l is a radical of an alcohol j For the time that unit j must be shut down.
Further, the step 5) is as follows according to the fluid mechanics law:
w ij,t +w ji,t ＝0
in the above formula, w ijt For the flow of natural gas pipeline node i to node j,
respectively injecting the gas quantity of the point j, the gas consumption of the gas unit and the gas load quantity of the point j into the gas source point; z (j) and v (j) are respectively a pipeline set taking the node j as a tail node and taking the node j as a head node; c ij Is the constant of the pipe i, j; psi i,t Is the air pressure at node i;ψ、
respectively an upper limit and a lower limit of air pressure; ij w、
respectively an upper limit and a lower limit of the pipeline flow; sgn (psi) i,t ,ψ j,t ) And a sign function representing the natural gas flow direction in the pipeline, wherein the gas flow flows from the node with high pressure to the node with low pressure.
Further, the step 6) is as follows:
Φ i ＝C p m i (T i ms -T i mr )
in the above formula, the first and second carbon atoms are,
respectively the outlet water temperature and the inlet water temperature of a pipeline b in the water supply system;
the outlet water temperature and the inlet water temperature of a pipeline b in the water return system are respectively; t is i ms 、T i mr Respectively the mixed temperature of the water supply system and the return water system at the node i; ms is b 、mr b Water flow in the water supply system pipeline b and the water return system pipeline b respectively;
respectively a pipeline set taking the node i as a head end and a pipeline set taking the node i as a tail end; phi i Being a node iA thermal load power; c p The specific heat capacity of water is 4200J/(kg x DEG C); m is i The amount of water injected into the node i;
the outlet water temperature and the inlet water temperature of the pipeline b are respectively; t is a Is ambient temperature; lambda [ alpha ] b Is the heat transfer coefficient of the pipe b; l is a radical of an alcohol b Is the length of the pipe b.
Further, the step 7) is expressed as:
H CHP,t +H EB,t ＝H source,t
H source,t ＝H load,t +H loss,t
H CHP,t ＝η chp,h *P CHP,t
H EB,t ＝η EB *P EB,t
C gas ＝C ng *Q gas
in the above formula, H CHP,t 、H EB,t The heat powers generated at t moments of the CHP waste heat boiler and the electric boiler are respectively; h source,t Providing the total thermal power for the system at the moment t; h load,t 、H loss,t Respectively the heat load and the heat loss of the system at the time t; p is CHP,t Electrical power generated for the CHP genset; p EB,t The power consumption of the electric boiler; eta EB Is the electric heat conversion coefficient of the electric boiler; eta chp,h CHP electric heating proportionality coefficient; eta chp,e Is the efficiency coefficient of the CHP unit; ζ represents a unit LHV Is natural gas with low heat value; q gas The gas consumption amount of the natural gas is increased; c ng Is the cost per unit of natural gas; c gas The total fuel cost for natural gas.
Further, the objective function in step 8) is:
in the formula, a i 、b i 、c i Coefficients that are a function of the coal cost, respectively;
the active power is the active power generated by the coal-fired unit where the node i is located at the time t; c ng,i,t A unit gas supply cost coefficient of the gas source point i at the moment t; q gas,i,t The air supply quantity of an air source point i at the moment t; gas,i Q、
respectively, the limit values of the gas supply at the gas source point i.
An apparatus, comprising: the system comprises one or more processors and a memory, wherein the memory is used for storing one or more programs, and the one or more programs are executed by the one or more processors to realize the SCUC optimized scheduling method of the gas-thermoelectric coupling system.
A storage medium, comprising: the computer executable instructions, when executed by a computer processor, are for performing a gas-thermoelectric coupling system SCUC optimized scheduling method.
The invention has the beneficial effects that:
1. the invention considers the safe operation constraint conditions of the generator set, including active power output constraint, climbing and downhill slope limitation, unit startup and shutdown constraint, unit operation time limitation and the like; in addition, two power generating sets of coal and gas are considered, and the use time period of the power generating sets can be reflected in the result due to the difference of fuel prices;
2. the invention considers the difficulty brought by the nonlinear function in the calculation process, and greatly improves the operation efficiency and accuracy by the piecewise linearization processing;
3. the invention can calculate the power flow distribution condition of each subsystem, and the output of the coupling equipment and the unit; the established SCUC optimization scheduling model of the gas-thermoelectric coupling system ensures that the system not only meets the requirements of gas, heat and electric loads, but also ensures that the fuel cost is minimum and meets the requirements of economy and safe operation under the condition of ensuring safe operation.
Drawings
The invention is further described below with reference to the accompanying drawings.
FIG. 1 is a flow chart of the steps of the present optimized scheduling method;
FIG. 2 is a gas load demand diagram of the present optimized scheduling method;
FIG. 3 is a graph of the thermal load requirements of the present optimized scheduling method;
FIG. 4 is a graph of the electrical load requirements of the present optimized scheduling method;
FIG. 5 is a block diagram of the present optimized scheduling method;
FIG. 6 is a diagram of an example of the optimized scheduling method;
fig. 7 is a device and equipment diagram of the implementation of the optimized scheduling method.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the description of the present invention, it is to be understood that the terms "opening," "upper," "lower," "thickness," "top," "middle," "length," "inner," "peripheral," and the like are used in an orientation or positional relationship that is merely for convenience in describing and simplifying the description, and do not indicate or imply that the referenced component or element must have a particular orientation, be constructed and operated in a particular orientation, and thus should not be considered as limiting the present invention.
An optimal scheduling method for a gas-thermoelectric coupling system SCUC is shown in FIGS. 1, 5 and 6, and comprises the following steps:
step 1: reading in original data of the gas-thermal-electric coupling system and carrying out initialization processing; the power system adopts an IEEE39 node system, the natural gas system adopts a 20-node network, the thermal system adopts a 6-node network, and the coupling link adopts CHP and electric boiler equipment. The original data comprises node data and branch data of each subsystem, the power system also comprises generator set data, 10 machine sets comprise a coal-fired machine set and a gas-fired machine set, and the gas-fired machine set comprises 1 CHP gas-fired machine set; load data, spinning reserve data, etc. The natural gas system also comprises gas source data, gas supply nodes, load data and the like; the thermodynamic system also includes ambient temperature, thermal load data, and the like.
Step 2: defining system state variables and decision variables; system state variables such as voltage and phase angle in the power system, node pressure of the natural gas system, flow and temperature of the thermodynamic system, CHP efficiency, low heating value, electric boiler efficiency and the like; the system decision variables include generator output in the power system, natural gas compressor transformation ratio, gas source point output, thermal system heat source output, temperature and the like. Calculating the power flow of the power system by using a direct current model, and setting power balance and power flow out-of-limit safety constraints; the branch and node active power flow balance and constraint are as follows:
in the formula, P ij Is the active power flow of the line i, j; theta.theta. i 、θ j The voltage phase angles of nodes at two ends of the line are respectively; x is a radical of a fluorine atom ij Is a line reactance;
for the line current to allow a maximum value, ij θ、
is a voltage phase angle limit value; m is a group of G 、M EB Coefficient matrixes of nodes where the generator set and the electric boiler are located are respectively set; the electricity consumption of the electric boiler is regarded as the load,
is the node load of the electric boiler,
is the node load;
the power is output by the generator set,
respectively the upper and lower limit values of the output of the generator; n is the node number, B ij A power system susceptance matrix.
And step 3: and establishing a quadratic function piecewise linearization model. The quadratic function curve piecewise linearization model is as follows:
in the formula, m is the number of the sections divided by the variable value interval; 1 is a segment node number;
variable values segmented for t period 1;
is the total length of the first 1 segment of the variable interval;
is the segment value of the 1 st subinterval, l is the {1,2, \8230 ∈, m };
1 segment of binary variable;
is the slope of segment 1; w is a ij,t Is the natural gas flow between nodes i, j.
And 4, step 4: the method comprises the following steps of setting unit running state parameters, setting unit climbing and start-stop time constraints, and thermal power unit output constraints, as follows:
in the above formula, the first and second carbon atoms are,
the generated power of the unit j at the time t is obtained;
the binary state variable of the unit is represented, the value of 1 represents that the unit j is in the running state at the moment of t, and the value of 0 represents that the unit j is in the shutdown state; j Pthe minimum value of the output of the unit j is shown,
representing the maximum value of the j output of the unit;
representing the maximum value that the unit j can obtain at the moment t;
the upper formula is climbing and start-up constraint, RU j For unit j ramp limit, SU j Starting a limit value for the unit j;
the above formula is shutdown constraint, SD j Stopping the limiting value for the unit j;
the upper form being downhill restraint, RD j A downhill limit value for the unit j;
the above formula is the minimum start time constraint of the unit, UT j The minimum startup time of the unit j;
the starting time of the unit j before scheduling;
the running state of the unit j before dispatching; g j The time for which the unit j must operate;
the above formula is the minimum shutdown time constraint of the unit, DT j The minimum shutdown time of the unit j;
the shutdown time of the unit j before scheduling; l is j For the time that unit j must be shut down.
And 5: calculating the load flow of a gas network according to a natural gas linear model, calculating a correlation matrix of related node branches, setting a natural gas consumption balance constraint, a natural gas unit output constraint and a flow and air pressure balance constraint of each node, and according to a fluid mechanics law, the method comprises the following steps:
w ij,t +w ji,t ＝0
in the above formula, w ij,t For the flow of natural gas pipeline node i to node j,
respectively injecting the gas quantity of the point j, the gas consumption of the gas unit and the gas load quantity of the point j into the gas source point; z (j) and v (j) are respectively a pipeline set taking the node j as a tail node and taking the node j as a head node; c ij Is the constant of the pipe i, j; psi i,t Is the air pressure at node i;ψ、
respectively an upper limit and a lower limit of air pressure; ij w、
respectively an upper limit and a lower limit of the pipeline flow; sgn (psi) i,t ,ψ j,t ) And a sign function representing the natural gas flow direction in the pipeline, wherein the gas flow flows from the node with high pressure to the node with low pressure.
Because the gas flow and the gas pressure are in a nonlinear relation, linearization is carried out according to the method in the step 3, and therefore the solution can be rapidly and accurately carried out.
And 6: calculating the flow of a heat supply network according to a hydraulic model and a thermodynamic model, calculating a related node branch incidence matrix, setting flow, heat and temperature balance constraints according to energy conservation and a thermodynamic theorem, and adding thermal load and thermal power balance constraints of each node as follows:
in the above formula, the first and second carbon atoms are,
the water temperature of an outlet and the water temperature of an inlet of a pipeline b in the water supply system are respectively;
the outlet water temperature and the inlet water temperature of a pipeline b in the water return system are respectively;
respectively the mixed temperature of the water supply system and the return water system at the node i; ms is b 、mr b Water flow in the water supply system and the backwater system pipeline b respectively;
respectively, a pipeline set taking the node i as a head end and a pipeline set taking the node i as a tail end; phi i Is the thermal load power of node i; c p The specific heat capacity of water is 4200J/(kg x DEG C); m is a unit of i The amount of water injected into the node i;
the outlet water temperature and the inlet water temperature of the pipeline b are respectively; t is a unit of a Is ambient temperature; lambda [ alpha ] b Is the heat transfer coefficient of the pipe b; l is b Is the length of the pipe bAnd (4) degree.
And 7: setting power balance of a coupling link and equipment output constraint; the coupling equipment comprises a CHP and an electric boiler, the electric boiler and a CHP waste heat boiler are used as heat sources to supply heat for a thermodynamic system, and electric power generated by the CHP unit is transmitted to a No. 30 node of the electric power system. According to the original data, the gas power generation is more expensive than the fire coal, if the environmental protection factor is not considered, the priority of the fire coal generator set is higher than that of the gas generator set in order to ensure the minimum total cost of the fuel, and the calculation result shows that the priority is higher than that of the gas generator set.
H CHP,t +H EB,t ＝H source,t
H source,t ＝H load,t +H loss,t
H CHP,t ＝η chp,h *P CHP,t
H EB,t ＝η EB *P EB,t
C gas ＝C ng *Q gas
In the above formula, H CHP,t 、H EB,t The heat powers generated at t moments of the CHP waste heat boiler and the electric boiler are respectively; h source,t The total thermal power provided for the system at the time t; h load,t 、H loss,t Respectively the heat load and the heat loss of the system at the time t; p is CHP,t Electrical power generated for the CHP genset; p EB,t The power consumption of the electric boiler; eta EB Is the electric heat conversion coefficient of the electric boiler; eta chp,h CHP electric heating proportionality coefficient; eta chp,e Is the efficiency coefficient of the CHP unit; zeta LHV Is natural gas with low heat value; q gas The gas consumption amount of the natural gas is increased; c ng Is the cost per unit of natural gas; c gas The total fuel cost for natural gas.
And 8: the method comprises the steps of calculating an objective function value by using a gas-thermoelectric coupling system SCUC optimization scheduling model, adopting a piecewise linearization processing method for a nonlinear function, and rapidly solving an optimal value by using a solver, wherein the fuel cost consumed by the system is 6.11X 106 dollars, the gas cost accounts for 23.3%, the coal cost accounts for 76.7%, and the objective function is as follows:
in the formula, a i 、b i 、c i Coefficients that are a function of the cost of the coal, respectively;
active power is generated for the coal-fired unit where the node i is located at the time t; c ng,i,t A unit gas supply cost coefficient of the gas source point i at the moment t; q gas,i,t The amount of air supplied to the air source point i at the time t; gas,i Q、
respectively, the limit values of the gas supply at the gas source point i.
The gas, thermal and electrical loads of the gas thermoelectric coupling system used in the present invention are shown in fig. 2-4.
The results calculated by the optimized scheduling method of the present invention are shown in tables 1-3. The power consumed by the electric boiler all day is 647MW, the power generated by the CHP unit is 3.89 x 103MW, and the waste heat boiler can provide 526.21MW heat power according to the heat supply proportion. The total daily air supply amount of a natural gas source point is 146.78Mm3, wherein 143.42Mm3 is used for supplying air load, and 3.35Mm3 is used for supplying power to a natural gas unit for generating power.
TABLE 1 heating and regeneration temperatures of nodes of thermodynamic system in different periods of time
Table 2 flow rates of partial branch lines of natural gas system at different time periods
TABLE 3 partial branch of electric power system different time period tidal current
Fig. 7 is a schematic structural diagram of an apparatus according to an embodiment of the present invention, where the embodiment of the present invention provides services for implementing the optimal scheduling method for the gas-thermoelectric coupling system SCUC according to the above embodiment of the present invention. Fig. 7 illustrates a block diagram of an exemplary device 12 suitable for use in implementing embodiments of the present invention. The device 12 shown in fig. 7 is only an example and should not bring any limitation to the function and scope of use of the embodiments of the present invention.
As shown in FIG. 7, device 12 is in the form of a general purpose computing device. The components of device 12 may include, but are not limited to: one or more processors or processing units 16, a system memory 28, and a bus 18 that couples various system components including the system memory 28 and the processing unit 16.
The system memory 28 may include computer system readable media in the form of volatile memory, such as Random Access Memory (RAM) 30 and/or cache memory 32. The device 12 may further include other removable/non-removable, volatile/nonvolatile computer system storage media. By way of example only, storage system 34 may be used to read from and write to non-removable, nonvolatile magnetic media (not shown in FIG. 5, and commonly referred to as a "hard drive"). Although not shown in FIG. 7, a magnetic disk drive for reading from and writing to a removable, nonvolatile magnetic disk (e.g., a "floppy disk") and an optical disk drive for reading from or writing to a removable, nonvolatile optical disk (e.g., a CD-ROM, DVD-ROM, or other optical media) may be provided. In these cases, each drive may be connected to bus 18 by one or more data media interfaces. Memory 28 may include at least one program product having a set (e.g., at least one) of program modules that are configured to carry out the functions of embodiments of the invention.
A program/utility 40 having a set (at least one) of program modules 42 may be stored, for example, in memory 28, such program modules 42 including but not limited to an operating system, one or more application programs, other program modules, and program data, each of which or some combination of which may comprise an implementation of a network environment. Program modules 42 generally carry out the functions and/or methodologies of embodiments of the invention as described.
The processing unit 16 executes various functional applications and data processing by running a program stored in the system memory 28, for example, implementing the gas thermoelectric coupling system SCUC optimized scheduling method provided by the embodiment of the present invention
In the description herein, references to the description of "one embodiment," "an example," "a specific example" or the like are intended to mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
The foregoing shows and describes the general principles, principal features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are given by way of illustration of the principles of the present invention, but that various changes and modifications may be made without departing from the spirit and scope of the invention, and such changes and modifications are within the scope of the invention as claimed.
Claims (10)
1. A gas-thermoelectric coupling system SCUC optimal scheduling method is characterized by comprising the following steps:
1) Reading in original data of the gas-thermoelectric coupling system and carrying out initialization processing; the power system adopts an IEEE39 node system, the natural gas system adopts a 20-node network, the thermal system adopts a 6-node network, and the coupling link adopts CHP and electric boiler equipment;
2) Defining system state variables and decision variables, calculating the power flow of the power system by using a direct current model, and setting power balance and power flow out-of-limit safety constraints;
3) Establishing a quadratic function piecewise linearization model;
4) Setting unit operation state parameters, and setting unit climbing and start-stop time constraints and thermal power unit output constraints;
5) Calculating the load flow of a gas network according to a natural gas linear model, calculating a correlation matrix of related node branches, and setting a natural gas consumption balance constraint, a natural gas unit output constraint and a flow and air pressure balance constraint of each node;
6) Calculating the flow of a heat supply network according to a hydraulic model and a thermodynamic model, calculating a related node branch incidence matrix, setting flow, heat and temperature balance constraints according to energy conservation and a thermodynamic theorem, and adding heat load and heat power balance constraints of each node;
7) The method comprises the following steps that power balance of a coupling link and equipment output constraint are set, the combined heat and power CHP comprises a gas turbine set and a waste heat boiler, power generation and heat supply are achieved, the electric boiler is used as a heat source, and the purpose of system optimization scheduling is achieved by adjusting the output of the coupling equipment;
8) Establishing an SCUC optimization scheduling model of the gas-thermoelectric coupling system to calculate an objective function value, adopting a piecewise linearization processing method for a nonlinear function, programming through Matlab software, and solving an optimal value by using a Cplex solver.
2. The gas-thermoelectric coupling system SCUC optimal scheduling method as recited in claim 1, wherein the system state variables in the step 2) include voltage and phase angle in the electric power system, node pressure of the natural gas system, flow rate and temperature of the thermodynamic system, CHP efficiency, low heating value and electric boiler efficiency; the system decision variables comprise the output of a generator in the power system, the transformation ratio of a natural gas compressor, the output of a gas source point, the output of a heat source of a thermodynamic system and the temperature; calculating the power flow of the power system by using a direct current model, and setting power balance and power flow out-of-limit safety constraints; the branch and node active power flow balance and constraint are as follows:
in the formula, P ij Is the active power flow of the line i, j; theta.theta. i 、θ j The voltage phase angles of nodes at two ends of the line are respectively; x is a radical of a fluorine atom ij Is a line reactance;
for the line current to allow a maximum value, ij θ、
is a voltage phase angle limit value; m is a group of G 、M EB Coefficient matrixes of nodes where the generator set and the electric boiler are located are respectively set; the electricity consumption of the electric boiler is regarded as the load,
is the node load of the electric boiler,
is the node load;
the power is output by the generator set,
respectively representing the upper limit value and the lower limit value of the output of the generator; n is node weaveNo. B ij A power system susceptance matrix.
3. The optimal scheduling method for the SCUC of the gas-thermoelectric coupling system according to claim 1, wherein the piecewise linearization model of the quadratic function curve in the step 3) is as follows:
In the formula, m is the number of the sections divided by the variable value interval; 1 is a segment node number;
variable values segmented for t period 1;
the total length of the first 1 section of the variable interval;
for the segmentation value of the 1 st subinterval, l ∈ {1,2, \8230;, m };
1 segment of binary variable;
is the slope of segment 1; w is a ij,t Is the natural gas flow between nodes i, j.
4. The optimal scheduling method for the gas-thermoelectric coupling system SCUC according to claim 1, wherein the set operation state parameters, the set ramp and start-stop time constraints, and the set output constraints are set in step 4), as follows:
in the above formula, the first and second carbon atoms are,
the generated power of the unit j at the time t is obtained;
the binary state variable of the unit is represented, the value of 1 represents that the unit j is in the running state at the moment of t, and the value of 0 represents that the unit j is in the shutdown state; j Pthe minimum value of the output of the unit j is shown,
representing the maximum value of the j output of the unit;
representing the maximum value that the unit j can obtain at the moment t;
the upper formula is climbing and start-up restriction, RU j For unit j ramp limit, SU j Starting a limit value for the unit j;
the above formula is shutdown constraint, SD j Stopping the limiting value for the unit j;
the upper form being downhill constraint, RD j A downhill limit value for the unit j;
the above formula is the minimum start time constraint of the unit, UT j The minimum startup time of the unit j;
the starting time of the unit j before scheduling;
the running state of the unit j before dispatching; g j Time for unit j to have to run;
the above formula is the minimum shutdown time constraint of the unit, DT j The minimum shutdown time of the unit j;
the shutdown time of the unit j before scheduling; l is j For the time that unit j must be shut down.
5. The optimal scheduling method for the gas-thermoelectric coupling system SCUC according to claim 1, wherein the step 5) is as follows according to the law of fluid mechanics:
w ij,t +w ji,t ＝0
in the above formula, w ij,t For the flow of natural gas pipeline node i to node j,
respectively injecting the gas quantity of the j point and the gas consumption of the gas unit into the gas source pointVolume, air load at point j; z (j) and v (j) are respectively a pipeline set taking the node j as a last node and taking the node j as a first node; c ij Is the constant of the pipe i, j; psi i,t Is the air pressure at node i;ψ、
respectively an upper limit and a lower limit of air pressure; ij w、
respectively an upper limit and a lower limit of the pipeline flow; sgn (psi) i,t ,ψ j,t ) And a sign function representing the natural gas flow direction in the pipeline, wherein the gas flow flows from the node with high pressure to the node with low pressure.
6. The optimal scheduling method for the gas-thermoelectric coupling system SCUC according to claim 1, wherein the step 6) is as follows:
Φ i ＝C p m i (T i ms -T i mr )
in the above-mentioned formula, the compound has the following structure,
respectively the outlet water temperature and the inlet water temperature of a pipeline b in the water supply system;
the outlet water temperature and the inlet water temperature of a pipeline b in a water return system are respectively; t is i ms 、T i mr Respectively the mixed temperature of the water supply system and the return water system at the node i; ms is b 、mr b Water flow in the water supply system and the backwater system pipeline b respectively;
respectively, a pipeline set taking the node i as a head end and a pipeline set taking the node i as a tail end; phi i Is the thermal load power of node i; c p The specific heat capacity of water is 4200J/(kg x DEG C); m is i The amount of water injected into the node i;
the outlet water temperature and the inlet water temperature of the pipeline b are respectively; t is a Is ambient temperature; lambda [ alpha ] b Is the heat transfer coefficient of the pipe b; l is b Is the length of the pipe b.
7. The optimal scheduling method for the SCUC of the gas-thermoelectric coupling system according to claim 1, wherein the step 7) is represented as:
H CHP,t +H EB,t ＝H source,t
H source,t ＝H load,t +H loss,t
H CHP,t ＝η chp,h *P CHP,t
H EB,t ＝η EB *P EB,t
C gas ＝C ng *Q gas
in the above formula, H CHP,t 、H EB,t Respectively generating thermal power at t moment for the CHP waste heat boiler and the electric boiler; h source,t Providing the total thermal power for the system at the moment t; h load,t 、H loss,t Respectively representing the heat load and the heat loss at the t moment of the system; p is CHP,t Electrical power generated for the CHP genset; p is EB,t The power consumption of the electric boiler; eta EB The electric heat conversion coefficient of the electric boiler; eta chp,h CHP electric heating proportionality coefficient; eta chp,e Is the efficiency coefficient of the CHP unit; zeta LHV Is natural gas with low heat value; q gas The gas consumption of natural gas is reduced; c ng The cost per unit of natural gas; c gas The total fuel cost for natural gas.
8. The optimal scheduling method for the gas-thermoelectric coupling system SCUC according to claim 1, wherein the objective function in step 8) is:
in the formula, a i 、b i 、c i Coefficients that are a function of the cost of the coal, respectively;
active power is generated for the coal-fired unit where the node i is located at the time t; c ng,i,t Is the unit gas supply cost coefficient of the gas source point i at the moment t; q gas,i,t The air supply quantity of an air source point i at the moment t; gas,i Q、
respectively, the limit values of the gas supply at the gas source point i.
9. An apparatus, comprising: one or more processors and memory for storing one or more programs, the one or more programs for execution by the one or more processors to implement the gas thermoelectric coupling system SCUC optimized scheduling method as recited in any of claims 1-8.
10. A storage medium, comprising: the computer executable instructions, when executed by a computer processor, are for performing the gas-thermoelectric coupling system SCUC optimized scheduling method of any of claims 1-8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011125210.5A CN112288592B (en) | 2020-10-20 | 2020-10-20 | Gas-thermal coupling system SCUC optimization scheduling method, device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011125210.5A CN112288592B (en) | 2020-10-20 | 2020-10-20 | Gas-thermal coupling system SCUC optimization scheduling method, device and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112288592A CN112288592A (en) | 2021-01-29 |
| CN112288592B true CN112288592B (en) | 2022-10-25 |
Family
ID=74423443
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011125210.5A Active CN112288592B (en) | 2020-10-20 | 2020-10-20 | Gas-thermal coupling system SCUC optimization scheduling method, device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112288592B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113222230A (en) * | 2021-04-29 | 2021-08-06 | 中国石油大学(北京) | Flow distribution method and device of natural gas pipe network under accident condition |
| CN113221358B (en) * | 2021-05-13 | 2022-09-27 | 浙江大学 | Standby output optimization method of electric-gas coupling system based on reliability parameters |
| CN113487342A (en) * | 2021-06-03 | 2021-10-08 | 华中科技大学 | Decision method for thermodynamic system to participate in electric power market as benefit whole |
| CN113659569A (en) * | 2021-08-09 | 2021-11-16 | 国网北京市电力公司 | Day-ahead optimal scheduling method and system for power system |
| CN114611823B (en) * | 2022-03-23 | 2022-11-08 | 特斯联科技集团有限公司 | Optimized dispatching method and system for electricity-cold-heat-gas multi-energy-demand typical park |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110007600A (en) * | 2019-04-24 | 2019-07-12 | 四川大学 | A kind of constrained multipotency stream coordinated scheduling aid decision-making system of tool |
| CN110322051A (en) * | 2019-06-06 | 2019-10-11 | 国网浙江省电力有限公司经济技术研究院 | Consider the integrated energy system Optimal Configuration Method of N-1 security constraint |
| CN111564873B (en) * | 2020-05-29 | 2021-10-19 | 江苏方天电力技术有限公司 | Integrated gas-electric energy system optimization method considering unit start and stop |
-
2020
- 2020-10-20 CN CN202011125210.5A patent/CN112288592B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN112288592A (en) | 2021-01-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112288592B (en) | Gas-thermal coupling system SCUC optimization scheduling method, device and storage medium | |
| CN109670730B (en) | Day-ahead economic dispatching method for comprehensive energy system | |
| Li et al. | Combined heat and power dispatch considering advanced adiabatic compressed air energy storage for wind power accommodation | |
| CN109919399B (en) | Day-ahead economic dispatching method and system for comprehensive energy system | |
| CN109740955A (en) | A kind of electric-gas integrated energy system planing method counted and improve staged carbon transaction | |
| Zhang et al. | Optimal clean heating mode of the integrated electricity and heat energy system considering the comprehensive energy-carbon price | |
| CN113256045B (en) | Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty | |
| CN107358345A (en) | The distributed triple-generation system optimizing operation method of meter and dsm | |
| Ding et al. | Simulation study on a novel solar aided combined heat and power system for heat-power decoupling | |
| Giap et al. | Waste heat utilization in reversible solid oxide fuel cell systems for electrical energy storage: fuel recirculation design and feasibility analysis | |
| Hu et al. | Optimal control of a hydrogen microgrid based on an experiment validated P2HH model | |
| Bo et al. | Hyperfine optimal dispatch for integrated energy microgrid considering uncertainty | |
| CN112836882A (en) | Regional comprehensive energy system operation optimization method considering equipment load rate change | |
| Liu | Optimal scheduling strategy of electricity‐heat‐hydrogen integrated energy system under different operating modes | |
| Cao et al. | Energy, exergy and economic (3E) analysis and multi-objective optimization of a combined cycle power system integrating compressed air energy storage and high-temperature thermal energy storage | |
| Yun et al. | Low-carbon optimal dispatch of integrated energy system considering the operation of oxy-fuel combustion coupled with power-to-gas and hydrogen-doped gas equipment | |
| CN111008760A (en) | Double-layer optimization method for combined cooling heating and power micro-grid system | |
| Mo et al. | Analysis of performance and suitable users of CCHP systems with active thermal energy storage | |
| CN112862158B (en) | Day-ahead optimal scheduling method and system for comprehensive energy system | |
| Man et al. | State estimation for integrated energy system containing electricity, heat and gas | |
| Zhang et al. | Low-carbon economic dispatch of integrated energy system based on liquid carbon dioxide energy storage | |
| CN110377973B (en) | Construction method of standard linear comprehensive energy system model | |
| Wang et al. | Configuration method for combined heat and power plants with flexible electricity regulation | |
| Carlson et al. | Nuclear power coupled with thermal energy storage: Impact of technical performance on economics in an exemplary electricity grid | |
| CN110716429A (en) | Control method and device of combined cooling heating and power system, computer and storage medium |
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 |