A simulation scheduling method and system for surface water-groundwater-mine water coupling

Abstract

The application discloses a kind of surface water-groundwater-mine water coupling simulation scheduling method and system, comprising: obtaining the catchment around target mine in research area, screening the nearest reference point with target mine;Calculate the water quality fluctuation coefficient of the catchment on the nearest reference point with target mine;Calculate the influence weight of the influence catchment on target mine;Build the water balance equation of each influence catchment, the groundwater control equation associated with target mine water and groundwater, the amount of water infiltration into groundwater in influence catchment as groundwater source and sink item, simulate the change state of mine water, and plan the amount of water that target mine can be scheduled at current simulation scheduling time.Such as water quality sensor, water level sensor, data acquisition module and data processing device.System is suitable for dynamic optimization management of complex water resource system in mining area and surrounding area, especially suitable for water resource safety scheduling and water disaster prevention and control of large water inflow mining area.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of water resources scheduling, in particular to a simulation scheduling method and system for surface water-groundwater-mine water coupling. BACKGROUND

[0002] In the utilization and management of water resources in a mining area, surface water, groundwater and mine water form a complex system that is interrelated and mutually influenced: surface water recharges groundwater through infiltration, groundwater forms mine water inflow through fissure penetration, and mine water can be reused as production water after treatment, and the dynamic coupling of water quantity and quality of the three is the core of sustainable utilization of water resources and water hazard prevention and control in the mining area.

[0003] However, the prior art has the following key defects: the simulation accuracy of mine water is insufficient, specifically, the existing surface water-groundwater coupling model simplifies mine water as a constant groundwater discharge term (such as a fixed mine drainage), without considering the dynamic influence of changes in groundwater and surface water on mine water level, and cannot accurately simulate the chain response of the coupling of the three, resulting in generally large simulation errors of mine water inflow, and the scheduling scheme for mine water utilization is prone to cause water hazards or waste of water resources. SUMMARY

[0004] In view of the above deficiencies of the prior art, the present application provides a simulation scheduling method and system for surface water-groundwater-mine water coupling.

[0005] To achieve the above-mentioned application purposes, the technical solution adopted by the present application is:

[0006] A simulation scheduling method for surface water-groundwater-mine water coupling is provided, comprising the following steps:

[0007] S1: obtaining the catchment area around the target mine in the study area, and taking a plurality of uniform reference points on the catchment area, and selecting the reference point closest to the target mine according to the distance between the catchment area and the target mine and the topographic height between the catchment area and the target mine;

[0008] S2: calculating the aeration zone diffusion flux between the location of the catchment area where the reference point closest to the target mine is located and the target mine, selecting the influence catchment area of the target mine, and calculating the water quality fluctuation coefficient of the catchment area where the reference point closest to the target mine is located on the influence catchment area;

[0009] S3: calculating the influence weight of the influence catchment area on the target mine according to the water quality fluctuation coefficient of the target mine and the water quality fluctuation coefficient of the catchment area where the reference point is located on the influence catchment area;

[0010] S4: constructing a water balance equation for each influence catchment area, a groundwater control equation related to the target mine water and groundwater, and discretizing the groundwater control equation to obtain a coupled equation set.

[0011] The amount of water infiltrating into groundwater from the basin is taken as a groundwater source and sink item, and is input into the coupling equation set to simulate the change state of the mine water under the change of the groundwater level and the basin water storage capacity, and to plan the water amount that can be dispatched by the target mine at the current simulation scheduling time.

[0012] Further, step S1 comprises:

[0013] S11: Obtain the basin around the target mine in the study area, and take a plurality of uniform reference points on the surrounding basin, obtain the coordinates of each reference point in the overhead plane of the study area , and the coordinates of the point where the target mine is located in the overhead plane of the study area , i is the number of the basin in the study area, e is the number of the reference point on the basin;

[0014] S12: Extract the elevation of the riverbed low point at the position of the reference point , and compare it with the elevation of the bottom of the target mine , respectively the elevation of the target mine mouth and the vertical depth of the target mine mouth;

[0015] If , the reference point is deleted from the basin, otherwise, the reference point is retained;

[0016] S13: According to the coordinates of the retained reference points , u the number of the retained reference points, calculate the straight line distance between the reference points and the target mine, and select the reference points closest to the target mine;

[0017] ;

[0018] Wherein, is the set of retained reference points on the basin.

[0019] Further, step S2 comprises:

[0020] S21: Calculate the air gap diffusion flux between the position of the selected reference point closest to the target mine and the target mine ;

[0021] ;

[0022] Wherein, is the air gap connectivity between the basin where the reference point is located and the target mine, Let be the permeability coefficient of the vadose zone. For reference point The water level in the basin For reference point The elevation of the top of the vadose zone between the watershed and the target mine. For the thickness of the vadose zone, For reference point The driving force of water level changes in the basin In time Internal reference point Water level changes in the basin For reference point The average annual water level change in the basin;

[0023] S22: Set the threshold for diffusion flux ;

[0024] like Then determine the reference point. The basin i It is the watershed affected by the target mine;

[0025] like Then determine the reference point. The basin i The watershed not affected by the target mine;

[0026] S23: Obtain all impact basins of the target mine within the study area and construct an impact basin set. Collect data on the impact of the watershed. Each affected watershed a upper reference point The watershed is based on time series. Water quality data , For the current simulated scheduling time, For the current simulated scheduling time Water quality data;

[0027] S24: Based on water quality data Calculate the impact basin a upper reference point The watershed in the time series Water quality fluctuation coefficient ;

[0028] ;

[0029] in, b Time series The time number in Time series At any time in the process.

[0030] Further, step S3 includes:

[0031] S31: Time series data acquisition Water quality data of internal target mines And calculate the water quality fluctuation coefficient of the target mine. ;

[0032] ;

[0033] in, For the current simulated scheduling time Water quality data for the target mine;

[0034] S32: Based on water quality fluctuation coefficient With water quality fluctuation coefficient Calculate the impact basin a Impact weight of the target mine ;

[0035] ;

[0036] in, A To influence the collection of watersheds The number of watersheds affected by the floods.

[0037] Further, step S4 includes:

[0038] S41: Constructing the set of watersheds with influence Each affected watershed a The water balance equation;

[0039] ;

[0040] in, To affect the watershed a Reference point The watershed is located at the current simulated scheduling time. Water storage capacity To affect the watershed a Reference point The watershed at the previous moment Water storage capacity To affect the watershed a Reference point Water level changes in the basin For reference point The area of ​​the watershed, Reference points The inflow and outflow of water in the basin Surface runoff, The amount of water within the watershed that seeps into groundwater;

[0041] S42: Considering the dynamic pooling of mine water inflow, construct the groundwater control equation relating the target mine water and groundwater;

[0042] ;

[0043] in, Groundwater at x , y Permeability coefficient in the z-direction Groundwater level, W For underground water source collection, For the collection of mine water inflow, The specific yield of groundwater;

[0044] S43: Discretize the groundwater control equations, with time variables... t Discretize into an implicit scheme to obtain the discretized groundwater control equations;

[0045] ;

[0046] in, For discrete space terms, For time step, These are the simulated scheduling times. , The groundwater level;

[0047] S44: By simultaneously setting the discretized groundwater control equations and water balance equations, a set of coupled equations is obtained;

[0048] ;

[0049] This will affect the amount of water infiltrating into groundwater within the watershed. As a source of groundwater W Take at this time ; Pre-simulation scheduling time Impact on the watershed a Reference point The amount of water within the watershed that seeps into groundwater;

[0050] The system collects the water storage of the reference point in each affected watershed and the groundwater level around the target mine at the current simulation scheduling moment. It then uses the coupled equation set to calculate the mine water inflow term and simulates the change state of mine water under the conditions of changes in groundwater level and watershed water storage.

[0051] S45: Based on the calculated mine water inflow, plan the amount of water that can be dispatched at the current simulated dispatch time of the target mine.

[0052] The application provides a simulation scheduling system applying the simulation scheduling method of surface water-groundwater-mine water coupling, and the simulation scheduling system comprises the following:

[0053] a water quality sensor for collecting water quality data in a target mine water affected by a basin;

[0054] a water level sensor for collecting water level data in the target mine water affected by the basin; a data collection module for collecting the water quality data and the water level data collected by the water quality sensor and the water level sensor and sending the data to a data processing device;

[0055] The data processing device is provided with a memory and a processor, the memory is used for storing computer program instructions, and the computer program instructions are executed by the processor to realize the simulation scheduling method.

[0056] The simulation scheduling method has the following beneficial effects:

[0057] The simulation scheduling method is used for simulating the change of mine water level of surface water-groundwater-mine water coupling, and the real-time mine water inflow item of mine water is obtained through the change of surface water and groundwater, so that accurate data support and reference are provided for the scheduling of mine water.

[0058] The coupling equation set of surface water-groundwater-mine water correlation is constructed, and the influence weight of surface water on groundwater and mine water is introduced, so that the change state of mine water is accurately simulated, a precise mine water scheduling scheme is formulated, water resource waste and water resource abuse are avoided, and the method is suitable for dynamic optimization management of complex water resource systems in a mining area and surrounding areas, and is especially suitable for water resource safety scheduling and water disaster prevention and control of a large mine water inflow mining area. BRIEF DESCRIPTION OF DRAWINGS

[0059] Figure 1 The simulation scheduling method of surface water-groundwater-mine water coupling is a flowchart. DETAILED DESCRIPTION

[0060] The specific embodiments of the application are described below, so that those skilled in the art can understand the application, but it should be clear that the application is not limited to the scope of the specific embodiments, and for those skilled in the art, it is obvious that various changes are within the spirit and scope of the application defined and determined by the appended claims, and all applications utilizing the concept of the application are included in the protection.

[0061] As shown in Figure 1 A simulation scheduling method of surface water-groundwater-mine water coupling comprises the following steps:

[0062] S1: Obtain the watershed surrounding the target mine within the study area, and select several uniform reference points within the watershed. Based on the distance between the watershed and the target mine, and the topographic elevation between the watershed and the target mine, select the reference point closest to the target mine. Step S1 specifically includes the following steps:

[0063] S11: Obtain the watershed surrounding the target mine within the study area, and take several uniform reference points on the surrounding watershed. Obtain the coordinates of each reference point in the top-view plane of the study area. and the coordinates of the target mine location in the top-down plane of the study area. , i The watershed designation within the study area, e The reference point number in the watershed;

[0064] In this embodiment, the top-view plane of the study area can be any map plane. By establishing a unified two-dimensional coordinate system on the map plane, as long as the coordinate system is consistent, the position coordinates of the reference point and the mine can be used to calculate the relative position coordinates of the two.

[0065] S12: Extract the elevation of the lowest point of the riverbed at the location of the reference point. and the elevation of the bottom of the target mine Comparison, These are the elevation of the target mine shaft and the vertical depth of the target mine shaft, respectively.

[0066] like If the condition is met, the reference point will be removed from the watershed; otherwise, the reference point will be retained.

[0067] S13: Based on the coordinates of the retained reference points , u To retain reference point numbers, calculate the straight-line distance between the reference point and the target mine, and select the reference point closest to the target mine;

[0068] ;

[0069] in, A set of reference points preserved within the watershed.

[0070] S2: Calculate the vadose zone diffusion flux between the watershed containing the nearest reference point to the target mine and the target mine; screen the watersheds affected by the target mine; and calculate the water quality fluctuation coefficient of the watershed containing the nearest reference point to the target mine within the affected watershed. Step S2 specifically includes the following steps:

[0071] S21: Calculate the reference point closest to the target mine selected from the screening. Vadose zone diffusion flux between the location of the watershed and the target mine ;

[0072] ;

[0073] in, For reference point The connectivity of the vadose zone between the watershed and the target mine is shown at a reference point in this embodiment. connectivity connectivity It is a dimensionless parameter; The permeability coefficient of the vadose zone is expressed in m / d (meters per day). It varies with the water content of the vadose zone; the higher the water content, the higher the permeability coefficient. The larger the value, the more it can be measured using a soil moisture sensor; For reference point Water level (m) in the basin; For reference point The elevation (m) of the top of the vadose zone between the watershed and the target mine was obtained from the borehole data of the study area; For the thickness of the vadose zone, For reference point The water level change driving coefficient (dimensionless parameter) of the watershed. In time Internal reference point Water level changes in the basin For reference point The average annual water level change in the basin;

[0074] S22: Set the threshold for diffusion flux ;

[0075] like Then determine the reference point. River Basin i It is the watershed affected by the target mine;

[0076] like Then determine the reference point. River Basin i The watershed not affected by the target mine;

[0077] S23: Obtain all impact basins of the target mine within the study area and construct an impact basin set. Collect data on the impact of the watershed. Each affected watershed a upper reference point The watershed is based on time series. Water quality data , This is the current simulated scheduling time. water quality data of the current simulation scheduling time ;

[0078] The length of each influence basin is too long, and the water quality data of different basin sections will have large differences. The application only collects the water quality data of the basin section closest to the target mine, and the water quality data of this basin section has the greatest impact on the target mine.

[0079] S24: According to the water quality data of the influence basin a , the water quality fluctuation coefficient of the reference point on the time series of the basin where the reference point is located is calculated. ;

[0080] ;

[0081] Wherein, b is the time number in the time series , and is any time in the time series .

[0082] S3: According to the water quality fluctuation coefficient of the target mine and the water quality fluctuation coefficient of the basin where the reference point of the influence basin is located, the influence weight of the influence basin on the target mine is calculated. Step S3 specifically includes the following steps:

[0083] S31: Collect the water quality data of the target mine in the time series , and calculate the water quality fluctuation coefficient of the target mine ; ;

[0084] ;

[0085] Wherein, is the water quality data of the target mine at the current simulation scheduling time ;

[0086] S32: According to the water quality fluctuation coefficient and the water quality fluctuation coefficient , the influence weight of the influence basin a on the target mine is calculated ;

[0087] ;

[0088] Wherein, A is the number of influence basins in the influence basin set .

[0089] S4: Construct the water balance equation for each affected watershed and the groundwater control equation relating the target mine water and groundwater. Discretize the groundwater control equation and solve it simultaneously to obtain a set of coupled equations.

[0090] The amount of water infiltrating into groundwater within the basin is used as the groundwater source term and input into the coupled equation set to simulate the changes in mine water under the conditions of changes in groundwater level and watershed water storage, and to plan the amount of water that can be dispatched at the current simulated dispatching time of the target mine.

[0091] Step S4 specifically includes the following steps:

[0092] S41: Constructing the set of watersheds with influence Each affected watershed a The water balance equation;

[0093] ;

[0094] in, To affect the watershed a Reference point The watershed is located at the current simulated scheduling time. Water storage capacity To affect the watershed a Reference point The watershed at the previous moment Water storage capacity To affect the watershed a Reference point Water level changes in the basin For reference point The area of ​​the watershed, Reference points The inflow and outflow of water in the basin Surface runoff, The amount of water within the watershed that seeps into groundwater;

[0095] S42: Considering the dynamic pooling of mine water inflow, construct the groundwater control equation relating the target mine water and groundwater;

[0096] ;

[0097] in, Groundwater at x , y Permeability coefficient in the z-direction Groundwater level, W It represents the inflow of groundwater, which is positive during recharge and negative during discharge (e.g., precipitation infiltration is recharge, artificial extraction is discharge). The negative value represents the mine water inflow, the positive value represents the mine drainage, and the positive value represents the mine water inflow. The specific yield of groundwater (dimensionless).

[0098] S43: Discretize the groundwater control equations and decompose the time variables. t Discretize into an implicit scheme to obtain the discretized groundwater control equations;

[0099] ;

[0100] in, For discrete space terms, The time step is typically one day. These are the simulated scheduling times. , The groundwater level;

[0101] S44: By simultaneously setting the discretized groundwater control equations and water balance equations, a set of coupled equations is obtained;

[0102] ;

[0103] This will affect the amount of water infiltrating into groundwater within the watershed. As a source of groundwater W Take at this time ; Pre-simulation scheduling time Impact on the watershed a Reference point The amount of water within the watershed that seeps into groundwater;

[0104] The system collects the water storage of the reference point in each affected watershed and the groundwater level around the target mine at the current simulation scheduling moment. It then uses the coupled equation set to calculate the mine water inflow term and simulates the change state of mine water under the conditions of changes in groundwater level and watershed water storage.

[0105] This method equates the diffusion effect of surface water from a watershed to the target mine to the diffusion effect of surface water to groundwater. Since the water quality data of the target mine is relatively stable, accurate influence weights can be calculated. As a reference standard for the diffusion of surface water into the ground, when calculating the amount of surface water infiltrating into groundwater, the influence weight of the watershed on the target mine can be directly used, because the influence of the watershed on the target mine is approximately equivalent to the influence of the watershed on groundwater, thus obtaining the total amount of groundwater infiltrating into each influencing watershed.

[0106] S45: Plan the water amount that can be dispatched at the current simulation scheduling moment of the target mine according to the calculated mine water inflow item. When the calculated mine water inflow item is negative, it means that the target mine can dispatch water outward (i.e. mine drainage), and when the mine water inflow item is positive, it means that the target mine cannot dispatch water outward.

[0107] A simulation scheduling system applying the above-mentioned simulation scheduling method of surface water-groundwater-mine water coupling, comprising:

[0108] A water quality sensor for collecting water quality data in the basin and the target mine;

[0109] A water level sensor for collecting water level data in the basin and the target mine;

[0110] A data collection module for collecting water quality data and water level data collected by the water quality sensor and the water level sensor, and sending them to the data processing device;

[0111] The data processing device is equipped with a memory and a processor, the memory is used to store computer program instructions, and the computer program instructions are executed by the processor to realize the above-mentioned simulation scheduling method.

[0112] The scheme is used to simulate the change of mine water level of surface water-groundwater-mine water coupling, and the real-time mine water inflow item of mine water is obtained through the change of surface water and groundwater, which provides accurate data support and reference for the scheduling of mine water.

[0113] By constructing the coupling equation set of surface water-groundwater-mine water correlation and introducing the influence weight of surface water on groundwater and mine water, the change state of mine water is accurately simulated, the accurate mine water scheduling scheme is facilitated to be made, and the waste and abuse of water resources are avoided. It is suitable for dynamic optimization management of complex water resource system in mining area and surrounding area, especially for water resource safety scheduling and water disaster prevention and control in large water inflow mining area.

Claims

1. A simulation scheduling method for surface water-groundwater-mine water coupling, characterized in that, The method comprises the following steps: S1: obtaining the catchment area around the target mine in the study area, and taking a plurality of uniform reference points on the catchment area, and screening the reference points closest to the target mine according to the distance between the catchment area and the target mine and the terrain height between the catchment area and the target mine; S2: calculating the aeration zone diffusion flux between the position of the catchment area where the reference point closest to the target mine is located and the target mine, screening the influence catchment area of the target mine, and calculating the water quality fluctuation coefficient of the catchment area where the reference point closest to the target mine is located on the influence catchment area; S3: calculating the influence weight of the influence catchment area on the target mine according to the water quality fluctuation coefficient of the target mine and the water quality fluctuation coefficient of the catchment area where the reference point is located on the influence catchment area; S4: constructing a water balance equation of each influence catchment area, a groundwater control equation related to the target mine water and the groundwater, and discretizing the groundwater control equation to obtain a coupled equation set; The amount of water infiltrated into the groundwater in the influence catchment area is taken as a groundwater source and sink item and is input into the coupled equation set, so as to simulate the change state of the mine water under the change of the groundwater level and the catchment storage capacity, and plan the water quantity that can be dispatched at the current simulation scheduling time of the target mine.

2. The method of claim 1, wherein, The step S1 comprises: S11: Obtain the catchment area around the target mine in the study area, and obtain the coordinates of each reference point in the overhead plane of the study area , and the coordinates of the point where the target mine is located in the overhead plane of the study area , i is the number of the catchment area in the study area, e is the number of the reference point on the catchment area; S12: Extract the elevation of the riverbed low point at the location of the reference point and compare it with the elevation of the target mine bottom , respectively the elevation of the target mine mouth, the vertical depth of the target mine mouth; If the reference point is removed from the catchment, otherwise it is retained; S13: according to the coordinates of the reserved reference points , u number the reserved reference points, and screen the reference points closest to the target mine according to the straight-line distance between the reserved reference points and the target mine; ； wherein, is a set of reference points retained on the basin.

3. The method of simulating scheduling of surface water-groundwater-mine water coupling according to claim 2, characterized in that, The step S2 comprises: S21: calculating the reference point closest to the target mine from the screened reference points aerobic zone diffusion flux between the location of the river basin and the target mine ; ； wherein, is a reference point aeration zone connectivity between the catchment where the reference point is located and the target mine, is a permeability coefficient of the aeration zone, is a reference point a water level of the catchment where the reference point is located, is a reference point an aeration zone top elevation between the catchment where the reference point is located and the target mine, is an aeration zone thickness, is a reference point a water level change driving coefficient of the catchment where the reference point is located, is a water level change amount of the catchment where the reference point is located within a time , is a water level change amount of the catchment where the reference point is located within a time , is an annual water level average change amount of the catchment where the reference point is located; S22: setting a threshold value of the diffusion flux ; like Then determine the reference point. River Basin i It is the watershed affected by the target mine; If , it is determined that the reference point is located in the basin i that is not the influence basin of the target mine; S23: Obtain all impact basins of the target mine within the study area and construct an impact basin set. Collect data on the impact of the watershed. Each affected watershed a upper reference point The watershed is based on time series. Water quality data , This is the current simulated scheduling time. For the current simulated scheduling time Water quality data; S24: according to water quality data computing impact basins a upper reference point where the basin is in time series coefficient of water quality fluctuation ; ； wherein b is a time series of time instant numbers, is any time instant in the time series of time instant numbers.

4. The simulated scheduling method of surface water-groundwater-mine water coupling according to claim 3, characterized in that, The step S3 comprises: S31: Collect time series Water quality data of the target mine And calculate the water quality fluctuation coefficient of the target mine ; ； wherein, is the current simulated dispatch time water quality data for the target mine; S32: According to the water quality fluctuation coefficient With the water quality fluctuation coefficient Calculate the impact basin a The impact weight of the target mine ; ； wherein, A is the number of impact basins in the basin set is the number of impact basins in the basin set 5. The method of simulated scheduling of surface water-groundwater-mine water coupling according to claim 4, characterized in that, The step S4 comprises: S41: constructing an impact watershed set the water balance equation of each impact watershed a ; ； in, To affect the watershed a Reference point The watershed is located at the current simulated scheduling time. Water storage capacity To affect the watershed a Reference point The watershed at the previous moment Water storage capacity To affect the watershed a Reference point Water level changes in the basin For reference point The area of ​​the watershed, Reference points The inflow and outflow of water in the basin Surface runoff, The amount of water within the watershed that seeps into groundwater; S42: considering the dynamic sink item of mine water inflow, constructing a groundwater control equation related to the target mine water and the groundwater; ； wherein, respectively, the permeability coefficient of groundwater in x , y , the permeability coefficient in z direction, is the groundwater level, W is the groundwater source and sink term, is the mine water inflow sink term, is the water yield of groundwater; S43: discretize the groundwater control equation, time variable t discretized into an implicit format to obtain a discretized groundwater control equation; ； wherein, is a discrete spatial term, is a time step, are simulated scheduled times , of the groundwater level; S44: obtaining a coupled equation set by simultaneously discretizing the groundwater control equation and the water balance equation; ； amount of water infiltrating into groundwater within a basin affected by the event as a groundwater source sink term W ; at this time ; previous simulation dispatch time a basin affected by the event a reference point amount of water infiltrating into groundwater within a basin affected by the event The water storage capacity of the catchment area where the reference point of each influence catchment area is located, the groundwater level around the target mine at the current simulation scheduling time are collected, the mine water inflow sink item is calculated by using the coupled equation set, and the change state of the mine water under the change of the groundwater level and the catchment storage capacity is simulated; S45: planning the water quantity that can be dispatched at the current simulation scheduling time of the target mine according to the calculated mine water inflow sink item.

6. A simulation scheduling system applying the simulation scheduling method of surface water-groundwater-mine water coupling according to claim 5, characterized in that, It comprises: A water quality sensor for collecting water quality data in the influence catchment area and the target mine water; A water level sensor for collecting water level data in the influence catchment area and the target mine water; A data processing device is provided with a memory and a processor, the memory is used for storing computer program instructions, and the computer program instructions are executed by the processor to realize the simulation scheduling method.

Images