Demand Response (DR) is one of the most cost-effective and unfailing techniques used by utilities for consumer load shifting. This research paper presents different DR programs in deregulated environments. The description and the classification of DR along with their potential benefits and associated cost components are presented. In addition, most DR measurement indices and their evaluation are also highlighted. Initially, the economic load model incorporated thermal, wind, and energy storage by considering the elasticity market price from its calculated locational marginal pricing (LMP). The various DR programs like direct load control, critical peak pricing, real-time pricing, time of use, and capacity market programs are considered during this study. The effect of demand response in electricity prices is highlighted using a simulated study on IEEE 30 bus system. Simulation is done by the Shuffled Frog Leap Algorithm (SFLA). Comprehensive performance comparison on voltage deviations, losses, and cost with and without considering DR is also presented in this paper.

Power utilities in the electrical system are undergoing a major rebuilding process and enter into the deregulated market operations throughout the world. In power systems, the challenges are made to build the proficiency of the mechanical segment and diminish the electrical power expenses of all the clients. To meet the power demand, electric power industries are moving towards renewable power generation sources such as wind [

In the deregulated market system, locational marginal pricing is a centralized process of market clearing, in which the Independent System Operator (ISO) is responsible to ascertain the power dispatch and energy prices. Network limits must be considered during scheduling loads, generators and bilateral transactions unlike in the case of system uniform pricing (i.e., unconstrained bidding) approach. Due to the consideration of network constraints in the market clearing process, it is not probable to ascertain the market equilibrium by simply analyzing the intersection of a cumulative demand curve and cumulative supply curve. Instead, the energy prices and the power dispatch schedules are calculated based on an optimization approach that consists of power flow and network-related constraints.

Electricity demand is primarily modeled as inelastic in power systems optimization problems generally. Though in real-life problems, a substantial amount of electricity demand is considered elastic. Electric loads such as district heating, PEV charging; HVAC systems are considered as flexible demands and comprise a significant percentage of the total demand. Distribution companies have good knowledge of the amount of flexibility in demand based on historic data. The majority of these demands are suspendable in which part of the demand is shifted concerning by considering both rate constraints and deadlines. Most of the distribution system operators employ some class of DR programs. DR is represented in the form of voluntary reduction of demand by customers or interruptible load contracts (ILC) based on electricity prices. These DR programs not only offer ancillary services to operate the power systems but also help in optimizing the generation from renewable sources. Dispatch models and deterministic security-constrained commitment are considered in traditional power system operation. These models utilize conservative forecasts of wind power generation [

To enhance the reliability of an electricity network and efficiency, the resource of demand response (DR) is used to describe the electricity demand. DR is usually typified as an incentive and price-based DR [

In this study, the generation allocation at an economic level (wind, ESS, and thermal) in the deregulated market system is mainly dealt with. By considering the price elasticity of the market based on the model is developed and reported. The developed load model consists of contracted price, penalty costs spot market price (LMP), and incentives. Under the previously developed ED model with wind, the load model is incorporated. Based on the above factors and fuel cost function, DR cost is calculated for every change in demand. The various DR programs such as time of use, direct load control, capacity market programs, real-time pricing, and critical peak pricing are mainly considered. By following these DR programs, peak demand is reduced and the load curve is shaped. Economic allocation for generators is done especially for shaped load curves that decrease the cost of generation and further the lines congestion. The simulation is done by SFLA and it is specifically tested on IEEE 30 bus system.

In a business sector clearing, locational marginal pricing is a concentrated procedure in which the Independent System Operator (ISO) is responsible for deciding the power dispatch plans along with energy costs. During the scheduling of loads, generators, and bilateral exchanges, system limits should be considered unlike the uniform pricing approach in the system. As the system limitations are considered during the business sector clearing process, it is not practical to analyze business sector harmony simply by considering the convergence of the total demand curve and the combined supply curve. Further, energy costs and power dispatch schedules are computed using an optimization approach that comprises power and system flow-related limitations.

Electricity LMP at the particular location (bus) is outlined as the minimum cost that is used to service the next increment of demand at a specific location along with all the constraints present during the operation of the power system. It is assumed that the minimized objective function F

Whereas, F and x indicate the total cost objective function and the solution vector comprises optimum values for the decision variables respectively. LMP at the real power and each bus “k” is generally Lagrange multiplier in addition to the constraints present in real power balance for the concerned or specific bus

DR is divided into two categories such as incentive and time-based programs. Mainly in time-based programs like Real Time Pricing (RTP), different periods have different electricity prices based on the electricity supply cost. For instance, in the case of peak period at a high price, the price is mostly medium for a peak at the off-stage and the price is less for a low load period. Hence, there is no penalty or incentive for these DR programs. Demand bidding (DB), direct load control (DLC), capacity market program (CAP), ancillary service (A/S) programs, and Interruptible/Curtailable (I/C) service are mainly considered as incentive-based programs. Direct load control is usually defined as a voluntary program and if the customers do not limit consumption, the inclusion of penalty charges is mandated. Whereas, CAP and I/C are compulsory programs, in which enrolled customers should pay incentives and they are limited during directed. There are large number of customers who are expectant to lessen the price of the load which they are much interested to curtail during demand bidding programs. The customers are offered reductions in load in electricity markets which are as operation reserves in case of ancillary services programs.

Interruptible/Curtailable service (I/C): The customers on I/C service tariffs have bill credit or discount during the exchange and also it is agreed to lessen the load in case of system emergencies. A penalty is included when the customers do not curtail. In the case of I/C tariffs, the size of the minimum customer must be higher than 200 KW.

The tariffs for customers are agreed either to decrease a definite block of electric load or to curtail consumption to a pre-specified level. These tariffs must be curtailed within 30–60 min after being notified by the advanced metering infrastructure (AMI) utility. The utility calls the interruption for nearly 200 h per year. All customers do not apply for interruptible programs. As in the silicon chip production industry, customers present in the continuous processes are not good candidates. In this case, incentive tariffs may vary in different markets.

The pre-specified load reductions are considered in the capacity market programs during the rise of system contingencies and are subjected to penalties when it is not curtailed if directed. The programs in the capacity market are mainly considered as insurance. The participants receive guaranteed payments when directed. In some years, load curtailments are not called though the participants are paid to be on call. The wholesale market providers like ISOs can take up these programs and operate the installed capacity markets which are the basic structural market analog of I/C tariffs.

Further, the obligation to curtail is to be agreed upon when the eligibility of the capacity market program is merely based on the reductions which are quite feasible and bearable. The provisions to entertain capacity payments are as follows. (i) The minimum reduction of around 100 KW, (ii) The minimum reduction for four hours, (iii) It is notified every two hours, and (iv) It is subjected to one inspection or test per capability period.

OPF problem formulation along with energy and wind storage in the deregulated environment is considered in this study. These studies detail the modeling of the load in the deregulated environment.

It is significant to include the development of a responsive load economic model to express the participation of the customers in DRPs on the characteristics of load profile. The problems present in DRP constraints and objective function are discussed in this section.

Load economic model describes the customer behavior towards electricity price change, incentives or penalties imposed. Elasticity is defined as the ratio of demand sensitivity concerning the price _{o} is the initial demand.

From ^{th} period following the kth period is calculated using

The difference in electricity consumption at the customer end by the change in electricity prices at different time periods are as follows. Some of the loads are quite firm and it is not shifted from one time period to another (for instance, lighting loads). Hence, it remains in two states either as off or on. Further, these types of loads are quite sensitive to a single period and it is described as ‘‘self-elasticity” which is constantly having a negative value. As like process loads, it is flexible and is transferrable from one period to the other, but it is quite sensitive to multi-period and is calculated using “cross-elasticity” and it constantly has a positive value.

In case there is a variation between initial load demand (d0 (t)) to new load demand d(t), the change in load demand is calculated using

If IC(t) is the incentive value remunerated to the customer in t^{th} hour for each MWh reduction in load, the total incentive for their participation in I/C, CAP, and DLC,

According to the contract, if the registered customers in DRPs have not achieved the concerned commitments, the system is penalized. If pen(t) and H(t) are the penalty for the identical period and contract level for the t^{th} hour respectively, the total penalty is accounted using

The consumption by the customers is accounted using

d(t) and d0 (t) are equal if the electricity price value is assumed to be the same for the employment of DRPs before and after consideration of the penalty and incentives.

The consumption by the customers after inclusion of cross elasticity and load change for 24 h period is calculated using

The combination of

The main objective of the DROPF problem is the total system cost minimization for different DRPs implementation to satisfy other operational constraints and load demands. The operational cost is the total fuel cost for every wind storage, thermal unit, and energy storage over a defined period. The implementation cost of DRPs is accounted for during the total system cost for every scheduling hour. Hence, the objective function is expressed as follows.

Subject to constraints

where PGi, t & QGi, t are the sum of the real power and reactive power injections (thermal, wind, and energy storage) at the bus i at time t.

Further state of charge of energy storage for a hour t + 1 is given as,_{gi} is Real power, V_{gi} is Bus voltage, W_{gi} is Wind generation, Q_{gi} is Reactive power, SOC is State of Charge, E_{Storage}-Energy storage, _{dis} is discharge efficiency and _{ch}is charging efficiency.

Hence,

The conventional generator constraints are subjected to consider current operating point deviation. Hence, the amount of change is limited during the generation based on the ramp rate of individual generators. The constraints are expressed during the increment in generation are as follows.

where _{i} _{i}^{th} unit in MW.

This algorithm is a memetic algorithm and is mostly used for the food hunting behavior of frogs. It is based on the evolution of memes conceded by the global exchange of information among themselves and the interactive frogs. This program is a combination of random and deterministic approaches. It combines the benefits of a genetic-based memetic algorithm and a social behavior-based PSO algorithm. This is mainly used to solve many complex optimization models which are non-differentiable and nonlinear.

In the initial step of the algorithm, it is used to generate initial population P of frogs arbitrarily in the search space. Further, the position of ith frog is indicated as Xi = [ Xi, 1, … Xi,D] whereas, D indicates the number of variables. Furthermore, the frogs are sorted in descending order by their fitness. In addition, the population is entirely partitioned into m subsets and is indicated as memeplexes in which each subset contains n frogs (P = m* n). The strategy of partitioning steps is as follows. (i) The first frog leaves into first memeplex, second to the second memeplex, the mth frog into mth memeplex, and (m + 1) frog goes to the first memeplex and hence so forth. In each memeplex, the position of frogs with the best and worst are identified as Xb and Xw respectively. Also, the position of the frog with global best is Xg. Then within each memeplex, a process similar to the PSO algorithm is applied to improve only the frog with the worst fitness in each cycle using the following equation

New position _{w} = _{w} + _{i}

where

Dimax is the maximum allowed variation in the position of frogs. Fitness function to be evaluated is given by,

SFLA algorithm is explained in the following flowchart

This section deals with the LMP and DR programs based simulation results. Simulations are carried out using SFLA in a 2.4 GHz processor using Matlab software.

At first, the effect of wind generators and ESSs on LMP is studied. As a by-product of ACOPF, LMP at each bus can be calculated. The location of the wind generator, as well as demands, play an important role in the calculation of LMP. IEEE 30 bus system is taken as the test system.

Initially, LMPs are calculated for the system with conventional generators. Then the wind generator is placed at different buses and the variations in LMPs are analyzed. Since the deregulated market is considered, LMP is the major factor affecting the market bids and contracts. Hence, more focus should be given on placing the wind generator. The wind generator is placed at different buses from bus 1 to bus 30(one at a time) and the obtained maximum LMP is shown in

Inferences:

Maximum LMP varies from 3.3142 to 56.1425 for different locations of wind generators.

From the LMPs obtained for 30 buses for each location of wind generator, it is noted that the LMP of the corresponding bus where the wind generators are installed is very minimum.

Wind generators installed at buses 2 to 7, resulted in very high LMP around 55 INR/MWh. Hence, these are not suitable locations for installing wind generators.

Except buses 8 and 28, the maximum LMPs obtained are around 14 to 15 INR/MWh.

Hence, from the results, the suitable locations of wind generators are bus no. 8 and 28 which resulted in a minimum LMP of around 3.5.

From the previous results, installing a wind generator either at bus 8 or 28 will result in minimum LMP. Further, the demand variations also affect the LMP. Here LMP calculation is performed for several values of total load demand with an equal interval of 10% load increase between 150 and 500 MW, considering the existence of all resources (thermal, wind, and ESS) to supply this load demand. An increase in load demand causes an increase in the LMP value as shown in

Different LMP curves are generated for different load curves. Each step of load increase and its equivalent LMP increase is because of the use of all available resources to meet the demand. For each of these steps, load curtailment using DRP is studied. Due to involved prices, DR programs are only useful for high demand. For higher demands, DR programs could be effective. LMP and load curve for a day ahead are shown in

Basic OPF problem incorporating different DRPs is implemented on IEEE 30 bus system. The 24-hour load demand is taken from

Program | Initial spot price (ρ) (INR/MWh) | Incentive (INR/MWh) | Penalty (INR/MWh) |
---|---|---|---|

Base | 11.6 flat rate (average value from LMP curve) | 0 | 0 |

DLC | 11.6 flat rate | 5 | 0 |

CAP | 11.6 flat rate | 5 | 4 |

I/C | 11.6 flat rate | 7 | 5 |

TOU | 8, 10, and 13 at the valley, off-peak and peak period | 0 | 0 |

CPP | 14 at 9, 10 and 11 h | 0 | 0 |

RTP | Individual LMP at each bus by solving ACOPF | 0 | 0 |

This case is solved with the actual load curve and without the implementation of any DR programs. The total load of the day is 23940 MW with a peak of 1305 MW. The average wind power is found by wind speed and power forecasting [

In this case, DRP is modeled such that, as the electricity price increases above the average price of the day during peak hours, customers will tend to curtail their load. The spot electricity price is considered 11.6 INR/MWh (from the average LMP curve). The incentive INR 15 is given to the customers in peak hours for their load reduction. Compared to the base case, the system total load is now reduced by 7.3% (22161 MW) with a peak value reduction of 7.5% (1207 MW) and the total cost is INR 606889.

This is an incentive-based mandatory type of DRP, so the penalty is also included in this program. The value of incentive and penalty is considered as

This is an incentive and penalty-based market program. Penalty of INR 15 is imposed on the customer if they don't curtail the load as per the contract (i.e., 20% of the total load). Compared to the baseload, the total load reduction is around 17.5% (4207 MW) and the incentive and penalty costs are INR 1160 and INR 471 respectively. The optimal cost for this case is INR 609843. Further, the peak demand is also reduced by 18%.

In this case, TOU rates are used with a very high price in critical peak periods i.e., INR 8, 10, and 13 at the valley, off-peak, and peak periods respectively. The demand reduction is achieved by about 1764 MW. Hence the total system cost is reduced to INR 597188, which is the lowest cost of all implemented DR cases.

In the CPP program with flat price rates (i.e., 11.6 INR/MWh), the price of INR 40 is considered for critical peak hours (10, 11and 12 h). In this case, the ED solution results in a total cost of INR 603378. Compare to the base case, during critical peak hour load profile is improved and peak reduction of 10.5% is achieved. Total power consumption is also reduced to 22824 MW. (vii) Program 6 (Real Time Pricing)

RTP is the time-based DRP in which no incentive or penalty will be imposed. The price structure considered for RTP is taken from ACOPF. In this case, the total demand for the system is reduced to 22101.

MW, with a system peak of 1267 MW. Here, an energy reduction of 7.68% is attained concerning the base case.

Results from these DRPs have been tabulated in

Program | Operating cost (INR) | Incentive cost (INR) | Penalty cost (INR) | Total cost (INR) |
---|---|---|---|---|

Base | 621378 | 0 | 0 | 621378 |

DLC | 605861 | 1028 | 0 | 606889 |

CAP | 604045 | 1160 | 461 | 604744 |

I/C | 609241 | 1760 | 1118 | 609883 |

TOU | 597188 | 0 | 0 | 597188 |

CPP | 603378 | 0 | 0 | 603378 |

RTP | 619652 | 0 | 0 | 619652 |

Program | Load reduction (MW) | Load reduction (%) | Peak reduction (MW) | Peak reduction (%) |
---|---|---|---|---|

Base | 23940 | -- | 1305 | --- |

DLC | 22161 | 7.43 | 1207 | 7.5 |

CAP | 20827 | 13 | 1128.8 | 13.5 |

I/C | 19733 | 17.57 | 1070.1 | 18 |

TOU | 22176 | 7.36 | 1203.8 | 7.816 |

CPP | 22824 | 4.67 | 1167 | 10.57 |

RTP | 22101 | 7.68 | 1267 | 2.93 |

The simulation results and the influence of considered DRPs on voltage profile and losses are discussed in this section. On application of different DRPs, the variations in load demand bring changes in the voltages at each bus and also the losses associated with it. The volunteered buses for load curtailment are buses 8, 9, 10, 21, 22, 24, 28, and 29. The uniform load curtailment at all volunteered buses (max of 20% of load) is carried out for all DR programs. From the results, it is seen that with load curtailment, voltage profile is improved with the decrease in the losses. The voltage deviations and losses for each DRP are presented in

Program | Maximum voltage deviation (%) | Losses (MW) |
---|---|---|

Base | 1.4 | 52.34 |

DLC | 1.32 | 40.32 |

CAP | 1.24 | 39.657 |

I/C | 1.01 | 38.127 |

TOU | 1.32 | 38.16 |

CPP | 1.22 | 40.57 |

RTP | 1.145 | 39.93 |

Parameters | Base case (MW) | I/C implemented case (MW) |
---|---|---|

P1 | 233.07 | 194.35 |

P2 | 220.82 | 180.85 |

P5 | 147.24 | 82.89 |

P8 | 225.13 | 204.89 |

P11 | 136.01 | 92.02 |

P13 | 130.39 | 132.54/ |

Pwind | 107.97 | 104.56 |

PESS | 111.26 | 83.56 |

Total demand | 1305 | 1070.1 |

Cost (INR/h) | 16937.97 | 16261.3 |

Losses (MW) | 6.86 | 5.85 |

Under the deregulated market different demand response programs have been considered. Demand response is one of the emerging fields in the deregulated market to achieve better load regulation and for the reduction of cost. The economic analysis of different DRPs has been carried out for IEEE 30 bus system. This study shows that the customer's demand is largely influenced by the price elasticity, electricity price, and the incentive and penalty values determined for the concerned DRPs. Results show the reductions in demand of customers under different DRPs. I/C program gives the best result among all DRPs. Incorporating SFLA to solve OPF gives the optimum scheduling to reduce the total cost of the system. Further voltage deviations and losses are also being reduced by the DRP implementation. It is confirmed that by implementing DRPs, a reduction in total cost can be achieved with an enhanced voltage profile.

We thank our management, beloved Principal, and our Head of the department.