CN108805445B - Grouping sequence scheduling method for providing rotary standby for air conditioner load group - Google Patents

Grouping sequence scheduling method for providing rotary standby for air conditioner load group Download PDF

Info

Publication number
CN108805445B
CN108805445B CN201810584434.9A CN201810584434A CN108805445B CN 108805445 B CN108805445 B CN 108805445B CN 201810584434 A CN201810584434 A CN 201810584434A CN 108805445 B CN108805445 B CN 108805445B
Authority
CN
China
Prior art keywords
conditioning load
air conditioning
air
load group
group
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810584434.9A
Other languages
Chinese (zh)
Other versions
CN108805445A (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.)
Zhejiang University ZJU
State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU, State Grid Jiangsu Electric Power Co Ltd filed Critical Zhejiang University ZJU
Priority to CN201810584434.9A priority Critical patent/CN108805445B/en
Publication of CN108805445A publication Critical patent/CN108805445A/en
Application granted granted Critical
Publication of CN108805445B publication Critical patent/CN108805445B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313—Resource planning in a project environment
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2310/00—The network for supplying or distributing electric power characterised by its spatial reach or by the load
    • H02J2310/10—The network having a local or delimited stationary reach
    • H02J2310/12—The local stationary network supplying a household or a building
    • H02J2310/14—The load or loads being home appliances
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
    • Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
    • Y02B70/3225—Demand response systems, e.g. load shedding, peak shaving
    • Y—GENERAL 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y—GENERAL 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
    • Y04S20/20—End-user application control systems
    • Y04S20/222—Demand response systems, e.g. load shedding, peak shaving

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Operations Research (AREA)
  • Physics & Mathematics (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Game Theory and Decision Science (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Air Conditioning Control Device (AREA)

Abstract

The invention discloses a grouping sequence scheduling method for providing rotary standby for an air conditioner load group. Establishing an aggregation model aiming at all air conditioner loads which can be scheduled, dividing all the air conditioner loads into a plurality of groups of air conditioner load groups, and obtaining aggregation power of the air conditioner load groups with subordinate states and subordinate state vectors of the air conditioner load groups; adjusting the set temperature of the air conditioning load at the scheduling moment by using the temperature adjustment value vector to provide rotary standby, and calculating the standby capacity and rebound capacity of the air conditioning load group in the rotary standby process by using the aggregation power; the method adopts a grouping sequence scheduling method to cooperatively optimize the scheduling time of the air conditioning load group and the selection of the air conditioning loads in the group, so that the air conditioning load group can provide rated spare capacity within rated duration. The invention has strong operability, enables the air-conditioning load group to simultaneously meet the requirements of different rated spare capacities and rated duration of the system, flexibly participates in the operation of the system, and has great popularization and application space.

Description

Grouping sequence scheduling method for providing rotary standby for air conditioner load group
Technical Field
The invention relates to a power load grouping scheduling method, in particular to a grouping sequence scheduling method for providing rotary standby for an air conditioner load group.
Background
The addition of new energy power generation brings more uncertainty and unpredictability to the power generation side of the system, and the requirement of the system for balance resources under different time scales is increased. In the traditional sense, balance resources are provided through power generation side resources, and great power generation resource waste is caused.
In recent years, due to the development of the smart home market, the popularization of a communication system and the popularization of a smart electric meter, demand side resources can respond to signals of the system and participate in assisting the balance of an electric power system. Because the coverage rate of the air conditioner in urban residents is high, the power of the air conditioner load can be flexibly and quickly controlled by changing the set temperature, and the air conditioner load group has the potential of greatly improving the rotary standby of the system. However, after sudden changes in the set temperature, the aggregate power of the air conditioning loads will quickly ramp back up after being reduced to a minimum value (power bounce phenomenon), limiting the maximum duration of time that the air conditioning loads provide spinning reserve. The control of the air conditioning load needs to consider the rebound phenomenon of the aggregated power, so that the air conditioning load group can reliably provide standby services with various durations, and the schedulable potential of the air conditioning load resource is fully exerted.
Disclosure of Invention
In view of the problems in the background art, the invention provides a method for scheduling a packet sequence for providing a rotary standby for an air conditioner load group, which considers that the packet sequence for providing the rotary standby under various duration times is scheduled after rebounding the load, can effectively control the rebounding load, and simultaneously meets the requirements of the system standby capacity and the duration time for the packet sequence scheduling.
The technical scheme of the invention is as follows:
1) establishing an aggregation model aiming at all air conditioner loads which can be scheduled, wherein in the aggregation model: dividing all the air-conditioning loads capable of being scheduled into a plurality of air-conditioning load groups, sequencing the air-conditioning load groups in sequence, and establishing the aggregate power P of the air-conditioning load groups with subordinate states through the working model of the single air-conditioning loadgAnd the subordinate state vector S of the air-conditioning load groupg(ii) a And using the temperature adjustment value vector gamma of all air conditioner loads at the scheduling time taugd adjusting the set temperature of the air conditioning load in the air conditioning load group to provide a rotary standby, and utilizing the aggregate power PgCalculating reserve capacity RC of air conditioning load group in process of providing rotary reservegAnd a rebound capacity PLg；
2) The scheduling time of each air-conditioning load group and the selection of the air-conditioning loads in the group are cooperatively optimized by adopting a grouping sequence scheduling method, so that the air-conditioning load groups can provide rated spare capacity within rated duration, and meanwhile, rebound capacity harmful to the system is restrained.
The polymerization power of the air conditioning load group in the step 1) is obtained by adopting the following mode:
1.1) the operating state of the individual air conditioning loads is controlled by a set temperature, the air conditioning loads are operated such that the room air temperature is maintained within an allowable temperature fluctuation range, and the on-off state of the air conditioning loads at time t is set for each air conditioning load using the following formula:
Figure BDA0001689134380000021
wherein m isi(T) represents the on-off state of the ith air conditioning load at time T, Tset,iIndicates the set temperature of the ith air conditioning load,. DELTA.TiRepresenting an allowable indoor temperature fluctuation interval when the ith air conditioner is in work;
taking the cooling mode in summer as an example, the indoor temperature theta is measured at the moment ti(t) when the temperature is higher than the upper limit of the set temperature, the air conditioner enters an on state, mi(t) ═ 1; indoor temperature theta at time ti(t) when the temperature is lower than the set temperature lower limit, the air conditioner enters an off state, mi(t)＝0。
1.2) establishing a working model of a single air-conditioning load of the following formula aiming at each air-conditioning load, and calculating the indoor air temperature corresponding to the air-conditioning load at the t +1 moment according to the switching state at the t moment:
Figure BDA0001689134380000022
wherein, thetai(t +1) represents the indoor air temperature corresponding to the ith air conditioning load at time t + 1; thetai(t) represents the indoor air temperature corresponding to the ith air conditioning load at time t; Δ t represents the simulation step; thetaa(t) represents the outdoor ambient temperature at time t; ciAnd RiRespectively corresponding heat capacity and heat resistance of the room to the ith air-conditioning load; m isi(t) represents the on-off state of the ith air conditioning load at time t; p is a radical ofiRated for the ith air conditioning loadPower; COPiThe energy efficiency coefficient of the ith air conditioner load;
1.3) after obtaining the indoor air temperature corresponding to the air-conditioning load at the next moment, repeating the iteration of the step 1.1) and then calculating the on-off state of the air-conditioning load at the next moment, and further continuously repeating the steps 1.1) -1.2) to perform the iteration calculation to obtain the on-off state of each air-conditioning load;
in specific implementation, for t th of air conditioning load0Time, t-th time is preset0Indoor air temperature theta corresponding to air conditioner load at momenti(t) an initial value in accordance with θi(t) obtaining an initial setting m using a formula of the working modeli(t0) And then the indoor air temperature at each moment t is generated through iterative calculation of the working model.
1.4) after obtaining the on-off state of each air-conditioning load, calculating the aggregate power of the g-th air-conditioning load group by adopting the following formula:
Figure BDA0001689134380000031
wherein, PgIs the aggregate power of the g-th air conditioning load group, g represents the ordinal number of the air conditioning load group, NmaxThe total number of all air-conditioning loads is represented, all the air-conditioning loads which can be scheduled are divided into a plurality of groups of air-conditioning load groups and are sequentially sorted, and the serial number is represented by g; v. ofiState, s, indicating whether the ith air conditioning load can be scheduled within the nominal standby durationg,iA subordinate state indicating whether or not the i-th air conditioning load belongs to the g-th air conditioning load group (s if the i-th air conditioning load belongs to the g-th air conditioning load group g,i1 if not, 0);
and constructing a subordinate state vector S of all following air conditioning loads to the g-th air conditioning load groupgExpressed as:
Figure BDA0001689134380000032
wherein s isg,1,sg,2,…,
Figure BDA0001689134380000033
The subordinate states of the 1 st to the last air conditioning loads to the g-th air conditioning load group are shown.
In step 1.4), the present invention defines the air conditioning load capable of being scheduled in the period of the rated standby as the air conditioning load which is always in the on state in the rated standby duration, and the status v of whether the ith air conditioning load can be scheduled in the rated standby durationiThe following formula is used for calculation:
Figure BDA0001689134380000034
wherein v isiA status indicating whether the ith air conditioning load can be scheduled within the rated standby duration,
Figure BDA0001689134380000035
indicating the customary off-time, t, of the ith air conditioning loadinsIndicating the starting time of the air conditioning load group to provide a spinning reserve, RC*Indicating rated reserve capacity, DT*Indicating nominal standby duration, Onoffi(tins) Indicates the start time t of the ith air conditioning load in providing the spinning reserveinsIf the ith air conditioning load is at the start time t of providing the spinning reserveinsIs in an on state, then Onoff is generatedi(tins) 1, and 0 if off).
As can be seen from the above equation, if the ith air conditioning load can be scheduled within the rated standby duration, then v i1, if the ith air conditioning load cannot be scheduled within the rated standby duration, viIs 0.
The g-th air-conditioning load group provides rotary standby by adjusting the set temperature of the air-conditioning load, if the set temperature is adjusted, the spare capacity RC of the air-conditioning load group in the rotary standby process is calculatedgAnd a rebound capacity PLgAs a backup characteristic and a rebound characteristic.
In step 1.1), the set temperature T of the ith air-conditioning loadset,iThe adjustment may occur during the provision of the spinning reserve, in particular by adjusting the set temperature T of the air conditioning load using the following formulaset,iTo provide rotational redundancy:
Figure BDA0001689134380000041
wherein, Tset,iIndicates the set temperature of the i-th air conditioning load,
Figure BDA0001689134380000042
scheduling time, s, for providing rotary standby for the g-th air conditioning load groupg,iA subordinate status, γ, indicating whether the ith air-conditioning load belongs to the g-th air-conditioning load groupiSetting a temperature adjustment value for the ith air conditioner load at the scheduling time;
if the temperature T is setset,iIf no adjustment is made, the reserve capacity RC of the g-th air conditioning load group is setgAnd a rebound capacity PLgIs zero.
If the temperature T is setset,iIf the adjustment occurs, the spare capacity RC of the g-th air conditioning load group is calculated by adopting the following formulagAnd a rebound capacity PLg；
In the invention, the g-th air conditioner load group is set
Figure BDA0001689134380000043
The reduction of the aggregate power after the change of the set temperature at each moment corresponds to the reserve capacity provided, the reserve capacity RC of the g-th air conditioning load groupgThe following formula is used for calculation:
Figure BDA0001689134380000044
Figure BDA0001689134380000045
Figure BDA0001689134380000046
wherein,
Figure BDA0001689134380000047
is shown in
Figure BDA0001689134380000048
The average value of the aggregated power of the g-th air conditioning load group before the scheduling time,
Figure BDA0001689134380000049
representing the aggregate power minimum, τ, of the g-th air conditioning load group during the provision of the spinning reserve periodg rIndicating the end time of the air conditioning load group for providing the rotation standby;
the aggregate power of the air conditioning load group reaches the lowest value
Figure BDA00016891343800000410
The operating characteristics of the air conditioning loads then cause aggregate power bounce, the bounce capacity PL of the g-th air conditioning load groupgThe following formula is used for calculation:
Figure BDA00016891343800000411
Figure BDA00016891343800000412
Figure BDA00016891343800000413
wherein,
Figure BDA00016891343800000414
represents the maximum value of the aggregation power which the air conditioning load group rebounds to after the aggregation power is reduced to the lowest point,
Figure BDA00016891343800000415
representing aggregate power P of air conditioning load groupsgAchieved by providing a spinning reserve period
Figure BDA00016891343800000416
The time of day.
In the step 2), the method specifically comprises the following steps:
2.1) aiming at the air conditioner load groups, adopting the standby scheduling control of a batch sequential scheduling strategy response system, and carrying out double-layer optimization aiming at the selection and scheduling time of the air conditioner load in each air conditioner load group, specifically constructing the following two-layer optimization objective function:
the rotary standby scheduling signal issued by the system comprises three parts: starting time t for providing spinning reserveinsRated spare capacity RC*Rated standby duration DT*. The invention optimizes the scheduling time of the g-th air conditioning load group as the rotary standby, so that the standby capacity and the rated standby capacity RC actually provided by the air conditioning load group*And (5) constructing an upper-layer optimization objective function expressed by the following formula when the difference is minimum:
Figure BDA0001689134380000051
wherein, RC*Indicating the rated reserve capacity actually provided by the air conditioning load group; pjRepresents the aggregate power of the jth air conditioning load group,
Figure BDA0001689134380000052
representing the average value of the aggregated power of the jth air-conditioning load group before the scheduling time;
Figure BDA0001689134380000053
representing aggregate power P of air conditioning load groupsgAchieved by providing a spinning reserve period
Figure BDA0001689134380000054
Time of (t)insA start time indicating that the air conditioning load needs to provide the rated reserve capacity;
and simultaneously constructing a first constraint condition of an upper-layer optimization objective function: the scheduling time meeting the g-th air conditioning load group is after the scheduling time of the g-1 st air conditioning load group, namely:
Figure BDA0001689134380000055
and simultaneously constructing a second constraint condition of the upper-layer optimization objective function: the scheduling time limit of the g-th air conditioning load group is met within the rated duration, namely:
Figure BDA0001689134380000056
the invention optimizes the subordinate state vector S of the air conditioning load in the g-th air conditioning load group to avoid the waste of the spare capacity of the air conditioning load groupgAnd the actual temperature adjustment value gamma of the air conditioning load so that
Figure BDA00016891343800000517
And the system spare capacity requirement can be met at the minimum time, and a lower-layer optimization objective function represented by the following formula is formed:
Figure BDA0001689134380000057
wherein,
Figure BDA0001689134380000058
indicating the maximum spare capacity of the air conditioning load group when the temperature adjustment value gamma of all the air conditioning loads reaches the maximum temperature adjustment amount
Figure BDA0001689134380000059
Timely g-th air conditionerSpare capacity RC reached by load groupgAs maximum spare capacity
Figure BDA00016891343800000510
Maximum temperature adjustment
Figure BDA00016891343800000511
Is a preset value;
and simultaneously constructing a first constraint condition of a lower-layer optimization objective function: spare capacity RC satisfying g-th air conditioner load groupgEqual to the bounce capacity of the g-1 st air conditioning load group, i.e.:
RCg＝PLg-1
and simultaneously constructing a second constraint condition of the lower-layer optimization objective function: the actual temperature adjustment range meeting the air conditioning load is within the preset maximum temperature adjustment amount, namely:
Figure BDA00016891343800000512
Figure BDA00016891343800000513
Figure BDA00016891343800000514
where γ represents a vector of temperature adjustment values for all air conditioning loads, γ1,γ2,…,
Figure BDA00016891343800000515
Represents a maximum temperature adjustment value corresponding to the 1 st air conditioning load to the last air conditioning load,
Figure BDA00016891343800000516
a maximum temperature adjustment value vector representing all air conditioning loads,
Figure BDA0001689134380000061
indicates a maximum temperature adjustment value corresponding to the ith air conditioning load,
Figure BDA0001689134380000062
the maximum temperature adjustment value corresponding to the 1 st air conditioning load to the last air conditioning load is represented;
and simultaneously constructing a third constraint condition of the lower-layer optimization objective function: the air conditioning load satisfying the g-th air conditioning load group is limited in the range of all schedulable air conditioning loads, namely:
Figure BDA0001689134380000063
wherein s isj,iA subordinate state v indicating whether or not the ith air conditioning load belongs to the jth air conditioning load groupiA status indicating whether an ith air conditioning load can be scheduled within a rated standby duration;
2.2) scheduling time of the g-th air conditioning load group for providing rotary standby
Figure BDA0001689134380000064
The set temperature T for adjusting the air conditioning load by providing a spinning reserve is affectedset,iAnd further influences the aggregate power P of the g-th air conditioning load groupgThereby forming influence on the upper and lower layer optimization objective function; dependent state vector SgWill affect the aggregate power P of the g-th air conditioning load groupgAnd further influences the reserve capacity RC of the g-th air conditioning load groupgThereby forming influence on the upper and lower layer optimization objective function; the temperature adjustment value vector gamma influences the set temperature T for providing a spinning reserve to adjust the air conditioning loadset,iAnd further influences the aggregate power P of the g-th air conditioning load groupgThereby having an influence on the upper and lower optimization objective functions.
Iterative solution is carried out on the upper-layer optimization objective function and the lower-layer optimization objective function in a feasible space by adopting a genetic algorithm to obtain an optimal dependent state vector SgVector of temperature adjustment value γ and individual air conditioning loadScheduling time
Figure BDA0001689134380000065
Optimizing the target functions of the upper layer and the lower layer simultaneously;
2.3) Using the dependency State vector SgTemperature adjustment value vector gamma and scheduling time of each air conditioner load
Figure BDA0001689134380000066
All air conditioning loads are distributed and scheduled such that the air conditioning load groups are able to provide a rated backup capacity for a rated duration.
The invention has the following beneficial effects:
1) the invention can effectively solve the limitation of the rebound load phenomenon on the duration time of the rotation reserve provided by the air conditioning load group, and realize the flexible control of the reserve duration time, thereby enabling the air conditioning load group to provide the rotation reserve under various duration times.
3) The invention optimizes the load selection and the temperature adjustment range of each air-conditioning load group, so that the potential of all air-conditioning loads for providing rotary reserve can be fully utilized, the air-conditioning load groups can simultaneously meet the requirements of different rated reserve capacities and rated duration of the system, the air-conditioning load groups can flexibly participate in the operation of the system, and the invention has great popularization and application space.
Drawings
FIG. 1 is a flow chart of the algorithm of the present invention.
Fig. 2 is a schematic diagram of an air conditioner load grouping sequential scheduling method.
Fig. 3 is a diagram of the control results of the air conditioning load groups providing spinning reserve according to the group sequential scheduling method.
Detailed Description
The following description is made with reference to the embodiments and the accompanying drawings.
As shown in fig. 1, the method of the present invention comprises:
1) establishing an aggregation model aiming at all air conditioner loads which can be scheduled, wherein in the aggregation model: dividing all the air-conditioning loads which can be dispatched into a plurality of air-conditioning load groups and carrying out smoothingSequence ordering, building up aggregate power P of air conditioning load group with subordinate state by working model of single air conditioning loadgAnd the subordinate state vector S of the air-conditioning load groupg(ii) a And using the temperature adjustment value vector gamma of all air conditioner loads at the scheduling time
Figure BDA0001689134380000071
Adjusting the set temperature of the air conditioning loads in the air conditioning load group to provide a rotary standby, using the aggregate power PgCalculating reserve capacity RC of air conditioning load group in process of providing rotary reservegAnd a rebound capacity PLg；
2) The scheduling time of each air-conditioning load group and the selection of the air-conditioning loads in the group are cooperatively optimized by adopting a grouping sequence scheduling method, so that the air-conditioning load groups can provide rated spare capacity within rated duration, and meanwhile, rebound capacity harmful to the system is restrained.
The principle of air conditioner load group sequential scheduling against bounce loads is shown in fig. 2. FIG. 2 shows
Figure BDA0001689134380000072
Power curve of g-1 air conditioning load group scheduled at time, an
Figure BDA0001689134380000073
The power curve of the g-th air conditioning load group scheduled at the moment. The g-1 air conditioning load is set in
Figure BDA0001689134380000074
After the set temperature is adjusted constantly, the polymerization power is reduced, and the spare capacity RC is provided for the systemg-1With spare capacity at RCg-1After a certain period of time, the polymerization power is quickly increased due to the rebound of the polymerization power, and after the polymerization power is increased to a certain degree, the polymerization power is stable and shows PLg-1The bounce capacity. Spare capacity RC of air conditioning load group of the latter group (g-th group)gRebound capacity PL of air conditioning load group of the preceding group (k-1 th group)g-1When equal, the rebound capacity of the former air conditioning load group can be controlled by the work of the latter air conditioning load groupThe rate reduction cancels out as shown on the right side of fig. 2. According to this principle, the air conditioning load groups are scheduled in sequence so that the bounce capacity can be cancelled in sequence until the duration reaches the rated value. Thus, the packet sequential scheduling enables the air conditioning load group to provide the rated reserve capacity for the rated duration.
The examples of the implementation of the complete process according to the invention are as follows:
the following description is made by using simulation example:
assuming that 60000 air-conditioning loads can be scheduled for a residential area, all operating in the summer cooling mode, the air-conditioning load parameters are set as follows:
TABLE 1 ith air conditioner load parameter settings
Figure BDA0001689134380000075
Figure BDA0001689134380000081
Assuming that the system standby demand occurs at 16:00 and the air conditioning loads are required to provide 5MW of standby capacity in a period of 16:00 to 18:00, the air conditioning load groups are sequentially scheduled according to the calculation result of the grouping sequential scheduling method, and the power curve of each air conditioning load group and the total power curve of the superposition of all the air conditioning load groups are shown in fig. 3.
Looking at the power curve (e.g., first) for each air conditioning load group alone, it can be seen that the aggregate power bounces off after being clipped to a minimum value, resulting in a reserve capacity of only about 10 minutes in duration. The packet sequential scheduling causes the power reduction of the next air conditioning load group to in turn compensate for the bounce capacity of the previous air conditioning load group. After scheduling to the fifth group of air conditioning loads, the bounce capacity has been reduced to an extremely low level, which is approximately negligible.
Therefore, the method can effectively reduce rebound load and realize long-time power reduction, thereby providing a rotary standby with longer duration and realizing the outstanding technical effect.

Claims (5)

1. A grouping sequence scheduling method for providing rotary standby for an air conditioning load group is characterized by comprising the following steps:
1) establishing an aggregation model aiming at all air conditioner loads which can be scheduled, wherein in the aggregation model: dividing all air conditioning loads capable of being scheduled into a plurality of air conditioning load groups, and establishing the aggregate power P of the air conditioning load groups with subordinate states through the working model of the single air conditioning loadgAnd the subordinate state vector S of the air-conditioning load groupg(ii) a And using the temperature adjustment value vector gamma of all air conditioner loads at the scheduling time
Figure FDA0003023664070000011
Adjusting the set temperature of the air conditioning loads in the air conditioning load group to provide a rotary standby, using the aggregate power PgCalculating reserve capacity RC of air conditioning load group in process of providing rotary reservegAnd a rebound capacity PLg；
2) And the scheduling time of each air-conditioning load group and the selection of the air-conditioning loads in the group are cooperatively optimized by adopting a grouping sequence scheduling method, so that the air-conditioning load groups can provide rated spare capacity within rated duration.
2. The method of claim 1, wherein the method comprises: the polymerization power of the air conditioning load group in the step 1) is obtained by adopting the following mode:
1.1) setting the on-off state of the air-conditioning load at time t for each air-conditioning load using the following formula:
Figure FDA0003023664070000012
wherein m isi(T) represents the on-off state of the ith air conditioning load at time T, Tset,iIndicating the set temperature, Δ T, of the ith air conditioning loadiIndicates the ith nullThe allowable indoor temperature fluctuation interval during load regulation work; thetai(t) represents the indoor air temperature corresponding to the ith air conditioning load at time t;
1.2) establishing a working model of a single air-conditioning load of the following formula aiming at each air-conditioning load, and calculating the indoor air temperature corresponding to the air-conditioning load at the t +1 moment according to the switching state at the t moment:
Figure FDA0003023664070000013
wherein, thetai(t +1) represents the indoor air temperature corresponding to the ith air conditioning load at time t + 1; Δ t represents the simulation step; thetaa(t) represents the outdoor ambient temperature at time t; ciAnd RiRespectively corresponding heat capacity and heat resistance of the room to the ith air-conditioning load; m isi(t) represents the on-off state of the ith air conditioning load at time t; p is a radical ofiRated power for the ith air conditioning load; COPiThe energy efficiency coefficient of the ith air conditioner load;
1.3) after obtaining the indoor air temperature corresponding to the air-conditioning load at the next moment, repeating the iteration of the step 1.1) and then calculating the on-off state of the air-conditioning load at the next moment, and further continuously repeating the steps 1.1) -1.2) to perform the iteration calculation to obtain the on-off state of each air-conditioning load;
1.4) after obtaining the on-off state of each air-conditioning load, calculating the aggregate power of the g-th air-conditioning load group by adopting the following formula:
Figure FDA0003023664070000021
wherein, PgIs the aggregate power of the g-th air conditioning load group, g represents the ordinal number of the air conditioning load group, NmaxIndicating the total number of all air conditioning loads, viState, s, indicating whether the ith air conditioning load can be scheduled within the nominal standby durationg,iA subordinate state indicating whether or not the ith air conditioning load belongs to the g-th air conditioning load group;
and constructing a subordinate state vector S of all following air conditioning loads to the g-th air conditioning load groupgExpressed as:
Figure FDA0003023664070000022
wherein,
Figure FDA0003023664070000023
the subordinate states of the 1 st to the last air conditioning loads to the g-th air conditioning load group are shown.
3. The method of claim 2, wherein the method comprises: state v in step 1.4) whether the ith air conditioning load can be scheduled within the rated standby durationiThe following formula is used for calculation:
Figure FDA0003023664070000024
wherein v isiA status indicating whether the ith air conditioning load can be scheduled within the rated standby duration,
Figure FDA0003023664070000025
indicating the customary off-time, t, of the ith air conditioning loadinsIndicating the starting time of the air conditioning load group to provide a spinning reserve, RC*Indicating rated reserve capacity, DT*Indicating nominal standby duration, Onoffi(tins) Indicates the start time t of the ith air conditioning load in providing the spinning reserveinsThe switch state of (1).
4. The method of claim 1, wherein the method comprises: in the step 1.1), the regulation is carried out by adopting the following formulaSet temperature T of whole air conditioner loadset,iTo provide rotational redundancy:
Figure FDA0003023664070000026
wherein, Tset,iIndicates the set temperature of the i-th air conditioning load,
Figure FDA0003023664070000027
scheduling time, s, for providing rotary standby for the g-th air conditioning load groupg,iA subordinate status, γ, indicating whether the ith air-conditioning load belongs to the g-th air-conditioning load groupiSetting a temperature adjustment value for the ith air conditioner load at the scheduling time;
if the temperature T is setset,iIf no adjustment is made, the reserve capacity RC of the g-th air conditioning load group is setgAnd a rebound capacity PLgIs zero;
if the temperature T is setset,iIf the adjustment occurs, the spare capacity RC of the g-th air conditioning load group is calculated by adopting the following formulagAnd a rebound capacity PLg；
Spare capacity RC of g-th air conditioning load groupgThe following formula is used for calculation:
Figure FDA0003023664070000031
Figure FDA0003023664070000032
Figure FDA0003023664070000033
wherein,
Figure FDA0003023664070000034
is shown in
Figure FDA0003023664070000035
The average value of the aggregated power of the g-th air conditioning load group before the scheduling time,
Figure FDA0003023664070000036
represents the aggregate power minimum value of the g-th air conditioning load group during the spinning standby period,
Figure FDA0003023664070000037
indicating the end time of the air conditioning load group for providing the rotation standby;
rebound capacity PL of the g-th air conditioning load groupgThe following formula is used for calculation:
Figure FDA0003023664070000038
Figure FDA0003023664070000039
Figure FDA00030236640700000310
wherein,
Figure FDA00030236640700000311
represents the maximum value of the aggregation power which the air conditioning load group rebounds to after the aggregation power is reduced to the lowest point,
Figure FDA00030236640700000312
representing aggregate power P of air conditioning load groupsgAchieved by providing a spinning reserve period
Figure FDA00030236640700000313
The time of day.
5. The method of claim 1, wherein the method comprises: in the step 2), the method specifically comprises the following steps:
2.1) carrying out double-layer optimization aiming at the selection and scheduling time of the air-conditioning loads in each air-conditioning load group, specifically, constructing the following two-layer optimization objective function:
constructing an upper-layer optimization objective function expressed by the following formula:
Figure FDA00030236640700000314
wherein, RC*Indicating the rated reserve capacity actually provided by the air conditioning load group; pjRepresents the aggregate power of the jth air conditioning load group,
Figure FDA00030236640700000315
representing the average value of the aggregated power of the jth air-conditioning load group before the scheduling time;
Figure FDA00030236640700000316
representing aggregate power P of air conditioning load groupsgAchieved by providing a spinning reserve period
Figure FDA00030236640700000317
Time of (t)insA start time indicating that the air conditioning load needs to provide the rated reserve capacity;
and simultaneously constructing a first constraint condition of an upper-layer optimization objective function: the scheduling time meeting the g-th air conditioning load group is after the scheduling time of the g-1 st air conditioning load group, namely:
Figure FDA00030236640700000318
and simultaneously constructing a second constraint condition of the upper-layer optimization objective function: the scheduling time limit of the g-th air conditioning load group is met within the rated duration, namely:
Figure FDA0003023664070000041
the lower layer optimization objective function represented by the following formula is formed:
Figure FDA0003023664070000042
wherein,
Figure FDA0003023664070000043
indicating the maximum spare capacity of the air conditioning load group when the temperature adjustment value gamma of all the air conditioning loads reaches the maximum temperature adjustment amount
Figure FDA0003023664070000044
Reserve capacity RC reached by the g-th air conditioning load groupgAs maximum spare capacity
Figure FDA0003023664070000045
And simultaneously constructing a first constraint condition of a lower-layer optimization objective function: spare capacity RC satisfying g-th air conditioner load groupgEqual to the bounce capacity of the g-1 st air conditioning load group, i.e.:
RCg＝PLg-1
and simultaneously constructing a second constraint condition of the lower-layer optimization objective function: the actual temperature adjustment range meeting the air conditioning load is within the preset maximum temperature adjustment amount, namely:
Figure FDA0003023664070000046
Figure FDA0003023664070000047
Figure FDA0003023664070000048
where γ represents a vector of temperature adjustment values for all air conditioning loads,
Figure FDA0003023664070000049
represents a maximum temperature adjustment value corresponding to the 1 st air conditioning load to the last air conditioning load,
Figure FDA00030236640700000410
a maximum temperature adjustment value vector representing all air conditioning loads,
Figure FDA00030236640700000411
indicates a maximum temperature adjustment value corresponding to the ith air conditioning load,
Figure FDA00030236640700000412
the maximum temperature adjustment value corresponding to the 1 st air conditioning load to the last air conditioning load is represented;
and simultaneously constructing a third constraint condition of the lower-layer optimization objective function: the air conditioning load satisfying the g-th air conditioning load group is limited in the range of all schedulable air conditioning loads, namely:
Figure FDA00030236640700000413
wherein s isj,iA subordinate state v indicating whether or not the ith air conditioning load belongs to the jth air conditioning load groupiA status indicating whether an ith air conditioning load can be scheduled within a rated standby duration;
2.2) adopting a genetic algorithm to carry out iterative solution on an upper-layer optimized objective function and a lower-layer optimized objective function in a feasible space to obtain the best solutionPreferred dependent state vector SgTemperature adjustment value vector gamma and scheduling time of each air conditioner load
Figure FDA00030236640700000414
Optimizing the target functions of the upper layer and the lower layer simultaneously;
2.3) Using the dependency State vector SgTemperature adjustment value vector gamma and scheduling time of each air conditioner load
Figure FDA00030236640700000415
All air conditioning loads are distributed and scheduled such that the air conditioning load groups are able to provide a rated backup capacity for a rated duration.
CN201810584434.9A 2018-06-08 2018-06-08 Grouping sequence scheduling method for providing rotary standby for air conditioner load group Active CN108805445B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810584434.9A CN108805445B (en) 2018-06-08 2018-06-08 Grouping sequence scheduling method for providing rotary standby for air conditioner load group

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810584434.9A CN108805445B (en) 2018-06-08 2018-06-08 Grouping sequence scheduling method for providing rotary standby for air conditioner load group

Publications (2)

Publication Number Publication Date
CN108805445A CN108805445A (en) 2018-11-13
CN108805445B true CN108805445B (en) 2021-07-06

Family

ID=64087677

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810584434.9A Active CN108805445B (en) 2018-06-08 2018-06-08 Grouping sequence scheduling method for providing rotary standby for air conditioner load group

Country Status (1)

Country Link
CN (1) CN108805445B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111271839B (en) * 2020-02-14 2021-05-11 河海大学 Method for adjusting short-term power of fixed-frequency air conditioner
CN115377963B (en) * 2022-07-18 2023-05-09 四川大学 Novel temperature set value control method for heterogeneous cluster air conditioner load

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104214912A (en) * 2014-09-24 2014-12-17 东南大学 Aggregation air conditioning load scheduling method based on temperature set value adjustment
CN106655200A (en) * 2016-09-13 2017-05-10 浙江大学 Method for calculating electric power system operation reserve response quantity provided by air-conditioner aggregation
CN106921159A (en) * 2017-02-21 2017-07-04 国网山东省电力公司电力科学研究院 A kind of active power modulator approach of many air conditioner load groups

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9632666B2 (en) * 2012-02-16 2017-04-25 GM Global Technology Operations LLC Team-oriented HVAC system
US20170018924A1 (en) * 2015-07-13 2017-01-19 Power Over Time, Inc. Systems and methods for reducing an electric utility reserve capacity using instrumented energy consuming devices

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104214912A (en) * 2014-09-24 2014-12-17 东南大学 Aggregation air conditioning load scheduling method based on temperature set value adjustment
CN106655200A (en) * 2016-09-13 2017-05-10 浙江大学 Method for calculating electric power system operation reserve response quantity provided by air-conditioner aggregation
CN106921159A (en) * 2017-02-21 2017-07-04 国网山东省电力公司电力科学研究院 A kind of active power modulator approach of many air conditioner load groups

Also Published As

Publication number Publication date
CN108805445A (en) 2018-11-13

Similar Documents

Publication Publication Date Title
CN104214912B (en) Aggregation air conditioning load scheduling method based on temperature set value adjustment
CN110492498B (en) Temperature control load participation power system primary frequency modulation method based on double-layer control
CN111555304B (en) Air conditioner load virtual energy storage scheduling method for power grid frequency modulation service
CN107732897B (en) Building microgrid model prediction regulation and control method fusing virtual energy storage system
CN103199555B (en) Control method of secondary frequency modulation of electrical power system with participation of load side resources
CN108036468B (en) Aggregation control method for air conditioning system
CN107631428A (en) A kind of air conditioner temperature/humidity control method and air conditioner
CN108494012B (en) Online optimization method for regional comprehensive energy system considering electricity-to-gas technology
CN110661253B (en) Multi-type electric power elastic load combined adjusting method in building group
CN108805445B (en) Grouping sequence scheduling method for providing rotary standby for air conditioner load group
CN109812946B (en) Control method suitable for large-scale resident air conditioner load group demand response
CN114512994A (en) Frequency modulation method, system, equipment and medium for cluster temperature control load system
CN107228454A (en) Method and device for airconditioning control
CN110165665A (en) A kind of source-lotus-storage dispatching method based on improvement multi-objective particle swarm algorithm
CN113757931A (en) Air conditioner control method and system
CN113028604A (en) Temperature regulation control and scheduling method based on aggregated air conditioner
WO2020055411A1 (en) Systems and methods for rapidly responding to commanded power profiles
CN116085953A (en) Energy-saving control method, system and medium based on dynamic air conditioner operation data
CN114498651B (en) Hybrid load cluster control method and system
CN112054505B (en) Building electric power demand elasticity rapid quantification method based on building basic information
CN114936529A (en) Temperature control load group aggregation model, modeling method and temperature control load group adjustable potential evaluation method
CN110912163B (en) Method for acquiring frequency control strategy of temperature-control load energy efficiency power plant
CN109190988B (en) Demand side response game method for realizing optimal coordination of temperature control load
CN112967152B (en) Water-loop heat pump load aggregation modeling and control method
CN113285463B (en) Controllable load multi-target control method based on sequencing

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