CN107425521A - A kind of dynamic aggregation modeling for considering load fluctuation and participation frequency modulation method - Google Patents

A kind of dynamic aggregation modeling for considering load fluctuation and participation frequency modulation method Download PDF

Info

Publication number
CN107425521A
CN107425521A CN201710508974.4A CN201710508974A CN107425521A CN 107425521 A CN107425521 A CN 107425521A CN 201710508974 A CN201710508974 A CN 201710508974A CN 107425521 A CN107425521 A CN 107425521A
Authority
CN
China
Prior art keywords
load
mrow
msub
msubsup
response
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.)
Pending
Application number
CN201710508974.4A
Other languages
Chinese (zh)
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.)
State Grid Corp of China SGCC
Southeast University
China Electric Power Research Institute Co Ltd CEPRI
State Grid Jiangsu Electric Power Co Ltd
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Southeast University
China Electric Power Research Institute Co Ltd CEPRI
State Grid Jiangsu Electric Power Co Ltd
Electric Power Research Institute of State Grid Jiangsu 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 State Grid Corp of China SGCC, Southeast University, China Electric Power Research Institute Co Ltd CEPRI, State Grid Jiangsu Electric Power Co Ltd, Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201710508974.4A priority Critical patent/CN107425521A/en
Publication of CN107425521A publication Critical patent/CN107425521A/en
Pending legal-status Critical Current

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
    • 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/28Arrangements for balancing of the load in a network by storage of energy
    • 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]
    • 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/003Load forecast, e.g. methods or systems for forecasting future load demand

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Cosmetics (AREA)

Abstract

The invention discloses a kind of dynamic aggregation modeling for considering load fluctuation and participate in frequency modulation method, there is provided various types of loads are docked in abstract synthesis polymerization model, model side, integrate load general character, ignore the load opposite sex;Opposite side docks electric power networks, aggregate load attribute and demand response characteristic under application scenarios needed for directly showing, such as respond capacity, response time, response speed, demand response is by grasping the load structure with analysis system, power consumer can be guided to select the rational electricity consumption time, or use rational way of energy storage, have the function that peak load shifting, cluster analysis can be applied to carry out part throttle characteristics research, constantly system frequency is changed using dynamic aggregation body and carries out dynamic corrections, effective help is provided for power system security stable operation, have a good application prospect.

Description

A kind of dynamic aggregation modeling for considering load fluctuation and participation frequency modulation method
Technical field
The present invention relates to intelligent grid demand response technical field, and in particular to a kind of dynamic for considering load fluctuation is gathered Build mould jointly and participate in frequency modulation method.
Background technology
Application demand response technology is one of core technology of intelligent grid, and application demand response technology can be excavated fully Load side resource, realize the integrated application of resource.At present, a considerable amount of potential demand responses in middle-size and small-size power consumer be present Resource is not developed, and typical Load Regulation resource includes air-conditioning, refrigerator, electric automobile, energy storage and accumulation of energy etc..These load kinds Class is numerous, and individual devices power consumption is small, as single air conditioner cooling in summer electrical power is no more than 5kWh.It is more and more huger in face of colony Greatly, the more and more flexible power consumer of electricity consumption, before Demand Side Response measure is implemented, it is not necessary that be also impossible to each user Modeling Research is all gone one by one, but they can be compared with the typical user's classification having been acknowledged by the method for cluster, Go out their part throttle characteristics by that analogy, so as to back-up needs response project implementation, greatly improve operating efficiency.
For power system, the burdened resource scheduling of middle-size and small-size power consumer is difficult, and part throttle characteristics is mixed and disorderly, in order to assist Adjust using these controllable resources, it is necessary to be polymerize to it and managed.It is necessary to find a kind of abstract synthesis polymerization model, mould Various types of loads are docked in type side, integrate load general character, ignore the load opposite sex;Opposite side docks electric power networks, directly opens up Aggregate load attribute and demand response characteristic under application scenarios needed for revealing, such as respond capacity, response time, response speed Deng.Demand response can guide power consumer to select rational electricity consumption time by grasping and the load structure of analysis system, or Using rational way of energy storage, have the function that peak load shifting, cluster analysis can be applied to carry out part throttle characteristics research, be current Need to solve the problems, such as.
The content of the invention
The purpose of the present invention is that the burdened resource for overcoming existing middle-size and small-size power consumer dispatches difficulty, and part throttle characteristics is miscellaneous Disorderly, the problem of how load is polymerize and managed.The dynamic aggregation modeling and participation of the consideration load fluctuation of the present invention Frequency modulation method, can apply cluster analysis carry out part throttle characteristics research, constantly system frequency is changed using dynamic aggregation body into Mobile state amendment, effective help is provided for power system security stable operation, had a good application prospect.
In order to achieve the above object, the technical solution adopted in the present invention is:
A kind of dynamic aggregation modeling for considering load fluctuation and participation frequency modulation method, it is characterised in that:Including following step Suddenly,
Step (A), based on the affiliated category of employment of small flexible load, interactive, controllability, the characteristic of randomness, formed Flexible load universal identification model;
Step (B), based on small flexible load responding time, response speed, the demand response characteristic for responding capacity, formed Flexible load demand response model;
Step (C), the influent factor based on external environment to Load aggregation, form demand response Context awareness model;
Step (D), consider dynamic aggregation requirement of the different target line for flexible load, establish the polymerization of flexible load general character Method;
Step (E), the polymerization formed based on different aggregate targets is overall, forms flexible load condensate synthesis mould Type;
Step (F), consider that condensate participates in power system assistant service requirement, form condensate response capacity-constrained bar Part;
Step (G), consider that condensate participates in stabilizing the output requirement of distributed power source, forming condensate regulating index will Ask;
Step (H), the fast response characteristic possessed based on condensate, form flexible load dynamic aggregation model and participate in system Frequency modulation method.
A kind of dynamic aggregation modeling of foregoing consideration load fluctuation and participation frequency modulation method, it is characterised in that:Step (A) flexible load universal identification model, is formed, as shown in formula (1) and formula (2),
Load={ loadi|1≤i≤I} (1)
loadi=(Num (t), YD, YC, Ys, YT) (2)
Wherein, Load is the set of various small flexible loads;loadiRepresent the i-th class small flexible load;I is area under one's jurisdiction Interior small-sized load sum, Num (t) is the number of load t;It is interactive that YD represents load, and load is shown as two-way when YD is 1 Interaction, load is shown as unidirectionally when YD is 0;YC represents load controllability, YCFor 1 when load show as can control, YCFor 0 when bear Lotus shows as uncontrollable;YSRepresent load randomness, YSFor 1 when load show as randomness, YSFor 0 when load show as really It is qualitative;YTRepresent load industry type, YTFor 1 when load be industrial load, YTFor 2 when be Commercial Load, YTFor 3 when load be Resident load, YTFor 4 when be other loads.
A kind of dynamic aggregation modeling of foregoing consideration load fluctuation and participation frequency modulation method, it is characterised in that:Step (B) flexible load demand response model, is formed, as shown in formula (3) and formula (4),
loadi={ lj|1≤j≤J} (3)
lj=(PN, f (t), g (t), T) and (4)
Wherein, ljFor i class controllable burdens loadiIn single model j, J is controllable burden loadiMiddle model sum, PN are The load rated capacity;F (t) is the modulating properties curve of the load;G (t) is the load characteristic curve of the load;T is that this is negative The fastest response time of lotus.
A kind of dynamic aggregation modeling of foregoing consideration load fluctuation and participation frequency modulation method, it is characterised in that:Step (C) demand response Context awareness model, is formed, as shown in formula (5) and formula (6),
Environ={ environh|1≤h≤H} (6)
environh=(Temp (t), Hum (t), I, XC) (7)
Wherein, Environ is the set of local demand response load local environment;environhFor H type environment category, Local common H kinds environment category;Temp (t) is this kind of environment environjUnder external temperature;Hum (t) is under this kind of environment Humidity;I indicates workaday division, and I=0 is nonworkdays, and I=1 is working day;XCFor under the environment for demand response Responding ability, be divided into 1,2,3, show faint response, medium response, complete response successively.
A kind of dynamic aggregation modeling of foregoing consideration load fluctuation and participation frequency modulation method, it is characterised in that:Step (D), consider dynamic aggregation requirement of the different target line for flexible load, establish flexible load general character polymerization, including with Lower step,
(D1), from n flexible load optional k object as initial center;
(D2) be, Euclidean distance according to initial center and the similarity of cluster centre, the similarity, assign them to Its most like cluster centre cluster;
(D3) cluster centre of obtained new class, is recalculated with average, constantly repeats this process until standard deviation is estimated Function convergence, according to the difference of application scenarios, the K-means canonical measure functions of selection are also different, as shown in formula (7),
Wherein, E is sample and cluster centre mean square deviation sum;P is a point in the space for represent object;miFor cluster CiAverage;
(D4), the sample based on time scale and cluster centre mean square deviation sum Et, as shown in formula (8),
Wherein, loadi.lj.T it is a point in the space for represent object;For the cluster based on time scale Value;ItFor the load load included under the clusteriSum;JtTo cluster body loadiIn single load ljSum;
(D5) sample and cluster centre mean square deviation sum E, based on response capacityP, as shown in formula (9),
Wherein, loadi.lj.f (t) is a point in the space for represent object;For the cluster based on time scale Average;IPFor the load load included under the clusteriSum;JPTo cluster body loadiIn single load ljSum.
A kind of dynamic aggregation modeling of foregoing consideration load fluctuation and participation frequency modulation method, it is characterised in that:Step (E) flexible load condensate collective model, is formed, as shown in formula (10) and formula (11),
Aggre={ aggk|1≤k≤K} (10)
aggk=(PN, R (t), T, v) and (11)
Wherein, aggre is the set after similar Load aggregation, and k is the classification number of polymerization, for J load loadi, aggregate number K≤2lnJ;aggkFor similar Load aggregation;PNFor aggkPolymerization total capacity;R (t) is aggkIt is adjustable Range of capacity, T are load maximum response time, and v is load average response speed.
A kind of dynamic aggregation modeling of foregoing consideration load fluctuation and participation frequency modulation method, it is characterised in that:It is if same Type load polymerize aggk, include load lj(j=1,2,3 ..., jk), for controllable burden (YC=1), each parameter of polymerization model with Load in condensate has following relation:
Wherein, regulation Potential curves and load characteristic curve after R (t), G (t) respectively polymerization, T=0, v=∞.
A kind of dynamic aggregation modeling of foregoing consideration load fluctuation and participation frequency modulation method, it is characterised in that:Step (F) condensate response capacity constraints, are formed, as shown in formula (17),
aggk.Pdown≤aggk.PDR≤aggk.PN-aggk.Pcontact (17)
Wherein, aggk.PDRIt is the real response capacity that condensate participates in power system assistant service requirement;aggk.Pcontact It is the contract capacity being had been sold out in condensate;aggk.PdownIt is by caused by external environment condition and self psychology constraint Respond capacity-constrained lower limit;aggk.PNFor condensate rated capacity.
A kind of dynamic aggregation modeling of foregoing consideration load fluctuation and participation frequency modulation method, it is characterised in that:Step (G) requirement of condensate regulating index, is formed, it is desirable to which wind power plant per minute active power fluctuation of the installed capacity less than 30MW should be small , can be grid-connected in 20%, it is specific as formula (18) and formula (19) are shown,
Wherein, Pw tRepresent the wind power output of t, unit kW, Pw t-1When representing t-1
The wind power output at quarter;aggk.R (t) is the Load Regulation total amount after Load aggregation, is had positive and negative
Number, it is negative when needing reduction plans, when needing to increase load for just.
A kind of dynamic aggregation modeling of foregoing consideration load fluctuation and participation frequency modulation method, it is characterised in that:Step (H) fast response characteristic, possessed based on condensate, form flexible load dynamic aggregation model and participate in system frequency modulation method, by Adjusted in real time in polymeric response characteristic, the changed power of response system primary frequency modulation second level is capable of from frequency modulation performance It is required that in this manner, the work(frequency characteristic of generator primary frequency modulation is simulated to condensate first, it is then automatic by condensate Detect the frequency departure of power system and change payload, be embodied in when system frequency declines and reduce or even close and drop Oligopolymer responds capacity;Increase condensate response capacity when system frequency rises.
The beneficial effects of the invention are as follows:The dynamic aggregation modeling of the consideration load fluctuation of the present invention and participation frequency modulation side Method, there is provided various types of loads are docked in abstract synthesis polymerization model, model side, integrate load general character, to ignore load different Property;Opposite side docks electric power networks, directly shows aggregate load attribute and demand response characteristic under required application scenarios, Capacity, response time, response speed are such as responded, demand response can guide electricity by grasping the load structure with analysis system Power user selects the rational electricity consumption time, or uses rational way of energy storage, has the function that peak load shifting, can apply cluster Analysis carries out part throttle characteristics research, and constantly system frequency is changed using dynamic aggregation body and carries out dynamic corrections, is power system Safe and stable operation provides effective help, has a good application prospect.
Brief description of the drawings
Fig. 1 is the dynamic aggregation modeling of the consideration load fluctuation of the present invention and the flow chart for participating in frequency modulation method;
Fig. 2 is the simulation model figure based on polymeric single region primary frequency modulation of the present invention;
Fig. 3 is the illustraton of model of the condensate FM module of the present invention.
Embodiment
Below in conjunction with Figure of description, the present invention is further illustrated.
As shown in figure 1, dynamic aggregation modeling and the participation frequency modulation method of the consideration load fluctuation of the present invention, including it is following Step,
Step (A), based on the affiliated category of employment of small flexible load, interactive, controllability, the characteristic of randomness, formed Flexible load universal identification model, as shown in formula (1) and formula (2),
Load={ loadi|1≤i≤I} (1)
loadi=(Num (t), YD, YC, Ys, YT) (2)
Wherein, Load is the set of various small flexible loads;loadiRepresent the i-th class small flexible load;I is area under one's jurisdiction Interior small-sized load sum;Num (t) is the number of load t;It is interactive that YD represents load, and load is shown as two-way when YD is 1 Interaction, load is shown as unidirectionally when YD is 0, such as electric automobile is two-way, and air conditioner load is unidirectional;It is controllable that YC represents load Property, load is shown as can control when YC is 1, and load is shown as uncontrollable when YC is 0, such as air-conditioning is controllable, and wind-powered electricity generation is It is uncontrollable;YS represents load randomness, and load shows as randomness when YS is 1, and load shows as certainty, example when YS is 0 If air-conditioning is randomness, distributed fuel cell is certainty;YT represents load industry type, and load is that industry is negative when YT is 1 Lotus, it is Commercial Load when YT is 2, load is resident load when YT is 3, is other loads when YT is 4, by taking air-conditioning as an example, air-conditioning Load identification model be expressed as
loadi=(0,1,1,1/2/3/4);
Step (B), based on small flexible load responding time, response speed, the demand response characteristic for responding capacity, formed Flexible load demand response model, as shown in formula (3) and formula (4),
loadi={ lj|1≤j≤J} (3)
lj=(PN, f (t), g (t), T) and (4)
Wherein, ljFor i class controllable burdens loadiIn single model j, J is controllable burden loadiMiddle model sum, PN are The load rated capacity;F (t) is the modulating properties curve of the load;G (t) is the load characteristic curve of the load;T is that this is negative The fastest response time of lotus;
Step (C), the influent factor based on external environment to Load aggregation, demand response Context awareness model is formed, born Polymerization of the environment to load residing for lotus also has a certain impact, for space scale, the close load group in geographical position, outside The probability density characteristicses of portion's ambient parameter are similar, and the probability density characteristicses of the load group's external environment condition parameter of geographical position farther out Differ greatly, moreover, the level of enforcement or feedback degree that the loads of different regions encourages for same demand response also have Difference;For time scale, working day and nonworkdays, the difference of the time such as festivals or holidays also have for the polymerization property of load Large effect, therefore, forms demand response Context awareness model, as shown in formula (5) and formula (6),
Environ={ environh|1≤h≤H} (6)
environh=(Temp (t), Hum (t), I, XC) (7)
Wherein, Environ is the set of local demand response load local environment;environhFor H type environment category, Local common H kinds environment category;Temp (t) is this kind of environment environjUnder external temperature;Hum (t) is under this kind of environment Humidity;I indicates workaday division, and I=0 is nonworkdays, and I=1 is working day;XCFor under the environment for demand response Responding ability, be divided into 1,2,3, show faint response, medium response, complete response successively;
Step (D), consider dynamic aggregation requirement of the different target line for flexible load, establish the polymerization of flexible load general character Method, by the similar Load aggregation of regulating characteristic curve together.Obtaining each load ljRegulation elasticity curve f (t) after, According to the demand of application scenarios or electricity market, automatic cluster is carried out to point of proximity using K-means algorithms, specifically included following Step,
(D1), from n flexible load optional k object as initial center;
(D2) be, Euclidean distance according to initial center and the similarity of cluster centre, the similarity, assign them to Its most like cluster centre cluster;
(D3) cluster centre of obtained new class, is recalculated with average, constantly repeats this process until standard deviation is estimated Function convergence, according to the difference of application scenarios, the K-means canonical measure functions of selection are also different, as shown in formula (7),
Wherein, E is sample and cluster centre mean square deviation sum;P is a point in the space for represent object;miFor cluster CiAverage;
(D4), the sample based on time scale and cluster centre mean square deviation sum Et, as shown in formula (8),
Wherein, loadi.lj.T it is a point in the space for represent object;For the cluster based on time scale Value;ItFor the load load included under the clusteriSum;JtTo cluster body loadiIn single load ljSum;
(D5) sample and cluster centre mean square deviation sum E, based on response capacityP, as shown in formula (9),
Wherein, loadi.lj.f (t) is a point in the space for represent object;For the cluster based on time scale Average;IPFor the load load included under the clusteriSum;JPTo cluster body loadiIn single load ljSum.
Step (E), the polymerization formed based on different aggregate targets is overall, forms flexible load condensate synthesis mould Type, as shown in formula (10) and formula (11),
Aggre={ aggk|1≤k≤K} (10)
aggk=(PN, R (t), T, v) and (11)
Wherein, aggre is the set after similar Load aggregation, and k is the classification number of polymerization, for J load loadi, aggregate number K≤2lnJ;aggkFor similar Load aggregation;PNFor aggkPolymerization total capacity;R (t) is aggkIt is adjustable Range of capacity, T are load maximum response time, and v is load average response speed.
A kind of dynamic aggregation modeling of foregoing consideration load fluctuation and participation frequency modulation method, it is characterised in that:It is if same Type load polymerize aggk, include load lj(j=1,2,3 ..., jk), for controllable burden (YC=1), each parameter of polymerization model with Load in condensate has following relation:
Wherein, regulation Potential curves and load characteristic curve after R (t), G (t) respectively polymerization, T=0, v=∞;
Step (F), consider that condensate participates in power system assistant service requirement, form condensate response capacity-constrained bar Part, as shown in formula (17),
aggk.Pdown≤aggk.PDR≤aggk.PN-aggk.Pcontact (17)
Wherein, aggk.PDRIt is the real response capacity that condensate participates in power system assistant service requirement;aggk.Pcontact It is the contract capacity being had been sold out in condensate;aggk.PdownIt is by caused by external environment condition and self psychology constraint Respond capacity-constrained lower limit;aggk.PNFor condensate rated capacity;
Step (G), consider that condensate participates in stabilizing the output requirement of distributed power source, forming condensate regulating index will Ask, according to State Grid Corporation of China's wind power plant access power system technology regulation, it is desirable to which installed capacity is every less than 30MW wind power plant Minute active power fluctuation should be less than 20%, can be grid-connected, specific as shown in formula (18) and formula (19),
Wherein, Pw tRepresent the wind power output of t, unit kW, Pw t-1Represent the wind power output at t-1 moment;aggk.R(t) For the Load Regulation total amount after Load aggregation, there is sign, be negative when needing reduction plans, be when needing and increasing load Just;
Step (H), the fast response characteristic possessed based on condensate, form flexible load dynamic aggregation model and participate in system Frequency modulation method, as shown in Figures 2 and 3, because polymeric response characteristic adjusts in real time, it can be responded from frequency modulation performance and be The changed power requirement for primary frequency modulation second level of uniting, in this manner, simulates generator primary frequency modulation to condensate first Work(frequency characteristic, then by condensate automatic detection power system frequency departure and change payload, be embodied in when being System frequency reduces when declining even closes and reduces condensate response capacity;Increase condensate response when system frequency rises to hold Amount.
The present invention consideration load fluctuation dynamic aggregation modeling and participate in frequency modulation method, first, middle-size and small-size user with Operator signs a contract, there is provided the response capacity of load, response time, response speed information.(it is difficult to directly obtain for some The user of response message, such as air conditioner load, it can be believed by gathering the correlation such as preference temperature section, architectural environment of user Breath, is obtained indirectly), then, in face of the load responding information of magnanimity, with reference to different optimization aims, establish different polymerization moulds Type f (x1,x2,x3), including using the response time it is polymerization site or to respond capacity as polymerization site, obtains corresponding poly- Gregarious body;Finally, when power system operator needs to be scheduled to meet a certain frequency modulation or active service load, root According to the response requirement needed for the assistant service, unified instruction is sent to a certain zoarium of birdsing of the same feather flock together.
In summary, the dynamic aggregation modeling of consideration load fluctuation of the invention and participation frequency modulation method, there is provided abstract Synthesis polymerization model, model side docks various types of loads, integrates load general character, ignores the load opposite sex;Opposite side pair Connect electric power networks, directly show needed for aggregate load attribute and demand response characteristic under application scenarios, such as respond capacity, Response time, response speed, demand response can guide power consumer selection to close by grasping the load structure with analysis system The electricity consumption time of reason, or rational way of energy storage is used, have the function that peak load shifting, cluster analysis can be applied to carry out load Characteristic research, constantly system frequency is changed using dynamic aggregation body and carries out dynamic corrections, is power system security stable operation Effective help is provided, had a good application prospect.
General principle, the main features and advantages of the present invention have been shown and described above.The technical staff of the industry should Understand, the present invention is not limited to the above embodiments, the original for simply illustrating the present invention described in above-described embodiment and specification Reason, without departing from the spirit and scope of the present invention, various changes and modifications of the present invention are possible, these changes and improvements It all fall within the protetion scope of the claimed invention.The claimed scope of the invention is by appended claims and its equivalent circle It is fixed.

Claims (10)

1. a kind of dynamic aggregation modeling for considering load fluctuation and participation frequency modulation method, it is characterised in that:Comprise the following steps,
Step (A), based on the affiliated category of employment of small flexible load, interactive, controllability, the characteristic of randomness, formed flexible Load universal identification model;
Step (B), based on small flexible load responding time, response speed, the demand response characteristic for responding capacity, formed flexible Workload demand response model;
Step (C), the influent factor based on external environment to Load aggregation, form demand response Context awareness model;
Step (D), consider dynamic aggregation requirement of the different target line for flexible load, establish flexible load general character polymerization side Method;
Step (E), the polymerization formed based on different aggregate targets is overall, forms flexible load condensate collective model;
Step (F), consider that condensate participates in power system assistant service requirement, form condensate response capacity constraints;
Step (G), consider that condensate participates in stabilizing the output requirement of distributed power source, form the requirement of condensate regulating index;
Step (H), the fast response characteristic possessed based on condensate, form flexible load dynamic aggregation model and participate in system frequency modulation Method.
2. a kind of dynamic aggregation modeling for considering load fluctuation according to claim 1 and participation frequency modulation method, it is special Sign is:Step (A), flexible load universal identification model is formed, as shown in formula (1) and formula (2),
Load={ loadi|1≤i≤I} (1)
loadi=(Num (t), YD, YC, YS, YT) (2)
Wherein, Load is the set of various small flexible loads;loadiRepresent the i-th class small flexible load;I is local small-sized Load sum;Num (t) is the number of load t;YDIt is interactive to represent load, YDFor 1 when load show as two-way interaction, YD For 0 when load show as unidirectionally;YCRepresent load controllability, YCFor 1 when load show as can control, YCFor 0 when load show as It is uncontrollable;YSRepresent load randomness, YSFor 1 when load show as randomness, YSFor 0 when load show as certainty;YTGeneration Table load industry type, YTFor 1 when load be industrial load, YTFor 2 when be Commercial Load, YTFor 3 when load be resident load, YTFor 4 when be other loads.
3. a kind of dynamic aggregation modeling for considering load fluctuation according to claim 1 and participation frequency modulation method, it is special Sign is:Step (B), flexible load demand response model is formed, as shown in formula (3) and formula (4),
loadi={ lj|1≤j≤J} (3)
lj=(PN, f (t), g (t), T) and (4)
Wherein, ljFor i class controllable burdens loadiIn single model j, J is controllable burden loadiMiddle model sum, PNIt is negative for this Lotus rated capacity;F (t) is the modulating properties curve of the load;G (t) is the load characteristic curve of the load;T is the load The fastest response time.
4. a kind of dynamic aggregation modeling for considering load fluctuation according to claim 1 and participation frequency modulation method, it is special Sign is:Step (C), demand response Context awareness model is formed, as shown in formula (5) and formula (6),
Environ={ environh|1≤h≤H} (6)
environh=(Temp (t), Hum (t), I, XC) (7)
Wherein, Environ is the set of local demand response load local environment;environhFor H type environment category, area under one's jurisdiction Interior common H kinds environment category;Temp (t) is this kind of environment environjUnder external temperature;Hum (t) is wet under this kind of environment Degree;I indicates workaday division, and I=0 is nonworkdays, and I=1 is working day;XCFor the sound under the environment for demand response Should be able to power, be divided into 1,2,3, show faint response, medium response, complete response successively.
5. a kind of dynamic aggregation modeling for considering load fluctuation according to claim 1 and participation frequency modulation method, it is special Sign is:Step (D), consider dynamic aggregation requirement of the different target line for flexible load, establish the polymerization of flexible load general character Method, comprise the following steps,
(D1), from n flexible load optional k object as initial center;
(D2) it is, Euclidean distance according to initial center and the similarity of cluster centre, the similarity, assigns them to it most Similar cluster centre cluster;
(D3) cluster centre of obtained new class, is recalculated with average, constantly repeats this process until standard deviation measure function Convergence, according to the difference of application scenarios, the K-means canonical measure functions of selection are also different, as shown in formula (7),
<mrow> <mi>E</mi> <mo>=</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </msubsup> <msub> <mi>&amp;Sigma;</mi> <mrow> <mi>p</mi> <mo>&amp;Element;</mo> <msub> <mi>C</mi> <mi>i</mi> </msub> </mrow> </msub> <mo>|</mo> <mi>p</mi> <mo>-</mo> <msub> <mi>m</mi> <mi>i</mi> </msub> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
Wherein, E is sample and cluster centre mean square deviation sum;P is a point in the space for represent object;miTo cluster Ci Value;
(D4), the sample based on time scale and cluster centre mean square deviation sum Et, as shown in formula (8),
<mrow> <msup> <mi>E</mi> <mi>t</mi> </msup> <mo>=</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>I</mi> <mi>t</mi> </msup> </msubsup> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>J</mi> <mi>t</mi> </msup> </msubsup> <msub> <mi>&amp;Sigma;</mi> <mrow> <msub> <mi>load</mi> <mi>i</mi> </msub> <mo>.</mo> <msub> <mi>l</mi> <mi>j</mi> </msub> <mo>.</mo> <mi>T</mi> <mo>&amp;Element;</mo> <msubsup> <mi>C</mi> <mi>i</mi> <mi>t</mi> </msubsup> </mrow> </msub> <mo>|</mo> <msub> <mi>load</mi> <mi>i</mi> </msub> <mo>.</mo> <msub> <mi>l</mi> <mi>j</mi> </msub> <mo>.</mo> <mi>T</mi> <mo>-</mo> <msubsup> <mi>m</mi> <mi>i</mi> <mi>t</mi> </msubsup> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
Wherein, loadi.lj.T it is a point in the space for represent object;For the cluster based on time scaleAverage;It For the load load included under the clusteriSum;JtTo cluster body loadiIn single load ljSum;
(D5) sample and cluster centre mean square deviation sum E, based on response capacityP, as shown in formula (9),
<mrow> <msup> <mi>E</mi> <mi>P</mi> </msup> <mo>=</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>I</mi> <mi>P</mi> </msup> </msubsup> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>J</mi> <mi>P</mi> </msup> </msubsup> <msub> <mi>&amp;Sigma;</mi> <mrow> <msub> <mi>load</mi> <mi>i</mi> </msub> <mo>.</mo> <msub> <mi>l</mi> <mi>j</mi> </msub> <mo>.</mo> <mi>T</mi> <mo>&amp;Element;</mo> <msubsup> <mi>C</mi> <mi>i</mi> <mi>P</mi> </msubsup> </mrow> </msub> <mo>|</mo> <msub> <mi>load</mi> <mi>i</mi> </msub> <mo>.</mo> <msub> <mi>l</mi> <mi>j</mi> </msub> <mo>.</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <msubsup> <mi>m</mi> <mi>i</mi> <mi>P</mi> </msubsup> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
Wherein, loadi.lj.f (t) is a point in the space for represent object;For the cluster based on time scale Value;IPFor the load load included under the clusteriSum;JPTo cluster body loadiIn single load ljSum.
6. a kind of dynamic aggregation modeling for considering load fluctuation according to claim 1 and participation frequency modulation method, it is special Sign is:Step (E), flexible load condensate collective model is formed, as shown in formula (10) and formula (11),
Aggre={ aggk|1≤k≤K} (10)
aggk=(PN, R (t), T, v) and (11)
Wherein, aggre is the set after similar Load aggregation, and k is the classification number of polymerization, for the load with J loadi, Aggregate number K≤2lnJ;aggkFor similar Load aggregation;PNFor aggkPolymerization total capacity;R (t) is aggkVariable capacity model Enclose, T is load maximum response time, and v is load average response speed.
7. a kind of dynamic aggregation modeling for considering load fluctuation according to claim 6 and participation frequency modulation method, it is special Sign is:If similar Load aggregation aggk, include load lj(j=1,2,3 ..., jk), for controllable burden (YC=1), it polymerize Each parameter of model has following relation with the load in condensate:
<mrow> <msub> <mi>agg</mi> <mi>k</mi> </msub> <mo>.</mo> <msub> <mi>P</mi> <mi>N</mi> </msub> <mo>=</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>j</mi> <mi>k</mi> </msub> </msubsup> <msub> <mi>l</mi> <mi>j</mi> </msub> <mo>.</mo> <msub> <mi>P</mi> <mi>N</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>T</mi> <mo>&amp;GreaterEqual;</mo> <msubsup> <mo>&amp;cap;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>j</mi> <mi>k</mi> </msub> </msubsup> <msub> <mi>l</mi> <mi>j</mi> </msub> <mo>.</mo> <mi>t</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>R</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>j</mi> <mi>k</mi> </msub> </msubsup> <msub> <mi>l</mi> <mi>j</mi> </msub> <mo>.</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>G</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>j</mi> <mi>k</mi> </msub> </msubsup> <msub> <mi>l</mi> <mi>j</mi> </msub> <mo>.</mo> <mi>g</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>v</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mi>T</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
Wherein, regulation Potential curves and load characteristic curve after R (t), G (t) respectively polymerization, T=0, v=∞.
8. a kind of dynamic aggregation modeling for considering load fluctuation according to claim 1 and participation frequency modulation method, it is special Sign is:Step (F), condensate response capacity constraints are formed, as shown in formula (17),
aggk.Pdown≤aggk.PDR≤aggk.PN-aggk.Pcontact (17)
Wherein, aggk.PDRIt is the real response capacity that condensate participates in power system assistant service requirement;aggk.PcontactIt is poly- The contract capacity being had been sold out in zoarium;aggk.PdownIt is to be responded by caused by external environment condition and self psychology constraint Capacity-constrained lower limit;aggk.PNFor condensate rated capacity.
9. a kind of dynamic aggregation modeling for considering load fluctuation according to claim 1 and participation frequency modulation method, it is special Sign is:Step (G), form the requirement of condensate regulating index, it is desirable to which installed capacity is per minute active less than 30MW wind power plant Power swing should be less than 20%, can be grid-connected, specific as shown in formula (18) and formula (19),
<mrow> <mfrac> <mrow> <mo>|</mo> <msubsup> <mi>P</mi> <mi>w</mi> <mi>t</mi> </msubsup> <mo>-</mo> <msubsup> <mi>P</mi> <mi>w</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>|</mo> </mrow> <msubsup> <mi>P</mi> <mi>w</mi> <mi>t</mi> </msubsup> </mfrac> <mo>&amp;le;</mo> <mn>0.2</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mrow> <mo>|</mo> <msubsup> <mi>P</mi> <mi>w</mi> <mi>t</mi> </msubsup> <mo>+</mo> <msub> <mi>agg</mi> <mi>k</mi> </msub> <mo>.</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <msubsup> <mi>P</mi> <mi>w</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>|</mo> </mrow> <msubsup> <mi>P</mi> <mi>w</mi> <mi>t</mi> </msubsup> </mfrac> <mo>&amp;le;</mo> <mn>0.2</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>
Wherein, Pw tRepresent the wind power output of t, unit kW, Pw t-1Represent the wind power output at t-1 moment;aggk.R (t) is negative Load Regulation total amount after lotus polymerization, there is sign, is negative when needing reduction plans, when needing to increase load for just.
10. a kind of dynamic aggregation modeling for considering load fluctuation according to claim 1 and participation frequency modulation method, it is special Sign is:Step (H), the fast response characteristic possessed based on condensate, form flexible load dynamic aggregation model and participate in system Frequency modulation method, being capable of the response system primary frequency modulation second from frequency modulation performance because polymeric response characteristic adjusts in real time The changed power requirement of level, in this manner, simulates the work(frequency characteristic of generator primary frequency modulation, then to condensate first By condensate automatic detection power system frequency departure and change payload, be embodied in and subtract when system frequency declines It is small in addition close and reduce condensate response capacity;Increase condensate response capacity when system frequency rises.
CN201710508974.4A 2017-06-28 2017-06-28 A kind of dynamic aggregation modeling for considering load fluctuation and participation frequency modulation method Pending CN107425521A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710508974.4A CN107425521A (en) 2017-06-28 2017-06-28 A kind of dynamic aggregation modeling for considering load fluctuation and participation frequency modulation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710508974.4A CN107425521A (en) 2017-06-28 2017-06-28 A kind of dynamic aggregation modeling for considering load fluctuation and participation frequency modulation method

Publications (1)

Publication Number Publication Date
CN107425521A true CN107425521A (en) 2017-12-01

Family

ID=60427760

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710508974.4A Pending CN107425521A (en) 2017-06-28 2017-06-28 A kind of dynamic aggregation modeling for considering load fluctuation and participation frequency modulation method

Country Status (1)

Country Link
CN (1) CN107425521A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108491983A (en) * 2018-04-11 2018-09-04 西南交通大学 A kind of extensive Unit Combination optimization method of more stack fuel cell electricity generation systems
CN109193719A (en) * 2018-08-03 2019-01-11 中国电力科学研究院有限公司 A kind of modeling method and system for assessing distributed energy storage systematic polymerization frequency modulation performance
CN109245133A (en) * 2018-11-02 2019-01-18 国网电子商务有限公司 Energy storage regulation method, electronic equipment and the computer readable storage medium of micro-capacitance sensor
CN110571793A (en) * 2019-08-23 2019-12-13 华北电力大学 multidimensional identification method for flexible load demand response effect
CN110880772A (en) * 2019-11-08 2020-03-13 武汉大学 Electricity selling company response power grid control method based on industrial park load aggregation
CN110991693A (en) * 2019-10-28 2020-04-10 重庆大学 Load price-quantity curve aggregation system considering uncertainty of electricity purchase cost
CN113297799A (en) * 2021-06-10 2021-08-24 国网综合能源服务集团有限公司 Air conditioner cluster load demand response potential evaluation method based on data driving
CN114268168A (en) * 2021-11-23 2022-04-01 贵州电网有限责任公司 Network security monitoring method applied to power monitoring system
CN115498661A (en) * 2022-11-18 2022-12-20 中国电力科学研究院有限公司 Power system frequency modulation parameter aggregation method and system considering new energy

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105720592A (en) * 2016-04-05 2016-06-29 东南大学 Method for aggregating air-conditioning loads to participate into auxiliary secondary frequency modulation of system
CN106372752A (en) * 2016-08-31 2017-02-01 东南大学 Variable frequency air conditioner thermal battery modeling and scheduling method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105720592A (en) * 2016-04-05 2016-06-29 东南大学 Method for aggregating air-conditioning loads to participate into auxiliary secondary frequency modulation of system
CN106372752A (en) * 2016-08-31 2017-02-01 东南大学 Variable frequency air conditioner thermal battery modeling and scheduling method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张良杰: "《电力系统柔性负荷的建模与控制策略研究》", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
赵丽: "《全局K_均值聚类算法研究与改进》", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108491983A (en) * 2018-04-11 2018-09-04 西南交通大学 A kind of extensive Unit Combination optimization method of more stack fuel cell electricity generation systems
CN108491983B (en) * 2018-04-11 2021-09-03 西南交通大学 Large-scale unit combination optimization method for multi-stack fuel cell power generation system
CN109193719A (en) * 2018-08-03 2019-01-11 中国电力科学研究院有限公司 A kind of modeling method and system for assessing distributed energy storage systematic polymerization frequency modulation performance
CN109193719B (en) * 2018-08-03 2021-11-26 中国电力科学研究院有限公司 Modeling method and system for evaluating aggregation frequency modulation performance of distributed energy storage system
CN109245133A (en) * 2018-11-02 2019-01-18 国网电子商务有限公司 Energy storage regulation method, electronic equipment and the computer readable storage medium of micro-capacitance sensor
CN110571793A (en) * 2019-08-23 2019-12-13 华北电力大学 multidimensional identification method for flexible load demand response effect
CN110991693A (en) * 2019-10-28 2020-04-10 重庆大学 Load price-quantity curve aggregation system considering uncertainty of electricity purchase cost
CN110880772B (en) * 2019-11-08 2023-04-18 武汉大学 Electricity selling company response power grid control method based on industrial park load aggregation
CN110880772A (en) * 2019-11-08 2020-03-13 武汉大学 Electricity selling company response power grid control method based on industrial park load aggregation
CN113297799A (en) * 2021-06-10 2021-08-24 国网综合能源服务集团有限公司 Air conditioner cluster load demand response potential evaluation method based on data driving
CN113297799B (en) * 2021-06-10 2024-02-06 国网综合能源服务集团有限公司 Air conditioner cluster load demand response potential evaluation method based on data driving
CN114268168A (en) * 2021-11-23 2022-04-01 贵州电网有限责任公司 Network security monitoring method applied to power monitoring system
CN115498661A (en) * 2022-11-18 2022-12-20 中国电力科学研究院有限公司 Power system frequency modulation parameter aggregation method and system considering new energy
CN115498661B (en) * 2022-11-18 2023-01-17 中国电力科学研究院有限公司 Power system frequency modulation parameter aggregation method and system considering new energy

Similar Documents

Publication Publication Date Title
CN107425521A (en) A kind of dynamic aggregation modeling for considering load fluctuation and participation frequency modulation method
CN107301472A (en) Distributed photovoltaic planing method based on scene analysis method and voltage-regulation strategy
CN104348188A (en) Distributed generation running and monitoring method
CN104318494A (en) Distributed generation intelligent monitoring system
CN114021390A (en) Random robust optimization method for urban comprehensive energy system and application thereof
CN107069776A (en) A kind of energy storage prediction distributed control method of smooth microgrid dominant eigenvalues
Nykamp et al. Storage operation for peak shaving of distributed PV and wind generation
CN107482618A (en) Electricity pricing method with the generating of wind fire and the power network of flexible load containing polymorphic type
CN114156951B (en) Control optimization method and device of source network load storage system
CN107273968A (en) A kind of Multiobjective Scheduling method and device based on dynamic fuzzy Chaos-Particle Swarm Optimization
CN115102953B (en) Power distribution network Yun Bianduan cooperative management and control system and method
CN109272258A (en) Region wind light generation stock assessment method based on K-means cluster
Ma et al. Optimal SOC control and rule-based energy management strategy for fuel-cell-based hybrid vessel including batteries and supercapacitors
Chen et al. Flexible dispatching method for park‐level integrated energy systems considering energy cascade utilization
CN118040652A (en) Wind-light-energy storage scheduling control method, system, equipment and medium
Li et al. A day-ahead peer to peer energy sharing strategy among energy hubs considering flexibility of energy storage and loads
CN115879715A (en) Response characteristic considered equivalent aggregation method for polymorphic flexible resources of virtual power plant
CN112202202A (en) Wind power plant group coordination control method under multilayer hierarchical structure
CN106786505A (en) A kind of dispersion charging pile coordinated dispatching method based on neighborhood information
Zhang et al. Coordinated dispatch of the wind-thermal power system by optimizing electric vehicle charging
CN105958513B (en) Regional power grid power dynamic balancing method based on HVAC system
CN115907117A (en) Hybrid unmanned aerial vehicle energy management method based on deep learning and dynamic planning
Nandhini et al. Green Energy Management and Environmental Monitoring for Smart Cities
CN106056317B (en) A kind of multiple flexible resource multimode Risk Calculation method
Zhang Research on multi-form flexible resource aggregation model of virtual power plant based on cloud model

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20171201

RJ01 Rejection of invention patent application after publication