CN112968445B - Multifunctional complementary demand response control method for residential building - Google Patents

Multifunctional complementary demand response control method for residential building Download PDF

Info

Publication number
CN112968445B
CN112968445B CN202110226687.0A CN202110226687A CN112968445B CN 112968445 B CN112968445 B CN 112968445B CN 202110226687 A CN202110226687 A CN 202110226687A CN 112968445 B CN112968445 B CN 112968445B
Authority
CN
China
Prior art keywords
day
load
temperature
demand response
time
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
CN202110226687.0A
Other languages
Chinese (zh)
Other versions
CN112968445A (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.)
Tongji University
Original Assignee
Tongji University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tongji University filed Critical Tongji University
Priority to CN202110226687.0A priority Critical patent/CN112968445B/en
Publication of CN112968445A publication Critical patent/CN112968445A/en
Application granted granted Critical
Publication of CN112968445B publication Critical patent/CN112968445B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/12Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/14Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
    • H02J3/144Demand-response operation of the power transmission or distribution network
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02BCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B70/00Technologies for an efficient end-user side electric power management and consumption
    • Y02B70/30Systems 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/3225Demand response systems, e.g. load shedding, peak shaving
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S20/00Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
    • Y04S20/20End-user application control systems
    • Y04S20/222Demand response systems, e.g. load shedding, peak shaving

Abstract

The invention relates to a multifunctional complementary demand response control method for residential buildings, which selects corresponding weights and reference loads according to different types of control targets of the residential buildings, such as maximum economy, maximum peak reduction and the like, and finally realizes the control of demand response through the optimization of a uniform objective function. Compared with the traditional demand response control strategy method, the method adopts a uniform optimization objective function, and can improve the universality of the demand response control method.

Description

Multifunctional complementary demand response control method for residential building
Technical Field
The invention relates to the technical field of demand response, in particular to a multifunctional complementary demand response control method for a residential building.
Background
In recent years, with the continuous promotion of domestic industrialization and urbanization construction, the problem of imbalance between the supply side and the demand side of a power grid becomes more and more serious. According to statistics, the generating and loading machines in 2018 of China reach 19 hundred million kW and are located in the first world, but the average utilization hours of generating equipment in the whole country in 2017 are only 3670 hours, which means that more than half of power plants are not fully utilized all the year round. In addition, widespread use of renewable energy sources further exacerbates the problem of imbalance on the supply side and the demand side.
Demand Response (DR) alleviates the problem of power supply and demand imbalance during peak periods of electricity by curtailing or shifting peak loads on the customer side. For the supply side, the demand response can slow down newly added peak regulation equipment, further improves the utilization hours of the existing equipment, and reduces the operation and maintenance cost; meanwhile, the land resources originally used for newly-built power plants and power grids are reserved. For the user on the demand side, a corresponding subsidy or reward may generally be available by participating in the demand response.
Previous demand response control strategies require different objective functions for different demand response objectives, such as maximum load shedding or minimum cost, which is not conducive to the popularization and implementation of demand response. Therefore, it is necessary to develop a set of general residential building demand response control strategies, and to perform optimization control by using the same optimization objective function for different demand response objectives.
Disclosure of Invention
The invention aims to solve the technical problem of overcoming the defects in the prior art and provides a demand response control method suitable for residential buildings.
In order to solve the technical problem, the solution of the invention is as follows:
a multifunctional complementary demand response control method for residential buildings can meet the maximum economic benefit (charging according to time-of-use electricity price), the maximum reduction at peak time (subsidy according to the maximum reduction at peak time), the most stable reduction all day, or the maximum reduction all day (namely the maximum energy-saving control). If the actual net electrical load in one day is X ═ X (X) 1 ,x 2 ,…,x n ) Then, thenThe optimal optimization objective function in the whole time of a day is as follows:
Figure BDA0002956527940000021
wherein x is i Is the actual net electrical load in a day, X ═ X 1 ,x 2 ,…,x n ) The value corresponding to the middle time i; y is i For reference, the electrical load Y ═ Y 1 ,y 2 ,…,y n ) The value corresponding to the middle time i; dist (x) i -y i ) Denotes x i And y i A distance function between; w is a i Is the weight corresponding to time i.
Specifically, the values of the terms in the objective function may be determined and solved by:
step 1: dividing the actual net electric load X into basic loads X according to the usage and the use property ba Weather sensitive load X we Transferable type load X sh And a convertible load X ch The method comprises the following steps of totally four types of loads, and determining values and constraint conditions through a calculation formula of each type of load:
X=X ba +X we +X sh +X ch
step 2: according to different demand response control purposes, selecting corresponding one-day reference electric load Y (Y) 1 ,y 2 ,…,y n );
And step 3: selecting corresponding distance function dist (x) according to different demand response control purposes i -y i ) And a weight w i
And 4, step 4: the parameters (start-stop state, temperature setting and the like) related to the equipment obtained by optimization are output as control signals.
The base load X in the step 1 ba =(x ba,1 ,x ba,2 ,...,x ba,n ) In order to maintain the electrical load necessary for residential building use, which is not flexible in the time dimension, such as room LED lighting, desktop outlet power, etc., the value of a day is available according to the daily usage habits of the user.
The weather sensitive load X in the step 1 we =(x we,1 ,x we,2 ,...,x we,n ) The air conditioner is mainly used for heating or refrigerating electric loads of residential buildings, and the air conditioner is influenced by outdoor weather and indoor set temperature.
X we =f 1 (T dry ,T dew ,/,P,T set ,S hvac )
Wherein f is 1 Is X we Load prediction model of (1), T dry =(t dry,1 ,t dry,2 ,…,t dry,n ) The temperature of the dry bulb outside the whole day room; t is dew =(t dew,1 ,t dew,2 ,…,t dew,n ) Is the whole day outdoor wet bulb temperature; j ═ J (J) 1 ,j 2 ,…,j n ) Direct solar radiation outside the whole sky; p ═ P (P) 1 ,p 2 ,…,p n ) Is the atmospheric pressure outside the whole day room; t is set =(t set,1 ,t set,2 ,…,t set,n ) Is the room temperature of the whole daySetting a value; s set =(s set,1 ,s set,2 ,…,s set,n ) The air conditioning system is in an open state all day long.
The load prediction model f 1 Modeling calculation can be performed on an air conditioning system such as a residential building through a physical modeling method such as EnergyPlus, TRNSYS, Modelica and the like, or prediction calculation can be performed on the basis of historical data through a data-driven model such as an Artificial Neural Network (ANN), a Support Vector Machine (SVM), an extreme learning machine (XGboost) and the like.
The temperature T of the whole day outdoor dry bulb dry All day outdoor wet bulb temperature T dew The total day outdoor direct solar radiation J and the total day outdoor atmospheric pressure P are meteorological parameter adjustments required for predicting the load of the air conditioning system through physical model calculation or a data driving model, and can be obtained through weather forecast data of a local meteorological station.
The set value T of the all-day room temperature set The following four constraints should be satisfied:
t set,min ≤t set,i ≤t set,max
Figure BDA0002956527940000031
Figure BDA0002956527940000041
Figure BDA0002956527940000042
the first constraint condition indicates that the indoor temperature setting value should be within the minimum setting value and the maximum setting value, and the setting range may be slightly larger than the human body comfort temperature range provided by ASHRAE due to the requirement of implementing the demand response. The second constraint indicates that the indoor temperature setpoint should be the minimum temperature step t step Integer multiple of this, the setting of this should match the control accuracy of the air conditioning system. Third andthe fourth constraint represents T set =(t set,1 ,t set,2 ,…,t set,n ) Of any arbitrary consecutive a indoor temperature settings, only up to 1 temperature set point adjustment is allowed, which is to prevent the set temperature from being adjusted too frequently.
All-day air conditioning system opening state S set The following two constraints should be satisfied:
s hvac,i ∈{0,1}
Figure BDA0002956527940000043
wherein, the first constraint condition indicates that the on state of the air conditioning system can be off (0) or off (1). The second constraint denotes S hvac =(s hvac,1 ,s hvac,2 ,…,s hvac,n ) In the open states of any continuous beta air conditioning systems, only the adjustment of the open states of the air conditioning systems for at most 1 time is allowed, and the adjustment is set to prevent the air conditioning systems from being started and stopped too frequently.
Transferable load X in said step 1 sh =(x sh,1 ,x sh,2 ,...,x sh,n ) The electric loads with elasticity in the time dimension during daily use of residential buildings, such as dryers, dishwashers and the like. X sh Can be seen as a superposition of multiple household electrical appliance transferable loads:
X sh =X sh1 +X sh2 +…
X sh1 =(x sh1,1 ,x sh1,2 ,...,x sh1,n )
X sh2 =(x sh2,1 ,x sh2,2 ,...,x sh2,n )
……
x sh,1 =x sh1,1 +x sh2,1 +...
x sh,2 =x sh1,2 +x sh2,2 +…
……
with one of the household appliances X sh1 For example, the operating constant power is p sh1 The starting working time is t start The required working time of a day is t span Hour, latest working time t ddl Then X sh1 =(x sh1,1 ,x sh1,2 ,...,x sh1,n ) Can be calculated by the following formula:
Figure BDA0002956527940000051
the constraints that need to be satisfied are as follows:
t start +t span ≤t ddl
convertible load X in said step 1 ch =(x ch,1 ,x ch,2 ,...,x ch,n ) The device is a device with a charging and discharging function in a residential building, such as an electric automobile, a storage battery and the like. X ch Can be seen as a superposition of multiple household electrical appliance transferable loads:
X ch =X ch1 +X ch2 +…
X ch1 =(x ch1,1 ,x ch1,2 ,...,x ch1,n )
X ch2 =(x ch2,1 ,x ch2,2 ,...,x ch2,n )
……
x ch,1 =x ch1,1 +x ch2,1 +...
x ch,2 =x ch1,2 +x ch2,2 +…
with one of the charging and discharging devices X ch1 For example, the full charge capacity is Q ch1 Input power p during charging ch1,ch The charging efficiency is eta ch1,ch With effective power output p at discharge ch1,ex Discharge efficiency of η ch1,ex Increasing penalty efficiency η due to battery discharge causing a decay in battery life ch1,pu At the beginning of workIs m between t start Initial electric quantity is soc start The lowest battery capacity set to soc after the end of the day end Then X ch1 =(x ch1,1 ,x ch1,2 ,...,x ch1,n ) Can be calculated by the following formula:
s bat,i ∈{-1,0,1}
Figure BDA0002956527940000061
Figure BDA0002956527940000062
wherein s is bat,i The charge and discharge states at time i are represented by-1, 0, and 1, respectively, discharge, power-off, and charge. x' ch1,i Is a strategy for actual charging and discharging of the battery. s bat,i The following three constraints should be satisfied:
Figure BDA0002956527940000063
Figure BDA0002956527940000064
Figure BDA0002956527940000065
the first constraint condition indicates that the charging and discharging equipment is required to ensure that the battery capacity is greater than or equal to Soc at each moment min The lowest battery capacity is considered to be defined.
The second constraint condition indicates that the capacity of the charging and discharging device must not be less than Soc after the end of a day end
The third constraint denotes S bat =(s bat,1 ,s bat,2 ,…,s bat,n ) In any continuous gamma charge-discharge stateIn the above, only the charge/discharge state adjustment is allowed at most 1 time, and this is provided to prevent the charge/discharge device from being switched over too frequently.
The reference electrical load Y in step 2 is (Y) 1 ,y 2 ,…,y n ) This can be chosen according to the following different needs response control objectives.
1) If the maximum economic benefit (according to the time-of-use electricity price), the maximum reduction amount at the peak moment (which can be considered to be subsidized according to the maximum reduction amount at the peak moment) or the maximum smooth reduction all day needs to be obtained, the x axis can be used as the reference electricity load, namely, Y is equal to (Y) 1 ,y 2 ,…,y n )=(0,0,…,0)。
2) If the total daily reduction is required to be maximum (namely, the maximum energy saving amount control), the reference power consumption load can be selected as a baseline load, namely, the predicted total daily power consumption load when no demand response is carried out
Figure BDA0002956527940000071
Figure BDA0002956527940000072
Distance function dist (x) of said step 3 i -y i ) The selection can be made according to different requirements and response control purposes:
1) if the maximum economic benefit (according to the time-of-use electricity price) needs to be obtained, the maximum reduction at the peak moment (according to the maximum reduction subsidy at the peak moment) or the maximum reduction all day (namely the maximum energy-saving control) needs to be obtained, the distance function can be expressed as an Euclidean distance:
Figure BDA0002956527940000073
2) if the smoothest reduction is to be achieved throughout the day, the distance function can be expressed as the squared Euclidean distance:
dist(x i -y i )=(x i -y i ) 2
the weight w in step 3 i The selection can be made according to different requirements and response control purposes:
1) if the maximum economic benefit (charged by time of use price) is to be obtained, the weight can be expressed as:
w i =q i
wherein q is i The electricity price corresponding to the time i.
2) The peak time maximum reduction (or subsidized by peak time maximum reduction) weight may be expressed as:
Figure BDA0002956527940000081
where t dr,start ,t dr,start +1,...,t dr,end The peak time (demand response) interval; t is t dr,start Is the demand response start time; t is t dr,end Is the end time of the demand response.
1) The total daily reduction is maximum (i.e. maximum energy saving control), and the distance function can be expressed as:
w i =1
2) if the smoothest reduction is to be achieved throughout the day, the distance function can be expressed as:
w i =1
drawings
FIG. 1 is a flow chart of a demand response control method for multi-energy complementation of residential buildings.
FIG. 2 illustrates load shedding in response to a residential building case demand.
Detailed Description
FIG. 1 is a flow chart of a demand response control method for multi-energy complementation of residential buildings according to the present invention.
The invention is further described in detail below with reference to the drawings and specific embodiments. The method specifically comprises the following steps:
step 1: and according to the demand response control purpose, sequentially determining the reference electric load, the distance function and the weight of the target function. In this embodiment, the control purpose selected by the residential building is the maximum reduction amount of the demand response time, demand responseThe response time is 19:00-21:00, and the reference electric load of the objective function is selected from Y ═ Y 1 ,y 2 ,…,y n ) (0, 0, …, 0), the distance function is chosen to be the euclidean distance
Figure BDA0002956527940000082
Figure BDA0002956527940000083
The weight can be expressed as:
Figure BDA0002956527940000084
step 2: calculating X as X according to the information of the residential building ba +X we +X sh +X ch The values of each item are shown in Table 1. Base load X ba Mainly comprises a refrigerator and a building lighting load, wherein the refrigerator is started all day long, the power is 100W, the lighting power is 660W, and the indoor use time of personnel is 7:00-22: 00. Weather sensitive load X we The hourly load needs to be calculated through EnergyPlus, and the ground source heat pump type air conditioner is adopted for the building. In the simulation, the upper limit t of the temperature set value set,max And a lower limit t set,min Set at 24 ℃ and 28 ℃ respectively. And the parameter alpha for limiting frequent temperature adjustment and the parameter beta for limiting frequent start and stop are 3. Transferable type load X sh Comprising a washing machine, a dish washer and a dryer, respectively using X sh1 ,X sh2 ,X sh3 Is shown, its rated power p sh1 ,p sh2 ,p sh3 1.34kW, 1.25kW and 3.65kW respectively, and the working time t span Are respectively 2h, 2h and 1h, t start 19:00, 20:00 and 21:00 respectively, and the latest working time t ddl 4:00, 5:00 and 6:00, respectively. Convertible load X ch The equipment is an electric automobile with a capacity Q ch1 30kWh, charging input power of 7kW, and charging efficiency eta ch1,ch 85%, the effective power output by battery discharge is 5kW, the efficiency is 80%, and the penalty efficiency eta is ch1,pu 70% of the initial charge, soc start Is 10 percent of the total weight of the composition,minimum battery capacity soc min 20%, the lowest capacity set for the battery capacity after the end of the day is soc end Is 100%.
And 3, step 3: the approximate solution is solved for the objective function through a genetic algorithm, the load reduction situation after the demand response is obtained is shown in figure 2, and the load is transferred and adjusted and controlled in a ratio of 19:00-21: the 12.87kW of the base line load in the 00 period is reduced to 2.60kW, and the reduction amplitude reaches 79.83%.
TABLE 1 residential architecture various kinds of electricity consumption information
Figure BDA0002956527940000091
Figure BDA0002956527940000101

Claims (1)

1. A multi-energy complementary demand response control method for residential buildings is characterized in that the method can meet the maximum economic benefit (charging according to time-of-use electricity price), the maximum reduction amount at peak time (subsidy according to the maximum reduction amount at peak time), the most stable reduction all day, or the maximum reduction amount all day; the actual net electrical load in a day is X ═ X 1 ,x 2 ,…,x n ),All the time of day The optimized objective function is:
Figure FDA0003656414430000011
wherein x is i Is the actual net electrical load in a day, X ═ X 1 ,x 2 ,…,x n ) The value corresponding to the middle time i; y is i For reference, the electrical load Y ═ Y 1 ,y 2 ,…,y n ) The value corresponding to the middle time i; dist (x) i -y i ) Denotes x i And y i A distance function between; w is a i The weight corresponding to the moment i; characterised in that, in particular, the determination of the objective function is carried out byThe values of the terms are solved:
step 1: dividing the actual net electric load X into basic loads X according to the usage and the use property ba Weather sensitive load X we Transferable type load X sh And a convertible load X ch The method comprises the following steps of totally four types of loads, and determining values and constraint conditions through a calculation formula of each type of load:
X=X ba +X we +X sh +X ch
step 2: according to different demand response control purposes, selecting corresponding one-day reference electric load Y (Y ═ Y) 1 ,y 2 ,…,y n );
And step 3: selecting corresponding distance function dist (x) according to different demand response control purposes i -y i ) And a weight w i
And 4, step 4: searching through a group optimization algorithm to obtain optimized equipment-related parameters (start-stop state and temperature setting) which are output as control signals;
the base load X in the step 1 ba =(x ba,1 ,x ba,2 ,...,x ba,n ) In order to maintain the electricity load necessary for the use of the residential building, the elasticity in the time dimension is not provided, such as room LED illumination, electricity consumption of a desktop socket and the like, and the value in one day can be obtained according to the daily use habit of a user;
the weather-sensitive load X in the step 1 we =(x we,1 ,x we,2 ,...,x we,n ) The electric load is mainly used for heating or refrigerating the residential building, and the electric load is influenced by outdoor weather and indoor set temperature;
X we =f 1 (T dry ,T dew ,J,P,T set ,S hvac )
wherein f is 1 Is X we Load prediction model of (1), T dry =(t dry,1, t dry,2 ,…,t dry,n ) The temperature of the dry bulb outside the whole day room; t is dew =(t dew,1 ,t dew,2 ,…,t dew,n ) Is made ofThe wet bulb temperature outside the sky; j ═ J (J) 1 ,j 2 ,…,j n ) Direct solar radiation outside the whole sky room; p ═ P (P) 1 ,p 2 ,…,p n ) Is the atmospheric pressure outside the whole day room; t is set =(t set ,1,t set,2 ,…,t set,n ) The room temperature set value throughout the day; s set =(s set,1, s set,2 ,…,s set,n ) The air conditioning system is in a starting state all day long;
the load prediction model f 1 Modeling calculation can be performed on air conditioning systems such as residential buildings through physical modeling methods such as EnergyPlus, TRNSYS, Modelica and the like, or prediction calculation can be performed on historical data through data driving models such as Artificial Neural Networks (ANN), Support Vector Machines (SVM), extreme learning machines (XGboost) and the like;
the temperature T of the whole day outdoor dry bulb dry All day outdoor wet bulb temperature T dew The all-weather outdoor direct solar radiation J and the all-weather outdoor atmospheric pressure P are meteorological parameter adjustment required by predicting the load of the air conditioning system through physical model calculation or a data driving model and can be obtained through weather forecast data of a local meteorological station;
the set value T of the all-day room temperature set The following four constraints should be satisfied:
t set,min ≤t set,i ≤t set,max
Figure FDA0003656414430000021
Figure FDA0003656414430000031
Figure FDA0003656414430000032
wherein the first constraint representsThe indoor temperature set value should be in the interval of the minimum set value and the maximum set value, and the set range can be slightly larger than the human body comfortable temperature interval provided by ASHRAE due to the requirement of implementing demand response; the second constraint indicates that the indoor temperature setpoint should be the minimum temperature step t step Integral multiple of the air conditioning system, the setting of which should match the control accuracy of the air conditioning system; the third and fourth constraints represent T set =(t set,1 ,t set,2 ,…,t set,n ) In any continuous alpha indoor temperature settings, only at most 1 time of temperature set point adjustment is allowed, and the setting is to prevent the set temperature from being adjusted too frequently;
the all-day air conditioning system is in an open state S set The following two constraints should be satisfied:
s hvac,i ∈{0,1}
Figure FDA0003656414430000033
wherein, the first constraint condition indicates that the on state of the air conditioning system can be off (0) or (1) off; the second constraint denotes S hvac =(s hvac,1 ,s hvac,2 ,…,s hvac,n ) In the opening states of any continuous beta air-conditioning systems, only the opening state adjustment of the air-conditioning system is allowed for at most 1 time, and the setting is to prevent the air-conditioning system from being started and stopped too frequently;
transferable load X in said step 1 sh =(x sh,1 ,x sh,2 ,...,x sh,n ) The electric loads with elasticity in the time dimension during daily use of residential buildings, such as dryers, dishwashers and the like; x sh Can be seen as a superposition of multiple household appliance transferable loads:
X sh =X sh1 +X sh2 +…
X sh1 =(x sh1,1 ,x sh1,2 ,...,x sh1,n )
X sh2 =(x sh2,1 ,x sh2,2 ,...,x sh2,n )
……
x sh,1 =x sh1,1 +x sh2,1 +...
x sh,2 =x sh1,2 +x sh2,2 +…
……
with one of the household appliances X sh1 For example, it operates at a constant power of p sh1 The starting working time is t start The required working time of a day is t span Hour, latest working time t ddl Then X sh1 =(x sh1,1 ,x sh1,2 ,...,x sh1,n ) Can be calculated by the following formula:
Figure FDA0003656414430000041
the constraints that need to be satisfied are as follows:
t start +t span ≤t ddl
convertible load X in said step 1 ch =(x ch,1 ,x ch,2 ,...,x ch,n ) The device is a device with the functions of charging and discharging in a residential building, such as an electric automobile, a storage battery and the like; x ch Can be seen as a superposition of multiple household appliance transferable loads:
X ch =X ch1 +X ch2 +…
X ch1 =(x ch1,1 ,x ch1,2 ,...,x ch1,n )
X ch2 =(x ch2,1 ,x ch2,2 ,...,x ch2,n )
……
x ch,1 =x ch1,1 +x ch2,1 +...
x ch,2 =x ch1,2 +x ch2,2 +…
……
with one of the charging and discharging devices X ch1 For example, the full charge capacity is Q ch1 Input power of p during charging ch1,ch Air chargerElectrical efficiency of η ch1,ch With effective power output p at discharge ch1,ex Discharge efficiency of η ch1,ex Increasing penalty efficiency η due to battery discharge causing a decay in battery life ch1,pu The starting working time is t start Initial electric quantity is soc start The lowest battery capacity set to soc after the end of the day end Then X ch1 =(x ch1,1 ,x ch1,2 ,...,x ch1,n ) Can be calculated by the following formula:
s bat,i ∈{-1,0,1}
Figure FDA0003656414430000051
Figure FDA0003656414430000052
wherein s is bat,i The charge-discharge state at the moment i is represented by-1, 0 and 1, which respectively represent discharge, power failure and charge; x' ch1,i The method is a strategy for actual charging and discharging of the battery; s is bat,i The following three constraints should be satisfied:
Figure FDA0003656414430000053
Figure FDA0003656414430000054
Figure FDA0003656414430000055
the first constraint condition indicates that the charging and discharging equipment is required to ensure that the battery capacity is greater than or equal to soc at each moment min To be considered to define the lowest battery capacity;
the second constraint condition indicates that the capacity of the charging and discharging device must not be less than soc after the end of a day end
The third constraint denotes S bat =(S bat,1 ,S bat,2 ,…,s bat,n ) In any gamma charge-discharge states, only up to 1 charge-discharge state adjustment is allowed, and the charge-discharge state adjustment is set to prevent the charge-discharge equipment from being switched over frequently.
CN202110226687.0A 2021-03-01 2021-03-01 Multifunctional complementary demand response control method for residential building Active CN112968445B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110226687.0A CN112968445B (en) 2021-03-01 2021-03-01 Multifunctional complementary demand response control method for residential building

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110226687.0A CN112968445B (en) 2021-03-01 2021-03-01 Multifunctional complementary demand response control method for residential building

Publications (2)

Publication Number Publication Date
CN112968445A CN112968445A (en) 2021-06-15
CN112968445B true CN112968445B (en) 2022-09-13

Family

ID=76276107

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110226687.0A Active CN112968445B (en) 2021-03-01 2021-03-01 Multifunctional complementary demand response control method for residential building

Country Status (1)

Country Link
CN (1) CN112968445B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106779365A (en) * 2016-12-02 2017-05-31 合肥工业大学 A kind of family's Optimization of Load Dispatching method based on demand response
CN108155639A (en) * 2018-03-13 2018-06-12 华北电力大学 The discrimination method of load is adjusted in a kind of wide area

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111884216B (en) * 2020-07-30 2022-01-14 上海理工大学 Multi-target control method based on building power demand response

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106779365A (en) * 2016-12-02 2017-05-31 合肥工业大学 A kind of family's Optimization of Load Dispatching method based on demand response
CN108155639A (en) * 2018-03-13 2018-06-12 华北电力大学 The discrimination method of load is adjusted in a kind of wide area

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"自适应多场景的单用户柔性负控效果评价算法";郑文明 等;《电力需求侧管理》;20220531;第21卷(第3期);第32-36页 *
中国公共建筑的建筑典型模型建立;陈智博等;《建筑节能》;20200225(第02期);全文 *
单体建筑需求响应研究;罗为等;《建筑节能》;20191125(第11期);第31-35,64页 *
面向电网调峰辅助服务的规模化空调负荷优化调度研究;李作锋等;《电器与能效管理技术》;20180130(第02期);全文 *

Also Published As

Publication number Publication date
CN112968445A (en) 2021-06-15

Similar Documents

Publication Publication Date Title
CN110866641B (en) Two-stage optimization scheduling method and system for multi-energy complementary system considering source storage load coordination
CN106527135B (en) Optimization control method and device for building energy management system
CN104850013A (en) Intelligent electricity utilization method of household appliances
CN102252454B (en) Combined type renewable building energy source supplying system
CN109340904B (en) Electric heating collaborative optimization operation method
CN109787262A (en) A kind of resident's building system active response micro-capacitance sensor Optimization Scheduling
GR1010085B (en) Method for improving the energy management of a nearly zero energy building
CN115857348A (en) Distributed energy system capacity optimization method considering comfortable energy supply of two-combined heat pump
CN103762611A (en) Method for smoothing microgrid call wire power
CN115241927B (en) Operation control method, device, equipment and medium for household intelligent energy system
CN110752629A (en) Energy optimization management method for AC/DC hybrid household micro-grid
CN112128945A (en) Method for providing active power compensation based on battery equivalent model
CN102103362B (en) Micro-grid system and micro-grid system-based fuzzy control method
CN103512151A (en) Method and device for having control over air conditioner installed in area
CN111915107A (en) Load clustering control method based on dynamic clustering
CN112883629A (en) Building power utilization load optimal scheduling method based on improved whale algorithm
CN112968445B (en) Multifunctional complementary demand response control method for residential building
CN113131519A (en) Family energy management optimization method based on mixed integer linear programming
CN110661253B (en) Multi-type electric power elastic load combined adjusting method in building group
CN112880133A (en) Flexible energy utilization control method for building air conditioning system
CN107563547B (en) Comprehensive energy management and control method for optimizing depth of energy consumption of user side
CN110729728A (en) Demand response method considering real-time comfort of user and load rate of power grid
CN113420413B (en) Flexible load adjustability quantification method and system based on load plasticity
CN105549389A (en) Household energy management algorithm based on building thermodynamic model
CN202101470U (en) Composite type renewable building energy supplying system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant