# A Solution to Prevent a Blackout Crisis: Determining the Behavioral Potential and Capacity of Solar Power

Authors: Saeed Vedadi Kalantar; Amir Ali Saifoddin; Ahmad Hajinezhad; Mohammad Hossein Ahmadi
Journal: International Journal of Photoenergy
 (2021)

---

## Abstract

Increasing the power grid peak in the summer time causes power outages in industries and residential areas in Iran. The most obvious example of this issue is the power outages in the summer of 2018. Management of the demand-side is the most important strategy to reduce the grid peak due to the high cost of the development of the power plant capacity (500$ per kilowatt). In the present study, the effect of behavioral parameters in decreasing the power grid peak was identified. The behavioral simulation was done as an agent-based model using the raw data of the time-use survey (TUS) of the Statistics Center of Iran. 4228 urban households were surveyed, and the quality of people’s behavior was determined in each time step of 15 minutes during the day and night with 2 deterministic and stochastic approaches. In the stochastic approach, the Markov chain method was used. It showed that the power grid peak can only be reduced by 10% with behavioral flexibility and up to 25% by upgrading technology. In addition, based on the power deficit in 2018, 2000 megawatts of solar power capacity must be added to the network at peak times to meet grid demand.


---

## Body

## 1. Introduction

Air conditioning) AC) is capable of saving lives during warm weather and has a major contribution to peak electricity demand during summer time. Throughout exceedingly hot days or heat waves, such an elevation in energy demand results in brownouts and dimouts [1]. The electric energy and electricity industry are also key factors for the economic system’s propulsion in order not to affect the economic system and public security setting [2]. In many countries with significantly demanding seasonal cooling, the role of AC in peak electricity needs is remarkably greater than overall year-round demands. For instance, AC demand in Madrid, Spain, encompassed one-third of overall peak consumption in June 2008 [3]. In such countries as Singapore or most of those in the Middle East, which have completely year-round cooling demands, the AC contribution to peak loads can be as large as 50% or higher [3]. In the plethora of countries, a sharp elevation is predicted in contribution of space cooling to peak electricity loading, and the elevations are uppermost in hot countries, including India, in which the contribution rises from only 10% currently to 45% in 2050 [4]. As many cooling demands are fulfilled with electricity-powered fans or ACs, power systems are greatly influenced by the growing need for cooling. Elevated AC loadings particularly increase not only total electricity requirement but also peak electricity loadings. Application of ACs and electric fans for staying cool comprises approximately 20% of overall electricity, which is currently utilized in buildings worldwide [4, 5]. Such a tendency starts to develop as the focus of global economic and demographic development is toward hotter countries. Oftentimes, demand-side management (DSM) programs are primarily aimed at lowering space-cooling loadings via demand response at periods of peak electricity demands [6]. The overall energy requirement of a building results from both its thermic setting and its inhabitants, which have a greater complication in the assessment and quantification than the building envelope and its thermic setting [7]. Therefore, the management of energy demands in buildings requires critical consideration of the behavior of occupants [8, 9]. In this research, the amount of cooling load in the power grid peak and its reduction potential by observing behavioral issues in the form of codified scenarios are investigated by studying the behavioral details of household energy consumption using the summertime-use survey of the Statistics Center of Iran in 2015 and by agent-based modeling of the residential sector in the Tehran. The bottlenecks and important goals pursued in this study are as follows: (1) development of bottom-up modeling tools in terms of behavioral functions; (2) accuracy of estimating the effective cooling load in the peak of Tehran electricity network by different building groups; (3) investigating the effect of increasing technology efficiency, comfort temperature, and a combination of both in the peak of the electricity network; (4) determining the behavioral potentials of reducing cooling energy consumption and its impact on the peak of the electricity network; and (5) determining the share of solar electricity to manage the power grid peak. Behavioral modeling methods in the field of energy with 2 stochastic and deterministic approaches and agent-based and bottom-up modeling methods are reviewed in the literature. The following organization is presented in this study. In Section 2, the background of blackouts in Iran is reported. In Section 3, the study methodological procedure is explained. In Section 4, the results and discussion is represented. And finally, in Section 5, the main research conclusions and discussions are synthesized.


## 2. Background

### 2.1. Blackouts and Brownouts in Iran

The power supply system includes the generation, transmission, and distribution sections, and if any of these sections are disrupted, blackouts will occur [10]. The trend of demand and production capacity at the peak moment from 2007 to 2020 was increasing (from about 35000 MW in 2007 to about 58000 MW in 2020) (Figure 1). During these years, except for 2019 and 2020, the amount of power demand exceeded the production amount at the peak moment.Figure 1
Maximum production capacity and consuming demand in network peak time in Iran. Extracted from the Iran Grid Management Company (IGMC) database.Table1 shows the hours of the power grid peak in the last 10 years. It is noteworthy that in 2020, the grid peak occurred at night, which indicates that the behavior of the energy consumption is changed and also the peak hours are changed from day to night.Table 1
Power grid peaks (hour and amount) in the last 10 years. Extracted from the Iran Grid Management Company (IGMC) and Tavanir Company databases.
YearPower grid peakPower grid peakDay and night peak difference (GW)Day (GW)HourNight (GW)HourJuly 21, 202058.24412:4058.25421:250.01July 22, 201957.63516:2356.56321:241.072July 11, 201857.09816:3854.94622:322.152July 30, 201755.44214:4552.45221:362.99July 20, 201652.69314:5050.321:302.393July 11, 201550.171548222.17July 26, 201448.52715:2446.3222.227July 17, 201345.69314:4344.8220.893August 14, 2012431642221August 3, 2011421541221In 2018, the capacity to generate electricity during peak hours was about 8000 MW less than the demand due to reduced rainfall and a sharp decline in the capacity of hydropower plants. If hydropower plants are taken out of operation, the lack of capacity of the country’s thermal power plants is about 8000 MW. The plan of the development of the capacity of thermal power plants in 2019 and 2020 was about 7000 MW (Table2) of which only 1500 MW of new capacity was built.Table 2
Ministry of Energy plans to increase capacity of power plants (MW). Extracted from the Thermal Power Plants Holding Company database.
Share of each section20172018201920202021Governmental23542652938760Private9585432099998555Private buyback agreement1604804801785640Cumulative sum13531065323137211955The share of different power plants in electricity generation and different sectors in electricity consumption is as shown in Figures2 and 3. To increase energy security, Iran must diversify its energy supply [11] and power plants. The highest electricity generation in Iran is in combined cycle power plants (36%) and the highest electricity consumption is in the industrial (36%) and residential (32%) sectors.Figure 2
The share of various power plants in electricity generation in Iran. Extracted from the Iran Grid Management Company (IGMC) database.Figure 3
The share of different sectors of electricity consumption in Iran. Extracted from the Iran Grid Management Company (IGMC) database.According to the power plant construction process in the country from 1997 to 2020 (Figure4), the capacity of power plants can be increased by an average of 2000 MW per year, which cannot compensate for the lack of capacity of thermal power plants, especially in years of water shortage. Therefore, to solve the blackout crisis, the solution must be sought along with increasing the capacity of power plants on the demand side [12].Figure 4
Construction of power plant trend (1995 to 2019). Extracted from Thermal Power Plants Holding Company database.As shown in Figure5, the industry’s electricity consumption in the peak hours of summer 2019 and 2020 is about 5000 MW, which is more than the storage power plants at this time (about 3000 MW). The high consumption of industrial electricity at the power grid peak shows the lack of a specific program in the field of the management of the load in this sector.Figure 5
Comparison of industrial electricity consumption and power plant reserves at the peak moment in 2019 and 2020. Extracted from the Iran Grid Management Company (IGMC) database.The peak of the power grid in most countries of the world, including Iran, is affected by the cooling load. As shown in Figure6, the share of the cooling load from the grid peak in the United States of America is nearly 30% and in China and South Korea is more than 15% [4].Figure 6
The share of cooling load from the peak of the electricity network in different countries.The temperature in the hot summer months in Tehran averages about 36 degrees Celsius, which indicates the high potential of cooling electricity consumption in this city (Figures7 and 8).Figure 7
Temperature in Tehran from June 22 to July 22.Figure 8
Temperature in Tehran from July 23 to August 22.


### 2.2. Importance of Solar Power in Meeting the Energy Demands

The demand for solar electricity in the world is increasing, so that by 2050, the cheapest type of electricity is projected to be solar electricity. The installed capacity of solar power plants in the world by 2019 is equal to 627 thousand MW. More than 29 countries had more than a thousand megawatts of solar electricity by 2019 [13]. The price of solar panels is falling with a significant slope [14]. Solar electricity has the largest share of renewable energy types connected to the electricity grid. Some countries, such as Algeria, have paid special attention to solar power to meet the demand of the people [15]. The growth rate of the capacity of solar power plants in the world compared to last year was 28%. Iran has a high potential in the field of solar energy with 1,648,000 km2 and 300 sunny days and 2200 kW/m2 of radiation capacity. The distribution of solar radiation capacity in Iran is as shown in Figure 9.Figure 9
The distribution of solar radiation capacity in Iran.


## 2.1. Blackouts and Brownouts in Iran

The power supply system includes the generation, transmission, and distribution sections, and if any of these sections are disrupted, blackouts will occur [10]. The trend of demand and production capacity at the peak moment from 2007 to 2020 was increasing (from about 35000 MW in 2007 to about 58000 MW in 2020) (Figure 1). During these years, except for 2019 and 2020, the amount of power demand exceeded the production amount at the peak moment.Figure 1
Maximum production capacity and consuming demand in network peak time in Iran. Extracted from the Iran Grid Management Company (IGMC) database.Table1 shows the hours of the power grid peak in the last 10 years. It is noteworthy that in 2020, the grid peak occurred at night, which indicates that the behavior of the energy consumption is changed and also the peak hours are changed from day to night.Table 1
Power grid peaks (hour and amount) in the last 10 years. Extracted from the Iran Grid Management Company (IGMC) and Tavanir Company databases.
YearPower grid peakPower grid peakDay and night peak difference (GW)Day (GW)HourNight (GW)HourJuly 21, 202058.24412:4058.25421:250.01July 22, 201957.63516:2356.56321:241.072July 11, 201857.09816:3854.94622:322.152July 30, 201755.44214:4552.45221:362.99July 20, 201652.69314:5050.321:302.393July 11, 201550.171548222.17July 26, 201448.52715:2446.3222.227July 17, 201345.69314:4344.8220.893August 14, 2012431642221August 3, 2011421541221In 2018, the capacity to generate electricity during peak hours was about 8000 MW less than the demand due to reduced rainfall and a sharp decline in the capacity of hydropower plants. If hydropower plants are taken out of operation, the lack of capacity of the country’s thermal power plants is about 8000 MW. The plan of the development of the capacity of thermal power plants in 2019 and 2020 was about 7000 MW (Table2) of which only 1500 MW of new capacity was built.Table 2
Ministry of Energy plans to increase capacity of power plants (MW). Extracted from the Thermal Power Plants Holding Company database.
Share of each section20172018201920202021Governmental23542652938760Private9585432099998555Private buyback agreement1604804801785640Cumulative sum13531065323137211955The share of different power plants in electricity generation and different sectors in electricity consumption is as shown in Figures2 and 3. To increase energy security, Iran must diversify its energy supply [11] and power plants. The highest electricity generation in Iran is in combined cycle power plants (36%) and the highest electricity consumption is in the industrial (36%) and residential (32%) sectors.Figure 2
The share of various power plants in electricity generation in Iran. Extracted from the Iran Grid Management Company (IGMC) database.Figure 3
The share of different sectors of electricity consumption in Iran. Extracted from the Iran Grid Management Company (IGMC) database.According to the power plant construction process in the country from 1997 to 2020 (Figure4), the capacity of power plants can be increased by an average of 2000 MW per year, which cannot compensate for the lack of capacity of thermal power plants, especially in years of water shortage. Therefore, to solve the blackout crisis, the solution must be sought along with increasing the capacity of power plants on the demand side [12].Figure 4
Construction of power plant trend (1995 to 2019). Extracted from Thermal Power Plants Holding Company database.As shown in Figure5, the industry’s electricity consumption in the peak hours of summer 2019 and 2020 is about 5000 MW, which is more than the storage power plants at this time (about 3000 MW). The high consumption of industrial electricity at the power grid peak shows the lack of a specific program in the field of the management of the load in this sector.Figure 5
Comparison of industrial electricity consumption and power plant reserves at the peak moment in 2019 and 2020. Extracted from the Iran Grid Management Company (IGMC) database.The peak of the power grid in most countries of the world, including Iran, is affected by the cooling load. As shown in Figure6, the share of the cooling load from the grid peak in the United States of America is nearly 30% and in China and South Korea is more than 15% [4].Figure 6
The share of cooling load from the peak of the electricity network in different countries.The temperature in the hot summer months in Tehran averages about 36 degrees Celsius, which indicates the high potential of cooling electricity consumption in this city (Figures7 and 8).Figure 7
Temperature in Tehran from June 22 to July 22.Figure 8
Temperature in Tehran from July 23 to August 22.


## 2.2. Importance of Solar Power in Meeting the Energy Demands

The demand for solar electricity in the world is increasing, so that by 2050, the cheapest type of electricity is projected to be solar electricity. The installed capacity of solar power plants in the world by 2019 is equal to 627 thousand MW. More than 29 countries had more than a thousand megawatts of solar electricity by 2019 [13]. The price of solar panels is falling with a significant slope [14]. Solar electricity has the largest share of renewable energy types connected to the electricity grid. Some countries, such as Algeria, have paid special attention to solar power to meet the demand of the people [15]. The growth rate of the capacity of solar power plants in the world compared to last year was 28%. Iran has a high potential in the field of solar energy with 1,648,000 km2 and 300 sunny days and 2200 kW/m2 of radiation capacity. The distribution of solar radiation capacity in Iran is as shown in Figure 9.Figure 9
The distribution of solar radiation capacity in Iran.


## 3. Method: Stochastic and Agent-Based Modeling

The process of stochastic and agent-based modeling of this study is shown in Figure10. In this research, the variable state is simulated based on the stochastic model and the basic state is deterministic.Figure 10
The process of stochastic and agent base modeling (behavioral and bottom-up).


### 3.1. Agent-Based Models

In agent-based simulations that consider behavioral issues, two deterministic and stochastic approaches are used. Deterministic agent-based models follow certain schedules. Schedules include load profiles that are numerically represented between zero and 1 [16]. Stochastic agent-based models use statistical distributions for determining residents’ conditions and behavior [17]. Most stochastic models are developed based on the first-order Markov chain techniques [18, 19]. Data from stochastic models are extracted from time-use surveys that describe consumer behavior [20]. The number of stochastic agent-based studies is four times more than deterministic approaches [21]. In this research, a stochastic agent-based model is used.


### 3.2. Behavioral Parameters

Density, presence of residents, and their desired comfort temperature are considered as behavioral parameters of the model. An agent-based cooling simulation that takes into account the presence of individuals differs by approximately 10% from the actual state [22]. This issue shows the high performance of behavioral models to correctly detect the presence of people in the building. It should be noted that most of the behavioral actions of cooling energy consumption in the building occur in the form of comfort temperature (Figure 11). For instance, Hoyt et al. [23] report that the comfort temperature of cooling and heating reflects the behavioral parameters of the household. The behavioral parameters of the research are presented in Tables 3 and 4.Figure 11
Influence of different behavioral characteristics of cooling energy consumption from the presence of people and comfort temperature.Table 3
Behavioral parameters of the model.
Specifications of the modelParameterDescriptionProbability of the presence of people at homeStochastic and deterministicThe cooling load is a function of the presence of people and the comfort temperature. The probability of the presence of people at home in the form of a Markov matrix and with the help of the 2015 TUS has been extracted.Number of occupants by building groups in Table11 person“Results of Statistics on Energy Consumption and Environmental Characteristics of Urban Households in 2016” by the Statistics Center of Iran2 persons3 persons4 persons5 persons6 persons and moreCooling comfort temperatures by building groups in Table118°C“Results of Statistics on Energy Consumption and Environmental Characteristics of Urban Households in 2016” by the Statistics Center of Iran22°C27°C30°CPopulation density per m2 of building based on building groups0.06550 m2 and lessDevelopment by the author0.0551-75 m20.04176-80 m20.03881-100 m20.028101-150 m20.02151-200 m20.0138201-300 m20.0115301 m2 and moreTable 4
Number and density of people in building groups.
Number of residentsBuilding groups1 person2 persons3 persons4 persons5 persons6 personsTotalTotal area (m2)Population density per m250 m2 and less1336150303115762179617784033302347047171594860/06551-75 m26334635087085219110394074411921719292918937584733690/0576-80 m21543310821819260628885812486574445804429194672760/04181-100 m2273972334745025027596943616972036762088443547060260/038101-150 m23076124961665208510187105216482851322757956974326810/028151-200 m26830558261137921838359107365444516802258572020/02201-300 m2470031812675391086984313238353294236213284040/0138301 m2 and more1804155112340119116189972526310409590451110/0115


### 3.3. Determining the Presence of People in the Building Using Time-Use Survey: A Stochastic Approach

A behavioral database of people with the appropriate resolution, which is created by the government in different countries in the form of a time-use survey, is required for stochastic behavioral modeling [24]. The model of Walker and Pokoski [25] is one of the first time-use-based stochastic models. The Richardson model, developed in the United Kingdom based on the 2000 time-use survey, was the inspiration for stochastic behavioral research in this field [18]. The census of the time-use survey conducted by the Statistics Center of Iran in 2014 and 2015 is according to Table 5. Since our study is related to the power grid peak and this peak occurs in the summer, the census data of summer 2015 was used.Table 5
Characteristics of the time-use survey data in this research.
RowTitleDescription1Target societyPeople over 15 years in urban families2Statistical timeThe statistical time of this plan, depending on the case, is last week or one continuous day and night during the census week.3Data collection methodIn this survey, the required household information is collected through face-to-face interviews with the most informed household member and the completion of a questionnaire. Personal information is also collected in a self-reported manner by completing a questionnaire by each person aged 15 years and older in the sample household.4Sample sizeThis study was conducted in urban areas of the country for 4 seasons. The sample size in the whole statistical year (autumn and winter 2014 and spring and summer 2015) was equal to 16912 households. In summer, 8248 people belonging to 4228 households were referred.5Time stepTime-use of people every 15 minutes is questioned.6ActivitiesAll activities that a person does during the day are divided into 15 main groups according to the ICATUS. These 15 groups eventually merge into 9 groups.7Place of activityThe location of the activity is given in 22 codes. Three of these 22 codes indicate the presence of people in the house, and 19 codes indicate the absence of people in the house.15 activity codes and their classification into 9 groups were performed according to Table6.Table 6
Code of activities and its classification into different groups.
CodeClassification ICATUSMerged classification01Working for companies or quasicompanies, nonprofit organizations, and government (work in the formal sector)Work and job activities (codes 01 to 05)02Working for the family in basic productive activities03Working for the household in nonprimary productive activities04Working for the family in construction activities05Working for the family to provide services for income06Providing unpaid domestic services for final consumption in your householdHousehold activities (codes 06 to 07)07Providing unpaid care services to family members08Providing services to the local community and helping other householdsVoluntary and charitable activities (code 08)09LearningLearning activities (code 09)10Social participation and sociabilitySocial participation (code 10)11Going to places or visiting cultural, recreational, and sports eventsEntertainment (codes 11 and 12)12Entertainment, games, and other recreational activities13Participating in indoor and outdoor sports and related coursesSports activities (code 13)14Using mass mediaUse of mass media (code 14)15Personal maintenance and carePersonal maintenance and care (code 15)The process of extracting the presence/absence statistics of individuals from the raw data of the time-use survey is shown in Figure12.Figure 12
A process model for extracting and monitoring data from a time-use survey.


### 3.4. Markov Chain

In the present study, the Markov chain method was used to understand the stochastic presence of people in the house. Most studies with a stochastic approach use the Markov chain method [26]. A Markov chain is a model for displaying a sequence of random variables in which the probability of each event occurring depends only on the previous event. So, the probability of occurrence of events in such a model depends only on the previous time, and other events do not interfere with the probability.
(1)PXt+1=x∣X1=x1,X2=x2,⋯,Xn=xt=PXt+1=xXt=xt.The Markov chain matrix was displayed as a transition probability matrix.Pij indicates the probability of transition from point i to j. In this way, the transition matrix represents the values of the state space set. In the present study, the state space has 2 elements (Equation (2)). The transition probability matrix is also in accordance with Equation (3). The Markov chain matrix can be represented as a transition probability matrix. Pij indicates the probability of transition from point i to j.
(2)S=1,2=Absenceinthebuilding,presenceinthebuilding,(3)Pij=P11P12P21P22.In this matrix, according to the principles of probability, it should be∑J=1KPij=1, which is the sum of the probabilities of each row in the matrix must be equal to 1. In general, in the Markov chain, pijn indicates the probability of reaching from state i to j in step.
(4)Pijn=∑rϵsPirKPrjn−k·∀K·0<K<n.The transition time step in this study is 15 minutes, and the transition probability matrix for each time step should be determined using the Markov chain. An example of a Markov chain step with and without considering the household dimension (Table4) is shown in Tables 7 and 8.Table 7
A step from the Markov chain regardless of the family dimension.
ResolutionModeNumberDeterministic probabilityNext modeNumberProbabilityStochastic probability00:15Presence819199.3089Presence818699.93890.9924Absence50.0610Absence570.6910Presence000.0075Absence57100Table 8
A step from the Markov chain considering the family dimension.
ResolutionProbability of presence at least one person at home (deterministic)Probability of presence at least one person at home (stochastic)00:150.9999670.999357


### 3.5. Bottom-Up Model

Table9 shows the specifications of the bottom-up model used to simulate cooling peak loads. According to the data of the Tehran Municipality website, there are more than 2,867,000 residential houses in eight-building groups in the capital of Iran, which in this respect has the largest share in the network peak [27]. In this simulation, the computational core of the EnergyPlus program is used through the user interface of the DesignBuilder software. The open-source and modular features of E+ facilitate the addition of new simulation modules for developers [28]. The computing core has been developed as open-source codes by the US Department of Energy [29]. The Ashrae 90.1 standard is used as the basic mode for simulating new versions of the EnergyPlus Computational Core [16]. The simulated cooling load was converted into electrical load by considering the common cooling equipment in the country according to data of the Statistics Center of Iran (Table 10). A 3-story apartment is simulated for all building groups (Figure 13). The bottom-up model is extended to all building groups with the archetype method.Table 9
Specifications of the bottom-up model.
Specifications of modelParameterDescriptionSimulation siteTehranAccording to the meteorological data of Tehran in the database of DesignBuilder softwareBuilding type3-story apartmentThe first and third floors correspond to the first and last floors of each building and the second floor as the middle floors of the buildingsComputational coreEnergyPlusUse DesignBuilder Graphics Interface 6.1.4.00.7Building direction20-degree angle with north directionThe predominant orientation of Tehran buildingsBuilding cooling coefficient96.1%96.1% of buildings in Tehran are cooled.Simulation time steps15 minutesBasic simulation modeAshrae 90.1According to the basic mode of DesignBuilder softwareBuilding materialsBrick, iron, and other common materials in the interior of the buildingAccording to the 2016 Population and Housing Census by the Statistical Centre of IranCooling technologyCommon technologies“Results of Statistics on Energy Consumption and Environmental Characteristics of Urban Households in 2016” by the Statistics Center of IranExtension methodArchetype methodThis method divides the building section into groups based on infrastructural, dimensional similarities, materials used, and more. The estimated energy for the representative of each group is developed for the whole building sector.Building groups50 m2 and lessSample building: 50 m251-75 m2Sample building: 65 m276-80 m2Sample building: 80 m281-100 m2Sample building: 90 m2101-150 m2Sample building: 125 m2151-200 m2Sample building: 175 m2201-300 m2Sample building: 250 m2301 m2 and moreSample building: 300 m2Table 10
Distribution of cooling equipment used in residential buildings—“Results of Statistics on Energy Consumption and Environmental Characteristics of Urban Households in 2016” by the Statistics Center of Iran.
Cooling equipment%Energy grade%Absorption chiller0/4A38.9B8.2Not stated52.9Compression chiller0.1Fan coil0.2A19.2C16.5E26.7None13.2Not stated24.3Refrigeration package0.2A21.5B35.4C5Not stated38Split cooler9.8A32.7B13.7C2.8D1.3E1F0.1G0.2None15.1Not stated33.1Water cooler56.8A9.4B5.3C4.7D1.9E2F1.6G1.3None31.8Not stated42Fan5.4Water cooler and split cooler0.8Water cooler and split cooler and fan0.7Water cooler and fan14Split cooler and fan11.5Other0.1100Figure 13
A sample simulated house for a building group.


### 3.6. Scenarios

In this section, seven scenarios about the comfort temperature and technology upgrade have been developed, and the effects of each on the electric cooling load of the power grid peak are investigated. The formulation of the scenarios is based on the main research model, namely, stochastic and agent-based model.


#### 3.6.1. Ideal Scenario

The ideal comfort temperature in Tehran is about 27 degrees Celsius (Heidari, 2009). 23% (percentage of dispersion of gas coolers, absorption, and compression chillers; water coolers cannot adjust the temperature) of buildings with comfort temperatures of 18, 22, and 30 degrees Celsius set their comfort temperature at 27 degrees. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.


#### 3.6.2. Cold Room Scenario

23% (percentage of dispersion of gas coolers, absorption, and compression chillers; water coolers cannot adjust the temperature) of buildings with comfort temperatures of 18, 27, and 30 degrees Celsius set their comfort temperature at 22 degrees. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.


#### 3.6.3. Refrigerator Scenario

23% (percentage of dispersion of gas coolers, absorption, and compression chillers; water coolers cannot adjust the temperature) of buildings with a comfort temperature of 22, 27, and 30 degrees Celsius set their comfort temperature at 18 degrees. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.


#### 3.6.4. Technology Upgrade Scenario

Considering that 65% of the country’s buildings use water coolers [30] and 22% use split coolers [31], one of the ways to control the electricity peak is increasing the efficiency of these two types of widely used cooling equipment in buildings. In this scenario, the dominant grade of cooling equipment in Iran, namely, grade G, has been upgraded to grade A. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.


#### 3.6.5. Combined Scenarios

A combination of the above three scenarios (idea, cold room, and refrigerator) and the technology upgrade scenario is examined. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.Figure14 shows the flowchart of the research simulation with details.Figure 14
Research modeling flowchart.


### 3.7. Limitations of Model

The limitations of modeling in this study are as follows:(i)
Markov chain development of the presence of people at home with a time step of 15 minutes has been done for the first time in Iran. Refining the raw data of the Statistics Center’s TUS and turning it into a Markov probability chain database is a difficult process with modeling complexities, which in turn has increased the implicit innovations in this research(ii)
The use of behavioral modeling methods with emphasis on human identity to study the peak of the electricity network has been done for the first time in this research. Lack of data and its integration is one of the limitations of this research(iii)
Some of the new concepts introduced in this study, such as behavioral flexibility, did not exist before. Therefore, its formulation was limited


## 3.1. Agent-Based Models

In agent-based simulations that consider behavioral issues, two deterministic and stochastic approaches are used. Deterministic agent-based models follow certain schedules. Schedules include load profiles that are numerically represented between zero and 1 [16]. Stochastic agent-based models use statistical distributions for determining residents’ conditions and behavior [17]. Most stochastic models are developed based on the first-order Markov chain techniques [18, 19]. Data from stochastic models are extracted from time-use surveys that describe consumer behavior [20]. The number of stochastic agent-based studies is four times more than deterministic approaches [21]. In this research, a stochastic agent-based model is used.


## 3.2. Behavioral Parameters

Density, presence of residents, and their desired comfort temperature are considered as behavioral parameters of the model. An agent-based cooling simulation that takes into account the presence of individuals differs by approximately 10% from the actual state [22]. This issue shows the high performance of behavioral models to correctly detect the presence of people in the building. It should be noted that most of the behavioral actions of cooling energy consumption in the building occur in the form of comfort temperature (Figure 11). For instance, Hoyt et al. [23] report that the comfort temperature of cooling and heating reflects the behavioral parameters of the household. The behavioral parameters of the research are presented in Tables 3 and 4.Figure 11
Influence of different behavioral characteristics of cooling energy consumption from the presence of people and comfort temperature.Table 3
Behavioral parameters of the model.
Specifications of the modelParameterDescriptionProbability of the presence of people at homeStochastic and deterministicThe cooling load is a function of the presence of people and the comfort temperature. The probability of the presence of people at home in the form of a Markov matrix and with the help of the 2015 TUS has been extracted.Number of occupants by building groups in Table11 person“Results of Statistics on Energy Consumption and Environmental Characteristics of Urban Households in 2016” by the Statistics Center of Iran2 persons3 persons4 persons5 persons6 persons and moreCooling comfort temperatures by building groups in Table118°C“Results of Statistics on Energy Consumption and Environmental Characteristics of Urban Households in 2016” by the Statistics Center of Iran22°C27°C30°CPopulation density per m2 of building based on building groups0.06550 m2 and lessDevelopment by the author0.0551-75 m20.04176-80 m20.03881-100 m20.028101-150 m20.02151-200 m20.0138201-300 m20.0115301 m2 and moreTable 4
Number and density of people in building groups.
Number of residentsBuilding groups1 person2 persons3 persons4 persons5 persons6 personsTotalTotal area (m2)Population density per m250 m2 and less1336150303115762179617784033302347047171594860/06551-75 m26334635087085219110394074411921719292918937584733690/0576-80 m21543310821819260628885812486574445804429194672760/04181-100 m2273972334745025027596943616972036762088443547060260/038101-150 m23076124961665208510187105216482851322757956974326810/028151-200 m26830558261137921838359107365444516802258572020/02201-300 m2470031812675391086984313238353294236213284040/0138301 m2 and more1804155112340119116189972526310409590451110/0115


## 3.3. Determining the Presence of People in the Building Using Time-Use Survey: A Stochastic Approach

A behavioral database of people with the appropriate resolution, which is created by the government in different countries in the form of a time-use survey, is required for stochastic behavioral modeling [24]. The model of Walker and Pokoski [25] is one of the first time-use-based stochastic models. The Richardson model, developed in the United Kingdom based on the 2000 time-use survey, was the inspiration for stochastic behavioral research in this field [18]. The census of the time-use survey conducted by the Statistics Center of Iran in 2014 and 2015 is according to Table 5. Since our study is related to the power grid peak and this peak occurs in the summer, the census data of summer 2015 was used.Table 5
Characteristics of the time-use survey data in this research.
RowTitleDescription1Target societyPeople over 15 years in urban families2Statistical timeThe statistical time of this plan, depending on the case, is last week or one continuous day and night during the census week.3Data collection methodIn this survey, the required household information is collected through face-to-face interviews with the most informed household member and the completion of a questionnaire. Personal information is also collected in a self-reported manner by completing a questionnaire by each person aged 15 years and older in the sample household.4Sample sizeThis study was conducted in urban areas of the country for 4 seasons. The sample size in the whole statistical year (autumn and winter 2014 and spring and summer 2015) was equal to 16912 households. In summer, 8248 people belonging to 4228 households were referred.5Time stepTime-use of people every 15 minutes is questioned.6ActivitiesAll activities that a person does during the day are divided into 15 main groups according to the ICATUS. These 15 groups eventually merge into 9 groups.7Place of activityThe location of the activity is given in 22 codes. Three of these 22 codes indicate the presence of people in the house, and 19 codes indicate the absence of people in the house.15 activity codes and their classification into 9 groups were performed according to Table6.Table 6
Code of activities and its classification into different groups.
CodeClassification ICATUSMerged classification01Working for companies or quasicompanies, nonprofit organizations, and government (work in the formal sector)Work and job activities (codes 01 to 05)02Working for the family in basic productive activities03Working for the household in nonprimary productive activities04Working for the family in construction activities05Working for the family to provide services for income06Providing unpaid domestic services for final consumption in your householdHousehold activities (codes 06 to 07)07Providing unpaid care services to family members08Providing services to the local community and helping other householdsVoluntary and charitable activities (code 08)09LearningLearning activities (code 09)10Social participation and sociabilitySocial participation (code 10)11Going to places or visiting cultural, recreational, and sports eventsEntertainment (codes 11 and 12)12Entertainment, games, and other recreational activities13Participating in indoor and outdoor sports and related coursesSports activities (code 13)14Using mass mediaUse of mass media (code 14)15Personal maintenance and carePersonal maintenance and care (code 15)The process of extracting the presence/absence statistics of individuals from the raw data of the time-use survey is shown in Figure12.Figure 12
A process model for extracting and monitoring data from a time-use survey.


## 3.4. Markov Chain

In the present study, the Markov chain method was used to understand the stochastic presence of people in the house. Most studies with a stochastic approach use the Markov chain method [26]. A Markov chain is a model for displaying a sequence of random variables in which the probability of each event occurring depends only on the previous event. So, the probability of occurrence of events in such a model depends only on the previous time, and other events do not interfere with the probability.
(1)PXt+1=x∣X1=x1,X2=x2,⋯,Xn=xt=PXt+1=xXt=xt.The Markov chain matrix was displayed as a transition probability matrix.Pij indicates the probability of transition from point i to j. In this way, the transition matrix represents the values of the state space set. In the present study, the state space has 2 elements (Equation (2)). The transition probability matrix is also in accordance with Equation (3). The Markov chain matrix can be represented as a transition probability matrix. Pij indicates the probability of transition from point i to j.
(2)S=1,2=Absenceinthebuilding,presenceinthebuilding,(3)Pij=P11P12P21P22.In this matrix, according to the principles of probability, it should be∑J=1KPij=1, which is the sum of the probabilities of each row in the matrix must be equal to 1. In general, in the Markov chain, pijn indicates the probability of reaching from state i to j in step.
(4)Pijn=∑rϵsPirKPrjn−k·∀K·0<K<n.The transition time step in this study is 15 minutes, and the transition probability matrix for each time step should be determined using the Markov chain. An example of a Markov chain step with and without considering the household dimension (Table4) is shown in Tables 7 and 8.Table 7
A step from the Markov chain regardless of the family dimension.
ResolutionModeNumberDeterministic probabilityNext modeNumberProbabilityStochastic probability00:15Presence819199.3089Presence818699.93890.9924Absence50.0610Absence570.6910Presence000.0075Absence57100Table 8
A step from the Markov chain considering the family dimension.
ResolutionProbability of presence at least one person at home (deterministic)Probability of presence at least one person at home (stochastic)00:150.9999670.999357


## 3.5. Bottom-Up Model

Table9 shows the specifications of the bottom-up model used to simulate cooling peak loads. According to the data of the Tehran Municipality website, there are more than 2,867,000 residential houses in eight-building groups in the capital of Iran, which in this respect has the largest share in the network peak [27]. In this simulation, the computational core of the EnergyPlus program is used through the user interface of the DesignBuilder software. The open-source and modular features of E+ facilitate the addition of new simulation modules for developers [28]. The computing core has been developed as open-source codes by the US Department of Energy [29]. The Ashrae 90.1 standard is used as the basic mode for simulating new versions of the EnergyPlus Computational Core [16]. The simulated cooling load was converted into electrical load by considering the common cooling equipment in the country according to data of the Statistics Center of Iran (Table 10). A 3-story apartment is simulated for all building groups (Figure 13). The bottom-up model is extended to all building groups with the archetype method.Table 9
Specifications of the bottom-up model.
Specifications of modelParameterDescriptionSimulation siteTehranAccording to the meteorological data of Tehran in the database of DesignBuilder softwareBuilding type3-story apartmentThe first and third floors correspond to the first and last floors of each building and the second floor as the middle floors of the buildingsComputational coreEnergyPlusUse DesignBuilder Graphics Interface 6.1.4.00.7Building direction20-degree angle with north directionThe predominant orientation of Tehran buildingsBuilding cooling coefficient96.1%96.1% of buildings in Tehran are cooled.Simulation time steps15 minutesBasic simulation modeAshrae 90.1According to the basic mode of DesignBuilder softwareBuilding materialsBrick, iron, and other common materials in the interior of the buildingAccording to the 2016 Population and Housing Census by the Statistical Centre of IranCooling technologyCommon technologies“Results of Statistics on Energy Consumption and Environmental Characteristics of Urban Households in 2016” by the Statistics Center of IranExtension methodArchetype methodThis method divides the building section into groups based on infrastructural, dimensional similarities, materials used, and more. The estimated energy for the representative of each group is developed for the whole building sector.Building groups50 m2 and lessSample building: 50 m251-75 m2Sample building: 65 m276-80 m2Sample building: 80 m281-100 m2Sample building: 90 m2101-150 m2Sample building: 125 m2151-200 m2Sample building: 175 m2201-300 m2Sample building: 250 m2301 m2 and moreSample building: 300 m2Table 10
Distribution of cooling equipment used in residential buildings—“Results of Statistics on Energy Consumption and Environmental Characteristics of Urban Households in 2016” by the Statistics Center of Iran.
Cooling equipment%Energy grade%Absorption chiller0/4A38.9B8.2Not stated52.9Compression chiller0.1Fan coil0.2A19.2C16.5E26.7None13.2Not stated24.3Refrigeration package0.2A21.5B35.4C5Not stated38Split cooler9.8A32.7B13.7C2.8D1.3E1F0.1G0.2None15.1Not stated33.1Water cooler56.8A9.4B5.3C4.7D1.9E2F1.6G1.3None31.8Not stated42Fan5.4Water cooler and split cooler0.8Water cooler and split cooler and fan0.7Water cooler and fan14Split cooler and fan11.5Other0.1100Figure 13
A sample simulated house for a building group.


## 3.6. Scenarios

In this section, seven scenarios about the comfort temperature and technology upgrade have been developed, and the effects of each on the electric cooling load of the power grid peak are investigated. The formulation of the scenarios is based on the main research model, namely, stochastic and agent-based model.


### 3.6.1. Ideal Scenario

The ideal comfort temperature in Tehran is about 27 degrees Celsius (Heidari, 2009). 23% (percentage of dispersion of gas coolers, absorption, and compression chillers; water coolers cannot adjust the temperature) of buildings with comfort temperatures of 18, 22, and 30 degrees Celsius set their comfort temperature at 27 degrees. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.


### 3.6.2. Cold Room Scenario

23% (percentage of dispersion of gas coolers, absorption, and compression chillers; water coolers cannot adjust the temperature) of buildings with comfort temperatures of 18, 27, and 30 degrees Celsius set their comfort temperature at 22 degrees. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.


### 3.6.3. Refrigerator Scenario

23% (percentage of dispersion of gas coolers, absorption, and compression chillers; water coolers cannot adjust the temperature) of buildings with a comfort temperature of 22, 27, and 30 degrees Celsius set their comfort temperature at 18 degrees. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.


### 3.6.4. Technology Upgrade Scenario

Considering that 65% of the country’s buildings use water coolers [30] and 22% use split coolers [31], one of the ways to control the electricity peak is increasing the efficiency of these two types of widely used cooling equipment in buildings. In this scenario, the dominant grade of cooling equipment in Iran, namely, grade G, has been upgraded to grade A. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.


### 3.6.5. Combined Scenarios

A combination of the above three scenarios (idea, cold room, and refrigerator) and the technology upgrade scenario is examined. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.Figure14 shows the flowchart of the research simulation with details.Figure 14
Research modeling flowchart.


## 3.6.1. Ideal Scenario

The ideal comfort temperature in Tehran is about 27 degrees Celsius (Heidari, 2009). 23% (percentage of dispersion of gas coolers, absorption, and compression chillers; water coolers cannot adjust the temperature) of buildings with comfort temperatures of 18, 22, and 30 degrees Celsius set their comfort temperature at 27 degrees. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.


## 3.6.2. Cold Room Scenario

23% (percentage of dispersion of gas coolers, absorption, and compression chillers; water coolers cannot adjust the temperature) of buildings with comfort temperatures of 18, 27, and 30 degrees Celsius set their comfort temperature at 22 degrees. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.


## 3.6.3. Refrigerator Scenario

23% (percentage of dispersion of gas coolers, absorption, and compression chillers; water coolers cannot adjust the temperature) of buildings with a comfort temperature of 22, 27, and 30 degrees Celsius set their comfort temperature at 18 degrees. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.


## 3.6.4. Technology Upgrade Scenario

Considering that 65% of the country’s buildings use water coolers [30] and 22% use split coolers [31], one of the ways to control the electricity peak is increasing the efficiency of these two types of widely used cooling equipment in buildings. In this scenario, the dominant grade of cooling equipment in Iran, namely, grade G, has been upgraded to grade A. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.


## 3.6.5. Combined Scenarios

A combination of the above three scenarios (idea, cold room, and refrigerator) and the technology upgrade scenario is examined. Behavioral characteristics of people during the day and night are determined based on the results of TUS. Grid stability is considered the same during the day and night.Figure14 shows the flowchart of the research simulation with details.Figure 14
Research modeling flowchart.


## 3.7. Limitations of Model

The limitations of modeling in this study are as follows:(i)
Markov chain development of the presence of people at home with a time step of 15 minutes has been done for the first time in Iran. Refining the raw data of the Statistics Center’s TUS and turning it into a Markov probability chain database is a difficult process with modeling complexities, which in turn has increased the implicit innovations in this research(ii)
The use of behavioral modeling methods with emphasis on human identity to study the peak of the electricity network has been done for the first time in this research. Lack of data and its integration is one of the limitations of this research(iii)
Some of the new concepts introduced in this study, such as behavioral flexibility, did not exist before. Therefore, its formulation was limited


## 4. Results

### 4.1. Model Validation

According to Figure15, the results of the basic state of the model showed that the electric load of the model on the peak day is 2 to 8% different from the real state, which is a good adaptation.Figure 15
Model validation, comparison of simulated electricity load, and real electricity load in Tehran.


### 4.2. Results of Basic State

The distribution of high electricity consumption hours during the day and night in the summer of 2019 and 2020 is shown in Figure16. The distribution of peak hours of the night in summer 2019 and 2020 is similar and between 20 and 22, especially at 21:30. The peak hours of the day in the summer of 2019 were between 14:00 and 16:50. However, most of the peak hours in summer 2020 were around 12:00 and 14:30. The peak of the day in summer 2020 compared to 2019 has 2 hours back.Figure 16
Distribution of day and night peak times in summer 2019 and 2020.Figure17 shows the stochastic and deterministic probabilities of the presence of people in the building. The maximum difference between the probabilities of at least one person at home in a stochastic state with a deterministic state is 5%, which occurs at 14:24. From 12 a.m. to 7:12 p.m., the difference between stochastic and deterministic probabilities is between 2 and 3%. Between 10 and 15 o’clock, the difference between the stochastic and deterministic probabilities is about 2%. From 16:48 to 24, the difference between the stochastic and deterministic probabilities of the presence of people decreases and goes to zero.Figure 17
Stochastic and deterministic probabilities of presence.The most excellent behavioral flexibility in reducing the peak is 12 to 14 (Figure18). Therefore, a strategy in peak shaving is to direct the peak time to these hours.Figure 18
Basic electricity consumption and comparison of stochastic and deterministic electric load (day).From 11 a.m. to 5 p.m., the behavioral flexibility potential of household cooling energy consumption is between 70 and 134 MW, which is approximately equal to 10% of the cooling peak of Tehran’s electricity grid. Table11 shows the potential for behavioral flexibility on the network peak day by building groups.Table 11
Behavioral flexibility potential of cooling energy consumption by different building groups (day).
Building groupsAverage potential for behavioral flexibility (MW)TimespanMaximum potential for behavioral flexibility (MW)Time point50 m2 and lessAbout 4.312 to 16:155.813:3051-75 m2About 25.411 to 16:30391276-80 m2About 5.310:30 to 178.751281-100 m2About 1410:30 to 16:452712101-150 m2About 1810:45 to 16:153412151-200 m2About 6.810:45 to 16:4512.512201-300 m2About 1.510:30 to 183.112301 m2 and moreAbout 0.810:45 to 16:451.412From 1:45 to 6 a.m., the behavioral flexibility potential of household cooling energy consumption is about 50 to 65 MW (Figure19 and Table 12), which is approximately equal to 5% of the cooling peak of the Tehran power grid.Figure 19
Basic electricity consumption and comparison of stochastic and deterministic electric load (night).Table 12
Behavioral flexibility potential of cooling energy consumption by different building groups (night).
Building groupsAverage potential for behavioral flexibility (MW)TimespanMaximum potential for behavioral flexibility (MW)Time point50 m2 and lessAbout 1.21:45 to 62.94:1551-75 m2About 15.22 to 617.24:3076-80 m2About 0.652 to 62.84:3081-100 m2About 51:45 to 611.24:15101-150 m2About 71:45 to 616.34:15151-200 m2About 1.12 to 62.84:15201-300 m2About 0.91:45 to 61.94:15301 m2 and moreAbout 0.22:45 to 60.64:15


### 4.3. Results of the Scenario

The cooling load peaks in the refrigerator (1.48 GW) and combined (ideal temperature and technology upgrade) (0.97 GW) scenarios are the maximum and minimum cooling loads in Tehran, respectively. The distance of 500 MW of these two scenarios indicates the high potential of reducing the peak of Tehran’s power grid by adjusting the comfort temperature and technology upgrade. As the efficiency of cooling technologies increases, the average cooling load in all comfort temperature scenarios decreases between 250 and 400 MW, and up to 25% of the network peak is reduced compared to the base scenario. The maximum difference between stochastic and deterministic electric loads in the refrigerator and cold scenarios is 154 MW and 139 MW, respectively. This shows the flexibility of cooling energy consumption based on the presence of people in the building, in different scenarios. As productivity decreases, the potential for behavioral flexibility in energy consumption increases. The same thing happened in night peak scenarios (Figures20–33 and Tables 13 and 14).Figure 20
Different modes in the ideal comfort temperature scenario (day).Figure 21
Different modes in the cold room scenario (day).Figure 22
Different modes in the refrigerator scenario (day).Figure 23
Different modes in the technology upgrade scenario (day).Figure 24
Different modes in the combined scenario (ideal temperature and technology upgrade) (day).Figure 25
Different modes in the combined scenario (cold room and technology upgrade) (day).Figure 26
Different modes in the combined scenario (refrigerator and technology upgrade) (day).Figure 27
Different modes in the ideal comfort temperature scenario (night).Figure 28
Different modes in the cold room scenario (night).Figure 29
Different modes in the refrigerator scenario (night).Figure 30
Different modes in the technology upgrade scenario (night).Figure 31
Different modes in the combined scenario (ideal temperature and technology upgrade) (night).Figure 32
Different modes in the combined scenario (cold room and technology upgrade) (night).Figure 33
Different modes in the combined scenario (refrigerator and technology upgrade) (night).Table 13
Comparison of different scenarios (day).
ScenariosStochasticDeterministicMaximum of behavioral flexibility (MW)Behavioral flexibility compared to base scenario (MW)Maximum load (GW)Maximum difference with the base scenario (MW)Decrease or increase compared to the base scenario in maximum load (%)Maximum load (MW)Maximum difference with the base scenario (MW)Decrease or increase compared to the base scenario in maximum load (%)Maximum difference with the base scenario (MW)Base scenario1.29——1.33————134.64Ideal scenario1.24-48.7-3.751.28-47.4-3.8-48.7-4.94129.7Cold room scenario1.3342.273.001.3741.123.142.274.28138.92Refrigerator scenario1.48195.914.281.52190.514.72195.919.83154.47Technology upgrade scenario (basic mode and technology upgrade)1.008-293.66-22.031.037-285.62-21.8-293.66-29.71104.93Combined scenario (ideal temperature and technology upgrade)0.97-332.4-24.960.998-323.3-24.8-332.4-33.63101.01Combined scenario (cold room and technology upgrade)1.03-261.53-19.621.069-254.36-20.15-261.53-26.46108.18Combined scenario (refrigerator and technology upgrade)1.15-141.89-11.271.18-138-10.85-141.89-14.34120.3Table 14
Comparison of different scenarios (night).
ScenariosDeterministicStochasticMaximum of behavioral flexibility (MW)Behavioral flexibility compared to base scenario (MW)Maximum difference with the base scenario (MW)Decrease or increase compared to the base scenario in maximum load (%)Maximum load (MW)Maximum difference with the base scenario (MW)Decrease or increase compared to the base scenario in maximum load (%)Maximum load (GW)—57.86——0.934——0.932Base scenario-3.2654.6-53.44-5.780.880-53.44-5.680.879Ideal scenario1.5859.4425.852.780.96025.852.780.958Cold room scenario7.6265.48122.9513.161.057122.937.831.055Refrigerator scenario-12.7745.09-206.15-22.050.728-206.12-22.100.726Technology upgrade scenario (basic mode and technology upgrade)-15.3342.53-248.33-26.650.685-248.30-26.600.684Combined scenario (ideal temperature and technology upgrade)-11.5846.28-186.57-20.020.747-186.55-20.060.745Combined scenario (cold room and technology upgrade)-6.8750.99-110.96-11.880.823-110.94-11.800.822Combined scenario (refrigerator and technology upgrade)


## 4.1. Model Validation

According to Figure15, the results of the basic state of the model showed that the electric load of the model on the peak day is 2 to 8% different from the real state, which is a good adaptation.Figure 15
Model validation, comparison of simulated electricity load, and real electricity load in Tehran.


## 4.2. Results of Basic State

The distribution of high electricity consumption hours during the day and night in the summer of 2019 and 2020 is shown in Figure16. The distribution of peak hours of the night in summer 2019 and 2020 is similar and between 20 and 22, especially at 21:30. The peak hours of the day in the summer of 2019 were between 14:00 and 16:50. However, most of the peak hours in summer 2020 were around 12:00 and 14:30. The peak of the day in summer 2020 compared to 2019 has 2 hours back.Figure 16
Distribution of day and night peak times in summer 2019 and 2020.Figure17 shows the stochastic and deterministic probabilities of the presence of people in the building. The maximum difference between the probabilities of at least one person at home in a stochastic state with a deterministic state is 5%, which occurs at 14:24. From 12 a.m. to 7:12 p.m., the difference between stochastic and deterministic probabilities is between 2 and 3%. Between 10 and 15 o’clock, the difference between the stochastic and deterministic probabilities is about 2%. From 16:48 to 24, the difference between the stochastic and deterministic probabilities of the presence of people decreases and goes to zero.Figure 17
Stochastic and deterministic probabilities of presence.The most excellent behavioral flexibility in reducing the peak is 12 to 14 (Figure18). Therefore, a strategy in peak shaving is to direct the peak time to these hours.Figure 18
Basic electricity consumption and comparison of stochastic and deterministic electric load (day).From 11 a.m. to 5 p.m., the behavioral flexibility potential of household cooling energy consumption is between 70 and 134 MW, which is approximately equal to 10% of the cooling peak of Tehran’s electricity grid. Table11 shows the potential for behavioral flexibility on the network peak day by building groups.Table 11
Behavioral flexibility potential of cooling energy consumption by different building groups (day).
Building groupsAverage potential for behavioral flexibility (MW)TimespanMaximum potential for behavioral flexibility (MW)Time point50 m2 and lessAbout 4.312 to 16:155.813:3051-75 m2About 25.411 to 16:30391276-80 m2About 5.310:30 to 178.751281-100 m2About 1410:30 to 16:452712101-150 m2About 1810:45 to 16:153412151-200 m2About 6.810:45 to 16:4512.512201-300 m2About 1.510:30 to 183.112301 m2 and moreAbout 0.810:45 to 16:451.412From 1:45 to 6 a.m., the behavioral flexibility potential of household cooling energy consumption is about 50 to 65 MW (Figure19 and Table 12), which is approximately equal to 5% of the cooling peak of the Tehran power grid.Figure 19
Basic electricity consumption and comparison of stochastic and deterministic electric load (night).Table 12
Behavioral flexibility potential of cooling energy consumption by different building groups (night).
Building groupsAverage potential for behavioral flexibility (MW)TimespanMaximum potential for behavioral flexibility (MW)Time point50 m2 and lessAbout 1.21:45 to 62.94:1551-75 m2About 15.22 to 617.24:3076-80 m2About 0.652 to 62.84:3081-100 m2About 51:45 to 611.24:15101-150 m2About 71:45 to 616.34:15151-200 m2About 1.12 to 62.84:15201-300 m2About 0.91:45 to 61.94:15301 m2 and moreAbout 0.22:45 to 60.64:15


## 4.3. Results of the Scenario

The cooling load peaks in the refrigerator (1.48 GW) and combined (ideal temperature and technology upgrade) (0.97 GW) scenarios are the maximum and minimum cooling loads in Tehran, respectively. The distance of 500 MW of these two scenarios indicates the high potential of reducing the peak of Tehran’s power grid by adjusting the comfort temperature and technology upgrade. As the efficiency of cooling technologies increases, the average cooling load in all comfort temperature scenarios decreases between 250 and 400 MW, and up to 25% of the network peak is reduced compared to the base scenario. The maximum difference between stochastic and deterministic electric loads in the refrigerator and cold scenarios is 154 MW and 139 MW, respectively. This shows the flexibility of cooling energy consumption based on the presence of people in the building, in different scenarios. As productivity decreases, the potential for behavioral flexibility in energy consumption increases. The same thing happened in night peak scenarios (Figures20–33 and Tables 13 and 14).Figure 20
Different modes in the ideal comfort temperature scenario (day).Figure 21
Different modes in the cold room scenario (day).Figure 22
Different modes in the refrigerator scenario (day).Figure 23
Different modes in the technology upgrade scenario (day).Figure 24
Different modes in the combined scenario (ideal temperature and technology upgrade) (day).Figure 25
Different modes in the combined scenario (cold room and technology upgrade) (day).Figure 26
Different modes in the combined scenario (refrigerator and technology upgrade) (day).Figure 27
Different modes in the ideal comfort temperature scenario (night).Figure 28
Different modes in the cold room scenario (night).Figure 29
Different modes in the refrigerator scenario (night).Figure 30
Different modes in the technology upgrade scenario (night).Figure 31
Different modes in the combined scenario (ideal temperature and technology upgrade) (night).Figure 32
Different modes in the combined scenario (cold room and technology upgrade) (night).Figure 33
Different modes in the combined scenario (refrigerator and technology upgrade) (night).Table 13
Comparison of different scenarios (day).
ScenariosStochasticDeterministicMaximum of behavioral flexibility (MW)Behavioral flexibility compared to base scenario (MW)Maximum load (GW)Maximum difference with the base scenario (MW)Decrease or increase compared to the base scenario in maximum load (%)Maximum load (MW)Maximum difference with the base scenario (MW)Decrease or increase compared to the base scenario in maximum load (%)Maximum difference with the base scenario (MW)Base scenario1.29——1.33————134.64Ideal scenario1.24-48.7-3.751.28-47.4-3.8-48.7-4.94129.7Cold room scenario1.3342.273.001.3741.123.142.274.28138.92Refrigerator scenario1.48195.914.281.52190.514.72195.919.83154.47Technology upgrade scenario (basic mode and technology upgrade)1.008-293.66-22.031.037-285.62-21.8-293.66-29.71104.93Combined scenario (ideal temperature and technology upgrade)0.97-332.4-24.960.998-323.3-24.8-332.4-33.63101.01Combined scenario (cold room and technology upgrade)1.03-261.53-19.621.069-254.36-20.15-261.53-26.46108.18Combined scenario (refrigerator and technology upgrade)1.15-141.89-11.271.18-138-10.85-141.89-14.34120.3Table 14
Comparison of different scenarios (night).
ScenariosDeterministicStochasticMaximum of behavioral flexibility (MW)Behavioral flexibility compared to base scenario (MW)Maximum difference with the base scenario (MW)Decrease or increase compared to the base scenario in maximum load (%)Maximum load (MW)Maximum difference with the base scenario (MW)Decrease or increase compared to the base scenario in maximum load (%)Maximum load (GW)—57.86——0.934——0.932Base scenario-3.2654.6-53.44-5.780.880-53.44-5.680.879Ideal scenario1.5859.4425.852.780.96025.852.780.958Cold room scenario7.6265.48122.9513.161.057122.937.831.055Refrigerator scenario-12.7745.09-206.15-22.050.728-206.12-22.100.726Technology upgrade scenario (basic mode and technology upgrade)-15.3342.53-248.33-26.650.685-248.30-26.600.684Combined scenario (ideal temperature and technology upgrade)-11.5846.28-186.57-20.020.747-186.55-20.060.745Combined scenario (cold room and technology upgrade)-6.8750.99-110.96-11.880.823-110.94-11.800.822Combined scenario (refrigerator and technology upgrade)


## 5. Conclusions

This study has estimated the electric cooling load of Tehran at the peak moment by stochastic and agent-based modeling. Density, presence of residents, and comfort temperature separately sorted by all building groups are considered as behavioral parameters of the model. Stochastic and deterministic presence of people based on the raw census data of the time-use survey of the Statistics Center of Iran in the summer of 2015, which was performed among 4228 households, was modeled using the Markov chain method for each time step of 15 minutes. The EnergyPlus computing core (DesignBuilder user interface) was used to develop the bottom-up model. In this study, scenarios based on household comfort temperature were developed. The results of this study showed that the most behavioral flexibility in reducing the peak is related to the 12 to 14 o’clock. So, a strategy in peak shaving is to direct the peak time to these hours. From 11 a.m. to 5 p.m., the behavioral flexibility potential of household cooling energy consumption is between 70 and 134 MW, which is approximately equal to 10% of the cooling peak of Tehran’s electricity grid. The cooling load peak in the refrigerator (1.48 GW) and combined (ideal temperature and technology upgrade) (0.97 GW) scenarios are the maximum and minimum cooling loads in Tehran, respectively. The distance of 500 MW of these two scenarios indicates the high potential of reducing the peak of Tehran’s power grid by adjusting the comfort temperature and technology upgrade. The maximum difference between stochastic and deterministic electric loads in the refrigerator and cold scenarios is 154 MW and 139 MW, respectively. As productivity decreases, the potential for behavioral flexibility in energy consumption increases. As the efficiency of cooling technologies increases, the average cooling load in all comfort temperature scenarios decreases between 250 and 400 MW, and up to 25% of the network peak is reduced compared to the base scenario. The same thing happened in night peak scenarios. This study, we developed possible scenarios of household comfort temperature and technology upgrade that provide a powerful tool for policymakers to manage and control the grid peak, while estimating the behavioral flexibility of reducing the power grid peak (almost 10%) and modeling energy consumption based on behavioral parameters with a resolution of 15 minutes. Based on the power deficit in 2018, 2000 megawatts of solar power capacity must be added to the network at peak times to meet grid demand.



---

*Source: 2092842-2021-11-09.xml*