CN111275251B - CCHP system cooling, heating and power combined optimization method containing sewage source heat pump - Google Patents

CCHP system cooling, heating and power combined optimization method containing sewage source heat pump Download PDF

Info

Publication number
CN111275251B
CN111275251B CN202010045449.5A CN202010045449A CN111275251B CN 111275251 B CN111275251 B CN 111275251B CN 202010045449 A CN202010045449 A CN 202010045449A CN 111275251 B CN111275251 B CN 111275251B
Authority
CN
China
Prior art keywords
power
gas
time
cost
gas turbine
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.)
Expired - Fee Related
Application number
CN202010045449.5A
Other languages
Chinese (zh)
Other versions
CN111275251A (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.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
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 Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN202010045449.5A priority Critical patent/CN111275251B/en
Publication of CN111275251A publication Critical patent/CN111275251A/en
Application granted granted Critical
Publication of CN111275251B publication Critical patent/CN111275251B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • 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
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Educational Administration (AREA)
  • General Health & Medical Sciences (AREA)
  • Engine Equipment That Uses Special Cycles (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

Aiming at the problem that the existing CCHP system of the traditional algorithm does not take the satisfaction rate of user requirements and the system operation cost as an optimized objective function, the invention provides a combined cooling, heating and power optimization method of the CCHP system of a heat pump containing a sewage source, which is characterized by comprising the following steps: s1, inputting original data of a combined cooling heating and power system; s2, constructing a CCHP micro-grid model; s3, constructing a target function; s4, constructing constraint conditions; and S5, solving an algorithm to obtain an optimization result. The method constructs the start-stop functions of the gas turbine and the waste heat boiler, calculates the start-stop cost of the two devices better, considers the satisfaction rate of user requirements, takes the satisfaction rate of the user requirements and the system operation cost as optimization objective functions, and optimizes the operation scheduling of the cooling, heating and power triple supply type micro-grid system.

Description

CCHP system cooling, heating and power combined optimization method containing sewage source heat pump
Technical Field
The invention belongs to the field of operation and control of power systems, and particularly relates to a method for solving the micro-grid dispatching output of a combined cooling heating and power system with a sewage source heat pump by using an improved crisscross algorithm.
Background
The CCHP system is a mature energy comprehensive utilization technology in many countries and regions, and is widely concerned and favored by governments and business industries of various countries due to the characteristics of being close to users, gradient utilization, high primary energy utilization efficiency, environmental friendliness, safe and reliable energy supply and the like. Many scholars both at home and abroad have conducted a great deal of research on the CCHP system. The heat pump is used as a low-investment and high-yield heat production device and is suitable for being used as a heat source in a CCHP system. Many scholars have also studied the energy complementation and coordination, the economy, etc. of the CCHP system with heat pump.
The sewage source heat pump is one of water source heat pumps, utilizes lower electric energy to extract heat in urban sewage or industrial sewage to recycle resources, not only saves primary energy, but also has higher economical efficiency, and the temperature of general sewage is higher than the environmental temperature, so that the heat producible by the sewage source heat pump is larger, and the sewage source heat pump also has excellent heat production capacity as heat source equipment. At present, the research on the sewage source heat pump in a CCHP system at home and abroad is less. The research on the optimization method of combined cooling heating and power of the CCHP system containing the sewage source heat pump is not sufficient.
The main purpose of the establishment of the CCHP system is to better meet the demands of users on the three loads of cold, heat and electricity, but the existing CCHP system does not take the satisfaction rate of the demands of the users and the operation cost of the system as an optimized objective function, so that the economy of the CCHP system is affected.
Disclosure of Invention
Aiming at the problem that the existing CCHP system of the traditional algorithm does not take the satisfaction rate of user requirements and the system operation cost as an optimization objective function, the invention provides a combined cooling heating and power optimization method of the CCHP system of the heat pump containing the sewage source, which can effectively solve the problems pointed out in the background technology.
The optimization method for combined cooling, heating and power of the CCHP system with the sewage source heat pump is characterized by comprising the following steps of:
s1, inputting original data of a combined cooling heating and power system, wherein the input of predicted cooling load, heat load, electric load, fan and photovoltaic cell output is included;
s2, constructing a CCHP micro-grid model, which comprises a gas turbine output power and heat production model, a waste heat boiler model, a fan and photovoltaic cell output model, an energy storage cell model, a gas boiler model, a sewage source heat pump model, an absorption refrigerator model and an electric refrigerator model;
s3, constructing a target function;
s4, constructing constraint conditions;
s5, solving an algorithm to obtain an optimization result;
wherein, the step S3 specifically comprises:
s3.1, taking the satisfaction rate of user requirements as an optimization objective function:
Figure BDA0002369228970000021
Figure BDA0002369228970000022
wherein, F 1.1 A rate of electrical load demand satisfaction; f 1.2 A thermal load satisfaction rate; f 1.3 A cold load satisfaction rate; f 1.4 Overall satisfaction rate for three types of users; mu, sigma, a,
Figure BDA0002369228970000023
Is a weight coefficient;
Figure BDA0002369228970000024
wherein, P MT.t The electric power provided by the micro gas turbine at the time t;
Figure BDA0002369228970000025
actual output of the wind turbine generator at the time t;
Figure BDA0002369228970000026
actual output of the photovoltaic cell at the moment t;
Figure BDA0002369228970000027
and
Figure BDA0002369228970000028
respectively the charging power and the discharging power at the time t; eta cha And η dis Respectively charge efficiency and discharge efficiency; p is EXB.t Purchasing power from the power distribution network for the combined supply system at the time t; p is EXS.t The power selling power from the combined supply system to the power distribution network; p EC.t The power consumption of the electric refrigerator in the t period is achieved; p is SE.t The power consumption of the sewage source heat pump in the time period t is determined;
Figure BDA0002369228970000029
an electrical load for a period of t;
since electricity purchase and electricity sale are not carried out simultaneously, the following regulations are provided:
P EXB.t P EXS.t =0
Figure BDA00023692289700000210
wherein H MT.t The heat production of the micro gas turbine is in a period t; h GF.t The heat production quantity of the gas boiler in the time period t is obtained; h SE.t The heat extraction quantity of the sewage source heat pump is t time period;
Figure BDA00023692289700000211
a thermal load for a period of t; alpha is alpha t The heat distribution coefficient of the waste heat boiler at the moment t is obtained; beta is a beta t The heat distribution coefficient of the sewage source heat pump at the moment t; eta WT The recovery efficiency of the waste heat boiler is obtained;
Figure BDA00023692289700000212
wherein Q is AC.t And Q EC.t Refrigerating capacities of the electric refrigerator and the absorption refrigerator in a t period respectively;
Figure BDA00023692289700000213
the time period t is the cold load;
s3.2, taking the running cost as an optimization objective function:
F 2 =C GAS +C EX +C PR +C OS
wherein, F 2 For operating costs; c GAS 、C EX 、C PR And C OS Respectively the system fuel cost, the difference between the system electricity purchasing cost and the electricity selling income, the system equipment operating cost and the equipment starting and stopping cost;
the fuel cost of the system:
Figure BDA00023692289700000214
wherein epsilon t Is the natural gas price; g MT.t The gas consumption of the gas turbine in the time period t; g GF.t The gas consumption of the gas boiler in the time period t;
the difference between the electricity purchasing cost and the electricity selling income of the system is as follows:
Figure BDA0002369228970000031
wherein, tau B.t And τ S.t The prices of electricity purchase and electricity sale from the combined cooling, heating and power supply system to the power distribution network are respectively;
electricity rate structure:
Figure BDA0002369228970000032
Figure BDA0002369228970000033
equipment operating cost:
Figure BDA0002369228970000034
Figure BDA0002369228970000035
wherein K is MT 、K GF 、K WT 、K PV 、K EC 、K SE The unit power running costs of a gas turbine, a gas boiler, a wind turbine generator, a photovoltaic cell, an electric refrigerator, a sewage source heat pump and an absorption refrigerator are respectively set; q AC.t The refrigerating power of the absorption refrigerator in the t period;
equipment start-stop cost:
Figure BDA0002369228970000036
Figure BDA0002369228970000037
Figure BDA0002369228970000038
Figure BDA0002369228970000039
Figure BDA00023692289700000310
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA00023692289700000311
the coefficient of cut for the gas generator;
Figure BDA00023692289700000312
the rated power of the gas generator; c OP.MT For gas turbine startup costs; c ST.MT Cost for gas turbine shutdown; c OP.GF The startup cost for the gas boiler; c ST.GF The shutdown cost of the gas boiler; SS MTO.t 、SS MTS.t 、SS GFO.t 、SS GFS.t For start-stop coefficient, SS MTO.t At 1, represents the micro gas turbine is turned on during the period t, SS MTS.t 1 represents the micro gas turbine is shut down during time t; SS GFO.t Is 1, the gas boiler is started in the time period t, SS GFS.t 1 represents the micro gas turbine is shut down during time t;
s3.3, synthesizing a satisfaction rate objective function and an operation cost objective function of user requirements to obtain a total objective function:
Figure BDA0002369228970000041
wherein the content of the first and second substances,
Figure BDA0002369228970000042
and δ are weight coefficients, respectively.
The specific method of the step S4 is as follows:
s4.1 energy balance constraints
Electric balance constraint:
the electric energy is mainly provided by gas turbine, distribution network, photovoltaic cell, wind turbine generator system, and the energy storage battery plays the undulant effect of exerting oneself of stabilizing wind turbine generator system and photovoltaic cell here:
Figure BDA0002369228970000043
and (3) thermal balance constraint:
the heat energy is mainly provided by the extracted heat energy generated by the sewage source heat pump and the waste heat boiler and the heat energy of the gas boiler:
Figure BDA0002369228970000044
cold balance constraint:
the cooling load demand is mainly provided by absorption chillers and electric chillers:
Figure BDA0002369228970000045
wherein, the heat of absorption formula refrigerator is provided by exhaust-heat boiler and sewage source heat pump:
Q EC.t =(H MT.t (1-α tWT +H SE.t (1-β t ))η EC
in the formula eta EC The refrigeration efficiency of the absorption refrigerator;
s4.2, restraining the operation working conditions of all devices:
the power generation power constraint of the micro gas turbine is as follows:
Figure BDA0002369228970000046
because the micro gas generator loses the natural gas to generate electricity and has certain constraint on the gas consumption, the climbing constraint exists on the corresponding generated power:
|P MT.t -P MT.t-1 |≤a P P MT.t
wherein, a P The power generation power climbing capacity coefficient of the micro gas turbine;
and (3) output force restraint of the gas boiler:
Figure BDA0002369228970000047
wherein the content of the first and second substances,
Figure BDA0002369228970000048
rated power for the gas boiler;
output restraint of the wind turbine generator:
Figure BDA0002369228970000049
wherein the content of the first and second substances,
Figure BDA00023692289700000410
the minimum cut coefficient of the fan is set;
Figure BDA00023692289700000411
rated power for the fan;
photovoltaic cell output restraint:
Figure BDA0002369228970000051
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0002369228970000052
is the rated power of the photovoltaic cell;
energy storage equipment restraint:
considering the problem of service life of the energy storage battery, the lowest storage capacity of the energy storage battery is limited, and the maximum charge capacity and the maximum discharge capacity of the energy storage battery are restricted in each time period, and the charging and the discharging cannot be simultaneously carried out in each time period, and the energy storage device is restricted as follows:
S BAT.min <S BAT.t <S BAT.max
Figure BDA0002369228970000053
Figure BDA0002369228970000054
Figure BDA0002369228970000055
wherein S is BAT.t The storage capacity of the energy storage battery at the moment t; s BAT.min The lowest allowable amount of stored electricity for the energy storage device; s BAT.max The rated power of the energy storage device; y is cha And y dis The maximum discharge rate and the minimum charge rate of the energy storage battery are respectively;
absorption chiller output restraint:
Figure BDA0002369228970000056
wherein the content of the first and second substances,
Figure BDA0002369228970000057
is the rated power of the absorption chiller;
electric refrigerator output restraint:
Figure BDA0002369228970000058
wherein the content of the first and second substances,
Figure BDA0002369228970000059
is the rated power of the electrical refrigerator.
The specific method of the step S5 is as follows:
s5.1: setting parameters such as population scale, iteration times and the like; setting a target function, and randomly generating particle positions; training a BP neural network according to original data;
s5.2: selecting a gas turbine to run in a mode of deciding electricity by heat or electricity by cold according to the water quantity and the water temperature of the sewage source;
s5.3: constructing a CCHP system model, an objective function and constraint conditions;
s5.4: carrying out transverse crossing on two non-repeated random pairings of particles in a population;
s5.5: comparing the fitness between the mediocre solution obtained after transverse crossing and the parent particles, and selecting the particles with high fitness as members of a new population;
s5.6: longitudinally crossing pairwise non-repeated random pairing between variables in each particle to obtain progeny particles;
s5.7: comparing the fitness between parent particles and child particles, and keeping the particles with high fitness;
s5.8: the iteration loop reaches a specified number of times or the precision reaches a convergence condition.
The invention has the beneficial effects that: the waste gas of the gas turbine is utilized to utilize waste heat through the waste heat boiler, meanwhile, the heat of urban sewage, industrial sewage and the like is extracted through the sewage source heat pump, the heat is converted into heat energy with higher quality through consuming a small amount of electricity, the heat energy and the heat obtained by extracting the waste heat in the waste gas of the gas turbine through the waste heat boiler are combined, meanwhile, start-stop functions of the gas turbine and the waste heat boiler are built, the start-stop cost of the two devices is calculated better, the satisfaction rate of user demands is considered, the satisfaction rate of the user demands and the system operation cost are taken as optimization objective functions, and the operation scheduling of the combined cooling, heating and power supply micro-grid system is optimized.
Drawings
Fig. 1 is a flowchart of an optimal scheduling method in a CCHP micro-grid according to the present invention.
Detailed Description
The invention will be further described with reference to the accompanying drawings.
1. And inputting the original data of the combined cooling heating and power system. Inputting the predicted cold load, heat load, electric load, fan and photovoltaic cell output, electricity rate structure and initial values of the operation of each device of the CCHP system.
And then constructing a CCHP microgrid model. The CCHP microgrid model comprises a gas turbine output power and heat production model, a waste heat boiler model, a fan and photovoltaic cell output model, an energy storage cell model, a gas boiler model, a sewage source heat pump model, an absorption refrigerator model and an electric refrigerator model.
2. Constructing an objective function:
201. taking the satisfaction rate of user requirements as an optimization objective function:
Figure BDA0002369228970000061
Figure BDA0002369228970000062
wherein, F 1.1 A rate of electrical load demand satisfaction; f 1.2 A thermal load satisfaction rate; f 1.3 A cold load satisfaction rate; f 1.4 Overall satisfaction rate for three types of users; mu, sigma,
Figure BDA0002369228970000063
Are weight coefficients.
Figure BDA0002369228970000064
Wherein, P MT.t The electric power provided by the micro gas turbine at the time t;
Figure BDA0002369228970000065
actual output of the wind turbine generator at the time t;
Figure BDA0002369228970000066
the actual output of the photovoltaic cell at the time t;
Figure BDA0002369228970000067
and
Figure BDA0002369228970000068
respectively the charging power and the discharging power at the time t; eta cha And η dis Respectively charge efficiency and discharge efficiency; p EXB.t Purchasing power from the power distribution network for the combined supply system at the time t; p EXS.t The power selling power from the combined supply system to the power distribution network; p EC.t The power consumption of the electric refrigerator in the t period is achieved; p SE.t The power consumption of the sewage source heat pump in the time period t is determined;
Figure BDA0002369228970000069
is the electrical load for the period t.
Since electricity purchase and electricity sale are not carried out simultaneously, the following regulations are provided:
P EXB.t P EXS.t =0
Figure BDA00023692289700000610
wherein H MT.t The heat production quantity of the micro gas turbine is t time period; h GF.t The heat production quantity of the gas boiler in the time period t; h SE.t The heat extraction quantity of the sewage source heat pump is t time period;
Figure BDA0002369228970000071
a thermal load for a period of t; alpha is alpha t The heat distribution coefficient of the waste heat boiler at the moment t is obtained; beta is a t Is the heat of the sewage source heat pump at the moment tA distribution coefficient; eta WT The recovery efficiency of the waste heat boiler is improved.
Figure BDA0002369228970000072
Wherein Q is AC.t And Q EC.t Refrigerating capacities of the electric refrigerator and the absorption refrigerator in a t period respectively;
Figure BDA0002369228970000073
is the cooling load for the period t.
202. Taking the running cost as an optimization objective function:
F 2 =C GAS +C EX +C PR +C OS
wherein, F 2 The cost of operation; c GAS 、C EX 、C PR And C OS Respectively the system fuel cost, the difference between the system electricity purchasing cost and the electricity selling profit, the system equipment operation cost and the equipment start-stop cost.
The fuel cost of the system:
Figure BDA0002369228970000074
wherein epsilon t Is the natural gas price; g MT.t The gas consumption of the gas turbine in the time period t; g GF.t Is the gas consumption of the gas boiler in the time period t.
The difference between the electricity purchasing cost and the electricity selling income of the system is as follows:
Figure BDA0002369228970000075
wherein, tau B.t And τ S.t The prices of electricity purchase and electricity sale from the combined cooling, heating and power supply system to the power distribution network are respectively;
electricity rate structure:
Figure BDA0002369228970000076
Figure BDA0002369228970000077
equipment operating cost:
Figure BDA0002369228970000078
Figure BDA0002369228970000079
wherein K is MT 、K GF 、K WT 、K PV 、K EC 、K SE The unit power operating costs of a gas turbine, a gas boiler, a wind turbine generator, a photovoltaic cell, an electric refrigerator, a sewage source heat pump and an absorption refrigerator are respectively; q AC.t The refrigerating power of the absorption refrigerator in the t period is shown.
Equipment start-stop cost:
Figure BDA00023692289700000710
Figure BDA0002369228970000081
Figure BDA0002369228970000082
Figure BDA0002369228970000083
Figure BDA0002369228970000084
wherein the content of the first and second substances,
Figure BDA0002369228970000085
the coefficient of cut for the gas generator;
Figure BDA0002369228970000086
the rated power of the gas generator; c OP.MT For gas turbine startup costs; c ST.MT Cost for gas turbine shutdown; c OP.GF The startup cost for the gas boiler; c ST.GF The shutdown cost of the gas boiler; SS MTO.t 、SS MTS.t 、SS GFO.t 、SS GFS.t For start-stop coefficient, SS MTO.t At 1, the micro gas turbine is turned on during the period t, SS MTS.t 1 represents the micro gas turbine is shut down during time t; SS GFO.t Is 1, the gas boiler is started in the time period t, SS GFS.t A time of 1 represents the micro gas turbine shutdown during time t.
203. And synthesizing a satisfaction rate objective function and an operation cost objective function of the user requirements to obtain a total objective function:
Figure BDA0002369228970000087
wherein the content of the first and second substances,
Figure BDA0002369228970000088
and δ are weight coefficients, respectively.
3. And (3) constructing a constraint condition:
301. energy balance constraint
Electric balance constraint:
the electric energy is mainly provided by gas turbine, distribution network, photovoltaic cell, wind turbine generator system, and the energy storage battery plays the undulant effect of exerting oneself of stabilizing wind turbine generator system and photovoltaic cell here:
Figure BDA0002369228970000089
and (3) thermal balance constraint:
the heat energy is mainly provided by the extracted heat energy generated by the sewage source heat pump and the waste heat boiler and the heat energy of the gas boiler:
Figure BDA00023692289700000810
cold balance constraint:
the cooling load demand is mainly provided by absorption chillers and electric chillers:
Figure BDA00023692289700000811
wherein, the heat of absorption formula refrigerator is provided by exhaust-heat boiler and sewage source heat pump:
Q EC.t =(H MT.t (1-α tWT +H SE.t (1-β t ))η EC
in the formula eta EC The refrigeration efficiency of the absorption refrigerator.
302. And (3) constraint of the operation condition of each device:
and (3) power generation power constraint of the micro gas turbine:
Figure BDA0002369228970000091
because the micro gas generator loses the natural gas to generate electricity and has certain constraint on the gas consumption, the climbing constraint exists on the corresponding generated power:
|P MT.t -P MT.t-1 |≤a P P MT.t
wherein, a P The power generation power climbing capacity coefficient of the micro gas turbine.
Output force restraint of the gas boiler:
Figure BDA0002369228970000092
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0002369228970000093
is the rated power of the gas boiler.
Output restraint of the wind turbine generator:
Figure BDA0002369228970000094
wherein the content of the first and second substances,
Figure BDA0002369228970000095
the minimum cut coefficient of the fan is set;
Figure BDA0002369228970000096
the rated power of the fan.
Photovoltaic cell output constraint:
Figure BDA0002369228970000097
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0002369228970000098
the rated power of the photovoltaic cell.
Energy storage equipment restraint:
considering the service life problem of the energy storage battery, the lowest storage capacity of the energy storage battery is limited, the maximum charging capacity and the maximum discharging capacity of the energy storage battery are restricted in each time period, the charging and the discharging cannot be simultaneously carried out in each time period, and the energy storage device is restricted as follows:
S BAT.min <S BAT.t <S BAT.max
Figure BDA0002369228970000099
Figure BDA00023692289700000910
Figure BDA00023692289700000911
wherein S is BAT.t The storage capacity of the energy storage battery at the moment t; s BAT.min The lowest allowable charge capacity for the energy storage device; s. the BAT.max The rated power of the energy storage device; y is cha And y dis Respectively, the maximum discharge rate and the minimum charge rate of the energy storage battery.
Absorption chiller output restraint:
Figure BDA00023692289700000912
wherein the content of the first and second substances,
Figure BDA00023692289700000913
is the rated power of the absorption chiller.
Electric refrigerator output restraint:
Figure BDA00023692289700000914
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0002369228970000101
is the rated power of the electrical refrigerator.
4. The solving method comprises the following steps:
the traditional group optimization algorithm has the advantages of high convergence rate and the like, so that the traditional group optimization algorithm is widely applied at present. But all have the phenomenon of "precocity" by using the traditional group optimization algorithm. The vertical and horizontal crossing algorithm has the advantages of high convergence speed and the like of the traditional group optimization algorithm, and the capability of jumping out of local optimum of the group is realized by utilizing the vertical crossing, so that the premature phenomenon is avoided. Therefore, an improved criss-cross algorithm is used as a combined cooling, heating and power optimization method for a CCHP system of a heat pump containing a sewage source, and the updating formula of the longitudinal cross position is as follows:
MS VC (i,d 1 )=rX(i,d 1 )+(1-r)X(i,d 2 )
i∈N(1,M),d 1 ,d 2 ∈N(1,D)
wherein r is [0,1 ]]Random numbers uniformly distributed thereon; MS (Mass Spectrometry) VC (i,d 1 ) Is particle i d 1 D a sum 2 Progeny resulting from the cross of dimensions.
The internal parts of the same particles are longitudinally crossed at a certain probability, so that the global search capability of the population can be improved, but the defect of low efficiency exists, and each crossing needs to be normalized and reproduced, so that the update speed of the population is influenced. On the basis, X (i, d) in the formula of longitudinal intersection is used 2 ) The improved CSO algorithm finds out potential linear relation between the longitudinal crossing scalars through a BP neural network, so that the longitudinal crossing process is more purposeful, and replaces two steps of normalization and recurrence. The improvement method comprises the following steps:
a. firstly, according to the operation experience, the original data of the optimal operation state is obtained.
b. And training N BP neural networks by using the obtained original data. Each neural network corresponds to two different variables of the same particle, and when the input of one neural network in the N BP neural networks is d 1 And X (i, d) 2 ) When output is as
Figure BDA0002369228970000102
Wherein X avg (i,d 1 ) Denotes the d-th 1 The average of the dimensional variables;
Figure BDA0002369228970000103
n is the number of variables.
c. Mixing X (i, d) 2 ) Is changed into
Figure BDA0002369228970000104
The significance of the method is that the cross direction is close to the optimal running state, so that the method has better global searching capability.
The method comprises the following steps:
s5.1: and setting parameters such as population scale and iteration times. And setting an objective function and randomly generating particle positions. And training the BP neural network according to the original data.
S5.2: and selecting the operation mode of the gas turbine by using heat to fix the power or using cold to fix the power according to the water quantity and the water temperature of the sewage source.
S5.3: and constructing a CCHP system model, an objective function and constraint conditions.
S5.4: and transversely crossing two non-repeated random pairings of the particles in the population.
S5.5: and comparing the fitness between the mediocre solution obtained after transverse crossing and the parent particles, and selecting the particles with high fitness as members of the new population.
S5.6: and longitudinally crossing pairwise non-repeated random pairing between variables in each particle to obtain progeny particles after crossing.
S5.7: and comparing the fitness between the parent particles and the child particles, and keeping the particles with high fitness.
S5.8: the iteration loop reaches a specified number of times or the precision reaches a convergence condition.
The invention constructs a CCHP system model containing a gas turbine output power and heat production model, a waste heat boiler model, a fan and photovoltaic cell output model, an energy storage cell model, a gas boiler model, a sewage source heat pump model, an absorption refrigerator model and an electric refrigerator model. The waste heat of the waste gas of the gas turbine is utilized by the waste heat boiler, meanwhile, the heat of urban sewage, industrial sewage and the like is extracted by the sewage source heat pump, the heat is converted into heat energy with higher quality by consuming a small amount of electricity, the heat energy is combined with the heat obtained by extracting the waste heat in the waste gas of the gas turbine by the waste heat boiler, when the water quantity of the sewage is large or the temperature is higher, the heat energy supplied by the heat pump is higher, the gas turbine adopts a mode of 'fixing heat by electricity', otherwise, when the water quantity of the sewage is lower or the temperature is lower, the gas turbine adopts a mode of 'fixing electricity by heat', and supplies the heat energy to the CCHP system. The method constructs the start-stop functions of the gas turbine and the waste heat boiler at the same time, calculates the start-stop cost of the two devices better, considers the satisfaction rate of user requirements, and optimizes the operation scheduling of the cooling, heating and power triple supply type micro-grid system by taking the satisfaction rate of the user requirements and the system operation cost as optimization objective functions.
Moreover, aiming at the problems that the traditional algorithm is easy to fall into local optimum, the global search capability is poor, the convergence speed is low and the like, the improved vertical and horizontal crossing algorithm which utilizes the BP neural network to replace irrelevant variables in the vertical crossing with potential relation operators is provided, and the capability of the vertical crossing to help a group to jump out of the local optimum is greatly improved.
The above description is only an embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art can easily change or replace the present invention within the technical scope of the present invention. Therefore, the protection scope of the present invention is subject to the protection scope of the claims.

Claims (4)

1. A combined cooling heating and power optimization method for a CCHP system with a sewage source heat pump is characterized by comprising the following steps:
s1, inputting original data of a combined cooling heating and power system, wherein the input of predicted cooling load, heat load, electric load, fan and photovoltaic cell output is included;
s2, constructing a CCHP micro-grid model, which comprises a gas turbine output power and heat production model, a waste heat boiler model, a fan and photovoltaic cell output model, an energy storage cell model, a gas boiler model, a sewage source heat pump model, an absorption refrigerator model and an electric refrigerator model;
s3, constructing a target function;
s4, constructing constraint conditions;
s5, solving the CCHP micro-grid model in the step S2 to obtain an optimization result;
wherein, the step S3 specifically comprises the following steps:
s3.1, taking the satisfaction rate of user requirements as an optimization objective function:
Figure FDA0003715432930000011
Figure FDA0003715432930000012
wherein, F 1.1 A rate of electrical load demand satisfaction; f 1.2 Heat load satisfaction rate; f 1.3 A cold load satisfaction rate; f 1.4 Overall satisfaction rate for three types of users; mu, sigma,
Figure FDA0003715432930000013
Is a weight coefficient;
Figure FDA0003715432930000014
wherein, P MT.t The electric power provided by the micro gas turbine at the time t;
Figure FDA0003715432930000015
actual output of the wind turbine generator at the time t;
Figure FDA0003715432930000016
the actual output of the photovoltaic cell at the time t;
Figure FDA0003715432930000017
and
Figure FDA0003715432930000018
respectively the charging power and the discharging power at the time t; eta cha And η dis Respectively, charging efficiency and discharging efficiency; p EXB.t Purchasing power from the power distribution network for the combined supply system at the time t; p EXS.t The power is the power sold by the combined supply system to the distribution network; p EC.t The electric power of the electric refrigerator at the time t is used; p is SE.t The power consumption of the sewage source heat pump at the time t is obtained;
Figure FDA0003715432930000019
is the electrical load at time t;
since electricity purchase and electricity sale are not carried out simultaneously, the following regulations are provided:
P EXB.t P EXS.t =0
Figure FDA00037154329300000110
wherein H MT.t The heat production capacity of the micro gas turbine at the moment t; h GF.t The heat production quantity of the gas boiler at the time t is obtained; h SE.t The heat extraction quantity of the sewage source heat pump at the moment t;
Figure FDA0003715432930000021
thermal load for time t; alpha is alpha t The heat distribution coefficient of the waste heat boiler at the moment t is obtained; beta is a t The heat distribution coefficient of the sewage source heat pump at the moment t is obtained; eta WT The recovery efficiency of the waste heat boiler is obtained;
Figure FDA0003715432930000022
wherein Q is AC.t And Q EC.t Refrigerating capacities of the electric refrigerator and the absorption refrigerator at the moment t respectively;
Figure FDA0003715432930000023
is the cold load at time t;
s3.2, taking the running cost as an optimization objective function:
F 2 =C GAS +C EX +C PR +C OS
wherein, F 2 The cost of operation; c GAS 、C EX 、C PR And C OS Respectively the system fuel cost, the difference between the system electricity purchasing cost and the electricity selling profit, the system equipment operation cost and the equipment start-stop cost;
the fuel cost of the system:
Figure FDA0003715432930000024
wherein epsilon t Is the natural gas price; g MT.t The gas consumption of the gas turbine at the time t is obtained; g GF.t The gas consumption of the gas boiler at the time t is measured;
the difference between the electricity purchasing cost and the electricity selling income of the system is as follows:
Figure FDA0003715432930000025
wherein, tau B.t And τ S.t The prices of electricity purchasing and electricity selling from the combined cooling, heating and power system to the power distribution network are respectively;
structure of electric rate:
Figure FDA0003715432930000026
Figure FDA0003715432930000027
equipment operating cost:
Figure FDA0003715432930000028
wherein K MT 、K GF 、K WT 、K PV 、K EC 、K SE The unit power running costs of a gas turbine, a gas boiler, a wind turbine generator, a photovoltaic cell, an electric refrigerator, a sewage source heat pump and an absorption refrigerator are respectively set; q AC.t The refrigeration power of the absorption refrigerator at the moment t;
equipment start-stop cost:
Figure FDA0003715432930000031
Figure FDA0003715432930000032
Figure FDA0003715432930000033
Figure FDA0003715432930000034
Figure FDA0003715432930000035
wherein, the first and the second end of the pipe are connected with each other,
Figure FDA0003715432930000036
the coefficient of cut for the gas generator;
Figure FDA0003715432930000037
the rated power of the gas generator; c OP.MT Startup costs for the gas turbine; c ST.MT Cost for gas turbine shutdown; c OP.GF The startup cost for the gas boiler; c ST.GF The shutdown cost of the gas boiler; SS MTO.t 、SS MTS.t 、SS GFO.t 、SS GFS.t For start-stop coefficient, SS MTO.t Time 1 represents the micro gas turbine being turned on at time t, SS MTS.t A time of 1 represents that the micro gas turbine is shut down at time t; SS GFO.t Is 1 time represents that the gas boiler is started at the time t, SS GFS.t 1 represents the shutdown of the micro gas turbine at time t;
s3.3, synthesizing a satisfaction rate objective function and an operation cost objective function of the user requirements to obtain a total objective function:
Figure FDA0003715432930000038
wherein the content of the first and second substances,
Figure FDA0003715432930000039
and δ are weight coefficients, respectively.
2. The optimization method for combined cooling, heating and power of the CCHP system of the sewage-source-containing heat pump according to claim 1, wherein the specific method of S4 is as follows:
s4.1 energy balance constraints
Electric balance constraint:
the electric energy is mainly provided by gas turbine, distribution network, photovoltaic cell, wind turbine generator system, and the energy storage battery plays the undulant effect of exerting oneself of stabilizing wind turbine generator system and photovoltaic cell here:
Figure FDA00037154329300000310
and (4) thermal balance constraint:
the heat energy is mainly provided by the extracted heat energy generated by the sewage source heat pump and the waste heat boiler and the heat energy of the gas boiler:
Figure FDA0003715432930000041
cold balance constraint:
the cooling load demand is mainly provided by absorption chillers and electric chillers:
Figure FDA0003715432930000042
wherein, the heat of absorption formula refrigerator is provided by exhaust-heat boiler and sewage source heat pump:
Q EC.t =(H MT.t (1-α tWT +H SE.t (1-β t ))η EC
in the formula eta EC The refrigeration efficiency of the absorption refrigerator;
s4.2, restraining the operation condition of each device:
and (3) power generation power constraint of the micro gas turbine:
Figure FDA0003715432930000043
because the micro gas generator loses the natural gas and generates electricity, certain constraint exists on the gas consumption amount of the natural gas, and therefore, climbing constraint exists on the corresponding generated power:
|P MT.t -P MT.t-1 |≤a P P MT.t
wherein, a P The power generation power climbing capacity coefficient of the micro gas turbine;
and (3) output force restraint of the gas boiler:
Figure FDA0003715432930000044
wherein, the first and the second end of the pipe are connected with each other,
Figure FDA0003715432930000045
rated power for the gas boiler;
output restraint of the wind turbine generator:
Figure FDA0003715432930000046
wherein the content of the first and second substances,
Figure FDA0003715432930000047
the minimum cutting coefficient of the fan is set;
Figure FDA0003715432930000048
rated power for the fan;
photovoltaic cell output constraint:
Figure FDA0003715432930000049
wherein, the first and the second end of the pipe are connected with each other,
Figure FDA00037154329300000410
is the rated power of the photovoltaic cell;
energy storage equipment restraint:
considering the service life problem of the energy storage battery, the lowest storage capacity of the energy storage battery is limited, the maximum charging capacity and the maximum discharging capacity of the energy storage battery are restricted in each time period, the charging and the discharging cannot be simultaneously carried out in each time period, and the energy storage device is restricted as follows:
S BAT,min <S BAT.t <S BAT.max
Figure FDA0003715432930000051
Figure FDA0003715432930000052
Figure FDA0003715432930000053
wherein S is BAT.t The storage capacity of the energy storage battery at the moment t; s. the BAT.min The lowest allowable amount of stored electricity for the energy storage device; s BAT.max The rated power of the energy storage device; y is cha And y dis Respectively the maximum discharge rate and the minimum charge rate of the energy storage battery;
absorption chiller output restraint:
Figure FDA0003715432930000054
wherein the content of the first and second substances,
Figure FDA0003715432930000055
is the rated power of the absorption chiller;
electric refrigerator output restraint:
Figure FDA0003715432930000056
wherein, the first and the second end of the pipe are connected with each other,
Figure FDA0003715432930000057
is the rated power of the electrical refrigerator.
3. The CCHP system cooling, heating and power cogeneration optimization method of a sewage-source heat pump as claimed in claim 1,
the specific method of S5 is as follows:
s5.1: setting the population scale, the iteration times and the objective function, and randomly generating the particle positions; training a BP neural network according to original data;
s5.2: selecting a gas turbine to run in a mode of fixing power by heat or fixing power by cold according to the water quantity and the water temperature of a sewage source;
s5.3: constructing a CCHP microgrid model, a target function and constraint conditions;
s5.4: carrying out transverse crossing on two non-repeated random pairings of particles in a population;
s5.5: comparing the fitness between the mediocre solution obtained after transverse crossing and the parent particles, and selecting the particles with the maximum fitness as members of the new population;
s5.6: longitudinally crossing pairwise non-repeated random pairing between variables in each particle to obtain progeny particles;
s5.7: comparing the fitness between parent particles and child particles, and keeping the particles with the maximum fitness;
s5.8: the iteration loop reaches a specified number of times or the precision reaches a convergence condition.
4. The CCHP system cooling, heating and power combined optimization method of the sewage-source heat pump according to claim 3, wherein for the S5.6 step, the following steps are adopted to improve the global search capacity:
a. firstly, obtaining original data of an optimal running state according to running experience;
b. training N BP neural networks by using the obtained original data; each neural network corresponds to two different variables of the same particle, and when the input of one neural network in the N BP neural networks is d 1 And X (i, d) 2 ) When output is
Figure FDA0003715432930000061
Wherein X avg (i,d 1 ) Denotes the d-th 1 The mean of the dimensional variables;
Figure FDA0003715432930000062
n is the number of optimization variables contained in each particle;
c. mixing X (i, d) 2 ) Instead of using
Figure FDA0003715432930000063
The meaning is to make the cross direction approach to the optimum operation state.
CN202010045449.5A 2020-01-16 2020-01-16 CCHP system cooling, heating and power combined optimization method containing sewage source heat pump Expired - Fee Related CN111275251B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010045449.5A CN111275251B (en) 2020-01-16 2020-01-16 CCHP system cooling, heating and power combined optimization method containing sewage source heat pump

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010045449.5A CN111275251B (en) 2020-01-16 2020-01-16 CCHP system cooling, heating and power combined optimization method containing sewage source heat pump

Publications (2)

Publication Number Publication Date
CN111275251A CN111275251A (en) 2020-06-12
CN111275251B true CN111275251B (en) 2022-10-04

Family

ID=71001641

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010045449.5A Expired - Fee Related CN111275251B (en) 2020-01-16 2020-01-16 CCHP system cooling, heating and power combined optimization method containing sewage source heat pump

Country Status (1)

Country Link
CN (1) CN111275251B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112257899A (en) * 2020-09-22 2021-01-22 国网河北省电力有限公司营销服务中心 CCHP system optimal scheduling method and terminal equipment
CN113239550B (en) * 2021-05-18 2022-03-04 华北电力大学 Energy storage capacity configuration method for combined cooling heating and power system and combined cooling and power system
CN114626616B (en) * 2022-03-21 2023-11-21 特斯联科技集团有限公司 Optimal operation method and system for recycling CCHP (CCHP) from high-temperature tail gas waste heat in park
CN116736780B (en) * 2023-08-15 2023-11-24 贵州汇通华城股份有限公司 Startup and shutdown control optimization method and system for regional energy station

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104392286A (en) * 2014-12-02 2015-03-04 山东大学 Microgrid operation optimizing method by considering combined supply of cooling, heating and power with stored energy operation strategy
CN105631557A (en) * 2016-03-17 2016-06-01 北京睿新科技有限公司 Micro-grid optimal scheduling method considering ice-storage air conditioning and having power and cooling coupled characteristic
CN109165788A (en) * 2018-09-04 2019-01-08 河北工业大学 A kind of optimization method of cooling heating and power generation system
CN109245093A (en) * 2018-09-28 2019-01-18 上海电力学院 A kind of supply of cooling, heating and electrical powers distributed busbar protection collaboration Optimization Scheduling
CN110598915A (en) * 2019-08-23 2019-12-20 广东工业大学 Economic optimization method for micro-grid for combined heat and power supply

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104392286A (en) * 2014-12-02 2015-03-04 山东大学 Microgrid operation optimizing method by considering combined supply of cooling, heating and power with stored energy operation strategy
CN105631557A (en) * 2016-03-17 2016-06-01 北京睿新科技有限公司 Micro-grid optimal scheduling method considering ice-storage air conditioning and having power and cooling coupled characteristic
CN109165788A (en) * 2018-09-04 2019-01-08 河北工业大学 A kind of optimization method of cooling heating and power generation system
CN109245093A (en) * 2018-09-28 2019-01-18 上海电力学院 A kind of supply of cooling, heating and electrical powers distributed busbar protection collaboration Optimization Scheduling
CN110598915A (en) * 2019-08-23 2019-12-20 广东工业大学 Economic optimization method for micro-grid for combined heat and power supply

Also Published As

Publication number Publication date
CN111275251A (en) 2020-06-12

Similar Documents

Publication Publication Date Title
CN111275251B (en) CCHP system cooling, heating and power combined optimization method containing sewage source heat pump
Li et al. A hybrid optimization-based scheduling strategy for combined cooling, heating, and power system with thermal energy storage
CN109165788B (en) Optimization method of combined cooling heating and power system
Chen et al. An optimization on an integrated energy system of combined heat and power, carbon capture system and power to gas by considering flexible load
CN105048516B (en) A kind of honourable extreme misery multi-source complementation Optimization Scheduling
Jalili et al. Stochastic optimal operation of a microgrid based on energy hub including a solar-powered compressed air energy storage system and an ice storage conditioner
CN103580061B (en) Microgrid operating method
CN109146182A (en) The economic load dispatching method of meter and the distributed triple-generation system of a variety of energy storage
CN108009693A (en) Grid-connected micro-capacitance sensor dual blank-holder based on two-stage demand response
CN109474025B (en) Optimized dispatching model of park level comprehensive energy system
Liu et al. Co-optimization of a novel distributed energy system integrated with hybrid energy storage in different nearly zero energy community scenarios
CN112270433B (en) Micro-grid optimization method considering renewable energy uncertainty and user satisfaction
CN111737884B (en) Multi-target random planning method for micro-energy network containing multiple clean energy sources
CN111160636B (en) CCHP type micro-grid scheduling optimization method
CN112365021A (en) Regional comprehensive energy system planning optimization method based on hybrid energy storage
CN111668878A (en) Optimal configuration method and system for renewable micro-energy network
CN111859683A (en) Optimal configuration method of park comprehensive energy system based on dynamic energy concentrator
Kang et al. Feed-forward active operation optimization for CCHP system considering thermal load forecasting
CN113255198A (en) Multi-objective optimization method for combined cooling, heating and power supply micro-grid with virtual energy storage
CN117081143A (en) Method for promoting coordination and optimization operation of park comprehensive energy system for distributed photovoltaic on-site digestion
Yang et al. Optimal scheduling of CCHP with distributed energy resources based on water cycle algorithm
CN112883630A (en) Day-ahead optimized economic dispatching method for multi-microgrid system for wind power consumption
CN111126675A (en) Multi-energy complementary microgrid system optimization method
CN113344357B (en) Design method of comprehensive energy system based on frequency domain dynamic index
CN114091917A (en) Dynamic environment-friendly economic dispatching method and system for cold-heat-power combined supply type micro-grid

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20221004

CF01 Termination of patent right due to non-payment of annual fee