Abstract
With the increasing penetration of new-type loads such as electric vehicles and hydrogen fuel vehicles in urban power grids, the peak-to-valley load difference increases sharply, and a multi-energy coordination model is desirable. This article proposes a day-ahead operation model of an urban energy system considering traffic flows and peak shaving, which can positively contribute to multi-energy complement and low-carbon emission. The proposed model minimizes the total cost of electricity and gas by optimizing the charging and discharging strategies of energy storage, in which the output of the wind turbine and energy management of the energy hub are adaptively adjusted. The urban energy system is represented by a second-order cone (SOC) energy flow model, and hence, the optimization problem is modeled as a mixed integer SOC programming (MISOCP). Finally, test results on an integrated urban energy network indicate that the energy storage and multi-energy coordination can alleviate the peak load cutting and valley filling. The relationship between urban grid operation cost and peak-valley difference is also discussed. The maximum utilization of renewable energy sources using gasoline vehicles has been presented in this study to illustrate cost and emission reductions for a sustainable integrated electricity and transportation infrastructure.
1 Introduction
1.1 Motivation
In recent years, energy consumption and global environmental problems have become increasingly serious, and the energy transition toward a low-carbon energy system is highly desirable (ZhaoĀ etĀ al.,Ā 2017). The energy internet provides great potential for reliable power supply and an improved energy efficiency by integrating various energy production/conversion components, including renewable energy sources (RESs), energy storage systems (ESSs), hydrogen fuel vehicles (HVs), and electric vehicles (EVs) (StrasserĀ etĀ al.,Ā 2015; Meng and Wang,Ā 2017; Yang and Fang,Ā 2017). At the same time, interactions between urban electrical, transportation, and natural gas systems have become increasingly common due to the utilization of coupling components such as microturbines, combined heat and power (CHP), and power to gas units (Mancarella,Ā 2014; Shabanpour-Haghighi and Seifi,Ā 2015). The coupling of urban multi-energy networks positively contributes to the improved energy utilization efficiency, multi-energy complement, and the construction of a low-carbon sustainable energy system (XuĀ etĀ al.,Ā 2020).
Traditionally, coupled energy systems have been operated and optimized independently (ChenĀ etĀ al.,Ā 2014). The urban distribution network that is inactive generally has limited ability toregulate power flow, which results in the problem of a large difference between the daily peak and valley of the net load (UddinĀ etĀ al.,Ā 2018). On the other hand, the traditional energy storage has relatively small storage capacity with a high investment cost, and hence, it is difficult to store energy effectively on a large scale and for a long time (ZhaoĀ etĀ al.,Ā 2020). Considering the fact that the traditional urban distribution network cannot fully meet the electricity demand of new-type loads such as EVs and HVs, the development of day-ahead co-optimization of multi-energy systems is of considerable interest.
Owing to the aforementioned considerations, the present work develops a day-ahead co-optimization strategy for urban energy networks considering traffic flows.
1.2 Literature Review
1.2.1 Research on Power-Transportation Coordination
EVs and HVs have become a promising alternative to gasoline vehicles (GVs) for decarbonizing the transportation sector and combating climate change (IEA,Ā GlobalĀ EVĀ Outlook,Ā 2020). An integrated demand response framework is designed to regulate flexible resources in both networks through optimal road tolls and electricity tariffs in LvĀ etĀ al.,Ā 2021a. In Rotering and IlicĀ (2011), the problem of growing peak load and grid overloading is addressed based on a forecast of future electricity prices, and the dynamic programming algorithm is employed to find the economically optimal solution for the vehicle owner. HuĀ etĀ al.Ā (2016) present an integrated optimization framework for battery sizing, charging, and on-road power management in plug-in hybrid EVs. The maximum utilization of RESs using GVs has been presented in Saber and VenayagamoorthyĀ (2011) to achieve cost and emission reductions, which contributes to a low-carbon integrated electricity and transportation infrastructure. We note that most works consider the coupling between the transportation system and the power grid. However, the coupling between the power, natural gas, and the transportation system has not been well considered.
1.2.2 Coordinated Optimization of Multi-Energy Systems
A day-ahead optimal scheduling of the urban energy system is proposed by JinĀ etĀ al.,Ā 2016, in which the flexible and reconfigurable topology of power distribution networks is considered. The joint and unified optimization decision of coupled power systems and natural gas systems is analyzed in MartĆnez-MaresĀ etĀ al.Ā (2011). Correa-Posada and SĆ”nchez-MartĆnĀ (2015) determine the dynamic optimal operation strategy of the electricityāgas integrated energy system, which includes a transient natural-gas flow model. An integrated electricityāgas system with steady-state energy flow models considering the uncertainty of wind power is proposed in Gao and LiĀ (2020). Additionally, a new Benders decomposition-based algorithm is proposed, which improves the solution efficiency for non-convex models (Gao and Li,Ā 2021). A new external dependency model based on the energy hub (EH) is introduced to consider the possible uncertainty in customer decision through stochastic model processing (NeyestaniĀ etĀ al.,Ā 2015). The stochastic optimal model is investigated with the comprehensive consideration of renewable generation and carbon-capture-based power-to-gas technology (LiĀ etĀ al.,Ā 2018). However, we note that most existing works on coordination optimization of integrated energy systems do not consider the existence of new-type loads such as EVs and HVs. Moreover, in most existing studies, a common assumption is that the driving mode of EVs is known. Arrival rates and times as well as the number of charging requests are specified in advance in a deterministic or stochastic manner, which can be determined from data-driven methods or queuing theory or modeled with probabilistic models through Markov decision processes. This assumption is reasonable for set-level research, such as residential parking management or a single CS management. However, it may not be appropriate if the system-level interdependency is under investigation, which is proposed in this article.
1.3 Contribution of This Study
To fill the research gaps summarized above, our model makes the following contributions:
1) We propose a mixed-UE traffic flow model based on an urban transportation system and divide the overall vehicular flows into five categories: GV traffic flow, EV traffic flow, HV traffic flow, HV charging flow, and EV charging flow. Our model, hence, accurately simulates the distribution of EVs and HVs in the urban transportation network.
2) We develop a day-ahead optimal dispatch of an urban energy system considering the transportation-network operation model, which includes the electricāgas co-optimization and the impact of traffic flows. Our model positively contributes to multi-energy complement and low-carbon operation. Moreover, the optimal energy management of EHs reduces the peak-valley difference caused by increasing penetrations of new-type loads such as EVs and HVs.
The remainder of this article is organized as follows: SectionĀ 2 elaborates on the mathematical model of day-ahead operation for an urban energy system considering traffic flows. SectionĀ 3 presents the performance of the proposed optimization model. Finally, SectionĀ 4 presents the conclusions of the work.
2 Day-Ahead Operation of an Urban Energy System Considering Peak Load Shaving
This section first presents the day-ahead electricityāgas co-optimization for urban energy systems, which includes the operating constraint with regard to peak load shaving. Then, the joint optimal dispatch model considers the electricity charging load and hydrogen charging load on the traffic side. Finally, through the optimal control of ESS and energy management of EH, the peak shaving and valley filling of the urban system with high penetrations of wind power are realized. The multi-energy systems are shown in FigureĀ 1. There are many nodes connected to EH in the urban distribution network. EH provides hydrogen for the hydrogen charging station on the traffic side through the internal P2H device, meets the heat load demand through the CHP unit, and provides power for the charging station on the traffic side simultaneously through the converter and CHP unit.
FIGURE 1
Five key assumptions are made in our work:
1) We use a nonlinear day-ahead optimal dispatch model of an urban energy system that considers the transportation-network operation model.
2) We use fixed natural gas price information to coordinate the operation of electricity and natural gas systems. Unfixed gas pricing often leads to market inefficiencies.
3) We assume that the wind power output is deterministic, without considering the randomness of wind power output.
4) We use time-of-use price to calculate the cost of purchasing electricity in an urban energy system.
2.1 Branch Flow Model of the Urban Power Distribution Network
Since AC power flow constraints are non-convex, the day-ahead dispatch of the urban energy system is a non-convex optimization problem, which is difficult to solve. The second-order cone (SOC) relaxation based on the branch flow model is implemented to improve its computational tractability. Therefore, we use the branch flow model to replace the nonlinear power flow model of the urban distribution network (Baran and Wu,Ā 1989):
EqsĀ 1, 2 replace the voltage and current quadratic terms in the original nonlinear power flow equation with new variables. The DistFlow model for radial electricity networks is formulated by 3ā7, in which variables of the squared voltage magnitude and squared branch current are employed. EqsĀ 3ā5 represent the active and reactive power flow balance. EqsĀ 6, 7 represent the injected active and reactive power of the node.
The operation of the urban AC distribution network should satisfy the following security constraints:
Constraint (8) limits the transmission power capacity of the root node. Constraint (9) represents the feeder capacity constraint. Then, curtailment constraints with respect to distributed generation are given in 10, 11. Constraint (12) limits the operating range of the squared voltage magnitude.
ESS operation constraints include charge and discharge state constraints, storage capacity constraints, charge and discharge power constraints, and daily allowable charge and discharge times constraints (Gabash and Li,Ā 2012), which are provided as follows:
Constraint (13) denotes the charging and discharging states of ESS, and constraint (14) represents upper limit of charge and discharge. EqsĀ 15, 16 are energy constraints of ESS and daily allowable charge and discharge times constraints, respectively. Then, the quadratic equality (5) is further relaxed into the following SOC inequality (Taylor and Hover,Ā 2012; Farivar and Low,Ā 2013):
where is the mathematical expression of the Euclid norm.
2.2 Model of the Energy Hub
The EH (ChenĀ etĀ al.,Ā 2019) depicted in FigureĀ 1A represents the coupling between electricity and gas networks. CHP represents the key component of this EH, which generates electricity and heat simultaneously, and the urban distribution network provides the charging service for EVs through the EH. The energy conversion efficiency between the input and output energy flows of the EH is formulated in 18
Eq.Ā 19 represents the dispatch factor of the EH input gas flow
Then, the EH also needs to satisfy the following operating constraints:
Constraint (20) gives the range of total gas purchases of EH. Constraints (20) and (21) represent the production capacity of CHP units and gas furnaces, respectively.
For a specific case, the nodes of the distribution network are connected to some components of EH, for example, only CHP units or gas furnaces, as shown in FigureĀ 2. On the other hand, the power distribution network converts electric power into hydrogen through P2H devices to provide hydrogen supply services for hydrogen-fueled vehicles on the traffic side (Korpas and Holen,Ā 2006; PanĀ etĀ al.,Ā 2021). The structure of P2H is shown in FigureĀ 3.
FIGURE 2
FIGURE 3
The energy conversion efficiency between the input and output energy flows of P2H is given in 23 and 24
Constraint (23) calculates the energy conversion efficiency between the input and output energy flows of P2H. Constraint (24) limits the hydrogen production.
2.3 Transportation System Modeling
The user equilibrium (UE) model has been widely used to describe the urban transportation system in transportation research. In the case of the given traffic network and traffic demand, the UE model can output the current traffic state, namely, the traffic flow distribution on the road. The transportation system model in this article is based on the mixed-UE traffic flow model (WeiĀ etĀ al.,Ā 2018). We divide the overall vehicular flows into five categories: GV traffic flow, EV traffic flow, HV traffic flow, HV charging flow, and EV charging flow. The difference between traffic flow and charging flow lies in that the former represents the route choice of EVs and HVs, while the latter represents the charging choice of EVs and HVs. Traffic link has regular link, charging link, and bypass link. The charging links correspond to EV and HV queuing and charging events in EVCSs and HVCSs, and bypass links denote bypass events at EVCSs and HVCSs (LvĀ etĀ al.,Ā 2021b; TengĀ etĀ al.,Ā 2021).
The improved mixed-UE transportation model is a nonlinear complementary problem (NCP) consisting of four constraint sets
where constraint sets are presented as follows:
Cons āĀ Flow:
Cons āĀ Time:
Cons āĀ Cost:
Cons āĀ Cp:
Constraint set Cons āĀ Flow describes flow conservation in the TN. Specifically, constraints (29)ā(31) indicate that the traffic flow on each link is equal to the sum of the flows on all paths that pass through the link. Constraints (26)ā(28) explain that the traffic flow on all paths connecting each O-D pair is equal to the corresponding traffic demand. Constraint set Cons āĀ Time gives link travel time and estimation of different link types, respectively. For regular link , constraint (33) uses the Bureau of Public Roads (BPR) function [BureauĀ ofĀ PublicĀ RoadsĀ (1,Ā 1964)] to represent the link travel time as a function of . In constraint (7), the first and second terms denote the charging time and queuing time, respectively. Constraint set Cons āĀ Cost describes the travel costs of GVs, EVs, and HVs that choose path Ļ. Note that the cost for EVs and HVs contains both travel cost and charging cost, respectively. Constraint set Cons āĀ Cp represents the complementarity condition of the Wardrop UE principle (Sheffi,Ā 1985), which indicates that no GV/EV/HV travelers could decrease their travel costs by unilaterally switching their route choices.
By changing the electricity price and hydrogen price for each period of time, we calculate the transportation network model 24 times. Then, we can get the load curve of each charging station and hydrogen charging station in a day according to EqsĀ 42, 43. In addition, the heat load can be obtained from the literature (Soroudi,Ā 2017)
2.4 The Day-Ahead Optimization Dispatch Model of an Urban Energy System
With the goal of minimizing the sum of gas purchase cost and electricity purchase cost, we established a day-ahead optimization dispatch model of an urban energy system, and the optimization model is summarized below (LuĀ etĀ al.,Ā 2020; Bakeer and Hossam,Ā 2021; ChengĀ etĀ al.,Ā 2021; LuĀ etĀ al.,Ā 2020; Bakeer and Hossam,Ā 2021; ChengĀ etĀ al.,Ā 2021).
The objective function is provided as
The objective function (44) is subject to Day-ahead operation constraints:
FigureĀ 4 shows the framework of the day-ahead optimization dispatch model of an urban energy system, which clarifies the interactions between power and transportation networks from the modeling perspective.
FIGURE 4
3 Case Study
The urban energy system shown in FigureĀ 5 is composed of a transportation network and a 33-node urban electricity network with five EHs located at nodes 6, 16, 18, 21, and 30 and two P2H units at nodes 8 and 10. CS1āCS8 are charging stations which correspond to nodes 6, 16, 12, 27, 18, 25, 21, and 30 in the power network, respectively. HFS1āHFS7 are hydrogen charging stations which correspond to nodes 6, 16, 21, 30, 18, 10, and 8 in the power network, respectively. In addition, two 0.4-MW capacity wind farms are connected at nodes 14 and 33, while two 1.8-MW capacity ESSs are connected at nodes 11 and 29, respectively. The computational tests were conducted on a laptop with Intel Core 1.0Ā GHz CPU and 16Ā GB memory, and the MISOCP model was programmed in GAMS and solved with the GUROBI solver.
FIGURE 5
3.1 Electric Vehicle Charging Load and Hydrogen Load
FigureĀ 6 shows the daily heat load, the charging load of EVs, and the hydrogen load of HVs on the traffic side. As can be observed from this figure, the charging load curve of EVs has the same trend with the hydrogen load curve of HVs, and both have the characteristic of double peaks. During the morning rush hour from 8:00 to 10:00 and the evening rush hour from 17:00 to 19:00, the electric load and hydrogen load on the traffic side are relatively high, while the heat load is higher during the day and lower at night (YeĀ etĀ al.,Ā 2021). The trend of these three loads is consistent with that of the total load.
3.2 The Anti-Peak-Shaving Characteristic of Wind Power Production
FigureĀ 7 shows the comparison of the wind power output curve, total original load curve, and total net load curve without considering energy storage. In this figure, load 0 represents the total original load curve, and load 1 represents the total net load curve without considering the energy storage. As can be observed from the figure, wind turbines produce more output at night, which results in a lower net load at the valley time. However, the output of wind turbines is less in the daytime, and hence, it cannot significantly reduce the peak load. Consequently, the daily wind power production results in a larger peak-to-valley difference of net load.
FIGURE 6
3.3 Coordination Between ESS and Wind Power Production
It is noteworthy that when the ESS participates in the power regulation of the power grid, the disadvantages brought by the anti-peak-shaving of wind power can be well compensated. As shown in FigureĀ 8, load 3 represents the total net load considering energy storage and multi-energy regulation, and load 2 represents the total net load without considering multi-energy regulation. ESS charges at night when the power load is low and the wind power output is large and discharges during the day when the load is relatively high and the wind power output is small. As a result, peak cutting and valley filling of the power load are achieved.
FIGURE 7
3.4 Complementary Operation of Multi-Energy Coupling
Natural gas is simultaneously converted into heat and electricity through CHP units and only into heat through the gas furnace. FigureĀ 9 shows the total purchase of natural gas at each time period, from which we can observe that the amount of gas purchase of gas is larger at the time period when the power load is high than that at the time period when the power load is low. As shown in FigureĀ 8, the difference between curve load 2 and load 3 is mainly due to the multi-energy complement. In addition, the EV load is relatively high at time periods from 7:00 to 17:00. To summarize, increasing the power supply of CHP units can effectively reduce the peak load, but this impact on valley filling is not obvious.
FIGURE 8
3.5 Impact of ESS Strategies on Urban Energy System Operation
As discussed above, ESS can play an important role in power systems. Reasonable ESS scheduling decisions can alleviate the impact of wind power on the distribution network and promote peak shaving and valley filling. FigureĀ 10 and FigureĀ 11 show the chargeādischarge curve of ESS and the SOC state of charge of ESS, respectively. We can note that for time periods from 0:00 to 7:00, the net load is relatively with low power consumption and high wind power production, and ESS charges to satisfy the load electricity demand. The SOC state of charge rises in preparation for the next stage of the discharge process. For time periods from 8:00 to 11:00, the electricity price is high and the net load is relatively high. ESS reduces the load peak by discharging. From 16:00 to 19:00, the load reaches the peak again, and ESS continues to discharge, which reduces the power purchase cost from the main grid. At the same time, in order to ensure the normal operation of the next charge and discharge cycle, the ESS charges and returns to the original state of charge at time periods from 20:00 to 24:00.
FIGURE 9
FIGURE 10
FIGURE 11
We introduce an index ΓPa to measure the average peak-valley difference:
TableĀ 1 lists mean peak-valley differences and the total operation costs under different ESS capacities. As shown in this table, the increase in ESS capacity positively reduces the average peak-valley difference of power demands and reduces the total operating costs. This is because a higher energy storage capacity can provide more sufficient charging power in load valley and discharging power in load peak, which also reduces the purchase of power at the peak time, thus reducing the operating cost.
TABLE 1
| Capacity of ESS/MW | 1.8 | 2.2 | 2.6 | 2.8 | 3.2 |
| Mean peak-valley difference/MW | 1.37 | 1.27 | 1.16 | 1.11 | 1.00 |
| Cost/$ | 25145.35 | 25029.00 | 24919.91 | 24868.68 | 24778.03 |
Average peak-valley difference under different ESS capacities.
3.6 Cost Comparison Under Different Values of Average Peak-Valley Difference
To control the load peak-valley difference to an ideal level, we add an additional linear constraint (46) into the original model. Meanwhile, we set the energy storage capacity to be fixed at the original valueĪ» denotes the load peak-valley difference which we want to get. By modifying the value of Ī», we can obtain the corresponding operating costs under different requirements of average peak-valley difference, and the results are summarized in TableĀ 2.
TABLE 2
| Ī»/MW | 1.37 | 1.27 | 1.15 | 1.10 | 1.05 |
| Cost/$ | 25145.35 | 25282.07 | 25731.15 | 25952.80 | 26179.76 |
Cost comparison under different constraint values of average peak-valley difference.
Generally, a more ideal peak-valley difference increases the total operating cost, which indicates that the operating cost and the requirement of peak shaving and valley filling are contradictory. The reason is that if we set the value of Ī» to be small, we need to increase the net load during the valley time. Therefore, part of the wind power will be curtailed during the load valley, and additional electricity needs to be purchased from the main grid, which reduces the peak-to-valley difference and inevitably increase the total operating cost.
4 Conclusion
We have presented a day-ahead co-optimization approach of unban energy systems considering traffic flows and peak shaving. Based on the numerical results obtained, we can draw the following conclusions:
⢠Optimal charging and discharging dispatch strategies of energy storages can alleviate the peak-valley difference caused by the anti-peak-shaving characteristic of wind power production.
⢠In addition, multi-energy complement can significantly reduce the peak-valley difference and contribute to the improved energy utilization efficiency. However, increasing the supply of natural gas can reduce the peak load, but the effect on valley filling is not obvious.
⢠When the ESS capacity is fixed, the higher the peak shaving and valley filling requirements, the greater the operating cost.
Our work, hence, might help the operators of urban energy systems to produce operating decisions of multi-energy production sources and satisfy the requirement from multi-energy demands. Our work might also help the design of a future low-carbon energy system. Meanwhile, in future research, we will consider the optimal dispatch model of the transportation system into the day-ahead dispatch model of the urban energy system to further improve social benefits. Reasonable traffic management strategies can give full play to the initiative of the urban distribution network to achieve peak shaving and valley filling and economic optimization.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.
Author contributions
YP: methodology, software, validation, formal analysis, writingāoriginal draft, and visualization. JM: methodology, investigation, and writingāreview editing. CL conceptualization, methodology, and writingāreview editing. SC conceptualization, resources, and supervision. LP methodology, validation, and writingāreview editing.
Funding
This work was supported by State Grid Zhejiang Electric Power Company Project āSteady-state models and parameter identification technology of energy storage stationsā.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisherās note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.
Glossary
Indices and sets
- i, j
Subscript indices of nodes in electricity networks
- t
Subscript indices of time periods
- min
Superscript index of minimum value
- max
Superscript index of maximum value
- Ωnode
Set of electrical nodes
- Ωline
Set of branches
- Ωsub
Set of root nodes
- ΩESS
Set of (energy storage system) ESS nodes
- ΩDGN
Set of DG nodes
- T
Operational cycle.
Set of regular links in a transportation network.
Set of charging links in the transportation network
- TRS
Set of originādestination (O-D) pairs
- TR
Set of origin nodes
- TS
Set of destination nodes.
Set of gasoline vehicles/electric vehicles/hydrogen fuel vehicles (GV/EV/HV) paths connecting an O-D pair r-s, where r āĀ TRs āĀ TS and rs āĀ TRS
Variables
- Ui,t
Voltage of electrical nodes
- Iij,t
Current magnitude through branch ij
- Pij,t, Qij,t
Active and reactive power flow through branch ij
The injected active and reactive power of node i
Output active and reactive power of the root node
Output active and reactive power of the distributed generation
Active and reactive power load
Active power charge and discharge of ESS
Input power flow of the energy hub (EH)
Input power flow of the gas furnace (GF)
Input power flow of combined heat and power (CHP)
Input power flow of power to hydrogen (P2H)
Binary variable
Energy level of the ESS
Dispatched electricity and heat demand
Input gas flow of the energy hub
- v
Percentage of gas fed into combined heat and power
- Hel,t
Hydrogen production
GV flow on path between O-D pair r-s
EV flow on path between O-D pair r-s
HV flow on path between O-D pair r-s
Aggregated traffic flow on link
Aggregated traffic flow on link
Aggregated traffic flow on link
Travel time on link
Average time that EVs spend on link
Average time that HVs spend on link
GV travel cost on path between O-D pair r-s
EV travel cost on path between O-D pair r-s
HV travel cost on path between O-D pair r-s
Minimal travel cost of GVs between O-D pair r-s
Minimal travel cost of EVs between O-D pair r-s
Minimal travel cost of HVs between O-D pair r-s
Peak and valley net load
Parameters
- Rij,t
The resistance of branch ij
- Xij,t
The reactance of branch ij
- tanāĻ
Power factor of DG output
Upper limits of ESS charge and discharge power
- α, β
Charge and discharge efficiency coefficient of ESS
Maximum storage capacity of ESS
Maximum charging and discharging times of ESS
- Ī·ee
Transformer efficiency of the energy hub
- Ī·CHP,e
Electrical efficiency of combined heat and power
- Ī·CHP,h
Thermal efficiency of combined heat and power
- Ī·gh
Gas furnace efficiency
- Ī·e2hy
Electrical efficiency of P2H
- Hc
The higher heating value of hydrogen
- Ļg, Ļe, Ļh
Ratio of GV/EV/HV traffic demand
If path passes link , ; otherwise, .
If path passes link , ; otherwise, .
If path passes link , ; otherwise, .
Free flow travel time on link
Traffic flow capacity of link
Maximum allowable vehicular flow of charging link
- pser,e, pser,h
Average service rate at charging stations
- Ee, Eh
Charging demand of unit traffic flow.
Maximum waiting time at charging stations
- qrs,t
Trip rate (traffic demand) between O-D pair r-s
- Ļ
Monetary cost of travel time.
Charging price at charging stations
- μe,t
Time of use price
- μg,t
Unit price of natural gas
- Npeak, Nvalley
Number of peak and valley periods
References
1
BakeerA.HossamS. (2021). Salama and Istvan Vokony. Integration of PV System with SMES Based on Model Predictive Control for Utility Grid Reliability Improvement. Prot. Control. Mod. Power Syst.6 (2), 191ā203. 10.1186/s41601-021-00191-1
2
BaranM. E.WuF. F. (1989). Network Reconfiguration in Distribution Systems for Loss Reduction and Load Balancing. IEEE Trans. Power Deliv.4 (2), 1401ā1407. 10.1109/61.25627
3
Bureau of Public Roads (1964). Traffic Assignment Manual. Washington, D.C.: U.S. Department of Commerce.
4
ChenS.WeiZ.SunG.CheungK. W.WangD.ZangH. (2019). Adaptive Robust Day-Ahead Dispatch for Urban Energy Systems. IEEE Trans. Ind. Electron.66 (2), 1379ā1390. 10.1109/tie.2017.2787605
5
ChenY.GribikP.GardnerJ. (2014). Incorporating Post Zonal Reserve Deployment Transmission Constraints into Energy and Ancillary Service Co-optimization. IEEE Trans. Power Syst.29 (2), 537ā549. 10.1109/tpwrs.2013.2284791
6
ChengYan.TangYi.DaiJianfeng. (2021). Uncertainty Modeling of Wind Power Frequency Regulation Potential Considering Distributed Characteristics of Forecast Errors. Prot. Control. Mod. Power Syst.6 (3), 276ā288.
7
Correa-PosadaC. M.SĆ”nchez-MartĆnP. (2015). Integrated Power and Natural Gas Model for Energy Adequacy in Short-Term Operation. IEEE Trans. Power Syst.30 (6), 3347ā3355. 10.1109/tpwrs.2014.2372013
8
FarivarM.LowS. H. (2013). Branch Flow Model: Relaxations and Convexification-Part I. IEEE Trans. Power Syst.28 (3), 2554ā2564. 10.1109/tpwrs.2013.2255317
9
GabashA.LiP. (2012). Active-Reactive Optimal Power Flow in Distribution Networks with Embedded Generation and Battery Storage. IEEE Trans. Power Syst.27 (4), 2026ā2035. 10.1109/tpwrs.2012.2187315
10
GaoH.LiZ. (2021). A Benders Decomposition Based Algorithm for Steady-State Dispatch Problem in an Integrated Electricity-Gas System. IEEE Trans. Power Syst.36 (4), 3817ā3820. 10.1109/tpwrs.2021.3067203
11
GaoH.LiZ. (2020). āEffect of Power to Gas on Integrated Electricity-Gas System with Uncertain Wind Generation,ā in IEEE 3rd Student Conference on Electrical Machines and Systems. SCEMS, 164ā168. 10.1109/scems48876.2020.9352402
12
HuX.MouraS. J.MurgovskiN.EgardtB.CaoD. (2016). Integrated Optimization of Battery Sizing, Charging, and Power Management in Plug-In Hybrid Electric Vehicles. IEEE Trans. Contr. Syst. Technol.24 (3), 1036ā1043. 10.1109/tcst.2015.2476799
13
IEA, Global EV Outlook (2020). Entering the Decade of Electric Drive?International Energy Agency.
14
JinX.MuY.JiaH.WuJ.XuX.YuX. (2016). Optimal Day-Ahead Scheduling of Integrated Urban Energy Systems. Appl. Energ.180, 1ā13. 10.1016/j.apenergy.2016.07.071
15
KorpasM.HolenA. T. (2006). Operation Planning of Hydrogen Storage Connected to Wind Power Operating in a Power Market. IEEE Trans. Energ. Convers.21 (3), 742ā749. 10.1109/tec.2006.878245
16
LiY.ZouY.TanY.CaoY.LiuX.ShahidehpourM.et al (2018). Optimal Stochastic Operation of Integrated Low-Carbon Electric Power, Natural Gas, and Heat Delivery System. IEEE Trans. Sustain. Energ.9 (1), 273ā283. 10.1109/tste.2017.2728098
17
LuS.GuW.MengK.YaoS.LiuB.DongZ. Y. (2020). Thermal Inertial Aggregation Model for Integrated Energy Systems. IEEE Trans. Power Syst.35 (3), 2374ā2387. 10.1109/tpwrs.2019.2951719
18
LvS.ChenS.WeiZ.ZhangH. (2021). Power Transportation Coordination: toward a Hybrid Economic-Emission Dispatch Model. IEEE Trans. Power Syst.10.1109/tpwrs.2021.3131306
19
LvS.WeiZ.ChenS.SunG.WangD. (2021). Integrated Demand Response for Congestion Alleviation in Coupled Power and Transportation Networks. Appl. Energ.283, 116206. 10.1016/j.apenergy.2020.116206
20
MancarellaP. (2014). MES (Multi-energy Systems): An Overview of Concepts and Evaluation Models. Energy65, 1ā17. 10.1016/j.energy.2013.10.041
21
MartĆnez-MaresA.Fuerte-EsquivelC. R.de IngenierĆaI. (2011). Integrated Energy Flow Analysis in Natural Gas and Electricity Coupled Systems. North American Power Symposium, 1ā7.
22
MengW.WangX. (2017). Distributed Energy Management in Smart Grid with Wind Power and Temporally Coupled Constraints. IEEE Trans. Ind. Electron.64 (8), 6052ā6062. 10.1109/tie.2017.2682001
23
NeyestaniN.Yazdani-DamavandiM.Shafie-khahM.ChiccoG.CatalaoJ. P. S. (2015). Stochastic Modeling of Multienergy Carriers Dependencies in Smart Local Networks with Distributed Energy Resources. IEEE Trans. Smart Grid6 (4), 1748ā1762. 10.1109/tsg.2015.2423552
24
PanG.HuQ.GuW.DingS.QiuH.LuY. (2021). Assessment of Plum Rain's Impact on Power System Emissions in Yangtze-Huaihe River basin of China. Nat. Commun.12, 6156. 10.1038/s41467-021-26358-w
25
RoteringN.IlicM. (2011). Optimal Charge Control of Plug-In Hybrid Electric Vehicles in Deregulated Electricity Markets. IEEE Trans. Power Syst.26 (3), 1021ā1029. 10.1109/tpwrs.2010.2086083
26
SaberA. Y.VenayagamoorthyG. K. (2011). Plug-in Vehicles and Renewable Energy Sources for Cost and Emission Reductions. IEEE Trans. Ind. Electron.58 (4), 1229ā1238. 10.1109/tie.2010.2047828
27
Shabanpour-HaghighiA.SeifiA. R. (2015). Energy Flow Optimization in Multicarrier Systems. IEEE Trans. Ind. Inf.11 (5), 1067ā1077. 10.1109/tii.2015.2462316
28
SheffiY. (1985). Urban Transportation Networks: Equilibrium Analysis with Mathematical Programming Methods. Englewood Cliffs, NJ, USA: Prentice-Hall.
29
SoroudiA. (2017). Power System Optimization Modeling in GAMS.
30
StrasserT.AndrenF.KathanJ.CecatiC.BuccellaC.SianoP.et al (2015). A Review of Architectures and Concepts for Intelligence in Future Electric Energy Systems. IEEE Trans. Ind. Electron.62 (4), 2424ā2438. 10.1109/tie.2014.2361486
31
TaylorJ. A.HoverF. S. (2012). Convex Models of Distribution System Reconfiguration. IEEE Trans. Power Syst.27 (3), 1407ā1413. 10.1109/tpwrs.2012.2184307
32
TengW.WangY.SunS. (2021). Robust Stability Control for Electric Vehicles Connected to DC Distribution Systems. Front. Energ. Res.9. 10.3389/fenrg.2021.740698
33
UddinM.RomlieM. F.AbdullahM. F.Abd HalimS.Abu BakarA. H.Chia KwangT. (2018). A Review on Peak Load Shaving Strategies. Renew. Sustain. Energ. Rev.82, 3323ā3332. 10.1016/j.rser.2017.10.056
34
WeiW.WuL.WangJ.MeiS. (2018). Network Equilibrium of Coupled Transportation and Power Distribution Systems. IEEE Trans. Smart Grid9 (6), 6764ā6779. 10.1109/tsg.2017.2723016
35
XuD.WuQ.ZhouB.LiC.BaiL.HuangS. (2020). Distributed Multi-Energy Operation of Coupled Electricity, Heating, and Natural Gas Networks. IEEE Trans. Sustain. Energ.11 (4), 2457ā2469. 10.1109/tste.2019.2961432
36
YangQ.FangX. (2017). Demand Response under Realātime Pricing for Domestic Households with Renewable DGs and Storage. IET Generation, Transm. Distribution11 (8), 1910ā1918. 10.1049/iet-gtd.2016.1066
37
YeJ.ZhaoD.ZhangL. (2021). Research on Combined Electricity and Heating System Scheduling Method Considering Multi-Source Ring Heating Network. Front. Energ. Res.9. 10.3389/fenrg.2021.800906
38
ZhaoG. Y.LiuZ. Y.HeY.CaoH. J.GuoY. B. (2017). Energy Consumption in Machining: Classification, Prediction, and Reduction Strategy. Energy133, 142ā157. 10.1016/j.energy.2017.05.110
39
ZhaoJ.MaY.LiuQ.WenL.JiaC.FangY. (2020). A Multi-Source Coordinated Optimal Operation Model Considering the Risk of Nuclear Power Peak Shaving and Wind Power Consumption. IEEE Access8, 189702ā189719. 10.1109/access.2020.3027705
Summary
Keywords
an urban energy system, peak shaving, traffic flows, multi-energy complement, low-carbon operation
Citation
Peng Y, Ma J, Lu C, Chen S and Pei L (2022) Day-Ahead Operation of an Urban Energy System Considering Traffic Flows and Peak Shaving. Front. Energy Res. 10:883598. doi: 10.3389/fenrg.2022.883598
Received
25 February 2022
Accepted
28 March 2022
Published
03 May 2022
Volume
10 - 2022
Edited by
Qinran Hu, Southeast University, China
Reviewed by
Hossam A. Gabbar, Ontario Tech University, Canada
Chuan He, Sichuan University, China
Shuai Lu, Southeast University, China
Updates
Copyright
Ā© 2022 Peng, Ma, Lu, Chen and Pei.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sheng Chen, chenshenghhu@163.com
This article was submitted to Smart Grids, a section of the journal Frontiers in Energy Research
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.