1. Introduction
Currently, China faces dual pressures from air pollution and greenhouse gas emissions. China’s regional air pollution problems are very serious: More than one-fifth of the country’s land has experienced extraordinarily high levels on the “smog index”. A “top-down”, mandatory commitment to the carbon emission reduction model dictated by the Kyoto Protocol is now gradually being transformed into a “bottom-up”, “Nationally Determined Contributions” (NDC) mechanism of the Paris Agreement [
1]. China expects the future of carbon emissions per unit of GDP to fall by 60–65% by 2030 compared to 2005. Due to the severity of domestic air pollution, combined with the commitment to join the international community in implementing an independent policy to reduce greenhouse gas emissions, both the “13th Five-Year Plan for Controlling Greenhouse Gas Emissions” and the “13th Five-Year Plan for Ecological Environmental Protection” issued by the State Council mention “enhanced coordinated control of carbon emissions and atmospheric pollutant emissions”. The collaborative governance of air pollutants and greenhouse gases is an excellent choice and policy outlet for China during a critical stage in its economic transformation.
According to “China statistical yearbook on environment 2017”, the total investment in treatment of environmental pollution reached 921.98 billion yuan in 2016. Among the three main areas of environmental pollution—air, water, and soil—the current investment related to air pollution is the highest. For example, in the industrial sector, 68.5% of the total investment in pollution treatment was spent on the treatment of waste gas (SO
2, NO
X, smoke, dust, etc.). In the short term, the cost of greenhouse gas reductions is derived from potential GDP loss caused by production cutbacks [
2]. Irrespective of both the financial constraints and the needs of economic development, only focusing on reducing emissions will not prove conducive to the long-term development of collaborative governance [
3,
4]. Therefore, it is necessary to consider the costs paid to reduce air pollutants and greenhouse gas emissions, namely abatement costs [
5]. This paper proposes a cooperative optimization-based approach to reduce emissions and abatement costs. Its underlying goal is to reduce both air pollutant and greenhouse gas emissions and cut down abatement costs simultaneously in order to achieve a coordinated balance.
While previous studies have been devoted to the collaborative governance of air pollutants and greenhouse gases in the context of technological improvements and macro control, in the long run, the energy structure, which is the ratio between one energy type consumption and the total energy consumption [
6], holds the key to achieving collaborative governance [
7]. The main air pollutants of China such as SO
2, and NO
X, and greenhouse gas CO
2 are produced from the consumption of fossil fuel-based energy sources, including coal, oil, and natural gas, which also offer the potential for collaborative governance [
8]. When equal amounts of coal, oil, and natural gas are consumed, the emissions of SO
2, NO
X, and CO
2 vary greatly. For example, 1 kg of oil is consumed to produce 5.77–36.25 grams of NO
x, while equal amounts of coal and natural gas yield 1.88–8 grams and 0.736–2.085 grams of nitrogen oxides, respectively. The difference in the method of use, even for the same amount of energy, will affect the emission ratio of SO
2, NO
X and CO
2. For example, due to advanced desulfurization facilities in industrial sectors, the amount of SO
2 produced by coal combustion is lower than that of other sectors in society. However, higher temperatures result in a larger amount [
9].
According to the differences in both the energy type and the combustion mechanism, the collaborative optimization method proposed in this research primarily concerns the arrangement of different energy sources (coal, oil, natural gas) in different sectors without disrupting production. This is an attempt to minimize the emissions of air pollutants and greenhouse gasses, while also reducing abatement costs and achieving both environmental and economic benefits. By taking Tianjin as an example, the effectiveness of the collaborative optimization method is also verified. The main contributions are as follows:
- (i)
The model for collaborative optimization of emissions and abatement costs is established. On the basis of analyzing the common constraints of different modes of energy consumption, the optimization of multiple targets, including reductions in air pollutant and greenhouse gas emissions, as well as abatement cost reductions, is realized by adopting a multi-objective optimization method.
- (ii)
A perspective of energy structure is proposed in order to solve the problem of collaborative optimization. According to the type and combustion mechanism of a given form of energy, limited energy (coal, oil, natural gas) is arranged reasonably into different social sectors, so as to achieve collaborative optimization on the premise of meeting social needs.
- (iii)
The collaborative optimization scheme of Tianjin is proposed. The relevant parameters of the collaborative optimization model are determined using the actual data, such as the yearbook of Tianjin. The collaborative optimization scheme, as well as the specific energy structure adjustment method in Tianjin, is also given in order to realize a reduction in both emissions and abatement costs.
The rest of this paper is structured as follows:
Section 2 consists of a literature review,
Section 3 introduces a collaborative optimization model, while
Section 4 offers empirical analysis of the collaborative optimization model in Tianjin.
Section 5 discusses the research results, and
Section 6 concludes the paper.
2. Literature Review
The problem of coordinated governance in the context of air pollutants and greenhouse gases has attracted the attention of scholars both at home and abroad. Some researchers have used LEAP, TRACE-P EI, CMAQ, and BenMAP models, as well as empirical research methods, to verify the synergistic effects of air pollutant and greenhouse gas emission reduction [
1,
10,
11]. The following research pays attention to the methods of collaborative reduction from a number of perspectives. For example, a self-management approach determined by a committee was established to facilitate the control of the total quantity of greenhouse gas and air pollutant emissions [
12]; structural emission reduction measures are easier to achieve with the coordinated control [
13]; factors such as discount rate and weight will affect the emissions reduction supply curve [
3]. A few papers have also focused on the correlation between energy consumption, and emissions of air pollutants and greenhouse gases. Meng et al. (2016) used multi-objective means to discuss collaborative emission reductions according to the common constraints of different forms of energy consumption, while also formulating a cleaner energy distribution plan [
9]. Yuan (2017), meanwhile, provided a comprehensive analytical framework in order to quantitatively evaluate the effects of energy measures in mitigating climate change and combating air pollution [
14]. Other studies have focused on the abatement costs of air pollutants and greenhouse gases. For example, Chae (2010) conducted cost effectiveness and synergy evaluations of air improvement and greenhouse gas control measures in Seoul, South Korea [
15], while Yang et al. (2013) estimated the environmental damage costs of air pollutants in various provinces and cities in China [
16]. Liu et al. (2019) have also conducted a cost–benefit analysis of measures taken by the transport sector in China’s Pearl River delta region to improve air quality and mitigate climate change [
4]. The above literature either discusses methods for the collaborative reduction of air pollutants and greenhouse gases from a technical and management level or analyzes abatement costs from an economic perspective. However, both emissions and abatement costs form important components of collaborative governance in China [
17]. The singular pursuit of emission reduction or abatement cost reduction is not conducive to supporting both environmental quality and economic development. In addition, due to the homologous nature of air pollutants and greenhouse gases, optimizing the energy structure is the key to achieve collaborative governance in the long-term [
7,
18]. However, only a few articles have discussed emission reduction schemes from the perspective of energy consumption, and these studies have not also considered the abatement costs.
Some scholars have combined energy structure with air quality and carbon emissions to conduct the research. Bilgen (2014) studied the environmental impacts of fuel types, industry, and global energy consumption, finding that the energy consumption structure and emission conditions of SO
2, NO
x, and CO
2 had a significant impact on acid rain and the greenhouse effect [
19]. Zhang et al. (2015) considered energy pollution reduction and other indicators to carry out a critical assessment of China’s existing energy model [
20], while Xing et al. (2017) examined various social factors that affect the emissions of greenhouse gases and air pollutants and emphasized the importance of energy transformation [
7]. Yang et al. (2018) adopted an LMDI model to quantify carbon emissions resulting from energy consumption in Wuhan, indicating that energy structure optimization contributed significantly towards the inhibition of carbon emissions growth [
8]. Existing research has investigated the relationship between energy structure, environmental quality, economic development, and the effectiveness of energy structure at a macro level. However, there has not been an in-depth systematic analysis of energy structures that can simultaneously ensure both environmental quality and economic development. In addition, most previous research results demonstrate that the proportion of new renewable energy schemes should be increased in order to optimize the energy structure. However, there are no renewable energy sources to completely replace fossil fuels at this time. Consequently, it is of practical significance for China to optimize its fossil fuel energy structure according to the pollution characteristics of these forms of energy in terms of the abatement costs related to air pollutants and greenhouse gases.
In recent years, multi-objective methods have been increasingly applied to social problems, such as environmental improvement and energy structure adjustment. Meng et al. (2016) used a multi-objective method to establish a collaborative emission reduction model for greenhouse gases and air pollutants, and also designed the optimal-pole algorithm subsequently [
20]. Mohamadpour and Hassanzadeh (2018) took reverse logistics as their research object, adopting a multi-objective mixed-integer linear planning model to alleviate environmental problems and solved it using the distance method and the ε-Constraint method [
21]. Li et al. (2019) proposed a new air quality index (AQI) analysis and prediction system, which uses a multi-objective multi-verse optimization algorithm to predict the hourly AQI series and to overcome randomness and non-stationarity [
22]. As the multi-objective function exists within a certain decision space, a solution that can combine multiple objectives and achieving optimality is urgently needed. As this kind of overall consideration is more suitable for the coordinated governance of air pollutants and greenhouse gases, this paper uses the mathematical method offered by the multi-objective approach to research.
This paper focuses on the collaborative optimization of emissions and abatement costs. The originality lies in the following two points: (1) considering emissions and abatement costs for air pollutants (SO2 and NOX) and greenhouse gases (CO2) simultaneously, a multi-objective optimization method is adopted to explore a clean and economically beneficial collaborative optimization scheme so as to achieve unification between environmental benefits and economic benefits; (2) from the perspective of energy structure, using the complementarity of coal, oil, and natural gas in terms of energy consumption, according to their combustion mechanism and pollution production characteristics, limited energy is reasonably arranged in different social sectors in order to achieve collaborative optimization without disrupting production.
3. Collaborative Optimization Model
We begin this section by determining the functional relationship between energy consumption and the emissions related to air pollutants and greenhouse gases. Following this, the abatement cost function is constructed in order to analyze the relationship between energy consumption and abatement costs. On this basis, the multi-objective method is adopted to establish a collaborative optimization model for emissions and abatement costs. Finally, in order to obtain solutions that are close to ideal values, the distance method is introduced to solve the problems related to a multi-objective approach.
3.1. Functional Relationship between Energy Consumption and Emissions
The energy in this paper refers to fossil fuels that release air pollutants (SO2 and NOX) or greenhouse gases (CO2).
According to the air pollutant and greenhouse gas emission prediction method proposed in “Resources (energy), Environment, Economic Prediction and Research Report of the National 13th Five-Year Plan”, society can be divided into multiple energy-consuming sectors. Following this, the quantitative relation between energy consumption and the emissions of air pollutants and greenhouse gases can be determined based on desulphurization efficiency, the denitration efficiency of different sectors, the carbon emission factors of different forms of energy, combustion loss rate, and the NOX emission factors of different energy in different sectors.
The SO
2 emission function is shown in Equation (1), where
is the amount of the
th energy consumed by the
th sector,
is the sulfur content per unit of the
th energy,
is the sulfur conversion rate of the
th energy, and
refers to the desulphurization efficiency of the
th sector.
The NO
X emission function is shown in Equation (2), where
is the NO
X emission factor of the
th sector of the
th energy consumption, and
refers to the denitration efficiency.
The CO
2 emission function is shown in Equation (3), where
refers to the carbon emission factor of the
th energy, and
refers to the combustion loss rate of the
th energy.
According to the above analysis, although the production mechanisms of SO2, NOX, and CO2 are different, they all depend on the energy consumed by the society and are all positively linearly correlated with energy consumption. In addition, due to differences in the technological means, the methods of use, and the amount of energy used by various sectors of society, the emissions of air pollutants and greenhouse gases also vary considerably. Even different sectors that use the same form of fuel will produce different emission amounts. For example, for NOx, this gas can be formed not only during combustion of the fuel, but also at the high temperature to oxidize nitrogen molecules in the air.
3.2. Functional Relationship between Energy Consumption and Abatement Costs
It can be seen from Equations (1) and (2) that exhaust gas with SO2 and NOX must be de-sulfurized and de-nitrated before it can be discharged into the air. Furthermore, the cost of performing these operations is classified as abatement costs; the CO2 emissions cannot be disposed of freely under constraints relating to known viable technologies. Carbon emission reduction, in the short term, can therefore only be achieved by reducing the production amount at the source of the emissions, thereby slowing economic growth. The potential loss of output is classified as the abatement costs.
3.2.1. Functional Relationship between Energy Consumption and Abatement Costs for Air Pollutants
It is generally believed that the abatement costs are related to the amount of air pollutants emissions and removals (Zhou et al., 2018) [
23]. Some studies suggest that abatement costs are also correlated with regional economic conditions (Tang & Chen, 2017) [
24]. The abatement cost function of air pollutants is usually obtained through regression analysis.
This paper comprehensively considers the amount of air pollutants emissions and removals, while also describing economic indicators in the target area in order to conduct the regression calculation of the cost of air pollutant governance in each region. The SO
2 abatement cost function is shown in Equation (4), and Equation (5) is obtained by logarithmic processing of Equation (4). Following this, linear regression will be performed by software in order to obtain the various parameter values.
where
is the abatement cost function of SO
2,
represents the amount of SO
2 removed, and
are the constants.
is a combined impact factor, including regional industrial structure, enterprise ownership structure, etc., which can generally be considered unchanged in the short term; and
are the elastic coefficients of emissions, removals and economic indicators respectively.
The parameter is normalized by a series of indicators that characterize the regional economic development level. is obtained through principal component analysis, where is the coefficient, represents the normalized value of the economic indicator, and .
Similarly, the abatement cost function of NO
X can be obtained by Equation (6), where
is the abatement cost function of NO
X,
is the amount of NO
X removed, and
are the constants.
The amount of air pollutants removed is defined as the difference between those produced and those emitted. That is, the product of the amount of air pollutants produced, and the rate of air pollutants removed. The quantitative relationship between energy consumption and the abatement costs of SO
2 and NO
X is obtained by substituting Equations (1) and (2) into Equations (4) and (6), as shown in Equations (7) and (8) respectively.
3.2.2. Functional Relationship between Energy Consumption and Abatement Costs for Greenhouse Gases
Some researchers have conducted a comprehensive analysis of the cost of greenhouse gas governance, concluding that the CO
2 marginal abatement costs of different regions, different energy structures, and industrial structures are quite different. In terms of the technical and structural constraints, carbon emissions cannot be disposed of freely. To achieve the short-term goal, it can only be accomplished by compressing the production scale of the carbon emission sector and slowing economic growth. Therefore, the governance of CO
2 comes at the cost of potential GDP losses [
25].
Wu et al. (2019) proposed an improved DDF-DEA method to effectively avoid the occurrence of marginal emission reduction of non-positive CO
2 and calculate the provincial CO
2 marginal abatement cost [
2]. Our analytical framework is based on this model. The greenhouse gas abatement cost function is shown in Equation (9), where
is the abatement costs of CO
2,
is the marginal abatement costs of reducing CO
2 by
units, and
is the amount of CO
2 removed.
Taking each year as a time scale, the marginal abatement cost of the year can be used as the unit cost of CO2 for the year. Since there is no centralized emission reduction measure for CO2, that is, where the amount of CO2 production is equal to its emissions, the amount of CO2 removed is regarded as the difference between the actual amount of CO2 emissions and that of target emissions in the th year. As this paper aims to reduce the amount of CO2 emissions, set .
The equation for the greenhouse gas abatement cost function of the
th year is as follows:
Both and are functions of , as shown in Equation (3).
3.3. Multi-Objective Model for Collaborative Optimization
Society can be divided into
sectors with
optional energy sources. The decision variable
is the
th energy consumption of the
th sector; and
are the emissions of SO
2, NO
X, and CO
2, respectively; while
denotes total abatement costs; and
,
denotes the abatement costs of SO
2, NO
X, and CO
2. The minimum value is adopted by the objective function, as shown below.
The energy supply must meet the social development needs, but it is also subject to some conditions, such as reserves and exploitation capacity. Therefore, energy demand constraints and supply constraints can be expressed as in Equations (11) and (12), where
refers to the coefficient of the
th energy converted into standard energy;
refers to the total energy demanded to develop the
th sector; and
refers to the upper limit of the
th energy supply.
In addition, energy consumption must be positive, i.e., .
In summary, the obtained collaborative optimization model is shown in Equation (13).
3.4. Distance Method and Solving Process
Using the distance method to solve the multi-objective problem is more conducive towards obtaining a solution that is close to the ideal value (Branke et al., 2008) [
26]. First, each objective function is solved independently under given constraints (Mirzapour Al-E-Hashem et al., 2011) [
27] where the minimum value of SO
2, NO
X, and CO
2 emissions are
,
, and
respectively, while the minimum cost of governance is
. The second stage involves taking the ideal values
,
,
, and
of the four objective functions of the target space and minimizing the Euclidean distance
between the four targets and their ideal values under the established constraints set by the model. The solution of the original multi-objective function can be obtained by solving the nonlinear programming problem in Equation (14).
Where
,
,
,
represents the relative weight of the four goals regarding SO
2 emissions, NO
x emissions, CO
2 emissions, and abatement costs respectively, while
.
4. Empirical Analysis of Collaborative Optimization Model in Tianjin
There are two reasons for choosing Tianjin as the sample of this empirical study. First, Tianjin is an important economic center and a major city in northern China. However, it is also one of the country’s most polluted areas where economic development and pollution control have become severe issues. Second, the energy consumption pattern of Tianjin is dominated by fossil fuels, which is typically representative of most industrial cities across China and, therefore, representative as an example for empirical analysis. This paper employs data from Tianjin from 2006 to 2016 to verify the model. Due to the lag in the publication of statistical yearbooks, the data for 2016 is the latest available data for this study.
4.1. Determination of Parameters in Emissions Function
Drawing on the work of Meng et al. (2016) [
7], in order to reduce the number of calculations according to the amount of energy consumption, Tianjin can be divided into five social sectors: large-scale agriculture (farming, forestry, herding, fishery, and water conservation), industry (including power generation and heating), transportation, retail and accommodation, and consumption of living. At the same time, based on the similarity of the pollutant discharge coefficient in energy consumption, the various forms of fossil fuel energy included in the energy yearbook are classified as coal, oil, and natural gas. Therefore, a collaborative optimization model with 15 variables is established.
As the SO
2 and NO
X generated in daily life are difficult to process centrally, the desulfurization and denitrification procedures considered in this paper mainly involve the industrial field. The desulfurization rates of coal, oil and, gas in Tianjin industrial departments are 0.3, 0.3, and 0 respectively, while the denitrification rates are 0.196, 0.08, and 0, respectively. According to both “Resources (energy), Environment, Economic Prediction and Research Report of The National 13th Five-Year Plan” and “China energy yearbook”, the parameters in the emission function, namely the discharge coefficients of various energies, are shown in
Table 1.
4.2. Determination of Parameters in Abatement Cost Function
4.2.1. Determination of Parameters in Air Pollutants Abatement Cost Function
In determining the abatement cost function of air pollutant in Tianjin, this paper makes the following assumptions and simplifications:
- (i)
The data of air pollutant emissions and removals in
Table A1 refers to those of the industrial sector. As the SO
2 and NO
X generated by other sectors are difficult to be processed centrally, only the industrial SO
2 and NO
X governance data are included in the statistical yearbook.
- (ii)
The data of air pollutant abatement costs in
Table A1 refers to the “Operating Cost of Air Pollutants Governance Facilities” in the “China Environmental Yearbook”, including energy consumption, equipment depreciation, equipment maintenance, personnel salaries, management expenses, process chemical costs, and other expenses related to the operation of the facility.
The “China Environment Yearbook” was used so as to obtain data about air pollutant governance in Tianjin on the basis of the above assumptions and simplifications. The indicators of the on-the-job worker’s average wage, per capita GDP, fixed assets investment, and industrial added value of employees were used to characterize Tianjin’s macroeconomic and industrial development level, while the “Tianjin Statistical Yearbook” was adopted to obtain Tianjin macroeconomic data (as listed in
Table A1).
The economic parameter
of each year was generated using principal component analysis. Due to the negative number relating to this, in order to facilitate the next dynamic measurement analysis, according to the principle
of statistics, the formula
was used to coordinate the translation to eliminate the negative influence, while the final economic parameters
of each year were also obtained (as listed in
Table A2).
Following this, the abnormal data was interpolated, while SPSS software was used to execute regression analysis on the variables involved in the SO
2 abatement cost function. The results show that the linear relationship between the dependent variable and all independent variables is not significant (as is described in
Table A3). It was also found that there were multiple collinearity problems, in which the tolerance of
was small and the VIF value was large. Especially when the data of 2006–2010 was processed separately, the tolerance was 0.067 and the VIF value reached 13.08, which means that multiple collinearity is more severe. Furthermore, through the Sobel test, it was found that the economic parameters play a mediating role between the amount of SO
2 removals and abatement costs (Sig < 0.05), while also playing a mediating role between that of SO
2 emissions and abatement costs (Sig < 0.05). Due to the existence of multiple collinearity and mediating effects, the economic parameter was removed from the abatement cost function.
The function of the SO
2 abatement costs for Tianjin was regressively determined (as outlined in
Table A4). For the F-test, with sig. 0.016, the linear relationship between the dependent variable and all independent variables was significant, producing a tolerance of 0.962 and a VIF of only 1.039. There was no collinearity problem in this instance. After removing the logarithm, according to Equation (3), the abatement cost function of SO
2 in Tianjin can be obtained in the following:
Similarly, after the abnormal data was interpolated, SPSS software was used to analyze the variables involved in the NOX abatement cost function. It was also found that the regression coefficient of to was not significant and that there is also a collinearity problem. The economic parameter was removed from the abatement cost function.
The function of the NO
2 abatement costs for Tianjin was regressively determined (as listed in
Table A5). For the F-test, with sig. 0.015, the linear relationship between the dependent variable and all independent variables was found to be significant, with a tolerance of 0.777, while for the VIF it was only 1.287. There is no collinearity problem. After removing the logarithm, according to Equation (3), the abatement cost function of NO
X in Tianjin can be obtained in the following:
4.2.2. Determination of Parameters in Greenhouse Gas Abatement Cost Function
This paper draws on Wu et al. (2019) to improve the DDF-DEA method so as to calculate the provincial marginal abatement cost of CO
2 [
2], where the marginal abatement cost of CO
2 in Tianjin in 2016 was 2208.83 yuan/t. After substituting into Equation (10), the abatement cost function of CO
2 in Tianjin in 2016 is as follows:
can be calculated by Equation (3) using the data in
Table 1 and
Table 2,
is function of
, as shown in Equation (3).
4.3. Determination of Constraint Parameters in a Collaborative Optimization Model
Corresponding to model 3.13, in terms of the demand constraint, is the coefficient of energy converted into standard coal, while the corresponding coefficients of coal, oil, and natural gas were 0.71, 1.4, and 14.3 respectively. constitutes the sector’s energy demand, which is expressed in this paper using the standard coal value of the th sector.
In the supply constraint, is the supply of energy. Coal, oil, and natural gas cannot be supplied indefinitely. However, if the energy supply is set too low, this will result in less space available for energy structure adjustment. Conversely, if the setting is too high, coal, and oil will be substituted by natural gas, which is inconsistent with the realities of natural gas supply. Therefore, this paper set the energy supply at 150% of the total energy consumption in Tianjin in the th year.
According to the terminal consumption of coal, oil, and natural gas by social sectors of Tianjin in 2016, as well as the coefficients of various forms of energy converted into standard coal, the final consumption amounts of Tianjin in 2016 are shown in
Table 2.
5. Results and Discussion
With the model parameters and distance method above, MATLAB can be used to solve the final optimization problem (including 8 constraints and 15 decision variables). Referring to the practice in Reference [
28], this section outlines the final results and discusses the importance of different objectives in collaborative optimization, how to adjust energy consumption in different sectors of society to achieve clean and economic collaborative optimization, and the impact of the addition of abatement cost objectives on the optimization scheme.
5.1. Optimization Scheme and Analysis
According to the above method, the emissions of SO
2, NO
X, and CO
2 were 233,000 tons, 261,000 tons, and 700,000 tons, respectively, while the aggregate abatement costs were 15.48 billion yuan (the specific optimization scheme is listed in
Table 3). In 2016, the actual emission amounts of SO
2, NO
X, and CO
2 in Tianjin were 260,000 tons, 294,000 tons, and 729,000 tons, respectively, while the abatement costs was found to be 16.10 billion yuan, indicating that the projected environmental and economic benefits obtained by the optimization scheme were much better than those of the current situation (as listed in
Table 3).
In terms of the importance of each target, the weight values of the four targets in the optimization scheme were 0.21, 0.11, 0.33, and 0.35, respectively. In order to achieve collaborative governance, as well as to improve environmental and economic benefits, it is necessary to focus on abatement costs (especially CO2 abatement costs) and CO2 emissions. On the one hand, Tianjin is not only the largest open coastal city in China, but also one of the country’s most important industrial cities. The large amount of energy consumption within the city means that its greenhouse gas emissions are high. Tianjin is the only municipality to have participated in low-carbon provinces and cities, greenhouse gas emission lists, and regional carbon emissions pilot, which all involve the challenge of reducing CO2 emissions. On the other hand, through upgrading industrial infrastructure and promoting the development of high-tech industries, Tianjin has made great contributions to reducing carbon emissions. This region has its reduced emission range and has a relatively high MAC compared to undeveloped regions. According to the calculations of Wu et al., in 2019, the Beijing–Tianjin–Hebei region had the highest shadow price of CO2 in China, making CO2 emissions and its abatement costs a priority.
In terms of the governance of air pollutants, the weight of SO2 (0.21) was higher than that of NOX (0.11), while SO2 should constitute the main control object of air pollutants. Tianjin is one of the SO2 pollution control areas regulated by the state, while the reduction emission of SO2 began at the start of this century and has since been rolled out fully. The government has taken active measures to promote the use of clean energy and eliminate backward enterprises, while also developing and promoting desulfurization technologies and other specific measures as well. However, due to the development and positioning of Tianjin as an industrial city, it is also necessary to continue to rely largely on fossil fuel-derived energy in order to maintain its economic development level. The energy characteristics of China—abundant coal, limited oil, and scarce gas—make coal an important commodity in Tianjin’s energy structure. The SO2 released by coal combustion accounts for over 90% of the region’s total SO2 emissions. Therefore, Tianjin SO2 has not yet been fundamentally governed.
Using the corresponding energy allocation scheme (as listed in
Table 4), the consumption ratio of coal, oil, and gas (converted to standard coal) was adjusted from the original 17.5%: 65.5%: 17.0% respectively, to 15.8%: 58.4%: 25.8%. In other words, the law that gas should be used instead of oil and coal should be followed in order to cope with the contradiction between the strong growth of energy consumption and the need to reduce environmental damage. Comparing the actual energy consumption data and optimization schemes of various social sectors, we can say that the optimized amount of natural gas used is equal to that of its supply. At this time, the energy allocation scheme is concentrated on increasing the proportion of natural gas in the industrial sector, while the amount of natural gas assigned to the industrial sector is dominant in the overall supply of natural gas. The consumption ratio of coal, oil, and natural gas (converted to standard coal) was also adjusted from the original 21.3%: 63.2%: 15.5% to 13.5%: 49.6%: 37.0%. In the energy structure, the proportion of natural gas used has been greatly increased, meaning that the main source of air pollutants in the industrial sector-SO
2 and NO
X, have been greatly reduced.
5.2. Impact of Abatement Costs Goals
On the basis of previous studies which have focused on collaborative reduction of air pollutants and greenhouse gases from a technical and management level, or unilaterally analyzing abatement costs from an economic perspective, this paper analyzes both emissions and abatement costs to ensure successful outcomes to collaborative governance. In order to analyze the impact of the addition of abatement costs goals on emission reduction effects and energy structure optimization with sufficient depth, the following comparisons are made to analyze the optimization schemes considering both emissions and governance costs and the schemes that only considering emissions.
If only environmental benefits are involved, the emissions of SO
2, NO
X, and CO
2 will be the main control targets, while the weight of abatement costs will be
. Using the same method, the emissions of SO
2, NO
X, and CO
2 are 220,000 tons, 175,000 tons, and 672,000 tons, respectively, compared with the actual emissions of Tianjin in 2016 (260,000 tons, 294,000 tons, 729,000 tons). At the same time, considering the environmental and economic benefits, gas emissions (233,000 tons, 261,000 tons, 700,000 tons, respectively) are lower, indicating that the emissions of SO
2, NO
X, and CO
2 in the environmental benefit optimization scheme are lower. Furthermore, the addition of economic benefits also limits the improvement level of environmental benefits. In other words, some environmental benefits must be sacrificed in order to improve economic efficiency (as described in
Table 3).
In terms of the importance of each target, the weights of the three objectives in the environmental benefit optimization scheme are 0.32, 0.17, and 0.51, respectively. In order to improve environmental benefits, the reduction of greenhouse gas emissions should be prioritized, while SO2 should be the main control object of air pollutants. Its weight should also be higher than that of NOX. This observation is similar to the conclusion when considering both environmental and economic benefits, which is also based on the unique circumstances of Tianjin’s urban development foundation and positioning.
In terms of energy allocation (as outlined in
Table 4), compared with the actual consumption data and environmental benefit optimization scheme, the consumption ratio of coal, oil, and gas (converted to standard coal) was adjusted from the original 17.5%: 65.5%: 17.0% to 0.1%: 74.3%: 25.6%. Put simply, natural gas is still the most desired form of fossil fuel-based energy, while coal should be almost completely replaced by oil and natural gas. When considering the environmental and economic benefits simultaneously, the consumption ratio of the three is 15.8%: 58.4%: 25.8%. In other words, in terms of balancing the economic benefits, it is necessary to retain part of the consumption ratio of coal. In view of the current situation in Tianjin, in order to reduce abatement costs, the status of coal within the overall energy structure cannot be completely replaced, which is more consistent with the reality.
At the same time, the environmental benefit optimization program still needs to meet the development needs. However, compared to the optimization scheme that involves both environmental and economic benefits, the energy allocation optimization scheme focuses on the use ratio of natural gas in the transportation sector, in which the natural gas use ratio accounts for 99.4% of its total energy consumption and the fuel consumption ratio is decreased from 91.2% to 0; at the same time, the proportion of consumption of living of coal dropped to almost zero, while the proportion of natural gas increased from 18.9% to 65.8%. This is in line with the fact that natural gas consumption has witnessed explosive growth in recent years, while it is again stated that the emission reduction targets concerning air pollutants are constrained by the supply of natural gas. Higher environmental benefits will be accomplished if the transportation sector’s project, “oil to gas” is promoted and the proportion of coal used in daily life is reduced. However, under the premise of limited costs concerning atmospheric control, it is better to satisfy the environmental and economic benefits of the industrial sector by prioritizing meeting the growing demand for natural gas.
6. Conclusions
Previous studies have typically discussed the collaborative emission or abatement costs for air pollutants and greenhouse gases unilaterally, or have instead explored the relationship between energy structure, environmental quality and economic development from the macro level and calculated the effectiveness of energy structure optimization. Considering the differences in the combustion mechanism and pollution production characteristics of fossil energy, such as coal, oil and natural gas, this research used a multi-objective method to construct a collaborative optimization model of the emissions and abatement costs for air pollutants and greenhouse gases. By rationally arranging the energy consumption of different sectors and adjusting the energy structure, the optimization of multiple objects, such as air pollutant emissions, greenhouse gas emissions, and abatement costs can be achieved. In order to verify the effectiveness of the proposed optimization method, this paper took Tianjin as an example, using yearbook data to determine the parameters in the emissions function, the abatement cost function, and the constraints. Finally, the Tianjin collaborative optimization model was found to be effective and the collaborative optimization scheme was also proposed.
Although the validity of the aforementioned model has been verified, there are still some factors that may affect its application in environmental management. Energy costs are often a key affecting factor. The energy characteristics of China—abundant coal, limited oil, and scarce gas—make oil and natural gas supply insufficient and costly relative to coal, which will affect the use of the model. Frequent energy price volatility exacerbates this effect. Energy conversion costs are another influencing factor. For example, in many cases, the switch from coal and oil to natural gas implies additional installation costs and requires the expansion of the pipeline network.
The proposal concerning "Collaborative Governance" in terms of air pollutants and greenhouse gases in developed countries is mainly used to evaluate the welfare effect of greenhouse gas emission reduction policies more comprehensively. China, as a developing country, was chosen due to the characteristics of its economic development and the economic structure of its energy sector. Collaborative governance, in order to realize air pollutants and greenhouse gas reduction emission, plays an essential role in preserving healthy living conditions and managing energy scarcity, while also continuing to succeed in the industrialization process. Atmospheric concerns, such as frequent extreme weather events and the destruction of our ecological balance, are in conflict with the process of industrialization unless processes such as collaborative governance evolve to counter their adverse effects. Collaborative governance based on China’s national interests is particularly important, while establishing a unified model framework and data specifications for initiating atmospheric collaborative governance in major regions of China is a topic worthy of further research.
A limitation of this research is only focusing on energy consumption, energy structure and abatement technology, but seldom involves state policy and strategy, and social and cultural factors due to the indirect effects of them on pollutant emissions. Therefore, further research could quantify these aspects in the described model. Another limitation is the algorithm. Although the distance method can effectively solve the multi-objective optimization problem, it also, to a certain extent, reduces the solution space. Furthermore, it also results in a strong dependence on the initial solution. Multi-objective intelligent algorithms (multi-objective genetics, multi-target particle swarms, etc.) and a number of other related methods can be applied to the future research of multi-objective collaborative optimization problems.