CN107919686A - A kind of wide area source lotus active optimization control method a few days ago - Google Patents

A kind of wide area source lotus active optimization control method a few days ago Download PDF

Info

Publication number
CN107919686A
CN107919686A CN201711036152.7A CN201711036152A CN107919686A CN 107919686 A CN107919686 A CN 107919686A CN 201711036152 A CN201711036152 A CN 201711036152A CN 107919686 A CN107919686 A CN 107919686A
Authority
CN
China
Prior art keywords
mrow
msubsup
load
power
few days
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201711036152.7A
Other languages
Chinese (zh)
Other versions
CN107919686B (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.)
North China Electric Power University
State Grid Gansu Electric Power Co Ltd
Electric Power Research Institute of State Grid Gansu Electric Power Co Ltd
Taiyuan Power Supply Co of State Grid Shanxi Electric Power Co Ltd
Original Assignee
North China Electric Power University
State Grid Gansu Electric Power Co Ltd
Electric Power Research Institute of State Grid Gansu Electric Power Co Ltd
Taiyuan Power Supply Co of State Grid Shanxi Electric Power Co Ltd
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 North China Electric Power University, State Grid Gansu Electric Power Co Ltd, Electric Power Research Institute of State Grid Gansu Electric Power Co Ltd, Taiyuan Power Supply Co of State Grid Shanxi Electric Power Co Ltd filed Critical North China Electric Power University
Priority to CN201711036152.7A priority Critical patent/CN107919686B/en
Publication of CN107919686A publication Critical patent/CN107919686A/en
Application granted granted Critical
Publication of CN107919686B publication Critical patent/CN107919686B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48Controlling the sharing of the in-phase component
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects
    • 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
    • Y02P80/00Climate change mitigation technologies for sector-wide applications
    • Y02P80/20Climate change mitigation technologies for sector-wide applications using renewable energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a kind of wide area source lotus active optimization control method a few days ago, including:Read the related data of all kinds of power supplys, centralized industrial load and distributed appliance load;Polymerization processing is grouped to distributed appliance load;Calculate the polymerization regulation performance parameter of each group aggregate load;Build and solve wide area source lotus active optimization model a few days ago, obtain all kinds of power supplys, centralized industrial load, the active power plan a few days ago for polymerizeing appliance load;The active power plan a few days ago of each distributed load in distribution polymerization appliance load.Active optimization control method has taken into full account the control characteristic of polymorphic elastic load to a kind of wide area source lotus provided by the invention a few days ago, solves the problems, such as the specificity analysis of magnanimity wide-area distribution type load using the mode of polymerization, realize the coordination optimization control of the polymorphic elastic load of wide area, normal power supplies and wind-powered electricity generation, improve digestion capability of the power grid to wind-powered electricity generation.

Description

A kind of wide area source-lotus active optimization control method a few days ago
Technical field
The invention belongs to New-energy power system active optimization control field, more particularly to a kind of wide area source-lotus to have a few days ago Work(optimal control method.
Background technology
With large-scale ten million multikilowatt wind power base installed capacity of wind-driven power sharp increase, the accounting of wind-powered electricity generation in systems increases Add.The characteristics of randomness and fluctuation of wind-powered electricity generation, peak load regulation network pressure is added, only rely on conventional power unit output and can not achieve wind The abundant consumption of electricity.Load side, centralized elastic load (such as high energy load) and distributed elastic load (such as thermal storage electric boiler Deng) good control characteristic is shown, dispatching of power netwoks can be participated in, promotes the consumption of wind-powered electricity generation.However, elastic load adjusts spy Property it is complicated, and distributed elastic load has the characteristics that distribution wide area, quantity magnanimity, adds the difficulty of dispatching of power netwoks.
Coordinate control aspect for source-lotus, both at home and abroad existing preliminary research.There is document to be made for large-scale wind power is grid-connected Into peaking problem, feasibility and validity that " lotus-source " coordinates control are participated in the centralized elastic load such as high energy load Studied.There is the theoretical research that document has carried out price type and stimulable type demand response participates in source-lotus coordination control, but simultaneously The adjustable characteristic of wide area elastic load is not deeply considered.Have a literature research thermal storage electric boiler distributed elastic load can Control characteristic, but source-lotus coordination control is not participated in it and is furtherd investigate.
In conclusion for the wide area elasticity such as centralized load (high energy load) and distributed load (thermal storage electric boiler) Load, it is common to participate in " lotus-source " optimal control aspect, still lack systematic Study at present.In order to enable wide area elastic load Preferably participate in power grid to adjust, fully dissolve wind-powered electricity generation, it is necessary to propose a kind of new wide area source-lotus active optimization controlling party a few days ago Method, realizes the active optimization control for considering wind-powered electricity generation, normal power supplies, wide area elastic load.
The content of the invention
The object of the present invention is to provide a kind of wide area source-lotus active optimization control method a few days ago, for solving extensive wind Under the grid-connected background of electricity, source-lotus active optimization control problem, reference is provided for operation of power networks.
A kind of wide area source-lotus active optimization control method a few days ago, comprises the following steps:
S1:Read the related data of all kinds of power supplys, centralized industrial load and distributed appliance load;
S2:Polymerization processing is grouped to distributed appliance load;
S3:Calculate the polymerization regulation performance parameter of each group aggregate load;
S4:Build and solve wide area source-lotus active optimization model a few days ago, obtain all kinds of power supplys, centralized industrial load, gather Close the active power plan a few days ago of appliance load;
S5:The active power plan a few days ago of each distributed load in distribution polymerization appliance load.
The S1 comprises the following steps:
S101:Obtain the quantity M of centralized industrial load in power gridH, the quantity N of distributed appliance load;
S102:Obtain next day wind-powered electricity generation active power output prediction dataConventional load power prediction dataEach collection Chinese style industrial load power prediction dataM=1 .., MH, each distributed appliance load power prediction dataM= 1,...,N;
S103:Obtain the performance parameter of normal power supplies:The peak power output P of conventional power unit iGi,max, minimum output power PGi,min, minimum run timeMinimum downtimeRise and power restriction PGi,up, decline and power restriction PGi,down
S104:Obtain centralized industrial load regulation performance parameter:The maximum load power of centralized industrial load m Minimum load powerMaximum emersion power restriction DPh +, minimum power restriction drops
S105:Obtain the regulation performance parameter of distributed appliance load:The minimum load power of distributed appliance load nMaximum load powerMinimum amount of stored heatMaximum amount of stored heatMaximum climbing rate DPx+ m,nClimbed with minimum Ratio of slope
The S2 comprises the following steps:
S201:N number of distributed appliance load is grouped, distributed civilian deferrable load only considers hot stored electric herein Boiler load;The thermal storage electric boiler load that model is identical, regulating power scope is identical is grouped into the same group, is divided into MxGroup, M group thermal storage electric boilers quantity is Nm
S202:To MxThe processing of distribution type thermal storage electric boiler Load aggregation, polymerization process such as attached drawing 1.
The S3 comprises the following steps:
S301:Calculate MxThe secondary daily load prediction power of group polymerization appliance load
S302:Calculate MxThe minimum load power of group polymerization appliance loadAnd maximum load power
S303:Calculate MxThe minimum amount of stored heat of group polymerization appliance loadWith maximum amount of stored heat
S304:Calculate MxThe maximum climbing rate DP of group polymerization appliance loadx+ mWith minimum climbing rate
The S4 comprises the following steps:
S401:Wide area source-lotus active optimization model a few days ago is built, sets target function is to abandon air quantity minimum, constraints Including power-balance constraint, wind power output constraint, normal power supplies operation constraint, aggregate load constraint, the constraint of high energy load.
S402:Solution obtains wind-powered electricity generation active power output plan a few days agoNormal power supplies active power output plan a few days agoConcentrate Formula industrial load day preload effective power meter drawIt polymerize appliance load day preload effective power meter to draw
The S5 comprises the following steps:
Calculate each distributed load active power plan a few days ago in polymerization appliance load
A kind of wide area source-lotus provided by the invention active optimization control method a few days ago, passes through:Read all kinds of power supplys, concentration The related data of formula industrial load and distributed appliance load;Polymerization processing is grouped to distributed appliance load;Calculate each The polymerization regulation performance parameter of group aggregate load;Build and solve wide area source-lotus active optimization model a few days ago, obtain all kinds of electricity Source, centralized industrial load, the active power plan a few days ago for polymerizeing appliance load;It is each distributed negative in distribution polymerization appliance load The active power plan a few days ago of lotus.The method has taken into full account that the adjusting of the grid-connected fluctuation of large-scale wind power and elastic load is special Property, solve the problems, such as distributed load quantity magnanimity by way of polymerization, can preferably play wide area elastic load, often Complementary adjustment effect of the power supply to wind-powered electricity generation is advised, promotes wind electricity digestion.
Brief description of the drawings
Below by drawings and examples, technical scheme is described in further detail.
Fig. 1 is a kind of wide area source-lotus provided by the invention active optimization control method flow chart a few days ago.
Fig. 2 is regional power grid schematic diagram in example 2 provided by the invention
Fig. 3 is that the next day active power of the wind-powered electricity generation of regional power grid and each type load predicts number in example 2 provided by the invention According to
Fig. 4 is the next day active power prediction data of the polymerization appliance load of regional power grid in example 2 provided by the invention
Fig. 5 is that the wind-powered electricity generation of regional power grid in example 2 provided by the invention, normal power supplies, high energy load, polymerization are civilian negative The active power plan a few days ago of lotus
Fig. 6 is the effective power meter a few days ago of each distributed appliance load of regional power grid in example 2 provided by the invention Draw
Embodiment
In order to have a clear understanding of technical scheme, its detailed structure will be set forth in the description that follows.Obviously, originally The specific simultaneously deficiency of implementing of inventive embodiments is limited to the specific details that those skilled in the art is familiar with.The typical case of the present invention is real Apply example to be described in detail as follows, in addition to these embodiments of detailed description, there can also be other embodiment.
The present invention is described in further details with reference to the accompanying drawings and examples.
Embodiment 1
Fig. 1 is a kind of flow chart of wide area source-lotus active optimization control method a few days ago.In Fig. 1, one kind provided by the invention Active optimization control method flow chart includes wide area source-lotus a few days ago:
S1:Read the related data of all kinds of power supplys, centralized industrial load and distributed appliance load;
S2:Polymerization processing is grouped to distributed appliance load;
S3:Calculate the polymerization regulation performance parameter of each group aggregate load;
S4:Build and solve wide area source-lotus active optimization model a few days ago, obtain all kinds of power supplys, centralized industrial load, gather Close the active power plan a few days ago of appliance load;
S5:The active power plan a few days ago of each distributed load in distribution polymerization appliance load.
The S1 comprises the following steps:
S101:Obtain the quantity M of centralized industrial load in power gridH, the quantity N of distributed appliance load;
S102:Obtain next day wind-powered electricity generation active power output prediction dataConventional load power prediction dataEach collection Chinese style industrial load power prediction dataM=1 .., MH, each distributed appliance load power prediction dataM= 1,...,N;
S103:Obtain the performance parameter of normal power supplies:The peak power output P of conventional power unit iGi,max, minimum output power PGi,min, minimum run timeMinimum downtimeRise and power restriction PGi,up, decline and power restriction PGi,down
S104:Obtain centralized industrial load regulation performance parameter:The maximum load power of centralized industrial load m Minimum load powerMaximum emersion power restriction DPh +, minimum power restriction drops
S105:Obtain the regulation performance parameter of distributed appliance load:The minimum load power of distributed appliance load nMaximum load powerMinimum amount of stored heatMaximum amount of stored heatMaximum climbing rate DPx+ m,n, minimum climbing Rate
The S2 comprises the following steps:
S201:N number of distributed appliance load is grouped, distributed civilian deferrable load only considers hot stored electric herein Boiler load;The identical thermal storage electric boiler load of thermal storage electric boiler model, amount of stored heat, rated power is grouped into the same group, is divided into For MxGroup, m group thermal storage electric boilers quantity are Nm
S202:To MxThe processing of distribution type thermal storage electric boiler Load aggregation, concretely comprises the following steps:M groups polymerize appliance load Power is equal to NmThe superposition of a thermal storage electric boiler load power, i.e.,M groups polymerization appliance load amount of stored heat etc. In NmThe superposition of a thermal storage electric boiler load amount of stored heat, i.e.,
The S3 comprises the following steps:
S301:Calculate Mx24 pre- power scales of point load of next day of group polymerization appliance load:M groups polymerization appliance load t The load prediction power of period is:
Wherein,The load prediction power of n-th of distributed appliance load t period
S302:Calculate MxThe minimum and maximum regulating power of group polymerization appliance load:
Wherein,For the minimum and maximum regulating power of n-th of distributed appliance load.
S303:Calculate MxThe minimum and maximum amount of stored heat of group polymerization appliance load:
Wherein,For the minimum and maximum amount of stored heat of n-th of distributed appliance load.
S304:Calculate MxThe maximum climbing rate and minimum climbing rate of group polymerization appliance load:
Wherein,For the minimum and maximum climbing rate of n-th of distributed appliance load.
The S4 comprises the following steps:
S401:Build wide area source-lotus active optimization model a few days ago
The essence of active optimization control problem is to meet that electric system is normal to source-lotus that polymorphic elastic load participates in a few days ago On the premise of stipulations beam, the adjustment effect of normal power supplies and elastic load is given full play to, more to dissolve large-scale wind power.Therefore To abandon the minimum target of air quantity, structure object function is as follows:
In formula,For t periods wind-powered electricity generation plan a few days ago power generating value;Predicted for t next day, active power for wind power period Value;Δ t is the duration of unit period, is set as 15 minutes;
Wherein constraints includes:
(1) power-balance constraint:
In formula, NGThe number of conventional power unit in expression system;Start and stop shapes of the conventional power unit i in period t in expression system State, is represented with 0 or 1;Represent active power outputs of the conventional power unit i in the t periods;For the pre- power scale of t period conventional loads;The unscheduled power a few days ago of respectively m-th of high energy load and polymerization appliance load in the t periods;MHFor centralized work The number of industry load (high energy load), MxFor the group number of polymerization appliance load (thermal storage electric boiler).
(2) wind power output constraints
(3) normal power supplies operation constraints
1) technology output bound constrains:
In formula, PGi,maxAnd PGi,minThe respectively output power bound of conventional power unit i.
3) Ramp Rate constrains:
In formula,For conventional power unit i period t-1 output;PGi,upAnd PGi,downThe respectively rising of conventional power unit i Go out power restriction and decline and power restriction.
(4) it polymerize appliance load (thermal storage electric boiler aggregate load) constraint
1) amount of stored heat constrains
Heat storage electric boiler design maximum water temperature is 95 DEG C, if for water temperature more than 95 DEG C, boiler will in storage heater Load operation is reduced, so amount of stored heat should be made in the case where meeting heat demand on daytime next day within prescribed limit:
In formula, amount of stored heatIt is represented by
2) load active power constrains
In formula:For the power load power of t period m classes heat storage electric boilers;WithRespectively m classes store The lower limit and upper limit value of hot type electric boiler power.
3) power swing constrains
The power adjustable of heat storage electric boiler is very high, but in order to ensure the safe and stable operation of electric boiler, its power Fluctuation should limit within limits:
In formula:DPx+ mWithFor the response speed limit of m class heat storage electric boiler ascending, descending power.
(5) centralized industrial load (high energy load) constraint
1) load bound constrains
Wherein,The respectively minimum and maximum capacity limit of high energy load m.
2) load climbing rate constrains
Wherein,And DPh +Power is contributed and rises from for the high energy load m maximum drops allowed from period t-1 to period t Value.
S402:Wind-powered electricity generation active power output plan a few days ago is obtained using model solutionT=1 ..., 96, normal power supplies a few days ago Active power output planT=1 ..., 9,6 centralized 24 point day preload effective power meters of industrial load are drawnT= 1 ..., 9,6 24 point day preload effective power meters of polymerization appliance load are drawnT=1 ..., 96.
The S5 comprises the following steps:
Calculate each distributed appliance load active power plan a few days ago in polymerization appliance load:According to each point of m groups Effective power meter of the cloth appliance load in the t periods divides into
Embodiment 2:
Fig. 2 is a regional power grid schematic diagram, and as example, a kind of wide area source-lotus provided by the invention is a few days ago active excellent Changing control method includes:
S1:Read the related data of all kinds of power supplys, centralized industrial load and distributed appliance load.
(1) in regional power grid, wind power plant rated power 300MW;Normal power supplies unit rated power G1:300MW、G2: 300MW;High energy load quantity MH=1, its rated capacity 50MW;Distributed appliance load N=10, its rated power are respectively 1MW*6,1.25MW*4;It is 0 to send power outside, i.e. Pin=Pout
(2) active prediction data a few days ago is obtained:Next day wind-powered electricity generation active power output prediction data, conventional load power prediction number According to, each centralized industrial load (i.e. high energy load) power prediction data, attached drawing 3 (a) is seen;Each distribution appliance load is (i.e. Thermal storage electric boiler load) power prediction data, see attached drawing 3 (b).
(3) performance parameter of normal power supplies is obtained:
Unit Gi PGi,max PGi,min PGi,up PGi,down
G1 300MW 150MW 30MW/h -30MW/h
G2 300MW 150MW 30MW/h -30MW/h
(4) centralized industrial load regulation performance parameter is obtained:
(5) distributed appliance load regulation performance parameter is obtained:
S2:Polymerization processing is grouped to distributed appliance load, it is as shown in the table.
S3:The polymerization regulation performance parameter of each group aggregate load is calculated, it is as a result as follows:
The secondary daily load prediction power of two groups of polymerization appliance loads is calculated, the result is shown in attached drawing 4;
S4:Build and solve wide area source-lotus active optimization model a few days ago, obtain wind-powered electricity generation, thermoelectricity, centralized industrial load Active power plan a few days ago, is shown in attached drawing 5 (a), polymerize the active power plan a few days ago of appliance load, sees attached drawing 5 (b).
S5:Each distributed load active power plan a few days ago in distribution polymerization appliance load, the result is shown in attached drawing 6.Finally It should be noted that:The above embodiments are merely illustrative of the technical scheme of the present invention and are not intended to be limiting thereof, although with reference to above-mentioned reality Apply example the present invention is described in detail, those of ordinary skill in the art still can be to the specific implementation of the present invention Mode technical scheme is modified or replaced equivalently, these are without departing from any modification of spirit and scope of the invention or equivalent substitution, Applying within pending claims.

Claims (8)

1. a kind of wide area source-lotus active optimization control method a few days ago, comprises the following steps:
S1:Read the related data of all kinds of power supplys, centralized industrial load and distributed appliance load;
S2:Polymerization processing is grouped to distributed appliance load;
S3:Calculate the polymerization regulation performance parameter of each group aggregate load;
S4:Build and solve wide area source-lotus active optimization model a few days ago, obtain all kinds of power supplys, centralized industrial load, the polymerization people With the active power plan a few days ago of load;
S5:The active power plan a few days ago of each distributed load in distribution polymerization appliance load.
A kind of 2. wide area source-lotus according to claim 1 active optimization control method a few days ago, it is characterised in that the S1 Comprise the following steps:
S101:Obtain the quantity M of centralized industrial load in power gridH, the quantity N of distributed appliance load;
S102:Obtain next day wind-powered electricity generation active power output prediction dataConventional load power prediction dataEach centralization Industrial load power prediction dataEach distribution appliance load power prediction data
S103:Obtain the performance parameter of normal power supplies:The peak power output P of conventional power unit iGi,max, minimum output power PGi,min, minimum run timeMinimum downtimeRise and power restriction PGi,up, decline and power restriction PGi,down
S104:Obtain centralized industrial load regulation performance parameter:The maximum load power of centralized industrial load mIt is minimum Load powerMaximum emersion power restriction DPh +, minimum power restriction drops
S105:Obtain the regulation performance parameter of distributed appliance load:The minimum load power of distributed appliance load n Maximum load powerMinimum amount of stored heatMaximum amount of stored heatMaximum climbing rate DPx+ m,n, minimum climbing rate
A kind of 3. wide area source-lotus according to claim 1 active optimization control method a few days ago, it is characterised in that the S2 Comprise the following steps:
S201:N number of distributed appliance load is grouped, thermal storage electric boiler model, specified amount of stored heat, rated power is identical Thermal storage electric boiler load be grouped into the same group, be divided into MxGroup, m group thermal storage electric boilers quantity are Nm
S202:To MxThe processing of distribution type thermal storage electric boiler Load aggregation, concretely comprises the following steps:M groups polymerization appliance load power etc. In NmThe superposition of a thermal storage electric boiler load power, i.e.,M groups polymerization appliance load amount of stored heat is equal to NmA storage The superposition of water-tube boiler load amount of stored heat, i.e.,
A kind of 4. wide area source-lotus according to claim 1 active optimization control method a few days ago, it is characterised in that the S3 Comprise the following steps:
S301:Calculate MxThe secondary daily load prediction power of group polymerization appliance load
S302:Calculate MxThe minimum load power of group polymerization appliance loadAnd maximum load power
S303:Calculate MxThe minimum amount of stored heat of group polymerization appliance loadWith maximum amount of stored heat
S304:Calculate MxThe maximum climbing rate DP of group polymerization appliance loadx+ mWith minimum climbing rate
A kind of 5. wide area source-lotus according to claim 1 active optimization control method a few days ago, it is characterised in that the S4 Comprise the following steps:
S401:Wide area source-lotus active optimization model a few days ago is built, to abandon air quantity minimum, constraints includes sets target function Power-balance constraint, wind power output constraint, normal power supplies operation constraint, polymerization appliance load constraint, the constraint of high energy load;
S402:Solution obtains wind-powered electricity generation active power output plan a few days agoNormal power supplies active power output plan a few days agoCentralized work Industry load day preload effective power meter drawIt polymerize appliance load day preload effective power meter to draw
A kind of 6. wide area source-lotus according to claim 1 active optimization control method a few days ago, it is characterised in that the S5 Comprise the following steps:
Calculate each distributed load active power plan a few days ago in polymerization appliance load
A kind of 7. wide area source-lotus according to claim 4 active optimization control method a few days ago, it is characterised in that the S3 Comprise the following steps:
S301:Calculate Mx24 pre- power scales of point load of next day of group polymerization appliance load:The m groups polymerization appliance load t periods Load prediction power is:
<mrow> <msubsup> <mi>P</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>m</mi> <mn>0</mn> </mrow> </msubsup> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>m</mi> </msub> </munderover> <msubsup> <mi>P</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>m</mi> <mn>0</mn> <mo>,</mo> <mi>n</mi> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
Wherein,The load prediction power of n-th of distributed appliance load t period
S302:Calculate MxThe minimum and maximum regulating power of group polymerization appliance load:
<mrow> <msubsup> <mi>P</mi> <mrow> <mi>x</mi> <mi>min</mi> </mrow> <mi>m</mi> </msubsup> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <msub> <mi>N</mi> <mi>m</mi> </msub> </munderover> <msubsup> <mi>P</mi> <mrow> <mi>x</mi> <mi>min</mi> </mrow> <mrow> <mi>m</mi> <mo>,</mo> <mi>n</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>P</mi> <mrow> <mi>x</mi> <mi>max</mi> </mrow> <mi>m</mi> </msubsup> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <msub> <mi>N</mi> <mi>m</mi> </msub> </munderover> <msubsup> <mi>P</mi> <mrow> <mi>x</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mrow> <mi>m</mi> <mo>,</mo> <mi>n</mi> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
Wherein,For the minimum and maximum regulating power of n-th of distributed appliance load;
S303:Calculate MxThe minimum and maximum amount of stored heat of group polymerization appliance load:
<mrow> <msubsup> <mi>Q</mi> <mrow> <mi>x</mi> <mi>min</mi> </mrow> <mi>m</mi> </msubsup> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <msub> <mi>N</mi> <mi>m</mi> </msub> </munderover> <msubsup> <mi>Q</mi> <mrow> <mi>x</mi> <mi>min</mi> </mrow> <mrow> <mi>m</mi> <mo>,</mo> <mi>n</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>Q</mi> <mrow> <mi>x</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mi>m</mi> </msubsup> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <msub> <mi>N</mi> <mi>m</mi> </msub> </munderover> <msubsup> <mi>Q</mi> <mrow> <mi>x</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mrow> <mi>m</mi> <mo>,</mo> <mi>n</mi> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
Wherein,For the minimum and maximum amount of stored heat of n-th of distributed appliance load;
S304:Calculate MxThe maximum climbing rate and minimum climbing rate of group polymerization appliance load:
DPx+ m=Nmin (DPx+ m,n),
Wherein, DPx+ m,nFor the minimum and maximum climbing rate of n-th of distributed appliance load.
A kind of 8. wide area source-lotus according to claim 5 active optimization control method a few days ago, it is characterised in that S401:Structure Build wide area source-lotus active optimization model a few days ago;
It is as follows to build object function:
<mrow> <mi>min</mi> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>0</mn> </mrow> <mn>24</mn> </munderover> <mrow> <mo>(</mo> <msubsup> <mi>P</mi> <mrow> <mi>W</mi> <mo>,</mo> <mi>f</mi> <mi>o</mi> <mi>r</mi> <mi>e</mi> <mi>c</mi> <mi>a</mi> <mi>s</mi> <mi>t</mi> </mrow> <mi>t</mi> </msubsup> <mo>-</mo> <msubsup> <mi>P</mi> <mi>W</mi> <mi>t</mi> </msubsup> <mo>)</mo> </mrow> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
In formula,For t periods wind-powered electricity generation plan a few days ago power generating value;For period next day t active power for wind power predicted value;Δt For the duration of unit period, it is set as 15 minutes;
Wherein constraints includes:
(1) power-balance constraint:
<mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>G</mi> </msub> </munderover> <msubsup> <mi>U</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <msubsup> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>+</mo> <msubsup> <mi>P</mi> <mi>W</mi> <mi>t</mi> </msubsup> <mo>=</mo> <msubsup> <mi>P</mi> <mi>D</mi> <mi>t</mi> </msubsup> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>H</mi> </msub> </munderover> <msubsup> <mi>P</mi> <mrow> <mi>h</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>m</mi> </msubsup> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>X</mi> </msub> </munderover> <msubsup> <mi>P</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>m</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
(2) wind power output constraints
<mrow> <mn>0</mn> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mi>W</mi> <mi>t</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>W</mi> <mo>,</mo> <mi>f</mi> <mi>o</mi> <mi>r</mi> <mi>e</mi> <mi>c</mi> <mi>a</mi> <mi>s</mi> <mi>t</mi> </mrow> <mi>t</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
(3) normal power supplies operation constraints
1) technology output bound constrains:
<mrow> <msubsup> <mi>U</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <msub> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>U</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <msub> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
In formula, PGi,maxAnd PGi,minThe respectively output power bound of conventional power unit i;
3) Ramp Rate constrains:
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>U</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <msubsup> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>-</mo> <msubsup> <mi>U</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <msubsup> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>&amp;le;</mo> <msub> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>u</mi> <mi>p</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>U</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <msubsup> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>U</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <msubsup> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>&amp;le;</mo> <msub> <mi>P</mi> <mrow> <mi>G</mi> <mi>i</mi> <mo>,</mo> <mi>d</mi> <mi>o</mi> <mi>w</mi> <mi>n</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
In formula,For conventional power unit i period t-1 output;PGi,upAnd PGi,downThe rising of respectively conventional power unit i is contributed Limit and decline and power restriction;
(4) it polymerize appliance load (thermal storage electric boiler aggregate load) constraint
1) amount of stored heat constrains
Heat storage electric boiler design maximum water temperature is 95 DEG C, and amount of stored heat is within prescribed limit:
<mrow> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>Q</mi> <mrow> <mi>x</mi> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> <mi>m</mi> </msubsup> <mo>#</mo> </mrow> </mtd> <mtd> <msubsup> <mi>Q</mi> <mrow> <mi>x</mi> <mi>t</mi> </mrow> <mi>m</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>Q</mi> <mrow> <mi>x</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mi>m</mi> </msubsup> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
In formula, amount of stored heatIt is represented by
2) load active power constrains
<mrow> <mtable> <mtr> <mtd> <msubsup> <mi>P</mi> <mrow> <mi>x</mi> <mi>min</mi> </mrow> <mi>m</mi> </msubsup> </mtd> <mtd> <mrow> <mo>#</mo> <msubsup> <mi>P</mi> <mrow> <mi>x</mi> <mi>t</mi> </mrow> <mi>m</mi> </msubsup> </mrow> </mtd> <mtd> <msubsup> <mi>P</mi> <mrow> <mi>x</mi> <mi>max</mi> </mrow> <mi>m</mi> </msubsup> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
In formula:For the power load power of t period m classes heat storage electric boilers;WithRespectively m classes heat accumulating type The lower limit and upper limit value of electric boiler power;
3) power swing constrains
The power adjustable of heat storage electric boiler is very high, but in order to ensure the safe and stable operation of electric boiler, the ripple of its power It is dynamic to limit within limits:
In formula:DPx+ mWithFor the response speed limit of m class heat storage electric boiler ascending, descending power;
(5) centralized industrial load (high energy load) constraint
1) load bound constrains
<mrow> <msubsup> <mi>P</mi> <mrow> <mi>h</mi> <mi>min</mi> </mrow> <mi>m</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>h</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>m</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>h</mi> <mi>max</mi> </mrow> <mi>m</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
Wherein,The respectively minimum and maximum capacity limit of high energy load m;
2) load climbing rate constrains
<mrow> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>DP</mi> <mi>h</mi> <mo>-</mo> </msubsup> <mo>?</mo> <msubsup> <mi>P</mi> <mrow> <mi>h</mi> <mo>,</mo> <mi>t</mi> </mrow> <mi>m</mi> </msubsup> </mrow> </mtd> <mtd> <msubsup> <mi>P</mi> <mrow> <mi>h</mi> <mo>,</mo> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> <mi>m</mi> </msubsup> </mtd> <mtd> <mrow> <msup> <msub> <mi>P</mi> <mi>h</mi> </msub> <mo>+</mo> </msup> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
Wherein,And DPh +Force value is contributed and rises from for the high energy load m maximum drops allowed from period t-1 to period t.
CN201711036152.7A 2017-10-30 2017-10-30 Wide-area source-load day-ahead active power optimization control method Expired - Fee Related CN107919686B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711036152.7A CN107919686B (en) 2017-10-30 2017-10-30 Wide-area source-load day-ahead active power optimization control method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711036152.7A CN107919686B (en) 2017-10-30 2017-10-30 Wide-area source-load day-ahead active power optimization control method

Publications (2)

Publication Number Publication Date
CN107919686A true CN107919686A (en) 2018-04-17
CN107919686B CN107919686B (en) 2021-08-13

Family

ID=61895001

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711036152.7A Expired - Fee Related CN107919686B (en) 2017-10-30 2017-10-30 Wide-area source-load day-ahead active power optimization control method

Country Status (1)

Country Link
CN (1) CN107919686B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108711889A (en) * 2017-11-07 2018-10-26 华北电力大学 A kind of wide area source-lotus active optimization control method a few days ago
CN109829595A (en) * 2018-09-05 2019-05-31 华北电力大学 A kind of response potentiality quantization method of polymorphic elastic load clustered control
CN110112767A (en) * 2019-03-19 2019-08-09 华北电力大学 The polymorphic Demand-side load of wide area participates in the lotus source optimization control method of peak-load regulating
CN111190355A (en) * 2020-01-14 2020-05-22 华北电力大学 Air conditioner and water heater combined control method and system based on cloud platform
CN111416339A (en) * 2019-03-21 2020-07-14 华北电力大学 Source-load day-ahead active power coordination control method based on double-layer planning model
CN113437764A (en) * 2021-05-21 2021-09-24 国网辽宁省电力有限公司鞍山供电公司 Source-load interaction peak regulation strategy based on active control of fused magnesium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103761680A (en) * 2014-01-16 2014-04-30 华南理工大学 Grid and provincial dispatching method and device for AC/DC interconnected large power grid with wind farm
CN105470987A (en) * 2015-12-23 2016-04-06 华北电力大学 Load-source coordination control based renewable energy day-ahead plan correction method
CN105490312A (en) * 2015-12-18 2016-04-13 国家电网公司 Multi-source reactive power optimization control method for power system
JP2016059126A (en) * 2014-09-08 2016-04-21 株式会社東芝 Power load estimation device, power load estimation method and power load estimation program

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103761680A (en) * 2014-01-16 2014-04-30 华南理工大学 Grid and provincial dispatching method and device for AC/DC interconnected large power grid with wind farm
JP2016059126A (en) * 2014-09-08 2016-04-21 株式会社東芝 Power load estimation device, power load estimation method and power load estimation program
CN105490312A (en) * 2015-12-18 2016-04-13 国家电网公司 Multi-source reactive power optimization control method for power system
CN105470987A (en) * 2015-12-23 2016-04-06 华北电力大学 Load-source coordination control based renewable energy day-ahead plan correction method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李建芳等: "基于自适应粒子群算法的主动配电网日前有功调度", 《南方电网技术》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108711889A (en) * 2017-11-07 2018-10-26 华北电力大学 A kind of wide area source-lotus active optimization control method a few days ago
CN109829595A (en) * 2018-09-05 2019-05-31 华北电力大学 A kind of response potentiality quantization method of polymorphic elastic load clustered control
CN109829595B (en) * 2018-09-05 2023-06-20 华北电力大学 Response potential quantification method for polymorphic elastic load cluster control
CN110112767A (en) * 2019-03-19 2019-08-09 华北电力大学 The polymorphic Demand-side load of wide area participates in the lotus source optimization control method of peak-load regulating
CN110112767B (en) * 2019-03-19 2023-05-30 华北电力大学 Load source optimization control method for peak regulation of wide-area polymorphic demand side load participation system
CN111416339A (en) * 2019-03-21 2020-07-14 华北电力大学 Source-load day-ahead active power coordination control method based on double-layer planning model
CN111190355A (en) * 2020-01-14 2020-05-22 华北电力大学 Air conditioner and water heater combined control method and system based on cloud platform
CN111190355B (en) * 2020-01-14 2021-09-17 华北电力大学 Air conditioner and water heater combined control method and system based on cloud platform
CN113437764A (en) * 2021-05-21 2021-09-24 国网辽宁省电力有限公司鞍山供电公司 Source-load interaction peak regulation strategy based on active control of fused magnesium

Also Published As

Publication number Publication date
CN107919686B (en) 2021-08-13

Similar Documents

Publication Publication Date Title
CN107919686A (en) A kind of wide area source lotus active optimization control method a few days ago
Wang et al. Short-term hydro-thermal-wind-photovoltaic complementary operation of interconnected power systems
JP2019525705A (en) Long-term federation peaking scheduling method for inter-provincial continuous hydropower plants
Wang et al. Renewable production of ammonia and nitric acid
CN107276096A (en) The distribution network voltage control method that a kind of photovoltaic is coordinated and optimized with air conditioner load
CN109936170A (en) Consider the honourable extreme misery complementation coordination optimization dispatching method of power supply flexibility nargin
CN103426032A (en) Method for economically and optimally dispatching cogeneration units
CN107332286A (en) A kind of cogeneration of heat and power containing heat accumulation and wind-powered electricity generation coordinated scheduling method
CN102623989A (en) Method for optimization and configuration of intermittent distributed generation (DG)
CN107634547B (en) Power output control method of electric combined system based on new energy output prediction error
CN104578176A (en) Method for making power generation plan in consideration of direct current interaction
CN104809545B (en) A kind of virtual plant runs modeling method
CN104143839B (en) Wind power plant cluster based on power prediction limits active power distribution method of exerting oneself
CN104617585A (en) Reactive compensation configuration method
CN112615370A (en) Wind power consumption coordination control method based on electric smelting magnesium load
CN103593711A (en) Distributed power source optimal configuration method
CN114330827A (en) Distributed robust self-scheduling optimization method for multi-energy flow virtual power plant and application thereof
CN115689166A (en) Method and system for aggregated utilization of regional distributed energy resources
CN108711889A (en) A kind of wide area source-lotus active optimization control method a few days ago
CN103280845A (en) Water-fire-electricity system trans-provincial peak-regulating load distribution method
CN106684931A (en) Active control method and system of new-energy power plant
CN102684228A (en) Method for optimizing configuration of intermittent distribution type power supply based on complementary
CN113240205A (en) Regional energy utilization system substitution optimization method based on multi-energy comprehensive utilization
CN109829595A (en) A kind of response potentiality quantization method of polymorphic elastic load clustered control
CN108062022A (en) A kind of optimal control method of thermoelectricity cooperative system

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

Granted publication date: 20210813