Optimal Operation of Multiple Energy System Based on Multi-Objective Theory and Grey Theory

: The manufacturing industry consumes electricity and natural gas to provide the power and heat required for manufacturing. Additionally, large amounts of electric energy and heat energy are used, and the electricity cost, amount of environmental pollution, and equipment maintenance cost are high. Thus, optimizing the management of equipment with new energy is important to satisfy the load demand from the system. This paper formulates the scheduling problem of these multiple energy systems as a multi-objective linear regression model (MLRM), and an energy management system is designed focusing on the economy and on greenhouse gas emissions. Furthermore, a variety of optimization objectives and constraints are proposed to make the energy management scheme more practical. Then, grey theory is combined with the common MLRM to accurately represent the uncertainty in the system and to make the model better reﬂect the actual situation. This paper takes load ﬂuctuation, total grid operation cost, and environmental pollution value as reference standards to measure the effect of the gray optimization algorithm. Lastly, the model is applied to optimize the energy supply plan and its performance is demonstrated using numerical examples. The veriﬁcation results meet the optimized operating conditions of the multi-energy microgrid system.


Introduction
With the acceleration of industrialization and urbanization, the energy crisis, air pollution, and other problems have become increasingly serious [1]. In order to reduce environmental pollution and to improve power quality, the traditional power grid is gradually transforming into a smart grid [2]. A microgrid can be developed as an effective means for distributed power sources to connect to a smart grid [3]. A microgrid is composed of a distributed power supply, energy storage, and load and has unique advantages in improving the utilization rate of renewable energy. It can reduce the power interaction with a superior power grid, can alleviate the impact on a superior power grid [4], and plays an important role in lowering carbon emissions and in improving economic benefits.
At present, many scholars have conducted preliminary research on optimizing the operation of a microgrid. The authors of Reference [5] researched the optimization of operations of a microgrid based on the chaotic particle swarm optimization algorithm. The authors of Reference [6] combined the optimization of the operation of a microgrid with the game analysis method and proposed an optimal configuration model of an intermediate microgrid in the distribution network based on game theory. The authors of Reference [7] optimized the operation of a microgrid by modifying the hypermutation operator, which is the key mechanism of the CLONALG algorithm. Based on the non-dominated sorting genetic algorithm (NSGA-II), an optimal solution set under all objectives was found in [8] and used in the operation planning of a microgrid. However, the above multi-objective optimization method has weak dynamic response ability, and it is difficult to solve the microgrid optimization problem with real-time parameters in practical application. Although literature [9] adopts grey theory to optimize the model and improve the dynamic response speed, its optimization only targets economy and does not consider environmental problems such as exhaust emissions during micro-grid operation.
The grey multi-objective programming algorithm is a dynamic programming algorithm. The parameters in the model contain grey numbers to make up for the deficiency in general linear programming. It not only knows the optimal configuration under established conditions and dynamic conditions but also makes the model better reflect the objective reality [10]. Therefore, an intelligent operation control strategy of a multi-energy microgrid based on the grey multi-objective programming algorithm is proposed in this paper to effectively improve the economic costs and to reduce environmental pollution. First, mathematical models of the economic operation cost, exhaust emission, and equipment maintenance cost of a multi-energy microgrid are established. Second, a mathematical model of the power supply, heat supply, and cold supply balance is established. Then, a grey multi-objective linear programming solution model is established based on the mathematical model to economically and environmentally optimize the operation of a multi-energy microgrid. Finally, the effectiveness of the proposed strategy is verified using simulations.
The main contributions of this article are as follows: (1) Combining wind power, photoelectric power, an energy storage system, and a gas system, the energy management system is designed with a focus on the economy and on greenhouse gas emissions. (2) This paper proposes a grey multi-target linear planning algorithm and optimizes multitarget multi-energy management using the grey multi-target linear planning algorithm.
As Figure 1 shows, this article studies a multiple energy microgrid system model. This model was established according to the MES system architecture of an area in Liaoning Province. The solid line shows the flow of electricity, the two-way arrow shows the two-way flow of electricity, the dotted line shows the flow of heat energy, and the thin dotted line is the signal transmission line. As shown in Figure 1, the power generation information of all power generation equipment and the power consumption information of power consumption equipment, the heat provided by gas, electricity price, and fuel cost are transmitted to the SCADA control system. By using the control algorithm proposed in this paper, the control instructions of each unit are obtained and transmitted to the power generation and heating equipment. Figure 1 shows that, in addition to the power grid, photovoltaic cells, wind turbines, energy storage batteries, and the micro gas turbine can also supply electricity to the system. The heat is supplied by a gas turbine engine, an electric boiler, and thermal storage as the main sources. The cost of the system includes the electricity and gas costs. The equipment used to detect environmental pollution mainly considers the micro gas turbine and electric boiler system as well as the equivalent environmental pollution produced when generating power for the power grid. Additionally, all power generation and energy storage devices in the system need to be in the best state [11,12]. The following section discusses the objective function and the constraint conditions of a multiple energy microgrid system model for optimizing the operation of its equipment.

Objective Function
As shown in Figure 1, the micro grid system model includes minimizing the system economic operation cost, maintenance cost, and considering the environmental benefits, minimizing emissions of the system. Set the formulation as follows: 1. The energy consumption costs 2. The waste gas emissions

The components maintenance costs
The output of the wind turbines and photovoltaic cells is without artificial restrictions, so the major variable of the operation maintenance cost function are battery, electric boiler, gas turbine, and regenerator. The formulation is

Objective Function
As shown in Figure 1, the micro grid system model includes minimizing the system economic operation cost, maintenance cost, and considering the environmental benefits, minimizing emissions of the system. Set the formulation as follows: 1.
The energy consumption costs 2.
The waste gas emissions 3.

The components maintenance costs
The output of the wind turbines and photovoltaic cells is without artificial restrictions, so the major variable of the operation maintenance cost function are battery, electric boiler, gas turbine, and regenerator. The formulation is

1.
The electricity supply and demand constraint This equation sets out the harming system disturbance rate at 50% as benchmark, only counting the disturbance rate that are higher than the benchmark.

6.
Cost of micro gas turbine maintenance 7.
Cost of thermal storage maintenance where is the number of electric heating boiler.

1.
Constraint of heat output of the electric boiler 2.

Constraint of heat output of the thermal storage
The heat capacity constraint equation of the thermal storage system is Time continuous function of heat storage of the thermal storage system is where µ h is the heat loss coefficient of the thermal storage. Set the initial state equation for the thermal storage system as follow: Heat supply constraint equation of the thermal storage system is 3. Constraint of heat output of the gas fired boiler 4. Constraint of electricity output of the battery SOC state constraint equation of the battery system is where S is the battery SOC lower limit index. Setting its calculating formula as the ratio of the micro grid system's important load and the battery rated capacity, in order to ensure that even if the power grid failed suddenly, the important load would not suffer the impact. Time continuous function of SOC state of the battery system is Set the initial state equation for the battery system as follows: Electricity supply constraint equation of the battery system is

5.
Constraint of output of the micro gas turbine The electricity output equation of the gas turbine system is Constraint of electricity output of the gas turbine system is The heat output equation of the gas turbine system is Simplified, considering the gas turbine used for heat and electricity production as a basic proportional relation, the ratio is a coefficient of thermal efficiency.
The fuel consumption equation of the gas turbine system is

Grey Multi-Objective Linear Programming Algorithm
In grey multi-objective linear programming algorithm, whiting the grey number first is often necessary to solve the equation [13,14]. Note grey number a(⊗) = [a, a], then the whitened value of a(⊗) is a(⊗). Note that grey information of a(⊗) is a(⊗) = [a, a], and the whitened value a(⊗) can set as a(⊗) ∈ [a, a] [15].
The grey multi-objective linear programming mathematics model of system is where C(⊗) is grey target matrix, A(⊗) and Aeq(⊗) are grey distribution matrix, while B(⊗) and Beq(⊗) are grey constraint matrix.
By consuming coefficient of the grey value to whiten, the expression is where , the corresponding grey multi-objective linear programming solution is S(θ c , θ b , θ a ), and the feasible solution domain is R(θ b , θ a ).

Model Solving Method
The whitened optimal solution matrix is an n-order zero matrix, and the whitened probability integral matrix is an m-order zero matrix. Then, the technical parameters of the power generation and energy storage modules are denoted A and Aeq, respectively, and the thermal and electric load values are denoted B and Beq, respectively. The corresponding coefficient matrix of the grey threshold can be used to present a variable range of the parameters. Set the corresponding whitened coefficients of Aeq and Beq as a and b, Aeq ∼ a ∈ [0, 1] and Beq ∼ b ∈ [0, 1] [16]. The whiting formulae are provided in (40) and (41). The whitened coefficient and corresponding coefficient are iteratively calculated. The algorithm for the flow of model substitution is shown in Figure 2. (41). The whitened coefficient and corresponding coefficient are iteratively calculated. The algorithm for the flow of model substitution is shown in Figure 2.
Output the global optimal solution min(min( )) ( , ) Read electric load value、model parameters、Initial value of model coefficients C, A, B (including gray field value) Output the optimal solution matrix of the corresponding uptake coefficient S1 Step 1: Initialization. Whiten the optimal solution matrix. Whiten the probability integral matrix, which is a [1/n, 1/m] order zero matrix. Obtain the technical parameters of the power generation and energy storage modules. Obtain the model constraint matrices A and Aeq, and the thermal and electric load values, and obtain the model constraint ma- Step 1: Initialization. Whiten the optimal solution matrix. Whiten the probability integral matrix, which is a [1/n, 1/m] order zero matrix. Obtain the technical parameters of the power generation and energy storage modules. Obtain the model constraint matrices A and Aeq, and the thermal and electric load values, and obtain the model constraint matrices B and Beq. If the parameters are in the variable range, the corresponding coefficient matrix is represented by the corresponding grey field value.
Step 3: Cycle the whitening uptake coefficient step by step to solve the optimal solution of multi-objective programming for the corresponding whitening uptake coefficient, and output the whitening optimal solution matrix S1 at the end of the cycle.
Step 4: Using the grey multi-objective linear programming model of theorem 1 and theorem 2, for any point A and point B of the corresponding solution plane, if 0 ≤ µ a1 < µ a2 ≤ 1 and 0 ≤ λ b2 < λ b1 ≤ 1, then S(θ b1 , θ a1 ) ⊂ S(θ b2 , θ a2 ). Therefore, the probability function is assigned to the whitening uptake coefficient, and the probability integral is solved to output the probability integral matrix.
Step 5: Weight the solution of the matrix to obtain the optimal solution for the probability integral that optimizes the global optimal solution under double constraints.

Example Analysis
First, consider the grey multi-objective linear programming model of the objective function values, for the purposes of this example, of the standard time-sharing electricity charge for industrial enterprises according to the rules and relevant regulations of electricity prices in the power supply business, according to the transformer capacity (including highpressure motors without any transformer capacity) or the maximum demand. Large data on the 10 kv industrial electricity prices from Liaoning province are shown in Table 1. As the gas price is not fixed, the online query gas prices are set according to the data from 24 September 2021 based on the Shenyang gas price at 3.48 yuan/cubic meter.
The statistics for the constraint function in the model, including the environmental impact coefficient and equipment capacity configuration, are shown in Table 2 [17][18][19]. Set the initial state of batteries as SOC(t 0 ) = K 0 = 1/30, and set the initial state of thermal storage as q sr (t 0 ) = −10. The micro gas turbine power supply constraint is [0, 80], the grey threshold of the heating coefficient is [0.4, 0.6], the boiler system power supply constraint is [01,20], and the grey threshold of the thermal storage's heat loss coefficient is [0.7, 0.7]. The data for wind turbines, photovoltaic cells, heat load, and electric load have normal distributions within the scope of the threshold. Thus, the function for the corresponding whitened coefficient has a normal distribution. They can be initialized as Use the grey multiple linear regression algorithms based on weather information to forecast the amount of energy generated by the wind turbines and PV battery. The forecasting result's error is under 10%. Therefore, the constraint matrix of its grey domain can be set as ± 10% of the forecasting result, which is shown in Table 3.
The result is as follows: Table 3. Grey domain settings of P w , P s , e load , q h . Due to the industrial enterprises' load being of a single wave form, during 21-24 PM, the load is in the low range, so the energy storage unit during this time period cannot self-prepare for the next 24 h and loses part of its energy storage function. Therefore, in order to achieve maximum energy-saving effects, the optimization time period is adjusted forward by 4 h. In order to conveniently read the results shown below, other than the energy storage equipment, the result is still shown with the x-axis starting at zero.

Simulation Result
First, in the condition m = n = 0.02, the optimal solution matrix of the whitened values and the corresponding probability distribution value of the whitened coefficient can be calculated, as shown three-dimensionally in Figure 3.

Simulation Result
First, in the condition m = n = 0.02, the optimal solution matrix of the whitened val and the corresponding probability distribution value of the whitened coefficient can calculated, as shown three-dimensionally in Figure 3.   To achieve the optimal solution after comprehensive consideration, the objecti functions need to be as economical as possible while meeting most energy constraints. T optimal solution of the resulting surface is shown in Figure 5. The optimal solution of t probability integral is 85.44%, and the corresponding points are S (0.02, 0.34) = 0.02. How ever, that point and the optimal solution value do not have physical meanings. The pla with all equipment running is shown as follows, and the electrical load requirements an their grey domain values are shown in Figure 5. To achieve the optimal solution after comprehensive consideration, the objective functions need to be as economical as possible while meeting most energy constraints. The optimal solution of the resulting surface is shown in Figure 5. The optimal solution of the probability integral is 85.44%, and the corresponding points are S (0.02, 0.34) = 0.02. However, that point and the optimal solution value do not have physical meanings. The plan with all equipment running is shown as follows, and the electrical load requirements and their grey domain values are shown in Figure 5. functions need to be as economical as possible while meeting most energy constraints. The optimal solution of the resulting surface is shown in Figure 5. The optimal solution of the probability integral is 85.44%, and the corresponding points are S (0.02, 0.34) = 0.02. However, that point and the optimal solution value do not have physical meanings. The plan with all equipment running is shown as follows, and the electrical load requirements and their grey domain values are shown in Figure 5.  Figure 6 shows an electricity supply plan for microgrid power. When the electric load is low, the period of the energy generated by wind power is high and the storage battery is charged to full capacity in preparation for a peak load. Due to the microgrid coordinating energy, a significant peak sharpening effect is seen, where 48% of the peak electricity  Figure 6 shows an electricity supply plan for microgrid power. When the electric load is low, the period of the energy generated by wind power is high and the storage battery is charged to full capacity in preparation for a peak load. Due to the microgrid coordinating energy, a significant peak sharpening effect is seen, where 48% of the peak electricity is provided by the microgrid and 50% of the flat electricity is provided by other distributed units of the system.            Figure 9 presents the micro gas turbine electricity supply plan, and the micro gas turbine's main power period is the peak load period for electricity and heat, all while running under maximum capacity. During this period, the cost of using power generated by the gas turbine is lower than that from the power grid. Thus, the gas turbine helps supply heat energy.
Energies 2022, 15, x FOR PEER REVIEW 13 of 23 Figure 9 presents the micro gas turbine electricity supply plan, and the micro gas turbine's main power period is the peak load period for electricity and heat, all while running under maximum capacity. During this period, the cost of using power generated by the gas turbine is lower than that from the power grid. Thus, the gas turbine helps supply heat energy.  Figure 10 shows the performance of the battery system, with the negative part indicating that the battery is in a charged state and the positive part indicating that the battery is in a discharged state. As seen, at 09:00, 10:00, 13:00, and other times when battery energy is required, the battery is in a discharged state, and the discharged power reaches a peak suppression effect when the battery reaches its peak, while the battery is in a charged or idle state during the remaining time.
Of note, the SOC of the battery is measured at the initial moment of each time period, so the SOC of each time period is determined according to the charge and discharge of the battery from the previous time. In contrast, the positive part of the black curve indicates that the battery is in a charged state while the negative part indicates that the battery is in a discharged state.   Figure 10 shows the performance of the battery system, with the negative part indicating that the battery is in a charged state and the positive part indicating that the battery is in a discharged state. As seen, at 09:00, 10:00, 13:00, and other times when battery energy is required, the battery is in a discharged state, and the discharged power reaches a peak suppression effect when the battery reaches its peak, while the battery is in a charged or idle state during the remaining time.
Of note, the SOC of the battery is measured at the initial moment of each time period, so the SOC of each time period is determined according to the charge and discharge of the battery from the previous time. In contrast, the positive part of the black curve indicates that the battery is in a charged state while the negative part indicates that the battery is in a discharged state.
The SOC of the battery is shown in Figure 11.
idle state during the remaining time. Of note, the SOC of the battery is measured at the initial moment of each time period, so the SOC of each time period is determined according to the charge and discharge of the battery from the previous time. In contrast, the positive part of the black curve indicates that the battery is in a charged state while the negative part indicates that the battery is in a discharged state. The SOC of the battery is shown in Figure 11.   Figure 13 shows the heat supply plan for the heating boiler system. The duration of heat supplied is 24 h a day and it is a little higher during the peak period, while no big fluctuations occur.   Figure 13 shows the heat supply plan for the heating boiler system. The duration of heat supplied is 24 h a day and it is a little higher during the peak period, while no big fluctuations occur.  Figure 13 shows the heat supply plan for the heating boiler system. The duration of heat supplied is 24 h a day and it is a little higher during the peak period, while no big fluctuations occur.  Figure 13 shows the heat supply plan for the heating boiler system. The duration o heat supplied is 24 h a day and it is a little higher during the peak period, while no big fluctuations occur.   Figure 14 shows the micro gas turbine for the system that provides heat, and the heating load for the grey value is shown in the figure by applying the global optimal to the heating load and by using the visible micro gas turbine as the main heat source.
Energies 2022, 15, x FOR PEER REVIEW 15 Figure 14 shows the micro gas turbine for the system that provides heat, and heating load for the grey value is shown in the figure by applying the global optima the heating load and by using the visible micro gas turbine as the main heat source. In Figure 15, the heat supply plan of the thermal storage system can be seen in hot trough period in the heat storage condition. When run, it instantaneously rele heat, and during the day, it charges and discharges energy many times. The remaining heat in thermal storage is shown in Figure 16. In Figure 15, the heat supply plan of the thermal storage system can be seen in the hot trough period in the heat storage condition. When run, it instantaneously releases heat, and during the day, it charges and discharges energy many times.
Energies 2022, 15, x FOR PEER REVIEW Figure 14 shows the micro gas turbine for the system that provides heat, a heating load for the grey value is shown in the figure by applying the global opt the heating load and by using the visible micro gas turbine as the main heat sourc In Figure 15, the heat supply plan of the thermal storage system can be seen hot trough period in the heat storage condition. When run, it instantaneously r heat, and during the day, it charges and discharges energy many times. The remaining heat in thermal storage is shown in Figure 16.  The remaining heat in thermal storage is shown in Figure 16. The remaining heat in thermal storage is shown in Figure 16.  Considering the corresponding equipment for energy supply and demand relations, and each period when the main equipment supplies the load power, a detailed analysis of the electricity supply and demand relationship is shown in Figure 17 and a detailed analysis of the heat supply and demand relationship is shown in Figure 18.

Energies 2022, 15, x FOR PEER REVIEW
Considering the corresponding equipment for energy supply and demand rel and each period when the main equipment supplies the load power, a detailed an of the electricity supply and demand relationship is shown in Figure 17 and a de analysis of the heat supply and demand relationship is shown in Figure 18.  In Figure 17, the system load includes electricity used for lighting, ventilation Considering the corresponding equipment for energy supply and demand relations and each period when the main equipment supplies the load power, a detailed analysi of the electricity supply and demand relationship is shown in Figure 17 and a detailed analysis of the heat supply and demand relationship is shown in Figure 18.  In Figure 17, the system load includes electricity used for lighting, ventilation, office buildings, and air conditioning; electricity production; and electricity from the electri  In Figure 17, the system load includes electricity used for lighting, ventilation, office buildings, and air conditioning; electricity production; and electricity from the electric heating boiler. The supply of electricity is from the grid, a miniature gas turbine, wind turbines, photovoltaic cells, as well as the battery. During 20:00-7:00, which is a trough in the electricity load period, the power is mainly provided by wind turbines, while during 8:00-19:00, which is the peak load period, power is mainly provided by the power grid and a gas turbine. Compared with the single-energy-source supply plan, the battery plays a very important role in coordinating energy.
As seen in Figure 18, thermal storage plays an important role in coordinating heat. During 20:00-6:00, the electric heating boiler is the only heat supply source. During 7:00-19:00, which is the peak heat load period, the heat load is mainly provided by the gas turbine, the electric boiler, and thermal storage. Therefore, the gas turbine supplies the majority of heat during the peak hot load period. The use of thermal storage avoids frequent alterations in the gas turbine and electric boiler, which can increase the heating efficiency and the service life of the equipment.
The running results of two optimization algorithms and the algorithm without distributed energy resource and optimization are shown in Figure 19. From these results, the calculated results are shown as being able to optimize the system. Since methods have advanced and the economic costs of this method are discussed, a comparison and an analysis of various kinds of optimization methods are shown in Table 4.
In common multiple linear regression algorithms, there are multiple objective functions and certain contradictions between objectives. In the case, there is no optimal solution of the usual meaning. The solution of the multi-objective planning problem can be divided into three situations: (1) If at least one solution enables all the objectives of the multiobjective planning problem to be achieved, the solution is called the optimal solution for the multi-objective planning; (2) If at least one solution can only meet part of the importance objectives in the multi-objective planning problem, the solution is called a suboptimal or feasible solution for multi-objective planning; (3) If solution that can meet any goal of the multi-target planning problem cannot be found, then the multi-target planning problem has no solution. At this time, it is shown that the existing constraints of the multi-target planning problem are contradictory, and reasonable and feasible goals need to be reset for solution.

Environment Pollution
Value (g/kWh) Figure 19. The running results of three algorithms. There are many solutions for common multiple linear regression algorithms, such as the ideal point method, sum of square weighting method, linear weighting method, priority method, minimum maximum method, etc. This paper takes the priority method as the solution method of the common multiple linear regression algorithms. The basic idea of the method is to divide into different priorities according to the importance between different targets. First seeking the optimal value of the high priority objective function, on the premise of ensuring that the high priority target, and then seeking the optimal value of the low priority target function, the optimal solution is the multi-objective planning problem (satisfactory solution).
From Table 4, forecasting the cost for one day as an example, when the system has no distributed energy resource and optimization, the total energy consumption cost is 223.6 USD, the environment pollution value is 1.97 × 10 6 , and the total cost of a day has a lot of room for improvement.
Using common multi-objective linear programming algorithms, because new energy resources are replacing the traditional power grid, the total cost is reduced by a lot and environmental pollution is also greatly reduced. However, for system uncertainty, this method is useful when the load of the system is low. Therefore, the method is adopted.
Using the proposed grey multi-objective linear regression optimization algorithms to obtain the global optimal solution, the new total cost is 27.1 USD, in which electricity makes up 120.9 USD, the gas fee is 58.6 USD, and the environment pollution value is 1.7 × 10 6 , compared with the cost when optimization is not used, 223.6 USD; the amount of savings from daily optimization, 50.3 USD; and the reduction in environmental pollution, 2.7 × 10 5 .
The optimization results in Table 4 mainly reflect the constraints between energy supply and demand. The optimization results are not only affected by the energy supply and demand, but also by the time-sharing electricity price. The rise in electricity price will increase the total cost of power grid operation, but after optimizing the system with a gray multi-target linear planning method, the cost-saving cost will increase. This is due to the optimized system, and its utilization rate of electricity is greatly improved, thus saving more power. The electric energy saved is calculated into the cost saving according to the electricity price at this time, which will ultimately increase the cost saving.
In order to highlight the characteristics of the method used in this paper, we compared the results with those from recent literature. Taking Reference [20] as an example, the following parts are included in the comparison: This paper proposes a distributed distribution feeder load management method based on an autonomous decision-making entity that not only relieves the system pressure, but also ensures customer satisfaction by optimizing the demand constraint limit to alleviate the demand rebound effect. Compared with the results of a method based on the water injection algorithm and the simple and fair distribution method, the method designed is verified as having a more efficient performance in alleviating demand rebound and in reducing transformer congestion.
However, in order to meet the requirements of the model solution or to facilitate calculation, the equipment parameters and load prediction data involved in the above article are regarded as ideal data, but this does not reflect reality. As can be seen from the optimization results of this paper, the changes in parameters and load requirements greatly affect the final optimization results. Therefore, the grey multi-target linear planning algo-rithm proposed in this paper can obtain relatively global optimal results when considering uncertainties or uncertain parameters in an actual situation.
In the process of optimizing the operation of a multi-energy microgrid in this paper, the corresponding energy management was conducted from the perspective of the economy and exhaust emissions. At the same time, based on relevant data from actual industrial enterprises, a multi-energy microgrid model integrating wind power, photoelectric power, an energy storage system, and a gas system is established, and the optimized constraint function and control target of a multi-energy microgrid are determined, which is more in line with the actual situation in a multi-energy microgrid.

Conclusions
In this paper, the management of a multiple energy system in the manufacturing industry is considered and balancing the optimization of optimal economic goals, optimal environmental targets, and use of optimal equipment is attempted. First, the scheduling problem is formulated as multiple linear regression algorithms. Second, grey theory is applied to the formulation to solve uncertainty factors. Third, improved grey multiple linear regression optimization algorithms are applied to optimize all of the power supply devices and energy storage equipment. Note that the scheduling problem is continuously solved and that only the equipment's capacity is implemented. This could be affected by adjusting the device's capacity, especially the storage device.
The following points are summaries of the innovations in this study: (1) Combining wind power, photoelectric power, an energy storage system, and a gas system, the energy management system is designed focused on the economy and on greenhouse gas emissions. Considering the actual operation of a multi-energy microgrid system, in order to make the energy management scheme more practical, a variety of optimization objectives and constraints are proposed in this paper. The optimization objectives and constraints are determined not only by taking into account the economic cost of micro-network operations, but also by mathematically modeling the exhaust emissions when the system is run and by setting the corresponding energy control targets in terms of economic benefits and environmental protection. (2) For the possible uncertainty resulting when constructing the multi-objective model and the constraints in practical engineering cases, the grey multi-objective linear planning algorithm is proposed. Using the grey multi-target linear planning algorithm, the multi-target multi-energy optimization management of a microgrid is realized. By comparisons with multiple optimization methods, the economic costs of the proposed method are verified and the amount of cost savings and the cost recovered using our method are analyzed. The verification results meet the optimized operating conditions of a multi-energy microgrid system after analysis.
The method proposed in this paper can realize the optimal operation of multi-energy microgrid systems, but some shortcomings are also present. The focus of the optimal operation of the multi-energy microgrid system proposed in this paper is on small and medium-sized industrial enterprises. Due to their own electric load and thermal load characteristics, a typical winter day is used as the data source for the load samples. This paper did not consider the use of a cooling load in high temperatures during summer, but considering recent trends in global warming, with the temperature increasing in summer, even industrial enterprises need to be equipped with corresponding refrigeration equipment to ensure that production activities progress as usual. Therefore, in a follow-up research, we should consider the impact of a summer cooling load on the optimal operation of a multi-energy microgrid system and should increase the accuracy and reliability of the microgrid optimization operation.
At the same time, the methods proposed in this paper are based on theoretical research and the focus is also limited to the power load, with obvious peak and trough periods for industrial enterprises. Certain practical engineering cases are also lacking to verify the practicability of our methods. In the future, further development of the software is needed to apply it to an actual energy efficiency management system, to constantly improve the flexibility and accuracy of the algorithm, and finally to invest in complex microgrids with various types of loads to contribute to green energy savings.