Loading [MathJax]/jax/element/mml/optable/GeneralPunctuation.js
Research article

A DEA model to sustainability improvement of the electricity supply chain in presence dual-role factors and undesirable outputs: A case on the power industry

  • Received: 28 March 2020 Accepted: 12 June 2020 Published: 09 July 2020
  • The energy and industrial sectors are the most attractive investment regions for enhancing efficiency in production processes. The power industry is one of the important investment targets for enhancing corporate sustainability. One of the most fundamental problems in the power industry is the control of wasted energy in oil and gas fields and power plant sectors and the power losses management in transmission and distribution lines. The investment to new technology innovation and environmental protection from pollution gases emission in energy and power plant sectors and the power losses management in transmission and distribution lines play an important role in the implementation progress of the power industry. The purpose of this study is to examine the effects of investment to flare gas and greenhouses gases reduction in energy and power plant sections and power losses control by equipping sections to improved engineering systems in transmission and distribution networks of the electricity supply chain. Indeed, the supply chain management needs information related to investment effect to activity level control as handling flare gas in energy sections and reducing harmful substance emissions and greenhouses gases in power plant sectors and harnessing power losses in transmission and distribution networks. The proposed approach evaluates the sustainability and efficiency of an electricity supply chain by a radial model in the presence of two categories of inputs under natural and managerial disposability, dual-role factors and undesirable produces. A real case on the Iran power industry is presented to demonstrate the applicability and practicability of the proposed method. Moreover, to demonstrate the capability of the proposed approach a supply chain identified by oil and gas companies, power plants, transmissions companies, dispatching companies and final consumers in the Iran power industry. One empirical implication has obtained from model performance in the electricity supply chain. The results indicate approximately, the oil and gas fields, the power plants and the distribution lines and the public divisions of power consumers have earned 100%, 90% and 90% efficiency of the total in supply chains, respectively. Particularly, this study recommends that transmission and distribution companies must have adequate decisional capacities regarding investment for transmitting power to industrial, agriculture divisions in the power industry.

    Citation: Mojgan Pouralizadeh. A DEA model to sustainability improvement of the electricity supply chain in presence dual-role factors and undesirable outputs: A case on the power industry[J]. AIMS Energy, 2020, 8(4): 580-614. doi: 10.3934/energy.2020.4.580

    Related Papers:

    [1] Mojgan Pouralizadeh, Aliraza Amirtaimoori, Rossana Riccardi, Mohsen Vaez-Ghasemi . Supply chain performance evaluation in the presence of undesirable products: A case on power industry. AIMS Energy, 2020, 8(1): 48-80. doi: 10.3934/energy.2020.1.48
    [2] Ladislas Mutunda Kangaji, Atanda Raji, Efe Orumwense . Cutting-edge progress in offshore wind and tidal stream power technology—State-of-the-Art. AIMS Energy, 2025, 13(1): 188-230. doi: 10.3934/energy.2025007
    [3] Chukwuebuka Okafor, Christian Madu, Charles Ajaero, Juliet Ibekwe, Festus Otunomo . Situating coupled circular economy and energy transition in an emerging economy. AIMS Energy, 2021, 9(4): 651-675. doi: 10.3934/energy.2021031
    [4] Surender Reddy Salkuti . Sustainable energy technologies for emerging renewable energy and electric vehicles. AIMS Energy, 2024, 12(6): 1264-1270. doi: 10.3934/energy.2024057
    [5] Kaovinath Appalasamy, R Mamat, Sudhakar Kumarasamy . Smart thermal management of photovoltaic systems: Innovative strategies. AIMS Energy, 2025, 13(2): 309-353. doi: 10.3934/energy.2025013
    [6] Mowffaq M. Oreijah, Mohammed Yunus . A parametric analysis to evaluate the performance metrics of power generation system involving Trilateral Flash Cycle using three different working fluids for low grade waste heat. AIMS Energy, 2019, 7(4): 483-492. doi: 10.3934/energy.2019.4.483
    [7] D. Marene Larruskain, Inmaculada Zamora, Oihane Abarrategui, Garikoitz Buigues, Víctor Valverde, Araitz Iturregi . Adapting AC Lines to DC Grids for Large-Scale Renewable Power Transmission. AIMS Energy, 2014, 2(4): 385-398. doi: 10.3934/energy.2014.4.385
    [8] Arben Gjukaj, Rexhep Shaqiri, Qamil Kabashi, Vezir Rexhepi . Renewable energy integration and distributed generation in Kosovo: Challenges and solutions for enhanced energy quality. AIMS Energy, 2024, 12(3): 686-705. doi: 10.3934/energy.2024032
    [9] Seif Bayoumi, Mohamed E.A. Ali, Islam Amin, Raneem El Torky, Selda Oterkus, Hosam Shawky, Erkan Oterkus . Environmentally-driven design of a floating desalination platform (Case study: reverse osmosis floating desalination platform of ras gharib, Egypt). AIMS Energy, 2021, 9(3): 623-650. doi: 10.3934/energy.2021030
    [10] Sebastian Auburger, Richard Wüstholz, Eckart Petig, Enno Bahrs . Biogas digestate and its economic impact on farms and biogas plants according to the upper limit for nitrogen spreading—the case of nutrient-burdened areas in north-west Germany. AIMS Energy, 2015, 3(4): 740-759. doi: 10.3934/energy.2015.4.740
  • The energy and industrial sectors are the most attractive investment regions for enhancing efficiency in production processes. The power industry is one of the important investment targets for enhancing corporate sustainability. One of the most fundamental problems in the power industry is the control of wasted energy in oil and gas fields and power plant sectors and the power losses management in transmission and distribution lines. The investment to new technology innovation and environmental protection from pollution gases emission in energy and power plant sectors and the power losses management in transmission and distribution lines play an important role in the implementation progress of the power industry. The purpose of this study is to examine the effects of investment to flare gas and greenhouses gases reduction in energy and power plant sections and power losses control by equipping sections to improved engineering systems in transmission and distribution networks of the electricity supply chain. Indeed, the supply chain management needs information related to investment effect to activity level control as handling flare gas in energy sections and reducing harmful substance emissions and greenhouses gases in power plant sectors and harnessing power losses in transmission and distribution networks. The proposed approach evaluates the sustainability and efficiency of an electricity supply chain by a radial model in the presence of two categories of inputs under natural and managerial disposability, dual-role factors and undesirable produces. A real case on the Iran power industry is presented to demonstrate the applicability and practicability of the proposed method. Moreover, to demonstrate the capability of the proposed approach a supply chain identified by oil and gas companies, power plants, transmissions companies, dispatching companies and final consumers in the Iran power industry. One empirical implication has obtained from model performance in the electricity supply chain. The results indicate approximately, the oil and gas fields, the power plants and the distribution lines and the public divisions of power consumers have earned 100%, 90% and 90% efficiency of the total in supply chains, respectively. Particularly, this study recommends that transmission and distribution companies must have adequate decisional capacities regarding investment for transmitting power to industrial, agriculture divisions in the power industry.


    An electricity supply chain is a network of suppliers, producers, transmitters and distributors in which raw materials are transformed into final products and delivered to the customers. The Energy sector is one of the most important types of developed infrastructures in any country. The fossil fuels are energy sources of incompatible with the environment so that they emissions various pollutions and greenhouse gases in economics activities. The flare gas emission is one of the most critical problems in oil and gas fields. According to published statistics in the year 2017, daily 4 million oil barrels produced in Iran oil fields and about 45 million cubic meters associated gas (gas in oil) have been burned to avoid from the possible explosion in oil and gas fields that burning fossil fuels not only a big thread for human health and the other organisms but also cause decrease economic return in industrial activities. The repairs problems respect to strengthening pressure systems of associated gas (gas in oil) lead burning the large quantities of gas in oil fields. The problems such as the lack of enough education of the workforce or the lack of pieces timely preparation cause wasting the amount considerable energy in the energy industry. Also, the burned gases release more than 250 toxic substances in the air. Moreover, the associated gas can be used in Liquid petrol gas (LPG) production for car and urban consumptions, power production and chemical and petrochemical derivatives production. Besides, the associate gas can be injected into oil reservoirs to the rehabilitation of thanks and prevention from drying oil reservoirs. Similarly, Power plants from production to consumption produce kinds of contaminations in the environment. Power plants are the largest fossil fuel consumers such as coal, fuel oil and gasoline and natural gas. The sixty–eight percent of the Iran power plants are non-renewable and they consume fossil fuels to power production. These fuels have been playing most of the key role in electricity production and release a huge amount of pollution substance in the production process. Hence, this is immediately necessary to enhance efficiency by the protection of the negative impacts of economic activities. The carbon dioxide gas (CO2) has the most contributions to pollution emissions in power plants. This gas cause climate changes and global warming also it is a threat to human health and other organisms. Therefore, we must reduce the number of greenhouse gases (GHG) by enhancing systems efficiently, otherwise, we will confront sever events such as heat waves, droughts, floods and other harmful factors to social and economics. Approximately, one percent of power plants' nominal capacity is devoted to power losses in transmission and distribution lines. The one percent of power plant capacity is equivalent to 2.5 billion kilowatt-hour power that to produce this amount of electricity releases about 1.8 ton Carbon dioxide (CO2) in air. Therefore, a profit solution to this problem is new ideas performance to investment opportunities and Technology innovation to harmful effects protection of environmental. In other words, if supply chain enterprises equipped with improved engineering capability and invest in improvement and repair of equipment in the divisions then undesirable outputs considerably decrease in production activities. Therefore, supply chain management should enable propose an appropriate approach to wasted energy harness to environmental efficiency enhancement in the power industry. Besides, power plants produce power also they need electricity to power regulations. The inner electricity consumption of power plants divided into technical and non-technical consumptions. The power station generators plant voltage regulators to control the output of power plants in an electric power distribution system. The voltages regulators install at power plants to power transmission with steady voltage and they may be installed at distribution lines to customers receive steady voltage. Therefore, the control of the electricity consumption of power plants can significantly enhance unified efficiency (operational and environmental) in power stations. Data envelopment analysis (DEA) is a profitable method to performance new ideas to investment opportunities and Technology innovation to harmful effects protection of environmental. Let us now suppose that supply chain divisions apply inputs to produce desirable and undesirable outputs as the inputs separate into two categories under natural and managerial disposability. Also, Let us consider undesirable outputs such as emissions of harmful substances in the air, water and ground and other detrimental variables of production activities. Besides, certain factors are considered to simultaneously play the role of both inputs and outputs in production processes. These factors are called dual-role factors. Also, the material flow is transferred from suppliers to manufacturers and from manufacturers to transmitters and from them to distributors and finally from distributors to customers in the production processes. Also, intermediate measures flow between divisions of two consecutive steps in the two inverse directions. Furthermore, the inverse intermediate measures exit from transmitter divisions and enter to manufacture divisions and exit from manufacture divisions and enter to supplier divisions.

    In this study, we are going to answer the following questions: how a decision-making unit or a supply chain enables decrease pollution gases emissions by investment on specialist workforce and flare gas recovery systems in oil and gas fields and the new technology innovation in non-renewable power plants and handling wasted energy by engineer workforce as power losses noticeable abatement in transmission and distribution lines? In this case, supply chain management should be able to identify whether inputs increase under managerial disposability to new technology innovation reduce undesirable productions in the electricity supply chain divisions or the increase inputs for investment ineffective for decrease a number of undesirable outputs. Also, it is immediately necessary to know whether the investment can effectively decrease the amount of undesirable outputs or increase the inputs under managerial disposability have a limited effect on decrease an amount of undesirable outputs. Moreover, the supply chain management needs information related to investment effect to the inputs level control under managerial disposability as handling flare gas in energy sections and reducing pollution emissions and greenhouses gases in power plant sectors and harnessing energy wasted in transmission and distribution networks. Furthermore, how the factors to simultaneously play the role of both inputs and outputs can be applied to the costs flare gas control in energy sections and the inner electricity consumption (technical and non-technical) management in power plant sectors and wasted energy harness in distribution lines. In this study, managerial disposability is accomplished by investment into the Energy section to flare gas reduction and environmental protection, construction and initiation of renewable power plants to pollution emissions prevention in the power plant section. Meanwhile, transmission and distribution lines are equipped with improved engineering capability to power losses reduction. Also, dual-role factors control the cost recovery of flare gas in energy sections and the inner electricity consumption of power plants (technical and non-technical) and increase the scientific level of staff to the power losses harness in transmission operation. In current paper applied two concepts of natural and managerial disposability to environmental assessment as the inputs separate into two categories under natural and managerial disposability. Furthermore, we consider natural and managerial disposability to inputs and free disposability of undesirable outputs and weak disposability of desirable outputs so that we calculate a supply chain inefficiency score in the presence of two categories of inputs, dual-role factors, and desirable and undesirable outputs. In the more details, the divisions of every pair of members belong to a consecutive stage are connected by two sets of intermediate measures in the two inverse ways. Moreover, the inverse intermediate measures enter to divisions are considered as non-discretionary inputs. To include the two concepts of natural and managerial disposability to operational and environmental assessment Fan et al. [1] proposed a radial model based on data envelopment analysis to study on eco-efficiency of industrial parks in china. Sueyoshi et al. [2] presented an environmental assessment on Energy and sustainability by data envelopment analysis. Wang et al [3] are calculated operational and environmental efficiency in China' thermal power industry by a global fractional model as taking effectiveness measure as a complement to an efficiency measure. Zhang et al. [4] proposed a three-stage model based on data envelopment analysis. They calculated industrial eco-efficiency of 30 provinces in china. Moreover, the other research studies are presented to the management of greenhouse gases emissions in the production chain [5] and the control of renewable energy [6] and energy management in hybrid electrical vehicle [7].

    The remainder of this paper is organized as follows: In Section 2, we present an appropriate literature review on how DEA has been used for research on investment opportunities and technology innovation. Also, it is indicated, the literature summary on the presence of a dual-factors role in Data Envelopment Analysis. Moreover, we present a DEA model for resource utilization and investment in technology innovation. We show how correctly specify natural and managerial disposability in a production processes model of supply chain performance evaluation problems. Section 3 is devoted to introducing a procedure to calculate supply chain efficiency in the presence of two categories of inputs, undesirable products, and dual-role factors and the two set intermediate measures. In Section 4, we present a case study to demonstrate the applicability of the proposed method to the Iran power industry. In Section 5, we present our conclusions.

    In the following subsections, various studies on Environmental and operational assessment and green supply chain management (GSCM), and dual-role factors are briefly summarized.

    To include the two concepts of natural and managerial disposability to environmental assessment in the technology and account for the harmful substances prevention and negative impact on productivity Sueyoshi and Golver [8] discussed the history of DEA from the contributions of Cooper who first invented DEA in the 19 century.

    Gotto et al. [9,10] proposed a description of the conventional uses of DEA for environmental assessment. Then the concept of natural and managerial disposability has applied as a conceptual basis for preceding research efforts, see for example [11].

    Sueyoshi et al. [12] proposed a stage DEA model to operational and environmental assessment of Japanese industrial sectors. They calculated a unified efficiency score under natural and managerial disposability of the decision-making unit by resource utilization and technology innovation.

    Kao [13] modified the conventional DEA model by taking into account the series relationship of the two sub-processes within the whole process. Ton and Tustusi [14] proposed a slacks-based network DEA model called network SBM.

    Khalili et al. [15] proposed the fuzzy model for measurement of efficiency in transformation process of supply chain agility.

    Toloo et al. [16] proposed the DEA approach with mixed integer programming model to determine the most efficiency supplier without imprecise data.

    Tavana et al. [17] extended the EBM model proposed by Ton et al. [14] and proposed a new Network EBM (NEMB).

    Mahdiloo et al. [18] used the DEA model and DMUs by employing better integration of environmental and technical efficiency objective. They measured environmental, technical and eco-efficiency for supplier selection. The researchers have showed that all previous models are computationally cannot measure eco-efficiency in the best way. They Proposed the new model provide a valid eco-efficiency indicator of DMUs by utilizing a better combination of the technical and environmental efficiency.

    Tajbakhsh et al. [19] proposed a multi-stage data envelopment analysis model to evaluate the sustainability of a chain of business partners. They assess supply chain sustainability in the banking sector and beverage case.

    Khodakerami et al. [20] proposed the DEA new two stages model of supply chain sustainability in resin producing companies. The authors considered performance measurement of some imprecise and uncertain problems related to in real life as this problem needs to use fuzzy set in DEA model.

    Devika et al. [21] applied DEA approach for measurement of the pareto frontier quality. They have considered the social impact with economics and environmental impacts on class producer, simultaneously.

    Nikfarjam et al. [22] propose the new method DEA for measuring the supply chain with integrated to approaches. They showed the proposed model can use for evaluating of performance for identify the benchmarking units for inefficiency supply chain.

    Babazadeh et al. [23] used DEA approach to evaluate the social and climate criteria in cultivation areas. They evaluated strategic design of biodiesel supply chain network by integration of DEA and mathematical programming. Besides, the authors believe there is lack in previous studies which did not focus on climatic and social criteria and proposed a new DEA model related to biodiesel supply chain planning.

    Pouralizadeh et al. [24] proposed a new DEA-based model to sustainability evaluate an electricity supply chain in presence undesirable outputs. They planned a supply chain by five stages and fifteen divisions from different districts in Iran. Also, the weak disposability assumption was adopted for activity level control in production activity. The proposed model enable determents the type and size of inputs to control undesirable outputs.

    Toolo [25] proposed a revision of proposed model in [19]. Hatefi et al. [26] proposed a new model based on distance function for classifying inputs and outputs.

    Farzipoor [27] proposed a model for selecting third–party reverse logistics providers in the presence of multiple dual role factors and proposed [28] a model for selecting 3PL providers in the presence of both dual-role factors and imprecise data. All of the references mentioned in this subsection do not used network DEA model for GSCM evaluation problem.

    Mirhedayrian et al. [29] presented a DEA-based model in the presence of undesirable outputs, dual-role factors, and fuzzy data to a supply chain. They indicated a method to improve environmental performance a green supply chain management and incorporate dual-role factor and undesirable output into (NSBM) model proposed by Tone and Tsutsui [14].

    In summary, all of the abovementioned references for environmental performance assessment of the supply chain do not consider network DEA model based on the new technology innovation and targeting investment for the reduction of undesirable products. Also, the aforementioned models to sustainability assessment of supply chain are not able to determine whether the investment effectively decrease the number of undesirable outputs or limited effect on decreasing an amount of undesirable outputs. In other word, the investment may be ineffective for some of supply chain divisions to undesirable outputs abatement.

    In this Section are reported fundamental concepts for environmental and operational assessment decision-maker unit and the approach to calculate the unified efficiency (operational and environmental) of the electricity supply chain.

    Let us suppose Xj=(x1j,x2j,...,xmj)T>0, Gj=(g1j,g2j,...,gsj)T>0,Bj=(b1,b2,...,bhj)T>0 presents column vectors of inputs, desirable and undesirable outputs in jth DMU (Decision maker unit), respectively. Sueyoshi and Gotto [5] have proposed a radial model to measure the unified efficiency (operational and environmental) of the kth DMU under natural and managerial disposability of inputs as follows.

    maxξ+ε[Rxidxi+Rxqdxq+Rxfdxf]nj=1x+ijλj+dxi=xiki=1,...,mnj=1x+iqλjdxq=x+qkq=1,...,m+nj=1grjλj+ξgrk=grkr=1,...,snj=1bfjλjdbf=bfkf=1,...,hnj=1λj=1λj0,j=1,...,n,ξURS,dxi0,i=1,...,mdxq0,q=1,...,m+,dbf0,f=1,...,h (1)

    In this model, the number of original m inputs are separated into two categories m (under natural disposability) and m+ (under managerial disposability), respectively. The model maintains m=m+m+. Also, ξ is an inefficiency score that measures the distance between efficiency frontier and one observed vector of the desirable outputs and dxi,dxq,dbf are slack variables belong to two categories input and undesirable output, respectively. In addition, ε is a small amount and it considered as 0.0001 for our computation convenience it is possible for model to use ε=0 in model (5). In this model Rxi,Rxq,Rbf are specified by decision maker as fallows.

    Rxi=(m+s+h)1(max{xij|j=1,...,n}min{xij|j=1,...,n})1Rxq=(m+s+h)1(max{xqj|j=1,...,n}min{xqj|j=1,...,n})1Rbf=(m+s+h)1(max{bfj|j=1,...,n}min{bfj|j=1,...,n})1 (2)

    The column vectors of structural variables (λ) are applied for connecting the input and output vectors by convex combination under variable return scale. A unified efficiency score under natural and managerial disposability is measured as follows

    UEMN=1[ξ+ε(mi=1Rxidxi+m+q=1Rxqdxq+hf=1Rbfdbf)] (3)

    where the inefficiency score and all slack variables are determined on the optimality of Model (1).

    The weak disposability concept has specified on two outputs vectors of hth division, (Gh,Bh) as follows:

    Phw(x)={(Gh,Bh):Ghnj=1Ghjλhj,Bh=nj=1Bhjλhj,Xhnj=1Xhjλhj,nj=1λhj=1,(j=1,...,n)} (4)

    Subscript, (j) shows jth (DMU) and λj indicates the jth intensity variable (j = 1, …, n). The inequality constraints (Xhnj=1Xhjλhj), (Gjnj=1Ghjλhj) indicates strong disposability on inputs and desirable outputs from hth division, respectively and Bh=nj=1Bhjλhj measures congestion on undesirable outputs from hth division. Similarity, strong disposability is specified on the two output vectors as follows.

    Phs(x)={(Gh,Bh):Ghnj=1Ghjλhj,Bhnj=1Bhjλhj,Xhnj=1Xhjλhj,nj=1λhj=1,(j=1,...,n)} (5)

    The inequality constraint Bnj=1Bjλj allow for strong disposability on undesirable outputs. The constraint nj=1λj=1 is incorporated into the two expressions which indicate variable return to scale in production processes. The production technology set to definition of natural and managerial Disposability is specified by the following two types of output vectors and an input vector for hth division of the supply chain as follows.

    PhN(x)={(Gh,Bh):Ghnj=1Ghjλhj,Bhnj=1Bhjλhj,Xhnj=1Xhjλhj,nj=1λhj=1,(j=1,...,n)} (6)
    PhM(x)={(Gh,Bh):Ghnj=1Ghjλhj,Bhnj=1Bhjλhj,Xhnj=1Xhjλhj,nj=1λhj=1,(j=1,...,n)} (7)

    Here PhN(x) is defined as a production possibility set under natural (N) disposability and PhM(x) managerial (M) disposability one from hth division. The production technology under natural disposability or PhN(x) has Xhnj=1Xhjλhj mentioned an organization can reduce a directional vector of input to attain to efficiency frontier. Likewise, The production technology under managerial disposability or PhM(x) has Xhnj=1Xhjλhj mentioned a directional vector of inputs from hth division of the supply chain can increase to attain to the efficient frontier. The both of production technology set have common constraints, Ghnj=1Ghjλhj, Bhnj=1Bhjλhj under natural and managerial disposability.

    Let us consider the general structure of the supply chain depicts in Figure 1. Let us consider, xhmj,ghrj,bhfj, whej indicate mth input (m=1,...,M), rth desirable outputs (r=1,...,S) and fth undesirable outputs (f=1,...,F) and eth dual-role factors (e=1,...,E) of h th division (h=1,...,H) in jth (j=1,...,n) supply chain, respectively. Also, ˉxhmj, ˜xhmj indicate original m inputs are separated into two categories m and m+, as M=m+m+. Furthermore, v(h,h)pj represent the intermediate measures between the h th division to the h th division of jth supply chain. The subscript (p, j) indicating pth intermediate measure (p=1,...,Ph) in jth supply chain (j=1,...,n) and z(h,h)aj represent invers intermediate measures exit from h th division and enter to h th division. The subscript (a, j) indicating ath intermediate measure (a=1,...,Ah) in jth supply chain. (j=1,...,n).

    The production technology set of h th division in the jth supply chain is defined as follows:

    Y={(vhj,zhj,ghj,bhj,whjxhj)|xhjcanproduce(vhj,zhj,yhj,whj)}. Thus, the outputs set of h th division in the jth supply chain can be indicated as follows:

    Phj(x)={(vhj,zhj,ghj,bhj,whj)|(vhk,zhj,ghj,bhj,whj,xhj)Y}
    Figure 1.  The general structure of supply chain.

    Let us now suppose a supply chain (DMU) is concluded from five-stage, supplier, Manufacture, transmitter, distributor, and customer. We treat each supply chain as a DMU. Let us consider hs,hm,ht,hd,hc the number of divisions in the supplier, manufacturer, transmitter, distributor and customer. Figure 2 shows an electricity supply chain structure in the power industry. The electricity supply chains are power suppliers in power production activities. They are comprised of fuel suppliers (oil and gas fields), power producers (power plants), electricity transmitters (transmission lines), power distributors (distribution lines) and final customers. These entities collaborate to power production and management in economic business.

    Figure 2.  The supply chain structure.

    In this study, the supply chains have been built in northern, southern, eastern, western and central districts in Iran. In this conformation Oil and gas fields and refineries provide demand fuels of power plants and district power plants Transfer produced power by regional power companies to the area distribution companies to dispatching to consumers or residents of their area. Other words, each supply chain or DMU is built of five stages and partners of each stage connected by intermediate measures to the successor stage. Supply chains are comparable and compete in the power industry. In Figure 2 is depicted intermediated measures sent from oil and gas fields to power plants, from power plants to transmissions companies, from transmissions companies to distributions companies and finally from them to customers. Furthermore, the inverse intermediate measures exit from transmitter divisions and enter to manufacture divisions and exit from manufacture divisions and enter to supplier divisions. These measures indicate entities' relationship in the supply chain. However, each division of entities operates independent from other divisions of per stage in production activities and supply chains compete to high efficiency earn in economic business (see Pouralizadeh et al. [24]).

    In this section, we propose a DEA model to sustainability assessment a supply chain. We suppose a supply chain contains an arbitrary number of suppliers, manufacturers, transmitters, distributors and customers. The model (1) be further developed as a network model by incorporate the two categories intermediate measures and dual-role factors for each supply chain division in order to efficiency assessment of the overall supply chain.

    We shall assume the inputs separate into two categories under natural and managerial disposability, weak disposability of good outputs reduction, free disposability of undesirable outputs and convexity and variable returns to scale in the production process to calculate inefficiency score. In this study we considered the different weights for partners of a particular stage of the network supply chain as Wh,(h=1,...,H) are weights for H divisions that are defined by decision-makers in production activities. In this method, the inefficiency performance evaluating of an overall supply chain can be formed by the inefficiency performance evaluating of all its divisions similar to model (1). The production factors of the jth supply chain (DMU) are summarized as follows:

    ˉXhj=(ˉxh1j,ˉxh2j,...,ˉxhij)T>0 : The input ith under natural disposability from hth division in jth supply chain, i=1,...,m,h=1,...,H, j=1,...,n.

    ˜Xhj=(˜xh1j,˜xh2j,...,˜xhqj)T>0 : The input qth under managerial disposability from hth division in jth supply chain, q=1,...,m+,h=1,...,H, j=1,...,n.

    Ghrj=(gh1j,gh2j,...,ghrj)T>0 : The desirable output rth from hth division in jth supply chain r=1,...,s,h=1,...,H,j=1,...,n.

    Bhj=(bhij,bh2j,...,bhfj)T>0 : The undesirable output fth from hth division in jth supply chain f=1,...,F,h=1,...,H,j=1,...,n.

    Whj=(wh1j,wh2j,...,whej)T>0 : The dual-role factor eth from hth division in jth supply chain e=1,...,E,h=1,...,H,j=1,...,n.

    V(h,h)j=(v(h,h)1j,v(h,h)2j,...,v(h,h)pj)T>0 : The pth Material flow or intermediate measure from division h to division h in jth supply chain, p=1,...,P,h=1,...,H,j=1,...,n.

    Z(h,h)j=(z(h,h)1j,z(h,h)2j,...,z(h,h)aj)T>0 : The ath invers intermediate measure from division h to division h in jth supply chain, a=1,...,A,h=1,...,H,j=1,...,n.

    s(h,h)pj : The slack variables of the pth intermediate measure from divisions h to divisions h in jth supply chain, (p = 1, …, P), (j = 1, …, n).

    s(h,h)aj0 : The input slack variables of the ath invers intermediate measure from division h to division h in jth supply chain (a = 1, …, A), (j = 1, …, n).

    s+(h,h)aj0 : The output slack variables of the ath intermediate measure or invers flow from division h to division h in jth supply chain (a = 1, …, A), (j = 1, …, n).

    λh=(λh1,λh2,...,λhn)T : An unknown column vector.

    Rhi=(Mh+Sh+Fh+Eh+Dh)1(max{xhij|j=1,...,n}min{xhij|j=1,...,n})1 : A data range related to i th input in hth division. i=1,...,m,h=1,...,H.

    Rhq=(Mh+Sh+Fh+Eh+Dh)1(max{xhqj|j=1,...,n}min{xhqj|j=1,...,n})1 : A data range related to qth input in hth division. h=1,...,H, q=1,...,m+.

    Rhf=(Mh+Sh+Fh+Eh+Dh)1(max{bhfj|j=1,...,n}min{bhfj|j=1,...,n})1 : A data range related to f th undesirable output input in hth division. h=1,...,H, f=1,...,F.

    Rp=(Mh+Sh+Fh+Eh+Dh)1(max{v(h,h)pj|j=1,...,n}min{v(h,h)pj|j=1,...,n})1 : A data range related to p th intermediate measure sent from h th division to h th divisions.

    p=1,...,P,hh,h,h{1,...,H}

    Ra=(Mh+Sh+Fh+Eh+Dh)1(max{z(h,h)aj|j=1,...,n}min{z(h,h)aj|j=1,...,n})1 : A data range related to a th invers intermediate measure sent from sent from h th division to h th divisions. a=1,...,A,hh,h,h{1,...,H}.

    ξh : Inefficiency score of hth division.

    ε : A small amount and it considered as 0.0001 for computation convenience.

    In proposed approach, the number of original m inputs of hth division are separated into two categories, mh (under natural disposability) and m+h (under managerial disposability), respectively. The model maintains Mh=mh+m+h. Also, s(h,h)p is slack variable of pth the intermediate measure (p = 1, ..., P) sent from hth division to h'th division and, s+(h,h)a is defined as slack variables of ath the inverse intermediate measures (a=1,...,A) sent from h th division to h th division. Also, the inverse intermediate measures enter to divisions are considered as non-discretionary inputs set and, the inverse intermediate measures exits from divisions are specified as desirable outputs set in model. The column vectors of structural variables (λh) are applied for connecting the input, desirable and undesirable output vectors, the dual-role factors and the set intermediate measures by convex combination under variable return scale in h th division. Mh,Sh,Fh,Eh,(h=1,...,H), indicate the total number of inputs, the desirable and undesirable outputs, the dual-role factors in h th division. Also, Ph,Ah show the total number of intermediate measures sent from h th division to the h th division and the inverse intermediate measures exit from h th division and enter to h th division (h, h : h,h=1,...,H), respectively.

    In proposed model Rhi,Rhq,Rhf,Rhp,Rha are specified by the decision maker for h th division as follows:

    Rhi=(Mh+Sh+Fh+Eh+Ph+Ah)1(max{ˉxhij|j=1,...,n}min{ˉxhij|j=1,...,n})1Rhq=(Mh+Sh+Fh+Eh+Ph+Ah)1(max{˜xhqj|j=1,...,n}min{˜xhqj|j=1,...,n})1Rhf=(Mh+Sh+Fh+Eh+Ph+Ah)1(max{bhfj|j=1,...,n}min{bhfj|j=1,...,n})1Rhp=(Mh+Sh+Fh+Eh+Ph+Ah)1(max{v(h,h)pj|j=1,...,n}min{v(h,h)pj|j=1,...,n})Rha=(Mh+Sh+Fh+Eh+Ph+Ah)1(max{z(h,h)aj|j=1,...,n}min{z(h,h)aj|j=1,...,n})1 (8)

    Moreover, slack variables correspond to inverse intermediate flows that are considered as non-discretionary inputs sets are not include in objective function and their corresponding constraints set is followed by the '*' symbol. Unified efficiency score is obtained by subtracting the level of inefficiency from unity. A unified efficiency score under natural and managerial disposability is measured from the supply chain as follows:

    UENM=1[ξ+ε(Rhidhi+Rhqdhq+Rhfdhf+Pp=1Rps(h,h)p+Aa=1RaS+(h,h)a)] (9)

    The objective function of DMU (supply chain) calculates by weighted average of optimal inefficiency of each division of the supply chain so the objective function weights could be obtained through an expert opinion process. Therefore, the inefficiency scores and all slack variables are determined on the optimality model as follows:

    θ=MaxHh=1Wh[ξh+ε(mi=1Rhidhi+m+q=1Rhqdhq+Ff=1Rhfdhf+Hh=1Pp=1Rps(h,h)p+Hh=1Aa=1Ras+(h,h)a]nj=1ˉxhijλhj+dhi=ˉxhiki=1,...,mh,h=1,...,Hnj=1˜xhqjλhjdhq=˜xhqkq=1,...,m+h,h=1,...,Hnj=1ghrjλhj+ξhghrk=ghrkr=1,...,Sh,h=1,...,Hnj=1bfjλhjdhf=bhfkf=1,...,Fh,h=1,...,Hnj=1whejλhj=wheke=1,...,Es,h=1,...,hsnj=1whejλhj=wheke=1,...,Em,h=1,...,hmnj=1whejλhj=wheke=1,...,Et,h=1,...,htnj=1λhjv(h,h)pj+s(h,h)p=nj=1λhjv(h,h)pjh=1,...,hs,p=1,...,Ps,h=1,...,hmnj=1λhjv(h,h)pj+s(h,h)p=Jj=1λhjv(h,h)pjh=1,...,hm,p=1,...,Pm,h=1,...,htnj=1λhjv(h,h)pj+s(h,h)p=nj=1λhjv(h,h)pjh=1,...,ht,p=1,...,Pt,h=1,...,hdnj=1λhjv(h,h)pj+s(h,h)p=nj=1λhjv(h,h)pjh=1,...,hd,p=1,...,Pd,h=1,...,hcnj=1λhjz(h,h)ajs+(h,h)a=z(h,h)akh=1,...,hs,h=1,...,hm,a=1,...,Asnj=1λhjz(h,h)aj+s+(h,h)a=z(h,h)akh=1,...,hs,h=1,...,hm,a=1,...,Amnj=1λhjz(h,h)ajs(h,h)a=z(h,h)akh=1,...,hm,h=1,...,ht,a=1,...,Amnj=1λhjz(h,h)aj+s(h,h)a=z(h,h)akh=1,...,hm,h=1,...,ht,a=1,...,Atnj=1λhj=1h=1,...,H,j=1,...,nλj0,s+(h,h)a0,s(h,h)a0,s(h,h)p,ξUR,j=1,...,n,h=1,...,H (10)

    Therefore, efficiency score on DMU is measured by β=1θ where the inefficiency score and all slack variables correspond to inputs under natural and managerial disposability and undesirable outputs and the two set intermediate measures are determined on the optimality of model (10). In a result, the inefficiency of the overall supply chain can be formed of weighted average of all of its partner's inefficiency in production processes as model (10). The first and second constraints categories correspond to inputs set under natural and managerial disposability, respectively. Also, the third and the fourth constraints categories related to desirable and undesirable outputs, respectively and the fifth, sixth, seventh, constraints categories are correspond to dual-role factors of supplier, manufacture and transmitter divisions. The eighth, ninth, tenth and eleventh the categories constraints correspond to intermediate measures sent from supplier divisions to manufacturer divisions, and manufacturer divisions to transmitter divisions, and from transmitter divisions to distributor divisions and from them to customer divisions, respectively.

    The twelfth and thirteenth the categories constraints related to inverse intermediated measures exit from manufacturer divisions and enter to supplier divisions. Also, the fourteenth and fifteen the categories constraints correspond to inverse intermediate measures exit from transmitter divisions and enter to manufacture divisions. The last constraints categories related to variable returns to scale in the production process. This model measures an investment opportunity for technology innovation for reducing the number of industrial pollutions (flaring gas) in oil and gas fields and power plants sectors and preventing from power losses in transmission and distribution lines. Moreover, these approaches examine the level of unified efficiency by a single inefficiency score that is assigned to desirable outputs. Meanwhile, constraints on the desirable output nj=1ghrjλhj+ξhghrk=ghrk do not have any slack so that they can be considered as equality, so belonging to weak disposability and other constraints relate to inputs and undesirable outputs maintain slacks in the model (10). Thus, these constraints on all inputs and undesirable output are considered as inequality, so implying the concept of strong disposability.

    Let us suppose, thi(i=1,...,m), lhq(i=1,...,m+), uhr(r=1,...,s), chf(f=1,...,F), yhe(e=1,...,E), present the dual variables correspond to the categories constraints of the inputs under natural and managerial disposability, desirable and undesirable outputs and dual-role factors from hth division, (h=1,...,H) in the model(10), respectively. Moreover, let us consider, ˉBp, Bp, ˜Bp ˆBp the dual variables correspond to eighth, ninth, tenth and eleventh of the categories constraints related to intermediate measures which are sent from the supplier divisions to manufacture divisions and from manufacture divisions to transmitter divisions and from them to distributor divisions finally from distributer divisions to customer divisions, respectively. Likewise, we suppose ˉIa, Ia, ˜Ia, ˆIa, present the dual variables related to the twelfth, thirteenth, fourteenth and fifteen of the categories constraints of inverse intermediate measures which exit from manufacture divisions, entire to supplier divisions and exit from transmitter divisions and entire to manufacture divisions, respectively. Furthermore, the dual variable σh is obtained from the last equation from the model (10) in hth division.

    The dual formulation of model (10) is as follows:

    minZ=Hh=1(mi=1thiˉxhikm+q=1lhq˜xhqk+Sr=1urghrkFf=1cfbhfk+Ee=1yhewhek+h=hmAa=1ˉIaz(h,h)ak+h=htAa=1ˆIaz(h,h)ak+σh)s.tmi=1thiˉxhijm+q=1lhq˜xhqj+Sr=1urghrjFf=1cfbhfj+Ee=1yhewhej+Pp=1ˉBpv(h,h)pj+Aa=1Iaz(h,h)aj+σh0h=1,...,hs,h=1,...,hm,j=1,...,nmi=1thiˉxhijm+q=1lhq˜xhqj+Sr=1urghrjFf=1cfbhfj+Ee=1yhewhej+Pp=1Bpv(h,h)pjPp=1ˉBpv(h,h)pj+Aa=1ˉIazaj(h,h (11)

    According to model (10) the supporting hyper plane is expressed for an arbitrary division from power customers as follows:

    t_{}^h\, \bar x_{}^h\, - l_{}^h\tilde x_{}^h + u_{}^h\, g_{}^h\, - c_{}^{{h_{}}}b_{}^{{h_{}}} + \, w_{}^{{h_{}}}y_{}^{{h_{}}} + \, {\sigma ^h}\, = 0\, \, \, \, \, \, \, \, \, \, \, \, \, \, \, \, \, h = 1, ..., {h_c} (12)

    In this case, all production factors have a single component. The concept of DTR defined as ({{db} \mathord{\left/ {\vphantom {{db} {dg{)_h}}}} \right. } {dg{)_h}}}/{({b \mathord{\left/ {\vphantom {b g}} \right. } g})_h} for h th division in the case a single component of the two production factors. Based upon the sign of ({{db} \mathord{\left/ {\vphantom {{db} {dg{)_h}}}} \right. } {dg{)_h}}}/{({b \mathord{\left/ {\vphantom {b g}} \right. } g})_h} the type of a supporting hyper plane is specified for an arbitrary division from power customers on desirable output (g) and undesirable output (b) as follows:

    (a) If ({{db} \mathord{\left/ {\vphantom {{db} {dg{)_h}}}} \right. } {dg{)_h}}}/{({b \mathord{\left/ {\vphantom {b g}} \right. } g})_h} > 0 then the DTR is as positive

    (b) If ({{db} \mathord{\left/ {\vphantom {{db} {dg{)_h}}}} \right. } {dg{)_h}}}/{({b \mathord{\left/ {\vphantom {b g}} \right. } g})_h} = 0 then the DTR is as zero.

    (c) If ({{db} \mathord{\left/ {\vphantom {{db} {dg{)_h}}}} \right. } {dg{)_h}}}/{({b \mathord{\left/ {\vphantom {b g}} \right. } g})_h} < 0 then the DTR is as negative.

    After solving the Model (11) the desirable outputs congestion or technology innovation for h th division is identified under assumption of a unique optimal solution by the dual variables correspond to desirable output constraints as follows:

    (a) If {(u_r^h)^ * } = 0 for some (at least one) r then the 'zero DTR' occurs on h th division from the supply chain under consideration.

    (b) If {(u_r^h)^*} < 0 for some (at least one) r then the 'negative DTR' occurs on h th division from the supply chain under consideration.

    (c) If {(u_r^h)^*} > 0 for all r then the 'positive DTR' occurs on h th division from the supply chain under consideration.

    Note, If {(u_r^h)^*} < 0 some r and {(u_r^h)^ * } = 0 for other r, then we consider that the negative DTR occurs on h th division from the supply chain under consideration. In other words, this case indicates a status of desirable congestion or technology innovation on undesirable outputs. Furthermore, if {(u_r^h)^*} < 0 for all r then, this case indicates the best status because technology innovation increase all of the desirable outputs and increase in any desirable output always abatement undesirable outputs. Furthermore, If {(u_r^h)^*} < 0 is defined for some r, then it indicates a case to abatement an amount of undesirable outputs. Therefore, the effect of investment is specified by the dual variable {(z_q^h)^ * } as if {(z_q^h)^ * } > W\varepsilon R_q^h then the qth the input for investment under managerial disposability able to decrease the amount of undesirable output in h th division and if {(z_q^h)^ * } = W\varepsilon R_q^h then the qth the input has a limited effect on reducing of undesirable output.

    In this section we apply the proposed model to the analysis of the power industry in Iran. In Subsection 4.1 we will describe the dataset and we will specify the inputs and outputs we will consider in our analysis, in Subsection 4.2 we will present the main results.

    The stylized supply chain in the power industry can be summarized in five main actors: gas and fuel suppliers, power generators, transmission networks, distribution facilities, and final users. Conventional power plants consume fuel oil, natural gas and diesel to produce electricity, while renewable ones are solar, wind and hydro plants. Conventional plants can be further divided depending on the kind of technology adopted, in thermal, gas and combined cycle plants. In general, thermal power plants operated by fossil fuels produce huge amounts of air pollutants. The pollutants which have been considered in the study are sulfur oxides (SOX), nitrogen oxides (NOX) and carbon dioxide (CO2).

    Our purpose is to highlight the theoretical and practical quality of the model, therefore each of the DMUs or the supply chain is built of five stages and each stage includes a set of partners connected to the predecessor stages members by some sustainable intermediate measures. In our application, we consider 10 supply chains (DMUs) including oil and gas fields (suppliers) that provide different fuels to power stations, power plants (manufacturers), regional power companies (transmitters), distribution companies (distributors) and customers. Per each supply chain, we consider two suppliers: oil and gas companies that satisfy the fuel demand of power plants (intermediate product) and that can also sell fuels as final output. Suppliers use one input (capital) under natural disposability and one input under managerial disposability (labor) and produce one desirable (oil or gas) and one undesirable output (flaring gas). The dual-role factor is considered as the cost of cleanup flare gas pollutions. Each manufacturer includes at least three power plants with different technologies (thermal, combined cycle, gas, hydro, wind and solar). They use fuels, capital and labor (under natural disposability) and labor of hydro power plant under managerial disposability to produce electricity and they sell it to regional power companies. To update and enlarge their capacity, manufacturers can substitute existing plants with more efficient ones or they can construct new plants. Three undesirable outputs are considered for manufacturers: CO2, Nox and SOX emissions. Also, we consider the dual-role factor as the Inner consumptions of power plants as technical and nontechnical consumptions. The transmitters transfer electricity from manufacturers to distributing companies and capacity and length of the lines are considered as inputs under natural disposability and the number employees of the department of programing and researches are used as input under managerial disposability. The dual-role factor is considered as specialist workforce in programming and researches. The loose in the transmission lines is considered as undesirable output while the construction of new lines is a desirable one. Distribution companies receive electricity from transmitters and dispatch them to the final consumers. They use two additional inputs capital estimated as capacity of the distribution lines and length of the distribution lines under natural disposability and the number of employees of engineering assistance department and programming as input under managerial disposability, one final desirable output as the meter of electricity and one undesirable output that is losses in the distribution lines. Finally, customers are classified as residential, agriculture, public and industrial. They use one input under natural disposability and one input under managerial disposability and produce two desirable outputs and one undesirable output. Table 1 indicates the production factors used for supply chain evaluation.

    Table 1.  Production factors in performance evaluation.
    Division Numerator Factors Definition
    Supplier {h_s} \bar x_{1j}^{h(s)} Capacity of oil (103 Barrels) and gas(106 m3)
    \tilde x_{1j}^{h(s)} Number of employees
    g_{1j}^{h(s)} Oil (103 Barrels) and gas (106 m3) sold
    b_{1j}^{h(s)} Flaring gas of oil field (103 barrels)and gas field(106 m3)
    w_{1j}^{h(s)} Cost of flaring gas recovery
    Manufacture {h_m} \bar x_{1j}^{h(m)} Power nominal of power plants
    \bar x_{2j}^{h(m)} Labor
    \tilde x_{1j}^{h(m)} Labor of hydro plant
    g_{1k}^{{h_m}} Percentage of new construction of power plant
    b_{1k}^{{h_m}} Emissions of Nox harmful substances(103 Kg/106 Kwh)
    b_{2k}^{{h_m}}
    b_{3k}^{{h_m}}
    w_{1j}^{h(m)}
    Emissions of Sox harmful substances(103 Kg/106 Kwh).
    Emission of Co2 harmful substances(103 Kg/106 Kwh)
    Inner consumption of power plant
    TransmitterDistribution {h_t}
    {h_d}
    \bar x_{1j}^{h(t)}
    \bar x_{2j}^{h(t)}
    \tilde x_{1j}^{h(t)}
    g_{1k}^{{h_t}}
    w_{1j}^{h(t)}
    \bar x_{1j}^{h(d)}
    \bar x_{2j}^{h(d)}
    \tilde x_{1j}^{h(d)}
    g_{1k}^{{h_d}}
    Capacity of regional
    company (Mwa)
    Length transmission line (Km circuit).
    Labor
    New construction of transmission lines (Km)
    Number of employees
    Capacity of distribution
    (Mwa)
    Length transmission line (Km).
    Labor
    New construction of distribution lines (Km).
    b_{1k}^{{h_d}} Percentage of losses of distribution line (%).
    Customer \bar x_{1j}^{h(c)} Average cost with fuel subsidy (Rial).
    \tilde x_{1j}^{h(c)}
    g_{1k}^{{h_c}}
    g_{2k}^{{h_c}}
    b_{1k}^{{h_t}}
    v_{mk}^{(h, h')}
    z_{a\, j}^{({h_m}, {h_s})}
    z_{a\, j}^{({h_t}, {h_m})}
    Direct selling of electricity (106 Kwa).
    Number of customer
    Sales of electricity
    (106 Kwh)
    Cut of power
    Material flow from division h to division h' (106 Kwa)
    Invers intermediate measures sent from manufactures divisions to supplier
    Invers intermediate measures sent from transmitters to manufacture

     | Show Table
    DownLoad: CSV

    More in detail, the parameters used to characterize this supply chain are defined as follows:

    {h_s} : Numerator of divisions in the supplier level ({h_s} : 1, 2).

    \bar x_{1j}^{h(s)} : Capacity of oil (103 Barrels) and gas (106 m3) fields of {h_s} th supplier in jth supply chain.

    \tilde x_{1j}^{h(s)} : Number of employees from {h_s} th supplier in jth supply chain.

    g_{1j}^{h(s)} : Oil (103 Barrels) and gas (106 m3) sold to other companies from the {h_s} th supplier in jth supply chain.

    b_{1j}^{h(s)} : Flaring gas of oil field (103 barrels) and gas field (106m3) of the {h_s} th supplier in the jth supply chain.

    w_{1j}^{h(s)} : The cost of cleanup of burned gas (flaring gas) of {h_s} th supplier in jth supply chain.

    {h_m} : Numerator of division in the manufacturer level ({h_m} : 3, 4, 5).

    \bar x_{1j}^{h(m)} : Power nominal of {h_m} th manufacturer in the jth supply chain (106 Kwh).

    \bar x_{2j}^{h(m)} : Number of employees of {h_m} th manufacturer in the jth supply chain.

    \tilde x_{1j}^{h(m)} : Number of hydropower employees of {h_m} th manufacturer in the jth supply chain.

    g_{1j}^{h(m)} : Percentage of new construction of power plant of the {h_m} th manufacturer in the jth supply chain.

    b_{1j}^{h(m)} : Emissions of Nox harmful substances of the {h_m} th manufacturer in the jth supply chain (103Kg/106Kwh).

    b_{2j}^{h(m)} : Emissions of SoX harmful substance of the {h_m} th manufacturer in the jth supply chain (103Kg/106Kwh).

    b_{3j}^{h(m)} : Emission of CO2 harmful substance of the {h_m} th manufacturer in the jth supply chain (103 Kg/106 Kwh).

    w_{1j}^{h(m)} : Inner consumption of power plants (technical and nontechnical consumptions) of the {h_m} th manufacturer in the jth supply chain (106 Kwh).

    {h_t} : Numerator of the divisions the level of the transmitters ({h_t} : 6, 7).

    \bar x_{1j}^{h(t)} : Capacity of transmission lines of the {h_t} th transmitter in the jth supply chain (Mwa).

    \bar x_{2j}^{h(t)} : Length transmission line of the {h_t} th transmitter in the jth supply chain (Km circuit).

    \tilde x_{1j}^{h(t)} : Number of employees department of programing and researches of the {h_t} th transmitter in the jth supply chain.

    g_{1j}^{h(t)} : New construction of transmission lines of the {h_t} th transmitter in the jth supply chain (Km circuit).

    b_{1j}^{h(t)} : Loose of transmission line of {h_t} th transmitter in the jth supply chain (%).

    w_{1j}^{h(t)} : Number of employees of deputy transfer and exploitation of {h_t} th transmitter in the jth supply chain.

    {h_d} : Numerator of division in the distributer level ({h_d} : 8, 9, 10, 11).

    \bar x_{1j}^{h(d)} : Capacity of distribution lines of {h_d} th distributer in the jth supply chain (Mwa).

    \bar x_{2j}^{h(d)} : Length distribution line of the {h_d} th distributer in the jth supply chain (Km).

    \tilde x_{1j}^{h(d)} : Number of employees of engineering assistance department and programming of the {h_d} th distributer in the jth supply chain.

    g_{1j}^{h(d)} : Meter of electricity of {h_d} th distributer in jth supply chain.

    b_{1j}^{h(d)} : Percentage of losses of distribution line of {h_d} th distributer in the jth supply chain.

    {h_c} : Numerator of division in the customer level ({h_c} : 12, 13, 14, 15).

    x_{1j}^{h(c)} : Average cost with fuel subsidy of the {h_c} th customer in the jth supply chain (Rial).

    \tilde x_{1j}^{h(c)} : Direct selling of electricity from transmitter Company to the {h_c} th customer in the jth supply chain (106 Kwh).

    g_{1j}^{h(c)} : Number of customers of {h_c} th customer in the jth supply chain.

    g_{2j}^{h(c)} : Sales of electricity of the {h_c} th customer in the jth supply chain (106 Kwh).

    b_{2j}^{h(c)} : Cut off power of the {h_c} th customer in the jth supply chain (minute/year).

    v_{p\, j}^{(h, h')} : Material flow from division h to division h' (106 Kwa).

    z_{a\, j}^{({h_m}, {h_s})} : Power flow sent from power plants to oil and gas fields (106 Kwa).

    z_{a\, j}^{({h_t}, {h_m})} : Labor sent from regional companies to power plants to repair and maintenance of systems.

    The dataset has been collected from the power industry company in Iran and the reference year is 2015 (see TAVANIR website for the detailed data). The total emissions due to electricity generation in Iran, the amount and type fuel used in all power plants have been considered in the computation of undesirable outputs. All the data of the two oil and gas fields (suppliers), power plants (manufacturers), regional power companies (transmitters), distribution companies (distributors) and customers (residential, public, agriculture, industrial) are available in the TAVANIR website [30]. Supplier inputs are obtained from oil and gas fields statistics of the energy industry in Iran. The desirable output is computed as the difference between the average annual production and the amount of oil and gas that are sent to power plants; undesirable output (flaring gas) is calculated with a 0.03% rate of the annual production of oil and gas. Information related to the demanding fuel of power plants is collected from TAVANIR Company [30] in the power industry and they are considered as intermediate measures from oil and gas fields to power plants. The capacity of power plants is a proxy of the input capital. Undesirable outputs for manufacturers are computed based on the amount of electricity produced by the different power plants using different technologies and fuels. Dataset of inputs and desirable output of regional power company are collected from the transmission Division of TAVANIR Company in power industry and losses of the transmission line (undesirable output) are estimated with a 3.02% factor based on the amount of loose of transmission in Iran. All of the data of distribution company are obtained from dispatch division of TAVANIR company in power industry likewise input of customer divisions are collected from TAVANIR company and desirable output of customers are computed as total sale of electricity to residential, public, agriculture and industry divisions but undesirable output is computed by time cut off of electricity in different divisions of consumers in 2015 (see Pouralizadeh et al. [24]).

    The data sets corresponding to the 10 supply chains (DMUs) under analysis are presented in Tables 217. Tables 2 and 3 shows inputs under natural and managerial disposability and desirable and undesirable outputs for suppliers 1 and 2. In Tables 47, we present the data of manufacturer (level 1, 2, 3). Tables 8 and 9 show the data of transmitters with two inputs under natural disposability and one input under managerial disposability, one good output and one undesirable output. Tables 1013 collect the data on distributors where two inputs under natural disposability and, one input under managerial disposability, one desirable and one undesirable output are considered. Finally, in Tables 1417 the data of customers are reported with one input under natural disposability and one input under managerial disposability, two desirable outputs and one undesirable output.

    Table 2.  The supplier level-inputs.
    DMU supplier 1 (division 1) supplier 2 (division 2)
    \bar x_{1k}^1 \tilde x_{1k}^1 \bar x_{1k}^2 \tilde x_{1k}^2
    1 2550 3200 7200 2500
    2 61200 1300 21600 2500
    3 21600 3200 10800 2400
    4 32400 3110 6480 1400
    5 12600 2800 19440 3000
    6 43200 2200 10800 2400
    7 46800 2400 10800 1380
    8 39600 1600 21600 2250
    9 9360 2150 19440 2180
    10 64800 2500 6480 2900
    Source: category: oil field of iran-wikipedia, https//en.wikipedia.org/wiki/category:oil fields of iran; https//en.wikipedia.org/wiki/category:ntural gas in iran

     | Show Table
    DownLoad: CSV
    Table 3.  The supplier level desirable and undesirable outputs and dual-role factors.
    DMU
    w_{1k}^1
    Supplier1
    (Division1)
    g_{1k}^1
    b_{1k}^1 w_{1k}^2 Supplier2
    (Division 2)
    g_{1k}^2
    b_{1k}^2
    1 0.011 1739.693 54 4.725 1186.216 151.2
    2 0.255 40572.996 1296 10.8 7203.230 345.6
    3 0.085 8995.883 432 5.738 3726.203 183.6
    4 0.191 26527.191 972 4.388 1930.025 140.4
    5 0.042 4552.857 216 11.475 10438.190 367.2
    6 0.149 23324.391 756 5.738 3350.675 183.6
    7 0.149 17080.471 756 5.4 2353.130 172.8
    8 0.127 15872.914 648 10.8 9455.104 345.6
    9 0.038 6062.772 194.4 11.475 9849.593 367.2
    10 0.255 25603.400 1296 4.388 2208.415 140.4
    Calculation Flaring gas and Sold oil and gas

     | Show Table
    DownLoad: CSV
    Table 4.  Manufacturers level inputs.
    DMU Manufacturer1 Manufacturer 2 Manufacturer3
    \bar x_{1k}^3 \bar x_{2k}^3 \tilde x_{1k}^3 \bar x_{1k}^4 \bar x_{2k}^4 \tilde x_{1k}^4 \bar x_{1k}^5 \bar x_{2k}^5 \tilde x_{1k}^5
    163224 4070610 154081600 511903 12000
    216200 226327 10400700 02626.952 260027
    310448 10000 5701.123300 016760 200526
    480224 10007 8622.43300 08344 20050
    55184 8900 1920.48900 716417.760 28230
    613672.88 23000 33122500 353936 8000
    7966.32 14500 83522700 3417844.8 8900
    81491.2 152021 103202260 916800 13000
    93872 15000 105903600 177072 4100106
    1011453.6 318040 6787.2760 02053.28 15900
    Source: http//amar.tavanir.org.ir//tolid and calculations million kilo watt hour

     | Show Table
    DownLoad: CSV
    Table 5.  Manufacturers level desirable and undesirable outputs and Dual-role factor.
    DMU Manufacturer 1 (Division 3)
    w_{1k}^3 g_{1k}^3 b_{1k}^3 b_{2k}^3 b_{3k}^3
    1598.234 12.2 454610.278 23891876.280 288025420.100
    292.234 12.2 302399.805 4207069.806 191952930.500
    3180.638 13 235104.740 195553.061 149621794
    4394.18 12.2 229464.218 12059407.75 145380628.200
    510.78 73.6 43498.708 38755.471 27536231.770
    625.768 100 256638.343 217529.667 163094448.800
    72.939 85.5 6683.633 5954.829 4230977.926
    881.863 85.5 15138.687 184259.151 9585079.623
    942.59 13 92035.892 76552.691 58572086.910
    10139.981 86.6 236364.062 196600.528 150423232.700
    Source: http//amar.tavanir.org.ir//tolid and calculations 1000kg/million kilo watt hour

     | Show Table
    DownLoad: CSV
    Table 6.  Manufacturers level desirable and undesirable outputs and Dual—role factor.
    DMU Manufacturer 2 (Division 4)
    w_{1k}^4 g_{1k}^4 b_{1k}^4 b_{2k}^4 b_{3k}^4
    10 85.5 5715.366 5092.145 3618030.390
    2541.271 0 283431.105 14895617.700 179572190
    3291.571 12.2 174773.192 9070013.802 110729096.200
    486.474 25.2 182851.984 152090.788 116367887.400
    596.326 12.2 49845.037 2619587.603 3158009.070
    610.299 85.5 27420.014 24430.049 17357845.530
    7424.975 12.2 273496.466 14373506.370 173277944.500
    80.063 12.2 311634.456 21776302.480 197440862.200
    9102.151 98.8 176752.534 147351.908 112467128.500
    10170.387 86.6 79593.197 66419.786 50641168.170
    Source: http//amar.tavanir.org.ir//tolid and calculations 1000kg/million kilo watt hour

     | Show Table
    DownLoad: CSV
    Table 7.  Manufacturers level desirable and undesirable outputs and Dual—role factor.
    DMU Manufacturer 3
    (Division 5)
    w_{1k}^5 g_{1k}^5 b_{1k}^5 b_{2k}^5 b_{3k}^5
    10 73.600 19603.894 17519.680 12447945.190
    26.325 73.600 27423877.76 24433491.25 17360291475
    3103.532 98.800 212448.268 690393.877 135090771.800
    492.426 13 140748.540 117070.408 89573051.780
    547.29 87 300157.654 9178172.226 190308335.200
    635.747 13 77463.980 64432.212 49298451.340
    7290.054 13 471751.939 21768344.370 299051808
    8782.679 13 510495.755 21776302.480 323709891.900
    945.519 13 94829.614 78876.425 60350025.180
    10138.404 1.200 59895.401 3147780.793 37947663.670
    Source: http//amar.tavanir.org.ir//tolid and calculations 1000kg/million kilo watt hour

     | Show Table
    DownLoad: CSV
    Table 8.  The Transmitter level inputs.
    DMU Transmitter 1 (division 6) Transmitter 2 (division 7)
    \bar x_{1k}^6\tilde x_{1k}^6 \bar x_{2k}^6 \bar x_{1k}^7\tilde x_{1k}^7 \bar x_{2k}^7
    1 2754274 8704 2508639 14697.700
    2 4101178 9127.800 493817 2244.500
    3 1365938 8643.400 4101178 9127.800
    4 1654525 10367.900 4101178 9127.800
    5 687126 2850.700 1365938 8643.400
    6 1406842 11166.400 493817 2244.500
    7 1417151 5780.500 876226 4480.400
    8 1081233 8273.300 1540723 6095.800
    9 2508639 14697.700 736735 3776.100
    10 1081233 8273.300 7716.422 1453.800
    Source: http//amar.tavanir.org.ir//entaghl

     | Show Table
    DownLoad: CSV
    Table 9.  Transmitter level desirable and undesirable outputs and Dual—role factor.
    DMU Transmitter 1 (division 6) Transmitter 2 (division 7)
    w_{1k}^6g_{1k}^6 b_{1k}^7 w_{1k}^7g_{1k}^7 b_{1k}^7
    1 1592990 508.845 8681541.4 51.880
    2 1151302.3 200.566 183110 301.829
    3 7291961.5 175.381 11551302.3 357.789
    4 5661596 328.197 11551302.3 117.468
    5 330324 67.759 7291961.5 263.987
    6 559431.3 254.862 183110 107.780
    7 6151576.2 447.605 330747 61.919
    8 88601.2 373.774 479386 202.020
    9 8681541.2 273.358 231110 84.462
    10 88601.2 294.146 4261453.8 38.828
    Source: http//amar.tavanir.org.ir//entaghal and calculations loose of electricity

     | Show Table
    DownLoad: CSV
    Table 10.  The distributor level inputs.
    DMU Distributor 1 (division 8) Distributor 2 (division 9)
    \bar x_{1k}^8\tilde x_{1k}^8 \bar x_{2k}^8 \bar x_{1k}^9\tilde x_{1k}^9 \bar x_{2k}^9
    1 779247 40437 406754 60332
    2 11349292 64702 233061 19739
    3 11349292 64702 306879 28043
    4 861255 12406 178742 8942
    5 90036 13383 2480122 26770
    6 11349292 64702 317529 15731
    7 3639109 37153 1444115 13785
    8 208430 51688 422169 24689
    9 779247 40437 189471 18162
    10 269026 35606 208430 51688
    Source: http//amar.tavanir.org.ir//tozee

     | Show Table
    DownLoad: CSV
    Table 11.  The distributor level inputs.
    DMU Distributor 3 (division 10) Distributor 4 (division 11)
    \bar x_{1k}^{10}\tilde x_{1k}^{10} \bar x_{2k}^{10} \bar x_{1k}^{11}\tilde x_{1k}^{11} \bar x_{2k}^{11}
    1 332536 13761 449258 10052
    2 178742 18122 132419 11101
    3 3651115 32533 90036 13383
    4 187438 12075 317547 56184
    5 3965115 32533 306879 28043
    6 132419 11101 189471 18162
    7 90036 13383 11349292 64702
    8 406754 60332 539565 52340
    9 332536 13761 406754 60332
    10 406754 60332 539565 52340
    Source: http//amar.tavanir.org.ir//tozee

     | Show Table
    DownLoad: CSV
    Table 12.  Distributor level desirable and undesirable outputs.
    DMU Distributor 1
    (Division 8)
    Distributor 2
    (Division 9)
    g_{1k}^8 b_{1k}^8 g_{1k}^9 b_{1k}^9
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    576253
    2046151
    2046151
    1288350
    265678
    2046151
    497281
    294579
    576253
    469733
    14.210
    7.200
    15.570
    15.570
    13.250
    15.57
    13.600
    11.230
    14.210
    12.540
    576253
    323920
    631924
    345484
    662102
    513660
    429044
    368658
    513660
    347768
    8.030
    10.400
    11.390
    10.730
    12.670
    11.510
    11.050
    13.330
    7.250
    11.230
    Source:http//amar.tavanir.org.ir//tozee

     | Show Table
    DownLoad: CSV
    Table 13.  Distributor level desirable and undesirable outputs.
    DMU Distributor 3
    (Division 10)
    g_{1k}^{10}
    b_{1k}^{10} Distributor 4
    (Division 11)
    g_{1k}^{11}
    b_{1k}^{11}
    1 248079 13.590 327034 14.200
    2 345484 10.730 208346 7.990
    3 429044 11.050 265678 13.250
    4 329071 7.670 309704 12.030
    5 429044 11.05 631924 11.390
    6 208346 7.990 333449 7.250
    7 265678 13.25 2046151 15.570
    8 550244 8.030 691491 8.100
    9 208346 13.590 631924 8.030
    10 550244 8.030 691491 8.100
    Source: http//amar.tavanir.org.ir//tozee

     | Show Table
    DownLoad: CSV
    Table 14.  The customer level inputs.
    DMU Customer 1
    (Division 12)
    \bar x_{1k}^{12}
    Customer 2
    (Division 13)
    \bar x_{1k}^{13}
    Customer 3
    (Division 14
    \bar x_{1k}^{14}
    Customer 4
    (Division 15)
    \bar x_{1k}^{15}
    1 1400 1094.800 1096.400 2802.500
    2 1400 1094.800 1096.800 2802.500
    3 1400 1094.800 1096.800 2802.500
    4 1400 1094.800 1096.800 2802.500
    5 1400 1094.800 1096.800 2802.500
    6 1400 1094.800 1096.800 2802.500
    7 1400 1094.800 1096.800 2802.500
    8 1400 1094.800 1096.800 2802.500
    9 1400 1094.800 1096.800 2802.500
    10 1400 1094.800 1096.800 2802.500
    Source: http//amar.tavanir.org.ir//tozee

     | Show Table
    DownLoad: CSV
    Table 15.  The customer level inputs under managerial disposability.
    DMU Customer1
    (Division 12)
    \tilde x_{1k}^{12}
    Customer2
    (Division 13)
    \tilde x_{1k}^{13}
    Customer 3
    (Division 14
    \tilde x_{1k}^{14}
    Customer 4
    (Division15)
    \tilde x_{1k}^{15}
    1 0.000 258.173 30.0710 7195.787
    2 0.000 4.89300 28.7950 68.90600
    3 0.000 38.7860 22.8310 3564.162
    4 0.000 4.89300 74.6070 6801.258
    5 0.000 0.00000 17.8910 2024.679
    6 0.000 0.00000 310.5440 2241.095
    7 0.000 0.00000 0000.000 1276.555
    8 0.000 112.4370 0000.000 6377.373
    9 0.000 258.1730 0000.000 4747.578
    10 0.000 61.16200 141.2120 218.9860
    Source: http//amar.tavanir.org.ir//tozee

     | Show Table
    DownLoad: CSV
    Table 16.  The Customer level desirable and undesirable outputs.
    DMU Customer 1 (division 12) Customer 2 (division 13)
    g_{1k}^{12} g_{2k}^{12} b_{1k}^{12} g_{1k}^{13} g_{2k}^{13} b_{2k}^{13}
    1 1830958 6122.147 778.277 347030 3241.136 147.510
    2 6441756 5485.296 725.081 1778416 2903.980 200.178
    3 7866277 5821.292 725.323 2168359 3081.860 199.937
    4 6560395 4865.888 727.327 1791210 2576.059 198.585
    5 3804176 3622.099 752.559 855850 1917.582 169.308
    6 8009286 3996.064 734.466 2078242 2115.563 190.588
    7 8271676 5563.775 693.427 2196721 2945.528 184.154
    8 3602333 6217.991 718.110 962150 3291.877 191.801
    9 3213868 3906.777 752.079 691239 2068.293 161.757
    10 3683518 3635.504 722.771 953080 1924.679 187.011
    Source: http//amar.tavanir.org.ir//tozee and calculations time cut off of electricity

     | Show Table
    DownLoad: CSV
    Table 17.  The Customer level desirable and undesirable outputs.
    DMU Customer 3 (division 14) Customer 4 (division 15)
    g_{1k}^{14} g_{2k}^{14} b_{1k}^{14} g_{1k}^{15} g_{2k}^{15} b_{1k}^{15}
    1 16364 2700.947 6.956 7663 5942.083 3.257
    2 37745 2419.983 4.249 57685 5323.964 6.492
    3 51444 2568.217 4.743 65030 5650.077 5.996
    4 37480 2146.715 4.155 53509 4722.774 5.932
    5 42460 1597.985 8.400 28981 3515.567 5.733
    6 45458 1762.970 4.169 73999 3878.533 6.786
    7 624532 2454.607 52.355 72330 5400.135 6.064
    8 106646 2743.231 21.259 24231 6035.109 4.830
    9 54540 1723.578 15.103 30174 3791.871 7.061
    10 110055 1603.899 21.595 23562 3528.578 4.623
    Source: http//amar.tavanir.org.ir//tozee and calculations time cut off of electricity

     | Show Table
    DownLoad: CSV

    The material flow or intermediate measures from suppliers divisions to manufacturers divisions, from manufactures divisions to the transmitters divisions and from transmitters divisions to distributors divisions and from them to the customers divisions are presented in appendix Tables 1823. Tables 24 and 25 of Appendix indicates inverse intermediate measures to exit from manufactures divisions and enter to suppliers divisions, exit from transmitters divisions and enter manufactures divisions. The division's weights and the overall weights of the 15 divisions are presented in Table 26.

    Table 18.  The inefficiency scores of supply chains (DMUs).
    DMU \, {\theta _o} \xi _k^{S1} \xi _k^{S2} \xi _k^{M2} \xi _k^{M2} \xi _k^{M3} \xi _k^{T1} \xi _k^{T2} \xi _k^{D1} \xi _k^{D2} \xi _k^{D3} \xi _k^{D4} \xi _k^{C1} \xi _k^{C2} \xi _k^{C3} \xi _k^{C4}
    1 0.005 0 0 0 0 0 0 0 0.006 0 0.16 0 0 0 0 0
    2 0.13 0 0 0 0 0 0.42 0 0 0.25 0.22 0 0.33 0 0.15 0.30
    3 0.20 0 0.41 0 0 0.92 0 0 0 0.50 0.19 0 0.36 0 0.28 0.36
    4 0.16 0 0 0 0 0.47 0.56 0 0 0 0.36 0 0.23 0 0.09 0.23
    5 0.10 0 0 0 0 0.63 0 0 0 0.58 0.19 0.58 0 0 0 0
    6 0.05 0 0.35 0 0 0 0 0 0 0 0 0.33 0 0 0.08 0
    7 0.15 0 0.21 0 0 0 0 0.4 0 0.6 0 0 0.32 0.32 0 0.32
    8 0.09 0 0 0 0 0 0 0 0 0 0 0.56 0.24 0 0 0.22
    9 0.04 0 0 0 0 0 0 0 0 0.58 0 0 0 0 0 0
    10 0.03 0 0 0 0 0 0 0 0.28 0 0 0.56 0 0 0 0

     | Show Table
    DownLoad: CSV
    Table 19.  The dual variables of inputs and undesirable output of residential division.
    DMU Dual variable of inputs under Natural disposability
    v
    Dual variable of inputs under Managerial disposability
    z
    Dual variable of desirable output
    U1
    Dual variable of desirable output
    U2
    DTR Effective of investment
    1 0.00000162 0.00000162 0.000000044239 0.00000000 P --
    2 0.00000162 0.00000162 −0.000000001466 0.000016489 N E
    3 0.00000162 0.00000162 −0.000000001407 0.000015815 N E
    4 0.00000162 0.00000162 −0.000000001682 0.000018914 N E
    5 0.00000162 0.00000162 0.0000000022301 0.0000020021 P --
    6 0.00000162 0.00000162 −0.000000002194 0.00000024667 N E
    7 0.00000162 0.00000162 −0.000000001492 0.00000016777 N E
    8 0.00000162 0.00000162 0.00000001082 0.00000067567 P --
    9 0.00000162 0.00000162 0.00000001432 0.00000089453 P --
    10 0.00000162 0.00000162 0.000000008172 0.0000014056 P --

     | Show Table
    DownLoad: CSV
    Table 20.  The dual variables of inputs and undesirable output of public division.
    DMU Dual variable of inputs under Natural disposability
    v
    Dual variable of inputs under Managerial disposability
    z
    Dual variable of desirable output
    U1
    Dual variable of desirable output
    U2
    DTR Effective of investment
    1 0.0000018 0.0000000055 0.0000002074 0.00000000 P --
    2 0.0000018 0.0000000055 0.00000002155 0.00001159 P --
    3 0.0000018 0.0000255 −0.00000001251 0.00003216 N E
    4 0.0000018 0.0000000055 0.00000003009 0.000024116 P --
    5 0.0000018 0.0000000055 0.00000003009 0.000024116 P --
    6 0.0000018 0.0000000055 −0.00000006928 0.00004084 N E
    7 0.0000018 0.0000000055 −0.00000004747 0.000027984 N L
    8 0.0000018 0.0000000055 0.0000065392 0.0000027592 P --
    9 0.0000018 0.0000000055 0.0000055462 0.000016276 P --
    10 0.0000018 0.000001523 0.00000007132 0.00003387 P --

     | Show Table
    DownLoad: CSV
    Table 21.  The dual variables of inputs and undesirable output of agriculture division.
    DMU Dual variable of inputs under Natural disposability
    v
    Dual variable of inputs under Managerial disposability
    z
    Dual variable of desirable output
    U1
    Dual variable of desirable output
    U2
    DTR Effective of investment
    1 0.0000015 0.0000000038642 0.0000036666 0.0000000 P --
    2 0.0000015 0.00014888 0.00000062738 0.000001500 P --
    3 0.0000015 0.000012614 0.0000005315 0.0000012715 P --
    4 0.0000015 0.000016028 0.0000006742 0.0000016157 P --
    5 0.0000015 0.000000074749 0.000017686 0.0000017686 P --
    6 0.0000015 0.00013692 0.0000013199 0.0000000000 P --
    7 0.0000015 0.00000003842 −0.00000006977 0.000042196 N L
    8 0.0000015 0.000000038642 0.00000065261 0.000000000 P --
    9 0.0000015 0.00000003760 0.00000061798 0.000011671 P --
    10 0.0000015 0.000084827 0.00000054518 0.000000000 P --

     | Show Table
    DownLoad: CSV
    Table 22.  The dual variables of inputs and undesirable output of industrial division.
    DMU Dual variable of inputs under Natural disposability
    v
    Dual variable of inputs under Managerial disposability
    z
    Dual variable of desirable output
    U1
    Dual variable of desirable output
    U2
    DTR Effective of investment
    1 0.00002017 0.00000000024 0.000011353 0.0000000 P --
    2 0.00002017 0.00000000024 −0.000000055 0.00001694 N L
    3 0.00002017 0.000016557 −0.00000001281 0.000016872 N E
    4 0.00002017 0.000010272 0.0000003085 0.000014925 P --
    5 0.00002017 0.000002201 0.000000777 0.00001833 P --
    6 0.00002017 0.00000000024 −0.0000002137 0.000026509 N L
    7 0.00002017 0.00000000024 −0.0000001456 0.000018061 N L
    8 0.00002017 0.00032940 0.000015013 −0.00004586 N E
    9 0.00002017 0.0000024261 0.0000010953 0.000014228 P --
    10 0.00002017 0.0000000002 0.0000018587 0.000012245 P --

     | Show Table
    DownLoad: CSV
    Table 23.  Effective and limited investment opportunity on supply chain 10 in the power.
    Effective investment Percent% Limitedinvestment Percent%
    Residential 5 0.5 0 0.0
    Public 2 0.2 1 0.10
    Agriculture 0 0.1 1 0.10
    Industrial 5 0.5 3 0.30

     | Show Table
    DownLoad: CSV

    We now describe the results obtained in the new proposed approach. The model (10) is applied to estimate the efficiency score of supply chain 10 (DMUS). The model(10) is solved by a linear programming solver using the GAMS software on a 8GB RAM, 2.0 GHz desktop computer, the runtime of the computation in this study is negligible in model. The results are listed in Table 27.

    The first column of Table 18 represents the global inefficiency score of the supply chains. It can be easily seen that no DMU can reach inefficiency equal to null. This implies that all the 10 supply chains can improve their performance in some of the divisions. Supply chain number 1 is the one that reaches the lowest inefficiency score (0.005) while supply chain number 3 is the worst performing one. Looking vertically in the tables, the more efficient divisions are divisions 1, 3, 4 and with efficient values (100% of the total). This implies that supplier 1, manufacturers 2 and 3 are the more efficient ones concerning the other divisions. Just one efficient unit (90%) is obtained in the case of divisions 7 (Transmitter 2).

    As an illustration, we consider the four divisions of residential, public, agriculture and industrial from power consumers to identify DTR measures and effective investment on customer divisions in supply chains 10 of the power industry. This study applies the proposed radial model to examine the sustainability performance of supply chains. In the first stage, the supply chain management enables according to the dual variables sign related to desirable output constraints determent if the investment to new technology innovation decrease undesirable output or no then in the next stage, if the dual variables sign was negative then the decision-maker define the type input under management disposability from customer divisions in supply chain as their increase has an effective or limited effect to abatement of undesirable outputs. In other word, if the dual variable sign is positive then the investment does not effective to decrease of undesirable outputs. Tables 19-22 indicate the dual variables of the optimal solution of model (11) related to inputs under natural and managerial disposability and desirable outputs in residential, public, agriculture and industrial divisions of power subscribers on electricity companies 10 in different regions of Iran.

    According to Table 19, u_1^* < 0 for supply chains number (2, 3, 4, 6, 7) of residential division and they belonged to negative DRT, so indicating technology innovation was essentially necessary for enhancing their efficiency and sustainability. Moreover, the amount of the dual variable of input under managerial disposability determines the type of investment on inputs. Other word, z_{}^* > W\varepsilon R_q^x for supply chains number (2, 3, 4, 6, 7) then input for investment under managerial disposability can effectively decrease the number of undesirable outputs and is rated as E (Effective investment). Therefore, they have a high potential for an investment opportunity to enhance the entire supply chain sustainability. Similarity, supply chains number (1, 5, 8, 9, 10) belonged to positive DTR so indicating that the technological innovation was not an essential necessary for increase their unified efficiency and sustainability improvements in performance assessment.

    Similarity, according to Table 20 supply chains number (3, 6, 7) of public division of power customer have u_1^* < 0 hence negative Damage to return (N) occurs on division public of supply chains so this case indicate technology innovation was essentially necessary for enhancing their efficiency and sustainability. In other word, the amount z_1^* > W\varepsilon R_1^x for supply chains number (3, 6) hence input for investment under managerial disposability can effectively decrease the number of undesirable outputs and is rated as E (Effective investment). Similarity, supply chains number 7 belonged to negative DTR and z_{}^* = W\varepsilon R_q^x so the input for investment under managerial disposability has a limited effect on decrease some of undesirable outputs because z_{}^* is a very small positive number so the investment has only a limited effect and it rated as L (limited investment). Finally, supply chains number (1, 2, 4, 5, 8, 9, 10) belonged to positive damage to return (P) so indicating that the technological innovation was not an essential necessary for increasing their unified Efficiency in performance assessment.

    According to Table 21 supply chains number 7 of agriculture division have u_1^* < 0 and belonged to negative DRT and z_1^* = W\varepsilon R_1^x then the input under managerial disposability for investment has a limited effect on decrease some of undesirable outputs and it rated as L (limited investment) and the other supply chains belonged to positive DTR so indicating that the technology innovation was not an essential necessary for increasing their unified Efficiency.

    Finally, according to Table 22 supply chains number (2, 3, 6, 7) of industrial division have u_1^* < 0 and supply chain number 8 has u_2^* < 0 and negative DRT, so indicating technology innovation was essentially necessary for enhancing their efficiency and sustainability. In other words, the amount z_{}^* > W\varepsilon R_q^x for supply chains number 3 and 8 hence the input for investment under managerial disposability can effectively decrease the number of undesirable outputs and is rated as E (Effective investment). Therefore, they have a high potential for an investment opportunity to enhance the entire supply chain sustainability. Similarity, supply chains numbers (2, 6, 7) belonged to negative DTR and z_{}^* = W\varepsilon R_q^x therefore, input for investment has a limited effect on decrease a number of undesirable outputs and it rated as L (limited investment). Finally, supply chains number (1, 4, 5, 9, 10) belonged to positive DTR so technological innovation was not an essential necessary for increase their unified efficiency and sustainability improvements in performance assessment.

    Table 23 summarizes effective and limited investment opportunity on ten supply chains of four division of consumers in the power industry, all of them are specified by DTR and they are classified into two investment categories (effective and limited investment).

    As summarized in the Table 23 the industrial division of the power consumers had depicted a high level of effective (0.50) and limited investment (0.30) opportunity. Therefore, the industrial sector may have a high potential for an investment opportunity to enhance the entire sustainability. Finally, it is worth noting the energy and industrial sectors are the most attractive investment regions for enhancing sustainability and efficiency in production processes.

    The electricity supply chain is a network of energy sectors, power production divisions, transmission and distribution lines, and power subscribers. The power industry is one of the important investment targets for reducing wasted energy in oil and gas fields and power plant sectors and enhancing corporate sustainability. Furthermore, the investment to decrease the power losses in transmitter and distributer lines is an immediately necessary to increasing operational and environmental efficiency. This study proposes a model radial to a supply chain sustainability assessment which measures an investment opportunity for technology innovation and decreasing the number of undesirable outputs in the different sectors of the supply chain. Also, technology innovation in the energy and industrial sectors not only prevents energy losses but also abatement global warming and climate changes. It is immediately necessary to know whether the investment to undesirable outputs abatement can effectively decrease a number of undesirable outputs or increase the inputs under managerial disposability have a limited effect on decrease a number of undesirable outputs. In other words, an important feature of the proposed approach is that it able to identify the investment on which division of supply chain has a major or minor impact in decrease a number of undesirable outputs. Besides, it is possible to increase an input under managerial disposability may be ineffective in undesirable produces reduction.

    This study has two empirical results of customer divisions. One of the two results is that the transmission and distribution companies must have adequate decisional capacities regarding investment for transmitting directly the power to industrial, agriculture divisions and sectors of high electricity consumption in the power industry. Particular, the residential and industrial divisions have significant capacities on investment and technology innovation for reducing undesirable outputs to achieve corporate sustainability in the supply chain. The other result is that the dual-role factors have an important key role in the handling of undesirable output in the energy sector and the abatement inner the electricity consumptions of power plants in electricity production sector and managing specialist workforce to decrease of losses power in transmitter lines. Moreover, they able to enhance the effectiveness of transmission and distribution lines in the network supply chain. In general, all of the studied researches about environmental performance assessment of supply chain do not consider network DEA model based on investment to new technology innovation and undesirable products reduction. Also, the proposed model is able to handling investment on capital assets to the pollution emissions abatement and the power losses in the electricity supply chain divisions. Indeed, the difference between the proposed model and other approaches is that the model is able to recognize increase which the categories inputs cause significant decrease in wasted energy and harmful emissions. The proposed approach has three methodological limitations in leading environmental performance assessment. First, the source energy is different among districts. Each region has its essential structure and different conditions for business activity. For instance, southern regions in Iran have noticeable energy sources and the high capacity of power plants respect to other regions. Such regional difference effects on the number of efficiency measures in each regional. Second, the proposed approach assumes that all unified efficiency measures are uniquely determined on optimality. If the uniqueness assumption of efficiency measures is dropped, the proposed model needs to incorporate strong complementary slackness conditions into the model to obtain a unique optimal solution. The assumption on uniqueness is appropriate to the measurement of the dual variable by the model (11). Third, this study has not considered many companies in the proposed DEA assessment that contain a negative value on production indexes. The problem considered in this study needs to further researches in future. Similarity, this study can be conducted for green supply chain management evaluation in a time horizon by Malmquist index computation on time-series data to examine the frontier shift among multiple periods.

    The author would like to thank the anonymous reviewer and editor whose constructive comments have improved the quality of this study.

    There is no conflict of interest in this paper.



    [1] Fan Y, Bai B, Qiao Q, et al. (2017) Study on eco-efficiency of industrial parks in chain based on data envelopment analysis. J Environ Manage 192: 107-115. doi: 10.1016/j.jenvman.2017.01.048
    [2] Sueyoshi T, Gotto M (2017) Measurement of return to scale on large photovoltaic power station in the United States and Germany. Energy Econ 64: 306-320. doi: 10.1016/j.eneco.2017.03.028
    [3] Wang Z, Zhang J, Wei YM (2017) Operational and environmental performance in China's thermal power industry: Taking an effectiveness measure as complement to an efficiency measure. J Environ Manage 192: 254-270. doi: 10.1016/j.jenvman.2017.01.066
    [4] Zhang J, Liu Y, Chang Y, et al. (2017) Industrial eco-efficiency in China: A provincial quantification using three stage data envelopment analysis. J Clean Prod 143: 238-249. doi: 10.1016/j.jclepro.2016.12.123
    [5] Raili K, Markku H, Pekka O (2018) Methanol-managing greenhouse gas emissions in the production chain by optimizing the resource base. AIMS Energy 6: 1074-1102.
    [6] Gereme GD, Gorge C, Jose DL, et al. (2020) Economic evaluation methodologies for renewable energy project. AIMS Energy 8: 339-364. doi: 10.3934/energy.2020.2.339
    [7] Ahmed N, Amira M (2020) Local energy management in hybrid electrical vehicle via fuzzy rules system. AIMS Energy 8: 421-437. doi: 10.3934/energy.2020.3.421
    [8] Colver F, Sueyoshi T (2009) Contribution of Professor William W Cooper in operation research and management science. Eur J Oper Res 196: 1-16. doi: 10.1016/j.ejor.2008.02.014
    [9] Sueyoshi T, Goto M (2010) Measurement of a linkage among environmental, operation al and financial performance in Japanese manufacturing firms: A use of data envelopment analysis with Strong complements slackness condition. Eur J Oper Res 207: 1742-1753. doi: 10.1016/j.ejor.2010.07.024
    [10] Sueyoshi T, Gotto M (2010) Should the US clean air act include CO2 emission control? Examination by data envelopment analysis. J Energy Policy 38: 5902-5911.
    [11] Sueyoshi T, Goto M (2012) Methodological comparison between two unified (operational and environmental) efficiency measurements for environmental assessment. Eur J Oper Res 210: 684-693.
    [12] Sueyoshi T, Gotto M (2014) Environmental assessment for corporate sustainability by resource utilization and technology innovation: DEA radial measurement on Japanese industrial sectors. Energy Econ 46: 295-307. doi: 10.1016/j.eneco.2014.09.021
    [13] Kao C (2009) Efficiency decomposition in network data envelopment analysis (A relational model). Eur J Oper Res 192: 949-962. doi: 10.1016/j.ejor.2007.10.008
    [14] Tone K, Tsutsui M (2010) An epsilon-based measure of efficiency in DEA-A third pole of technical efficiency. Eur J Oper Res 207: 1554-156. doi: 10.1016/j.ejor.2010.07.014
    [15] Khalili D, Taghavefard M, Olfat L, et al. (2011) A hybrid approach based on fuzzy DEA and simulation to measure the efficiency of agility in supply chain: Real case of dairy industry. Int J Manage Sci Eng Manage 6: 163-172.
    [16] Toloo M (2014) Selecting and full ranking supplier with imprecise data: A new DEA method. Int J Adv Manuf Technol 74: 1141-1148. doi: 10.1007/s00170-014-6035-9
    [17] Tavana M, Mirzagoltabar H, Mirhedayatian SM, et al. (2013) A New Network Epsilon-based DEA model for supply chain performance evaluation. J Comput Ind Eng 66: 501-513. doi: 10.1016/j.cie.2013.07.016
    [18] Mhdiloo M, Sean RF, Lee KH (2015) Technical environmental and eco-efficiency measurement for supplier selection: An extension and application of data envelopment analysis. Int J Prod Econ 168: 279-289. doi: 10.1016/j.ijpe.2015.07.010
    [19] Tajbakhsh A, Hassini E (2015) A Data Envelopment Analysis approach to evaluate sustainability in supply chain networks. J Clean Prod 105: 74-85. doi: 10.1016/j.jclepro.2014.07.054
    [20] Khodakaremi M, Shabani A, Farzipor sean R, et al. (2015) Developing distinctive two-Stage data envelopment analysis models: An application evaluating the sustainability of supply chain management. Measurement 70: 62-74. doi: 10.1016/j.measurement.2015.03.024
    [21] Devika K, Jafarian A, Norbakhsh V (2014) Designing a Sustainability Closed- loop supply chain network based on triple bottom line approach: A comparison of meta heuristics hybridization technique. Eur J Oper Res 235: 594-615. doi: 10.1016/j.ejor.2013.12.032
    [22] Nikfarjam H, Rostami MM, Mamezadeh CS (2015) Measuring supply chain efficiency based on a hybrid approach. Transportation Research part D: Transp Environ 39: 141-150. doi: 10.1016/j.trd.2015.06.004
    [23] Babazadeh R, Razmi I, Rabbani M, et al. (2017) An integrated data envelopment analysis mathematical programming approach to strategic biodiesel supply chain network design problem. J Cleaner Prod 147: 694-707. doi: 10.1016/j.jclepro.2015.09.038
    [24] Pouralizadeh M, Amirteimoori A, Riccardi R, et al. (2020) Supply chain performance evaluation in the presence of undesirable outputs: A case on power industry. AIMS Energy 8: 48-80. doi: 10.3934/energy.2020.1.48
    [25] Toolo M (2009) On classifying inputs and outputs in DEA: A revised model. Eur J Oper Res 198: 358-360. doi: 10.1016/j.ejor.2008.08.017
    [26] Hatefi SM, Jolia F (2010) A new model for classifying inputs and outputs and evaluating the performance of DMUs based on trans log outputs distance function. Appl Math Modell 34: 1439-1449. doi: 10.1016/j.apm.2009.08.017
    [27] Farzipoor Saen R (2010) Developing a new data envelopment analysis methodology for supplier selection in the presence of both undesirable outputs and imprecise data. J Int Adv Manuf Technol 51: 1243-1250. doi: 10.1007/s00170-010-2694-3
    [28] Farzipoor Saen R (2010) Restricting weights in supplier selection decision in the presence of dual-role factors. J Appl Math Modell 34: 2820-2830. doi: 10.1016/j.apm.2009.12.016
    [29] Mirhedayatian SM, Farzipoor Sean R, Azadi M (2014) A Novel Network Data Envelopment Analysis model for evaluation green supply chain management. Int J Prod Econ 147: 544-554. doi: 10.1016/j.ijpe.2013.02.009
    [30] Available from: http://amar.tavanir.org.ir/pages/report/stat94/sanatebargh/sanatebargh
  • energy-08-04-580-S01.pdf
  • This article has been cited by:

    1. Mojgan Pouralizadeh, Sustainability Assessment of Electricity Supply Chain via Resource Waste Reduction and Pollution Emissions Management: A Case Study of the Power Industry, 2021, 5, 25206478, 20210004, 10.1520/SSMS20210004
    2. Vítor João Pereira Domingues Martinho, 2021, Chapter 3, 978-3-030-75838-7, 85, 10.1007/978-3-030-75839-4_3
    3. Mohammad Tavassoli, Saeedeh Ketabi, Mahsa Ghandehari, A novel fuzzy network DEA model to evaluate efficiency of Iran’s electricity distribution network with sustainability considerations, 2022, 52, 22131388, 102269, 10.1016/j.seta.2022.102269
    4. Mohammad Tavassoli, Mahsa Ghandehari, Classification and forecasting of sustainable-resilience suppliers via developing a novel fuzzy MIP model and DEA in the presence of zero data, 2023, 1936-9735, 10.1007/s12063-023-00401-z
    5. Mojgan Pouralizadeh, Marginal profit maximization estimation of supply chains by waste energy decrement: a case study of the power industry, 2024, 58, 0399-0559, 3143, 10.1051/ro/2024103
    6. Shanshan Chen, Ruchuan Zhang, Peiwen Li, Aijun Li, How to improve the performance of China's energy-transport-economy-environment system: An analysis based on new strategy parallel-series input-output data envelopment analysis models, 2023, 281, 03605442, 128235, 10.1016/j.energy.2023.128235
    7. Jian Zhang, Jianzhou Wen, Zhen Lu, Jiang Qian, Ning Wei, Investment Allocation Method for Distribution Networks Based on a Panel Data Model and an Incentive–Penalty Mechanism, 2025, 27, 15829146, 656, 10.24818/EA/2025/69/656
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3560) PDF downloads(277) Cited by(7)

Figures and Tables

Figures(2)  /  Tables(23)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog