Economic, Energetic, and Environmental Impact Evaluation of the Water Discharge Networks from Mining Works

: Sustainable development represents an optimistic scenario for the evolution of contemporary civilization. The object of this paper is to deﬁne certain evaluation criteria regarding the performances of water discharge networks from mining works, and propose a method for aggregating the speciﬁc indicators for monetary costs, energetic costs, and environmental impact-related costs. The global pollution index (GPI) represents a method for assessing environment health status or pollution levels. The GPI quantitatively expresses this status based on its index, which results from a ratio between the ideal value and the given value of certain quality indices that are considered speciﬁc for the analyzed environmental factors at certain moments. The proposed method in this paper tries to perform a synergistic aggregation of the balance sheet of harmfulness and classic balance sheets for matter and energy for an industrial process.


Introduction
The sustainable development postulates refer to economic efficiency, social solidarity, and ecological responsibility. The pursuit of the three postulates ensures the harmonization of the basic components of human community: natural capital, built capital, and social capital.
Economic efficiency, energetic efficiency, and reduced environmental impact are the desired results of the three principles. In the general economic and energetic context regarding environmental issues, optimization of the water transport systems for the purpose of reducing power consumption has become, currently, a very big problem. Besides the functional and energy optimization of the hydraulic system, the problem requires the elaboration of top-notch algorithms that are capable of solving complex problems regarding water transport and distribution networks.
Starting from Le Goff and Giuletti's question, "Thermodynamics and Economics-Is there a conflict?" [1], we aim in this paper to develop the CAREC diagram (energy-cost characteristic) and demonstrate the correlation between economic and energy efficiency as a basis for the sustainable development of the energy field. In this respect, we add a third axis to the CAREC diagram-pollution-as a dimension that will integrate the economic point of view (profit driven) with an energetic point of view (efficiency driven) and an environmental impact point of view (social responsibility). The algorithm proposed in this paper will provide a way to assess the sustainability of new development projects in mining (specifically for water networks). Based on our approach, various industries can determine their appropriate values for this algorithm in order to adapt it to their specific to adapt it to their specific needs. Our paper expands the concerns of the two authors to the ecological component, introducing the GPI method for the adequate quantification of the environmental impact.

About Energy and Sustainable Development
In the context of today's world-characterized by continuous changes in the economic, social, and natural environment (fierce competition, resource crisis, conflict at the level of the individual, the self against armed conflicts, natural disasters)-a complex approach to sustainable energy production and use is required.
In over 30 years, sustainable development has evolved from a mere concept to a widespread change in attitude, behaviors, and mentalities. What was once extrapolated only in theory, academic reports, diplomatic conventions, norms, and scientific works has turned into practical, everyday realities of life [2][3][4][5][6][7][8][9][10][11][12]. This is an important topic that addresses life from the perspective of future generations. The concept of sustainability has become a strategic part of the economy, social endeavors, and environmental projects, as well as an intrinsic part of local and regional policies [13].
The evolution of this "process"-as Rosalie Gardiner (2002) described it-turned to a steep ascent curve in recent years. The process concerned the heads of state and governments, as well specialists from different domains or non-profit organizations, who were all trying to incorporate in the development policies the definition of sustainable development: "meets the needs of present generations without compromising the ability of future generations to meet their own needs" [14].
Thus, the "Europe 2030" project [15] has energy efficiency as one of its main objectives, transposed into a series of national documents, such as: the National Sustainable Development Strategy of Romania 2013-2020-2030 [16], and the Romanian Energy Strategy 2016-2030, with perspectives for 2050 [17] (Figure 1). "The energy's trilemma"-as energy problems have been called lately-represents a challenge for the development of energy supply that must be solved in close interdependence and correlation ( Figure 2) [18]. We consider that in some industrial areas-such as mining, for instance-improvements can be made. "The energy's trilemma"-as energy problems have been called lately-represents a challenge for the development of energy supply that must be solved in close interdependence and correlation ( Figure 2) [18]. We consider that in some industrial areas-such as mining, for instance-improvements can be made.

Specific Costs and Energetic Indices for the Execution of Hydraulic Networks
First, we must take into account that hydraulic networks use materials and technologies that require a significant quantity of energy, and for the entire duration of the exploitation of the network, energy is required for the pumping of water. As a result of the conditions of current energy penury and its extremely high cost, the optimum network solution cannot be obtained unless optimization criteria are met, within which a significant role is held by the specific cost and energy indices for the execution and exploitation of these networks. The experimental base was constituted by the mine evacuation networks in Jiu Valley, which consists of networks of different sizes and different ages. Experimental determinations have been achieved with the help of specific devices applied to hydraulics: an ultrasonic flowmeter, manometric sensors, and environmental quality monitoring systems. The research method was based on analytical relations, experimental measurements, and the use of qualitative performance indicators, according to the latest information from the literature of the analyzed domain.

Cost Index
Based on statistical data and the comparative situation of the cost indices obtained during the execution of water discharge networks from mining works, the average values were determined for the specific investment costs in the networks' pipes, whose variation with diameter is graphically presented in Figure 3.

Specific Costs and Energetic Indices for the Execution of Hydraulic Networks
First, we must take into account that hydraulic networks use materials and technologies that require a significant quantity of energy, and for the entire duration of the exploitation of the network, energy is required for the pumping of water. As a result of the conditions of current energy penury and its extremely high cost, the optimum network solution cannot be obtained unless optimization criteria are met, within which a significant role is held by the specific cost and energy indices for the execution and exploitation of these networks. The experimental base was constituted by the mine evacuation networks in Jiu Valley, which consists of networks of different sizes and different ages. Experimental determinations have been achieved with the help of specific devices applied to hydraulics: an ultrasonic flowmeter, manometric sensors, and environmental quality monitoring systems. The research method was based on analytical relations, experimental measurements, and the use of qualitative performance indicators, according to the latest information from the literature of the analyzed domain.

Cost Index
Based on statistical data and the comparative situation of the cost indices obtained during the execution of water discharge networks from mining works, the average values were determined for the specific investment costs in the networks' pipes, whose variation with diameter is graphically presented in Figure 3.

Specific Costs and Energetic Indices for the Execution of Hydraulic Networks
First, we must take into account that hydraulic networks use materials and technologies that require a significant quantity of energy, and for the entire duration of the exploitation of the network, energy is required for the pumping of water. As a result of the conditions of current energy penury and its extremely high cost, the optimum network solution cannot be obtained unless optimization criteria are met, within which a significant role is held by the specific cost and energy indices for the execution and exploitation of these networks. The experimental base was constituted by the mine evacuation networks in Jiu Valley, which consists of networks of different sizes and different ages. Experimental determinations have been achieved with the help of specific devices applied to hydraulics: an ultrasonic flowmeter, manometric sensors, and environmental quality monitoring systems. The research method was based on analytical relations, experimental measurements, and the use of qualitative performance indicators, according to the latest information from the literature of the analyzed domain.

Cost Index
Based on statistical data and the comparative situation of the cost indices obtained during the execution of water discharge networks from mining works, the average values were determined for the specific investment costs in the networks' pipes, whose variation with diameter is graphically presented in Figure 3.  Steel pipes were taken into account, as they have the highest level of operational safety. However, they also have the disadvantages of a shorter lifespan and additional cost: with the internal and external corrosion protections these pipes require, they are not cheaper than cast iron pipes.
Using the integral rational regression method, based on the chart, a general relationship results for the establishment of an average value for the specific investment cost for a certain diameter:

Energy Index for the Included Energy (Energy)
The water discharge networks imply important consumptions of fuel and other forms of energy, on one hand for the production and mounting of pipes, and on the other hand, for the operation of the pumping installations for water circulation in the pipes.
That is why it is necessary to both rationalize and reduce the power consumption within the existing water discharge systems-especially in the new ones, by design-as well as the power that is consumed for their execution.
The chart in Figure 4 compares the data of specific equivalent average power and fuel consumptions that occur during the manufacturing and mounting of steel pipes used for the water discharge, in accordance with their diameter.
The specific value for embedded energy "e", in relation with the pipes diameter "D", can be approximated using the integral rational regression, in a function with the following form: (2) Sustainability 2017, 9, 240 4 of 14 Steel pipes were taken into account, as they have the highest level of operational safety. However, they also have the disadvantages of a shorter lifespan and additional cost: with the internal and external corrosion protections these pipes require, they are not cheaper than cast iron pipes.
Using the integral rational regression method, based on the chart, a general relationship results for the establishment of an average value for the specific investment cost for a certain diameter:

Energy Index for the Included Energy (Energy)
The water discharge networks imply important consumptions of fuel and other forms of energy, on one hand for the production and mounting of pipes, and on the other hand, for the operation of the pumping installations for water circulation in the pipes.
That is why it is necessary to both rationalize and reduce the power consumption within the existing water discharge systems-especially in the new ones, by design-as well as the power that is consumed for their execution.
The chart in Figure 4 compares the data of specific equivalent average power and fuel consumptions that occur during the manufacturing and mounting of steel pipes used for the water discharge, in accordance with their diameter.
The specific value for embedded energy "e", in relation with the pipes diameter "D", can be approximated using the integral rational regression, in a function with the following form: (2) Figure 4. The influence of the diameter on the variation of the specific energy consumption included in the water discharge network's pipes [19].

The Effect of Partial Obstruction of the Flow Section on the Pressure Losses and of the Pipe's Transport Capacity Reduction
Operational energy losses in the water pipes are calculated based on the pressure losses.
The following are included in the total pressure losses (Δpf): linear losses (Δpl), and local losses (Δpξ): The linear (longitudinal) losses are determined with the Weisbach-Darcy equation: where: Figure 4. The influence of the diameter on the variation of the specific energy consumption included in the water discharge network's pipes [19].

The Effect of Partial Obstruction of the Flow Section on the Pressure Losses and of the Pipe's Transport Capacity Reduction
Operational energy losses in the water pipes are calculated based on the pressure losses.
The following are included in the total pressure losses (∆p f ): linear losses (∆p l ), and local losses (∆p ξ ): The linear (longitudinal) losses are determined with the Weisbach-Darcy equation: where: • λ represents the fluid dynamics resistance coefficient (Table 1) For the calculation of local losses the Borda-Carnot-type equations are used: where ξ is the coefficient for local loss. Pressure losses are influenced by the impurities deposited on the pipes' walls. By disconnecting certain pipes on sections of various diameters, the obstructed areas were measured (deposits of a mixture of water, oil, dust, and rust on a pipe wall).
For the purpose of introducing the obstruction effect on the longitudinal pressure losses, we propose to use a new term in the Weisbach-Darcy equation called "the thickness of the uniform obstruction", which is defined by the equation: where: • D i is the pipe's interior diameter, (m); • S tot is the area of the flow's cross-section to the new pipe (m 2 ); • S obt is the area of the area obstructed by sediments (m 2 ).
Although the space of water and obstructing deposits are in the shape of a circular section in the pipe's cross-section, in order to allow a quantitative assessment of the pipe's diameter reduction due to obstruction, the obstructed section shall be considered an annular shape, with the obstruction material being uniformly located on the pipe's wall and having the thickness δ.
The values of the uniform obstruction's thickness for various relative obstructions are given in Table 2. The obstructed area was calculated based on experimental measurements using the equation: where r is the pipe radius (m), and φ is the angle in the centre (degrees). For regular quality steel pipes, the status of the material in accordance with the k retardance is: The effect of the section reduction increases parabolically along with the decrease of the diameter. In the case of mining water discharge networks that have been in operation for five years, an absolute obstruction is found between 100-300 mm 2 , as shown in Figure 5. The obstructed area was calculated based on experimental measurements using the equation: where r is the pipe radius (m), and φ is the angle in the centre (degrees).
For regular quality steel pipes, the status of the material in accordance with the k retardance is: The effect of the section reduction increases parabolically along with the decrease of the diameter. In the case of mining water discharge networks that have been in operation for five years, an absolute obstruction is found between 100-300 mm 2 , as shown in Figure 5. Experiments highlight that an equal degree of pipe obstruction is translated into a relative percentage decrease, differentially, of the pipe's transport capacity, and respectively, in a relative variable increase of the pressure loss determined by the relative differential obstruction of the section for a series of diameters. The values obtained for an absolute obstruction of 200 mm 2 are presented in Figure 4, and centralized in Table 3. Experiments highlight that an equal degree of pipe obstruction is translated into a relative percentage decrease, differentially, of the pipe's transport capacity, and respectively, in a relative variable increase of the pressure loss determined by the relative differential obstruction of the section for a series of diameters. The values obtained for an absolute obstruction of 200 mm 2 are presented in Figure 4, and centralized in Table 3. Based on Figure 5 and Table 3, and after considering several real cases of water discharge pipes, and calculating the weighted averages according to the length of the majority of the pressure losses and transport capacity reduction due to partial obstruction of the flow section, I noticed that the operational energy loss due to pipe obstruction represented 3.5-5.5 % of the energy used for water circulation.

Energy Index for Operational Energy
The following are included under the name operational energy: the energy used from the power network, energy losses, and energy capitalization at the end user. In the case of water discharge networks from mining works, the operational energy represents the sum of the energy of discharged water (current), longitudinal and local energy losses in the network, additional losses due to obstruction, pump losses, drive engine losses, and losses on the power line that supplies the drive engine.
Taking into account the components mentioned previously, and correlating the experimental data with information in the literature in this field, an average estimated value was calculated for the specific operational energy consumption e sop (kWh/m 2 ) ( Figure 6).  Based on Figure 5 and Table 3, and after considering several real cases of water discharge pipes, and calculating the weighted averages according to the length of the majority of the pressure losses and transport capacity reduction due to partial obstruction of the flow section, I noticed that the operational energy loss due to pipe obstruction represented 3.5-5.5 % of the energy used for water circulation.

Energy Index for Operational Energy
The following are included under the name operational energy: the energy used from the power network, energy losses, and energy capitalization at the end user. In the case of water discharge networks from mining works, the operational energy represents the sum of the energy of discharged water (current), longitudinal and local energy losses in the network, additional losses due to obstruction, pump losses, drive engine losses, and losses on the power line that supplies the drive engine.
Taking into account the components mentioned previously, and correlating the experimental data with information in the literature in this field, an average estimated value was calculated for the specific operational energy consumption esop (kWh/m 2 ) ( Figure 6). The specific operational energy curve esop, in relation with the pipes' D diameter, can be approximated using integral rational regression, in a function of the form:

Environment Quality Index
The global pollution index represents a method of assessing the environment's health, or pollution status. It expresses this status from the quantity point of view, which is based on an index that has resulted from a ratio between the ideal value and the given value-at a specific moment-of certain quality indices that are considered specific for the analyzed environmental factors.
The GPI method involves going through several stages of synthetic assessment based on quality indices that might reflect the general status of one of the analyzed environmental factors, and then their correlation using a graphic method.
In this respect, it is proposed to enlist quality at a given moment of every environmental factor in a reliability scale, by granting marks that would express its closeness to, or distance from, the ideal status.
The reliability scale is expressed in grades from 1 to 10, with 10 representing a natural status unaltered by human activity, and 1 representing an irreversible and extremely serious deteriorated situation of the analyzed environmental factor.
Generally, it is considered that it is possible to assess the environment health in a certain area and at a certain given moment by: • Air quality; • Water quality; • Soil quality; • Population's health status; • Recorded deficit of species, plants, and animals (biodiversity).
Each of these factors can be characterized by several qualitative indices that are representative of the assessment of the pollution degree, and for which admissible limits are established. A reliability grade is granted according to the inclusion within the specified limits. The ideal status is graphically represented by a regular geometrical shape with radiuses equal to each other, and with the grade 10 of the reliability unity. By joining the points that resulted from the location of the values expressing the real status, an irregular geometrical shape is obtained with a smaller surface, inscribed in the regular geometrical shape of the ideal status.
The index of a global pollution status of an ecosystem-GPI-results from the ratio between the surfaces representing the ideal status; S I -and the surface representing the real status-S R .
When there are no changes in the quality of the environmental factors-i.e., when there is no pollution-this index is 1. Graphically, the geometrical shape illustrating the real status of the environment overlaps, and the figure illustrates the ideal status. When there are changes in the quality of the environmental factors, the index shall get higher, and show higher over unity values as the surface of the real triangle, square, or pentagon diminishes.
The reliability grades that are obtained for every environmental factor in the analyzed area contribute to the graphic execution of a chart as a method to simulate the synergic effect. The geometrical shape is an equilateral triangle when the three environmental factors are analyzed, a square when we have data for four environmental factors, and it can be a regular pentagon when five environmental factors are considered.
The alternative approach in this paper-the evaluation of impact on the environment-is carried out using a matrix method (quality indices) and the V. Rojanschi method (global pollution index) [21].
Applying the correlations between the quality indices and the reliability grades [21], the following values are obtained for the targeted environmental factors:

•
The value of the reliability grade for maximum air pollution (reliability grade 5), for water (reliability grade 4), for soil (reliability grade 4), for biodiversity (reliability grade 4), for population health (reliability grade 6); • The value of the reliability grade for minimum air pollution (reliability grade 9), for water (reliability grade 7), for soil-subsoil (reliability grade 8), for biodiversity (reliability grade 7), and for population health (reliability grade 8).
Based on the reliability grades granted to the considered environmental factors, a global pollution index chart was built where the following are represented: ideal status, minimum pollution status, and maximum pollution status (Figure 7). Sustainability 2017, 9,240 10 of 14 Applying the correlations between the quality indices and the reliability grades [21], the following values are obtained for the targeted environmental factors: • The value of the reliability grade for maximum air pollution (reliability grade 5), for water (reliability grade 4), for soil (reliability grade 4), for biodiversity (reliability grade 4), for population health (reliability grade 6); • The value of the reliability grade for minimum air pollution (reliability grade 9), for water (reliability grade 7), for soil-subsoil (reliability grade 8), for biodiversity (reliability grade 7), and for population health (reliability grade 8).
Based on the reliability grades granted to the considered environmental factors, a global pollution index chart was built where the following are represented: ideal status, minimum pollution status, and maximum pollution status (Figure 7). Calculating the afferent areas to the environment's status, the following values are obtained: • Area afferent to ideal status-S I = 240 cm 2 ; • Area afferent to maximum pollution status-S pol.max = 46.07 cm 2 ; • Area afferent to minimum pollution status-S pol.min = 136.965 cm 2 .
The value of the global pollution index for the two statuses overlapping the comparison:

Aggregation Method for the Economic, Energy and Ecologic Indices
In order to have the synergistic image of the economic, energy, and ecologic effects afferent to an industrial objective, the specialty literature uses various aggregation methods for performance indices.
Based on the CAREC chart (energy-cost characteristic) [23], which we adapted to the studied case, processed, and completed with the component regarding the environmental impact, we built the three-dimensional CAREEC chart (energy-ecology-cost characteristic), as shown in Figure 8.
The coordinates of the proposed chart are represented by: C M -monetary cost; E-total energy used; and P ol -environmental impact. The values of the specific costs, which correspond to coordinates, were expressed as a percentage of the maximum value in order to have a common measurement unit.
The values E Cmin , C Emin , and P olCmin correspond to 100%. The CAREEC chart enables the highlighting of various costs that are afferent to optimization measures based on economical, energy, and ecological criteria.

Aggregation Method for the Economic, Energy and Ecologic Indices
In order to have the synergistic image of the economic, energy, and ecologic effects afferent to an industrial objective, the specialty literature uses various aggregation methods for performance indices.
Based on the CAREC chart (energy-cost characteristic) [23], which we adapted to the studied case, processed, and completed with the component regarding the environmental impact, we built the three-dimensional CAREEC chart (energy-ecology-cost characteristic), as shown in Figure 8.
The coordinates of the proposed chart are represented by: CM-monetary cost; E-total energy used; and Pol-environmental impact. The values of the specific costs, which correspond to coordinates, were expressed as a percentage of the maximum value in order to have a common measurement unit.
The values ECmin, CEmin, and PolCmin correspond to 100%. The CAREEC chart enables the highlighting of various costs that are afferent to optimization measures based on economical, energy, and ecological criteria. As can be observed in Figure 8, the following correlations can be done:

•
The difference between ECmin-Emin shows the additional energy costs in order to reduce monetary costs; • The difference between EPolmin-Emin shows the additional energy costs in order to reduce environmental costs (afferent to pollution); • The difference between CEmin-Cmin shows the additional monetary costs in order to reduce energy costs; • The difference between CPolmin-Cmin shows the additional monetary costs in order to reduce environmental costs; As can be observed in Figure 8, the following correlations can be done: • The difference between E C min -E min shows the additional energy costs in order to reduce monetary costs; • The difference between E Polmin -E min shows the additional energy costs in order to reduce environmental costs (afferent to pollution); • The difference between C Emin -C min shows the additional monetary costs in order to reduce energy costs; • The difference between C Polmin -C min shows the additional monetary costs in order to reduce environmental costs; • The difference between P olEmin -P olmin shows the additional environmental costs in order to reduce energy costs; • The difference between P olCmin -P olmin shows the additional environmental costs in order to reduce monetary costs.
Analyzing the water discharge networks from three mines, the following average percentage values were obtained:

•
To achieve an energy saving of E Cmin − E min = 60%, the monetary costs are increased by 65%, and the environmental costs are increased by 35%; • To achieve an energy saving of E polmin − E min = 31%, the monetary costs are increased by 42%, and the environmental costs are increased by 28%; • To reduce environmental costs by P olCmin − P olmin = 46.5%, the monetary costs are increased by 49%, and the energy costs are increased by 31%.

Conclusions
Our proposed method is a synergistic aggregation of the harmfulness balance sheet and the classic balance sheets for matter and energy for an industrial procedure.
The difficulty of the method resides in the: • Collection of experimental data; • Definition of the system frontier to which the balance sheet applies. The conclusions differed for the analysis of an objective when the location is in a depression, a hilly area or a field area; • Definition of maximum allowed harmfulness; • Consideration of local environmental characteristics; • Harmonization of national and international regulations; • Achieving a compromise between regional and global interest.
The advantages offered by this approach are materialized via: • Using concepts and analytical tools from exergetic balance sheets; • Obtaining synthetic values that represent the global ecological load exercised by the analyzed system on its own environment; • The method allowing the comparison of different solutions regarding the execution of a product; • Providing information regarding the ecological carrying capacity of an ecosystem and attaining certain alarm thresholds; • With significant additions, it can be fitted into an analytical-experimental apparatus designed for the operationalizing of the sustainable development concept; • It provides the possibility of a periodical control of the real measurements' results regarding environmental protection, economic efficiency, and energy efficiency, which also facilitate their quantification.
The mining industry is also a primary energy producer (coal) and a large energy consumer of various types of energy, including: primary (hydraulic) and secondary (electric energy, compressed air). Therefore, for the sustainable development of this industry, various measures and interventions are necessary that use methods and tools to reduce and rationalize the consumption and minimize environmental impact.