CN112731873B - Coordination control method for blast furnace combined cycle system and post-combustion carbon capture system - Google Patents

Coordination control method for blast furnace combined cycle system and post-combustion carbon capture system Download PDF

Info

Publication number
CN112731873B
CN112731873B CN202011513399.5A CN202011513399A CN112731873B CN 112731873 B CN112731873 B CN 112731873B CN 202011513399 A CN202011513399 A CN 202011513399A CN 112731873 B CN112731873 B CN 112731873B
Authority
CN
China
Prior art keywords
blast furnace
controller
combined cycle
parameter
reboiler
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
Application number
CN202011513399.5A
Other languages
Chinese (zh)
Other versions
CN112731873A (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.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN202011513399.5A priority Critical patent/CN112731873B/en
Publication of CN112731873A publication Critical patent/CN112731873A/en
Application granted granted Critical
Publication of CN112731873B publication Critical patent/CN112731873B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/418Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
    • G05B19/41865Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by job scheduling, process planning, material flow
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/32Operator till task planning
    • G05B2219/32252Scheduling production, machining, job shop
    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Quality & Reliability (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

The invention provides a blast furnace combined cycle system and burnt carbonThe coordination control method of the capture system can enhance the load tracking capability of the blast furnace combined cycle system on the premise of generating less influence on the carbon capture system after combustion. The coordinated control method of the invention adjusts the steam extraction flow of the reboiler in advance through the deviation of the system output power and the set value, and adjusts the opening of the blast furnace gas bypass valve and CO through the deviation of the steam extraction flow of the reboiler and the set value2The flow of the barren solution is adjusted by the deviation of the capture rate and a set value, so that the output power and CO are controlled by the blast furnace gas combined circulation system and the carbon capture system after combustion in a mutually coordinated and matched manner2The collection rate. The method of the invention improves the load response speed of the blast furnace combined cycle system and enhances the operation flexibility on the premise of having little influence on the carbon capture system after combustion.

Description

Coordination control method for blast furnace combined cycle system and post-combustion carbon capture system
Technical Field
The invention belongs to the technical field of thermal control, and particularly relates to a coordination control method of a blast furnace combined cycle system and a carbon capture system after combustion.
Background
Under the situation that climate change situation is getting more severe, reducing carbon emission has become a key problem for sustainable development of human society and economy. The steel industry, as a typical energy-intensive industry, is the key industry for greenhouse gas emissions, accounting for about 6% of global greenhouse gas emissions. The carbon emission of iron and steel enterprises in China accounts for 13% -15% of the total national emission, and the emission reduction of the industry is reluctant. Under the background that renewable energy sources are not fully deployed, a Carbon Capture and Storage (CCS) technology is the only effective way for realizing decarburization in the steel industry at present.
Blast Furnace Gas (BFG) is the byproduct Gas with the highest yield in steel plants, and is also the largest carbon emission source, and how to reasonably utilize BFG and reduce CO2 emission is the current research focus. The Blast Furnace Gas-based Gas-steam Combined Cycle (BFG-CCGT) is a way of efficiently utilizing Gas resources, avoids heat exchange loss of a boiler compared with the traditional Gas-fired boiler power generation, and has the thermoelectric conversion efficiency of 40-45% under the condition of no external heat supply. And post-combustion carbon capture based on chemisorption in terms of carbon abatement (Post-comfort Carbon Capture, PCC) is the most mature technology to develop due to its advantages of low cost and easy modification. Kevin et al, "carbon capture and utilization in the iron and steel industry: in the chemical engineering challenge and opportunity (Carbon capture and optimization in the steel industry: charles and opportunities for chemical engineering) "it is pointed out that the integration of a PCC system into a BFG-CCGT is more favorable for CO in consideration of the dispersed arrangement of a blast furnace and the secondary Carbon emission of the smoke discharged by a generator set2Large scale concentrated capture.
The integrated system has a complex coupling relationship between variables, and has a high requirement on the peak regulation capability of the system under the condition that renewable energy sources such as wind, light and the like are greatly accessed in the future, so that a coordination control method which effectively enhances the load tracking performance of a power plant and realizes the flexible coordination of the relationship between electricity and carbon of the system on the premise of generating less influence on a carbon capture system is necessary.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the coordination control method of the blast furnace combined cycle system and the burnt carbon capture system is provided, and the load tracking capability of the blast furnace combined cycle system can be enhanced on the premise of generating small influence on the burnt carbon capture system.
In order to solve the technical problem, the invention provides a coordination control method of a blast furnace combined cycle system and a carbon capture system after combustion, which comprises the following steps:
step 10), acquiring a set value of steam extraction flow of a reboiler;
step 30), adjusting the opening of a blast furnace gas bypass valve by adopting a first controller according to the deviation between the steam extraction flow of the reboiler and the set value of the steam extraction flow of the reboiler; adjusting the opening degree of the inlet guide vane by adopting a second controller according to the deviation between the gas turbine exhaust temperature of the blast furnace combined cycle system and the gas turbine exhaust temperature set value; adjusting the steam extraction flow of the reboiler by adopting a third controller according to the deviation of the system output power of the blast furnace combined cycle system and a power set value; CO from post-combustion carbon capture system2The deviation between the collection rate and the set collection rate is adjusted by a fourth controllerLean solution flow rate; the steam extraction flow of the reboiler finally reaches the steam extraction flow set value of the reboiler, so that the temperature of the reboiler is in an allowable range.
As a further improvement of the embodiment of the present invention, the step 10) specifically includes:
step 101) acquiring all-working-condition steady-state operation data of a blast furnace combined cycle system and a carbon capture system, wherein the operation data comprises blast furnace gas components of the blast furnace combined cycle system, system output power of the blast furnace combined cycle system, gas turbine exhaust temperature of the blast furnace combined cycle system, and CO of the carbon capture system after combustion2The capture rate, the reboiler temperature and the reboiler extraction flow rate; the collected operation data is divided into training data and test data;
102) mixing blast furnace gas components of the blast furnace combined cycle system, system output power of the blast furnace combined cycle system, gas turbine exhaust temperature of the blast furnace combined cycle system and CO of a carbon capture system after combustion2The capture rate and the reboiler temperature are used as input variables of the BP neural network model, the steam extraction flow of the reboiler is used as an output variable of the BP neural network model, the BP neural network model is trained by adopting training data, and the BP neural network model is verified by adopting test data;
and 103) collecting blast furnace gas components of the current blast furnace combined circulation system, and calculating by using the BP neural network model to obtain a reboiler steam extraction flow set value by combining a power set value, a gas turbine exhaust temperature set value, a capture rate set value and a reboiler optimal temperature value.
As a further improvement of the embodiment of the present invention, the first controller, the second controller, and the third controller are all PI controllers.
As a further improvement of the embodiment of the present invention, between the step 10) and the step 30), there is further included:
and 20) optimizing parameters of the first controller, the second controller and the third controller based on a multi-objective genetic algorithm to meet the preset control quantity constraint, so as to obtain the optimal parameters of the first controller, the second controller and the third controller.
As a further improvement of the embodiment of the present invention, the preset control amount constraint is as shown in formula (1):
Figure GDA0003472301560000041
in the formula u1Indicates the opening degree, Deltau, of the blast furnace gas bypass valve1Represents the rate of change of the opening of the blast furnace gas bypass valve u2Representing inlet guide vane opening, Δ u2Representing the rate of change of inlet guide vane opening, u3Denotes the reboiler extraction flow, Δ u3Indicating the rate of change of the reboiler extraction flow, u4Denotes the lean flow rate, Δ u4Indicating a rate of change of the lean liquid flow; u. of1minIndicates the control quantity u1Lower limit value u1maxIndicates the control quantity u1Upper limit value, Δ u1minIndicates the control quantity u1Lower limit of rate of change, Δ u1maxIndicates the control quantity u1Upper limit of rate of change, u2minIndicates the control quantity u2Lower limit value u2maxIndicates the control quantity u2Upper limit value, Δ u2minIndicates the control quantity u2Lower limit of rate of change, Δ u2maxIndicates the control quantity u2Upper limit of rate of change, u3minIndicates the control quantity u3Lower limit value u3maxIndicates the control quantity u3Upper limit value, Δ u3minIndicates the control quantity u3Lower limit of rate of change, Δ u3maxIndicates the control quantity u3Upper limit of rate of change, u4minIndicates the control quantity u4Lower limit value u4maxIndicates the control quantity u4Upper limit value, Δ u4minIndicates the control quantity u4Lower limit of rate of change, Δ u4maxIndicates the control quantity u4An upper limit value of the rate of change.
As a further improvement of the embodiment of the present invention, the step 20) specifically includes:
step 201) generating an initial population:
generating initial populations of n individuals which are uniformly distributed by taking a preset parameter range as a constraint, wherein the initial population of each individual represents a group of initial parameters; the initial parameters comprise parameters of a first controller, parameters of a second controller and parameters of a third controller; n represents an integer greater than 0;
step 202) selecting a parent individual:
setting an objective function as shown in equations (2) and (3):
Figure GDA0003472301560000042
Figure GDA0003472301560000051
expression (2) represents the tracking performance for the set value, and expression (3) represents the control overshoot;
in the formula, r1Indicates the power set point, r2Indicating the set value of the exhaust gas temperature of the combustion engine, y1Representing the system output power, y, of the combined cycle system of the blast furnace2The method comprises the steps of (1) representing the gas turbine exhaust temperature of a blast furnace combined cycle system, wherein n represents the number of sampling points; r is1bIndicating the power set point after change, r2bIndicating a changed exhaust gas temperature set value, alpha1Representing the proportion of a system output power control target in an objective function, alpha, of a blast furnace combined cycle system2Representing the proportion of a combustion engine exhaust temperature control target of a blast furnace combined cycle system in a target function;
selecting m parent individuals in an initial population by adopting a championship competition ordering method based on non-dominated ordering and crowded distance calculation; m represents an integer greater than 0;
step 203), crossing and varying parent individuals to generate new individuals to obtain a new population;
wherein, the cross formula is shown as formula (4):
child parent1+ rand Ratio (parent2-parent1) formula (4)
In the formula, child represents a child individual, parent1 and parent2 represent a parent individual, rand represents a random number between 0 and 1, and Ratio represents a crossover operator;
step 204) combining the new population and the current population into an expanded population, and calculating the grades and crowding distances of all individuals;
step 205) eliminating the individuals with high level and small crowding distance in the expanded population by a pruning algorithm, and pruning the individuals into a new population of n individuals, wherein the parameters corresponding to the individuals have better load tracking performance compared with the previous generation;
step 206) judging whether the termination condition is met, if not, repeating the steps 32 to 36;
step 207) obtaining preferred parameters of the first controller, the second controller and the third controller.
As a further improvement of the embodiment of the present invention, the preset parameter range is as shown in formula (5):
Figure GDA0003472301560000061
in the formula, P1Representing the P parameter, I, of the first controller1Denotes the I parameter, P, of the first controller2Representing the P parameter, I, of the second controller2Denotes the I parameter, P, of the second controller3P parameter, I, representing the third controller3An I parameter representing a third controller; p1minDenotes the lower limit of the first parameter, P1maxDenotes the upper limit value of the first parameter, I1minDenotes the lower limit of the second parameter, I1maxDenotes the upper limit of the second parameter, P2minDenotes the lower limit of the third parameter, P2maxRepresents the upper limit value of the third parameter, I2minDenotes the lower limit of the fourth parameter, I2maxDenotes the upper limit value of the fourth parameter, P3minDenotes the lower limit of the fifth parameter, P3maxRepresents the upper limit value of the fifth parameter, I3minDenotes the lower limit value of the sixth parameter, I3maxAnd represents the sixth parameter upper limit value.
As a further improvement of the embodiment of the present invention, the termination condition is that a preset maximum algebra is reached or an error satisfies a condition.
Compared with the prior art, the technical scheme of the invention has the following beneficial effects:
the embodiment of the invention providesThe method for coordinately controlling the combined blast furnace circulation system and the post-combustion carbon capture system comprises the steps of adjusting the steam extraction flow of a reboiler in advance through the deviation of the output power of the combined blast furnace circulation system and a set value, adjusting the opening of a gas bypass valve of the blast furnace and the opening of a CO bypass valve through the deviation of the steam extraction flow of the reboiler and the set value2The flow of the barren solution is adjusted by the deviation of the capture rate and the set value, so that the output power of the blast furnace gas combined cycle system and the CO of the carbon capture system after combustion are controlled by the blast furnace gas combined cycle and the carbon capture system after combustion in a mutually coordinated and matched manner2The collection rate. Because the time scale difference between the blast furnace combined cycle system and the burnt carbon capture system is large, the steam extraction flow of the reboiler has a quick influence on the power of the blast furnace combined cycle system, the influence on the burnt carbon capture system is slow, a part of steam is lent from the burnt carbon capture system to preferentially adjust the power of the unit, the load is returned after the load adjustment is finished, and the CO of the burnt carbon capture system is subjected to CO return2The capture rate has less influence. Compared with the conventional method that the opening of the blast furnace gas bypass valve is adjusted through the deviation of the output total power and the set value, and the steam extraction flow of the reboiler is adjusted through the deviation of the reboiler temperature and the set value, the method provided by the embodiment of the invention can improve the load response speed of the blast furnace combined circulation system and enhance the operation flexibility on the premise of having little influence on the post-combustion carbon capture system.
Drawings
FIG. 1(a) is a block flow diagram of a method for coordinated control of a combined blast furnace circulation system and a post-combustion carbon capture system in accordance with an embodiment of the present invention;
FIG. 1(b) is a schematic structural diagram of a combined blast furnace circulation system and a post-combustion carbon capture system according to an embodiment of the present invention;
FIG. 2 is a diagram of a BP neural network model structure in the method according to the embodiment of the present invention;
FIG. 3 is a Pareto frontier chart of an optimal solution set obtained by optimizing PI parameters through a multi-objective genetic algorithm in the method of the embodiment of the invention;
FIG. 4(a) is a graph illustrating the variation of the output power of the controlled quantity system during the variation of the power setpoint and the exhaust temperature setpoint of the combustion engine according to the exemplary method and the conventional control method of the present invention;
FIG. 4(b) is a graph of the change in engine exhaust temperature during a change in a power set point and an engine exhaust temperature set point for an embodiment of the present invention and a conventional control method;
FIG. 4(c) is a graph illustrating CO during a change in power setpoint and engine exhaust temperature setpoint, according to an exemplary method of the present invention and a conventional control method2A change curve of the trapping rate;
FIG. 4(d) is a graph showing the reboiler temperature profile during a change in the power setpoint and the turbine exhaust temperature setpoint for an exemplary method and conventional control method of the present invention;
FIG. 5(a) is a graph showing the variation of the steam extraction rate of the reboiler during the variation of the power setpoint and the exhaust temperature setpoint of the gas turbine according to the exemplary method and the conventional control method of the present invention;
FIG. 5(b) is a graph showing the change in lean flow rate during a change in power set point and engine exhaust temperature set point for an exemplary method of the present invention and a conventional control method;
FIG. 5(c) is a graph showing the variation of the opening of the blast furnace gas bypass valve during the variation of the power set point and the set point of the exhaust temperature of the gas turbine according to the method of the embodiment of the present invention and the conventional control method;
FIG. 5(d) is a graph showing the variation of the inlet guide vane opening during the variation of the power set point and the exhaust temperature set point for a combustion engine according to the method of the present invention and the conventional control method.
Detailed Description
The technical solution of the present invention will be described in detail below with reference to the accompanying drawings.
An embodiment of the present invention provides a method for coordinated control of a combined cycle system of a blast furnace and a post-combustion carbon capture system, as shown in fig. 1(a) and 1(b), including:
and step 10) acquiring a set value of the steam extraction flow of the reboiler.
Step 30) according to the steam extraction flow u of the reboiler3Adjusting the opening u of a blast furnace gas bypass valve (BPV) by a first controller according to the deviation of the set value of the steam extraction flow of the reboiler1(ii) a According to the exhaust temperature y of the combustion engine of the blast furnace combined cycle system2Deviation from set value of exhaust gas temperature of combustion engine, using second controllerAdjusting Inlet Guide Vane (IGV) opening u2(ii) a According to the system output power y of the blast furnace combined cycle system1Deviation from the power set point, a third controller is adopted to adjust the steam extraction flow u of the reboiler3(ii) a CO from post-combustion carbon capture system2Trapping rate y3Adjusting lean solution flow u using a fourth controller in accordance with the deviation from the collection rate set value4. Extraction flow u through reboiler3And blast furnace gas bypass valve (BPV) opening u1And controlling the system output power of the blast furnace combined cycle system. By Inlet Guide Vane (IGV) opening u2Controlling combustion engine exhaust temperature y of blast furnace combined cycle system2. Through lean solution flow u4Controlling capture rate y of post-combustion carbon capture system3. The steam extraction of the reboiler participates in load regulation and then gradually returns to the steam extraction flow set value of the reboiler, so that the temperature y of the reboiler is ensured4Within the allowable range.
The method of the embodiment of the invention adjusts the steam extraction flow of the reboiler in advance through the deviation of the output power of the blast furnace combined cycle system and the set value, and adjusts the BPV opening and the CO opening degree through the deviation of the steam extraction flow of the reboiler and the set value2The deviation of the capture rate and a set value adjusts the flow rate of a barren solution, the steam extraction flow rate of a reboiler and the opening degree of a blast furnace gas bypass valve control the system output power of the blast furnace combined cycle system, the flow rate of the barren solution controls the capture rate of the burnt carbon capture system, so that the blast furnace gas combined cycle system and the burnt carbon capture system are mutually coordinated and matched to control the output power of the blast furnace combined cycle system and the CO of the burnt carbon capture system2The collection rate. Because the time scale difference between the blast furnace combined cycle system and the burnt carbon capture system is large, the steam extraction flow of the reboiler has a quick influence on the power of the blast furnace combined cycle system, the influence on the burnt carbon capture system is slow, a part of steam is lent from the burnt carbon capture system to preferentially adjust the power of the unit, the load is returned after the load adjustment is finished, and the CO of the burnt carbon capture system is subjected to CO return2The capture rate has less influence. Compared with the conventional method that the BPV opening is adjusted through the deviation of the output total power and the set value, and the reboiler temperature and the set value deviation are used for adjusting the reboiler steam extraction flow, the method provided by the embodiment of the invention can be used for combustionOn the premise that the influence of the rear carbon capture system is small, the load response speed of the blast furnace combined cycle system is increased, and the operation flexibility is enhanced.
As a preferred example, the reboiler extraction flow rate set point is obtained by the following steps:
step 101) acquiring all-working-condition steady-state operation data of a blast furnace combined cycle system and a burnt carbon capture system, wherein the operation data comprises blast furnace gas components of the blast furnace combined cycle system, system output power of the blast furnace combined cycle system, gas turbine exhaust temperature of the blast furnace combined cycle system and CO of the burnt carbon capture system2The capture rate, the temperature of the reboiler and the steam extraction flow rate of the reboiler; the collected operation data is divided into training data and test data;
102) mixing blast furnace gas components of the blast furnace combined cycle system, system output power of the blast furnace combined cycle system, gas turbine exhaust temperature of the blast furnace combined cycle system and CO of a carbon capture system after combustion2The capture rate and the temperature of the reboiler serve as input variables of the BP neural network model, the steam extraction flow of the reboiler serves as output variables of the BP neural network model, as shown in FIG. 2, the BP neural network model is trained by adopting training data, and the BP neural network model is verified by adopting test data;
and 103) collecting blast furnace gas components of the current blast furnace combined circulation system, and calculating by using a BP (back propagation) neural network model to obtain a reboiler steam extraction flow set value by combining a power set value, a gas turbine exhaust temperature set value, a capture rate set value and a reboiler optimal temperature value.
According to the method, according to the blast furnace gas components of the current blast furnace combined cycle system, the power set value, the gas turbine exhaust temperature set value, the capture rate set value and the reboiler optimal temperature value are combined, the reboiler steam extraction flow set value calculated by the established BP neural network model is closer to the actual system operation condition, and the reboiler temperature can be effectively ensured to be within the allowable range.
Preferably, the first controller, the second controller and the third controller are all PI controllers. The PI controller is simple in structure and control method and easy to realize in engineering.
As a preferred example, the method of this embodiment further includes: and optimizing parameters of the first controller, the second controller and the third controller based on a multi-objective genetic algorithm so as to meet the preset control quantity constraint and obtain the optimal parameters of the first controller, the second controller and the third controller.
Further, the preset control amount constraint is as shown in formula (1):
Figure GDA0003472301560000101
in the formula u1Indicates the opening degree, Deltau, of the blast furnace gas bypass valve1Represents the rate of change of the opening of the blast furnace gas bypass valve u2Representing inlet guide vane opening, Δ u2Representing the rate of change of inlet guide vane opening, u3Denotes the reboiler extraction flow, Δ u3Indicating the rate of change of the reboiler extraction flow, u4Denotes the lean flow rate, Δ u4Indicating a rate of change of the lean liquid flow; u. of1minIndicates the control quantity u1Lower limit value u1maxIndicates the control quantity u1Upper limit value, Δ u1minIndicates the control quantity u1Lower limit of rate of change, Δ u1maxIndicates the control quantity u1Upper limit of rate of change, u2minIndicates the control quantity u2Lower limit value u2maxIndicates the control quantity u2Upper limit value, Δ u2minIndicates the control quantity u2Lower limit of rate of change, Δ u2maxIndicates the control quantity u2Upper limit of rate of change, u3minIndicates the control quantity u3Lower limit value u3maxIndicates the control quantity u3Upper limit value, Δ u3minIndicates the control quantity u3Lower limit of rate of change, Δ u3maxIndicates the control quantity u3Upper limit of rate of change, u4minIndicates the control quantity u4Lower limit value u4maxIndicates the control quantity u4Upper limit value, Δ u4minIndicates the control quantity u4Lower limit of rate of change, Δ u4maxIndicates the control quantity u4An upper limit value of the rate of change.
In the method, the physical limits and the actual operating condition range of control elements such as a valve, a feed pump and the like in an actual system are considered, and the rate constraint and the upper and lower limit constraints of the control variable are set, so that the actual operating condition is better met, and the optimized controller parameters can be used in the actual system.
As a preferred example, the method for optimizing the parameters of the first controller, the second controller and the third controller based on the multi-objective genetic algorithm to satisfy the preset control constraints to obtain the preferred parameters of the first controller, the second controller and the third controller specifically includes:
step 201) generating an initial population:
generating initial populations of n individuals which are uniformly distributed by taking a preset parameter range as a constraint, wherein the initial population of each individual represents a group of initial parameters; the initial parameters comprise parameters of a first controller, parameters of a second controller and parameters of a third controller;
step 202) selecting a parent individual:
setting an objective function as shown in equations (2) and (3):
Figure GDA0003472301560000111
Figure GDA0003472301560000112
expression (2) represents the tracking performance for the set value, and expression (3) represents the control overshoot;
in the formula, r1Indicates the power set point, r2Indicating the set value of the exhaust gas temperature of the combustion engine, y1Representing the system output power, y, of the combined cycle system of the blast furnace2The method comprises the steps of (1) representing the gas turbine exhaust temperature of a blast furnace combined cycle system, wherein n represents the number of sampling points; r is1bIndicating the power set point after change, r2bIndicating a changed exhaust gas temperature set value, alpha1Representing the proportion of a system output power control target in an objective function, alpha, of a blast furnace combined cycle system2Indicating combustion in a combined cycle system of a blast furnaceThe engine exhaust temperature control target accounts for the proportion of a target function;
selecting m parent individuals in an initial population by adopting a championship competition ordering method based on non-dominated ordering and crowded distance calculation;
step 203), crossing and varying parent individuals to generate new individuals to obtain a new population;
wherein, the cross formula is shown as formula (4):
child parent1+ rand Ratio (parent2-parent1) formula (4)
In the formula, child represents a child individual, parent1 and parent2 represent a parent individual, rand represents a random number between 0 and 1, and Ratio represents a crossover operator;
step 204) combining the new population and the current population into an expanded population, and calculating the grades and crowding distances of all individuals;
step 205) eliminating the individuals with high level and small crowding distance in the expanded population by a pruning algorithm, and pruning the individuals into a new population of n individuals, wherein the parameters corresponding to the individuals have better load tracking performance compared with the previous generation;
step 206) judging whether the termination condition is met, if not, repeating the step 202 to the step 206;
step 207) obtaining preferred parameters of the first controller, the second controller and the third controller.
Further, the preset parameter range is as shown in formula (5):
Figure GDA0003472301560000121
in the formula, P1Representing the P parameter, I, of the first controller1Denotes the I parameter, P, of the first controller2Representing the P parameter, I, of the second controller2Denotes the I parameter, P, of the second controller3P parameter, I, representing the third controller3An I parameter representing a third controller; p1minDenotes the lower limit of the first parameter, P1maxDenotes the upper limit value of the first parameter, I1minDenotes the lower limit of the second parameter, I1maxWhich represents the upper value of the second parameter,P2mindenotes the lower limit of the third parameter, P2maxRepresents the upper limit value of the third parameter, I2minDenotes the lower limit of the fourth parameter, I2maxDenotes the upper limit value of the fourth parameter, P3minDenotes the lower limit of the fifth parameter, P3maxRepresents the upper limit value of the fifth parameter, I3minDenotes the lower limit value of the sixth parameter, I3maxAnd represents the sixth parameter upper limit value.
The embodiment of the invention sets the parameter range, can accelerate the optimization speed and obtain the optimal controller parameter as soon as possible.
Further, the termination condition is that a preset maximum algebra is reached or the error meets the condition.
A specific example is provided below to verify the performance of the method of the present invention.
Example 1
And step 10) acquiring a set value of the steam extraction flow of the reboiler. The method specifically comprises the following steps:
step 101) input variable is BFG component u11System output power y1Exhaust gas temperature y of the gas turbine2、CO2Trapping rate y3And reboiler temperature y4The output variable is the steam extraction flow u of the reboiler3Acquiring all-condition steady-state operation data of a blast furnace combined circulation system and a burnt carbon capture system;
step 102), selecting 70% of data for training a BP neural network model, and performing model verification on the rest 30% of data;
step 103) according to the current BFG component, the system output set value is 30.72MW, the gas turbine exhaust temperature set value is 600.04 ℃, and CO is carried out2The set value of the trapping rate is 55 percent, the optimum temperature of the reboiler is 390.4 ℃, and the set value of the steam extraction flow of the reboiler is calculated to be 17.14kg/s through a BP neural network model.
Step 20) optimizing the controller parameters based on the multi-objective genetic algorithm. The method specifically comprises the following steps:
step 201) with a preset parameter range as shown in the following formula as a constraint,
Figure GDA0003472301560000141
generating an initial random population of evenly distributed n-80 individuals, representing a set of initial parameter solutions, wherein each individual comprises an initial value comprising a PI parameter of a first PI controller, a PI parameter of a second PI controller and a PI parameter of a third PI controller;
step 202) the objective function is set to the formula (2) and the formula (3), in the formula (3), alpha1=0.8,α20.2; selecting individuals in a population by adopting a championship sorting method based on non-dominated sorting and crowded distance calculation, namely randomly selecting m to 30 individuals from a mother population, optimally selecting the individuals with low grade and large crowded distance, and continuously repeating the steps until a new population reaches the original population scale, wherein the grade is determined based on the non-dominated sorting, and the smaller the target value is, the lower the grade is; the crowding distance represents the distance between each point and the adjacent point in the same level, and the bigger the crowding distance is, the better the diversity of the algorithm is;
step 203) according to formula (4), where Ratio is 0.9; crossing and varying parent individuals, randomly generating a self-adaptive direction relative to the success or failure of the previous generation by a variation algorithm, selecting the direction and the step length to meet a preset parameter range, and generating new individuals;
step 204) combining the new population and the current population into an expanded population, and calculating the grades and crowding distances of all individuals;
step 205) eliminating the individuals with high grade and small crowding distance in the expanded population by a pruning algorithm, and pruning the individuals into a new population with n equal to 80 individuals, wherein PI parameters corresponding to the individuals generally have better load tracking performance compared with the PI parameters of the previous generation;
step 206) judging whether the preset maximum algebra is reached or the error meets the condition, if not, repeating the step 202) -the step 206).
Step 207) of final normalization to obtain a Pareto (Pareto) front map, as shown in fig. 3. Each point in the pareto frontier graph can be an optimal point, the tracking performance and the overshoot are comprehensively considered, and the point with the minimum sum of the two indexes is selected to be a punctuation point in the graph. And satisfies a predetermined control constraint as follows:
Figure GDA0003472301560000151
obtaining parameters of a first PI controller: p1=-1.79,I1-0.0004; parameters of the second PI controller: p2=-0.22,I2-0.0064; parameters of the third PI controller: p3=-1.95,I3-0.11. Setting parameters of a fourth PI controller: p4=23.98,I4=4.80。
Step 30) according to the steam extraction flow u of the reboiler3And adjusting the BPV opening u by adopting a first PI controller according to the deviation of the set value of the steam extraction flow of the reboiler1(ii) a According to the exhaust temperature y of the combustion engine of the blast furnace combined cycle system2And the deviation of the set value of the exhaust temperature of the gas turbine, and adjusting the opening u of the IGV by adopting a second PI controller2(ii) a According to the system output power y of the blast furnace combined cycle system1Regulating the steam extraction flow u of the reboiler by a third PI controller according to the deviation from the power set value3(ii) a CO from post-combustion carbon capture system2Trapping rate y3Adjusting lean solution flow u using a fourth PI controller in response to a deviation from the capture rate set point4
Comparative example 1
The system output power of the blast furnace combined cycle system is controlled through the BPV opening degree, the gas turbine exhaust temperature of the blast furnace combined cycle system is controlled through the IGV opening degree, and the CO of the carbon capture system after combustion is controlled through the lean solution flow2And the capture rate is controlled by the reboiler steam extraction flow to control the reboiler temperature of the carbon capture system after combustion.
The following simulation experiments were performed for example 1 and comparative example 1: the carbon capture system of the gas turbine of the blast furnace gas combustion engine is stabilized at the system output power of 35.72MW, the combustion engine exhaust temperature of 595.04 ℃ and CO2At the working point with the capture rate of 50%, when t is 100s, the set value of the system output power is changed to 30.72MW, the set value of the exhaust temperature of the gas turbine is changed to 600.04 ℃ and the CO is changed2The collection rate set value was 55% and the reboiler extraction flow rate set value was 17.14kg/s, and a control simulation result was obtained.
As shown in fig. 4(a), 4(b), 4(c) and 4(d), in the figure, the improved control is example 1 and the conventional control is comparative example. It can be seen from the graph that the method of the embodiment of the invention has almost perfect tracking performance for the power set point and the exhaust temperature of the combustion engine rises faster than the conventional control method. And the temperature change of the reboiler is very small and can be almost ignored, thereby enhancing the load tracking capability and the operation flexibility of the system.
As shown in fig. 5(a), 5(b), 5(c) and 5(d), in the figure, the improved control is example 1 and the conventional control is comparative example. As can be seen from the figure, the method of the embodiment of the invention rapidly increases the steam extraction flow of the reboiler when the power set value is reduced so as to accelerate the load regulation; and simultaneously increasing the BPV opening to decrease the fuel amount, thereby decreasing the system output power. The IGV opening is decreased to decrease the compressor inlet air amount as the engine exhaust temperature set point increases, thereby increasing the exhaust temperature. Since the reboiler extraction flow rate "replaces" a portion of the BPV opening, the method of this embodiment controls a more gradual change in BPV opening than conventional control. After the load adjustment is finished, the steam extraction flow of the reboiler gradually returns to a set value to finish the appointed CO2The collection rate is required.
In summary, the method provided by the embodiment of the invention introduces auxiliary control of the steam extraction flow of the reboiler on the basis of the conventional blast furnace combined cycle control, provides a coordinated control method, optimizes the PI parameter of the controller through a multi-objective genetic algorithm, and enhances the load tracking capability of the system while finishing the set emission reduction target compared with the conventional control method.
The foregoing illustrates and describes the principles, general features, and advantages of the present 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 intended to further illustrate the principles of the invention, and that various changes and modifications may be made without departing from the spirit and scope of the invention, which is also intended to be covered by the appended claims. The scope of the invention is defined by the claims and their equivalents.

Claims (7)

1. A coordination control method for a blast furnace combined cycle system and a carbon capture system after combustion is characterized by comprising the following steps:
step 10), acquiring a set value of steam extraction flow of a reboiler;
step 30), adjusting the opening of a blast furnace gas bypass valve by adopting a first controller according to the deviation between the steam extraction flow of the reboiler and the set value of the steam extraction flow of the reboiler; adjusting the opening degree of the inlet guide vane by adopting a second controller according to the deviation between the gas turbine exhaust temperature of the blast furnace combined cycle system and the gas turbine exhaust temperature set value; adjusting the steam extraction flow of the reboiler by adopting a third controller according to the deviation of the system output power of the blast furnace combined cycle system and a power set value; CO from post-combustion carbon capture system2Adjusting the flow of the barren solution by adopting a fourth controller according to the deviation of the capture rate and the capture rate set value; the steam extraction flow of the reboiler finally reaches the steam extraction flow set value of the reboiler, so that the temperature of the reboiler is in an allowable range;
the step 10) specifically comprises:
step 101) acquiring all-working-condition steady-state operation data of a blast furnace combined cycle system and a carbon capture system, wherein the operation data comprises blast furnace gas components of the blast furnace combined cycle system, system output power of the blast furnace combined cycle system, gas turbine exhaust temperature of the blast furnace combined cycle system, and CO of the carbon capture system after combustion2The capture rate, the reboiler temperature and the reboiler extraction flow rate; the collected operation data is divided into training data and test data;
102) mixing blast furnace gas components of the blast furnace combined cycle system, system output power of the blast furnace combined cycle system, gas turbine exhaust temperature of the blast furnace combined cycle system and CO of a carbon capture system after combustion2The capture rate and the reboiler temperature are used as input variables of the BP neural network model, the steam extraction flow of the reboiler is used as an output variable of the BP neural network model, the BP neural network model is trained by adopting training data, and the BP neural network model is verified by adopting test data;
and 103) collecting blast furnace gas components of the current blast furnace combined circulation system, and calculating by using the BP neural network model to obtain a reboiler steam extraction flow set value by combining a power set value, a gas turbine exhaust temperature set value, a capture rate set value and a reboiler optimal temperature value.
2. The method of claim 1, wherein the first controller, the second controller, and the third controller are PI controllers.
3. The method of claim 2, further comprising, between the steps 10) and 30):
and 20) optimizing parameters of the first controller, the second controller and the third controller based on a multi-objective genetic algorithm to meet the preset control quantity constraint, so as to obtain the optimal parameters of the first controller, the second controller and the third controller.
4. The coordinated control method of the combined cycle system of a blast furnace and the post combustion carbon capture system as set forth in claim 3, wherein said preset control amount is constrained as shown in formula (1):
Figure FDA0003472301550000021
in the formula u1Indicates the opening degree, Deltau, of the blast furnace gas bypass valve1Represents the rate of change of the opening of the blast furnace gas bypass valve u2Representing inlet guide vane opening, Δ u2Representing the rate of change of inlet guide vane opening, u3Denotes the reboiler extraction flow, Δ u3Indicating the rate of change of the reboiler extraction flow, u4Denotes the lean flow rate, Δ u4Indicating a rate of change of the lean liquid flow; u. of1minIndicates the control quantity u1Lower limit value u1maxIndicates the control quantity u1Upper limit value, Δ u1minIndicates the control quantity u1Lower limit of rate of change, Δ u1maxIndicates the control quantity u1Upper limit of rate of change, u2minIndicates the control quantity u2Lower limit value u2maxIndicates the control quantity u2Upper limit value, Δ u2minIndicates the control quantity u2Lower limit of rate of change, Δ u2maxIndicates the control quantity u2Upper limit of rate of change, u3minIndicates the control quantity u3Lower limit value u3maxIndicates the control quantity u3Upper limit value, Δ u3minIndicates the control quantity u3Lower limit of rate of change, Δ u3maxIndicates the control quantity u3Upper limit of rate of change, u4minIndicates the control quantity u4Lower limit value u4maxIndicates the control quantity u4Upper limit value, Δ u4minIndicates the control quantity u4Lower limit of rate of change, Δ u4maxIndicates the control quantity u4An upper limit value of the rate of change.
5. The coordinated control method of the combined cycle system of a blast furnace and the carbon capture system after combustion as set forth in claim 3, wherein the step 20) specifically comprises:
step 201) generating an initial population:
generating initial populations of n individuals which are uniformly distributed by taking a preset parameter range as a constraint, wherein the initial population of each individual represents a group of initial parameters; the initial parameters comprise parameters of a first controller, parameters of a second controller and parameters of a third controller; n represents an integer greater than 0;
step 202) selecting a parent individual:
setting an objective function as shown in equations (2) and (3):
Figure FDA0003472301550000031
Figure FDA0003472301550000032
expression (2) represents the tracking performance for the set value, and expression (3) represents the control overshoot;
in the formula, r1Indicates the power set point, r2Indicating the set value of the exhaust gas temperature of the combustion engine, y1Representing the system output power, y, of the combined cycle system of the blast furnace2The method comprises the steps of (1) representing the gas turbine exhaust temperature of a blast furnace combined cycle system, wherein n represents the number of sampling points; r is1bIndicating the power set point after change, r2bIndicating a changed exhaust gas temperature set value, alpha1Representing the proportion of a system output power control target in an objective function, alpha, of a blast furnace combined cycle system2Representing the proportion of a combustion engine exhaust temperature control target of a blast furnace combined cycle system in a target function;
selecting m parent individuals in an initial population by adopting a championship competition ordering method based on non-dominated ordering and crowded distance calculation; m represents an integer greater than 0;
step 203), crossing and varying parent individuals to generate new individuals to obtain a new population;
wherein, the cross formula is shown as formula (4):
child parent1+ rand Ratio (parent2-parent1) formula (4)
In the formula, child represents a child individual, parent1 and parent2 represent a parent individual, rand represents a random number between 0 and 1, and Ratio represents a crossover operator;
step 204) combining the new population and the current population into an expanded population, and calculating the grades and crowding distances of all individuals;
step 205) eliminating the individuals with high level and small crowding distance in the expanded population by a pruning algorithm, and pruning the individuals into a new population of n individuals, wherein the parameters corresponding to the individuals have better load tracking performance compared with the previous generation;
step 206) judging whether the termination condition is met, if not, repeating the step 202) to the step 206);
step 207) obtaining preferred parameters of the first controller, the second controller and the third controller.
6. The method of claim 5, wherein the predetermined parameter range is as shown in equation (5):
Figure FDA0003472301550000041
in the formula, P1Representing the P parameter, I, of the first controller1Denotes the I parameter, P, of the first controller2Representing the P parameter, I, of the second controller2Denotes the I parameter, P, of the second controller3P parameter, I, representing the third controller3An I parameter representing a third controller; p1minDenotes the lower limit of the first parameter, P1maxDenotes the upper limit value of the first parameter, I1minDenotes the lower limit of the second parameter, I1maxDenotes the upper limit of the second parameter, P2minDenotes the lower limit of the third parameter, P2maxRepresents the upper limit value of the third parameter, I2minDenotes the lower limit of the fourth parameter, I2maxDenotes the upper limit value of the fourth parameter, P3minDenotes the lower limit of the fifth parameter, P3maxRepresents the upper limit value of the fifth parameter, I3minDenotes the lower limit value of the sixth parameter, I3maxAnd represents the sixth parameter upper limit value.
7. The coordinated control method of the combined cycle system of a blast furnace and the carbon capture system after combustion as set forth in claim 5, wherein the termination condition is that a preset maximum algebra is reached or that an error is satisfied.
CN202011513399.5A 2020-12-18 2020-12-18 Coordination control method for blast furnace combined cycle system and post-combustion carbon capture system Active CN112731873B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011513399.5A CN112731873B (en) 2020-12-18 2020-12-18 Coordination control method for blast furnace combined cycle system and post-combustion carbon capture system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011513399.5A CN112731873B (en) 2020-12-18 2020-12-18 Coordination control method for blast furnace combined cycle system and post-combustion carbon capture system

Publications (2)

Publication Number Publication Date
CN112731873A CN112731873A (en) 2021-04-30
CN112731873B true CN112731873B (en) 2022-03-08

Family

ID=75603527

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011513399.5A Active CN112731873B (en) 2020-12-18 2020-12-18 Coordination control method for blast furnace combined cycle system and post-combustion carbon capture system

Country Status (1)

Country Link
CN (1) CN112731873B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115770469A (en) * 2022-10-13 2023-03-10 昆明理工大学 Method for purifying blast furnace gas and collecting hot blast furnace carbon
CN116679572B (en) * 2023-08-03 2023-09-29 北京绿能碳宝科技发展有限公司 Carbon dioxide trapping self-learning method based on deep Q learning network

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107081040A (en) * 2017-03-17 2017-08-22 北京交通大学 One kind is based on CO after chemical chain making oxygen by air separation and burning2The electricity generation system of trapping
CN108710356A (en) * 2018-06-19 2018-10-26 东南大学 CO after the effective burning of one kind2Trap coal generating system progress control method
CN111552175A (en) * 2020-05-14 2020-08-18 东南大学 Overall optimization scheduling and rapid variable load control method for supercritical coal-fired power plant-carbon capture system after chemical adsorption combustion

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4577280A (en) * 1983-11-03 1986-03-18 Westinghouse Electric Corp. Control system for fluid flow distribution
EP1701006B1 (en) * 2005-02-22 2016-10-05 Kabushiki Kaisha Toshiba Electric power-generating and desalination combined plant and operation method of the same
EP2248999A1 (en) * 2008-12-24 2010-11-10 Alstom Technology Ltd Power plant with CO2 capture
EP2383522B1 (en) * 2010-04-28 2016-11-02 General Electric Technology GmbH Thermal integration of a carbon dioxide capture and compression unit with a steam or combined cycle plant
US9411326B2 (en) * 2012-08-21 2016-08-09 General Electric Company Plant control optimization system including visual risk display
WO2015035134A1 (en) * 2013-09-06 2015-03-12 Neuco, Inc. Combined cycle power generation optimization system
CN104049539B (en) * 2014-05-29 2016-08-17 华中科技大学 A kind of full working scope analogue system of Steam Turbine Regenerative System
CN110764419B (en) * 2019-11-15 2022-06-10 江苏方天电力技术有限公司 CO of large coal-fired power plant2Capture global scheduling and predictive control system and method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107081040A (en) * 2017-03-17 2017-08-22 北京交通大学 One kind is based on CO after chemical chain making oxygen by air separation and burning2The electricity generation system of trapping
CN108710356A (en) * 2018-06-19 2018-10-26 东南大学 CO after the effective burning of one kind2Trap coal generating system progress control method
CN111552175A (en) * 2020-05-14 2020-08-18 东南大学 Overall optimization scheduling and rapid variable load control method for supercritical coal-fired power plant-carbon capture system after chemical adsorption combustion

Also Published As

Publication number Publication date
CN112731873A (en) 2021-04-30

Similar Documents

Publication Publication Date Title
CN112633560B (en) Power station optimal scheduling method containing coal-fired cogeneration unit
CN112731873B (en) Coordination control method for blast furnace combined cycle system and post-combustion carbon capture system
CN103133147B (en) Method and apparatus for changing the power output in combustion gas turbine systems
CN107482688B (en) Scheduling optimization method for carbon capture virtual power plant
CN105787211B (en) For the Combined Cycle Heat Recovery Boiler pressure method of adjustment of combustion gas turbine deterioration
CN110764419B (en) CO of large coal-fired power plant2Capture global scheduling and predictive control system and method
CN106019935B (en) Multiple target boiler combustion optimization based on constraint fuzzy association rules
CN111552175B (en) Overall optimization scheduling and rapid variable load control method for supercritical coal-fired power plant-carbon capture system after chemical adsorption combustion
CN108710356B (en) Effective post-combustion CO2Operation control method for trapping coal-fired power generation system
CN102914966A (en) Method for dynamically setting parameters of coordinated control system on basis of preliminary coal supply control model
CN111637442B (en) Configuration self-adaptive efficient flexible clean coal-fired power generation system and operation method
CN113489024B (en) Multi-mode steam extraction auxiliary peak regulation and frequency modulation control system and method for cogeneration unit
CN115111601B (en) Multi-target boiler combustion optimization control method based on embedded algorithm fusion under variable load
CN107102545B (en) Load optimization distribution and residual gas amount correction calculation method and system for shutdown process of combined cycle unit
CN113657033A (en) Model-based multi-objective optimization method for gas turbine
CN117498393A (en) Peak regulation-heat supply-heat storage association system and cooperative control method thereof
CN108196452B (en) Heat supply unit optimal initial pressure operation control system and method based on heat balance
CN117028141A (en) Variable pitch coordinated optimization control method for wind turbine generator
CN113110316B (en) Primary frequency modulation control method for steam turbine of combined cycle unit
Fu et al. Using molten-salt energy storage to decrease the minimum operation load of the coal-fired power plant
CN113341716B (en) Large-scale coal-fired power plant CO based on artificial intelligence 2 Optimized scheduling method for trapping system
CN110970936B (en) Method for calculating primary frequency modulation performance of deep peak shaving unit
CN113587208A (en) Online quantitative calculation method and system for energy storage of steam pipe network
CN103629691B (en) Boiler combustion optimization
Flores et al. Optimizing natural gas combined cycle part load operation

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