Next Article in Journal
Optimising Frequency-Based Railway Services with a Limited Fleet Endowment: An Energy-Efficient Perspective
Previous Article in Journal
Optimal Allocation of Distributed Generation Considering Protection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

What Does Cost Structure Have to Say about Thermal Plant Energy Efficiency? The Case from Angola

1
Department of Mathematics, Isfahan (Khorasgan) Branch, Islamic Azad University, Isfahan 81595-158, Iran
2
COPPEAD Graduate Business School; Federal University of Rio de Janeiro, Rua Paschoal Lemme, Rio de Janeiro 355. 21949-900, Brazil
*
Author to whom correspondence should be addressed.
Energies 2020, 13(9), 2404; https://doi.org/10.3390/en13092404
Submission received: 27 March 2020 / Revised: 4 May 2020 / Accepted: 7 May 2020 / Published: 11 May 2020

Abstract

:
This paper analyzes the efficiency of thermal power plants in Angola by means of a two-stage Data Envelopment Analysis (DEA) approach. In the first stage, a novel super-efficiency DEA model for undesirable outputs (CO2 emission levels and discharge of polluted water) is initially used to measure their efficiency levels. Then, in the second stage, relevant cost structure variables frequently used to describe a productive technology are employed as analytical thresholds for assessing energy production performance either in terms of capital or labor-intensity levels. Precisely, bootstrapped regression trees are used to discriminate super-efficiency scores yielding an energy production performance predictive model based on the technology type as proxied by its cost structure and their respective thresholds, since Angolan thermal plants are heterogeneous. Findings suggest that Angolan power plants are old and labor intensive, as some of them date back to the colonial era, and that lack of capital investment should be revised in favor of installing carbon capture devices. The approach developed here consists of a valuable approach for identifying priorities when technologically updating a heterogeneous thermal industry to face pollutant concerns.

Graphical Abstract

1. Introduction

This research focuses on a relatively understudied topic in thermal plants, which is the relationship between its productive efficiency and cost structure in a developing country such as Angola, a former Portuguese colony and actually an important oil producer and exporter in Africa. Although the efficiency of thermal power plants at the machine or equipment level is a well-defined and parameterized research stream based on the laws of thermodynamics ([1,2], here we are interested in assessing the efficiency of thermal energy production at the industry level where different plants consist of the units of analysis and represent a reflex of how a set of individual equipment interacts with other factors such as labor, managerial style, capital, etc.
Despite several previous studies focus on power plant efficiency in developed countries, only a handful of them have addressed this issue in developing nations, despite their growing relevance with respect to energy generation [3,4]. Scant previous research on thermal power plant energetic efficiency has employed distinct methods, but the most common rely on non-parametric approaches such as Data Envelopment Analysis (DEA). See for instance [5] and [6] for comprehensive reviews on the subject. This research analyzes the technical efficiency of energy production in thermal plants located in Angola using a new super-efficiency DEA model to account for undesirable outputs. The major underlying idea is to assess whether cost structure variables are related to technical super-efficiency scores in thermal energy production, thus providing a means for describing how different technological patterns can be employed to control for the emission of undesirable outputs (pollutants) such as CO2 emissions and discharged water. While scores are computed using a super-efficiency model, the cost structure variables are tested as efficiency thresholds in bootstrapped regression trees, thus allowing for the discrimination of different technological patterns in thermal energy production.
It is worth mentioning why the set of Angolan thermal power plants represents a controlled environment for conducting this search. First, these thermal power plants are all state owned and operate in a monopolistic market; therefore, the impact or influence of different managerial practices on efficiency scores and cost structures is quite limited. Additionally, these 32 plants were built over the course of more than forty years after Angola’s independence from Portugal and therefore different production technologies used in thermal energy generation and their evolution are reflected within this set of plants.
Another distinctive aspect of this research concerns the analysis of the cost structure as a means to identify technological patterns and propose improvement paths within the ambit of an industry or sector. In fact, there is an emerging research stream that can be observed in a few efficiency studies on energy and infrastructure areas that tries to correlate efficiency levels to the cost structure of the productive process. Wang et al. [5], Lin and Yang [6], and Barros and Wanke [7] analyzed possible improvement paths for energy production efficiency in the thermal power industry based on cost structure variables. Wanke and Barros [8] conducted a similar analysis within the ambit of airport operations. All of these studies focused on the emerging economies of Asia and Africa where issues regarding a better balance between capital intensity and labor expenditure are deemed relevant to achieve higher productivity levels.
This paper contributes to the literature body on energy efficiency planning in three distinct ways. Firstly, and for the first time ever, this research unveils the relationships among capital, labor, pollutant emissions, and energy efficiency with generation thermal power plants. While using Angolan data, the results presented here can be generalized to other contexts since capital and labor cost structure variables are the key fundamentals for energy generation performance. Secondly, this paper adds to this beginning research stream by applying a novel super-efficiency DEA model for undesirable outputs altogether with statistical learning techniques such as bootstrapped regression trees, thus allowing to explore the impact of different productive technologies by proxying them to the cost structure. In fact, this is the first time that DEA and machine learning are complementarily used to explore productive technologies by means of their underlying labor and capital cost structure. Thirdly, this paper contributes to a more general research stream beyond energy generation by digging into the cost structure variables and the respective outcomes of the productive process. While technological patterns are most discussed when different energy sources are compared (e.g., thermal vs. hydro), a deeper analysis within each type is deemed necessary. In this research, the term productive technology describes the technical means by which fuel, labor, and capital are combined to produce thermal energy. While the concept of productive technology transcends the modeling of a functional relationship between inputs and outputs, some of their intrinsic structural concepts could be captured and related to each other such as efficiency and cost structure for decision-making purposes.
This text is structured in five more sections. Section 2 presents the background of thermal power plants in Angola, then Section 3 presents the literature review focusing on the previous efficiency studies on thermal power plants using non-parametric methods such as DEA and its variants as cornerstones. Section 4 presents the new super-efficiency DEA model capable of handling undesirable or bad outputs developed for this research. Section 5 presents the dataset used and the bootstrapped regression trees method, which are followed in Section 6 by the results and discussion. Section 7 makes some concluding remarks.

2. Background on Angolan Thermal Power Plants

This search for different technological patterns in thermal plant energy efficiency and their underlying impact on pollutant emissions is conducted in 32 Angolan thermal plants during the period from 2010 to 2016. It is worth mentioning why this set of plants represents a controlled environment for conducting this search. First of all, these thermal power plants are all state-owned and operate in a monopolistic market, therefore the impact or influence of different managerial practices on efficiency scores and cost structures is quite limited. Besides, these 32 plants have been built over the course of more than forty years since the independence of this country from Portugal. Therefore, different productive technologies used in thermal energy generation and their evolution are reflected within this set of plants.
Precisely, Angolan thermal power plants are state-owned facilities that burn fuel—oil, coal, and natural gas—assuring energy supply for domestic use. The charge of generating electricity has been given to a state firm, ENE-EP (Empresa Nacional de Electricidade—Empresa Pública), which controls energy production in every province of Angola. As regards the planning of thermal energy production, the required fuel resources and other technical capacity issues of each plant are decided in a centralized fashion by ENE-EP. However, another state company, EDEL-EP (Empresa de Distribuição de Electricidade—Empresa Pública), controls electricity distribution in Angola. Specifically in the Luanda region, power distribution is controlled by a subsidiary of EDEL created specifically for this purpose. This arrangement harkens to colonial times when the Portuguese government performed infrastructure planning, which was highly concentrated in the Luanda region.
In fact, the technological choices embedded in the construction of the thermal power plant industry in Angola over the course of the years yielded a unique cost structure that reflects the type of fuel burnt, the age of the plant along with its relatively small size, and the managerial style imposed by the state control. Older technologies, small-scaled operations, and state control are elements (as corroborated by the literature review in Section 3) that are often related to inefficient operations.
It is reasonable to expect that all these elements are not only reflected in the efficiency levels, but also in the different cost structure variables and ratios that characterize the productive technology in thermal plants such as capacity cost, labor and capital costs, besides cost–asset and capital–labor ratios. All these variables are analyzed in this study and are further discussed and operationalized in Section 4. Thus, this study seeks to identify the subset of cost structure variables and ratios that best explain efficiency levels and pollutant emissions in Angolan thermal plants.

3. Literature Review on Thermal Power Plant Efficiency

Thermal power plant efficiency is an emerging research area. More than 20 papers have been published over the course of 18 years with most being published from 2010 on. Not only are these papers characterized by a diversity of methodological approaches, although the non-parametric related DEA methods prevail, but also their research questions are very diverse. They range from efficiency comparison between different technologies or the evolution of efficiency levels over the course of time to the impact of regulatory issues and ownership, besides emerging trends as regards carbon dioxide emissions and other sustainability issues.
For example, seminal papers in the area encompass Park and Lesourd [9] who measured the performance of 64 fuel power plants in South Korea. Lam and Shiu [10] computed the efficiency of China’s thermal power generation based on data for 1995 and 1996. Chien et al. [3] applied a DEA–Malmquist model to measure the efficiency of eight thermal power plants in Taiwan. Sarica and Or [11] analyzed and compared the performance of electricity generation plants in Turkey. Barros and Peypoch [12] analyzed the technical efficiency of Portuguese thermoelectric power generating plants with a two-stage procedure. Nakano and Managi [13] assessed efficiency in the Japanese steam-power sector and scrutinized the impact of reforms on the relative performance of its firms over the course of almost two decades. Sozen et al. [14] studied the performance of 11 lignite-fired, 1 coal-fired, and 3 natural gas-fired Turkish state companies using DEA. Liu et al. [15] assessed the efficiency levels of major thermal and combined cycle power plants in Taiwan for the period 2004–2006 using the DEA approach. Rezaee et al. [16] and Rezaee [17] presented a novel model based on game theory and DEA to measure efficiency of Iranian thermal power plants in light of diverse objectives. Shrivastava et al. [18] focused on the relative productivity of 60 Indian coal-fired plants by means of classic DEA models. Sueyoshi and Goto [19] used DEA to conduct an efficiency assessment of coal-fired plants in the US. Du and Mao [20] measured the environmental efficiency and costs related to carbon dioxide emissions in Chinese coal plants by means of a novel DEA model.
Recent papers still maintain this focus. Ghosh and Kathuria [21] investigated the impact of institutional quality typified as regulatory governance on the performance of thermal power plants in India. They estimated a translog stochastic frontier model using an index of state-level independent regulation as one of the determinants of inefficiency. Their findings show that technical efficiency is sensitive to both unbundling of state utilities and regulatory experience. Barros and Wanke [7] utilized a two-stage approach for efficiency evaluation in thermal power plants. First, a DEA Slacks Based Model (SBM) was employed to assess the relative efficiency of thermal power plants. In the second stage, beta regression models were combined with DEA-SBM efficiency scores to produce a model for predicting energy production performance. Yan et al. [22] evaluated the carbon emission efficiency using the Undesirable-SBM model and data from China’s power industry in 30 provinces from 2003 to 2014. They performed a spatial autocorrelation analysis that is based on Moran’s index to confirm the non-equilibrium spatial distribution of the carbon emission efficiency for the power industry. Wu et al. [23] employed an improved two-stage analysis model to analyze eco-efficiency of 58 Chinese coal-fired power plants. Firstly, the principal component analysis was selected for pre-treatment of variables in order to reduce dimensionality and distinguish prioritized factors. Secondly, the super-efficiency DEA was chosen to assess eco-efficiency with overall discriminatory rankings.
Table 1 presents a literature review synthesis of the previous studies on efficiency in thermal energy generation. These previous studies not only suggest that environmental issues are increasingly growing in importance over the years as more and more plants focus on reducing carbon dioxide emissions, but they also reveal that the production technology of each plant may be a relevant study field to understand efficiency in thermal energy generation. Precisely, plant age, fuel type, scale size, frontier shift, and catching-up effects appear to be the descriptors used the most in previous papers of the technological patterns that lie within thermal energy generation. Besides, ownership and regulatory impacts are also relevant issues addressed by some of these researches undertaken in different countries worldwide, especially in Asia and Europe, while research on African countries is still negligible. Besides, although the use of non-parametric methods such as DEA and SBM (Slacks Based Model) prevail, it is worth noting that the use of super-efficiency models is still scarce and focused on addressing environmental impacts (Wu et al., [23]).
Therefore, this research fills a literature gap by not only addressing the issue of thermal energy production efficiency in an important African country, but also by proposing a novel super-efficiency model to better discriminate efficiency scores in a sector where productivity variations are subtle, especially within the ambit of the Angolan state-owned plants, as discussed in Section 2. Additionally, the review synthesized in Table 1 also sheds some light, by contrast, on the nature of the contributions of this paper since most previous studies did not consider the impact of cost structure to diagnose technological patterns and their relationship with efficiency levels and pollutant emissions (Barros and Wanke, [7]). Some of them, however, use fuel prices to trace their impact on eco-efficiency levels (Wu et al., [23]). In fact, the joint use of super-efficiency DEA models and bootstrapped regression trees is an additional innovative feature of this research when compared to previous research described in Table 1 since for the first time statistical learning methods are employed in the second stage of analysis when efficiency scores are usually correlated or regressed onto contextual variables.
Hence, the distinctive aspect of the current research in comparison to the other previous researches is the joint use of novel super-efficiency models and statistical learning techniques to unveil the relationship between cost structure (as a proxy of the productive technology) and efficiency levels in thermal power plant energy generation. As explored in Section 5.2., statistical learning methods constitute a useful approach for unveiling hidden relationships between efficiency scores and cost structure variables, as long as they do not rely on the traditional parametric assumptions found in the typical regression approach (Tobit, truncated bootstrapped regression, beta, and Generalized Method of Moments (GMM), cf. Table 1), yielding higher analytical flexibility and explanatory power.

4. The Proposed Super-Efficiency DEA Model with Undesirable Outputs

A comprehensive review of previous studies shows that most DEA applications have considered primary cross sectional data and evaluated relative efficiencies in a single period, usually one year (Emrouznejad and Yang, [28]; Fernández et al., [29]). Exceptions are found in window analyses (Charnes et al., [30]) and other models based on the Malmquist productivity index (Färe and Grosskopf, [31]; Yao et al., [4]). Looking beyond the inherent differences between these models, we posit that their main objective is to address the changing patterns of efficiency scores along distinct time periods. While these approaches may be useful for decision making, they do not take into account an aggregated measure of efficiency that represents and provides a synthesis of multiperiod productive systems.
In this sense, a clear exception is found within the ambit of dynamic DEA models. Nemoto and Goto [32]s proposed a dynamic DEA model to assess the overall efficiency of a multiperiod production system. This overall efficiency can be viewed as price or economic efficiency. However, even in a particular period, the assumption of exact costs of individual inputs is unrealistic (Thompson et al., [33]). Moreover, the true monetary value (or exact discount factor) of an input in the time horizon remains unknown in practice (for more details see Jahanshahloo et al., [34]; Silva and Stefanou, [35]; Soleimani-damaneh, [36]).
Here, in this paper, we focused on Multiperiod Data Envelopment Analysis (MDEA) in which the aggregative efficiency in the context of time serial data is measured. Any previous information on prices or input and output weights across multiple periods are required. Hence, a Multiperiod Aggregative Efficiency (MAE) that corresponds conceptually to a technical (but not price or economic) efficiency of multiperiod production units can be delivered. This multiperiod technical efficiency, economic-free, is deemed necessary to proxy in an unbiased way for the technological pattern of a thermal plant over the course of time in terms of its cost structure, which is treated here in the ambit of the contextual variable set.
Besides, a distinctive feature of the model here developed is regards super-efficiency. The super-efficiency concept is traditionally a method used in DEA to break the ties between fully efficient units. Putting it into other words, the super-efficiency approach is an alternative to make a better discrimination for each Decision Making Unit (DMU). Suppose we have n DMUs and there are L periods, t = 1 , 2 , , L , and in each period D M U j , j = 1 , 2 , , n , consume m inputs, x i j t , i = 1 , 2 , , m , to produce s outputs, y r j t , r = 1 , 2 , , s . In order to compute the multiperiod aggregative efficiency, abbreviated as MAE in the context of time serial data, Sam Park and Park [37] proposed a two-phase DEA model (PP-model). The phase-I of the PP-model related to D M U o ( o { 1 , 2 , , n } ) is as follows:
Phase   ( I )   max   ψ o s . t . j = 1 n μ j t y r j t ψ o y r o t         r = 1 , 2 , , s         t = 1 , 2 , , L j = 1 n μ j t x i j t x i o t         i = 1 , 2 , , m         t = 1 , 2 , , L μ j t 0 ,   j = 1 , 2 , , n       t = 1 , 2 , , L
Indeed, the above model is an aggregated output-oriented CCR model that evaluates D M U o within all L periods simultaneously [38]. Let ψ o * be the optimal value of Equation (1), if ψ o * = 1 we say D M U o is weakly efficient. The following model is solved in phase-II of the PP-model.
Phase   ( II )   max t = 1 L r = 1 s s r + t + t = 1 L i = 1 m s i t s . t . j = 1 n μ j t y r j t s r + t = ψ o * y r o t         r = 1 , 2 , , s         t = 1 , 2 , , L j = 1 n μ j t x i j t + s i t = x i o t ,     i = 1 , 2 , , m ,     t = 1 , 2 , , L μ j t 0 ,     j = 1 , 2 , , n ,     t = 1 , 2 , , L s r + t 0 ,     r = 1 , 2 , , s ,     t = 1 , 2 , , L s i t 0 ,     i = 1 , 2 , , m ,     t = 1 , 2 , , L
Generally there is more than one DMU with ψ o * = 1 per Equation (1); therefore, it is necessary to propose a ranking model to rank efficient DMUs. In this study inspired by the AP-model (Andersen and Petersen, [39]) in which D M U o is removed from the production possibility set, Equations (1) and (2) are converted to the following models, respectively.
Phase   ( I )   max   ψ o s . t . j = 1 , j o n μ j t y r j t ψ o y r o t               r = 1 , 2 , , s             t = 1 , 2 , , L j = 1 , j o n μ j t x i j t x i o t i = 1 , 2 , , m t = 1 , 2 , , L μ j t 0 j o = 1 , 2 , , L
Phase   ( I I )     max t = 1 L r = 1 s s r + t + t = 1 L i = 1 m s i t s . t . j = 1 , j o n μ j t y r j t s r + t = ψ o * y r o t r = 1 , 2 , , s t = 1 , 2 , , L j = 1 , j o n μ j t x i j t + s i t = x i o t i = 1 , 2 , , m t = 1 , 2 , , L μ j t 0 j o t = 1 , 2 , , L s r + t 0 r = 1 , 2 , , s t = 1 , 2 , , L s i t 0 i = 1 , 2 , , m t = 1 , 2 , , L
The drawback of the above models is that they consider an overall efficiency score for each DMU, whereas the behavior of a DMU may change from one period to another. Consequently, it is reasonable to consider a different efficiency score for each time period; therefore, the following models are proposed.
Phase   ( I )     max     Ψ o = 1 L t = 1 L ψ o t s . t . j = 1 , j o n μ j t y r j t ψ o t y r o t         r = 1 , 2 , , s         t = 1 , 2 , , L j = 1 , j o n μ j t x i j t x i o t         i = 1 , 2 , , m         t = 1 , 2 , , L μ j t 0         j o         t = 1 , 2 , , L
Phase   ( I I )     max     ξ o = t = 1 L r = 1 s s r + t + t = 1 L i = 1 m s i t s . t . j = 1 , j o n μ j t y r j t s r + t = ψ o * t y r o t         r = 1 , 2 , , s         t = 1 , 2 , , L j = 1 , j o n μ j t x i j t + s i t = x i o t         i = 1 , 2 , , m         t = 1 , 2 , , L μ j t 0         j o         t = 1 , 2 , , L s r + t 0         r = 1 , 2 , , s         t = 1 , 2 , , L s i t 0         i = 1 , 2 , , m         t = 1 , 2 , , L
Here ψ o t represents the efficiency score of D M U o in t -th, t = 1 , 2 , , L , period. The objective function of Equation (5) is the average efficiency of D M U o within all L periods. The advantage of Equation (5) is that not only does it maximize the average efficiency of D M U o , it also computes the ψ o t . The objective of Equation (6) is to find a solution that maximizes the sum of input excesses and output shortfalls while keeping ψ o t = ψ o * t . Let ξ o * be the optimal value of Equation (6), then we have the following definitions.
Definition 1.
D M U o is called efficient if Ψ o * 1 and max 1 t L { ψ o * t } 1 . If Ψ o * > 1 or max 1 t L { ψ o * t } > 1 , we say D M U o is inefficient.
Definition 2.
D M U o is called weakly efficient if Ψ o * = 1 and ξ o * 0 .
Conventional DEA models rely on the assumption that outputs have to be maximized. However, it was mentioned already in the literature that the production process might also produce undesirable outputs (Mariano et al., [40]; Ozkan and Ulutas, [41]; Scheel, [42]). In Equation (1), the outputs of the DMUs are all desirable outputs and Equation (1) cannot be applied when some of the outputs are undesirable. Now, assume that each DMU j uses m inputs to produce s desirable outputs and k undesirable outputs. The inputs, desirable outputs, and undesirable outputs of each D M U j , j = 1 , 2 , , n in period t are defined as x i j t , i = 1 , 2 , , m , y r j g t , r = 1 , 2 , , s , and y p j b t , p = 1 , 2 , , k , respectively. Similar to Seiford and Zhu [43] and Hadi-Vencheh et al. [44], we assume strong disposability of the undesirable outputs. The data of the undesirable outputs are then transformed using the following Equation (7).
y p b t ^ = y p b t + w p t > 0 ,   p = 1 , 2 , , k ,   t = 1 , 2 , , L
In Equation (7), w t is a positive vector for period t , which can be used to let all the negative undesirable outputs be positive. Considering the transformed data, Equations (5) and (6) are converted to Equations (8) and (9), respectively.
Phase   ( I )     max     Ψ o = 1 L t = 1 L ψ o t s . t . j = 1 , j o n μ j t y r j t ψ o t y r o t         r = 1 , 2 , , s         t = 1 , 2 , , L j = 1 , j o n μ j t x i j t x i o t         i = 1 , 2 , , m         t = 1 , 2 , , L j = 1 , j o n μ j t y p j b t ^ y p o b t ^         p = 1 , 2 , , k         t = 1 , 2 , , L μ j t 0         j o         t = 1 , 2 , , L
Phase   ( I I )     max   ξ o = t = 1 L r = 1 s s r + g t + t = 1 L p = 1 k s p b t + t = 1 L i = 1 m s i t s . t . j = 1 , j o n μ j t y r j t s r + g t = ψ o * t y r o t         r = 1 , 2 , , s           t = 1 , 2 , , L j = 1 , j o n μ j t x i j t + s i t = x i o t         i = 1 , 2 , , m           t = 1 , 2 , , L j = 1 , j o n μ j t y p j b t ^ + s p b t = y p o b t ^         p = 1 , 2 , , k           t = 1 , 2 , , L μ j t 0         j o         t = 1 , 2 , , L s r + g t 0         r = 1 , 2 , , s         t = 1 , 2 , , L s p b t 0         p = 1 , 2 , , k         t = 1 , 2 , , L s i t 0         i = 1 , 2 , , m         t = 1 , 2 , , L
As can be seen from (8) and (9), the undesirable outputs of the DMUs are considered as inputs when evaluating the performance of the DMUs.
As regards the model fitting the Angolan thermal plant industry analyzed here, it is worth noting the justification for some methodological choices with respect to the returns-to-scale assumption, the model productive orientation, and the strong disposability of pollutant assumption. First, a constant returns-to-scale assumption was adopted here because when compared to international standards, even the largest Angolan plant is small when compared to their USA, UE, or even South African counterparts. Besides, as shown by results analyzed and discussed in Section 5, these thermal plants are labor intensive and in such cases the returns-to-scale effect tends to be behave mostly linearly with respect to the number of employees. Second, an output orientation was chosen due to the socioeconomic characteristics of the country and, why not, for the African region as a whole since energy shortages are frequent and interruptions in energy supply cannot be sourced by alternative transmission lines given that transmission networks within Angola and between African countries barely exist. Third, a strong disposability of undesirable pollutants was considered as an adequate assumption due to the existence of technological devices for reducing CO2 emissions and water discharges without affecting energy production levels.

5. Data and Bootstrapped Regression Trees

5.1. Data

Data for 2010–2016 encompassing the operation of 32 Angolan thermal plants were obtained from ENE-EP. The inputs of these thermal plants are those commonly used in previous papers and encompassed, besides the costs of investment, fuel, and labor, the productive capacity and the number of employees. It is worth mentioning that productive capacity was adopted here to proxy a relevant resource constraint, thus modeling that energy production cannot be indefinitely expanded in the short and medium terms. In fact, productive capacity is considered here as a fixed input required to produce energy despite severe limitations with respect to its short/medium term variations.
As regards the outputs, energy production is used as the sole desirable output, while on the other hand carbon dioxide emissions and discharged polluted water are considered as the undesirable outputs. It is worth mentioning that, due to a lack of completeness of the dataset, fuel costs had to be estimated as a proportion of the investment costs and the labor costs, observing accepted technical standards (for a conventional coal power plant, capital costs lie around 65%, whereas fuels costs are about 30% (EIA, [45,46]). For conventional gas power plants, capital costs are about 32%, whereas fuels costs are about 61% (EIA, [45,46])). For modern combined cycle gas power plants, the capital costs are 23% and fuel costs increase to over 70% with labor costs tending to be negligible (EIA, [45,46])). Table 2 presents their descriptive statistics and the respective dataset used is given in Appendix A.
In addition, a number of contextual variables are used to proxy the productive technology. They are also described in Table 2 and involve major cost structure elements. These variables are fourfold: the capacity cost per MW (calculated as the ratio of the logarithms between the total cost of the plant and its productive capacity), the labor cost per employee (calculated as the ratio of the logarithms between total salaries paid and number of employees), the capital cost (calculated as the ratio of the logarithms between amortizations and total assets), and the cost–asset ratio (calculated as the ratio of the logarithms between total costs and total assets). Lastly, the capital–labor ratio is calculated as the ratio of the logarithms between total assets and total salaries paid to employees. It is also important to mention that principles of accrual accounting were adopted here to compute the capital cost as the amortizations to total assets ratio. In fact, amortization, depreciation, and depletion are methods that are used to prorate the cost of a specific type of asset over the asset’s life. This prorated cost yields, therefore, an accounting proxy for the capital cost.

5.2. Bootstrapped Regression Trees

Tree methods were first used by researchers [2,47,48] and have gained popularity through the major theoretical and practical contributions of Breiman et al. [5]. They involve stratifying and segmenting the predictor space into a number of simple regions (James et al., [49]). These features are particularly useful since (a) efficiency scores reflect uncertainty derived not only from vagueness in input/output collection and (b) explanatory variables may be endogenous or exogenous. Readers should refer to Faraway [2], Opitz and Maclin [50], Polikar [51], and Torgo [52] on how to resample (ensemble) trees using bagging.
Bagging (“bootstrapping aggregation”) is a bootstrap ensemble method introduced by Breiman [53] and combines predictors across different subsets of the training data. The R functions used to perform such bootstrapped regression tree analysis and their respective default values used for the analysis are presented in Table 3 (Ledolter, [54]).

6. Results and Discussion

The levels of super-efficiency computed for the Angolan thermal power plant sample using the proposed super-efficiency model for undesirable outputs are presented in Figure 1, Figure 2, Figure 3 and Figure 4. The full rank of DMU scores is given in Appendix B. Readers should recall that a DMU is considered efficient (inefficient) if super-efficiency scores are below (above) 1. If super-efficiency score is exactly 1, the DMU is classified as weakly efficient. As displayed in Figure 1, pooled efficiency scores are strongly concentrated around 1 and inefficiency prevails in Angolan thermal power plants. Super-efficiency levels also appear to be stagnant over the course of the years, as suggested by Figure 2, even though a slight decrease in efficiency is seen from 2014 on. Although this stagnant behavior may be justified by the fact that Angola thermal power generation plants are publicly owned and controlled by ENE-EP, a state company (cf. Section 2), the wide dispersion of super-efficiency scores suggest that inefficiency may be driven by different technological patterns in thermal energy production reflected in cost-structure variables, and that they may also impact the emission of pollutants. Based on the literature review, the evidence implies that the smaller, older, and coal thermal plants are less efficient and more polluting than the larger, newer plants that burn gas. The former would, therefore, be located above one in Figure 1 while the latter below one. One is the super-efficiency threshold that divides efficient from non-efficient plants.
Figure 3 presents a scatterplot panel of efficiency levels against the two undesirable outputs. As expected, super-efficiency scores fluctuate around 1, although not symmetrically. Readers should recall that inefficient plants present scores higher than 1.0 and that a score of 0.6, for instance, indicates a less efficient plant than another one with a score of 0.5. Figure 3 reveals that inefficient thermal power plants tended to emit more CO2, although this direct relationship did not appear to be so strong in terms of polluted water. In fact, the following results could be observed when comparing the groups between efficient and inefficient plants. While the efficient plants present a CO2 emission level that was 1.08% lower than those observed in the inefficient plants, the discharge of polluted water between both groups presented even a smaller variation (around 0.33% less).
Similarly, Figure 4 also presents a scatterplot panel of efficiency levels, but now against cost structure variables. Figure 4 clearly shows that inefficient energy generation prevailed when Angolan thermal plants presented a lower KL ratio, meaning that they are labor intensive, and had a higher cost–asset ratio, which means that they have higher operating costs in comparison to the value of their assets. As regards the lower KL ratios, the state-controlled operation of Angolan thermal plants may explain the excessive number of highly paid employees in relation to the size/scale of the plant. On the other hand, higher cost–asset ratios may reflect the operation of older, smaller, poorly maintained plants. With this cost structure in mind, the adoption of costly technologies for controlling CO2 emissions may be enhanced by opportunities that emerge from rebalancing KL and cost–asset ratios. This is necessary because adopting carbon capture devices would certainly increase the cost to asset ratio. Therefore, labor expenses should be rightsized first in order to open room in cost expenditure before adopting such an antipollutant measure.
Firstly, an adequate equilibrium between a rightsized labor force and the intensity of capital seems to be the cornerstone for improving efficiency levels while simultaneously controlling the discharge of undesirable outputs as a byproduct of the energy generation process. The drop in oil prices imposed strong budgetary restrictions upon the Angolan economy causing the classical conflict between labor and capital to increase with respect to scarce resource allocation. Secondly, with respect to the capital investments required for controlling carbon dioxide emissions, it is deemed necessary to apprehend how economically feasible these investments are for a thermal plant and how they impact relevant ratios such as the capital–labor and the cost–asset ratios. As a matter of fact, thermal energy generation is one of the biggest causers of the greenhouse effect on a worldwide basis. Except for CO2, all other emissions from a thermal plant can be mitigated with the technology available at a feasible cost (http://www.brighthubengineering.com/power-plants/57788-power-plant-emissions/). This happens because carbon dioxide is an unavoidable part of the thermal generation process. Therefore, systems for capturing carbon dioxide emissions consists of a costly alternative for reducing pollutant emissions in the context of the investments required for building up and/or renovating a thermal plant (U.S. Department of Energy (DOE) and U.S. National Energy Technology Laboratory (NETL). 2010. DOE/NETL Carbon Dioxide Capture and Storage RD&D Roadmap. http://www.netl.doe.gov/File%20Library/Research/Carbon%20Seq/Reference%20Shelf/CCSRoadmap.pdf), but they are often economically feasible in newer plants.
These results are confirmed by the bootstrapped regression trees presented in Figure 5 and Figure 6. While Figure 5 depicts the bootstrapped regression tree structure and its thresholds at each tree node, Figure 6 shows the most impacting variables in terms of overall increase in the Mean Squared Errors (MSE). The interpretation of a regression tree is very straightforward in terms of allowing the decision-maker to segment the results. For example, the first branch of the tree presented in Figure 5 states that “if labor cost is lower than 0.33 and capacity cost is lower than 0.80, and KL ratio is lower than 1.37, then the average plant efficiency is 0.962 (almost weak-efficient)”. The other branches read similarly. Both figures suggest that capacity cost and labor cost are the most impacting variables on Angolan thermal power plant efficiency.
Policy implications for Angolan thermal plants suggest the need for a better training of the workforce, downsizing of personnel in order to keep labor costs under control, and investing in carbon capture equipment. Since it is expensive and impacts the plant’s capacity cost, such equipment can be acquired by rebalancing the KL equilibrium in Angolan plants so that total operating costs are kept under control. There are, however, technological limits to adopting such measures, which are revealed by the capital cost. Lower values of capacity cost, often related to older amortized plants, may present physical constraints to the deployment of newer carbon capture technologies. In such cases, efficiency improvements may be confined to the traditional conversion from diesel to combined cycle gas, which has already occurred in some Angolan thermal plants.

7. Conclusions

This research assessed the thermal plants in Angola in terms of their energetic efficiency by jointly applying a novel super-efficiency DEA model that handles undesirable outputs and bootstrapped regression trees that discriminate productive technologies based on cost structure variables. This novel super-efficiency model makes it possible to capture subtle productivity variations in a state-owned industry where managerial practices tend to be quite similar over the course of time. This paper gives insights on how the technological patterns or productive technologies of thermal energy generation are reflected in the cost structure variables of each plant by means of efficiency levels, which constitutes a relatively novel approach not only in the energy efficiency strand, but also on infrastructure efficiency.
Efficiency levels were computed based on three outputs (energy production, carbon dioxide emissions, and discharged polluted water) and on five inputs (fuel, investment, labor costs, plant capacity, and number of employees). The findings suggest that efficiency levels of fuel consumption and undesirable emissions included are mostly affected by the capacity cost and the labor cost, which are the reflex of rightsizing and training the workforce in parallel with adopting expensive carbon capture devices. Specifically with respect to the pattern of pollutant emissions/discharges, CO2 emissions appear to be more impacted by the technological pattern of the power plant than the level of discharged polluted water, which may suggest that carbon capture technologies have evolved and can be deployed faster than technologies for recycling water in the energy generation productive processes.
Limitations of this research are related to the very nature of case studies built on the evidence of single countries. Although these results cannot be generalized to other countries with different regulatory regimes, some useful lessons for conducting similar research in other countries have been learned. It may be advisable to focus on capital, labor, and operating expenses and their countervailing forces while seeking opportunities for adopting antipollutant technologies. Further research should confirm these results in other environments.

Author Contributions

The authors have the equal contributions in this study. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by FAPERJ—Cientista do Nosso Estado.

Acknowledgments

The authors would like to tank two anonymous reviewers for their helpful comments and suggestions which improved the first draft of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Angolan Thermal Plants Dataset

DMUYearEnergy Production (MWh)CO2 Emission (Tons Per Year)Polluted Water (Liters Per Year)Capacity (MW)Number of EmployeesEmployee Cost Per Year (Log)Fuel Cost Per Year (Log)Investment Cost (Log)KL Ratio (Log)Capacity Cost per MW (Log)Labor Cost Per Year (Log)Capital Cost Per Year (Log)Cost–Asset Ratio (Log)Fuel_Type
Central Termica de Xitoto Novo2016106.000 2.0001232.614 70.00045.0002.7689.71416.6603.1000.3400.3540.7005.230Diesel
Central Termica de Xitoto Antigo201610.0002.4121245.0068.00040.0002.7737.90113.0303.6302.6700.4000.6102.070Diesel
Central Termica do Aeroporto201610.0002.3921270.1628.00023.0002.71410.15617.5972.7302.6760.6560.8224.838Diesel
Central Termica de Arimba201640.6002.3721257.51635.00015.0002.7049.96217.2212.2000.6740.9960.7306.664Diesel
Central Termica de Anexo SE201630.0002.4321254.75820.00030.0002.77410.24117.7090.7901.1620.5340.7626.982Diesel
Central Termica de Kileva201645.6002.3381250.71630.00069.0002.7646.73710.71011.4080.7000.2300.510−2.420Diesel
Central Termical do Lobito201620.0002.3181228.6705.00031.0002.7217.71312.7064.8103.6200.4900.702−1.748Diesel
Central Termica de Biopio2016117.0002.2961238.28215.00046.0002.75012.12821.5060.6301.6480.3400.87010.592Diesel
Central Termica do Cunene201615.0002.2961234.1445.00020.0002.7298.15113.5721.1803.7700.7660.7202.540Diesel
Central Termica Do Kuando Kubango201630.0002.2561270.16210.00016.0002.72111.89421.067−0.5902.3100.9500.9129.340Diesel
Central Termica do Benfica201650.4002.2321282.92635.00034.0002.8018.44914.0981.1540.5300.4840.7602.250Diesel
Central Termica do Cuito201610.0002.2121268.8187.00042.0002.74311.03919.3341.5662.7620.3701.0002.250Diesel
Central Termica do Moxico201640.0002.1761250.72630.00023.0002.75912.35221.9450.9200.9500.6860.77012.058Diesel
Central Termica de Cazenga201640.0002.1561246.38830.00050.0002.77341.40117.9282.1600.8300.3200.7207.280Combined Cycle Gas
Central Termica do Caminho de Ferro de Luanda201696.8002.1361231.27250.00065.0002.78810.67418.5601.9020.5800.2500.64010.784Diesel
Central Termica do Morro Bento201630.0002.1121232.22415.00032.0002.7529.15215.5521.2101.4400.4900.7204.630Diesel
Central Termica do Morro da Luz201634.8002.0921245.00810.00025.0002.7477.87212.9965.7702.1660.6240.6001.940Diesel
Central Termica KM9201615.0002.0921246.39810.00035.0002.74811.58220.4162.3302.3740.4460.8607.300Diesel
Central Termica Benfica201620.0001.1361247.77817.00018.0002.72311.27219.8211.2701.3880.8460.8407.850Diesel
Central Termica Praia do Bispo201635.0002.0741249.16820.00015.0002.67110.93919.207−0.3301.0600.9640.9067.162Diesel
Central Termica dos Quarteis201630.0002.0521246.54830.00037.0002.7819.44516.1102.9360.7160.4360.7504.340Diesel
Central Termica Boavista I201632.0002.2521255.88015.00030.0002.7289.08415.4401.3801.4100.5100.7304.356Diesel
Central Termica Boavista II201620.0002.1921237.99610.00045.0002.7858.88114.9761.7042.0800.3600.7203.180Diesel
Central Termica Boavista III201622.6002.1481259.34016.00042.0002.7968.70014.6051.1601.3040.3900.7003.404Diesel
Central Termica Kassaki201610.0002.1281290.8288.00028.0002.7527.87813.0035.2142.5800.5600.6300.930Diesel
Central Termica CEEF201615.0002.1061396.92410.00020.0002.68110.06417.4471.6402.0820.7300.8386.434Diesel
Central Termica do Bengo201632.2002.0861472.14018.00025.0002.7478.93715.1282.0081.2180.6240.6904.418Diesel
Central Termica do Uige201620.0002.1281402.33215.00020.0002.7278.17913.6325.7101.4200.7640.6401.890Diesel
Central Termica Banza Congo201620.0002.0861368.17010.00035.0002.7347.44412.1544.8202.0600.4400.5900.680Diesel
Central Termica Soyo2016110.6002.4021898.21270.00015.0002.70442.21518.4043.7600.4020.9960.65410.232Combined Cycle Gas
Central Termica Landana201630.2002.3861462.91220.00030.0002.70810.40618.1052.4901.0600.5000.8544.204Diesel
Central Termica Malongo201645.0002.7201333.35040.00015.0002.6989.18315.6671.4200.5440.9900.7205.854Diesel
Central Termica de Xitoto Novo201593.0001.9401197.22270.00045.0002.7639.71116.6603.1000.3400.3520.7005.230Diesel
Central Termica de Xitoto Antigo201510.0002.3411209.2588.00040.0002.7737.90113.0303.4302.6700.4000.6102.070Diesel
Central Termica do Aeroporto201510.0002.3211233.6918.00023.0002.7099.70516.7012.7302.6730.6530.7814.319Diesel
Central Termica de Arimba201540.3002.3011221.40835.00015.0002.7019.80816.9142.2000.6620.9930.7306.127Diesel
Central Termica de Anexo SE201530.0002.3611218.72920.00030.0002.77010.12417.4790.7901.1560.5320.7566.721Diesel
Central Termica de Kileva201542.8002.2691214.80330.00069.0002.7646.73710.7109.3840.7000.2300.510-2.420Diesel
Central Termical do Lobito201520.0002.2491193.3905.00031.0002.7217.38812.0554.8103.6200.4900.666-1.874Diesel
Central Termica de Biopio2015103.5002.2281202.72615.00046.0002.75011.05219.3540.6301.5640.3400.8258.231Diesel
Central Termica do Cunene201515.0002.2281198.7075.00020.0002.7258.14913.5721.1803.7700.7630.7202.540Diesel
Central Termica Do Kuando Kubango201530.0002.1881233.69110.00016.0002.72111.12819.535-0.5902.2300.9500.8768.065Diesel
Central Termica do Benfica201545.2002.1661246.08835.00034.0002.7978.44714.0980.6420.5300.4820.7602.250Diesel
Central Termica do Cuito201510.0002.1461232.3847.00042.0002.74310.53418.3241.4532.7410.3700.9552.250Diesel
Central Termica do Moxico201540.0002.1131214.81330.00023.0002.75411.65720.5590.9200.8900.6830.77010.219Diesel
Central Termica de Cazenga201540.0002.0931210.59930.00050.0002.77341.40117.9281.8400.8300.3200.7207.280Combined Cycle Gas
Central Termica do Caminho de Ferro de Luanda201588.4002.0731195.91650.00065.0002.78810.43418.0801.7560.5650.2500.64010.127Diesel
Central Termica do Morro Bento201530.0002.0511196.84215.00032.0002.7528.99015.2281.2101.4400.4900.7054.630Diesel
Central Termica do Morro da Luz201532.4002.0311209.25910.00025.0002.7447.86112.9784.9702.1630.6220.6001.940Diesel
Central Termica KPM9201515.0002.0311210.60910.00035.0002.74110.81518.8892.3302.3320.4430.8106.295Diesel
Central Termica Benfica201520.0001.1031211.94917.00018.0002.72011.03019.3401.2701.3790.8430.8257.370Diesel
Central Termica Praia do Bispo201535.0002.0121213.29920.00015.0002.66910.74718.826-0.3301.0600.9620.8887.111Diesel
Central Termica dos Quarteis201530.0001.9911210.75430.00037.0002.7749.40816.0432.0730.7130.4330.7504.340Diesel
Central Termica Boavista I201531.0002.1861219.82015.00030.0002.7289.08415.4401.3801.4100.5100.7304.353Diesel
Central Termica Boavista II201520.0002.1261202.44810.00045.0002.7858.88114.9761.4622.0800.3600.7203.180Diesel
Central Termica Boavista III201521.3002.0841223.18016.00042.0002.7968.68914.5821.1601.3020.3900.7003.362Diesel
Central Termica Kassaki201510.0002.0641253.7648.00028.0002.7527.87813.0034.5672.5800.5600.6300.930Diesel
Central Termica CEEF201515.0002.0431356.81210.00020.0002.6819.63416.5871.6402.0760.7300.7995.402Diesel
Central Termica do Bengo201531.1002.0231429.87018.00025.0002.7448.84914.9541.8141.2040.6220.6904.409Diesel
Central Termica do Uige201520.0002.0641362.06615.00020.0002.7248.17813.6325.0001.4200.7620.6401.890Diesel
Central Termica Banza Congo201520.0002.0231328.88510.00035.0002.7347.44412.1544.6652.0600.4400.5900.680Diesel
Central Termica Soyo2015100.3002.3311843.70670.00015.0002.70139.64617.1223.7600.3810.9930.6428.441Combined Cycle Gas
Central Termica Landana201530.1002.3181420.90620.00030.0002.7089.96117.2142.4901.0600.5000.8123.907Diesel
Central Termica Malongo201545.0002.5101295.06540.00015.0002.6989.15415.6101.4200.5420.9900.7205.467Diesel
Central Termica de Xitoto Novo201480.0001.8801161.83070.00045.0002.7579.70816.6603.1000.3400.3500.7005.230Diesel
Central Termica de Xitoto Antigo201410.0002.2701173.5108.00040.0002.7737.90113.0303.2302.6700.4000.6102.070Diesel
Central Termica do Aeroporto201410.0002.2501197.2208.00023.0002.7059.25615.8062.7302.6700.6500.7403.800Diesel
Central Termica de Arimba201440.0002.2301185.30035.00015.0002.6989.65316.6082.2000.6500.9900.7305.590Diesel
Central Termica de Anexo SE201430.0002.2901182.70020.00030.0002.76610.00817.2500.7901.1500.5300.7506.460Diesel
Central Termica de Kileva201440.0002.2001178.89030.00069.0002.7646.73710.7107.3600.7000.2300.510-2.420Diesel
Central Termical do Lobito201420.0002.1801158.1105.00031.0002.7217.06211.4034.8103.6200.4900.630-2.000Diesel
Central Termica de Biopio201490.0002.1601167.17015.00046.0002.75010.03317.3160.6301.4800.3400.7805.870Diesel
Central Termica do Cunene201415.0002.1601163.2705.00020.0002.7218.14713.5721.1803.7700.7600.7202.540Diesel
Central Termica Do Kuando Kubango201430.0002.1201197.22010.00016.0002.72110.39118.060-0.5902.1500.9500.8406.790Diesel
Central Termica do Benfica201440.0002.100109.25035.00034.0002.7928.44514.0980.1300.5300.4800.7602.250Diesel
Central Termica do Cuito201410.0002.080195.9507.00042.0002.74310.03517.3261.3402.7200.3700.9102.250Diesel
Central Termica do Moxico201440.0002.0501178.90030.00023.0002.75010.96119.1730.9200.8300.6800.7708.380Diesel
Central Termica de Cazenga201440.0002.0301174.81030.00050.0002.77341.40117.9281.5200.8300.3200.7207.280Combined Cycle Gas
Central Termica do Caminho de Ferro de Luanda201480.0002.0101160.56050.00065.0002.78810.19417.6001.6100.5500.2500.6409.470Diesel
Central Termica do Morro Bento201430.0001.9901161.46015.00032.0002.7528.82814.9041.2101.4400.4900.6904.630Diesel
Central Termica do Morro da Luz201430.0001.9701173.51010.00025.0002.7417.85012.9604.1702.1600.6200.6001.940Diesel
Central Termica KM9201415.0001.9701174.82010.00035.0002.73410.06917.4042.3302.2900.4400.7605.290Diesel
Central Termica Benfica201420.0001.0701176.12017.00018.0002.71610.79018.8651.2701.3700.8400.8106.890Diesel
Central Termica Praia do Bispo201435.0001.9501177.43020.00015.0002.66710.55618.444-0.3301.0600.9600.8707.060Diesel
Central Termica dos Quarteis201430.0001.9301174.96030.00037.0002.7679.37115.9751.2100.7100.4300.7504.340Diesel
Central Termica Boavista I201430.0002.1201183.76015.00030.0002.7289.08415.4401.3801.4100.5100.7304.350Diesel
Central Termica Boavista II201420.0002.0601166.90010.00045.0002.7858.88114.9761.2202.0800.3600.7203.180Diesel
Central Termica Boavista III201420.0002.0201187.02016.00042.0002.7968.67814.5601.1601.3000.3900.7003.320Diesel
Central Termica Kassaki201410.0002.0001216.7008.00028.0002.7527.87813.0033.9202.5800.5600.6300.930Diesel
Central Termica CEEF201415.0001.9801316.70010.00020.0002.6819.20715.7321.6402.0700.7300.7604.370Diesel
Central Termica do Bengo201430.0001.9601387.60018.00025.0002.7418.76014.7801.6201.1900.6200.6904.400Diesel
Central Termica do Uige201420.0002.0001321.80015.00020.0002.7218.17713.6324.2901.4200.7600.6401.890Diesel
Central Termica Banza Congo201420.0001.9601289.60010.00035.0002.7347.44412.1544.5102.0600.4400.5900.680Diesel
Central Termica Soyo201490.0002.2601789.20070.00015.0002.69837.14815.8763.7600.3600.9900.6306.650Combined Cycle Gas
Central Termica Landana201430.0002.2501378.90020.00030.0002.7089.51616.3242.4901.0600.5000.7703.610Diesel
Central Termica Malongo201445.0002.3001256.78040.00015.0002.6989.12515.5521.4200.5400.9900.7205.080Diesel
Central Termica de Xitoto Novo201381.0001.8201138.60070.00045.0002.7289.69416.6603.1000.3400.3400.7005.230Diesel
Central Termica de Xitoto Antigo201311.0002.2001150.0408.00040.0002.7737.90113.0303.2302.6700.4000.6102.070Diesel
Central Termica do Aeroporto201311.0002.1801173.2708.00023.0002.7059.25615.8062.7302.6700.6500.7403.800Diesel
Central Termica de Arimba201341.0002.1601161.60035.00015.0002.6989.65316.6082.2000.6500.9900.7305.590Diesel
Central Termica de Anexo SE201331.0002.2201159.05020.00030.0002.7479.99917.2500.7901.1500.5200.7506.460Diesel
Central Termica de Kileva201340.0002.1401155.31030.00069.0002.7646.73710.7107.3600.7000.2300.510-2.420Diesel
Central Termical do Lobito201321.0002.1201134.9505.00031.0002.7217.06211.4034.8103.6200.4900.630-2.000Diesel
Central Termica de Biopio201391.0002.0901143.82015.00046.0002.75010.03317.3160.6301.4800.3400.7805.870Diesel
Central Termica do Cunene201316.0002.0901140.0005.00020.0002.7218.14713.5721.1803.7700.7600.7202.540Diesel
Central Termica Do Kuando Kubango201331.0002.0501173.27010.00016.0002.72110.34917.976-0.5902.1400.9500.8406.790Diesel
Central Termica do Benfica201341.0002.0301185.06035.00034.0002.7928.44514.0980.1300.5300.4800.7602.250Diesel
Central Termica do Cuito201311.0002.0101172.0307.00042.0002.74310.03517.3261.3402.7200.3700.9102.250Diesel
Central Termica do Moxico201341.0001.9901155.32030.00023.0002.75010.96119.1730.9200.8300.6800.7708.380Diesel
Central Termica de Cazenga201341.0001.9701151.32030.00050.0002.77341.40117.9281.5200.8300.3200.7207.280Combined Cycle Gas
Central Termica do Caminho de Ferro de Luanda201381.0001.9501137.35050.00065.0002.78810.19417.6001.6100.5500.2500.6409.470Diesel
Central Termica do Morro Bento201331.0001.9301138.23015.00032.0002.7528.82814.9041.2101.4400.4900.6904.630Diesel
Central Termica do Morro da Luz201331.0001.9101150.04010.00025.0002.7417.85012.9604.1702.1600.6200.6001.940Diesel
Central Termica KPM9201315.0001.9101151.32010.00035.0002.73410.06917.4042.3302.2900.4400.7605.290Diesel
Central Termica Benfica201321.0001.0401152.60017.00018.0002.71610.79018.8651.2701.3700.8400.8106.890Diesel
Central Termica Praia do Bispo201335.0001.9001153.88020.00015.0002.66710.55618.444-0.3301.0600.9600.8707.060Diesel
Central Termica dos Quarteis201331.0001.8801151.46030.00037.0002.7679.37115.9751.2100.7100.4300.7504.340Diesel
Central Termica Boavista I201331.0002.0601160.09015.00030.0002.7289.08415.4401.3801.4100.5100.7304.350Diesel
Central Termica Boavista II201321.0002.0001143.56010.00045.0002.7858.88114.9761.2202.0800.3600.7203.180Diesel
Central Termica Boavista III201321.0001.9601163.28016.00042.0002.7968.67814.5601.1601.3000.3900.7003.320Diesel
Central Termica Kassaki201311.0001.9401192.3708.00028.0002.7527.85312.9533.9202.5700.5600.6300.930Diesel
Central Termica CEEF201316.0001.9201290.37010.00020.0002.6819.20715.7321.6402.0700.7300.7604.370Diesel
Central Termica do Bengo201331.0001.9001359.85018.00025.0002.7258.75214.7801.6201.1900.6100.6904.400Diesel
Central Termica do Uige201321.0001.9401295.36015.00020.0002.7088.17013.6324.2901.4200.7500.6401.890Diesel
Central Termica Banza Congo201321.0001.9001263.81010.00035.0002.7347.44412.1544.5102.0600.4400.5900.680Diesel
Central Termica Soyo201391.0002.1901753.42070.00015.0002.69837.14815.8763.7600.3600.9900.6306.650Combined Cycle Gas
Central Termica Landana201331.0002.1801351.32020.00030.0002.7089.51616.3242.4901.0600.5000.7703.610Diesel
Central Termica Malongo201346.0002.2301231.64040.00015.0002.6988.98115.2641.4200.5300.9900.7205.080Diesel
Central Termica de Xitoto Novo201281.0002.2601115.83070.00045.0002.7289.69416.6603.1000.3400.3400.7005.240Diesel
Central Termica de Xitoto Antigo201211.0001.7801127.0408.00040.0002.7739.73016.6882.5402.9800.4000.7005.230Diesel
Central Termica do Aeroporto201211.0002.1601149.8108.00023.0002.6897.85913.0304.4602.6700.6400.6102.070Diesel
Central Termica de Arimba201241.0002.1401138.37035.00015.0002.6889.24315.7992.7800.6100.9800.7403.800Diesel
Central Termica de Anexo SE201231.0002.1201135.87020.00030.0002.7479.62316.4981.2501.1300.5200.7305.590Diesel
Central Termica de Kileva201240.0002.1801132.20030.00069.0002.76410.04517.3250.5900.7700.2300.7506.460Diesel
Central Termical do Lobito201221.0002.0901112.2505.00031.0002.7216.67710.6348.0204.1700.4900.510-2.420Diesel
Central Termica de Biopio201291.0002.0701120.95015.00046.0002.7507.04511.3404.4001.2000.3400.630-2.000Diesel
Central Termica do Cunene201216.0002.0501117.2005.00020.0002.7219.99917.2771.1904.4300.7600.7805.870Diesel
Central Termica Do Kuando Kubango201231.0002.0501149.81010.00016.0002.7218.12913.5361.0701.8800.9500.7202.540Diesel
Central Termica do Benfica201241.0002.0101161.36035.00034.0002.79210.36317.934-1.6000.6100.4800.8406.790Diesel
Central Termica do Cuito201211.0001.9901148.5907.00042.0002.7438.44714.1510.9402.6600.3700.7602.250Diesel
Central Termica do Moxico201241.0001.9701132.22030.00023.0002.7509.97417.1991.1800.6300.6800.9102.250Diesel
Central Termica de Cazenga201241.0001.9501128.29030.00050.0002.77343.89119.1730.4300.8300.3200.7708.380Combined Cycle Gas
Central Termica do Caminho de Ferro de Luanda201281.0001.9301114.60050.00065.0002.78810.39418.0001.1000.5000.2500.7207.280Diesel
Central Termica do Morro Bento201231.0001.9101115.47015.00032.0002.75210.20817.6642.4001.8400.4900.6409.470Diesel
Central Termica do Morro da Luz201231.0001.9001127.04010.00025.0002.7418.78814.8351.4802.1500.6200.6904.630Diesel
Central Termica KPM9201215.0001.8801128.29010.00035.0002.7347.81712.9004.3602.1500.4400.6001.940Diesel
Central Termica Benfica201221.0001.8801129.55017.00018.0002.71610.07917.4422.5401.3500.8400.7605.290Diesel
Central Termica Praia do Bispo201235.0001.0201130.80020.00015.0002.66710.73018.7921.9101.1600.9600.8106.890Diesel
Central Termica dos Quarteis201231.0001.8601128.43030.00037.0002.76710.51818.270-1.7000.7000.4300.8707.060Diesel
Central Termica Boavista I201231.0001.8401136.89015.00030.0002.7289.35115.9751.5901.4200.5100.7504.340Diesel
Central Termica Boavista II201221.0002.0201120.69010.00045.0002.7859.09415.4030.3902.1100.3600.7304.350Diesel
Central Termica Boavista III201221.0001.9601140.02016.00042.0002.7968.88614.9761.2901.3000.3900.7203.180Diesel
Central Termica Kassaki201211.0001.9201168.5208.00028.0002.7528.65614.5601.7302.6000.5600.7003.320Diesel
Central Termica CEEF201216.0001.9001264.56010.00020.0002.6817.83012.9784.9902.0600.7300.6300.930Diesel
Central Termica do Bengo201231.0001.8801332.65018.00025.0002.7259.22815.7320.9501.1500.6100.7604.370Diesel
Central Termica do Uige201221.0001.8601269.46015.00020.0002.7088.70314.6971.9001.4200.7500.6904.400Diesel
Central Termica Banza Congo201221.0001.9001238.53010.00035.0002.7348.15113.5683.9702.1200.4400.6401.890Diesel
Central Termica Soyo201291.0001.8601718.35070.00015.0002.69829.35011.9775.0600.2900.9900.5900.680Combined Cycle Gas
Central Termica Landana201231.0002.1401324.30020.00030.0002.7088.03213.3563.4801.0600.5000.6302.620Diesel
Central Termica Malongo201246.0002.1401207.01040.00015.0002.6989.51116.3242.7300.5300.9900.7703.790Diesel
Central Termica de Xitoto Novo201185.0001.7001064.50070.00045.0002.7289.69416.6603.1300.3400.3400.7005.240Diesel
Central Termica de Xitoto Antigo201113.0002.0501075.2008.00040.0002.7739.70216.6322.5702.9700.4000.7005.230Diesel
Central Termica do Aeroporto201113.0002.0301096.9208.00023.0002.6897.83512.9814.4902.6600.6400.6102.070Diesel
Central Termica de Arimba201146.0002.0101086.00035.00015.0002.6889.24315.7992.8100.6100.9800.7403.800Diesel
Central Termica de Anexo SE201136.0002.0701083.62020.00030.0002.7479.62316.4981.2801.1300.5200.7305.590Diesel
Central Termica de Kileva201144.0001.9901080.12030.00069.0002.7649.93217.1000.6200.7600.2300.7506.460Diesel
Central Termical do Lobito201124.0001.9701061.0905.00031.0002.7216.67710.6348.0504.1700.4900.510-2.420Diesel
Central Termica de Biopio2011100.0001.9501069.38015.00046.0002.7507.04511.3404.4301.2000.3400.630-2.000Diesel
Central Termica do Cunene201118.0001.9501065.8105.00020.0002.7089.99317.2771.2204.4300.7500.7805.870Diesel
Central Termica Do Kuando Kubango201134.0001.9101096.92010.00016.0002.7218.12913.5361.1001.8800.9500.7202.540Diesel
Central Termica do Benfica201145.0001.8901107.94035.00034.0002.79210.36317.934-1.5700.6100.4800.8406.790Diesel
Central Termica do Cuito201112.0001.8701095.7507.00042.0002.7438.42114.0980.9702.6500.3700.7602.250Diesel
Central Termica do Moxico201145.0001.8601080.13030.00023.0002.7359.96717.1991.2100.6300.6700.9102.250Diesel
Central Termica de Cazenga201145.0001.8401076.39030.00050.0002.77343.89119.1730.4600.8300.3200.7708.380Combined Cycle Gas
Central Termica do Caminho de Ferro de Luanda201189.0001.8201063.33050.00065.0002.78810.39418.0001.1300.5000.2500.7207.280Diesel
Central Termica do Morro Bento201134.0001.8001064.15015.00032.0002.75210.20817.6642.4301.8400.4900.6409.470Diesel
Central Termica do Morro da Luz201135.0001.7801075.19010.00025.0002.7418.78814.8351.5102.1500.6200.6904.630Diesel
Central Termica KPM9201116.0001.7801076.39010.00035.0002.7117.80612.9004.3902.1500.4300.6001.940Diesel
Central Termica Benfica201123.0000.9701077.59017.00018.0002.70410.00817.3132.5701.3400.8300.7605.290Diesel
Central Termica Praia do Bispo201138.0001.7601078.79020.00015.0002.66710.73018.7921.9401.1600.9600.8106.890Diesel
Central Termica dos Quarteis201134.0001.7501076.52030.00037.0002.74310.50718.270-1.6600.7000.4200.8707.060Diesel
Central Termica Boavista I201135.0001.9201084.59015.00030.0002.7289.35115.9751.6201.4200.5100.7504.340Diesel
Central Termica Boavista II201123.0001.8601069.14010.00045.0002.7859.09415.4030.4202.1100.3600.7304.350Diesel
Central Termica Boavista III201124.0001.8201087.58016.00042.0002.7968.88614.9761.3201.3000.3900.7203.180Diesel
Central Termica Kassaki201112.0001.8001114.7708.00028.0002.7528.62814.5041.7702.5900.5600.7003.320Diesel
Central Termica CEEF201117.0001.7901206.39010.00020.0002.6817.79812.9155.0202.0500.7300.6300.930Diesel
Central Termica do Bengo201134.0001.7701271.35018.00025.0002.7259.22815.7320.9801.1500.6100.7604.370Diesel
Central Termica do Uige201123.0001.8001211.06015.00020.0002.7088.70314.6971.9301.4200.7500.6904.400Diesel
Central Termica Banza Congo201124.0001.7701181.56010.00035.0002.7348.15113.5684.0002.1200.4400.6401.890Diesel
Central Termica Soyo2011100.0002.0401639.30070.00015.0002.68829.33011.9775.0900.2900.9800.5900.680Combined Cycle Gas
Central Termica Landana201134.0002.0401263.38020.00030.0002.7088.03213.3563.5101.0600.5000.6302.620Diesel
Central Termica Malongo201185.0001.6101151.49040.00015.0002.6989.51116.3242.7600.5300.9900.7703.790Diesel
Central Termica de Xitoto Novo201013.0001.6401021.92070.00045.0002.7289.69416.6603.1200.3400.3400.7005.240Diesel
Central Termica de Xitoto Antigo201013.0001.9901032.1908.00040.0002.7739.70216.6322.5602.9700.4000.7005.230Diesel
Central Termica do Aeroporto201046.0001.9701053.0408.00023.0002.6897.72912.7684.4802.6600.6400.6002.070Diesel
Central Termica de Arimba201036.0001.9501042.56035.00015.0002.6889.24315.7992.8000.6100.9800.7403.800Diesel
Central Termica de Anexo SE201039.6002.0101040.27020.00030.0002.7479.62316.4981.2701.1300.5200.7305.590Diesel
Central Termica de Kileva201024.0001.9301036.92030.00069.0002.7649.81816.8720.6100.7600.2300.7406.460Diesel
Central Termical do Lobito2010100.0001.9101018.6405.00031.0002.7216.66410.6088.0404.1600.4900.510-2.420Diesel
Central Termica de Biopio201018.0001.8901026.61015.00046.0002.7507.04511.3404.4201.2000.3400.630-2.000Diesel
Central Termica do Cunene201034.0001.8901023.1805.00020.0002.7089.97317.2381.2104.4200.7500.7805.870Diesel
Central Termica Do Kuando Kubango201045.0001.8501053.04010.00016.0002.7218.12913.5361.0901.8800.9500.7202.540Diesel
Central Termica do Benfica201012.0001.8401063.62035.00034.0002.77110.35317.934-1.5800.6100.4700.8406.790Diesel
Central Termica do Cuito201045.0001.8201051.9207.00042.0002.7438.42114.0980.9602.6500.3700.7602.250Diesel
Central Termica do Moxico201045.0001.8001036.93030.00023.0002.7359.96717.1991.2000.6300.6700.9102.250Diesel
Central Termica de Cazenga201089.0001.7801033.33030.00050.0002.77343.89119.1730.4500.8300.3200.7708.380Combined Cycle Gas
Central Termica do Caminho de Ferro de Luanda201034.0001.7601020.79050.00065.0002.78810.39418.0001.1200.5000.2500.7207.280Diesel
Central Termica do Morro Bento201035.0001.7501021.59015.00032.0002.75210.20817.6642.4201.8400.4900.6409.470Diesel
Central Termica do Morro da Luz201016.0001.7301032.19010.00025.0002.7258.78014.8351.5002.1500.6100.6904.630Diesel
Central Termica KPM9201023.0001.7301033.34010.00035.0002.7117.69812.6854.3802.1500.4300.5901.940Diesel
Central Termica Benfica201038.0000.9401034.48017.00018.0002.70410.00817.3132.5601.3400.8300.7605.290Diesel
Central Termica Praia do Bispo201034.0001.7101035.63020.00015.0002.65710.72418.7921.9301.1600.9500.8106.890Diesel
Central Termica dos Quarteis201035.0001.6901033.46030.00037.0002.74310.50718.270-1.6700.7000.4200.8707.060Diesel
Central Termica Boavista I201023.0001.8601041.21015.00030.0002.7289.35115.9751.6101.4200.5100.7504.340Diesel
Central Termica Boavista II201024.0001.8001026.37010.00045.0002.7859.09415.4030.4102.1100.3600.7304.350Diesel
Central Termica Boavista III201012.0001.7701044.07016.00042.0002.7968.82814.8611.3101.2900.3900.7203.180Diesel
Central Termica Kassaki201017.0001.7501070.1808.00028.0002.7528.62814.5041.7602.5900.5600.7003.320Diesel
Central Termica CEEF201034.0001.7301158.13010.00020.0002.6817.79812.9155.0102.0500.7300.6300.930Diesel
Central Termica do Bengo201023.0001.7101220.50018.00025.0002.7259.16015.5950.9701.1400.6100.7604.370Diesel
Central Termica do Uige201024.0001.7501162.62015.00020.0002.7088.70314.6971.9201.4200.7500.6904.400Diesel
Central Termica Banza Congo2010100.0001.7101134.30010.00035.0002.7348.04513.3563.9902.1200.4400.6301.890Diesel
Central Termica Soyo201034.0001.9801573.73070.00015.0002.68829.33011.9775.0800.2900.9800.5900.680Combined Cycle Gas
Central Termica Landana201028.0001.9801212.84020.00030.0002.7088.03213.3563.5001.0600.5000.6302.620Diesel
Central Termica Malongo201045.0001.5601105.43040.00015.0002.6989.51116.3242.7500.5300.9900.7703.790Diesel

Appendix B. Complete Ranking of DMUs

DMUPhase I ψ o 2010 ψ o 2011 ψ o 2012 ψ o 2013 ψ o 2014 ψ o 2015 ψ o 2016 Phase IIRankDMU
Location
DMU
Name
DMU Capacity (MW)
11.08250.850.91690.9910.99041.17731.24531.4068-15NamibeXitoto Novo70
21.2261.05241.09151.13611.13451.34691.41041.4101-23NamibeXitoto Antigo8
31.11931.2121.18841.1551.15461.04441.0431.0376-16NamibeAeroporto8
40.95460.9510.95170.95190.95690.94180.96010.9691-8LubangoArimba35
51.28011.29811.28051.26271.26141.28191.27621.3001-29LubangoAnexo SE20
61.12080.83380.84270.85850.86321.49971.47071.4769-17BenguelaKileva30
70.77570.90640.88480.85210.85240.72320.72450.4868-4BenguelaLobito5
80.52080.4480.45140.45340.46220.37260.40351.0545-1BenguelaBiópio15
90.75470.75340.75390.75440.75470.7550.75550.756-3CuneneCunene5
100.81430.82740.81810.8060.80020.81080.81040.8275-5Cuando-CubangoKuando Kubango10
111.19041.05091.05341.0571.05571.35171.35661.4074-19HuamboBenfica35
121.32481.38251.38231.38271.38121.24931.24561.2496-31BiéCuito7
131.20691.25221.21491.17361.17421.18581.20141.2459-21MoxicoMoxico30
142.39542.58892.53062.45042.4532.45372.49371.7977-32LuandaCazenga30
151.25331.08511.07621.07791.07691.46171.46391.5313-26LuandaCaminho de Ferro de Luanda50
161.231.16081.15171.13621.13661.33311.3371.3543-24LuandaMorro Bento15
171.01050.89560.91490.93760.93951.12631.1211.1387-9LuandaMorro da Luz10
181.25661.43241.36481.30781.30841.1291.13111.1228-27LuandaKM910
191.01570.97420.97690.97860.98221.12771.03331.0371-10LuandaBenfica17
200.92660.93550.93640.93740.93990.84880.94350.9448-7LuandaPraia do Bispo20
211.29791.21241.21431.21351.21431.40041.40491.4252-30LuandaQuarteis30
221.20561.17731.181.1761.17531.21831.24741.2648-20LuandaBoavista I15
231.26931.15841.19721.24481.2451.351.34361.3463-28LuandaBoavista II10
241.2361.14931.16751.20841.20861.30931.30641.3022-25LuandaBoavista III16
251.06911.00571.01811.03351.03151.13121.13091.133-14LuandaKassaki8
261.03211.08181.07261.05131.05080.98910.98920.9897-12LuandaCEEF10
271.12481.09361.08761.08111.08221.16971.17431.1852-18BengoBengo18
281.02360.9730.9750.97720.97911.07961.091.0914-11UigeUige15
291.06570.95780.98681.02781.0281.16661.16571.1274-13ZaireBanza Congo10
300.60980.40690.44870.50.50550.50550.851.0524-2ZaireSoyo70
311.21691.34051.2961.2451.24381.11951.12521.1483-22CabindaLandana20
320.85120.91060.89220.87180.8660.92130.58170.9149-6CabindaMalongo40

References

  1. Eser, P.; Singh, A.; Chokani, N.; Abhari, R. Effect of increased renewables generation on operation of thermal power plants. Appl. Energy 2016, 164, 723–732. [Google Scholar] [CrossRef]
  2. Faraway, J.J. Extending the Linear Model with R: Generalized Linear, Mixed Effects and Nonparametric Regression Models; Chapman & Hall/CRC: Boca Raton, FL, USA, 2006. [Google Scholar]
  3. Chien, C.F.; Chen, W.C.; Lo, F.Y.; Lin, Y.C. A case study to evaluate the productivity change of the thermal power plants of the Taiwan power company. IEEE Trans. Energy Conserv. 2007, 22, 680–688. [Google Scholar] [CrossRef]
  4. Yao, X.; Guo, C.; Shao, S.; Jiang, Z. Total-factor CO2 emission performance of China’s provincial industrial sector: A meta-frontier non-radial Malmquist index approach. Appl. Energy 2016, 184, 1142–1153. [Google Scholar] [CrossRef]
  5. Wang, Y.-S.; Xie, B.-C.; Shang, L.-F.; Li, W.-H. Measures to improve the performance of China’s thermal power industry in view of cost efficiency. Appl. Energy 2013, 112, 1078–1086. [Google Scholar] [CrossRef]
  6. Lin, B.; Yang, L. Efficiency effect of changing investment structure on China’s power industry. Renew. Sustain. Energy Rev. 2014, 39, 403–411. [Google Scholar] [CrossRef]
  7. Barros, C.P.; Wanke, P. Efficiency in Angolan thermal power plants: Evidence from cost structure and pollutant emissions. Energy 2017, 130, 129–143. [Google Scholar] [CrossRef]
  8. Wanke, P.; Barros, C.P. Efficiency thresholds and cost structure in Senegal airports. J. Air Transp. Manag. 2017, 58, 100–112. [Google Scholar] [CrossRef]
  9. Park, S.U.; Lesourd, J.B. The efficiency of conventional fuel power plants in South Korea: A comparison of parametric and non-parametric approaches. Int. J. Prod. Econ. 2000, 63, 59–67. [Google Scholar] [CrossRef]
  10. Lam, P.L.; Shiu, A. A data envelopment analysis of the efficiency of China’s thermal power generation. Util. Policy 2001, 10, 75–83. [Google Scholar] [CrossRef]
  11. Sarica, K.; Or, I. Efficiency assessment of Turkish power plants using data envelopment analysis. Energy 2007, 32, 1484–1499. [Google Scholar] [CrossRef]
  12. Barros, C.P.; Peypoch, N. Technical efficiency of thermoelectric power plants. Energy Econ. 2008, 30, 3118–3127. [Google Scholar] [CrossRef]
  13. Nakano, M.; Managi, S. Regulatory reforms and productivity: An empirical analysis of the Japanese electricity industry. Energy Policy 2008, 36, 201–209. [Google Scholar] [CrossRef]
  14. Sozen, A.; Alp, I.; Özdemir, A. Assessment of operational and environmental performance of the thermal power plants in Turkey by using data envelopment analysis. Energy Policy 2010, 38, 6194–6203. [Google Scholar] [CrossRef]
  15. Liu, C.H.; Lin, S.J.; Lewis, C. Evaluation of thermal power plant operational performance in Taiwan by data envelopment analysis. Energy Policy 2010, 38, 1049–1058. [Google Scholar] [CrossRef]
  16. Rezaee, M.J.; Moini, A.; Makui, A. Operational and non-operational performance evaluation of thermal power plants in Iran: A game theory approach. Energy 2012, 38, 96–103. [Google Scholar] [CrossRef]
  17. Rezaee, M.J. Using Shapley value in multi-objective data envelopment analysis: Power plants evaluation with multiple frontiers. Int. J. Electr. Power Energy Syst. 2015, 69, 141–149. [Google Scholar] [CrossRef]
  18. Shrivastava, N.; Sharma, S.; Chauhan, K. Efficiency assessment and benchmarking of thermal power plants in India. Energy Policy 2012, 40, 159–176. [Google Scholar] [CrossRef]
  19. Sueyoshi, T.; Goto, M. DEA environmental assessment of coal fired power plants: Methodological comparison between radial and non-radial models. Energy Econ. 2012, 34, 1854–1863. [Google Scholar] [CrossRef]
  20. Du, L.; Mao, J. Estimating the environmental efficiency and marginal CO2 abatement cost of coal-fired power plants in China. Energy Policy 2015, 85, 347–356. [Google Scholar] [CrossRef]
  21. Ghosh, R.; Kathuria, V. The effect of regulatory governance on efficiency of thermal power generation in India: A stochastic frontier analysis. Energy Policy 2016, 89, 11–24. [Google Scholar] [CrossRef]
  22. Yan, D.; Lei, Y.; Li, L.; Song, W. Carbon emission efficiency and spatial clustering analyses in China’s thermal power industry: Evidence from the provincial level. J. Clean. Prod. 2017, 156, 518–527. [Google Scholar] [CrossRef]
  23. Wu, Y.; Ke, Y.; Xu, C.; Xiao, X.; Hu, Y. Eco-efficiency measurement of coal-fired power plants in China using super efficiency data envelopment analysis. Sustain. Cities Soc. 2018, 36, 157–168. [Google Scholar] [CrossRef]
  24. See, K.F.; Coelli, T. An analysis of factors that influence the technical efficiency of Malaysian thermal power plants. Energy Econ. 2012, 34, 677–685. [Google Scholar] [CrossRef]
  25. Munisamy, S.; Arabi, B. Eco-efficiency change in power plants: Using a slacks-based measure for the meta-frontier Malmquist–Luenberger productivity index. J. Clean Prod. 2015, 105, 218–232. [Google Scholar] [CrossRef]
  26. Mahmoudi, R.; Emrouznejad, A.; Khosroshahi, H.; Khashei, M.; Rajabi, P. Performance evaluation of thermal power plants considering CO2 emission: A multistage PCA, clustering, game theory and data envelopment analysis. J. Clean. Prod. 2019, 223, 641–650. [Google Scholar] [CrossRef] [Green Version]
  27. Wei, X.; Zhang, N. The shadow prices of CO2 and SO2 for Chinese Coal-fired Power Plants: A partial frontier approach. Energy Econ. 2020, 85, 104576. [Google Scholar] [CrossRef]
  28. Emrouznejad, A.; Yang, G. A survey and analysis of the first 40 years of scholarly literature in DEA: 1978–2016. Socio-Econ. Plan. Sci. 2018, 61, 4–8. [Google Scholar] [CrossRef]
  29. Fernández, D.; Pozo, C.; Folgado, R.; Jiménez, L.; Guillén-Gosálbez, G. Productivity and energy efficiency assessment of existing industrial gases facilities via data envelopment analysis and the Malmquist index. Appl. Energy 2018, 212, 1563–1577. [Google Scholar] [CrossRef]
  30. Charnes, A.; Clark, C.T.; Cooper, W.W.; Golany, B. A developmental study of data envelopment analysis in measuring the efficiency of maintenance units in the U.S. air forces. Ann. Oper. Res. 1985, 2, 95–112. [Google Scholar] [CrossRef]
  31. Färe, R.; Grosskopf, S. Intertemporal Production Functions with Dynamic DEA; Kluwer Academic Publishers: Boston, MA, USA, 1996. [Google Scholar]
  32. Nemoto, J.; Goto, M. Dynamic data envelopment analysis: modeling intertemporal behavior of a firm in the presence of productive inefficiencies. Econ. Lett. 1999, 64, 51–56. [Google Scholar] [CrossRef]
  33. Thompson, R.G.; Dharmapala, P.S.; Thrall, R.M. Linked-cone DEA profit ratios and technical efficiency with application to Illinois coal mines. Int. J. Prod. Econ. 1995, 39, 99–115. [Google Scholar] [CrossRef]
  34. Jahanshahloo, G.R.; Soleimani-damaneh, M.; Reshadi, M. On the Pareto (dynamically) efficient paths. Int. J. Comput. Math. 2006, 83, 629–633. [Google Scholar] [CrossRef]
  35. Silva, E.; Stefanou, S.E. Nonparametric dynamic efficiency measurement: Theory and application. Am. J. Agric. Econ. 2007, 89, 398–419. [Google Scholar] [CrossRef]
  36. Soleimani-damaneh, M. An effective computational attempt in DDEA. Appl. Math. Model. 2009, 33, 3943–3948. [Google Scholar] [CrossRef]
  37. Sam Park, K.; Park, K. Measurement of multi-period aggregative efficiency. Eur. J Oper. Res. 2009, 193, 567–580. [Google Scholar] [CrossRef]
  38. Charnes, A.; Cooper, W.; Rhodes, E. Measuring the efficiency of decision making units. Eur. J. Oper. Res. 1978, 2, 429–444. [Google Scholar] [CrossRef]
  39. Andersen, P.; Petersen, N.C. A procedure for ranking efficient units in data envelopment analysis. Manag. Sci. 1993, 39, 1261–1264. [Google Scholar] [CrossRef]
  40. Mariano, E.B.; Gobbo, J.A., Jr.; Camioto, F.C.; Rebelatto, D.A.N. CO2 emissions and logistics performance: A composite index proposal. J. Clean Prod. 2016, 163, 166–178. [Google Scholar] [CrossRef]
  41. Ozkan, N.F.; Ulutas, B.H. Efficiency analysis of cement manufacturing facilities in Turkey considering undesirable outputs. J. Clean Prod. 2017, 156, 932–938. [Google Scholar] [CrossRef]
  42. Scheel, H. Undesirable outputs in efficiency valuations. Eur. J. Oper. Res. 2001, 132, 400–410. [Google Scholar] [CrossRef]
  43. Seiford, L.M.; Zhu, J. Modeling undesirable factors in efficiency evaluation. Eur. J. Oper. Res. 2002, 142, 16–20. [Google Scholar] [CrossRef]
  44. Hadi-Vencheh, A.; Kazemi Matin, R.; Tavassoli Kajani, M. Undesirable factors in efficiency measurement. Appl. Math. Comput. 2005, 63, 547–552. [Google Scholar] [CrossRef]
  45. EIA. How Much Does It Cost to Generate Electricity with Different Types of Power Plants? 2017. Available online: https://www.eia.gov/tools/faqs/faq.php?id=19&t=3 (accessed on 1 March 2018).
  46. EIA. How Much Does It Cost to Build Different Types of Power Plants? 2018. Available online: https://www.eia.gov/tools/faqs/faq.php?id=487&t=3 (accessed on 7 July 2018).
  47. Morgan, J.; Sonquist, J. Problems in the analysis of survey data and a proposal. J. Am. Stat. Assoc. 1963, 58, 415–434. [Google Scholar] [CrossRef]
  48. Breiman, L.; Friedman, J.; Olshen, R.; Stone, C. Classification and Regression Trees; Wadsworth: Belmont, CA, USA, 1984. [Google Scholar]
  49. James, G.; Witten, D.; Hastie, T.; Tibshirani, R. An Introduction to Statistical Learning with Applications in R; Springer: New York, NY, USA, 2013. [Google Scholar]
  50. Opitz, D.; Maclin, R. Popular ensemble methods: An empirical study. J. Artif. Intell. Res. 1999, 11, 169–198. [Google Scholar] [CrossRef]
  51. Polikar, R. Ensemble based systems in decision making. IEEE Circuits Syst. Mag. 2006, 6, 21–45. [Google Scholar] [CrossRef]
  52. Torgo, L. Data Mining with R: Learning with Case Studies; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
  53. Breiman, L. Bagging predictors. Mach. Learn. 1996, 24, 123–140. [Google Scholar] [CrossRef] [Green Version]
  54. Ledolter, J. Data Mining and Business Analytics with R; Wiley: Hoboken, NJ, USA, 2013. [Google Scholar]
Figure 1. Pooled super-efficiency scores boxplot.
Figure 1. Pooled super-efficiency scores boxplot.
Energies 13 02404 g001
Figure 2. Distribution of the super efficiency scores per year.
Figure 2. Distribution of the super efficiency scores per year.
Energies 13 02404 g002
Figure 3. Scatter plots for super-efficiency levels and pollutant emissions.
Figure 3. Scatter plots for super-efficiency levels and pollutant emissions.
Energies 13 02404 g003
Figure 4. Scatter plots for super-efficiency levels and cost structure variables.
Figure 4. Scatter plots for super-efficiency levels and cost structure variables.
Energies 13 02404 g004
Figure 5. Bootstrapped regression tree thresholds.
Figure 5. Bootstrapped regression tree thresholds.
Energies 13 02404 g005
Figure 6. Relative importance of the cost structure variables for efficiency generation in thermal plants.
Figure 6. Relative importance of the cost structure variables for efficiency generation in thermal plants.
Energies 13 02404 g006
Table 1. Literature review synthesis.
Table 1. Literature review synthesis.
Author (Year)CountrySample SizeMethods
Park and Lesourd [9]South Korea64 fuel power plantsDEA and Stochastic-Frontier
Lam and Shiu [10] ChinaThermal Power Generation in 30 provinces, autonomous regions, and municipalitiesDEA and Tobit Regression
Chien et al. [3]Taiwan8 thermal power plantsDEA–Malmquist
Sarica and Or [3] Turkey65 thermal, hydro, and wind power plantsConstant returns to scale (CRS) DEA, variable returns to scale (VRS) DEA, and Assurance region DEA
Barros and Peypoch [11] Portugal7 thermoelectric plants data from 1996 to 2004DEA and Simar and Wilson bootstrapped procedure
Nakano and Managi [13] Japan10 local monopoly companies data from 1965 to 2003DEA, Luenberger productivity indicator and Generalized Method of Moments (GMM) estimation
Sozen et al. [14] Turkey15 thermal power plantsCRS DEA and VRS DEA
Liu et al. [15]Taiwan9 thermal power plantsDEA
Rezaee et al. [16]Iran24 power plantsDEA and game theory
See and Coelli [24]Malaysia14 thermal plantsStochastic Frontier Analysis (SFA)
Shrivastava et al. [18]India60 coal-fired power plantsCRS (Constant Returns-to-Scale) DEA and VRS (Variable Returns-to-Scale) DEA
Sueyoshi and Goto [19] United States20 U.S. coal-fired power plantsNon-radial DEA model
Wang et al. [5]China30 thermal powerDEA, Malmquist–Luenberger productivity index
Lin and Yang [6]ChinaPower industry for 31 provinces from 2005 to 2010Slacks-Based Measure (SBM) in dynamic DEA model
Rezaee [17]Iran20 power plantsShapley value and multiobjective DEA
Du and Mao [20]China1158 power plants dataParametric Linear Programing
Munisamy and Arabi [25]Iran48 thermal power plantsSBM, Malmquist–Luenberger index
Ghosh and Kathuria [21]India77 coal-based thermal power plantsSFA
Barros and Wanke [7]Angola32 Angolan thermal power plants from 2010 to 2014SBM-Undesirable and Beta regression
Yan et al. [22]ChinaPower industry in 30 provincesSBM-Undesirable, Malmquist index
Wu et al. [23]China58 coal-fired power plantsSuper Efficiency DEA, Kruskal–Wallis rank, Tobit regression
Xie et al. [2]ChinaThermal power plants in 30 provinces Nonparametric weighted Russell directional distance method
Mahmoudi et al. [26]Iran52 thermal power plantsMultivariate data analysis techniques, game theory, and Shannon entropy combined with DEA
Wei and Zhang [27]China93 coal-fired power plantsPartial parametric environmental production frontier
Table 2. Descriptive statistics for the sample.
Table 2. Descriptive statistics for the sample.
VariablesMinMaxMeanSD
InputsPlant capacity in MW57021.53316.663
Investment costs * (log)10.60821.94515.68452.3792
Employee costs per year * (log)2.6572.8012.736810.0332
Fuel costs per year * (log)6.66443.89111.02357.3653
Number of employees156931.61313.671
OutputsProduction in 1000 MWh1011735.29224.206
Carbon dioxide emissions in tons per year0.9402.7201.9890.251
Discharge of polluted water in liters per year1,018.641,898.211,190.07136.460
TrendTrend173.9872.005
Trend squared14919.91616.407
Cost StructureCapital–labor ratio−1.70011.4082.2571.923
Capacity cost per MW0.2904.4301.5710.959
Labor cost per employee0.2300.9960.5760.226
Capital cost0.51010.7210.089
Cost–asset ratio−2.42012.0584.2372.869
* Investment costs, employee costs per year, and fuel costs per year were originally measured in current Kwanzas and subsequently converted to 2016 USD.
Table 3. R packages and respective default values for bootstrapped regression trees.
Table 3. R packages and respective default values for bootstrapped regression trees.
MethodR LibraryR FunctionValues
RF-CARTrandomForestrandomForestx, where x represents a data frame or matrix of predictors,
y = NULL,
xtest = NULL,
ytest = NULL,
ntree = 500,
mtry = max(floor(ncol(x)/3), 1),
replace = TRUE,
classwt = NULL,
sampsize = nrow(x),
nodesize = 5,
maxnodes = NULL,
importance = FALSE,
localImp = FALSE,
nPerm = 1,
norm.votes = TRUE,
do.trace = FALSE,
keep.forest = TRUE,
corr.bias = FALSE,
keep.inbag = FALSE
CARTtreeTreena.action = na.pass,
method = "recursive.partition",
split = c("deviance", "gini"),
model = FALSE,
x = FALSE,
y = TRUE,
wts = TRUE,
mincut = 5,
minsize = 10,
mindev = 0.01

Share and Cite

MDPI and ACS Style

Hadi-Vencheh, A.; Wanke, P.; Jamshidi, A. What Does Cost Structure Have to Say about Thermal Plant Energy Efficiency? The Case from Angola. Energies 2020, 13, 2404. https://doi.org/10.3390/en13092404

AMA Style

Hadi-Vencheh A, Wanke P, Jamshidi A. What Does Cost Structure Have to Say about Thermal Plant Energy Efficiency? The Case from Angola. Energies. 2020; 13(9):2404. https://doi.org/10.3390/en13092404

Chicago/Turabian Style

Hadi-Vencheh, Abdollah, Peter Wanke, and Ali Jamshidi. 2020. "What Does Cost Structure Have to Say about Thermal Plant Energy Efficiency? The Case from Angola" Energies 13, no. 9: 2404. https://doi.org/10.3390/en13092404

APA Style

Hadi-Vencheh, A., Wanke, P., & Jamshidi, A. (2020). What Does Cost Structure Have to Say about Thermal Plant Energy Efficiency? The Case from Angola. Energies, 13(9), 2404. https://doi.org/10.3390/en13092404

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