当好一个吃瓜群众 《范伟打天下》家族门徒观战攻略

Abstract

To fully accommodate renewable and derivative energy sources in mine energy systems under supply and demand uncertainties, this paper proposes an optimized electricity–heat scheduling method for mining areas that incorporates Power-to-Gas (P2G) technology and Conditional Value-at-Risk (CVaR). First, to address uncertainties on both the supply and demand sides, a P2G unit is introduced, and a Latin hypercube sampling technique based on Cholesky decomposition is employed to generate wind–solar-load sample matrices that capture source–load correlations, which are subsequently used to construct representative scenarios. Second, a stochastic optimization scheduling model is developed for the mine electricity–heat energy system, aiming to minimize the total scheduling cost comprising day-ahead scheduling cost, expected reserve adjustment cost, and CVaR. Finally, a case study on a typical mine electricity–heat energy system is conducted to validate the effectiveness of the proposed method in terms of operational cost reduction and system reliability. The results demonstrate a 1.4% reduction in the total operating cost, achieving a balance between economic efficiency and system security.

1. Introduction

With the accelerated global transition towards a clean and low-carbon energy structure and the deepening implementation of China’s “dual carbon” strategic goals, the traditional mining industry, characterized by high energy consumption and high emissions, is confronted with unprecedented pressure to achieve energy conservation, emission reduction, and sustainable development [1,2]. Mining production activities, particularly large-scale open-pit mining and deep underground mining, are quintessential energy-intensive industries. These activities not only consume substantial amounts of electricity but also involve significant demands for thermal energy, such as for heating in mining areas and heat used in washing and sorting processes [3]. The satisfaction of energy demands within mining areas is highly reliant on fossil fuels, such as coal and natural gas. This dependency not only results in substantial operational costs but also constitutes a significant source of greenhouse gas emissions and environmental pollution. Consequently, the transition towards a cleaner, more efficient, and intelligent energy system in mining areas has become an imperative choice for the sustainable development of the industry.

During the production process, a substantial amount of derivative energy is generated within mining areas. In terms of their physical states, these derivative energies can be categorized into gaseous, liquid, and solid forms. The gaseous derivatives include exhausted air and coal mine methane, among others. The liquid derivatives consist of mine water inflow, etc. The solid derivatives include coal gangue, among others [4]. For the energy system of mining areas, the full utilization of derivative energy sources not only effectively reduces environmental pollution but also enhances overall energy utilization efficiency. Additionally, mining enterprises typically possess vast tracts of idle land (such as waste rock dumps, tailing ponds, and industrial plazas) and rooftop resources, which provide favorable conditions for the deployment of distributed photovoltaic (PV) and wind power generation systems. However, the inherent intermittency and fluctuation of renewable energy sources pose severe challenges to the stable and reliable energy supply in mines [5,6]. Mining production is characterized by its continuity and safety sensitivity, and thus, it demands extremely high stability and quality of power supply [7]. The effective mitigation of fluctuations in renewable power generation output, along with the high-proportion integration of renewable energy and its derivatives while achieving precise alignment with mining production loads, constitutes a fundamental challenge in establishing a green mining area’s electricity–heat integrated energy system (MIES).

In this context, Power-to-Gas (P2G) technology, as a highly promising flexibility resource and energy storage solution, offers a novel approach to addressing the aforementioned challenges [8,9]. P2G technology comprises two distinct phases. The first phase involves electrolyzing water to produce hydrogen using surplus electricity, particularly during off-peak periods or times of abundant wind/solar generation. In the second phase, the produced hydrogen is combined with carbon dioxide (CO₂) through methanation, thereby not only consuming electrical energy but also reducing CO₂ emissions. Consequently, P2G technology has gained widespread adoption.

Currently, scholars both domestically and internationally have conducted relevant research on MIES. Liang R et al. [10], considering the demand response capability of the entire production process in mining areas and the energy balance constraints of the system, proposed a high-quality coal mine energy system optimization model to achieve the full integration of renewable energy. Miao Q et al. [11] proposed a self-correcting optimization strategy for mining area energy systems with asymmetric prediction error fusion coefficients to enhance load forecasting accuracy, thereby effectively reducing system scheduling costs. Liu J et al. [12] proposed a multi-objective optimal coordination dispatch model for abandoned coal mine energy systems, incorporating P2G and power-to-hydrogen (P2H) technologies, aiming to minimize operational costs and carbon emissions. Wu F et al. [13] conducted a comprehensive review on utilizing underground coal mine spaces for energy storage, analyzing and discussing various energy storage technologies tailored for subsurface environments, along with their associated risks and challenges, with the aim of advancing the development of underground energy storage systems in abandoned mines. Xiong Y et al. [14] conducted a systematic review of advancements in green coal mining, encompassing intelligent green mining technologies and the green co-extraction of coal and its derivative resources, with a particular emphasis on the rational utilization of derivative resources as a critical research focus in mining area energy systems.

On the other hand, the inherent uncertainty in renewable energy output and load fluctuations poses significant risks to the dispatch decision-making and secure operation of MIES. Consequently, developing optimal dispatch strategies to mitigate or quantify the impact of these uncertainties on mining area energy systems, thereby ensuring both economic efficiency and operational security, has emerged as a critical research focus in academia. Currently, the predominant approaches for addressing uncertainties primarily include robust optimization [15,16], stochastic optimization [17,18], and fuzzy optimization [19,20], among others. Each method exhibits distinct characteristics and applicability in handling different types of uncertainty in mining area energy systems. Robust optimization accounts for the most extreme operating scenarios of the system, often yielding overly conservative dispatch results that compromise operational economic efficiency. Meanwhile, the selection of membership degrees in fuzzy optimization exhibits strong subjectivity. Stochastic optimization transforms uncertain problems into deterministic ones by generating discrete scenarios based on probability distributions of source–load variations. However, the computational burden increases significantly with a large number of discrete scenarios. Typically, scenario reduction techniques are required to cluster and retain representative scenarios for dispatch decision-making. It should be noted that the performance of traditional K-means algorithms depends heavily on the pre-specified number of clusters [21]. Moreover, multiple stochastic factors within the same region—including wind power, photovoltaic generation, and load demand—typically exhibit significant correlations. Xiang Y et al. [22] investigated the spatiotemporal correlation between wind and photovoltaic power to enhance PV utilization efficiency in complex mountainous regions. Ru Y et al. [23] established a probabilistic model of joint PV-wind power generation using copula functions. Subsequently, they employed Latin hypercube sampling combined with an improved K-means clustering technique to derive typical output scenarios. Zhong M et al. [24] utilized a MemoryFormer architecture to mine latent temporal correlations between reconstructed and enhanced wind-PV power output data.

However, the existing research exhibits several limitations: (1) None of the aforementioned studies quantitatively analyzed the economic risks induced by uncertainty factors. The Conditional Value-at-Risk (CVaR) [25], as a risk measurement metric, demonstrates superior advantages over Value at Risk (VaR) [26], including monotonicity, subadditivity, and accurate tail risk estimation. These properties enable CVaR to effectively evaluate operational risks in mining-area energy systems at specified confidence levels, thereby enhancing system security and reliability. (2) Existing studies on MIES have largely overlooked the utilization of derivative energy sources within the system, such as mine ventilation air methane (VAM) and mine water discharge. This oversight leads to significant resource wastage. Specifically, VAM refers to coal mine methane with a concentration below 0.75%, also known as coal mine ventilation air methane. Its long-term direct emission not only represents energy loss but also exacerbates the greenhouse effect due to methane’s high global warming potential. Therefore, MIES should incorporate technologies for VAM recovery and utilization to achieve both energy efficiency and environmental benefits.

Building upon existing research, this study develops an optimal dispatch strategy for mining-area energy systems by comprehensively considering P2G technology and source–load correlations. The proposed framework incorporates CVaR to quantitatively assess operational risks, thereby establishing a risk-aware optimization model for electricity–heat integrated energy systems in mining regions. The main contributions of this work are as follows:

• A Cholesky decomposition technique is employed to analyze source–load correlations in the mining-area integrated energy system (MIES). Building upon this correlation structure, Latin hypercube sampling (LHS) is implemented to generate a source–load sample matrix with predetermined correlation coefficients, followed by Affinity Propagation clustering to produce representative scenarios.
• With the objective of minimizing the total cost, which is constituted by dispatching cost, expected adjustment cost, and CVaR, a stochastic optimization dispatch model for the mining area energy system incorporating P2G technology is formulated. This model enables the reduction in operational costs of the mining area energy system and the enhancement of the renewable energy utilization rate while effectively addressing the impact of uncertain power outputs on both the supply and demand sides.
• A simulation study of the proposed optimal dispatch model was conducted on a typical mining area energy system to verify the model’s performance in terms of reducing operational costs and enhancing renewable energy accommodation, as well as the necessity of considering the correlation between supply and demand sides.

The rest of this paper is organized as follows: The structure of MIES containing P2G is analyzed in Section 2. In Section 3, the operation model of the unit in MIES is introduced. The model is simulated from several angles in Section 4, and the main conclusions are presented in Section 5.

2. The MIES Structure and Source–Load Risk Measurement

2.1. Analysis of the Architecture of the MIES Containing P2G

The architecture of the mining area’s electricity–heat integrated energy system incorporating P2G technology, as constructed in this study, is depicted in Figure 1.

In the MIES, the energy supply side consists of the external power grid, gas network, wind energy, photovoltaic (PV), mine ventilation air methane, and mine water discharge. The utilization equipment for derivative energy includes devices such as exhausted air heat storage oxidation units and water source heat pumps. The energy coupling equipment primarily comprises a micro gas turbine (MT) unit and a P2G system. On the demand side, the system caters to both electrical and thermal loads. In this system, the primary sources of uncertainty stem from the renewable energy sources on the supply side, namely wind energy and PV, as well as the loads on the demand side. Neglecting these uncertainties may increase the operational risks of the system, especially under extreme weather conditions, where the output of renewable energy sources may significantly decrease, thereby heightening the risk to energy supply security. Moreover, this can lead to the underutilization of renewable and derivative energy sources, thereby increasing the operational costs and CO₂ emissions of the system.

2.2. Generation of Typical Scenes on Both Sides of the Source and Load

The operation of the MIES is confronted with multiple uncertainties, including the fluctuating output of renewable energy sources and the variability of electrical load. Affected by natural environmental factors such as temperature and wind speed, these uncertain variables within the same region usually exhibit a certain degree of correlation. Based on this, this section employs the Spearman rank correlation coefficient to analyze historical data of wind speed, solar irradiance, and electrical load in a certain region of China. The resulting correlation heatmap is shown in Figure 2.

This study employs stochastic optimization to characterize the uncertainties of supply and demand. Stochastic optimization involves sampling based on the probability distributions of supply and demand to obtain discrete scenarios for handling uncertainties. Specifically, the probability models for the prediction errors of solar irradiance, wind speed, and electrical load are described using the beta distribution, Weibull distribution, and normal distribution, respectively. In terms of correlation control, the Cholesky decomposition method is a simple, practical, and widely adopted technique for managing the correlation structure among samples of random variables [27]. The Cholesky-based LHS method approximates the correlation matrix of the generated sample matrix to that of the historical data through procedures such as sampling and sorting. The detailed sampling steps can be found in reference [28]. Finally, 1000 samples of solar irradiance, wind speed, and electrical load are generated.

To avoid excessive computational burden caused by generating a large number of samples, clustering algorithms are commonly employed to reduce the sample set. The selected clustering algorithm should strive to preserve the correlation structure among the random variables within the generated samples. Affinity Propagation (AP) clustering is an unsupervised clustering algorithm based on message-passing mechanisms for cluster formation [29]. It is capable of adaptively determining the final number of clusters, and its resulting cluster centers correspond to actual data points within the sample set. Compared with traditional clustering algorithms such as K-means, the AP algorithm exhibits greater robustness to outliers and anomalies, resulting in more stable clustering outcomes. The detailed procedure of the AP clustering algorithm can be found in Ref. [29] and is therefore not reiterated in this paper.

2.3. Risk Measurement of the MIES

Due to the uncertainties associated with energy sources and loads, the system’s operational cost is exposed to risk. CVaR, a widely used risk measure in the field of economics, provides an effective means of capturing tail risk. A brief overview of the CVaR theory is presented below.

Let $x$ denote the decision variable and $y$ the uncertainty variable, with $p\left(y\right)$ representing the probability density function of $y$. Given a fixed decision, the associated loss function is defined as $L(x,y)$. At a specified confidence level $\alpha$, the VaR is defined as follows:

$$f_{VaR}\left(x\right)=min\left\{\theta \left|{\int }_{L(x,y)}^{\infty }p\left(y\right)\ge \alpha \right.\right\}$$

(1)

where $f_{VaR}\left(x\right)$ is the VaR value at the confidence level $\alpha$, and $\theta$ is the boundary value of the loss function $L(x,y)$. In this paper, the CVaR is evaluated at a 95% confidence level ($\alpha =0.95$) to quantify the expected tail risk of operational costs under worst-case scenarios.

The VaR represents the maximum potential loss over a given future period. However, due to its limitations—such as difficulty in capturing tail risk, non-convexity, and failure to satisfy subadditivity—this paper adopts CVaR as the risk measure. CVaR is defined as the expected loss exceeding the VaR threshold, as expressed in Equation (2).

$$F_{CVaR}={\displaystyle \frac{1}{1?\alpha }}{\int }_{L(x,y)\ge \theta }^{\infty }L(x,y)p\left(y\right)dx$$

(2)

Typically, CVaR cannot be directly computed using Equation (2); instead, an auxiliary function is introduced to facilitate its calculation, as shown in Equation (3).

$$F_{CVaR}=f_{VaR}+{\displaystyle \frac{1}{1?\alpha }}{\left[L\left(x,y\right)?\theta \right]}^{+}p(y)dy$$

(3)

where ${\left[L\left(x,y\right)?\theta \right]}^{+}=max\left\{L\left(x,y\right)?\theta ,0\right\}$.

When the probability distribution of the random variable y is continuous, it can be discretized by sampling N scenarios. In this case, Equation (3) can be reformulated as follows:

$$F_{CVaR}=\theta +{\displaystyle \frac{1}{1?\alpha }}{\sum }_{s=1}^N{\left[L\left(x,y\right)?\theta \right]}^{+}{\pi }_s$$

(4)

where $N$ is the number of discrete scenes, and ${\pi }_s$ is the probability of scenes occurring; it reflects the empirical frequency of observed source–load patterns, ensuring statistically representative sampling.

3. Optimization Scheduling Model for MIES

3.1. Objective Function

The objective function of the MIES primarily comprises two components: the operational cost and the risk-related cost measure. The total operational cost includes external energy transaction cost $F_{buy}$, equipment operation and maintenance cost $F_{om}$, system carbon emission cost $F_{co}$, renewable and derivative energy curtailment cost $F_{cut}$, and adjustment expectation cost $F_{ad}$.

$$\left\{\begin{array}{c}\min F=F_1+F_2+\beta F_{CVaR}\\ F_1=F_{buy}+F_{om}+F_{co}+F_{cut}\\ F_2=F_{ad}\end{array}\right.$$

(5)

where $F_{CVaR}$ is the risk measurement cost of the system, and $\beta$ is the risk coefficient. CVaR is mainly used to measure the risk of changes in system operating costs caused by uncertain factors. $F_1$ is the system operation scheduling cost, and $F_2$ is the adjustment expected cost. $t$ represents the scheduling time period, and T is the scheduling cycle $t\in T$.

(1) External energy transaction cost $F_{buy}$

$$F_{buy}={\sum }_{t=1}^T\left(P_{grid,t}{c}_{grid}+V_{gas,t}{c}_{gas}\right)$$

(6)

where $P_{grid,t}$ and $c_{grid}$ are the purchased electricity from the grid and the corresponding electricity prices, respectively. $V_{gas,t}$ and $c_{gas}$ are the natural gas volume purchased from the gas market and the corresponding gas price, respectively.

(2) Equipment operation and maintenance cost $F_{om}$

$$F_{om}={\sum }_{t=1}^T{P}_{k,t}{C}_k$$

(7)

where $P_{k,t}$ is the power consumed by the $k$-th class device, and $C_k$ is the operation and maintenance cost of the $k$-th class device.

(3) Carbon emission cost $F_{co}$

$$F_{co}=\left[{\sum }_{t=1}^T\left(P_{k?1,t}{\epsilon }_{k?1}?M_{co,t}\right)?N_{co}\right]{\lambda }_{co}$$

(8)

where ${\epsilon }_{k?1}$ is the carbon emission coefficient per unit power of the $k?1$st equipment in the mining area’s electric thermal energy system, excluding P2H. $M_{co,t}$ is the amount of carbon dioxide consumed during the hydrogen methanation process, while $N_{co}$ is the system’s carbon quota. ${\lambda }_{co}$ is the average carbon trading price.

(4) Renewable and derivative energy curtailment cost $F_{cut}$

$$F_{cut}={\sum }_{t=1}^T\left[\left(P_{cut,t}^{WT}+P_{cut,t}^{PV}\right)c_{cut}^R+\left(P_{cut,t}^{RT}+P_{cut,t}^{WS}\right)c_{cut}^O\right]$$

(9)

where $P_{cut}^{WT}$ and $P_{cut}^{PV}$ are the abandoned energy of wind and photovoltaic power, respectively. $c_{cut}^R$ represents the abandoned energy cost of renewable energy. $P_{cut}^{RT}$ and $P_{cut}^{WS}$ represent the abandoned energy of mine ventilation air methane and mine water discharge, respectively, while $c_{cut}^O$ represents the abandoned energy cost of derivative energy.

(5) Adjustment expectation cost $F_{ad}$

$$F_{ad}={\sum }_{s=1}^N\left[{\pi }_s{\sum }_{t=1}^T\left(F_{t,s}^{r,MT}+F_{t,s}^{r,grid}\right)\right]$$

(10)

where $F_{t,s}^{r,MT}$ and $F_{t,s}^{r,grid}$ represent the cost of calling gas turbines and external grid backup resources in scenario s, respectively. The detailed calculation formulas for the scheduling and dispatch costs of various types of reserve resources can be found in Ref. [30] and are therefore omitted in this paper.

3.2. Operating Constraints

3.2.1. Two-Stage P2G System

In the P2G system, the first stage involves the electrolysis process in which water is decomposed into hydrogen and oxygen through the electrolyzer.

$$\left\{\begin{array}{c}{P}_t^{EL?H}={\eta }_{EL}{P}_t^{EL}\\ P_{EL}^{min}\le P_t^{EL}\le P_{EL}^{max}\\ P_{down}^{EL}\le P_t^{EL}?P_{t?1}^{EL}\le P_{up}^{EL}\\ V_t^{H2}=P_t^{EL}{\eta }_{EL}/H_{HV}\end{array}\right.$$

(11)

where $P_{up}^{EL}$ and $P_{down}^{EL}$ are the up and down ramp rates of the electrolyzer, respectively. $P_t^{EL?H}$ is the hydrogen production power of the electrolyzer, and $H_{HV}$ is the high calorific value of hydrogen gas. $P_{EL}^{min}$ and $P_{EL}^{max}$ are the minimum and maximum hydrogen production power of the electrolyzer, respectively. $V_t^{H2}$ is the gas production of the P2G system. ${\eta }_{EL}$ is the efficiency of the electrolyzer. The first line reflects the conversion efficiency of electrical energy to hydrogen energy, the second line ensures that the equipment operates under safe conditions, the third line represents the climbing constraints of the equipment, and the fourth line represents the calculation of hydrogen production.

The efficiency of the electrolyzer varies with its power output, exhibiting an initial increase followed by a decrease. Within the operating range considered in this study, the electrolyzer efficiency is approximately linear and is expressed by the following formula:

$${\eta }_{EL}=?{\displaystyle \frac{0.18P_t^{EL}}{P_{rate}^{EL}}}+0.926$$

(12)

where $P_{rate}^{EL}$ is the rated power input that the electrolytic cell can maintain during continuous operation.

After the electrolysis process is completed, the produced hydrogen can be transported to a reactor where it reacts with carbon dioxide to generate ${CH}_4$ through methanation. Methanation enables the storage of excess renewable energy in the form of ${CH}_4$, thereby facilitating the accommodation of renewable energy output within the system. The reaction process can be represented as follows:

$$\left\{\begin{array}{c}{V}_{MR,t}^{C{H}_4}=\xi {\eta }_{MR}{P}_t^{EL?H}/L_{C{H}_4}\\ V_{MR,t}^{C{H}_4}=V_{MR,t}^{C{O}_2}\\ M_{co,t}={\rho }_{co2}{V}_{MR,t}^{C{O}_2}\end{array}\right.$$

(13)

where $V_{MR,t}^{C{H}_4}$ is the volume of natural gas produced by the methane reactor. The energy conversion efficiency of the methane reactor is represented by ${\eta }_{MR}$. $M_{co,t}$ is the consumption of carbon dioxide. ${\rho }_{co2}$ is the density of carbon dioxide. $V_{MR,t}^{C{O}_2}$ is the hourly volumetric consumption of CO₂ by the MR during methanation. $\xi$ is a practical adjustment factor accounting for non-ideal conversion efficiency in the methanation reaction.

The cost of P2G primarily originates from electricity expenses and raw material costs, with CO₂ being the predominant contributor. Considering the costs associated with CO₂ procurement, transportation, and storage, the overall system cost of CO₂ typically exceeds the base price observed in the carbon trading market. This is because the effective CO₂ cost for P2G (0.502 DKK/kg) includes purification and transportation premiums, exceeding the carbon market price (0.252 DKK/kg) due to technical requirements for methanation [8].

3.2.2. Micro Gas Turbine Unit (MT)

The gas turbine obtains natural gas from the external gas network and the P2G system. While generating electricity through combustion, it simultaneously recovers heat from the high-temperature exhaust gases to meet the thermal load demands of the mining site. The operational constraints of the MT unit are as follows:

$$\left\{\begin{array}{c}?s_{MT,t}^{sd}\le s_{MT,t}?s_{MT,t?1}\le s_{MT,t}^{su}\\ 0\le s_{MT,t}^{su}+s_{MT,t}^{sd}\le 1\\ s_{MT,t}{P}_{MT}^{min}\le P_{MT,t}\le s_{MT,t}{P}_{MT}^{max}\\ ?r_{MT}\le P_{MT,t}?P_{MT,t?1}\le r_{MT}\\ P_{MT,t}={\displaystyle \frac{\left(V_{gas,t}+V_{MR,t}^{CH4}\right){\eta }_{MT}{L}_{NG}}{3600}}\\ Q_{MT,t}=P_{MT,t}{\displaystyle \frac{1?{\eta }_{MT}?{\eta }_{MT,l}}{{\eta }_{MT}}}\\ P_{MT,t}^r\le min\left(P_{MT}^{max}?P_{MT,t},r_{MT}\right)\end{array}\right.$$

(14)

where $Q_{MT,t}$ is the residual heat of the MT during the period $t$, and ${\eta }_{MT,l}$ is the heat loss coefficient of the MT. $s_{MT,t}$ is a 0–1 variable, which means MT is on or off during the period $t$. $P_{MT,t}$ is the output of MT during the period $t$. $P_{MT}^{min}$ and $P_{MT}^{max}$ represent the minimum and maximum output of MT, respectively.$?r_{MT}$ represents the ramping rate of MT. $L_{NG}$ is the low calorific value of natural gas. $P_{MT,t}^r$ represents the reserve power from MT. $s_{MT,t}^{su}$ and $s_{MT,t}^{sd}$ is a 0–1 variable, which means MT starts up or not during the period $t$. ${\eta }_{MT}$ is the power generation efficiency of MT.

3.2.3. Exhausted Air Heat Storage Oxidation Unit (RT)

During coal mining operations, low-concentration methane gas (commonly referred to as “exhaust gas”) is generated. Direct discharge of this gas can lead to significant environmental pollution. The operational constraints of the RT are as follows:

$$\left\{\begin{array}{c}{Q}_{RT,t}=V_{fw,t}{d}_{fw,t}{L}_{NG}{\eta }_{RT}\\ Q_{RT}^{min}\le Q_{RT,t}\le Q_{RT}^{max}\\ ?r_{RT}\le Q_{RT,t}?Q_{RT,t?1}\le r_{RT}\end{array}\right.$$

(15)

where $Q_{RT,t}$ is the heat generation of the RT unit. $Q_{RT}^{min}$ and $Q_{RT}^{max}$ are the minimum and maximum values of heat generation, respectively. $V_{fw,t}$ is the flow rate of the exhaust mixed gas, $d_{fw,t}$ is the concentration of the exhaust mixed gas, and ${\eta }_{RT}$ is the efficiency of the RT unit. $r_{RT}$ is the climbing efficiency of the RT unit.

3.2.4. Water Source Heat Pump Unit (WS) and Air Source Heat Pump Unit (AS)

Water source heat pumps (WS) provide heating services to mining areas by extracting low-temperature thermal energy from mine drainage. Air source heat pump (AS) serves as a backup heat source to supplement heating when the thermal energy of mine drainage is insufficient. The operational constraints of WS and AS are as follows:

$$\left\{\begin{array}{c}{Q}_{WS,t}=P_{WS,t}{\eta }_{WS}\\ Q_{WS}^{min}\le Q_{WS,t}\le Q_{WS}^{max}\\ ?r_{WS}\le Q_{WS,t}?Q_{WS,t?1}\le r_{WS}\end{array}\right.$$

(16)

where $Q_{WS,t}$ is the heat generation of the water source heat pump unit, and $P_{WS,t}$ is the corresponding power consumption. ${\eta }_{WS}$ is the heating energy efficiency ratio.

$$\left\{\begin{array}{c}{Q}_{AS,t}=P_{AS,t}{\eta }_{AS}\\ Q_{AS}^{min}\le Q_{AS,t}\le Q_{AS}^{max}\\ ?r_{AS}\le Q_{AS,t}?Q_{AS,t?1}\le r_{AS}\end{array}\right.$$

(17)

where $Q_{AS,t}$ is the heat generation of the air source heat pump unit, and $P_{AS,t}$ is the corresponding power consumption. ${\eta }_{AS}$ is the heating energy efficiency ratio.

3.2.5. Battery and Thermal Energy Storage Unit (BT, TT)

Batteries and thermal energy storage can provide flexible regulation capabilities for system operation by temporarily storing and releasing energy, with the following operational constraints:

$$\begin{array}{c}{P}_t^{BT}=P_0^{BT}+{\sum }_{t=1}^T{E}_t^{ch}{\eta }_{ch}?t?{\sum }_{t=1}^T{\displaystyle \frac{E_t^{dis}}{{\eta }_{dis}}}?t\\ P_T^{BT}=P_0^{BT}\\ P_{BT,t}^{min}\le P_t^{BT}\le P_{BT,t}^{max}\\ 0\le E_t^{ch}\le s_{ch,t}{E}_t^{ch,max},?0\le E_t^{dis}\le s_{dis,t}{E}_t^{dis,max}\\ 0\le s_{ch,t}+s_{dis,t}\le 1\end{array}$$

(18)

$$\begin{array}{c}{Q}_t^{TT}=Q_0^{TT}+{\sum }_{t=1}^T{H}_t^{ch}{\eta }_{TT,ch}?t?{\sum }_{t=1}^T{\displaystyle \frac{H_t^{dis}}{{\eta }_{TT,dis}}}?t\\ Q_T^{TT}=Q_0^{TT}\\ Q_{TT,t}^{min}\le Q_t^{TT}\le Q_{TT,t}^{max}\\ 0\le H_t^{ch}\le u_{ch,t}{H}_t^{ch,max},?0\le H_t^{dis}\le u_{dis,t}{H}_t^{dis,max}\\ 0\le u_{ch,t}+u_{dis,t}\le 1\end{array}$$

(19)

where $P_t^{BT}$ is the energy stored in BT, and $P_0^{BT}$ is the initial state. $E_t^{ch}$ and $E_t^{dis}$ are the charging and discharging power, respectively. ${\eta }_{ch}$ and ${\eta }_{dis}$ are the charging and discharging efficiency, respectively. $P_{BT,t}^{min}$ and $P_{BT,t}^{max}$ are the minimum and maximum power, respectively. $E_t^{ch,max}$ and $E_t^{dis,max}$ are the maximum charging and discharging power, respectively. $s_{ch,t}$ and $s_{dis,t}$ are the 0–1 variables. In addition, TT have similar constraints.

3.2.6. Renewable Energy and Purchased Energy Constraint

$$\left\{\begin{array}{c}0\le P_{cut,t}^{WT}\le P_{WT,t}^f\\ 0\le P_{cut,t}^{PV}\le P_{PV,t}^f\end{array}\right.$$

(20)

where $P_{WT,t}^f$ is the predicted WT output, and $P_{PV,t}^f$ is the predicted PV output.

$$\left\{\begin{array}{c}0\le P_{grid,t}\le P_{grid}^{max}\\ P_{grid,t}^r\le P_{grid}^{max}?P_{grid,t}\\ 0\le V_{gas,t}\le V_{gas}^{max}\end{array}\right.$$

(21)

where $P_{grid}^{max}$ is the maximum amount of electricity purchased from the external power grid, and $V_{gas}^{max}$ is the maximum amount of gas purchased from the external gas grid.

3.2.7. Power Balance Constraint

$$\left\{\begin{array}{c}{P}_{MT,t}+P_{PV,t}+P_{WT,t}+P_{grid,t}+E_t^{dis}=E_t^{ch}+P_{cut,t}^{WT}+P_{cut,t}^{PV}+P_{AS,t}+P_{WS,t}+P_t^{EL}+P_t^{load}\\ Q_{AS,t}+Q_{WS,t}+Q_{RT,t}+Q_{MT,t}+H_t^{dis}=H_t^{ch}+Q_t^{load}\end{array}\right.$$

(22)

where $P_t^{load}$ and $Q_t^{load}$ are the system electrical load and thermal load, respectively.

3.2.8. CVaR Constraint

For discrete distributions, based on the relevant introduction in Section 2.3, CVaR can be transformed into the following linear programming problem for solution.

$$F_{CVaR}=\underset{\xi ,{\theta }_s}{\mathit{min}}\xi +{\displaystyle \frac{1}{1?\alpha }}{\sum }_{s\in \Omega }{\pi }_s{\theta }_s$$

(23)

s.t.

$${\sum }_{t=1}^T\left(F_{buy}+F_{om}+F_{co}+F_{cut}+F_{ad}\right)?\xi \le {\theta }_s$$

(24)

where $\xi$ is the VaR value, $\mathsf{\Omega }$ is the set composed of all scenes $s$. ${\theta }_s$ represents the portion of the scheduling cost of the electric thermal energy system in the mining area that exceeds $\xi$ in scenario s, ${\theta }_s\ge 0$.

3.3. Model Solving

The decision variables in the model can be categorized into the scheduling plan variables $X$, the scenario-specific reserve deployment variables $Y$, and the auxiliary variables ($\xi$ and ${\theta }_s$) introduced in the CVaR constraint formulation. Since both the variable $Y$ and the objective function $F_{CVaR}$ are scenario-dependent, and the feasible domain of variable $Y$ is constrained by the values of variable $X$, it is necessary to determine variable $X$ prior to solving for variable $Y$. Accordingly, the original problem can be reformulated as a bi-level optimization model, as shown in Equations (25) and (26), with the upper-level optimization problem defined in Equation (25).

$$\left\{\begin{array}{c}minF\left(X,Y\right)=F_1\left(X\right)+{\sum }_{s=1}^N{\pi }_S{F}_{2,s}\left(X,Y_s\right)+\beta F_{CVaR}\\ s.t.?\\ h\left(X\right)=0\\ g\left(X\right)\le 0\end{array}\right.$$

(25)

where $F_1\left(X\right)$ and $F_{2,s}\left(X,Y_s\right)$ are the operating and scheduling costs of the MIES and the adjustment costs under scenario $s$, respectively. $h\left(X\right)$ and $g\left(X\right)$ are the equality constraint and inequality constraint related to variable $X$, respectively.

After the variable $X$ is determined, the lower level solves for the backup adjustment amount $Y_s$ in scenarios.

$$\left\{\begin{array}{c}minG\left(X,Y_s\right)={\sum }_{s=1}^N{\pi }_S{F}_{2,s}\left(X,Y_s\right)+\beta F_{CVaR}\\ s.t.?\\ h_s\left(X,Y_s\right)=0,?s=\mathrm{1,2},\dots ,N\\ g_s(X,Y_s)\le 0,?s=\mathrm{1,2},\dots ,N\end{array}\right.$$

(26)

where $G\left(X,Y_s\right)$ is a lower-level optimization model that can be transformed into $N$ independent subproblems for solution.

This article uses an improved genetic algorithm to iteratively optimize the upper variable $X$. After $X$ is determined, the lower objective function and constraints can be solved using mature commercial software CPLEX.

4. System Simulation

4.1. Parameter Settings

In this study, a mine located in the Inner Mongolia Autonomous Region of China is used as a case study to validate the proposed model. The forecasted wind and solar power outputs, as well as various types of loads within the mine’s electricity–heat energy system, are illustrated in Figure 3 and Figure 4. The parameters of different units within the energy system are provided in

Table 1. System equipment parameters.

Parameters	Value	Parameters	Value
$c_{gas}$	$3.14?\mathrm{D}\mathrm{K}\mathrm{K}/{\mathrm{m}}^3$	${\lambda }_{co}$	$0.252?\mathrm{D}\mathrm{K}\mathrm{K}/\mathrm{k}\mathrm{g}$
$P_{EL}^{min}$	50 kW	$P_{EL}^{max}$	400 kW
${\eta }_{MR}$	1.08	$r_{MT}$	30 kW
$P_{MT}^{max}$	500 kW	$L_{NG}$	$38.97?\mathrm{MJ}/{\mathrm{m}}^3$
${\eta }_{RT}$	0.82	${\eta }_{WS}$	3
${\eta }_{ch}$	0.8	${\eta }_{dis}$	0.8
$P_{BT,t}^{min}$	50 kWh	$E_t^{ch,max}$	50 kW
${\eta }_{TT,ch}$	0.75	${\eta }_{TT,dis}$	0.75
$Q_{TT,t}^{min}$	30 kWh	$H_t^{ch,max}$	30 kW
$P_{grid}^{max}$	1000 kW	${\eta }_{MT}$	0.65
${\rho }_{co2}$	$1.977?\mathrm{kg}/{\mathrm{m}}^3$	$C_{PV}$	0.351 DKK/kWh
$C_{MT}$	0.085 DKK/kWh	$C_{RT}$	0.127 DKK/kWh
$H_{HV}$	$35?\mathrm{MJ}/{\mathrm{m}}^3$	$P_{MT}^{min}$	50 kW
$P_{up}^{EL}/P_{down}^{EL}$	50 kW	$r_{RT}$	30 kW
${\eta }_{AS}$	2.7	$E_t^{dis,max}$	50 kW
$P_{BT,t}^{max}$	300 kWh	$Q_{TT,t}^{max}$	200 kWh
$H_t^{dis,max}$	30 kW	$C_{WT}$	0.462 DKK/kWh
${\eta }_{MT,l}$	0.1	$C_{WS}$	0.032 DKK/kWh

. The electricity market price profile is depicted in Figure 5. The system operation model is programmed using MATLAB R2019 and the YALMIP toolbox, with the CPLEX 12.8.0 solver employed for optimization. The computer used features a Quad-Core Intel Core i5 processor and 16 GB of memory.

4.2. Analysis of Source–Load Bilateral Parameter Clustering Results

In this section, the Silhouette index is employed to evaluate the quality of clustering results. This metric measures the similarity of each data point to other points within the same cluster and its dissimilarity to points in other clusters. The Silhouette value ranges from [?1, 1], where a higher value indicates a more appropriate and well-separated clustering outcome. Assume that a given clustering method partitions the original dataset into $n$ clusters. For any data point $i$, the average distance to all other points within the same cluster is defined as $a\left(i\right)$. A smaller value of $a\left(i\right)$ indicates that data point $i$ is more similar to the other points in its own cluster. The minimum of the average distances between data point $i$ and all points in the other clusters is defined as $b\left(i\right)$. A larger value of $b\left(i\right)$ indicates that data point $i$ is less similar to the points in the other clusters. Therefore, an effective clustering method corresponds to a smaller $a\left(i\right)$ and a larger $b\left(i\right)$. Accordingly, the Silhouette index $S\left(i\right)$ is defined as follows:

$$S\left(i\right)={\displaystyle \frac{b\left(i\right)?a\left(i\right)}{max\left\{a\left(i\right),b\left(i\right)\right\}}}$$

(27)

Table 2. Comparison of profile coefficient between AP and K-means.

Number of Categories	K-Means Clustering Silhouette Coefficient	AP Clustering Silhouette Coefficient
3	0.374	0.598
6	0.275	0.418
10	0.151	0.275

presents a comparison of the Silhouette coefficients between the Affinity Propagation (AP) clustering algorithm and the K-means algorithm under different numbers of clusters. As shown in Table 2, for the same number of clusters, the Silhouette index values obtained by the AP clustering algorithm are consistently higher than those of the K-means algorithm, demonstrating the superior clustering performance of AP. Furthermore, unlike K-means, the AP algorithm does not rely on predefined initial values, leading to more stable clustering results.

In this section, both clustering performance and computational efficiency are comprehensively considered, leading to the selection of six clusters. A comparison between the representative scenarios generated through stochastic optimization and the actual historical output curves is illustrated in Figure 6. The probabilities associated with each representative scenario are presented in

Table 3. Probability distribution of typical scenarios.

Scenario	1	2	3	4	5	6	7	8
Probability Value	0.080	0.075	0.070	0.165	0.080	0.145	0.215	0.170

.

As shown in Figure 6, the generated typical scenarios effectively preserve the variation trends of the source–load power curves while exhibiting fluctuations within a reasonable range of the actual power values. This reflects the inherent stochastic characteristics of the source–load power and indicates that the typical scenarios obtained through clustering are representative.

4.3. Results of Optimized Dispatch of Mining Area Energy System

Figure 7 and Figure 8, respectively, show the scheduling results of electrical and thermal energy in the mining area energy system.

As illustrated in Figure 7, during the time periods from 00:00 to 05:00 and from 20:00 to 24:00, solar irradiance is zero. During these hours, electrical power is primarily supplied by gas turbines and wind power units, with the shortfall compensated by electricity purchased from the external grid. At 08:00, solar irradiance gradually increases, and photovoltaic units begin to supply power. As the output from renewable energy sources becomes surplus, the amount of electricity purchased from the external grid decreases significantly. Simultaneously, the power consumption of the P2G unit increases noticeably, contributing to the absorption of excess wind and solar energy. Specifically, when renewable energy generation exceeds load demand during periods of low electricity prices, P2G systems prioritize the conversion of electrical energy into hydrogen through electrolyzers, followed by methanation to produce synthetic natural gas. As exemplified by the operational data at 13:00, the P2G system consumes 168 kW of power while the surplus renewable generation reaches 201 kW, achieving a conversion efficiency of 83.6%. Consequently, the remaining 23.6% of unconverted renewable energy is either stored in battery systems or strategically curtailed based on system operational requirements.

During MIES system operation, the BT unit primarily performs charging and discharging in response to electricity market price fluctuations, thereby achieving the objective of economic arbitrage. As observed from the charging and discharging profile of the BT unit, charging primarily occurs during periods such as 04:00–08:00 and 23:00–24:00. In conjunction with Figure 4, it can be seen that these time intervals correspond to relatively low electricity market prices and high wind power output. The BT unit discharges during periods with relatively high electricity prices, such as 07:00–09:00. During these hours, the electricity market price is comparatively elevated, enabling the ES unit to achieve the goal of “charging at low prices and discharging at high prices,” thereby reducing the overall system cost.

As shown in Figure 8, the MIES system’s thermal load is primarily supplied by the MT unit and the RT unit. Due to the significant environmental harm caused by exhaust gas (low-concentration methane), while the environmental impact of drainage water is relatively minimal, this study assigns a higher penalty cost to the curtailment of exhaust gas. As a result, the output of the exhaust gas oxidation device is comparatively higher. In addition, due to the thermoelectric coupling characteristics of the MT unit, electricity generation is inevitably accompanied by the production of residual heat. This may lead to a mismatch between thermal supply and demand, resulting in potential thermal energy waste. Therefore, the TT unit ensures the efficient utilization of residual heat resources through thermal charging and discharging. For instance, at 03:00, if the TT unit does not perform thermal charging, it may result in a thermal energy waste of approximately 34 kW.

To analyze the impact of scenario probabilities on the types of reserve resources deployed, the reserve utilization under typical Scenario 3 (probability 0.070) and Scenario 7 (probability 0.215) is compared. The corresponding reserve power dispatch is illustrated in Figure 9.

In Scenario 3, both the external power grid and gas turbines are utilized to meet the system’s reserve requirements. In contrast, Scenario 7 involves limited reliance on the external grid. This distinction arises from the lower occurrence probability of Scenario 3, which results in a smaller expected reserve deployment cost. Consequently, minimizing operational scheduling costs becomes the dominant objective, making the use of external grid resources for reserves more economically favorable in this scenario. Conversely, in Scenario 7, which has a higher probability of occurrence, the high cost of grid-based reserve deployment leads to a greater expected adjustment cost. Under these conditions, relying solely on gas turbines for reserve capacity is more conducive to the system’s economic operation.

4.4. Model Comparison and Analysis

In order to verify the effectiveness of the model in this paper, the following scenarios are designed for comparative analysis:

Scheme 1: Deterministic operation model of the MIES without P2G equipment. In this scheme, the predicted (average) values of renewable energy generation and load are used, ignoring uncertainty.

Scheme 2: Stochastic optimization model of the MIES without P2G equipment, but without considering CVaR.

Scheme 3: Robust optimization model of the MIES with P2G equipment.

Scheme 4: Stochastic optimization model of the MIES with P2G equipment, but without considering CVaR.

Scheme 5: Stochastic optimization model of the MIES with P2G equipment, and considering CVaR, that is, the model in this paper.

(1) Operating costs under different solutions

The operating costs of the MIES under different operating schemes are shown in

Table 4. Operating costs under different schemes.

	Operation and Scheduling Costs (DKK)	Adjusting Expected Costs (DKK)	Total Cost (DKK)
Scheme 1	21,047.5	2451.9	23,499.4
Scheme 2	25,578.4	2955.4	28,533.8
Scheme 3	29,304.9	3369.4	32,674.3
Scheme 4	23,067.5	2676.4	25,743.9
Scheme 5	27,585.1	3178.3	30,763.4

.

In terms of total cost, Scheme 1 yields the lowest operational cost at 23,499.4 DKK, while Scheme 3 incurs the highest cost at 32,674.3 DKK. The cost levels of the other schemes fall between those of Scheme 1 and Scheme 3, and all employ stochastic optimization methods to address uncertainties in the mine energy system. The primary reason is that Scheme 1 relies on the predicted values of uncertainties on both the supply and demand sides, without accounting for their stochastic fluctuations, and performs deterministic optimization scheduling for the mine energy system. Although Scheme 1 achieves a lower operational cost, its disregard for the uncertainties on both the supply and demand sides may compromise system operational stability. In contrast, Scheme 3 adopts a robust optimization approach, which schedules the MIES system based on the worst-case scenarios of supply and demand uncertainties. However, since such extreme cases have a low probability of occurrence, the solution tends to be overly conservative, resulting in the highest operational cost.

A comparison between Schemes 2, 4, and 5 reveals that the operational cost of Scheme 5 lies between those of Schemes 2 and 4, with Scheme 4 outperforming Scheme 2 in terms of cost efficiency. The primary reason is that, after integrating the P2G equipment in Scheme 4, the system can better accommodate renewable energy output and reduce natural gas purchases from the external gas network, thereby lowering the overall energy system cost in the MIES. Scheme 5 builds upon Scheme 4 by further incorporating CVaR, which leads to a moderate increase in operational costs. This is because the CVaR-based stochastic optimization method not only considers the expected costs under typical scenarios but also accounts for operational risks caused by supply and demand uncertainties, thus achieving a better balance between economic efficiency and system robustness.

To quantify the operational risks of the system in different scenarios, two evaluation indicators are introduced, and the results are shown in

Table 5. Risk indicators for all schemes.

	EUD (kW)	LOLP (%)
Scheme 1	14.2	12.3
Scheme 3	0	0
Scheme 5	2.7	1.8

.

$$EUD={\sum }_{s=1}^N{\pi }_s\left({\displaystyle \frac{1}{24}}{\sum }_{t=1}^{24}max\left(P_{t,s}^{load}?P_{t,s}^{supply},0\right)\right)$$

(28)

where $EUD$ is the expected unmet demand. That is, the average hourly power shortfall across all scenarios. $P_{t,s}^{supply}$ includes MT, grid, and renewable generation.

$$LOLP={\displaystyle \frac{Number?of?scenarios?with?P_{t,s}^{load}>P_{t,s}^{supply}}{Total?scenarios}}\times 100\%$$

(29)

where $LOLP$ is the percentage of scenarios with any unmet demand.

Deterministic scheduling (Scheme 1) ignores uncertainties, leading to recurrent load shortfalls (EUD > 14 kW in 12.3% of scenarios). While robust optimization (Scheme 3) eliminates these risks, our stochastic-CVaR method (Scheme 5) achieves a 3.4-fold reduction in EUD versus Scheme 1 at only 31% higher cost, demonstrating superior cost–risk balance.

The optimization results illustrating the variation of the system’s total operational cost with respect to the number of generated source–load scenarios under different schemes are presented in Figure 10.

As shown in Figure 10, an increase in the number of scenarios leads to a rise in the total operational cost obtained via robust optimization, which is attributed to the consideration of higher-cost scenarios in scheduling decisions. In contrast, the operational costs obtained under Schemes 4 and 5 remain relatively stable despite changes in scenario quantity. This stability arises because the objective function in stochastic optimization is based on the expected value; therefore, when source–load scenarios follow the same distribution, the number of scenarios has a limited impact on the optimization outcome.

(2) Analysis of the impact of the risk preference coefficient on results

Different risk preference coefficients affect the system’s total operational cost and CVaR, with the results presented in Figure 11.

It can be seen from Figure 11 that with the increase in risk coefficient $\beta$, the total operating cost of the mining area energy system gradually increases, while the value of CVaR gradually decreases, indicating that decision-makers’ aversion to risk deepens and the scheduling strategy tends to be conservative. During the increase in $\beta$, the total operating cost increases from 19,764.1 DKK to 31,543.2 DKK, and CVaR decreases from 17,532.3 DKK to 13,763.4 DKK; that is, a 27.4% decrease in CVaR leads to a 37.3% increase in total operating costs.

When the risk coefficient $\beta$ is low, the operation strategy focuses on minimizing expected costs while neglecting tail risk—e.g., reducing MT reserve capacity and purchasing electricity during low-price periods—resulting in high CVaR and vulnerability to extreme events. At a medium $\beta$, the strategy balances cost and risk by maintaining moderate MT reserves and charging batteries during low-risk periods, thereby lowering CVaR at the expense of higher costs. When $\beta$ is high, the strategy prioritizes worst-case avoidance, such as maintaining high MT reserves, which further reduces CVaR but significantly increases operational costs.

(3) Analysis of the impact of source–load correlation

In order to explore the impact of source–load correlation on the operation of the mining area energy system, the following operation schemes are set.

Scheme 6: Consider source–load correlation (that is, consider the correlation between source loads in Section 2.2).

Scheme 7: Do not consider source–load correlation.

The operation results of the mining area energy system considering the source–load correlation are shown in Section 4.2. This section further analyzes the operation results of the mining area energy system without considering the source–load correlation, as shown in Figure 12.

A comparison of Figure 8 and Figure 12 reveals that, for most time periods, considering the correlation between supply and load leads to reductions in both gas turbine output and electricity procurement from the grid. However, during the midday peak period from 12:00 to 15:00, due to an upward fluctuation in electrical load and a simultaneous downward fluctuation in combined wind and solar output, the gas turbine output increases to meet the demand. Overall, the complementary fluctuations between wind, solar, and load effectively reduce the system’s scheduling costs. The impact of source–load correlation on system operating costs is shown in

Table 6. Impact of source–load correlation on operating costs.

	Operation and Scheduling Costs (DKK)	Adjusting Expected Costs (DKK)	Total Cost (DKK)
Scheme 6	27,585.1	3178.3	30,763.4
Scheme 7	27,999.2	3224.4	31,223.5
Cost reduction rate	1.5%	1.4%	1.4%

.

As shown in Table 6, when correlation is taken into account, both the operational scheduling cost and the expected adjustment cost of the mine energy system decrease by 1.5% and 1.4%, respectively. In terms of total operational cost, considering correlation results in a 1.4% reduction in system cost. Therefore, incorporating source–load correlation in the mine energy system is beneficial for optimizing system operating costs.

5. Conclusions

This paper comprehensively considers the correlations on both the supply and demand sides of the mine energy system, as well as the characteristics of renewable energy surplus in mining areas. A stochastic optimization model for a mine electricity–heat energy system incorporating P2G technology is proposed. CVaR is introduced to quantify the risk posed by multiple uncertainties to the system’s operation. This approach enhances the economic performance of system operation while constraining operational risks within acceptable limits. The main conclusions drawn from the case study analysis are summarized as follows.

(1)

The Latin hypercube sampling (LHS) method based on the correlation coefficient matrix can more accurately characterize the uncertainty and interdependencies of multiple random variables in the mine energy system. The Affinity Propagation (AP) clustering algorithm adaptively determines the number of cluster centers and selects actual scenarios from the original dataset as representative scenarios, resulting in clustering outcomes with strong representativeness.

(2)

By considering the correlation between supply and load and optimizing the operational power of adjustment devices based on the complementary characteristics of multiple uncertainties, the total operational cost of the mine energy system is reduced by 1.4%, thereby enhancing the system’s economic performance. The integration of the P2G unit into the traditional mine energy system effectively addresses issues such as renewable energy surplus and high carbon emissions.

(3)

The economic dispatch model of the mine energy system incorporating CVaR enables decision-makers to formulate reasonable operating strategies according to their risk preferences, achieving a balance between economic efficiency and operational robustness.