Wind-light-water tile integrated resource optimization configuration method with adjustable gas power generation

Abstract

The invention relates to a wind-light-water tile integrated resource optimization configuration method with adjustable gas power generation. The method specifically comprises the following steps: 1, constructing a coal mine-power grid system topological structure; 2, according to the coal mine-power grid system topology established in the step 1, a refined gas extraction model is established based on the actual production process of the coal mine; 3, on the basis of the refined gas extraction model, a gas virtual energy storage model and a gas power generation model containing adjustable gas are constructed in combination with the adjustable characteristic of gas storage; 4, on the basis of the gas virtual energy storage model and the gas power generation model containing the adjustable gas, constructing a wind-light-water-tile integrated resource optimization configuration model by taking the minimum annual total economic cost as a target, and constructing a power distribution network tide constraint, a coal mine safety constraint and an equipment operation constraint; and 5, solving the wind-light-water tile integrated resource optimization configuration model in the step 4 to obtain an optimal wind-light-water tile resource configuration scheme. The method provides powerful support for fusion development of new energy and traditional energy.

Description

Technical Field

[0001] The present invention belongs to the technical field of energy planning, and in particular relates to a method for optimizing the configuration of wind, solar, water and gas integrated resources including adjustable gas power generation. Background Art

[0002] As an important part of my country's energy structure adjustment, gas power generation plays a key role in reducing greenhouse gas emissions and promoting the development of clean energy. With the increase in the proportion of renewable energy, gas power generation has become an important means to balance the load of the power grid and ensure the stability of power supply with its flexible regulation ability. At the same time, gas power generation also provides a way to effectively utilize coal mine gas and avoids the pollution caused by direct gas emissions to the environment. However, how to give full play to the regulatory role of gas power generation in the context of the rapid development of new energy has become a hot issue in the energy field.

[0003] The traditional configuration methods of gas power generation and renewable energy sources such as wind, light, and water mainly rely on their own independent optimization strategies, and fail to fully tap the synergy potential between different energy systems. In these traditional methods, the scheduling and planning of energy systems are often based on historical experience data and static statistical models, ignoring the dynamic characteristics of different energy systems in time and space dimensions. This method shows great limitations when facing energy demand fluctuations and resource supply uncertainties. Especially in terms of gas resource supply, since gas production is affected by multiple factors such as coal mining progress and geological conditions, there is significant uncertainty and volatility, and traditional configuration methods often cannot respond to these changes in a timely and accurate manner. This leads to the fact that when resource supply is tight or fluctuating greatly, it is difficult for the gas power generation system to play its due regulatory role, which in turn affects the economy, reliability and environmental benefits of the overall energy system.

[0004] At present, traditional methods still have the problem of insufficient attention to the synergy between multiple energy systems. For example, renewable energy sources such as wind, light, and water are intermittent and volatile, while gas power generation has strong regulation capabilities. In theory, the balance and optimization of the energy system can be achieved through reasonable configuration. However, most existing methods treat various types of energy as independent units for optimization, and fail to fully utilize the regulation capabilities of gas power generation to compensate for the instability of other renewable energy sources, resulting in waste of energy resources and inefficient system operation. Overall, traditional configuration methods have significant deficiencies in the comprehensive scheduling and optimization of multi-energy systems. They do not fully utilize the adjustability of gas power generation and are difficult to meet the current needs of energy transformation and sustainable development. Summary of the invention

[0005] Purpose of the invention: In order to solve the problems existing in the above-mentioned prior art, the present invention provides a method for optimizing the configuration of wind, solar, water and gas integrated resources including adjustable gas power generation.

[0006] Technical solution: The present invention discloses a method for optimizing the configuration of wind, solar, water and gas integrated resources with adjustable gas power generation, which specifically includes the following steps:

[0007] Step 1: Based on the coal mine power distribution system and the urban power distribution system, construct the coal mine-grid system topology;

[0008] Step 2: Based on the coal mine-grid system topology established in step 1, a refined gas extraction model is constructed based on the actual production process of the coal mine.

[0009] Step 3: Based on a refined gas extraction model Combined with the adjustable characteristics of gas storage, a gas virtual energy storage model and a gas power generation model with adjustable gas are constructed;

[0010] Step 4: Based on the gas virtual energy storage model and the gas power generation model with adjustable gas, a wind, solar, water and gas integrated resource optimization configuration model is constructed with the goal of minimizing the annual total economic cost, and distribution network flow constraints, coal mine safety constraints and equipment operation constraints are constructed;

[0011] Step 5: Solve the wind, solar, water and tile integrated resource optimization configuration model in step 4 to obtain the optimal wind, solar, water and tile resource configuration plan.

[0012] Furthermore, the coal mine-grid system topology structure includes: a coal mine distribution network, an urban distribution network, a gas extraction system and a gas power generation system; the power generation equipment in the coal mine distribution network and the urban distribution network is used as a power supply node, the voltage conversion equipment is used as a substation node, and the load is used as a load node; each node is connected by a branch line; the logical relationship between the coal mine distribution network, the urban distribution network, the gas extraction system and the gas power generation system is that the gas extraction system extracts coal mine gas from the gas extraction station node of the coal mine distribution network, and the extracted gas is converted into electrical energy through the gas power generation system and transmitted to the gas power generation node in the urban distribution network.

[0013] Furthermore, the expression of the refined gas extraction model is as follows:

[0014]

[0015] Where t represents the time period, is the total amount of gas extraction, is the gas extraction volume during tunnel excavation, is the gas extraction volume during mining, and The expression is as follows:

[0016]

[0017]

[0018] Among them, b t is the tunnel excavation task variable at time t. If there is an excavation task, then b t =1, otherwise b t =0; E0 is the gas release per unit time during tunnel excavation, are the gas extraction volumes of the mined coal seam, adjacent layer, and goaf at time t, respectively; a1, a2, a3, b1, and b3 are gas extraction correlation coefficients, is the coal mining volume of the coal mining machine in the I time period.

[0019] Furthermore, the expression of the gas virtual energy storage model is as follows:

[0020]

[0021] Where t represents the time period, are the inlet mass flow rate and outlet mass flow rate of the gas storage tank, C d is the flow coefficient, are the volume fraction of gas before injection into the gas tank, the volume fraction of gas in the gas tank, and the volume fraction of gas released from the gas tank, respectively. in is the volume fraction of gas before pressurization in the booster pump room at time t, is the volume of gas pressurized by the booster pump room at time t; They are respectively the gas pressure before injection into the gas tank, the gas pressure in the gas tank, and the gas pressure released from the gas tank; and M is the gas pressure before being pressurized by the booster pump room at time t, meth is the molar mass of gas, R meth is the ideal gas constant; T meth is the absolute temperature of the gas; is the mass of gas in the gas tank, Δt is the time difference between the gas inlet and gas outlet of the gas tank; It is the gas charge state of the gas storage tank; The maximum mass in the gas tank, is the volume of gas released from the gas storage tank at time t;

[0022] The gas power generation model with adjustable gas is expressed as follows:

[0023]

[0024] in, They are respectively the gas pressure after pressure regulation, the air pressure involved in the mixture, the gas pressure after the mixture, and the gas pressure injected into the gas generator set; The volume of gas after pressure regulation, the volume of gas after mixing, and the volume of air involved in mixing; are the gas volume fraction after pressure regulation at time t, the gas volume fraction after mixing, and the gas volume fraction injected into the gas generator set, T out , T in They are respectively the gas temperature of the mixed gas after being output through the intercooler and the gas temperature before being output; is the power generation capacity of the gas generator set; is the power generation efficiency of the gas generator set; C meth is the calorific value of gas, It is the volume of gas injected into the gas generator set.

[0025] Furthermore, the expression of the wind, solar, water and tile integrated resource optimization configuration model is:

[0026] minC tota =C inve +C oper +C ab n +C poll

[0027] Among them, C tota , C inve , C oper , C aban , C poll They represent the total annual economic cost of the system, investment cost, operating cost, penalty cost for abandoned wind and solar power, and environmental governance cost respectively; C tota , C inve , C oper , C aban , C poll The expression is as follows:

[0028]

[0029] Among them, r is the discount rate; y is the planned life of the equipment, N equ is the total number of equipment types, C k is the unit capacity investment cost of the kth equipment; N k is the construction capacity of the kth type of equipment, t represents the time period, and T is the total number of time periods in a scheduling cycle; The amount of electricity purchased from the large power grid for the system; is the time-of-use electricity price in time period t, μ k is the unit operation and maintenance cost of the kth type of equipment, is the power of the kth device, k PV , k WT They are the penalty costs for abandoning solar power and abandoning wind power respectively; are the abandoned solar power and abandoned wind power respectively; N pol The total amount of pollutant gases emitted by gas units; is the nth pollutant gas emission coefficient, The environmental penalty cost for the nth pollutant gas.

[0030] Furthermore, the coal mine safety constraints include the following constraints:

[0031]

[0032]

[0033] in, They are photovoltaic, wind turbine, and hydroelectric power respectively; They are photovoltaic, wind turbine, hydroelectric and gas power generation at time t-1 respectively; They are the maximum power generation of photovoltaic, wind turbine, hydroelectric and gas power generation respectively; They are the lower limits of ramp power for photovoltaic, wind turbine, hydroelectric and gas power generation respectively; They are the upper limits of the climbing power for photovoltaic, wind turbine, hydropower and gas power generation respectively.

[0034] Furthermore, the power grid flow constraints include the following constraints:

[0035]

[0036] U min ≤U t,i ≤U max

[0037]

[0038] Among them, i, j, and h are all nodes in the topological structure of the coal mine-power grid system, and node i is the upstream node of node j, and node j is the upstream node of node h. bc is the boundary node set, t represents the time period, P t,ij , Q t,j are respectively the active power and reactive power from node i to node j, r ij is the resistance between node i and node j, α t,i , β t,ij , α t,j are slack variables, I t,ij is the current from node i to node j at time t, is the voltage amplitude at node i, is the voltage amplitude of node j; P t,j , Q t,j are the injected active power and reactive power of node j respectively; ij is the reactance between nodes i and j, P t,jh , Q t,jh are respectively the active power and reactive power from node j to node h; U min , U max are the lower and upper limits of the node voltage amplitude, are the maximum and minimum capacities of the reactive power compensation device at node i in the coal mine distribution network, is the reactive compensation amount at node i of the coal mine distribution network at time t.

[0039] Furthermore, the coal mine safety constraints are:

[0040]

[0041] Where t represents the time period, κ t is the operating status of the coal mining machine, 0 indicates maintenance, and 1 indicates operation; t1 and t2 are the start and end times of the maintenance period, respectively, and T is the total number of time periods in a scheduling cycle.

[0042] Beneficial effects: The present invention aims to solve the problems of insufficient utilization of gas resources and insufficient coordinated optimization of multi-energy systems in traditional methods. By establishing a gas virtual energy storage model, a gas power generation model, and an integrated resource optimization configuration model, the method of this article can give full play to the power regulation capability of gas power generation, and combine the characteristics of various renewable energy sources such as wind, solar, and water to achieve the optimal configuration of wind, solar, and water gas units while ensuring the safety and economy of the power grid and coal mines, thereby improving resource utilization efficiency, reducing operating costs, and providing strong support for the integrated development of new energy and traditional energy. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] Figure 1 The present invention is a flow chart of the method.

[0044] Figure 2 This is a topological diagram of the coal mine-grid structure containing adjustable gas power generation according to the present invention.

[0045] Figure 3 The figure is a flow chart of the gas state change of the present invention.

[0046] Figure 4 The photovoltaic, wind turbine and hydroelectric power sources of the present invention are predicted to output power diagrams in the present invention.

[0047] Figure 5The relaxation error diagrams of each branch of the distribution network of the present invention at each time period, wherein Figure (a) is the relaxation error diagram of the power flow of the urban distribution network, and Figure (b) is the relaxation error diagram of the power flow of the coal mine distribution network. DETAILED DESCRIPTION

[0048] The accompanying drawings, which constitute a part of the present invention, are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations on the present invention.

[0049] like Figure 1 The present invention provides a method for optimizing the configuration of wind, solar, water and gas integrated resources with adjustable gas power generation, comprising the following steps:

[0050] (1) Based on the selected coal mine power distribution system and urban power distribution system, the coal mine-grid system topology is constructed, and the system node and branch information, photovoltaic power generation prediction value, wind power generation prediction value, and hydropower generation prediction value are input.

[0051] First, based on the selected coal mine power distribution system and urban power distribution system, the topological structure of the coal mine-grid system is constructed as follows: Figure 2As shown. This topology can be regarded as a comprehensive power network model, which includes the connection relationship between the distribution network inside the coal mine and the urban power grid. Specifically, this topology consists of four parts: coal mine distribution network, urban distribution network, gas extraction system and gas power generation system. The coal mine distribution network consists of 34 nodes and 33 branches, of which node 9 is the gas extraction station node; the urban distribution system consists of 33 nodes and 32 branches, of which node 5 is connected to photovoltaic power generation, node 21 is connected to wind power generation, node 14 is connected to hydropower generation, and node 4 is connected to gas power generation. The logical relationship between the subsystems is that the gas extraction station at node 9 of the coal mine distribution network extracts coal mine gas through the gas extraction system, and the extracted gas is converted into electrical energy through the gas power generation system and transmitted to node 4 in the urban distribution network. In this embodiment, the 34 nodes of the coal mine distribution network include power generation nodes (power supply nodes), substation nodes, and load nodes. The power generation nodes (power supply nodes) are: This is the source of electric energy in the power system, usually referring to the location of a power plant or other type of power generation equipment. At this node, electric energy is output from the generator and injected into the power grid (such as node 1). Substation node: This is an important node in the power system. The substation is responsible for converting high voltage to lower voltage, or vice versa, to meet the needs of different power transmission stages. The substation is a key location for converting electric energy from a high-voltage transmission network to a medium-voltage and low-voltage distribution network (such as nodes 2, 4, 5, 6, 11, 14, 17, 27, and 30). Load node: This is where electricity is consumed, usually referring to the terminal of the distribution network, connected to user-side loads such as coal mining machines, crushers, and transfer machines. The load node is the final use point of electric energy (such as nodes 17, 24, 25, and 34, etc.).

[0052] After building the topology, detailed system node and branch information needs to be entered. This information usually includes but is not limited to the node's voltage level, connection type, capacity limit, and the impedance, transmission capacity, distance and other parameters of the branch between nodes. These parameters are the basic data for power flow calculation, network optimization, fault analysis and other work. Specifically, node information includes the geographical location of each node, electrical characteristics (such as voltage level, phase angle, etc.), and the type of equipment connected to the node (such as generators, transformers, loads, etc.). Branch information includes the length of the power line connecting the two nodes, electrical parameters (such as resistance, reactance, etc.), and maximum transmission capacity. By entering these detailed data, the accuracy and reliability of the model can be ensured, providing a solid foundation for subsequent system analysis.

[0053] In the coal mine-grid system, photovoltaic power generation, wind power generation and hydropower generation are all important renewable energy inputs. In order to achieve efficient scheduling and optimized operation of the system, it is necessary to predict the power generation of these renewable energy sources in advance. The predicted value of photovoltaic power generation is usually based on historical meteorological data, solar radiation intensity, panel inclination and other factors, and the power generation in a specific time period in the future is calculated through a model; the predicted value of wind power generation is based on parameters such as wind speed, wind direction, and wind turbine power curve; the predicted value of hydropower generation usually depends on information such as reservoir water level, flow forecast and generator unit operating characteristics. After these predicted values ​​are input into the system, the power supply in different time periods can be reasonably scheduled in combination with the actual load demand of the system to ensure the balanced operation of the system and maximize the utilization efficiency of renewable energy.

[0054] (2) Based on the coal mine-grid system topology and parameters provided in step 1, a refined gas extraction model is constructed considering the actual production process of the coal mine. Combined with the adjustable characteristics of gas storage, a gas virtual energy storage model and a gas power generation model containing adjustable gas are constructed.

[0055] The gas extraction subsystem mainly consists of two parts: gas extraction and gas storage. Gas extraction is the process of extracting gas from underground coal seams, upper and lower adjacent layers, goafs and other areas to the ground through a pipeline network. It is mainly composed of pipeline systems, control valves, water ring vacuum pumps, auxiliary devices, etc. Gas storage is to store the extracted gas in ground gas storage tanks for temporary storage, which is mainly composed of gas storage tanks.

[0056] 2.1: The refined gas extraction model is constructed by considering the actual production process of coal mines, where the actual production process of coal mines can be expressed as:

[0057]

[0058] G goaf =K g (G goaf1 -G mine ) (4)

[0059]

[0060] Among them, G mine , G adja , G goaf1 are the gas extraction volumes of the mined coal seam, adjacent layer, and goaf respectively; G goaf1 is the total amount of gas that can be extracted from the coal mine; K1, K c They are respectively the gas extraction imbalance coefficient, the ratio of the pre-extraction volume of the coal seam during the working face recovery period to the extraction volume during the pre-extraction period; L1, L2, L3, L4 are respectively the width of the working face recovery area, the length of the recovery working face, the width of the recovery working face and the annual advancement length of the recovery working face; m J, η J are the thickness of the Jth adjacent coal layer and the gas emission rate; K y is the gas extraction rate of the adjacent layer; M is the pre-extraction coal seam, γ is the apparent density of coal; W0, W c and W1 are respectively the original gas content of the mined coal seam, the residual gas content of the coal seam, and the pre-extraction gas content that meets the standard; W 0J , W cJ are the original gas content and residual gas content of the Jth adjacent coal seam, respectively; M C Coal reserves for mining coal seams; is the coal mining volume of the coal mining machine at the first moment; k1 is the negative pressure extraction coefficient; k2 is the expected extraction rate of coal seam gas; k3 is the underground emission rate of coal seam gas; where formula (5) is the relationship between the coal reserves and coal mining volume of the mined coal seam, Δt is the pre-extraction time, K g is the correlation coefficient of gas extraction in goaf.

[0061] In addition, when there is excavation in the coal mine, additional gas is generated during the excavation process:

[0062]

[0063] Where t represents the time period, is the gas extraction volume during tunnel excavation; b t is a 0 / 1 variable indicating whether there is a tunnel excavation task at time t, b t =1 means there is excavation task, b t =0 means there is no excavation task; E0 is the amount of gas released per unit time during tunnel excavation.

[0064] 2.2: The refined gas extraction model is constructed by considering the actual production process of coal mines. The refined gas extraction model can be obtained by substituting formula (5) into formulas (1) to (4), and can be expressed as:

[0065]

[0066] in, They are the total amount of gas extraction, the amount of gas extraction during tunnel excavation, and the amount of gas extraction during mining. are the gas extraction volumes of the mined coal seam, adjacent layer, and goaf at time t, respectively. a1, a2, a3, b1, and b3 are gas extraction correlation coefficients. f1(t), f2(t), and f3(t) are the gas extraction functions of the mined coal seam, adjacent layer, and goaf, respectively.

[0067] When establishing a virtual energy storage model for gas, the following assumptions or restrictions need to be made: the gas is an ideal gas; the temperature and pressure in the gas storage tank are evenly distributed; the dynamic process of the gas storage tank before reaching a stable state; and the temperature is constant during the compression process.

[0068] The extracted gas enters the gas storage tank after passing through the pressure pump room. The changes in gas volume and pressure before and after are as follows:

[0069]

[0070] In the formula, is the volume of gas pressurized by the booster pump room at time t; They are respectively the gas pressure before and after pressurization in the booster pump room at time t; They are respectively the volume fraction of gas before and after pressurization in the booster pump room at time t.

[0071] 2.3: If Figure 3 As shown, the virtual gas energy storage model and the gas power generation model containing adjustable gas are constructed by combining the adjustable characteristics of gas storage, wherein the virtual gas energy storage model can be expressed as:

[0072]

[0073] Among them, C d is the flow coefficient, are the inlet mass flow rate and outlet mass flow rate of the gas storage tank respectively; They are the volume fraction of gas before injection into the gas storage tank, the volume fraction of gas in the gas storage tank, and the volume fraction of gas released from the gas storage tank; They are respectively the gas pressure before injection into the gas tank, the gas pressure in the gas tank, and the gas pressure released from the gas tank; is the mass of gas in the gas storage tank at time t; M meth is the molar mass of gas; V tank is the volume of the gas tank; R meth is the ideal gas constant; T meth is the absolute temperature of the gas; is the gas charge state of the gas tank at time t, indicating the current remaining gas capacity of the gas tank; It is the maximum value of gas pressure and mass in the gas storage tank; is the gas output volume of the gas storage tank at time t.

[0074] 2.4: If Figure 3As shown, the virtual gas energy storage model and the gas power generation model containing adjustable gas are constructed in combination with the adjustable gas storage characteristics, wherein the gas power generation model containing adjustable gas can be expressed as:

[0075]

[0076] in, They are respectively the gas pressure released from the gas storage tank at time t, the gas pressure after pressure regulation, the air pressure involved in the mixture, the gas pressure after mixing, and the gas pressure injected into the gas generator set; They are respectively the volume of gas released from the gas storage tank at time t, the volume of gas after pressure regulation, the volume of gas after mixing, and the volume of air involved in mixing; are the volume fraction of gas released from the gas storage tank at time t, the volume fraction of gas after pressure regulation, the volume fraction of gas after mixing, and the volume fraction of gas injected into the gas generator set; T out , T in They are respectively the gas temperature of the mixed gas after being output through the intercooler and the gas temperature before being output; is the power generation of the gas generator set at time t; is the power generation efficiency of the gas generator set; C meth is the calorific value of gas. It is the volume of gas injected into the gas generator set.

[0077] (3) Based on the coal mine-grid system topology and parameters provided in step 1, and based on the gas extraction model, gas virtual energy storage model, and gas power generation model with adjustable gas provided in step 2, a wind, solar, water, and gas integrated resource optimization allocation model is constructed with the goal of minimizing the annual total economic cost, taking into account the distribution network flow constraints, equipment operation constraints, and coal mine safety constraints.

[0078] 3.1: The wind, solar, water and tile integrated resource optimization allocation model is constructed with the goal of minimizing the annual total economic cost, which can be expressed as:

[0079] minC tota =C inve +C oper +C ab n +C poll (29)

[0080]

[0081] Among them, C tota , C inve , C oper , C aban , Cpoll They represent the total annual economic cost, investment cost, operating cost, penalty cost for wind and solar power abandonment, and environmental governance cost of the system respectively; r is the discount rate; y is the planned life of the equipment; N equ is the total number of equipment types; C k is the unit capacity investment cost of the kth equipment; N k is the construction capacity of the kth type of equipment; T is the total number of time periods in a scheduling cycle, is the amount of electricity purchased by the system from the large power grid at time t; is the time-of-use electricity price at time t; μ k is the unit operation and maintenance cost of the kth type of equipment; is the power of the kth device at time t; k PV , k WT They are the penalty costs for abandoning solar power and abandoning wind power respectively; are the abandoned solar power and abandoned wind power at time t; N pol The total amount of pollutant gases emitted by gas units; is the nth pollutant gas emission coefficient; The environmental penalty cost for the nth pollutant gas.

[0082] 3.2: Considering the power flow constraints of the distribution network, it can be expressed as:

[0083]

[0084] U min ≤U t,i ≤U max (38)

[0085]

[0086] In the formula, i, j, and h are all nodes in the topological structure of the coal mine-power grid system, and node i is the upstream node of node j, and node j is the upstream node of node h. t,ij , Q t,ij are the active and reactive power from node i to node j at time t; r ij 、x ij are the resistance and reactance between node i and node j respectively; I t,ij is the current from node i to node j at time t; P t,j , Q t,j are the injected active and reactive power of node j at time t respectively; Ω bc is the boundary node set; U t,i , U t,j are the voltage amplitudes of node i and node j at time t respectively; U min , U max are the lower and upper limits of the node voltage amplitude respectively; is the reactive compensation amount at node i of the coal mine distribution network at time t; are the maximum and minimum capacities of the reactive power compensation device at node i of the coal mine distribution network; P t,jh , Q t,jh are the active power and reactive power from node j to node h respectively.

[0087] 3.3: The equipment operation constraints mentioned above can be expressed as:

[0088]

[0089]

[0090] In the formula, They are photovoltaic, wind turbine, hydroelectric and gas power generation at time t respectively; They are photovoltaic, wind turbine, hydroelectric and gas power generation at time t-1 respectively; They are the maximum power generation of photovoltaic, wind turbine, hydroelectric and gas power generation respectively; They are the lower limits of ramp power for photovoltaic, wind turbine, hydroelectric and gas power generation respectively; They are the upper limits of the climbing power for photovoltaic, wind turbine, hydropower and gas power generation respectively.

[0091] 3.4: The coal mine safety constraints can be expressed as:

[0092]

[0093] In the formula, κ t is the operating status of the coal mining machine at time t, 0 represents maintenance and 1 represents operation; t1 and t2 are the start and end times of the maintenance period respectively.

[0094] In addition, for the quadratic term in the power flow constraints of the distribution network, a slack variable α is introduced t,i , β t,ij , α t,j ; Using second-order cone relaxation for convex transformation, it can be expressed as:

[0095]

[0096] Therefore, equations (34) to (37) can be transformed into:

[0097]

[0098]

[0099] (4) Based on the wind, solar, water and tile integrated resource optimization configuration model provided in step 3, call the solver to solve and output the wind, solar, water and tile resource configuration capacity.

[0100] In order to effectively solve this optimization configuration model, it is usually necessary to select a suitable solver. According to the nature of the model (with integer variables), the Gurobi solver is selected to solve it. Finally, the optimal wind, solar, hydropower and gas resource configuration plan is output. These output results mainly include: the configuration capacity of photovoltaic, wind power, hydropower, and gas, and the total system cost.

[0101] For this embodiment, simulation verification is performed. The system configuration and parameters involved are shown in Table 1; the time-of-use electricity price is shown in Table 2; Figure 4 Predict the output of photovoltaic, wind turbines and hydropower sources in the day ahead; set the discount rate for all equipment in the system to 5%; and set the system planning period to 15 years.

[0102] Table 1

[0103] equipment Unit investment cost / (yuan / kWh) Unit maintenance cost (yuan / kW) Investment capacity limit Photovoltaic 3400 0.08 10MW Fan 3900 0.11 10MW Hydropower 8000 0.09 20MW Gas power generation 9000 0.05 20MW Gas tank <![CDATA[5000 yuan / m 3 > <![CDATA[0.01 yuan / m 3 > <![CDATA[10000m 3 ]]>

[0104] Table 2

[0105]

[0106] This embodiment uses the Gurobi solver to solve the model. The final output of the optimal wind, solar, water and tile resource configuration results is shown in Table 3, and the total system cost is shown in Table 4;

[0107] Table 3

[0108] equipment Planning capacity Photovoltaic 500kW Fan 1000kW Hydropower 500kW Gas power generation 1700kW Gas tank <![CDATA[9000m 3 ]]>

[0109] Table 4

[0110]

[0111]

[0112] The hardware system configuration of this embodiment is: Inter Core i7-11370H processor, 64-bit Windows 11 operating system; the software system configuration is: MATLAB_R2021b software, calling the GUROBI solver (version 9.5.0) through the YALMIP toolkit for example simulation.

[0113] From the system configuration results in Table 3, it can be seen that it is necessary to configure 500kW photovoltaic, 1000kW wind turbine, 500kW hydropower, 1700kW gas, and 9000m gas storage tank. 3It can meet the system operation requirements. And from Table 4, we can see that the total annual cost of the system is 27.9599 million yuan, of which the investment cost is 2.4033 million yuan, the operation cost is 25.4228 million yuan, the penalty cost for wind and solar abandonment is 0 yuan, and the environmental governance cost is 133,800 yuan.

[0114] In addition, the error of the method adopted in this embodiment during planning is as follows Figure 5 As shown, from Figure 5 It can be seen from (a) in the figure that the maximum magnitude of the power relaxation error of the urban distribution network branch is 1e-3. Figure 5 It can be seen from (b) that the maximum magnitude of the power relaxation error of the coal mine distribution network branch is 1e-6, which fully demonstrates that the method adopted in this embodiment is feasible for performing second-order cone relaxation without losing accuracy.

[0115] To sum up, the method proposed in the present invention can give full play to the power regulation capability of gas power generation, and combine the characteristics of various renewable energy sources such as wind, solar, and water, to achieve the optimal configuration of wind, solar, and water-gas units while ensuring the safety and economy of the power grid and coal mines, thereby improving resource utilization efficiency, reducing operating costs, and providing strong support for the integrated development of new and traditional energy.

[0116] It should also be noted that the various specific technical features described in the above specific embodiments can be combined in any suitable manner without contradiction. In order to avoid unnecessary repetition, the present invention will not further describe various possible combinations.

Claims

1. A method for optimizing the configuration of wind, solar, water and gas integrated resources with adjustable gas power generation, characterized in that: The specific steps include: Step 1: Based on the coal mine power distribution system and the urban power distribution system, construct the coal mine-grid system topology; Step 2: Based on the coal mine-grid system topology established in step 1, a refined gas extraction model is constructed based on the actual production process of the coal mine. Step 3: Based on a refined gas extraction model Combined with the adjustable characteristics of gas storage, a gas virtual energy storage model and a gas power generation model with adjustable gas are constructed; Step 4: Based on the gas virtual energy storage model and the gas power generation model with adjustable gas, a wind, solar, water and gas integrated resource optimization configuration model is constructed with the goal of minimizing the annual total economic cost, and distribution network flow constraints, coal mine safety constraints and equipment operation constraints are constructed; Step 5: Solve the wind, solar, water and tile integrated resource optimization configuration model in step 4 to obtain the optimal wind, solar, water and tile resource configuration plan.

2. The method for optimizing the configuration of wind, solar, water and gas integrated resources with adjustable gas power generation according to claim 1 is characterized in that: The coal mine-grid system topology structure includes: a coal mine distribution network, an urban distribution network, a gas extraction system and a gas power generation system; the power generation equipment in the coal mine distribution network and the urban distribution network is used as a power supply node, the voltage conversion equipment is used as a substation node, and the load is used as a load node; each node is connected by a branch line; the logical relationship between the coal mine distribution network, the urban distribution network, the gas extraction system and the gas power generation system is that the gas extraction system extracts coal mine gas from the gas extraction station node of the coal mine distribution network, and the extracted gas is converted into electrical energy through the gas power generation system and transmitted to the gas power generation node in the urban distribution network.

3. The method for optimizing the configuration of wind, solar, water and gas integrated resources with adjustable gas power generation according to claim 1 is characterized in that: The expression of the refined gas extraction model is as follows: Where t represents the time period, is the total amount of gas extraction, is the gas extraction volume during tunnel excavation, is the gas extraction volume during mining, and The expression is as follows: Among them, b t is the tunnel excavation task variable in time period t. If there is an excavation task, then b t =1, otherwise b t =0; E0 is the gas release per unit time during tunnel excavation, are the gas extraction volumes of the mined coal seam, adjacent layer, and goaf in time period t, respectively; a1, a2, a3, b1, and b3 are gas extraction correlation coefficients, is the coal mining volume of the coal mining machine in the I time period.

4. The method for optimizing the configuration of wind, solar, water and gas integrated resources with adjustable gas power generation according to claim 1 is characterized in that: The expression of the gas virtual energy storage model is as follows: Where t represents the time period, are the inlet mass flow rate and outlet mass flow rate of the gas storage tank, C d is the flow coefficient, are the volume fraction of gas before injection into the gas tank, the volume fraction of gas in the gas tank, and the volume fraction of gas released from the gas tank, respectively. in is the gas volume fraction before pressurization in the booster pump room during time period t, is the volume of gas pressurized by the booster pump room during time period t; They are respectively the gas pressure before injection into the gas tank, the gas pressure in the gas tank, and the gas pressure released from the gas tank; and is the gas pressure before the booster pump room pressurizes the gas in time period t, M meth is the molar mass of gas, R meth is the ideal gas constant; T meth is the absolute temperature of the gas; is the mass of gas in the gas tank, Δt is the time difference between the gas inlet and gas outlet of the gas tank; It is the gas charge state of the gas storage tank; is the maximum mass in the gas tank, is the volume of gas released from the gas storage tank during time period t; The gas power generation model with adjustable gas is expressed as follows: in, They are respectively the gas pressure after pressure regulation, the air pressure involved in the mixture, the gas pressure after the mixture, and the gas pressure injected into the gas generator set; The volume of gas after pressure regulation, the volume of gas after mixing, and the volume of air involved in mixing; are the gas volume fraction after pressure regulation, the gas volume fraction after mixing, and the gas volume fraction injected into the gas generator set at time period t, respectively. out , T in They are respectively the gas temperature of the mixed gas after being output through the intercooler and the gas temperature before being output; is the power generation capacity of the gas generator set; is the power generation efficiency of the gas generator set; C meth is the calorific value of gas, It is the volume of gas injected into the gas generator set.

5. The method for optimizing the configuration of wind, solar, water and gas integrated resources with adjustable gas power generation according to claim 4 is characterized in that: The expression of the wind, solar, water and tile integrated resource optimization configuration model is: minC tota =C inve +C oper +C abn +C poll Among them, C tota , C inve , C oper , C aban , C poll They represent the total annual economic cost of the system, investment cost, operating cost, penalty cost for abandoned wind and solar power, and environmental governance cost respectively; C tota , C inve , C oper , C aban , C poll The expression is as follows: Among them, r is the discount rate; y is the planned life of the equipment, N equ is the total number of equipment types, C k is the unit capacity investment cost of the kth equipment; N k is the construction capacity of the kth type of equipment, T is the total number of time periods in a scheduling cycle; The amount of electricity purchased from the large power grid for the system; is the time-of-use electricity price in time period t, μ k is the unit operation and maintenance cost of the kth type of equipment, is the power of the kth device, k PV , k WT They are the penalty costs for abandoning solar power and abandoning wind power respectively; are the abandoned solar power and abandoned wind power respectively; N pol The total amount of pollutant gases emitted by gas units; is the nth pollutant gas emission coefficient, The environmental penalty cost for the nth pollutant gas.

6. The method for optimizing the configuration of wind, solar, water and gas integrated resources with adjustable gas power generation according to claim 4 is characterized in that: The coal mine safety constraints include the following constraints: in, They are photovoltaic, wind turbine, and hydroelectric power respectively; They are photovoltaic, wind turbine, hydroelectric and gas power generation at time t-1 respectively; They are the maximum power generation of photovoltaic, wind turbine, hydroelectric and gas power generation respectively; They are the lower limits of ramp power for photovoltaic, wind turbine, hydroelectric and gas power generation respectively; They are the upper limits of the climbing power for photovoltaic, wind turbine, hydropower and gas power generation respectively.

7. The method for optimizing the configuration of wind, solar, water and gas integrated resources with adjustable gas power generation according to claim 1 is characterized in that: The power grid flow constraints include the following constraints: IN min ≤U t,i ≤U max Among them, i, j, and h are all nodes in the topological structure of the coal mine-power grid system, and node i is the upstream node of node j, and node j is the upstream node of node h. bc is the boundary node set, t represents the time period, P t,ij , Q t,j are respectively the active power and reactive power from node i to node j, r ij is the resistance between node i and node j, α t,i , β t,ij , α t,j are slack variables, I t,ij is the current from node i to node j during time period t, is the voltage amplitude at node i, is the voltage amplitude of node j; P t,j , Q t,j are the injected active power and reactive power of node j respectively; ij is the reactance between nodes i and j, P t,jh , Q t,jh are respectively the active power and reactive power from node j to node h; U min , U max are the lower and upper limits of the node voltage amplitude, are the maximum and minimum capacities of the reactive power compensation device at node i in the coal mine distribution network, is the reactive compensation amount at node i of the coal mine distribution network at time t.

8. The method for optimizing the configuration of wind, solar, water and gas integrated resources with adjustable gas power generation according to claim 1 is characterized in that: The coal mine safety constraints are: Where t represents the time period, κ t is the operating status of the coal mining machine, 0 indicates maintenance, and 1 indicates operation; t1 and t2 are the start and end times of the maintenance period, respectively, and T is the total number of time periods in a scheduling cycle.

Images