Price-Based Demand Response: A Three-Stage Monthly Time-of-Use Tariff Optimization Model
Abstract
:1. Introduction
1.1. Background and Motivation
- How could we divide different time periods?
- How could we set different tariffs for each time period?
- How could we set different tariffs for each season/month?
1.2. Implementation Status of the TOU Tariff
- To improve the TOU tariff pricing mechanism via scientific time period division and establishing a reasonable tariff difference.
- To establish a critical peak tariff mechanism.
- To establish a sound seasonal electricity pricing mechanism.
- Unreasonable peak–valley tariff ratio
- 2.
- Widely implemented critical peak and unpopular deep valley tariffs
- 3.
- Seasonal tariff with lack of narrow time scales
1.3. Literature Review
- Time period division
- 2.
- Price optimization
- 3.
- Execution cycle setting
- The influence of initial clustering center selection bias is harder to mitigate using most applied clustering models of the typical daily load curve, which leads to inaccuracy in the subsequent procedures. Furthermore, the above research only looked at the load curve on the demand side, and insufficient consideration was given to renewable energy generation on the supply side. This is not sufficiently applicable in light of the increasing penetration of renewable energy.
- The current literature does not give enough thought to more refined/high-frequency TOU tariff sets (such as monthly TOU tariffs), which may become increasingly crucial with the continuous development of the power industry. Furthermore, the typical three time periods—peak, flat, and valley—which are the foundation of most current research, need to be further extended for areas with significant renewable energy potential (such as the critical peak and deep valley).
1.4. Contribution and Article Organization
- Rather than using a basic load curve, time period division and price optimization are based on the net load curve, which accomplishes joint consideration for the supply and demand.
- The K-means++ algorithm is adopted to cluster the typical daily net load curve, which can enhance the clustering effect by 4.27–26.70% compared to that of the conventional K-means technique.
- Compared to the earlier annual/seasonal models, the monthly TOU tariff optimization model designed in this paper has a finer particle size. Furthermore, critical peak and deep valley time periods have also been taken into account in order to adapt to the new power system with ever-higher penetration rates of renewable energy than those of the peak, flat, and valley time periods considered in the majority of the previous research.
2. Methodology
2.1. K-Means++ Algorithm
- Step 1: Randomly choose the first cluster center from the net load data .
- Step 2: Calculate the distance from each net load curve to the initial cluster center .
- Step 3: Set other cluster centers based on the probability that each net load curve will be selected as the next cluster center (the net load curve with the maximum probability value). The calculation of probability is as follows:
- Step 4: Repeat the above steps until cluster centers are obtained.
2.2. MBT Approach
- Step 1: Reconstruct the original typical daily net load curve in ascending order as .
- Step 2: Set initial boundaries as and calculate the objective function .
- Step 3: Update the boundaries by moving forward a 1-h step, and calculate the objective function again.
- Step 4: Stop iteration and output the final result when .
- The divided quantity of time periods. The previous studies mostly focused on three periods, i.e., the peak, flat, and valley periods, while this paper has expanded this to five periods.
- The criteria for splitting time periods. Unlike earlier time period division based on the load data, we also use the net load as the primary foundation for time period division, which is similar to typical daily curve clustering.
2.3. PEED Matrix
2.4. TOU Tariff Optimization Model
2.4.1. Objective Function
2.4.2. Constraints
- Cost constraints
- Tariff constraints
- Power constraints
2.5. The Effectiveness Evaluation Index System of the TOU Optimization
- Economic characteristic indicators
- Load characteristic indicators
- Generation characteristic indicators
3. Empirical Analysis
3.1. The Basic Information
3.2. Clustering of the Typical Daily Net Load Curve
3.3. Division of the Time Periods
3.4. TOU Optimization Results
3.5. Effectiveness Evaluation
4. Discussion
4.1. Clustering Effect Comparison
4.2. Time Period Division Comparison
4.3. TOU Tariff Optimization Comparison
4.3.1. Comparison with the Effectiveness before Optimization
4.3.2. Comparison with the Effectiveness Based on Three Traditional Time Periods
5. Conclusions
- As it is different from the previous clustering models, the adopted K-means++ algorithm can overcome the impact of improper selection of initial clustering centers, with an error reduction of 4.27–26.70%. Furthermore, the consideration of the renewable energy generation introduces a more reasonable time period division standard.
- In addition to considering the net load, more extensive and adaptable time division criteria have been proposed to adapt to the increasing penetration rate of renewable energy. Comparing to the three traditional time periods, this paper further introduces the critical peak and deep valley time periods, which enabled a cost saving of 10.36% of the net load cost reduction.
- Using the three-stage monthly TOU tariff optimization model, the monthly net load cost decreased by 2.03–17.27%, the peak–valley difference decreased by 0.31–53.94%, and the maximum monthly abandonment of renewable energy decreased by 1900 MW.
- We must develop more rigorous and extensive time period division methods. With the increase in the installed capacity and penetration rate of renewable energy, the applicability of research based on peak–flat–valley time periods has significantly diminished. In addition, the time periods should not be fixed and evenly divided, which is inconsistent with the actual characteristics of the users and power sources.
- We must develop a reasonable price difference between each time period. At present, the unreasonable TOU tariff has, to some extent, led to the problem of users’ unwillingness to respond. Therefore, it is necessary to further design a users’ response willingness representation model (such as the consumer psychology and PEED matrix ones) and propose more scientific and reasonable electricity price optimization methods.
- We must develop dynamic and comprehensive price adjustment and feedback mechanisms. There is a significant time-varying characteristic in the electricity load and power generation. The traditional electricity pricing mechanisms mostly have annual adjustment cycles, with only a few with quarterly adjustment cycles. In the future, the dynamic adjustment of electricity prices should be more frequent and adapt to the actual demand.
- As analyzed in Section 1.2, different regions implement different TOU tariffs. The time period and price setting are quite different due to resource endowment and user composition. Considering the data limitations, we only took a northwestern province as an example, which should be further expanded to include more regions in the future.
- Although most regions design a TOU tariff for the entire industry, some regions also attempt to break down some typical industries. In the future, we will screen some typical industries (such as those with a high load proportion and strong adjustability) for separate TOU tariff designs.
Author Contributions
Funding
Data Availability Statement
Acknowledgments
Conflicts of Interest
Appendix A
References
- Zhao, Y.; Su, Q.; Li, B.; Zhang, Y.; Wang, X.; Zhao, H.; Guo, S. Have those countries declaring “zero carbon” or “carbon neutral” climate goals achieved carbon emissions-economic growth decoupling? J. Clean. Prod. 2022, 363, 132450. [Google Scholar] [CrossRef]
- Wang, S.; Li, F.; Zhang, G.; Yin, C. Analysis of energy storage demand for peak shaving and frequency regulation of power systems with high penetration of renewable energy. Energy 2023, 267, 126586. [Google Scholar] [CrossRef]
- Liu, Z.; Zhao, S.; Zhang, X.; Tang, Y.; You, G.; Li, J.; Zhao, S.; Hou, X. Renewable energy utilizing and fluctuation stabilizing using optimal dynamic grid connection factor strategy and artificial intelligence-based solution method. Renew. Energy 2023, 219, 119379. [Google Scholar] [CrossRef]
- Wu, Z.; Chen, L.; Wang, J.; Zhou, M.; Li, G.; Xia, Q. Incentivizing the Spatiotemporal Flexibility of Data Centers Toward Power System Coordination. IEEE Trans. Netw. Sci. Eng. 2023, 10, 1766–1778. [Google Scholar] [CrossRef]
- Chen, Y.; Li, J.; Li, H. Tradeoffs in cost competitiveness and emission reduction within microgrid sustainable development considering price-based demand response. Sci. Total Environ. 2020, 73, 135545. [Google Scholar] [CrossRef] [PubMed]
- Zheng, S.; Sun, Y.; Li, B.; Qi, B.; Shi, K.; Li, Y.; Du, Y. Bargaining-based cooperative game among multi-aggregators with overlapping consumers in incentive-based demand response. IET Gener. Transm. Distrib. 2020, 14, 1077–1090. [Google Scholar] [CrossRef]
- McKenna, K.; Keane, A. Residential Load Modeling of Price-Based Demand Response for Network Impact Studies. IEEE Trans. Smart Grid 2016, 7, 2285–2294. [Google Scholar] [CrossRef]
- Grimm, V.; Orlinskaya, G.; Schewe, L.; Schmidt, M.; Zöttl, G. Optimal design of retailer-prosumer electricity tariffs using bilevel optimization. Omega 2021, 102, 102327. [Google Scholar] [CrossRef]
- Choi, D.; Murali, K. The impact of heterogeneity in consumer characteristics on the design of optimal time-of-use tariffs. Energy 2022, 254, 124248. [Google Scholar] [CrossRef]
- Li, R.; Wang, Z.; Gu, C.; Li, F.; Wu, H. A novel time-of-use tariff design based on Gaussian Mixture Model. Appl. Energy 2016, 162, 1530–1536. [Google Scholar] [CrossRef]
- Yuan, X.; Cai, Q.; Song, D. Power Consumption Behavior Analysis Based on Cluster Analysis. In Proceedings of the International Symposium on Artificial Intelligence and Robotics 2021, Shenzhen, China, 25–26 April 2021; Volume 11884, p. 2605823. [Google Scholar]
- Liu, F.; Dong, T.; Hou, T.; Liu, Y. A Hybrid Short-Term Load Forecasting Model Based on Improved Fuzzy C-Means Clustering, Random Forest and Deep Neural Networks. IEEE Access 2021, 9, 59754–59765. [Google Scholar] [CrossRef]
- Li, M.; Xu, D.; Zhang, D.; Zou, J. The seeding algorithms for spherical k-means clustering. J. Glob. Optim. 2020, 76, 695–708. [Google Scholar] [CrossRef]
- He, Y.; Ye, Z.; Wei, W. Research on time division of peak-valley tou tariff considering price elasticity of user demand. Prices Mon. 2022, 12, 1–9. [Google Scholar]
- Bai, X.; Zhang, R.; Xu, L. Research on optimization strategy of time-of-use electricity price for e-commerce sales under the background of electricity reform—Analysis considering user-side and price-based demand response. Price Theory Pract. 2021, 10, 54–57+193. [Google Scholar]
- Gou, J.; Lin, C.; Li, J.; Geng, B.; Li, Z.; Cao, Y.; Bao, Y. Charging Strategy for Electric Vehicles Considering Consumer Psychology and Trip Chain. Front. Energy Res. 2021, 9, 723370. [Google Scholar] [CrossRef]
- Andruszkiewicz, J.; Lorenc, J.; Weychan, A. Seasonal variability of price elasticity of demand of households using zonal tariffs and its impact on hourly load of the power system. Energy 2020, 196, 117175. [Google Scholar] [CrossRef]
- Gholizadeh-Roshanagh, R.; Zare, K. Electric power distribution system expansion planning considering cost elasticity of demand. IET Gener. Transm. Distrib. 2019, 13, 5229–5238. [Google Scholar] [CrossRef]
- Yang, H.; Wang, L.; Ma, Y.; Zhang, D.; Wu, H. Optimization strategy of price-based demand response considering the bidirectional feedback effect. IET Gener. Transm. Distrib. 2021, 15, 1752–1762. [Google Scholar] [CrossRef]
- Zhang, J.; Zhang, Q.; Li, G.; Wu, J.; Wang, C.; Li, Z. Hybrid Model for Renewable Energy and Load Forecasting Based on Data Mining and EWT. J. Electr. Eng. Technol. 2022, 17, 1517–1532. [Google Scholar] [CrossRef]
- de Sá Ferreira, R.; Barroso, L.A.; Lino, P.R.; Carvalho, M.M.; Valenzuela, P. Time-of-Use Tariff Design Under Uncertainty in Price-Elasticities of Electricity Demand: A Stochastic Optimization Approach. IEEE Trans. Smart Grid 2013, 4, 2285–2295. [Google Scholar] [CrossRef]
- Taik, S.; Kiss, B. Clustering-based Optimal Dynamic Pricing for Residential Electricity Consumers. In Proceedings of the 23rd International Conference on Process Control (PC), Strbske Pleso, Slovakia, 1–4 June 2021; pp. 143–148. [Google Scholar]
- Bian, S.; Wang, Z.; Song, W.; Zhou, X. Feature extraction and classification of time-varying power load characteristics based on PCANet and CNN+Bi-LSTM algorithms. Electr. Power Syst. Res. 2023, 217, 109149. [Google Scholar] [CrossRef]
- Venizelou, V.; Philippou, N.; Hadjipanayi, M.; Makrides, G.; Efthymiou, V.; Georghiou, G.E. Development of a novel time-of-use tariff algorithm for residential prosumer price-based demand side management. Energy 2018, 142, 633–646. [Google Scholar] [CrossRef]
- Samadi, M.; Aghamohamadi, M.; Mahmoudi, A. Optimal time period clustering of time-of-use schemes based on elastic loads’ responsiveness. Int. Trans. Electr. Energy Syst. 2020, 30, e12275. [Google Scholar] [CrossRef]
- Xue, W.; Zhao, X.; Li, Y.; Mu, Y.; Tan, H.; Jia, Y.; Zhao, Y. Research on the Optimal Design of Seasonal Time-of-Use Tariff Based on the Price Elasticity of Electricity Demand. Energies 2023, 16, 1625. [Google Scholar] [CrossRef]
- Ahmed, M.; Seraj, R.; Islam, S. The k-means Algorithm: A Comprehensive Survey and Performance Evaluation. Electronics 2020, 9, 1295. [Google Scholar] [CrossRef]
- Dedić, A.; Konjić, T.; Ćalasan, M.; Dedić, Z. Fuzzy C-Means Clustering Applied to Load Profiling of Industrial Customers. Electr. Power Compon. Syst. 2022, 49, 1068–1084. [Google Scholar] [CrossRef]
- Tsekouras, G.; Hatzilau, I.; Prousalidis, J. A new pattern recognition methodology for classification of load profiles for ships electric consumers. J. Mar. Eng. Technol. 2009, 14, 45–58. [Google Scholar] [CrossRef]
- Tsekouras, G.J.; Kotoulas, P.B.; Tsirekis, C.D.; Dialynas, E.N.; Hatziargyriou, N.D. A pattern recognition methodology for evaluation of load profiles and typical days of large electricity customers. Electr. Power Syst. Res. 2008, 78, 1494–1510. [Google Scholar] [CrossRef]
- Sinaga, K.; Yang, M. Unsupervised K-Means Clustering Algorithm. IEEE Access 2020, 8, 80716–80727. [Google Scholar] [CrossRef]
- Arthur, D.; Vassilvitskii, S. K-Means++: The Advantages of Careful Seeding. In Proceedings of the Eighteenth Annual ACM-SIAM Symposium on Discrete Algorithms, SODA 2007, New Orleans, LA, USA, 7–9 January 2007. [Google Scholar]
- Yang, H.; Wang, L.; Zhang, Y.; Tai, H.M.; Ma, Y.; Zhou, M. Reliability Evaluation of Power System Considering Time of Use Electricity Pricing. IEEE Trans. Power Syst. 2019, 34, 1991–2002. [Google Scholar] [CrossRef]
- Yang, H.; Wang, L.; Ma, Y. Optimal Time of Use Electricity Pricing Model and Its Application to Electrical Distribution System. IEEE Access 2019, 7, 123558–123568. [Google Scholar] [CrossRef]
- Andruszkiewicz, J.; Lorenc, J.; Weychan, A. Demand Price Elasticity of Residential Electricity Consumers with Zonal Tariff Settlement Based on Their Load Profiles. Energies 2019, 12, 4317. [Google Scholar] [CrossRef]
- Zhai, Z.; Zhang, L.; Hou, X.; Yang, Q.; Huang, Z. Price elasticity of electricity demand in China A new instrument variable based on marketization policy. Energy Sustain. Dev. 2023, 76, 101275. [Google Scholar] [CrossRef]
Region | Season | Peak | Flat | Valley |
---|---|---|---|---|
California | Summer (1 May to 31 October) | 12:00–18:00 | 8:30–12:00, 18:00–21:30 | 21:30–8:30 |
Winter (1 November to 30 April) | / | 8:30–9:30 | 21:30–8:30 | |
New South Wales | Summer (1 June to 31 August) | 17:00–21:00 | 7:00–17:00, 21:00–22:00 | 22:00–7:00 |
Winter (1 November to 31 March) | 14:00–20:00 | 7:00–14:00, 20:00–22:0 | 22:00–7:00 | |
Other seasons | 14:00–20:00 | 7:00–14:00, 20:00–22:0 | 22:00–7:00 |
Region | Peak–Valley Tariff Ratio | Critical Peak Tariff | Deep Valley Tariff | Seasonal Tariff |
---|---|---|---|---|
Beijing | 4:1 | √ | √ | |
Tianjin | 3.26:1 | √ | √ | |
Northern Hebei | 3:1 | √ | √ | |
Southern Hebei | 5.67:1 | √ | √ | |
Shanxi | 3.56:1 | √ | √ | |
Shandong | 5.67:1 | √ | √ | √ |
Henan | 2.78:1 | √ | √ | |
Hubei | 2.87:1 | √ | √ | |
Hunan | 4:1 | √ | √ | |
Sichuan | 4:1 | √ | √ | |
Chongqing | 4.21:1 | √ | √ | |
Jiangxi | 3:1 | √ | √ | |
Jiangsu | 4.11:1 | √ | √ | √ |
Shanghai | 3.2:1 | √ | √ | |
Zhejiang | 3.57:1 | √ | √ | |
Anhui | 4.15:1 | √ | √ | |
Fujian | 4:1 | √ | √ | |
Heilongjiang | 3:1 | √ | √ | |
Jilin | 3:1 | √ | √ | |
Inner Mongolia | 3:1 | √ | √ | |
Liaoning | 3:1 | √ | √ | |
Shannxi | 3:1 | √ | √ | |
Gansu | 3:1 | |||
Ningxia | 3:1 | |||
Qinghai | 4.21:1 | √ | ||
Xinjiang | 4.71:1 | √ | √ | |
Guangdong | 4.47:1 | √ | √ | |
Guangxi | 3:1 | √ | √ | |
Guizhou | 3:1 | |||
Hainan | 4.25:1 | √ | √ | |
Yunnan | 3:1 | √ | √ |
Peak | Flat | Valley | |
---|---|---|---|
Time period division (hours) | 8:00–10:00 18:00–22:00 | 11:00–17:00 23:00 | 0:00–7:00 |
Tariff (CNY/MWh) | 0.7618 | 0.4304 | 0.1808 |
Critical Peak | Peak | Flat | Valley | Deep Valley | |
---|---|---|---|---|---|
January | 0.6876 | 0.4399 | 0.2808 | 0.234 | 0.1872 |
February | 0.647 | 0.5068 | 0.2854 | 0.2378 | 0.1902 |
March | 0.8551 | 0.7177 | 0.3774 | 0.3145 | 0.2516 |
April | 0.6827 | 0.6605 | 0.5632 | 0.3266 | 0.2612 |
May | 0.7928 | 0.6607 | 0.3482 | 0.2902 | 0.2321 |
June | 0.6968 | 0.5807 | 0.3202 | 0.2668 | 0.2134 |
July | 0.6847 | 0.5753 | 0.3534 | 0.2945 | 0.2356 |
August | 0.8086 | 0.7649 | 0.7587 | 0.3106 | 0.1887 |
September | 0.681 | 0.5675 | 0.3232 | 0.2693 | 0.2154 |
October | 0.6885 | 0.5212 | 0.4931 | 0.2349 | 0.1879 |
November | 0.6885 | 0.5936 | 0.4445 | 0.2951 | 0.2361 |
December | 0.698 | 0.663 | 0.3112 | 0.2593 | 0.2075 |
Month | I1 (CNY) | I2 (MW) | I3 (MW) | I4 (MW) | I5 | I6 (MW) |
---|---|---|---|---|---|---|
January | 96,452 | 9223 | 97 | 9126 | 98.95% | 0 |
February | 78,675 | 8685 | 0 | 8685 | 100.00% | 0 |
March | 74,804 | 7595 | −132 | 7727 | 101.74% | 282 |
April | 65,580 | 7234 | 0 | 7234 | 100.00% | 0 |
May | 65,207 | 7502 | 3153 | 4349 | 57.97% | 0 |
June | 76,140 | 7599 | 2628 | 4971 | 65.42% | 0 |
July | 78,093 | 8152 | 536 | 7616 | 93.42% | 0 |
August | 78,350 | 6171 | 3163 | 3008 | 48.75% | 0 |
September | 84,759 | 8567 | 1959 | 6608 | 77.13% | 0 |
October | 100,206 | 9295 | 2188 | 7107 | 76.46% | 0 |
November | 97,851 | 9159 | 638 | 8521 | 93.04% | 0 |
December | 101,128 | 9258 | 675 | 8583 | 92.70% | 0 |
K-Means++ | K-Means | |
---|---|---|
5.07 × 108 | 5.92 × 108 | |
5.65 × 108 | 6.70 × 108 | |
6.52 × 108 | 7.13 × 108 | |
5.83 × 108 | 6.43 × 108 | |
3.32 × 108 | 4.11 × 108 | |
4.11 × 108 | 4.38 × 108 | |
6.27 × 108 | 7.18 × 108 | |
6.06 × 108 | 6.33 × 108 | |
5.74 × 108 | 6.08 × 108 | |
6.44 × 108 | 6.97 × 108 | |
5.40 × 108 | 6.41 × 108 | |
5.71 × 108 | 7.79 × 108 |
Month | I1 (CNY) | I2 (MW) | I3 (MW) | I4 (MW) | I5 | I6 (MW) |
---|---|---|---|---|---|---|
January | 109,418 | 9465 | 120 | 9345 | 98.73% | 0 |
February | 95,094 | 9007 | −465 | 9472 | 105.16% | 1622 |
March | 80,146 | 8872 | −189 | 9061 | 102.13% | 404 |
April | 69,192 | 8196 | −487 | 8683 | 105.94% | 1900 |
May | 78,400 | 8469 | 2808 | 5661 | 66.84% | 0 |
June | 86,626 | 8441 | 1952 | 6489 | 76.87% | 0 |
July | 86,263 | 8421 | 782 | 7640 | 90.72% | 0 |
August | 81,977 | 8149 | 1618 | 6530 | 80.14% | 0 |
September | 96,053 | 8942 | 1674 | 7268 | 81.28% | 0 |
October | 113,189 | 9735 | 1267 | 8468 | 86.99% | 0 |
November | 107,729 | 9510 | 800 | 8710 | 91.59% | 0 |
December | 103,219 | 9856 | 563 | 9293 | 94.29% | 0 |
Month | I1 (CNY) in Five Time Periods | I1 (CNY) in Three Time Periods |
---|---|---|
January | 96,452 | 110,712 |
February | 78,675 | 97,320 |
March | 74,804 | 81,609 |
April | 65,580 | 74,369 |
May | 65,207 | 78,837 |
June | 76,140 | 80,067 |
July | 78,093 | 82,875 |
August | 78,350 | 75,566 |
September | 84,759 | 97,592 |
October | 100,206 | 106,547 |
November | 97,851 | 112,245 |
December | 101,128 | 114,737 |
Sum | 997,245 | 1,112,476 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
You, P.; Li, S.; Li, C.; Zhang, C.; Zhou, H.; Wang, H.; Zhao, H.; Zhao, Y. Price-Based Demand Response: A Three-Stage Monthly Time-of-Use Tariff Optimization Model. Energies 2023, 16, 7858. https://doi.org/10.3390/en16237858
You P, Li S, Li C, Zhang C, Zhou H, Wang H, Zhao H, Zhao Y. Price-Based Demand Response: A Three-Stage Monthly Time-of-Use Tariff Optimization Model. Energies. 2023; 16(23):7858. https://doi.org/10.3390/en16237858
Chicago/Turabian StyleYou, Peipei, Sitao Li, Chengren Li, Chao Zhang, Hailang Zhou, Huicai Wang, Huiru Zhao, and Yihang Zhao. 2023. "Price-Based Demand Response: A Three-Stage Monthly Time-of-Use Tariff Optimization Model" Energies 16, no. 23: 7858. https://doi.org/10.3390/en16237858
APA StyleYou, P., Li, S., Li, C., Zhang, C., Zhou, H., Wang, H., Zhao, H., & Zhao, Y. (2023). Price-Based Demand Response: A Three-Stage Monthly Time-of-Use Tariff Optimization Model. Energies, 16(23), 7858. https://doi.org/10.3390/en16237858