Charge Pricing Optimization Model for Private Charging Piles in Beijing

This paper develops a charge pricing model for private charging piles (PCPs) by considering the environmental and economic effects of private electric vehicle (PEV) charging energy sources and the impact of PCP charging load on the total load. This model simulates users’ responses to different combinations of peak-valley prices based on the charging power of PCPs and user charging transfer rate. According to the regional power structure, it calculates the real-time coal consumption, carbon dioxide emissions reduction, and power generation costs of PEVs on the power generation side. The empirical results demonstrate that the proposed peak-valley time-of-use charging price can not only minimize the peak-valley difference of the total load but also improve the environmental effects of PEVs and the economic income of the power system. The sensitivity analysis shows that the load-shifting effect of PCPs will be more obvious when magnifying the number of PEVs by using the proposed charging price. The case study indicates that the proposed peak, average, and valley price in Beijing should be 1.8, 1, and 0.4 yuan/kWh, which can promote the large-scale adoption of PEVs.


Introduction
In order to alleviate the problem of environmental pollution and the greenhouse gas effect, many countries, including China, have formulated relevant policies to promote the development of electric vehicles (EVs). However, the uncontrolled charging behavior of increasing numbers of EVs will not only affect the stability of the power grid, but also the coal consumption and carbon emissions of EV charging energy sources. The charge pricing of EVs is the key to guide the orderly charging of EV users and, importantly, to ensure the benefits of charging facility operators. The research on charge pricing of EVs can be divided into three aspects.
From the perspective of stakeholders' economic benefits: Finn et al. [1] determined the charging price of EVs to increase the economic benefits of participants; Anderson [2] proposed a two-stage charge pricing approach to balance the interests between public institutions and consumers; Lv [3] designed an optimal time-of-use (TOU) tariff scheme to balance the interests of the users and the power grid by using a multi-objective genetic algorithm; to minimize the users' charging costs, Zhang et al. [4] set a charging price to standardize their charging behavior, Yu et al. [5] introduced an incentive factor into the charging/discharging power limitation and charging/discharge price, and Zhang and Li [6] proposed a regime-switching based risk management scheme against the electricity peak prices.
From the perspective of smoothing the power grid's load: Ge et al. [7] ordered the charging of EVs to achieve the peak-filling effect of EVs by optimizing the valley price; Xiang et al. [8] established the vehicle-to-grid charging and discharging period model and the peak-valley electricity price model of the price limitation, generator sets' installed capacity, PEV energy sources' coal consumption and emissions reduction, and the power system's electricity-sale incomes, the peak-valley TOU charging price is proposed to minimize the peak-valley difference of the total load.

Charge Pricing Factors of PEVs
In the process of PEV charge pricing, this paper considers factors including the charging probability of PEVs, price response of PEV users, charging load forecast of PCPs, energy saving and emissions reduction of PEVs, and economic benefits of the power system. The charging load of PCPs is predicted according to the PEV charging probability and users' price response. Considering the line loss and charge loss, the sum of PEV load and conventional load equals the total load. Based on the total load and the local power structure, the environmental and economic effects of PCPs are calculated.

Charging Probability of PEVs
This paper divides a day into 24 h (t = 0, 1, 2, . . . , 23). PEV users are assumed to immediately charge when they return home. According to the total number of PEVs (N) and the travel final return probability of PEVs (f (t)), the number of PEVs beginning to charge (N r (t)) can be expressed as: where [] is the rounding symbol.

Price Response of PEV Users
The equations of user price response are provided by Chang et al. [22]: the dead zone thresholds (yuan/kWh) of peak-valley, peak-average, and valley-average price differences mean the peak-valley, peak-average, and valley-average price differences which users begin to respond to, and their values are 0.2, 0.2, and 0.1 respectively; the saturated zone thresholds (yuan/kWh) of peak-valley, peak-average, and valley-average price differences represent the peak-valley, peak-average, and valley-average price differences which users no longer response to, and their values are 1.4, 0.8, and 0.6 respectively; 0.25, 1, and 0.8 are price response curve linear region slopes in peak-valley, peak-average, and valley-average periods; 30%, 60%, and 40% are the saturation values of user transfer percentages in peak-valley, peak-average, and valley-average periods. Furthermore, this paper introduces the existing charging price (P 0 , yuan/kWh) into Equations (2)-(4) by considering its influence on user price response. Thus, peak-valley (a pv ), peak-average (a pa ), and valley-average (a av ) user response transfer percentage can be adjusted as where P p , P a , and P v are the peak, average, and valley charging prices, yuan/kWh. After users respond to the TOU tariff, the number of PEVs beginning to charge at each moment (N c (t)) should be an non-negative integer. Thus, the equation of N c (t) in the research of Chang et al. [22] is adjusted as: where N p and N a are the sum of the number of PEVs beginning to charge at peak and average period respectively; Values of v(t) can be 1, 2, and 3 which represents the peak, average, and valley period, respectively.

Charging Load Forecast of PCPs
This paper assumes that the daily travel mileage of private vehicles are the same, which is determined by the annual travel mileage of private vehicles. The charging time (T c , hour) in the research of Su and Liang [23] is determined as: where p is the charging power of a single PCP, MW; M e is the 100-kilometer electric consumption of fuel vehicles, kWh; L is the annual travel mileage of private vehicles, km. According to Equation (7), the charging time is 2 h. Thus, the charging power of PCPs (p c (t), MW) at each moment is expressed as: The load of PCPs (F e (t), MW) with the line loss (ω) at each moment is calculated as: Considering the charging loss (η) of PCPs, the purchased electric quantity of EV users (Q(t), MWh) is: The total load (F(t), MW) can be expressed as where F t (t) is the conventional load, MW.

Energy-Saving and Emission-Reduction Effects of PCPs
Considering the intermittency of wind power and hydroelectric power, this paper introduces the utilization rate of wind power (α(t)) in 24 h and the generating guarantee rate of hydroelectric power (β) to determine their available generation capacities (G w (t) and G h (t), MW) in each moment. where G w and G h are the installed capacities of wind power and hydroelectric power, MW. Coal power generating units are assumed to be peaking units. Coal consumption per electricity unit (s(t), t/MWh) is determined as: where G n is the installed capacity of nuclear power, MW; c is the coal consumption of coal power per electricity unit, t/MWh. Coal consumption of PEVs (S(t), t) is calculated as: Carbon emissions per electricity unit (e e (t), t/MWh) is defined as: where e w , e h , e n , and e c are carbon emissions of wind, hydroelectric, nuclear, and coal power respectively, t/MWh. Carbon emissions of PEVs (E e (t), t) is calculated as Carbon emissions of fuel vehicles (E g (t), t) is determined as where W g is 100kilometer gasoline consumption of fuel vehicles, L; e g is carbon emissions of gasoline, t/L.

Economic Benefits of the Power System
The power plants generate the electricity from the wind, hydroelectric, nuclear, and thermal energy, and then the electricity is transferred by the power grid to PCPs for PEVs' charging. Thus, this paper takes power plants and the power grid as an integrated system to calculate the economic benefits that the power system supplies energy to PEVs. The equation of generation costs per electricity unit (C(t), yuan/MWh) is expressed as where C w , C h , C n , and C c are generation costs of wind, hydroelectric, nuclear, and coal power respectively, yuan/MWh. Net incomes that the power system sells to PEV users (I(t), yuan) is calculated as: Sustainability 2017, 9, 2075 6 of 15

Charge Pricing Model of PCPs
Based on the work of Section 2, the charge pricing model aims at proposing an optimized peak-valley TOU charging price which could minimize the peak-valley difference of the total load and achieve better environmental and economic effects of PCPs.

Constraints
(1) Price constraint The order in magnitude of charging prices should be valley price, average price and peak price, and the charging prices should be in the range limit of charging price.
where the upper limit of charging price is P up , yuan/kWh; the lower limit of charging price is P low , yuan/kWh.
(2) Energy-saving effect constraint Coal consumption of PEVs (S(t)) should not be higher than the initial one (S 0 (t)). (

3) Emission-reduction effect constraint
Carbon emissions reduction of PEVs equals carbon emissions of fuel vehicles minus carbon emissions of PEVs, which should not be less than the initial one.
(4) Economic benefit constraint Net incomes of the power system cannot be lower than initial ones.
The peak-valley difference of the total load should not be higher than the initial one.

Optimization Goals
Minimize the peak-valley difference of the total load: Figure 1 shows the design of the peak-valley TOU charge pricing optimization model for PCPs. The coal consumption of PEVs, carbon emissions reduction of PEVs, and net incomes of the power system are determined by the charging regularity of PEVs, the price response of users and the local power structure. Under the constraints of price range limit, generation set installed capacity, energy-saving and emission-reduction effects, economic benefit, and peak-valley load difference constraints, the peak-valley TOU charging price is proposed to minimize the peak-valley difference of the total load. The number of PEVs begining to charge at each moment (Nc(t)) after users respond to the TOU tariff is generated by Nr(t), peak, average and valley price (Pp, Pa and Pv) and user price response formulas (Eqs. (2)(3)(4)). (Eqs. (5 and 6))

Flow Chart of the Optimized Peak-Valley TOU Charging Price
The number of PEVs beginning to charge at each moment (Nr(t)) is generated by the total number of PEVs (N) and the travel final return probability of PEVs (f(t)). (Eq. (1)) The charging power of PCPs at each moment (pc(t)) is determined by Nc(t), p and Tc. (Eq. (8)) The total load (F(t)) is generated by Fe(t) and the typical load (Ft). (Eq. (11)) Coal consumption per electricity unit (s(t)) is calculated by The load of PCPs (Fe(t)) is calculated by pc(t) and the line loss (ω). (Eq. (9)) The purchased electric quantity of EV users (Q(t)) is calculated by pc(t) and the charge loss (η).
(Eq. (10)) Net incomes of the power system (I(t)) is generated by

Empirical Results
This paper takes Beijing as an empirical case, and uses the energy consumption and carbon emissions coefficients of BAIC EV200 and Beijing Hyundai Elantra.

Initial Data of the Model
According to the investigated data of United States Transportation Department on household vehicles of United States, Su and Liang [23] approximated PEV travel final return time to the normal distribution by using the maximum likelihood estimation method. This paper transfers the probability density distribution function of PEV travel final return time [23] into a segmentation function by dividing a day into 24 h periods, and uses the rectangle method to integrate the segmentation function with MATLAB. The calculated travel final return probability of PEVs (f(t)) of each moment is shown as Table 1. Beijing plans to construct 360,000 PCPs and basically achieve equipping each PEV with a PCP by 2020 [21]. Thus, the total number of PEVs should be 360,000. Currently, the charging price of PCP (P 0 ) is 0.4733 yuan/kWh in Beijing [24]. The charging power of a single PCP (p) is 0.0025 MW [23]. The 100-kilometer electric consumption of BAIC EV200 (M) is 15 kWh [22]. The annual travel mileage of private vehicles (L) is 12,566 km [25]. The line loss (ω) and the charge loss (η) are 6.68% and 17%, respectively [26]. According to the research of Zhang et al. [26], the conventional load and the power structure of Beijing are shown as Tables 2 and 3. It can be seen that the peak load concentrates on the period from 10:00 to 18:00 and the valley one distributed from 23:00 to 7:00 in Table 2. Table 3 shows that coal power is the main energy source to supply the charging of PEVs. Unlike other generation technology, coal power generates more carbon emissions. This means that more electricity from the coal power will cause PEVs produce more carbon emissions. Table 4 shows the utilization rate of wind power in 24 h which is provided by Song [27]. According to Table 3, hydroelectric-power installed capacity as a proportion of all power installed capacity is 21.05%. Thus, the hydroelectric power's generating guarantee rate should be 80-90% [28]. This paper uses the middle value, 85%, to estimate the utilization rate of hydroelectric power. Unlike the research of Zhang et al., this paper introduces the utilization rate of wind power and hydroelectric power, and assumes to use coal power to modulate the peak load of the power grid. This is to say that the generating capacity of coal power can be more than 10,000 MW when the other power is not enough to fit the electricity demand in the power grid.  According to the energy-saving power generation dispatching method, the order of the power generation should be wind power, hydroelectric power, nuclear power, and coal power [31]. The coal consumption of coal power per electricity unit (c) is 4.47 t/MWh [32]. Carbon emissions of gasoline (e g ) are 0.002135 t/L [33]. According to the research of Chang et al. [22], the 100kilometer gasoline consumption of Elantra (s g ) is 9L, the upper limit of charging price (P up ) is 3.576 yuan/kWh and the lower limit of charging price (P low ) is 0.396 yuan/kWh. Due to the existing charging price being 0.4733 yuan/kWh and the acceptable highest charging price for users being 3.576 yuan/kWh, there is a lot consumer surplus. It is feasible to adopt the peak-valley TOU charging price to guide PEV users' charging behavior and make PCP operators more profitable.

Optimization Results
In order to effectively achieve the load-shifting effect of PEVs, 24 h are divided into three periods-peak period, average period, and valley period-according to the distribution time of Beijing grid peak, average and valley load in Table 2. Based on the charge pricing model in Section 3, the optimized peak-valley TOU charging price is proposed to minimize the peak-valley difference of the total load in Table 5. It can be seen that the peak and average prices are more than double and quadruple of the valley one, and the valley one is lower than the existing charging price 0.4733 yuan/kWh. Comparing with the initial PCP load, the black part indicates the increased load of optimized PCP load and the white part indicates the decreased one in Figure 2. It can be seen that the charging load of PCPs transfers from peak and average period to the valley one. Due to the low valley charging price in Table 5, the charging loads of PCPs at 8:00 and 9:00 in average period decrease to 0 MW and 1.78 MW, and the ones at 10:00, 11:00 and 12:00 in peak period reduced to 1.78 MW, 0 MW and 0 MW, respectively. The price difference between peak period and valley period is much more than that between average period and valley period. Thus, the charging load of PCPs from 10:00 to 18:00 decreases much more than that of 7:00-10:00 and 18:00-23:00. Apart from dead zone thresholds and saturated zone thresholds, the determination of user price response is based on price response curve linear region slopes, saturation values of user transfer percentages, and the comparison of the existing charging price in Equations (2)-(4). It means that only the large price difference can not cause a great number of users to change charging behavior. Thus, the decreased charging load of PCPs between 17:00 and 21:00 is still high. Sustainability 2017, 9,2075 10 of 15 the optimized peak-valley TOU charging price is proposed to minimize the peak-valley difference of the total load in Table 5. It can be seen that the peak and average prices are more than double and quadruple of the valley one, and the valley one is lower than the existing charging price 0.4733 yuan/kWh. Comparing with the initial PCP load, the black part indicates the increased load of optimized PCP load and the white part indicates the decreased one in Figure 2. It can be seen that the charging load of PCPs transfers from peak and average period to the valley one. Due to the low valley charging price in Table 5, the charging loads of PCPs at 8:00 and 9:00 in average period decrease to 0 MW and 1.78 MW, and the ones at 10:00, 11:00 and 12:00 in peak period reduced to 1.78 MW, 0 MW and 0 MW, respectively. The price difference between peak period and valley period is much more than that between average period and valley period. Thus, the charging load of PCPs from 10:00 to 18:00 decreases much more than that of 7:00-10:00 and 18:00-23:00. Apart from dead zone thresholds and saturated zone thresholds, the determination of user price response is based on price response curve linear region slopes, saturation values of user transfer percentages, and the comparison of the existing charging price in Equations 2-4. It means that only the large price difference can not cause a great number of users to change charging behavior. Thus, the decreased charging load of PCPs between 17:00 and 21:00 is still high.  Figure 3 shows that the coal consumption per electricity unit of the valley period is extremely low, and the one between 0:00 and 6:00 is much less than 1.5 t/MWh. While the coal consumption per electricity unit of the peak and average period is particularly high, especially the one between 9:00 and 20:00 which is higher than 2.5 t/MWh. Although the load of PCPs increases during the valley period and decreases during the peak and average period in Figure 2, the increased area of PEV coal consumption during the valley period are much less than the decreased one during the peak and average period in Figure 3. It results in decreasing coal consumption of PEVs in Figure 5.  Figure 3 shows that the coal consumption per electricity unit of the valley period is extremely low, and the one between 0:00 and 6:00 is much less than 1.5 t/MWh. While the coal consumption per electricity unit of the peak and average period is particularly high, especially the one between 9:00 and 20:00 which is higher than 2.5 t/MWh. Although the load of PCPs increases during the valley period and decreases during the peak and average period in Figure 2, the increased area of PEV coal consumption during the valley period are much less than the decreased one during the peak and average period in Figure 3. It results in decreasing coal consumption of PEVs in Figure 5.  Figure 4 shows that carbon dioxide emissions per electricity during the valley period are lower than 0.35 t/MWh, while the ones during the peak period are much higher than 0.55 t/MWh. Like coal consumption of EVs in Figure 3, the increased carbon emissions of PEVs during the valley period are much lower than the decreased ones during the peak and average period. Due to TOU stimulation, some users transfer their charging from the average and peak period to the valley one in Figure 2. Moreover, the low valley price has encouraged the more charging of users during the valley period, and make the optimized charging power of PCPs more than the initial one. According to Equations (10) and (18), it can be known that carbon emissions of fuels vehicles will rise in Figure 4. Thus, the decreased carbon emissions of PEVs and the increased carbon emissions of fuels vehicles result in the increased carbon emissions reduction in Figure 5.  Figure 5 shows the initial and optimized values of the total load's peak-valley difference, PEVs' coal consumption, PEVs' carbon emissions reduction, and the power system's net incomes. The optimized peak-valley difference of the total load reaches a minimum value of 8,741.15 MW, which decreases 2% of the initial total load's peak-valley difference. The coal consumption of PEVs at the power generation side reduces 20%. Carbon emissions reduction at the power generation side increases 28%. Net incomes of the power system that sell to PEV users enhances more than two times. Combining with Figures 3, 4, and 5, it can be known that minimizing the peak-valley difference of   Figure 4 shows that carbon dioxide emissions per electricity during the valley period are lower than 0.35 t/MWh, while the ones during the peak period are much higher than 0.55 t/MWh. Like coal consumption of EVs in Figure 3, the increased carbon emissions of PEVs during the valley period are much lower than the decreased ones during the peak and average period. Due to TOU stimulation, some users transfer their charging from the average and peak period to the valley one in Figure 2. Moreover, the low valley price has encouraged the more charging of users during the valley period, and make the optimized charging power of PCPs more than the initial one. According to Equations (10) and (18), it can be known that carbon emissions of fuels vehicles will rise in Figure 4. Thus, the decreased carbon emissions of PEVs and the increased carbon emissions of fuels vehicles result in the increased carbon emissions reduction in Figure 5.   Figure 4 shows that carbon dioxide emissions per electricity during the valley period are lower than 0.35 t/MWh, while the ones during the peak period are much higher than 0.55 t/MWh. Like coal consumption of EVs in Figure 3, the increased carbon emissions of PEVs during the valley period are much lower than the decreased ones during the peak and average period. Due to TOU stimulation, some users transfer their charging from the average and peak period to the valley one in Figure 2. Moreover, the low valley price has encouraged the more charging of users during the valley period, and make the optimized charging power of PCPs more than the initial one. According to Equations (10) and (18), it can be known that carbon emissions of fuels vehicles will rise in Figure 4. Thus, the decreased carbon emissions of PEVs and the increased carbon emissions of fuels vehicles result in the increased carbon emissions reduction in Figure 5.  Figure 5 shows the initial and optimized values of the total load's peak-valley difference, PEVs' coal consumption, PEVs' carbon emissions reduction, and the power system's net incomes. The optimized peak-valley difference of the total load reaches a minimum value of 8,741.15 MW, which decreases 2% of the initial total load's peak-valley difference. The coal consumption of PEVs at the power generation side reduces 20%. Carbon emissions reduction at the power generation side increases 28%. Net incomes of the power system that sell to PEV users enhances more than two times. Combining with Figures 3, 4, and 5, it can be known that minimizing the peak-valley difference of   Figure 5 shows the initial and optimized values of the total load's peak-valley difference, PEVs' coal consumption, PEVs' carbon emissions reduction, and the power system's net incomes. The optimized peak-valley difference of the total load reaches a minimum value of 8741.15 MW, which decreases 2% of the initial total load's peak-valley difference. The coal consumption of PEVs at the power generation side reduces 20%. Carbon emissions reduction at the power generation side increases 28%. Net incomes of the power system that sell to PEV users enhances more than two times. Combining with Figures 3, 4, and 5, it can be known that minimizing the peak-valley difference of the total load is beneficial to reduce the coal consumption and carbon emissions of EV charging energy sources. Sustainability 2017, 9,2075 12 of 15 the total load is beneficial to reduce the coal consumption and carbon emissions of EV charging energy sources.

Sensitivity Analysis
The total number of private vehicles is 4,523,000 [18] in Beijing, which is 12 times more than that of PEVs. Thus, we magnify the number of PEVs from 1 to 12 times to analyze its impact on the peakvalley difference of the total load. The peak-valley difference reduction rate means the decreased peak-valley difference as a proportion of initial peak-valley difference, which can be expressed as: In Figure 6, the value of peak-valley difference reduction rate rises with the magnifying number of PEVs, which means that the optimized peak-valley TOU charging price can achieve the loadshifting effect of PCPs. In addition, the load-shifting effect is enhanced with the increasing number of PEVs. When all the private vehicles in Beijing are replaced by PEVs, the peak-valley difference reduction rate will reach 20%.  Figure 7 shows the total load with magnifying the number of PEVs 10-fold. At that time, PEVs' number will account for 80% of private vehicles. The initial scenario (S0) means the scenario that the initial number of PEVs which is 360,000 [1]. The magnifying number of PEVs will rise the peak load and enlarge the peak-valley difference of the total load. It can be seen that the load-shifting effect of PCPs in S0 is not obvious. However, the proposed charging price

Sensitivity Analysis
The total number of private vehicles is 4,523,000 [18] in Beijing, which is 12 times more than that of PEVs. Thus, we magnify the number of PEVs from 1 to 12 times to analyze its impact on the peak-valley difference of the total load. The peak-valley difference reduction rate means the decreased peak-valley difference as a proportion of initial peak-valley difference, which can be expressed as: r = ((MAX(F 0 (t)) − MI N(F 0 (t))) − (MAX(F(t)) − MI N(F(t))))/(MAX(F 0 (t)) − MI N(F 0 (t))) (27) In Figure 6, the value of peak-valley difference reduction rate rises with the magnifying number of PEVs, which means that the optimized peak-valley TOU charging price can achieve the load-shifting effect of PCPs. In addition, the load-shifting effect is enhanced with the increasing number of PEVs. When all the private vehicles in Beijing are replaced by PEVs, the peak-valley difference reduction rate will reach 20%. Sustainability 2017, 9,2075 12 of 15 the total load is beneficial to reduce the coal consumption and carbon emissions of EV charging energy sources.

Sensitivity Analysis
The total number of private vehicles is 4,523,000 [18] in Beijing, which is 12 times more than that of PEVs. Thus, we magnify the number of PEVs from 1 to 12 times to analyze its impact on the peakvalley difference of the total load. The peak-valley difference reduction rate means the decreased peak-valley difference as a proportion of initial peak-valley difference, which can be expressed as: In Figure 6, the value of peak-valley difference reduction rate rises with the magnifying number of PEVs, which means that the optimized peak-valley TOU charging price can achieve the loadshifting effect of PCPs. In addition, the load-shifting effect is enhanced with the increasing number of PEVs. When all the private vehicles in Beijing are replaced by PEVs, the peak-valley difference reduction rate will reach 20%.  Figure 7 shows the total load with magnifying the number of PEVs 10-fold. At that time, PEVs' number will account for 80% of private vehicles. The initial scenario (S0) means the scenario that the initial number of PEVs which is 360,000 [1]. The magnifying number of PEVs will rise the peak load and enlarge the peak-valley difference of the total load. It can be seen that the load-shifting effect of PCPs in S0 is not obvious. However, the proposed charging price 10 Figure 7 shows the total load with magnifying the number of PEVs 10-fold. At that time, PEVs' number will account for 80% of private vehicles. The initial scenario (S0) means the scenario that the initial number of PEVs which is 360,000 [1]. The magnifying number of PEVs will rise the peak load and enlarge the peak-valley difference of the total load. It can be seen that the load-shifting effect of PCPs in S0 is not obvious. However, the proposed charging price displays a better load-shifting effect with the magnifying number of PEVs, which can decrease 16.6% of the total load's peak-valley difference. It means that the proposed peak-valley TOU charging price can effectively achieve the load-shifting effect of PCPs with the magnifying number of PEVs.
Sustainability 2017, 9,2075 13 of 15 displays a better load-shifting effect with the magnifying number of PEVs, which can decrease 16.6% of the total load's peak-valley difference. It means that the proposed peak-valley TOU charging price can effectively achieve the load-shifting effect of PCPs with the magnifying number of PEVs.

Conclusions
To minimize the peak-valley difference of the total load, this paper proposes a peak-valley TOU charge pricing model for PCPs with five constraints including charging price range limit, energysaving and emission-reduction effects, economic benefit and peak-valley load difference constraints. This model analyzes factors that influence the charge pricing and simulates the charging load of PEVs and price response of users to determine the coal consumption and carbon emissions of PEVs and net incomes of the power system. The charging load of PCPs is predicted according to the charging probability of PCPs and PCP constructing number of Beijing in 2020. The analysis results show that the optimized peak, average and valley charging price in Beijing should be 1.8 yuan/kWh from 10:00 to 17:00, 1 yuan/kWh and 0.4 yuan/kWh from 23:00 to 6:00. The proposed peak-valley TOU charging price can minimize the peak-valley difference of the total load. Moreover, the load-shifting effect of PCPs will be improved with magnifying the number of PEVs. When the proportion of PEVs in total private vehicles numbers reaches 80%, the proposed charging price can decrease 16.6% of the total load's peak-valley difference. Therefore, the proposed charging price is beneficial for the government to promote the large-scale development of PEVs. The proposed charging price can reduce 20% coal consumption and increase 28% carbon emissions reduction of PEVs at the power generation side. Additionally, it enhances more than 2 times net incomes from PEV electricity sales of the power system.
Author Contributions: In this paper, Yakun Zhang sorted out the studies on charge pricing in Section 1; Yanni Liang analyzed the factors that influence the charge pricing, designed the charge pricing model for PCPs, and wrote Sections 1, 2, and 3; Xingping Zhang analyzed the results of the optimization model and sensitivity analysis, drew the conclusions for this paper, and wrote Sections 4 and 5; Yinhe Bo and Hongyang Zhang collected the related data and basic parameters for the charge pricing model in Section 4.

Conclusions
To minimize the peak-valley difference of the total load, this paper proposes a peak-valley TOU charge pricing model for PCPs with five constraints including charging price range limit, energy-saving and emission-reduction effects, economic benefit and peak-valley load difference constraints. This model analyzes factors that influence the charge pricing and simulates the charging load of PEVs and price response of users to determine the coal consumption and carbon emissions of PEVs and net incomes of the power system. The charging load of PCPs is predicted according to the charging probability of PCPs and PCP constructing number of Beijing in 2020. The analysis results show that the optimized peak, average and valley charging price in Beijing should be 1.8 yuan/kWh from 10:00 to 17:00, 1 yuan/kWh and 0.4 yuan/kWh from 23:00 to 6:00. The proposed peak-valley TOU charging price can minimize the peak-valley difference of the total load. Moreover, the load-shifting effect of PCPs will be improved with magnifying the number of PEVs. When the proportion of PEVs in total private vehicles numbers reaches 80%, the proposed charging price can decrease 16.6% of the total load's peak-valley difference. Therefore, the proposed charging price is beneficial for the government to promote the large-scale development of PEVs. The proposed charging price can reduce 20% coal consumption and increase 28% carbon emissions reduction of PEVs at the power generation side. Additionally, it enhances more than 2 times net incomes from PEV electricity sales of the power system.