Next Article in Journal
Towards Place-Based Research to Support Social–Ecological Stewardship
Next Article in Special Issue
Methodological Issues Regarding Biofuels and Carbon Uptake
Previous Article in Journal
Determinants of Agricultural Diversification in a Hotspot Area: Evidence from Colonist and Indigenous Communities in the Sumaco Biosphere Reserve, Ecuadorian Amazon
Previous Article in Special Issue
A Numerical Study of the Temperature Reduction by Water Spray Systems within Urban Street Canyons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Security-Constrained Unit Commitment Considering Differentiated Regional Air Pollutant Intensity

School of Electrical Engineering and Automation, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Sustainability 2018, 10(5), 1433; https://doi.org/10.3390/su10051433
Submission received: 6 April 2018 / Revised: 28 April 2018 / Accepted: 30 April 2018 / Published: 4 May 2018
(This article belongs to the Special Issue Sustainable Energy Development under Climate Change)

Abstract

:
Conventional environmental-economic power dispatch methods constrain the total amount of emissions of power plants, and they succeed in reducing emissions from the power sector. However, they fail to address the mismatch between emission reductions and the resulting changes in regional air quality. This paper proposes an ecology- and security-constrained unit commitment (Eco-SCUC) model considering the differentiated impacts of generation-associated emissions on regional air quality. A Gaussian puff dispersion model is applied to capture the temporal-spatial transport of air pollutants. Additionally, an air pollutant intensity (API) index is defined for assessing the impacts of emissions on the air quality in regions with differentiated atmospheric environmental capacities. Then the API constraints are formulated based on air quality forecast and included in SCUC model. Moreover, the stochastic optimization is employed to accommodate wind power uncertainty, and the Benders decomposition technique is used to solve the formulated mixed-integer quadratic programming (MIQP) problem. Case studies demonstrate that the Eco-SCUC can cost-effectively improve air quality for densely-populated regions via shifting generation among units and can significantly reduce the person-hours exposed to severe air pollution. Furthermore, the benefits of wind power for air quality control are investigated.

1. Introduction

Air quality issues arouse wide concerns all over the world. Electricity production is responsible for large shares of atmospheric pollutants including sulfur dioxide (SO2), nitrogen oxide (NOx), and fine particulate matter (PM2.5) [1,2,3]. According to [4,5,6], coal-fired power plants are primary stationary sources of air pollutants that adversely harm human health. Population exposure to ambient air pollutant concentration that exceeds the World Health Organization (WHO) guidelines would cause significant mortality and morbidity [7]. During the past decades, air quality control, especially the control of environmental and health damages attributable to coal-fired generation, has received much attention in the power sector [8,9,10,11].
In China, although renewable energy resources are playing significant roles in improving air quality, they cannot fully substitute for traditionally fueled generation at the moment. Coal-fired power plants still account for over 72% of China’s electricity generation in 2015, and the nation is expected to maintain a coal-dominant generation-mix through 2040 [12]. Moreover, renewable generation is regarded as long-term solution due to its uncertain nature, such as wind speed and solar radiation [13,14], and the costly investment [15,16]. Therefore, in the short-run, it is necessary to develop new dispatch model to implement more timely and refined air quality control.
Traditional environmental economic dispatch (ED) or unit commitment (UC) constrains the total emissions of the power sector [17], and various environment-aware technologies and policies are considered [18,19]. Wei et al. [20] proposed a bi-objective ED considering carbon capture plants and wind power uncertainty. Zhang et al. [21] incorporated carbon emission trading in UC problem, which is solved by an improved particle swarm optimization. Abdollahi et al. [22] introduced a cost-emission-based UC taking into account demand response and emission regulation policies. Other critical techniques include: energy storage [23,24,25], electric vehicle [26], and power-to-gas [27,28]. It is noted that most environmental ED/UC models enforce undifferentiated control on emissions or impose uniform incentives to regulate air pollutants.
It is true that the undifferentiated approach can result in a net reduction of air pollutants on the individual generating unit level. However, on the level of regional air quality management, a more refined model is in need because of the trans-regional dispersion characteristics of air pollutants. The impact of generation-associated emissions on regional air quality is a complex problem. It involves meteorology (wind direction, wind speed, atmospheric stability), geography (distances between polluting units and polluted regions), ecology (self-purification capability of the atmosphere), and electricity (system load, network constraints). To address this issue, an air pollutant dispersion model is needed to describe the dynamic transport of air pollutants under variant meteorological conditions [29,30,31]. Sullivan et al. in 1973 [32] incorporated an air pollutant dispersion model into ED to regulate ground-level SO2 concentrations. Schwezer et al. [33] integrated meteorology-aware constraints into a fuel-mix problem. Chu et al. in 1977 [34] proposed an optimal ED implementing marginal environmental taxes under different weather conditions. However, related work has stagnated for decades due to the technology limitations. Nowadays, the air quality modeling and advances in weather forecasting enable us to continue this kind of research [35,36].
On the one hand, the geographical distribution of generating units influences the population exposure to the generation-associated air pollution. For example, the near-city units may severely pollute its downwind human habitats, while the remote units near coal-mines and the wind farms in unmanned regions hardly impact residential air quality. Standing from the people-oriented point of view, these ‘low-impact’ units are encouraged to generate more power within the system limits. On the other hand, regional air quality stress is impacted by local atmospheric environmental capacity (AEC). AEC is for recognizing atmospheric stress limits above which human and ecosystem health will face severe challenges. Generally, AEC is studied in terms of critical load [37,38] or ecological footprint [39,40], and its value is impacted by regional population, atmospheric ventilation, dry/wet deposition, humidity [41,42]. Densely-populated regions have lower AEC and sparsely populated region is with higher AEC [43]. This inspires us to formulate a refined power dispatch to ‘transfer’ the potential air pollution away from densely-populated regions to sparsely populated regions and unmanned regions, thereby reducing the total person-hours exposed to severe air pollution.
Thus, this paper proposes an ecology- and security-constrained unit commitment (Eco-SCUC) model integrating air quality intensity (API) constraints and stochastic wind power. The API index is defined for assessing regional air pollutant saturation level considering the dispersion processes, weather conditions and geographical distribution of system and AEC. When the API ≥ 100%, the region and its population will suffer air pollutant saturated time. The aim of the Eco-SCUC is to balance the regional API distribution and reduce the person-hours exposed to saturated pollution via shifting generation among units. For example, in Figure 1, under the predefined wind direction, the API constraints will help curtail generation of power plant A to prevent air quality deterioration in densely populated region and encourage generation of power plant B and other remote power plants to meet the system load, because their emissions are transported to sparsely-populated and unmanned regions with high AECs. Thus, the health impacts in densely populated region are alleviated and the total person-hours exposed to severe air pollution for the entire region are reduced. As depicted in the figure, the dispatching center is the controller; the coal-fired power plants and wind farm are the controlled objects; the air quality ground stations are the data collectors, and the collected data are sent to the dispatching center to formulate API constraints; the solutions to the Eco-SCUC are the control signals to scheduling generation.
The main contributions and novelty of this paper include:
  • Gaussian puff dispersion is used to model the transport of emissions, and the contribution of emissions to air pollutant concentration increments in targeted regions is calculated. The API is calculated as the ratio of local air pollutant concentration to local AEC-based threshold concentration. The API constraints are formulated and integrated into typical SCUC.
  • The Monte Carlo simulation with Latin hypercube sampling (LHS) and scenario reduction is used to generate wind power volatility scenarios. Additionally, Benders decomposition is used to solve the formulated stochastic mixed-integer quadratic programming (MIQP) problem.
  • Case studies prove the effectiveness of the proposed model in improving air quality for densely populated region and reducing the person-hours exposed to severe air pollution for the entire region via shifting generation among units.
  • Compared to undifferentiated dispatches, the Eco-SCUC considers trans-regional pollutant dispersion effects and differentiated population exposure to harmful pollutants. Instead of enforcing uniform regulations on emission, API constraints are formulated according to the differentiated impacts of coal-fired units on regional air quality. The model succeeds in mitigating the air pollution in densely-populated regions via shifting generation away from high-impact units to low-impact units and wind power.
The environmental benefits that the Eco-SCUC model can deliver are as follows: firstly, it can balance air pollution between densely-populated regions and sparsely-populated regions, thereby alleviating urban air pollution pressure; secondly, as the model schedules generation with respect to the differentiated impacts of units on regional air quality, it can reduce the average and peak concentrations of air pollutant in populated regions and help improve the level of air quality assessment of the city; thirdly, it significantly reduces the person-hours exposed to harmful air pollutants, such as SO2, and abates the associated human health impacts. Note that the proposed model does not include chemical conversion from primary pollutants to secondary pollutants for now, and integrating secondary pollutants into the Eco-SCUC model requires greater research effort.
The rest of this paper is organized as follows: Section 2 formulates the differentiated API constraints; Section 3 describes the wind power generation scenarios; Section 4 formulates the Eco-SCUC and provides the solution methodology; and Section 5 presents the case studies while the conclusion is drawn in Section 6.

2. Formulating Differentiated API Constraints

In order to formulate regional differentiated API constraints, two critical parameters need to be ascertained using applicable methods: generation-associated concentration and regional threshold concentration. The former generally refers to the contribution of emissions of power plants to the pollutant concentration at predefined location. Therefore, it is necessary to include a deterministic and accurate dispersion model which can accommodate the varying meteorological conditions. In this context, a Gaussian puff dispersion model is applied to capture the dynamic transport of air pollutants from the emitting units to the concentration monitoring point. Secondly, the threshold concentration at the predefined location reflects the air pollutant carrying capacity in the concerned region, which depends on its AEC, population, area, etc. In this paper, an AEC method is employed to determine the threshold concentration of air pollutants with respect to geographical variations. The methodology of API constraints formulation is shown in Figure 2, and the modeling of the Gaussian puff model and AEC model is presented in the following sub-sections accordingly.

2.1. Gaussian Puff Dispersion

There are two commonly used air dispersion models: the steady-state Gaussian plume model and the integrated Gaussian puff model. The former is based on the assumption of constant wind speed and unchanging wind directions [34]. In comparison, the latter can accommodate fluctuating wind directions, wind speeds, and atmospheric stabilities [30]. Thus, the Gaussian puff model is used to describe the link between emissions at stacks (source) and air pollutant concentrations in polluted regions (receptor). The diagram of the Gaussian puff model is shown in Figure 3.
The Gaussian puff dispersion function is expressed as follows [34]:
G i ( τ , t ; x , y , z ) = 1 ( 2 π ) 3 / 2 σ x ( τ , t ) σ y ( τ , t ) σ z ( τ , t ) e Ψ / 2
The intermediary parameter Ψ is calculated using Equations (A1) and (A2). When the atmospheric condition is steady, standard deviation σx/y/z(·) is in log-linear relation with respect to dispersion period (tτ), which is expressed in (A3); when atmospheric condition is unsteady, the relationship is expressed in Equation (A4). The atmospheric stability class ϖ can be estimated by Pasquill grading method based on cloud coverage, solar radiation, wind speed and actual measurements [44]. The solar radiation level (−2 to +3) with respect to cloud coverage and solar altitude is given in Table A1; stability class ϖ (A–F) with respect to wind speed and solar radiation is given in Table A2; relation between standard deviation and dispersion period with respect to ϖ is shown in Figure A1.
Based on Equation (1), the contribution of emissions to regional air pollutant concentration Ci is:
C i ( t ; j ) = τ t τ Q i ( τ ) G i ( τ , t ; j )
Q i ( t ) = s i F i ( p i t )
F i ( p i t ) = ( a i p i t 2 + b i p i t + c i ) Δ T
where j represents the location (xj,yj) in the targeted polluted region, in this paper we only care the ground-level concentration (assuming z ≈ 0).

2.2. Atmoshperic Environmental Capacity (AEC)

AEC represents the atmospheric stress limit within which the environmental system can return to its original condition [43]. Widely applied approaches to estimating AEC are: (1) A-value method; and (2) multisource modeling [45]. The latter requires accurate air quality model, such as the community multi-scale air quality (CMAQ) modeling. Its disadvantage is the neglect of accumulation-removal process of pollutants in the entire region. Thus, A-value method is applied in this paper. According to [43], for a given volume in atmosphere and time interval, AEC is defined as the removal rate of pollutant at the threshold concentration where the air pollutant emission rate (source) and removal rate (sink) reach the source-sink equilibrium. In the equilibrium, regional threshold concentration can be ascertained according to local AEC. The equilibrium is expressed in Equations (A5) and (A6), and the derived concentration and A-value is expressed in Equations (A7) and (A8).
When the generation and removal of air pollutant reach the equilibrium at location j and time t ( c ¯ = C T , t j ,   t ), the threshold concentration at location j and time t can be expressed as:
C T , t j = η t j A t j S j
To sum up, the steps of calculating regional threshold concentration is as follows:
Step 1: Determine the removal rate of air pollutant ηtj.
The removal rate of air pollutant can be obtained via official environmental statistic Yearbooks or Bulletins. In this paper, climate-average values of AEC during the air pollution-prone seasons (November-February) are used, referring to [41].
Step 2: Determine A-value Atj and examined area Sj.
Atj can be calculated using Equation (A8), where U, H, vd, and vw are obtained through observation data or modeling data. This paper mainly concerns atmospheric ventilation and considers A t j = π U t j H / 2 , where Utj is hourly predicted wind speed and H = 276.85 m referring to winter-average H value in [41]. The resolution of Sj should be in consistency with that of ηtj and Atj.
Step 3: Calculate the threshold concentration CT,tj.
Based on the obtained ηtj, Atj and Sj, regional CT,tj are calculated by Equation (5).

2.3. Air Pollutant Intensity (API)

API is defined as the ratio of local air pollutant concentration to local threshold concentration:
A P I i ( t ; j ) = C i ( t ; j ) C T , t j × 100 %
A P I 0 , t j = C 0 , t j C T , t j × 100 %
A P I ( t ; j ) = i = 1 N G A P I i ( t ; j ) + A P I 0 , t j
where APIi(t;j) is calculated by Equations (2) and (6); API0,tj is calculated by Equation (7); C0,tj comes from 24-h air quality forecast from ground stations; CT,tj is calculated by Equation (5). According to Equation (8), the API at the predefined location is the sum of all generation-associated APIs and background API.
Figure 4 shows the distribution of threshold concentration of SO2 and API of SO2 of one city in Northeast China, respectively. It illustrates that densely populated regions (low-AEC regions) have lower threshold concentrations and higher APIs.

2.4. Differentiated API Constraints

The differentiated API constraints can be expressed as:
A P I ( t ; j ) max { A P I E , t j ,   A P I 0 , t j }         j
APIE,tj refers to the expected API value at location j and time t, which is normally set to 100% (air pollutant saturated state). Additionally, the value of APIE,tj can be differentiated according to the local pollution level. For example, it can be set to higher value (>100%) when the air quality control is loosened in the targeted region; reversely, it can be set to lower value (<100%) when the region suffers severe air pollution and the air quality control is tightened.
When the background API value API0,tjAPIE,tj, the impacting units are expected to shut down or curtail their generation without violating system security constraints; when API0,tj < APIE,tj, the impacting units are allowed to generate more to meet the load.

3. Formulating Wind Power Scenarios

The volatile and intermittent nature of wind leads to the inaccuracy of wind power prediction, challenging the secure operation of power system. According to [47], the intermittency refers to the unavailability of wind power within a relatively long period (≥1 h) and can be addressed by hourly wind forecast; the volatility refers to small fluctuations and can be mitigated by simulating possible wind power scenarios (<1 h). In this paper, the Weibull distribution is used to simulate the wind speed patterns [15], based on which the Monte Carlo simulation is used with LHS and scenario reduction techniques to generate possible wind power scenarios.

3.1. Distribution of Wind Speed

Field data suggest that wind speed at a given location follows Weibull distribution [48], whose cumulative distribution function (CDF) and probability density function (PDF) are:
F V ( v ) = 0 v f v ( λ ) d λ = 1 exp [ ( v c ) k ]
f V ( v ) = k v k 1 c k exp [ ( v c ) k ] ,         v 0
The mean and standard variance of Weibull function is expressed by Equations (A9)–(A11). When k = 1, the distribution of V follows exponential distribution; when k = 2, the distribution of V follows Rayleigh distribution. To embrace the variety of wind speed patterns, our derivation of wind power distribution is based the general assumption that k > 0. Actual wind measurement data are used to estimate parameters c and k, as in Equations (A12) and (A13), and the estimation values are as in Equations (A14)–(A16).

3.2. Distribution of Wind Power

The relationship between wind power and wind speed depends on factors including generator efficiencies, wind rotor, inverter, gearbox and other wind turbine characteristics [49]. As shown in Figure A2, the simplified piecewise relationship can be expressed using Equation (A17). Two discrete probabilities for PW = 0 and PW = pwr [50] are expressed in Equations (A18) and (A19). According to the theory for the function of random variables [15], the PDF and CDF of wind power can be expressed as:
f P W ( p w ) = k h v in p w r c [ ( 1 + h p w p w r ) v in / c ] k 1 exp { [ ( 1 + h p w p w r ) v in / c ] k }
F P W ( p w ) = Pr ( P W p w ) = { 0 ( p w < 0 ) 1 exp { [ ( 1 + h p w p w r ) v in / c ] k } + exp [ ( v out c ) k ] ( 0 p w < p w r ) 1 ( p w p w r )
h = v r v in 1
Examples of the CDF and PDF of wind power are plotted in Figure 5, respectively.

3.3. Scenario Generation and Reduction

To capture the stochastic availability of wind power, first, the Monte Carlo simulation is used to generate a large number of scenarios based on the wind power distribution and forecasted wind power. Second, LHS [47], which can be easily included in the Monte Carlo method, is applied to reduce the number of runs and achieve more appropriate random distribution of wind power. Simulation tests show that the Monte Carlo method with LHS can better approximate the required PDF of wind power than the simple Monte Carlo. Third, scenario reduction technique [51] is used to remain representing scenarios by eliminating scenarios of low probabilities and aggregating scenarios of close probabilities within an acceptable level of variation of the simulation precision.

4. Modeling Eco-SCUC

As shown in Figure 6, the proposed Eco-SCUC integrates differentiated API constraints into typical SCUC model for a hybrid power system consisting of coal-fired units and wind power generation. To accommodate wind power uncertainty, the day-ahead forecasted wind power and system load profile are required as electrical inputs for the Eco-SCUC. To differentiate the impacts of coal-fired units on regional air quality stress, the environmental inputs including air quality and weather forecasts are needed. The goal of the Eco-SCUC is to balance regional air quality stress among urban, suburban, and unmanned regions via shifting generation among units.
In this paper, SO2 is assumed to be the major air pollutant in the targeted region, thus, it is taken as the single pollutant under study. However, the model can be extended to accommodate other single or multiple pollutants (e.g., NOx, PM2.5). This study does not consider chemical conversions from primary pollutants to secondary pollutants for now.

4.1. Integrating Wind Power Scenarios

SCUC with volatile wind power is formulated as a stochastic optimization [52] problem with the objective function of Equation (15) and constraints (Equations (16)–(36)). To consider the uncertainties of wind power, Monte Carlo simulation with LHS and scenario reduction techniques are used to generate possible scenarios based on the wind power distribution described in Equation (12). In each scenario, the simulated hourly random wind power will replace the forecasted hourly wind power in re-dispatch, thus, the hourly constraints (Equations (16)–(27)) are re-used as the scenario constraints (Equations (28)–(36)).

4.1.1. Objective Function

The objective function of typical SCUC is the sum of fuel, startup and shutdown costs of coal units, while assuming that the operating cost of wind-powered units is zero:
min   t = 1 N T i = 1 N G [ γ i F i ( p i t ) I i t + U i t + D i t ]
where Fi(pit) is calculated using Equation (4); and Uit/Dit is normally regarded as step function with respect to time t—in this paper we take it as a constant value for simplification.

4.1.2. Hourly Constraints

Hourly SCUC constraints include system power balance constraints (Equation (16)), wind power limit (Equation (17)), unit minimum on/off limits (Equations (18) and (19)), unit ramping up/down limits (Equations (20) and (21)), unit generation limits (Equation (22)), system regulation up/down limits (Equations (23) and (24)), system spinning/operating reserve limits (Equations (25) and(26)), and generalized DC power flow constraints (Equation (27)):
i = 1 N G p i t I i t + i = 1 N W p w i t = P D t + P L t       t
0 p w i t p w i t f       i , t
[ X i , t 1 on T i on ] [ I i , t 1 I i t ] 0       i , t
[ X i , t 1 off T i off ] [ I i , t 1 I i t ] 0       i , t
p i t p i , t 1 [ 1 I i t ( 1 I i , t 1 ) ] P i RU + I i t ( 1 I i , t 1 ) P i . min
p i , t 1 p i t [ 1 I i , t 1 ( 1 I i t ) ] P i RD + I i , t 1 ( 1 I i t ) P i . min
P i , min I i t p i t P i , max I i t       i , t
i = 1 N G R U i t I i t R U t       t
  i = 1 N G R D i t I i t R D t       t
i = 1 N G R S i t I i t R S t       t
  i = 1 N G R O i t I i t R O t       t
| S F ( K P P + K W P W K D P D ) | P L

4.1.3. Scenario Constraints

For each simulated scenario, the scenario constraints include system power balance constraints (Equation (28)), wind power limit (Equation (29)), unit generation limits (Equation (30)), system regulation up/down limits (Equations (31) and (32)), system spinning/operating reserve requirements (Equations (33) and (34)), generalized DC power flow constraints (Equation (35)) and permissible adjustment of real power generation limits (Equation (36)):
i = 1 N G p i t s I i t + i = 1 N W p w i t s = P D t + P L t       t , s
0 p w i t s p w i t f , s       i , t
P i , min I i t p i t s P i , max I i t       i , t , s
i = 1 N G R U i t s I i t R U t       t , s
i = 1 N G R D i t s I i t R D t       t , s
i = 1 N G R S i t s I i t R S t       t , s
i = 1 N G R O i t s I i t R O t       t , s
| S F ( K P P s + K W P W s K D P D ) | P L       s
| p i t s p i t | Δ i       i , t , s
Equation (36) ensures that the re-dispatch of coal-fired units can accommodate the volatility of wind power. For each scenario and time interval, Equation (36) is applied for mitigating the violations based on which corrective and preventive actions are formulated. Δi represents the adjustable magnitude of real power of coal-fired units in ten minutes (i.e., 10/60 of hourly ramping rate of units) [47].

4.2. Integrating API Constraints

The differentiated API constraints are expressed in Equation (9). As the violation of API constraints may fail the optimal solution to the formulated MIQP problem, the big M-slack method is employed. The objective function of the proposed model can be expressed as:
t = 1 N T i = 1 N G [ γ i F i ( p i t ) I i t + U i t + D i t ] + t = 1 N T j = 1 N J M t j max { ζ t j , 0 }
A P I ( t ; j ) max { A P I E , t j , A P I 0 , t j } = ζ t j
where Mtj is the penalty function determined by different environmental-justice-associated factors, such as geographical population, regional GDP, ecological GDP, etc., in this paper, Mtj is taken as a large enough constant for simplicity; ζtj is the artificial slack variable which is calculated using Equation (38).
The objective function of the Eco-SCUC is the minimization of total cost including the penalties caused by violation of regional API constraints as expressed in Equation (37), and the constraints include hourly SCUC constraints (Equations (16)–(27)) and wind power scenario constraints (Equations (28)–(36)). The Eco-SCUC incorporating API constraints and volatile wind power can be expressed as:
min (37)
s.t. (16)–(36), (38)

4.3. Solution Methodology

To reduce the computation burden of the formulated MIQP, Benders decomposition technique is applied to unbundle the problem into a master UC problem followed by a network security check sub-problem. Furthermore, re-dispatch will be conducted for accommodating the volatility of wind power in the scenario feasibility check sub-problems. If any violation appears in the sub-problems, the Benders cuts will be created and added to master UC problem. The flowchart of the two-stage algorithm is plotted in Figure 7.
Step 1: the input parameters including unit data, system data, 24-h load forecast, hourly wind power forecast, predicted meteorological conditions, regional background API data are collected;
Step 2: Monte Carlo simulation with LHS is applied to generate hourly wind power volatility scenarios; the number of possible scenarios is reduced by scenario reduction technique;
Step 3: the initial dispatch is checked in the hourly network security check sub-problem 1 in the first stage (part I); if any violation occurs, feasibility Benders cut will constrain the violating unit outputs and will be added to the master UC problem for updating new commitment;
Step 4: the re-dispatch for accommodating wind power volatility is considered in the scenario feasibility check sub-problem 2 in the second stage (part II); if any violation occurs, optimality Benders cut will prevent and revise the commitment solution and will be added to the master UC problem; the scenario check sub-problem is optimized by Equation (36); and
Step 5: the iterations between the master UC problem and feasibility check sub-problem will continue until the generation re-dispatch can satisfy the simulated wind power scenarios.
Detailed descriptions about Benders decomposition can be referred to [13,47]. The formulated MIQP can be solved by multiple commercial solvers [53]. In this paper, the problem is solved using CPLEX 12.1 optimizer in MATLAB (MathWorks, Natick, MA, USA).

5. Case Studies

5.1. Description

The Eco-SCUC is tested for a 24-h period in a 14-bus power system consisting of six6 coal units and wind power. Figure 8 shows the geographical distribution of the units (G1–G6) and regions (R1–R5) with differentiated threshold concentrations. The wind turbines are integrated at Bus 9.
Parameters of system are listed as follows: Table A3 lists the fuel consumption coefficients and stack parameters of coal units; Table A4 shows the generation costs, on/off limits, generation limits, ramping limits of coal units; and Table A5 lists the transmission line impedances and power flow limits. The population in R1–R5 is 2.65, 1.13, 0.56, 0.27, and 0.08 million, respectively.
Forecasted parameters are listed as follows: Table A6 gives the hourly forecasted wind power and system load; Table A7 shows the simulated wind power scenarios (1–10) based on a standard deviation of 10% of the predicted wind power as listed in Table A6; Table A8 provides the hourly forecasted wind speed, wind direction and atmospheric stability class near coal units; and Table A9 lists the forecasted hourly background concentration, threshold concentration, and calculated regional background API of SO2.

5.2. Results

Four cases are studied in this paper:
Case 1: typical SCUC minimizing the total costs, neither wind generation nor API constraints are included; Case 2: typical SCUC minimizing the total costs, considering wind generation, but the API constraints are still absent; Case 3: typical SCUC minimizing the total emissions, considering wind generation, but the API constraints are still absent; and Case 4: the Eco-SCUC, considering both wind generation and API constraints, the API expectation value is set to 100% (APIE,tj = 100%).
Comparisons are made among Cases 1–4. Figure 9 shows the generation results of G1-G6; Figure 10 compares their hourly total emissions; Figure 11 shows the geographical distribution of API increases (ΔAPI) at 6:00, 12:00, 18:00, and 24:00, respectively.
As shown in Figure 9, the hourly generation results of G1-G6 in Cases 1–4 are different. As Case 1 does not consider wind power, it has higher hourly generation compared to Cases 2–4. Cases 2 and 3 have different results because they schedule generation according to the minimization of the fuel cost and total emission, respectively. In comparison, the hourly generation of G3 and G5 in Case 4 (Eco-SCUC) is obviously lower. This is because, according to geographical locations (see Figure 8) and wind directions (see Table A8), G3 and G5 are situated upwind of densely-populated R1 and R2, thus, wind brings pollutants from the units downwind to densely-populated regions. Taking into account differentiated API constraints, Case 4 will curtail generation of G3 and G5 to reduce the resulting API increments in R1 and R2. Conversely, generation of other units is increased to meet the load.
As shown in Figure 10, Case 1 has the highest hourly total emission among the four cases during most times of the day, because it does not contain wind power. Generally, Case 3 has the lowest hourly emissions, because the objective of Case 3 is to minimize total emissions of system. Finally, Case 4 sometimes has higher hourly emissions because it schedules generation according to the differentiated impacts of units on regional air quality, instead of the merely low emission rate of the unit. In other words, the unit that has a low emission rate may have a high impact on the regional API.
As shown in Figure 11, Case 1 without wind power has the worst ΔAPI distribution results compared with Cases 2–4, for example, its ΔAPI in densely-populated R1 is the highest among the four cases. Compared to Cases 1–3, Case 4 (Eco-SCUC) has the lowest ΔAPI in densely-populated regions because, as stated above, Eco-SCUC schedules the generation of units according to their differentiated impacts on regional API. In order to justify the effectiveness of the Eco-SCUC in mitigating regional air pollution, comparisons on different cases, units, and regions are made in Table 1, Table 2 and Table 3.
As shown in Table 1, Case 1 suffers the most total cost, total fuel, total emission, and regional SO2 concentration compared to other cases, this is mainly because there are no wind power and API constraints in Case 1. Furthermore, comparisons among Cases 24 show that Case 4 suffers 2.9% more total cost than Case 2 and 9.7% more total emissions than Case 3. However, the resulting concentrations in Case 4 are much lower: in densely-populated R1, Cav reduces 16.2% from Case 2, and Cpk reduces 12.8% from Case 3. In sparsely populated regions, API constraints in the Eco-SCUC may lead to an increase of pollutant concentration, explaining that in Case 4, Cav and Cpk in R5 increase to 39.8 μg/m3 and 49.6 μg/m3, respectively. This proves that, via shifting generation among units, the Eco-SCUC can cost-effectively mitigate pollution in densely-populated regions while the concentrations in sparsely-populated regions do not exceed local threshold concentrations.
As shown in Table 2, the hourly generation of G1–G6 varies significantly between Case 3 and Case 4. Accordingly, the resulting regional ΔAPI is different (see Figure 11). It can be observed that Case 4 (Eco-SCUC) can effectively mitigate total ΔAPI for densely populated R1–R4 compared to Case 3, while ΔAPI in less densely populated R5 augments from 3.8% to 18.1%. However, Cav and Cpk in R5 in Case 4 are still below its threshold concentration. This phenomenon stems from the integration of differentiated API constraints in the Eco-SCUC. To be specific, G3 and G5 (underlined) are taken as high-impact units and their generation is curtailed, while generation of other units is augmented to meet the load. This justifies that via shifting generation away from high-impact units to low-impact ones, the Eco-SCUC can drive the potential air pollution away from densely-populated regions to sparsely populated ones.
To illustrate the effects in a larger scale, 48 historical samples of meteorological data are used to represent an annual case. Table 3 compares the annual pollutant saturated hours (Ts) in R1–R5 and the annual sum of person-hours (P × Ts) exposed to saturated pollutant for the whole region.
As shown in Table 3, Ts in R1 in Case 4 is nearly half of the value in Cases 13, which indicates that under the Eco-SCUC, densely populated regions face much less air pollution. In return, Ts in R5 in Case 4 is 4 h per year more than Case 1, which has little impact on sparsely populated R5. Once the system includes other regions with much lower population density, Ts in R1–R5 will all decrease. In the whole, the annual sum of person-hour in Case 4 in R1–R5 is 54.0% of that in Case 1, which proves that the Eco-SCUC can protect much more residents from severe air pollution.
From the above results, several conclusions can be drawn:
  • Even with integration of wind power, the API constraints are still necessary in the dispatch model to mitigate air quality deterioration. For example, in Table 1 in R1, Cav in Case 2 is only 0.4% lower than Case 1, while Cav in Case 4 is 16.6% lower than Case 1. Although Case 2 decreases its fuel consumption because of less power demand from coal-fired units, it fails to effectively mitigate regional air quality degradation.
  • The dispatch minimizing total emissions may encourage emission reductions at wrong locations and times regarding regional air quality control. In Table 1 in R1, Cav and Cpk in Case 4 are 9.5% and 12.7% lower than Case 3, respectively. This justifies that both weather conditions and geographical distribution of system are critical in environmental dispatch.
  • The proposed Eco-SCUC takes advantage of wind generation and API constraints, and it can cost-effectively alleviate air quality deterioration for densely populated regions. Via shifting generation among units, the Eco-SCUC succeeds in balancing regional air quality saturation level. Although the Eco-SCUC improves regional air quality at a cost of slight increase of total costs and emissions, it can significantly reduce the person-hour exposure to severe air pollution, which is beneficial to human and environmental health.
We examine the effectiveness of the model with respect to differentiated API expectations with and without wind power. API expectation (APIE,tj) refers to the expected value of API as expressed in Equation (9). When the city tightens its air quality control, APIE,tj can be set to lower values; when the city loosens its control, APIE,tj can be set to higher values. We adjust APIE,tj from 50% to 200%.
Figure 12 shows the resulting Cav and Cpk in R1 under the Eco-SCUC, respectively. As APIE,tj decreases, the model tightens control on emissions, both Cav and Cpk drop. This demonstrates that the decrease of APIE,tj helps improve regional air quality. Especially during extreme haze-fog periods, APIE,tj could be set to relatively low value to cope with bad weather. Moreover, Figure 12 proves that the Eco-SCUC with wind power is more effective than the counterpart without wind power in reducing regional pollutant concentrations.

5.3. Discussion

This paper links regional API with electricity generation scheduling. The proposed Eco-SCUC is a useful step towards alleviating the environmental dilemma in the power sector, as it enforces differentiated API constraints above which human and ecosystem health faces severe challenges. Based on the case studies, several potential policy implementations can be drawn:
  • Environmental or ecological taxes can be imposed to power plants based on the violation of API constraints. In the Eco-SCUC, the violation of API constraint can be directly linked to the ecosystem or human health damages. This solves the problem that the impact of coal-fired power plants on regional air quality cannot be quantitatively analyzed.
  • Air quality compensation mechanism can be implemented by the government to press the power plants to adjust their generation in response to regional air quality. The Eco-SCUC helps calculate the compensated amount based on regional API changes.
  • The developing cyber-physical system (CPS) technologies enable power grid corporations to expand their ancillary services to include API-associated regulations, thereby improving regional air quality.

6. Conclusions

This paper proposes an Eco-SCUC model addressing regional air quality management via shifting generation among units. An API index is defined for evaluating regional air pollutant saturation levels and is calculated based on a Gaussian puff dispersion model, geographically-differentiated AEC, and forecasted weather conditions. API constraints are formulated according to the differentiated impacts of coal-fired units on regional air quality and integrated into a SCUC model. Stochastic optimization is applied to accommodate volatile wind power; as such the benefits of wind power for air quality control are estimated. The test results demonstrate that the model can cost-effectively prevent air quality deterioration in densely-populated regions and reduce the person-hours exposed to severe air pollution.
Future studies will focus on the investigation of chemical conversion from primary pollutants to secondary pollutants and its integration into the Eco-SCUC model.

Author Contributions

J.Y. conceived the idea of this paper and guided the research; D.G. designed the Eco-SCUC and wrote the paper; and M.B. proposed the solution to the optimization problem and debugged the program. All of the authors revised the manuscript.

Acknowledgments

This research was funded by National Natural Science Foundation of China (51377035) and the National Natural Science Foundation of China (51361130153). The authors also would like to thank the editor and the anonymous reviewers for their valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

A. Acronyms
EDEconomic dispatch
UCUnit commitment
SCUCSecurity-constrained unit commitment
Eco-SCUCEcology- and security-constrained unit commitment
AECAtmospheric environmental capacity
APIAir pollutant intensity
LHSLatin hypercube sampling
CDFCumulative distribution function
PDFProbability density function
MIQPMixed-integer quadratic programming
B. Sets and indices
tIndex for time
τIndex for pollutant emitting time
iIndex for unit
jIndex for monitoring point
sIndex for wind power generation scenario
C. Parameters
(xp,yp,zp)Coordinates of dispersion puff center
(xs,ys,zs)Coordinates of the emitting stack
(xj, yj, zj)Coordinates of concentration monitoring location
Ψ, λ, hIntermediary parameter
σx/y/z(·)Standard deviation in x/y/z direction of the puff center
ΔnEqual time interval of dispersion period (t-τ), n = 1, 2, …N
vx/y/z(n)Wind velocity in x/y/z direction during time interval Δn
ϖAtmospheric stability class A–F
a(·)/b(·)/α(·)/β(·)Atmospheric stability coefficients
siPollutant emission per unit fuel of unit i, ton/ton
γiFuel cost of unit i, $/ton
FiFuel consumption function of unit i, ton
ai, bi, ciFuel consumption coefficients of unit i, ton/MW2h, ton/mwh, ton/h
c ¯ Volume-average value of concentration c, μg/m3
c ¯ 0 Initial volume-average concentration, μg/m3
ρAtmospheric volume, m3
V / V t Advection/pseudo-diffusion speed, m/s
ηρRemoval rate of air pollutant in the examined volume ρ, ton/a
ηtjRemoval rate of air pollutant at location j and time t, ton/a
δρDirac function for the examined volume ρ
v d / v w Dry/wet deposition speed, m/s
AAEC coefficient, m2/s
UPredicted wind speed, m/s
HHeight of atmospheric boundary layer, m
SArea of the examined region, km2
vValue of random wind speed variable V, m/s
vmMth sampled wind speed value, m/s
v ¯ /sAverage/sample variance of measured wind speed, m/s
FV(v)/fV(v)CDF/PDF of random V
c/kScale/shape factor at the predefined region
μ/ σ v 2 Mean/standard variance of wind speed
Γ(·)Gamma function
pwValue of random wind power variable PW, MW
pwrRated wind power, MW
vr, vin, voutRated/cut-in/cut-out wind speed, m/s
NTNumber of periods under study (24-h)
NGNumber of coal units
NWNumber of wind-powered units
NMSample number of measured wind speed
NJNumber of air quality monitoring points
ΔTTime interval of dispatch (1 h)
p w i t f Forecasted generation of wind-powered unit i at time t, MW
Uit/DitStartup/shutdown cost of unit i at time t, $
PDt/PLtSystem demand/losses at time t, MW
RSt/ROtSystem spinning/operating reserve requirement at time t, MW
T i o n / o f f Minimum on/off time of unit i
Pi,min/Pi,maxLower/upper limit of real power generation of unit i, MW
P i R U / P i R D Ramping up/down limit of unit i, MW
RUt/RDtSystem regulation up/down limit at time t, MW
SFShift factor matrix
KP/W/DBus-unit/bus-wind/bus-load incidence matrix
PD/PLSystem demand/losses vector
MtjPenalty function for location j at time t
ζtjArtificial slack variable for location j at time t
D. Variables
Gi(·)Gaussian puff dispersion function
Ci(t;j)Change of pollutant concentration caused by unit i at time t at location j, μg/m3
CT,tjThreshold concentration at location j and time t, μg/m3
C0,tjBackground concentration at location j and time t, μg/m3
Cav/CpkDaily average/peak concentration of pollutant, μg/m3
Qi(t)Pollutant emission of unit i at time t, ton
pitGeneration of coal-fired unit i at time t, MW
pwitGeneration of wind-powered unit i at time t, MW
API(t;j)API at location j and time t, %
APIi(t;j)Change of API caused by unit i at location j and time t, %
API0,tjBackground API at location j and time t, %
APIE,tjExpected API at location j and time t, %
IitBinary commitment of unit i at time t
X i t o n / o f f On/off time of unit i at time t
RUit/RDitRegulation up/down capacity of unit i at time t, MW
RSit/ROitSpinning/operating reserve of unit i at time t, MW
PVector of generation of coal-fired units
PWVector of generation of wind-powered units

Appendix A

Appendix A.1

The Gaussian puff model is detailed as:
Ψ = [ x x p ( τ , t ) σ x ( τ , t ) ] 2 + [ y y p ( τ , t ) σ y ( τ , t ) ] 2 + [ z z p ( τ , t ) σ z ( τ , t ) ] 2
[ x p ( τ , t ) y p ( τ , t ) z p ( τ , t ) ] = [ x s y s z s ] + n = 1 N [ v x ( n ) v y ( n ) v z ( n ) ] Δ n
{ σ x ( τ , t ) = σ y ( τ , t ) = α ( ϖ ) ( t τ ) a ( ϖ ) σ z ( τ , t ) = β ( ϖ ) ( t τ ) b ( ϖ )
{ σ x ( τ , n ) = [ a ( n , ϖ ) 1 / α ( n , ϖ ) Δ n + σ x ( τ , n 1 ) 1 / α ( n , ϖ ) ] α ( n , ϖ ) σ y ( τ , n ) = [ a ( n , ϖ ) 1 / α ( n , ϖ ) Δ n + σ y ( τ , n 1 ) 1 / α ( n , ϖ ) ] α ( n , ϖ ) σ z ( τ , n ) = [ b ( n , ϖ ) 1 / β ( n , ϖ ) Δ n + σ z ( τ , n 1 ) 1 / β ( n , ϖ ) ] β ( n , ϖ )

Appendix A.2

The atmospheric environmental capacity (AEC) model is detailed as:
c ¯ t + 1 ρ ρ V c d ρ = 1 ρ ρ η ρ δ ρ d ρ 1 ρ ρ ( c v d + c v w ) d ρ 1 ρ ρ V t c d ρ
c ¯ = 1 ρ ρ c d ρ
c ¯ = η ρ A S [ 1 exp ( A H S t ) ] + c ¯ 0 exp ( A H S t )
A = π 2 U H + ( v d + v w ) S

Appendix A.3

The wind speed distribution model is detailed as:
μ = c Γ ( 1 + 1 k )
σ v 2 = c 2 Γ ( 1 + 2 k ) μ 2
Γ ( α ) = 0 y α 1 exp ( y ) d y , y > 0
μ = v ¯ = 1 N M m = 1 N M v m
σ v 2 = s 2 = 1 N M m = 1 N M ( v m v ¯ ) 2
k = 1 λ 1
c = μ Γ 1 λ
λ = 1 Γ ( σ v 2 μ 2 + 1 )

Appendix A.4

The wind power distribution model is detailed as:
P W = { 0 ( V < v in   o r   V v out ) p w r   ( v r V < v out ) ( V v in ) p w r / ( v r v in ) ( v in V < v r )
Pr ( P W = 0 ) = Pr ( V < v in ) + Pr ( V v out ) = 1 exp [ ( v in c ) k ] + exp [ ( v out c ) k ]
Pr ( P W = p w r ) = Pr ( v r V < v out ) = exp [ ( v r c ) k ] exp [ ( v out c ) k ]

Appendix B

Figure A1. Relationship between standard deviation of puff center σx/y/z(τ,t) and dispersion time (tτ) according to Pasquill atmospheric stability classes (A–F): (a) in the x and y directions; and (b) in the z direction.
Figure A1. Relationship between standard deviation of puff center σx/y/z(τ,t) and dispersion time (tτ) according to Pasquill atmospheric stability classes (A–F): (a) in the x and y directions; and (b) in the z direction.
Sustainability 10 01433 g0a1
Figure A2. Curve of wind power with respect to wind speed.
Figure A2. Curve of wind power with respect to wind speed.
Sustainability 10 01433 g0a2
Table A1. Solar radiation classification according to total/low cloud coverage and solar altitude.
Table A1. Solar radiation classification according to total/low cloud coverage and solar altitude.
Total/Low Cloud Coverage (1/10)Solar Altitude (°)
Night≤1515–3535–65>65
≤4/≤4−2−1+1+2+3
5–7/≤4−10+1+2+3
≥8/≤4−100+1+1
≥5/5–70000+1
≥8/≥80000+1
Table A2. Atmospheric stability classification according to wind speed and solar radiation level.
Table A2. Atmospheric stability classification according to wind speed and solar radiation level.
Ground-Level Wind Speed(m/s)Solar Radiation Level
+3+2+10−1−2
≤1.9AA–BBDEF
2–2.9A–BBCDEF
3–4.9BB–CCDDE
5–5.9CC–DDDDD
≥6CDDDDD
Table A3. Fuel consumption coefficients and stack parameters of coal-fired units.
Table A3. Fuel consumption coefficients and stack parameters of coal-fired units.
UnitFuel Consumption CoefficientsStack Parameters
ai (ton/MW2h)bi (ton/MWh)ci (ton/h)si (ton/ton)xs (km)ys (km)zs (m)
G10.0010760.05057221.20.00487314150
G20.0009600.05005224.80.00465239150
G30.0010520.05607216.00.0045572135
G40.0011320.06394016.80.004952224105
G50.0011760.06886013.20.00425155135
G60.0015280.06986412.40.00413818135
Table A4. Generation costs, on/off limits, generation limits, and ramping limits of units.
Table A4. Generation costs, on/off limits, generation limits, and ramping limits of units.
UnitGeneration CostsOn/Off LimitsGeneration LimitsRamp Limits
γi ($/ton)Uit ($)Dit ($)Ton (h)Toff (h)Pi,min (MW)Pi,max (MW) P i R U / P i R D (MW)
G155.572623.6374.15536080120/120
G255.884501.2213.055380100125/125
G360.360497.6186.6442907580/80
G460.624428.7264.344280110115/115
G560.836423.5236.733295120120/120
G663.712202.5113.6221202560/60
Table A5. Transmission line impedances and power flow limits.
Table A5. Transmission line impedances and power flow limits.
LineInitial BusTerminal BusImpedance (p.u.)Power Flow Limit (MW)
1120.170250
2150.258180
3230.150175
4240.197185
5250.140225
6340.018320
7450.037175
8470.037195
9490.150200
10560.187225
116110.197450
126120.197150
136130.150180
14780.140380
15790.019260
169100.039325
179140.037255
1810110.152250
1912130.183150
2013140.192260
Table A6. Hourly forecasted wind power and system load.
Table A6. Hourly forecasted wind power and system load.
t p w i t f (MW)Bus-Load (MW)
Bus 4Bus 5Bus 7Bus 9Bus 10Bus 11Bus 12Bus 13
1142. 6190.2113.183.798.493.8177.1107.155.2
2131.0204.2121.489.9105.6100.8190.1115.059.1
3123.8203.6121.189.6105.3100.5189.5114.758.8
4142. 6205.4122.190.4106.2101.3191.2115.759.4
5123.8207.4123.391.3107.3102.3193.1116.859.9
6105.1212.1126.193.3119.7104.7197.5119.561.3
796.5237.2141.0104.4122.7117.0220.8133.668.5
889.3252.0149.8110.9130.3124.3215.9151.982.8
9110.9236.0146.3108.2127.3121.4225.3138.671.1
10112.3234.0135.099.9117.5111.4227.8137.470.5
11128.2267.7145.1122.2133.6137.0221.3156.473.8
12148.3254.7149.2125.2137.2140.5247.1160.369.4
13207.4263.0148.2124.5126.4139.6235.5154.172.6
14198.7260.8136.9123.5135.2138.5234.2158.174.7
15193.0263.6128.6124.8136.7139.9236.1154.466.4
16185.7249.7118.4109.9129.2123.2223.1140.672.1
17184.3225.6114.199.3116.7121.3220.0127.165.2
18174.2214.1127.394.2129.7105.6199.3120.661.9
1963.5246.2136.4108.3127.3121.5229.2165.885.1
2086.5206.0122.490.6106.5101.6191.7116.059.5
21110.9211.1125.592.9109.2104.2196.5118.961.0
22119.5245.7146.0108.1127.1121.2228.7138.471.0
23138.2197.4117.486.9102.197.4183.8111.257.0
24146.9191.7114.084.499.294.6181.3108.155.4
Table A7. Hourly simulated and reduced wind power volatility scenarios.
Table A7. Hourly simulated and reduced wind power volatility scenarios.
tWind Power Volatility Scenario (MW)
12345678910
1144.9146.8148.7146.8146.6146.9143.7137.9141.3140.1
2133.7135.4133.9136.8133.5128.6126.9128.3129.8130.1
3121.5120.9122.0121.1120.2119.7118.9114.3113.5114.4
4145.6146.8148.0145.8143.4143.4145.8148.7149.4148.8
5124.7122.7120.7119.2121.0121.7118.8125.2122.8124.8
6103.5104.1102.2102.8101.999.497.292.494.395.4
795.593.393.393.293.491.189.590.892.694.8
889.688.790.890.288.789.590.289.891.395.5
9113.5111.7110.2112.1112.7117.0118.4124.6123.2122.2
10115.4112.7113.2113.7112.4112.1113.0121.8122.4117.0
11126.1128.2126.5126.9127.9130.6130.3135.1131.9132.7
12151.8153.3155.3158.6160.2156.5156.9147.7147.2149.3
13212.2210.3207.8205.6208.2219.2197.9212.4210.5209.7
14198.6203.1203.2205.9205.4200.5203.1202.3199.9195.9
15195.9191.4193.4195.9191.9184.9182.1177.3170.3168.4
16182.4181.9185.5174.4182.0181.8189.8193.6192.9189.6
17183.6182.6186.8187.5191.4191.2188.2187.9184.1182.1
18177.9180.3180.8177.0174.0172.3171.2167.2168.1167.6
1964.560.564.469.364.162.763.164.664.662.7
2088.587.285.785.986.180.287.382.983.883.7
21111.8111.8110.5112.1114.9113.6111.3109.5110.6115.7
22116.8116.5118.5121.1118.6120.5123.1131.3132.3122.5
23140.7141.8140.1137.6137.3140.6142.6138.8135.6128.7
24150.1151.7154.1154.7151.7159.6151.5147.4144.3145.3
Table A8. Hourly forecasted wind speed, wind direction, and atmospheric stability near units.
Table A8. Hourly forecasted wind speed, wind direction, and atmospheric stability near units.
tG1G2G3G4G5G6
v1φ1ϖ1vφϖvφϖvφϖvφϖvφϖ
17.6292D15.2290E3.2310E13.1290D13.220E8.6337F
28.3285D15.4295E4.3333E13.4284D14.218E8.4343F
37.8281D15.0301E4.3323E12.9280D13.515E8.5349F
46.9275D14.2304E4.6320E12.4283D13.77E7.9356F
56.4276D14.9310E4.2317E12.1286D13.09E7.5348F
66.2271D14.0312E5.3304E11.5291D12.24E6.5355F
75.3270F14.1318F5.6302D10.5283F11.0354F5.5357D
85.0270F14.1317F4.9305D10.4276F11.0350F5.3355D
95.1274F14.3316F6.3311D10.2273F11.2359F5355D
105.2277F14.5315F5.3313D10.7269F11.32F5.15D
115.4275F14.6315F4.3315D10.3268F11.18F5.34D
125.0269F14.3320F5.1305D10.2271F11.5355F5.4358D
135.3268F14.2310E6.3298F10.2270C11.6357E5.93F
145.3270F14.2319E6.4295F10.3275C10.8354E5.7356F
154.9274F14.5316E7.3290F10.4276C10.97E5.12F
165.1275F14.3315E6.3285F10.1273C10.83E5.23F
175.0272F14.3315E5.3294F10.1275C11.1355E5355F
185.5274F14.1319D4.2303E10.2276C11.0355D5.2360D
195.5279F15.5310D6.3300E12.4273C11.8358D5.6355D
206.1281F15.9308D5.9322E12.6276C11.9360D6.3352D
216.4284F16.5304D5.5318E11.9284C12.5360D6.5356D
227.4295F17.5310D5.9323E11.0297C13.4345D7.0354D
237.9292D17.9296D6.0325E12.6292C14.5339D7.5331D
248.0296D16.5292D6.3316E12.9302C14.6331D7.2345D
1 In this table, v represents the wind speed, in m/s; φ represents the wind direction, in °; and ϖ represents the Pasquill atmospheric stability class, in A–F.
Table A9. Hourly forecasted regional background and threshold concentration of SO2.
Table A9. Hourly forecasted regional background and threshold concentration of SO2.
tBackground Concentration C0,tj (μg/m3)Threshold Concentration CT,tj (μg/m3)Background API API0.tj (%)
R1 1R2R3R4R5R1R2R3R4R5R1R2R3R4R5
134.532.231.928.825.335.738.240.342.560.296.684.379.167.842.0
233.231.629.826.723.235.738.240.342.560.293.082.773.962.838.5
331.530.228.725.521.935.738.240.342.560.288.279.171.260.036.4
430.829.828.825.721.535.337.539.842.360.287.379.572.460.835.7
531.530.328.825.822.535.337.539.842.360.289.280.872.461.037.4
633.632.429.326.123.635.337.539.842.360.295.286.473.661.739.2
736.433.530.726.924.335.137.539.542.259.5103.789.377.763.740.8
839.435.832.528.525.135.137.539.542.259.5112.395.582.367.542.2
943.939.134.829.826.835.137.539.542.259.5125.1104.388.170.645.0
1048.542.738.932.427.934.436.939.541.659.5141.0115.798.577.946.9
1152.946.842.634.531.434.436.939.541.659.5153.8126.8107.882.952.8
1253.148.343.836.834.234.436.939.541.659.5154.4130.9110.988.557.5
1352.948.143.536.533.634.636.939.441.659.3152.9130.4110.487.756.7
1451.845.242.736.433.534.636.939.441.659.3149.7122.5108.487.556.5
1548.245.141.635.932.834.636.939.441.659.3139.3122.2105.686.355.3
1648.144.141.735.231.934.536.640.241.559.3139.4120.5103.784.853.8
1750.143.240.735.032.434.536.640.241.559.3145.2118.0101.284.354.6
1854.247.243.836.634.734.536.640.241.559.3157.1129.0109.088.258.5
1953.946.943.235.435.134.537.240.241.560.1156.2126.1107.585.358.4
2050.244.341.534.233.234.537.240.241.560.1145.5119.1103.282.455.2
2147.042.138.932.432.434.537.240.241.560.1136.2113.196.878.153.9
2243.138.235.130.130.135.237.540.442.460.1122.4101.986.971.050.1
2340.736.232.527.927.535.237.540.442.460.1115.696.580.465.845.8
2438.134.632.328.326.635.237.540.442.460.1108.292.380.066.744.3
1 In this table, R1–R5 represents the region 1–5 as shown in Figure 8.

References

  1. Zhao, Y.; Wang, S.; Duan, L.; Lei, Y.; Cao, P.; Hao, J. Primary air pollutant emissions of coal-fired power plants in China: Current status and future prediction. Atmos. Environ. 2008, 42, 8442–8452. [Google Scholar] [CrossRef]
  2. Li, X.; Chen, X.; Yuan, X.; Zeng, G.; León, T.; Liang, J.; Chen, G.; Yuan, X. Characteristics of particulate pollution (PM2.5 and PM10) and their spacescale-dependent relationships with meteorological elements in China. Sustainability 2017, 9, 2330. [Google Scholar] [CrossRef]
  3. Zhao, B.; Wang, S.; Dong, X.; Wang, J.; Duan, L.; Fu, X.; Hao, J.; Fu, J. Environmental effects of the recent emission changes in China: Implications for particulate matter pollution and soil acidification. Environ. Res. Lett. 2013, 8, 024031. [Google Scholar] [CrossRef]
  4. Zhou, Q.; Yabar, H.; Mizunoya, T.; Higano, Y. Evaluation of integrated air pollution and climate change policies: Case study in the thermal power sector in Chongqing City, China. Sustainability 2017, 9, 1741. [Google Scholar] [CrossRef]
  5. Chen, L.; Shi, M.; Li, S.; Gao, S.; Zhang, H.; Sun, Y.; Mao, J.; Bai, Z.; Wang, Z.; Zhou, J. Quantifying public health benefits of environmental strategy of PM2.5 air quality management in Beijing-Tianjin-Hebei region, China. J. Environ. Sci. 2017, 57, 33–40. [Google Scholar] [CrossRef] [PubMed]
  6. Pretorious, I.; Piketh, S.; Burger, R. Emissions management and health exposure: Should all power stations be treated equal? Air Qual. Atmos. Health 2017, 10, 509–520. [Google Scholar] [CrossRef]
  7. Song, C.; He, J.; Wu, L.; Jin, T.; Chen, X.; Li, R.; Ren, P.; Zhang, L.; Mao, H. Health burden attributable to ambient PM2.5 in China. Environ. Pollut. 2017, 223, 575–586. [Google Scholar] [CrossRef] [PubMed]
  8. Li, L.; Liu, D. Study on an air quality evaluation model for Beijing city under haze-fog pollution based on new ambient air quality standards. Int. J. Environ. Res. Public Health 2014, 11, 8909–8923. [Google Scholar] [CrossRef] [PubMed]
  9. Li, L.; Tang, D.; Kong, Y.; Yang, Y.; Liu, D. Spatial analysis of haze–fog pollution in China. Energy Environ. 2016, 27, 726–740. [Google Scholar] [CrossRef]
  10. Zhao, H.; Ma, W.; Dong, H.; Jiang, P. Analysis of co-effects on air pollutants and CO2 emissions generated by end-of-pipe measures of pollution control in China’s coal-fired power plants. Sustainability 2017, 9, 499. [Google Scholar] [CrossRef]
  11. Mollahassani-Pour, M.; Rashidinejad, M.; Abdollahi, A. Appraisal of eco-friendly preventive maintenance scheduling strategy impacts on GHG emissions mitigation in smart grids. J. Clean. Prod. 2017, 143, 212–223. [Google Scholar] [CrossRef]
  12. Chinese Coal-Fired Electricity Generation Expected to Flatten as Mix Shifts to Renewables. Available online: https://www.eia.gov/todayinenergy/detail.php?id=33092 (accessed on 29 September 2017).
  13. Wang, R.; Wang, J.; Guan, Y. Robust unit commitment with wind power and pumped storage hydro. IEEE Trans. Power Syst. 2012, 27, 800–810. [Google Scholar] [CrossRef]
  14. Xu, Q.; Ding, Y.; Zheng, A. An optimal dispatch model of wind-integrated power system considering demand response and reliability. Sustainability 2017, 9, 758. [Google Scholar] [CrossRef]
  15. Liu, X.; Xu, W. Minimum emission dispatch constrained by stochastic wind power availability and cost. IEEE Trans. Power Syst. 2010, 25, 1705–1713. [Google Scholar] [CrossRef]
  16. Lyu, X.; Shi, A. Research on the renewable energy industry financing efficiency assessment and mode selection. Sustainability 2018, 10, 222. [Google Scholar] [CrossRef]
  17. Wang, X.; Gong, Y.; Jiang, C. Regional carbon emission management based on probabilistic power flow with correlated stochastic variables. IEEE Trans. Power Syst. 2015, 30, 1094–1103. [Google Scholar] [CrossRef]
  18. Wang, Y.; Huang, Y.; Wang, Y.; Li, F.; Zhang, Y.; Tian, C. Operation optimization in a smart micro-grid in the presence of distributed generation and demand response. Sustainability 2018, 10, 847. [Google Scholar] [CrossRef]
  19. Zhu, Z.; Liu, D.; Liao, Q.; Tang, F.; Zhang, J.; Jiang, H. Optimal power scheduling for a medium voltage AC/DC hybrid distribution network. Sustainability 2018, 10, 318. [Google Scholar] [CrossRef]
  20. Wei, W.; Liu, F.; Wang, J.; Chen, L.; Mei, S.; Yuan, T. Robust environmental-economic dispatch incorporating wind power generation and carbon capture plants. Appl. Energy 2016, 183, 674–684. [Google Scholar] [CrossRef]
  21. Zhang, N.; Hu, Z.; Dai, D.; Dang, S.; Yao, M.; Zhou, Y. Unit commitment model in smart grid environment considering carbon emissions trading. IEEE Trans. Smart Grid 2016, 7, 420–427. [Google Scholar] [CrossRef]
  22. Abdollahi, A.; Moghaddam, M.P.; Rashidinejad, M.; Sheikh-El-Eslami, M.K. Investigation of economic and environmental-driven demand response measures incorporating UC. IEEE Trans. Smart Grid 2012, 3, 12–25. [Google Scholar] [CrossRef]
  23. Zhao, C.; Wang, J.; Watson, J.; Guan, Y. Multi-stage robust unit commitment considering wind and demand response uncertainties. IEEE Trans. Power Syst. 2013, 28, 2708–2717. [Google Scholar] [CrossRef]
  24. Shi, N.; Luo, Y. Energy storage system sizing based on a reliability assessment of power systems integrated with wind power. Sustainability 2017, 9, 395. [Google Scholar] [CrossRef]
  25. Mou, D. Wind power development and energy storage under China’s electricity market reform—A case study of Fujian Province. Sustainability 2018, 10, 298. [Google Scholar] [CrossRef]
  26. Ban, M.; Yu, J. Procedural simulation method for aggregating charging load model of private electric vehicle cluster. J. Mod. Power Syst. Clean Energy 2015, 3, 170–179. [Google Scholar] [CrossRef]
  27. Ye, J.; Yuan, R. Integrated natural gas, heat, and power dispatch considering wind power and power-to-gas. Sustainability 2017, 9, 602. [Google Scholar] [CrossRef]
  28. Ban, M.; Yu, J.; Shalidehpour, M.; Yao, Y. Integration of power-to-hydrogen in day-ahead security-constrained unit commitment with high wind penetration. J. Mod. Power Syst. Clean Energy 2017, 5, 337–349. [Google Scholar] [CrossRef]
  29. Zhou, Y.; Levy, J.I.; Evans, J.S.; Hammitt, J.K. The influence of geographic location on population exposure to emissions from power plants throughout China. Environ. Int. 2006, 32, 365–373. [Google Scholar] [CrossRef] [PubMed]
  30. Cao, X.; Roy, G.; Hurley, W.J.; Andrews, W.S. Dispersion coefficients for Gaussian puff models. Bound.-Layer Meteorol. 2011, 139, 487–500. [Google Scholar] [CrossRef]
  31. Fann, N.; Fulcher, C.M.; Hubbell, B.J. The influence of location, source, and emission type in estimates of the human health benefits of reducing a ton of air pollution. Air Qual. Atmos. Health 2009, 2, 169–176. [Google Scholar] [CrossRef] [PubMed]
  32. Sullivan, R.L.; Hackett, D.F. Air quality control using a minimum pollution-dispatching algorithm. Environ. Sci. Technol. 1973, 7, 1019–1022. [Google Scholar] [CrossRef] [PubMed]
  33. Schweizer, P.F. Determining optimal fuel mix for environmental dispatch. IEEE Trans. Autom. Control 1974, 19, 534–537. [Google Scholar] [CrossRef]
  34. Chu, K.; Jamshidi, M.; Levitan, R.E. An approach to on-line power dispatch with ambient air pollution constraints. IEEE Trans. Autom. Control 1977, 22, 385–396. [Google Scholar] [CrossRef]
  35. Levy, J.I.; Spengler, J.D.; Hlinka, D.; Sullivan, D.; Moon, D. Using CALPUFF to evaluate the impacts of power plant emissions in Illinois: Model sensitivity and implications. Atmos. Environ. 2002, 36, 1063–1075. [Google Scholar] [CrossRef]
  36. Napelenok, S.L.; Cohan, D.S.; Hu, Y.; Russell, A.G. Decoupled direct 3D sensitivity analysis for particulate matter (DDM-3D/PM). Atmos. Environ. 2006, 40, 6112–6121. [Google Scholar] [CrossRef]
  37. Duan, L.; Xie, S.; Zhou, Z.; Ye, X.; Hao, J. Calculation and mapping of critical loads for S, N and Acidity in China. Water Air Soil Pollut. 2001, 130, 1199–1204. [Google Scholar] [CrossRef]
  38. Porter, E.; Blett, T.; Potter, D.U.; Huber, C. Protecting resources on Federal Lands: Implications of critical loads for atmospheric deposition of nitrogen and sulfur. Bioscience 2005, 55, 603–612. [Google Scholar] [CrossRef]
  39. Pingxing, L.; Jie, F. Regional ecological vulnerability assessment of the Guangxi Xijiang River Economic Belt in Southwest China with VSD model. J. Resour. Ecol. 2014, 5, 163–170. [Google Scholar] [CrossRef]
  40. Wang, M.; Liu, J.; Wang, J.; Zhao, G. Ecological footprint and major driving forces in West Jilin Province, Northeast China. Chin. Geogr. Sci. 2010, 20, 434–441. [Google Scholar] [CrossRef]
  41. Yanxi, C. Study of the Atmospheric Environment Capacity in Harbin City. Master’s Thesis, Harbin Institute of Technology, Harbin, China, 2015. (In Chinese). [Google Scholar]
  42. Zhou, Y.; Zhou, J. Urban atmospheric environmental capacity and atmospheric environmental carrying capacity constrained by GDP–PM2.5. Ecol. Indic. 2017, 73, 637–652. [Google Scholar] [CrossRef]
  43. Xu, D.; Wang, Y.; Zhu, R. Atmospheric environmental capacity and urban atmospheric load in mainland China. Sci. China Earth Sci. 2018, 61, 33–46. [Google Scholar] [CrossRef]
  44. Zoras, S.; Triantafyllou, A.G.; Deligiorgi, D. Atmospheric stability and PM10 concentrations at far distance from elevated point sources in complex terrain: Worst-case episode study. J. Environ. Manag. 2006, 80, 295–302. [Google Scholar] [CrossRef] [PubMed]
  45. An, X.; Zuo, H.; Chen, L. Atmospheric environmental capacity of SO2 in winter over Lanzhou in China: A case study. Adv. Atmos. Sci. 2007, 24, 688–699. [Google Scholar] [CrossRef]
  46. Air Quality Index for Harbin City. Available online: http://www.pm25.in/haerbin (accessed on 24 January 2018).
  47. Wang, J.; Shahidehpour, M.; Li, Z. Security-constrained unit commitment with volatile wind power generation. IEEE Trans. Power Syst. 2008, 23, 1319–1327. [Google Scholar] [CrossRef]
  48. Masters, G.M. Renewable and Efficient Electric Power Systems; Wiley: New York, NY, USA, 2004. [Google Scholar]
  49. Karki, R.; Patel, J. Reliability assessment of a wind power delivery system. Proc. Inst. Mech. Eng. Part O J. Risk Reliab. 2008, 223, 51–58. [Google Scholar] [CrossRef]
  50. Hetzer, J.; Yu, D.C.; Bhattarai, K. An economic dispatch model incorporating wind power. IEEE Trans. Energy Convers. 2008, 23, 603–611. [Google Scholar] [CrossRef]
  51. Hajipour, E.; Bozorg, M.; Fotuhi-Firuzabad, M. Stochastic capacity expansion planning of remote microgrids with wind farms and energy storage. IEEE Trans. Sustain. Energy 2015, 6, 491–498. [Google Scholar] [CrossRef]
  52. Wu, L.; Shahidehpour, M.; Li, T. Stochastic security-constrained unit commitment. IEEE Trans. Power Syst. 2007, 22, 800–811. [Google Scholar] [CrossRef]
  53. Lopez, J.A.; Ceciliano-Meza, J.L.; Moya, I.G.; Gomez, R.N. A MIQCP formulation to solve the unit commitment problem for large-scale power systems. Int. J. Electr. Power 2012, 36, 68–75. [Google Scholar] [CrossRef]
Figure 1. Cyber-physical framework of the proposed Eco-SCUC. Coal-fired power plant A and B are assumed to be identical near-city power plants. Remote power plants represent power plants near coal-mines which are located in sparsely-populated or unmanned regions. Blue dots represent air quality ground stations; and blue lines and green lines represent inputs and outputs for the dispatching center.
Figure 1. Cyber-physical framework of the proposed Eco-SCUC. Coal-fired power plant A and B are assumed to be identical near-city power plants. Remote power plants represent power plants near coal-mines which are located in sparsely-populated or unmanned regions. Blue dots represent air quality ground stations; and blue lines and green lines represent inputs and outputs for the dispatching center.
Sustainability 10 01433 g001
Figure 2. Methodology of formulating differentiated API constraints. Meteorological data are inputs for Gaussian puff model, and the generation-associated concentrations Ci are calculated; the regional threshold concentrations CT are derived according to AEC model. Based on Ci and CT, regional APIs are obtained and differentiated constraints are formed according to regional background API value.
Figure 2. Methodology of formulating differentiated API constraints. Meteorological data are inputs for Gaussian puff model, and the generation-associated concentrations Ci are calculated; the regional threshold concentrations CT are derived according to AEC model. Based on Ci and CT, regional APIs are obtained and differentiated constraints are formed according to regional background API value.
Sustainability 10 01433 g002
Figure 3. Top view of the Gaussian puff dispersion model.
Figure 3. Top view of the Gaussian puff dispersion model.
Sustainability 10 01433 g003
Figure 4. (a) SO2 threshold concentration and (b) API of SO2 by administrative-district-level of one city in Northeast China. Daily-average concentration data are from [46].
Figure 4. (a) SO2 threshold concentration and (b) API of SO2 by administrative-district-level of one city in Northeast China. Daily-average concentration data are from [46].
Sustainability 10 01433 g004
Figure 5. CDF and PDF of wind generation PW: (a) CDF of PW; (b) PDF of PW. Assumptions: c = 15, vin = 5 m/s, vout = 45 m/s, vr = 15 m/s, k = 1.2, 1.7, or 2.2. CDF includes both discrete probabilities (marked as colorful dots) and continuous probabilities, while PDF only contains continuous portions.
Figure 5. CDF and PDF of wind generation PW: (a) CDF of PW; (b) PDF of PW. Assumptions: c = 15, vin = 5 m/s, vout = 45 m/s, vr = 15 m/s, k = 1.2, 1.7, or 2.2. CDF includes both discrete probabilities (marked as colorful dots) and continuous probabilities, while PDF only contains continuous portions.
Sustainability 10 01433 g005
Figure 6. Diagram of the Eco-SCUC model. Blue represents the electricity network; red represents the environment network; the full lines represent physical structure; and the dotted lines represent data and information.
Figure 6. Diagram of the Eco-SCUC model. Blue represents the electricity network; red represents the environment network; the full lines represent physical structure; and the dotted lines represent data and information.
Sustainability 10 01433 g006
Figure 7. Eco-SCUC solution methodology. Part I includes the master UC problem and hourly network security check sub-problem; Part II includes the wind power scenario feasibility check sub-problem.
Figure 7. Eco-SCUC solution methodology. Part I includes the master UC problem and hourly network security check sub-problem; Part II includes the wind power scenario feasibility check sub-problem.
Sustainability 10 01433 g007
Figure 8. Geographical distribution of unit G1–G6 and region R1–R5. R1 represents the city center; R1–R4 represent the densely-populated regions; and R5 represents the sparsely-populated region.
Figure 8. Geographical distribution of unit G1–G6 and region R1–R5. R1 represents the city center; R1–R4 represent the densely-populated regions; and R5 represents the sparsely-populated region.
Sustainability 10 01433 g008
Figure 9. Hourly generation results for G1–G6 in Cases 1–4.
Figure 9. Hourly generation results for G1–G6 in Cases 1–4.
Sustainability 10 01433 g009
Figure 10. Hourly system emissions in Cases 1–4.
Figure 10. Hourly system emissions in Cases 1–4.
Sustainability 10 01433 g010
Figure 11. Distribution of API increments (ΔAPI) in Cases 1–4. Four representative hours of the day are chosen: 6:00, 12:00, 18:00, and 24:00. ΔAPI refers to the contribution of coal units on regional API. White dotted lines represent the five regions as shown in Figure 8.
Figure 11. Distribution of API increments (ΔAPI) in Cases 1–4. Four representative hours of the day are chosen: 6:00, 12:00, 18:00, and 24:00. ΔAPI refers to the contribution of coal units on regional API. White dotted lines represent the five regions as shown in Figure 8.
Sustainability 10 01433 g011
Figure 12. (a) Cav and (b) Cpk of SO2 in R1 under the Eco-SCUC with respect to differentiated APIE,tj, with and without wind power.
Figure 12. (a) Cav and (b) Cpk of SO2 in R1 under the Eco-SCUC with respect to differentiated APIE,tj, with and without wind power.
Sustainability 10 01433 g012
Table 1. Comparisons of total cost, total fuel, total emission, and SO2 concentrations at R1–R5.
Table 1. Comparisons of total cost, total fuel, total emission, and SO2 concentrations at R1–R5.
CaseCost (103$)Fuel (ton)Emission (ton)Cav (μg/m3)Cpk (μg/m3)
R1R2R3R4R5R1R2R3R4R5
1595.9310138.055.460.257.652.941.338.170.768.265.153.249.0
2512.738683.850.559.955.750.440.136.568.863.957.348.445.4
3523.468957.248.555.551.247.138.231.266.460.153.345.143.1
4527.629174.253.250.247.043.335.339.857.954.348.441.449.6
Table 2. Comparisons of total generation and daily average generation-associated API in Cases 3 and 4.
Table 2. Comparisons of total generation and daily average generation-associated API in Cases 3 and 4.
UnitTraditional SCUC (Case 3)Eco-SCUC (Case 4)
Generation (MW)ΔAPIi,av 1Generation (MW)ΔAPIi,av 1
R1R2R3R4R5R1R2R3R4R5
G14763.70.0%0.3%0.0%0.0%0.7%5572.50.0%0.5%0.0%0.0%3.4%
G25334.90.5%0.0%0.0%5.6%0.9%6015.80.8%0.0%0.0%6.8%5.1%
G35518.51.1%15.6%12.3%0.0%1.1%5044.20.8%9.4%9.1%0.0%0.9%
G44239.10.6%0.2%0.0%0.4%0.6%4763.90.9%0.3%0.0%0.6%1.2%
G53678.431.4%14.4%13.5%9.8%0.5%1622.115.9%7.1%7.2%1.5%0.4%
G60.00.0%0.0%0.0%0.0%0.0%516.10.0%0.0%0.0%0.0%7.1%
Total23534.633.6%30.5%25.8%15.8%3.8%23534.618.4%17.3%16.3%8.9%18.1%
1 ΔAPIi,av represents daily average value of regional API increment caused by coal-fired unit.
Table 3. Comparisons of pollutant saturated hours and sum of person-hour exposed to pollutant saturation in an annual case.
Table 3. Comparisons of pollutant saturated hours and sum of person-hour exposed to pollutant saturation in an annual case.
CaseTs (Hour)P × Ts (109 Person·Hour)
R1R2R3R4R5Sum of R1–R5
11101733411194594.078
21056685398181553.754
31063719379178573.894
4561437316153632.204

Share and Cite

MDPI and ACS Style

Guo, D.; Yu, J.; Ban, M. Security-Constrained Unit Commitment Considering Differentiated Regional Air Pollutant Intensity. Sustainability 2018, 10, 1433. https://doi.org/10.3390/su10051433

AMA Style

Guo D, Yu J, Ban M. Security-Constrained Unit Commitment Considering Differentiated Regional Air Pollutant Intensity. Sustainability. 2018; 10(5):1433. https://doi.org/10.3390/su10051433

Chicago/Turabian Style

Guo, Danyang, Jilai Yu, and Mingfei Ban. 2018. "Security-Constrained Unit Commitment Considering Differentiated Regional Air Pollutant Intensity" Sustainability 10, no. 5: 1433. https://doi.org/10.3390/su10051433

APA Style

Guo, D., Yu, J., & Ban, M. (2018). Security-Constrained Unit Commitment Considering Differentiated Regional Air Pollutant Intensity. Sustainability, 10(5), 1433. https://doi.org/10.3390/su10051433

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop