Analysis method of factors affecting water inrush and sand burst in high water pressure loose aquifer

Abstract

The invention discloses an analysis method for factors affecting water inrush and sand burst during mining in a loose aquifer with high water pressure, relates to the field of geological assessment, and solves the problem that the existing analysis method for factors affecting water inrush and sand burst has poor effect when conducting geological assessment. The method comprises: step S1: acquiring mining stage division data, regional real-time operation coefficients, and high-density foundation pit areas to obtain regional mining operation data; step S2: analyzing the regional mining operation data to obtain characteristic foundation pit quantity-area ratios and regional hydrological risk coefficients to obtain regional hydrological analysis data; step S3: analyzing the regional mining operation data and the regional hydrological analysis data to obtain the regional geological analysis coefficient; step S4: conducting an operation risk assessment on a target mining operation area according to the regional mining operation data, the regional hydrological analysis data, and the regional geological risk coefficient. The invention can improve the accuracy and pertinence of the analysis method for factors affecting water inrush and sand burst during geological assessment.

Description

Technical Field

[0001] The invention belongs to the field of geological assessment and relates to data analysis technology, specifically to an analysis method for factors affecting water inrush and sand burst in mining under high water pressure loose aquifers. Background Art

[0002] The existing analysis methods of factors affecting water inrush and sand burst have the following defects when evaluating water inrush and sand burst accidents:

[0003] 1. When evaluating water inrush and sand burst accidents, the existing analysis methods of factors affecting water inrush and sand burst usually only consider predictions before engineering operations, and cannot make real-time predictions based on the progress of the project, which results in a lack of real-time and accuracy in the prediction results.

[0004] 2. When evaluating water inrush and sand burst accidents, the existing analysis methods of factors affecting water inrush and sand burst usually require comprehensive predictions for all areas involved in the engineering operation, and fail to select accident-prone areas for targeted predictions based on the progress of the project. This results in a lack of pertinence in the prediction process, which can easily lead to a waste of exploration resources.

[0005] To this end, we proposed an analysis method for factors affecting water inrush and sand burst during mining in loose aquifers with high water pressure. Summary of the invention

[0006] In view of the shortcomings of the prior art, the purpose of the present invention is to provide an analysis method for factors affecting water inrush and sand burst in mining in loose aquifers with high water pressure. The present invention is based on respectively obtaining mining stage division data, regional historical operation coefficient, regional real-time operation coefficient and high-density foundation pit area, and obtaining regional mining risk coefficient by analyzing regional historical operation coefficient and regional real-time operation coefficient. The mining stage division data, regional real-time operation coefficient and high-density foundation pit area are defined as regional mining operation data. By analyzing the regional mining operation data, the foundation pits in the high-density foundation pit area are divided into first type foundation pits and second type foundation pits, and a hydrological analysis is performed on the second type foundation pits to obtain the characteristic foundation pit quantity area ratio, the average sudden inrush depth value of the foundation pit and The average sudden surge pressure value of the foundation pit is obtained by analyzing the area ratio of the number of characteristic foundation pits, the average sudden surge depth value of the foundation pit and the average sudden surge pressure value of the foundation pit, and the regional hydrological risk coefficient is defined as the regional hydrological analysis data. The average rock porosity, foundation pit depth value and foundation pit air humidity value of the high-density foundation pit area are obtained by analyzing the regional mining operation data and the regional hydrological analysis data. The regional geological analysis coefficient is obtained by calculating the average rock porosity, foundation pit depth value and foundation pit air humidity value of the foundation pit. According to the regional mining operation data, the regional hydrological analysis data and the regional geological risk coefficient, the target mining operation area is evaluated for operational risk.

[0007] In order to achieve the above-mentioned purpose, the present invention adopts the following technical scheme: The method for analyzing factors affecting water inrush and sand burst in high-pressure loose aquifers also includes the following specific steps:

[0008] Step S1: respectively obtaining the mining stage division data, the regional historical operation coefficient, the regional real-time operation coefficient and the high-density foundation pit area, and obtaining the regional mining risk coefficient by analyzing the regional historical operation coefficient and the regional real-time operation coefficient, and defining the mining stage division data, the regional real-time operation coefficient and the high-density foundation pit area as the regional mining operation data;

[0009] Step S2: by analyzing the regional mining operation data, the foundation pits in the high-density foundation pit area are divided into the first type of foundation pits and the second type of foundation pits, and the second type of foundation pits are subjected to hydrological analysis to obtain the characteristic foundation pit quantity area ratio, the average sudden gushing depth value of the foundation pits, and the average sudden gushing pressure value of the foundation pits; by analyzing the characteristic foundation pit quantity area ratio, the average sudden gushing depth value of the foundation pits, and the average sudden gushing pressure value of the foundation pits, the regional hydrological risk coefficient is obtained, and the characteristic foundation pit quantity area ratio and the regional hydrological risk coefficient are defined as regional hydrological analysis data;

[0010] Step S3: by analyzing the regional mining operation data and the regional hydrological analysis data, the average rock porosity, the depth value and the air humidity value of the foundation pit corresponding to the high-density foundation pit area are obtained, and the average rock porosity, the depth value and the air humidity value of the foundation pit are calculated to obtain the regional geological analysis coefficient;

[0011] Step S4: Conducting an operation risk assessment on the target mining operation area based on the regional mining operation data, regional hydrological analysis data and regional geological risk factor.

[0012] Furthermore, the step S1 further includes the following specific steps:

[0013] Step S11: Acquire regional mining operation engineering data, divide the regional mining operation engineering into b mining operation stages in chronological order according to the regional mining operation engineering data, and name them as the first mining operation stage to the bth mining operation stage, respectively, to obtain mining stage division data;

[0014] Step S12: when the mining operation is in the ath mining operation stage, historical operation data of the target mining operation area is analyzed to obtain the regional historical operation coefficient;

[0015] Step S13: when the mining operation is in the ath mining operation stage, real-time operation data analysis is performed on the target mining operation area to obtain the regional real-time operation coefficient;

[0016] Step S14: summing the regional historical operation coefficient and the regional real-time operation coefficient to obtain the regional mining risk coefficient;

[0017] Step S15: defining the regional mining risk coefficient and the high-density foundation pit area as regional mining operation data.

[0018] Furthermore, the step S12 further includes the following specific steps:

[0019] Step S121: obtaining the rock stratum mining depth corresponding to the target mining operation area before the start of the ath mining operation stage to obtain the first rock stratum mining depth, obtaining the rock stratum target mining depth corresponding to the target mining operation area in the ath mining operation stage to obtain the second rock stratum mining depth, calculating the difference between the second rock stratum mining depth and the first rock stratum mining depth, and taking the absolute value of the obtained difference to obtain the stage rock stratum mining depth;

[0020] Step S122: obtaining the rock stratum mining area corresponding to the target mining operation area before the start of the ath mining operation stage to obtain the first rock stratum mining area, obtaining the rock stratum target mining area corresponding to the target mining operation area in the ath mining operation stage to obtain the second rock stratum mining area, calculating the difference between the first rock stratum mining area and the second rock stratum mining area, and taking the absolute value of the obtained difference to obtain the stage rock stratum mining area;

[0021] Step S123: obtaining the number of mining foundation pits corresponding to the target mining operation area before the start of the ath mining operation stage, obtaining the first foundation pit number value, obtaining the target mining foundation pit number corresponding to the target mining operation area in the ath mining operation stage, obtaining the second foundation pit number;

[0022] Step S124: Calculate the ratio of the first foundation pit quantity value to the first rock stratum mining area to obtain the first foundation pit quantity area ratio, calculate the ratio of the second foundation pit quantity value to the second rock stratum mining area to obtain the second foundation pit quantity area ratio, calculate the difference between the first foundation pit quantity area ratio and the second foundation pit area quantity ratio, and take the absolute value of the obtained difference to obtain the stage foundation pit area quantity ratio;

[0023] Step S125: Calculate the regional historical operation coefficient corresponding to the a-th mining operation stage by using the ratio of the stage rock stratum mining depth, the stage rock stratum mining area and the stage foundation pit area;

[0024] The regional historical operation coefficient corresponding to the a-th mining operation stage is calculated as follows:

[0025] ;

[0026] Among them, Ygz is the regional historical operation coefficient, Jsd is the stage rock stratum mining depth, Jmj is the stage rock stratum mining area, and Mjb is the stage foundation pit area quantity ratio.

[0027] Furthermore, the step S13 further includes the following specific steps:

[0028] Step S131: when performing mining operations on the target operation area, a plan view of the target operation area is obtained in real time, and characteristic circles are used to fill the foundation pit areas in the plan view of the target operation area to obtain a plurality of characteristic circles;

[0029] Step S132: Acquire the high-density foundation pit area;

[0030] The step S132 further includes the following specific steps:

[0031] Step S1321: Take the lower left corner vertex of the target operation area plan as the coordinate origin, take the lower frame of the target operation area plan as the coordinate X axis, and draw a straight line perpendicular to the coordinate X axis through the coordinate origin as the coordinate Y axis to obtain the coordinate system of the lower mining operation area;

[0032] Step S1322: Mark each characteristic circle in the lower mining operation area coordinate system, and obtain the center coordinates of each characteristic circle in the lower mining operation area coordinate system to obtain a plurality of center coordinates;

[0033] Step S1323: randomly selecting a coordinate from a plurality of circle center coordinates as a sample circle center coordinate;

[0034] Step S1324: Calculate the density value corresponding to the sample circle center coordinates through the sample circle center coordinates and multiple circle center coordinates;

[0035] Calculate the density value corresponding to the sample center coordinates. The specific formula is as follows:

[0036] ;

[0037] Among them, Mz is the density value corresponding to the coordinates of the sample center, (x i ,y i ) is the sample center coordinate, K() is the characteristic Gaussian kernel, n is the number value corresponding to the characteristic circle in the coordinate system of the down-sampling operation area, x is the horizontal coordinate of the arbitrary circle center coordinate, y is the vertical coordinate of the arbitrary circle center coordinate, and h is the bandwidth parameter set in the characteristic Gaussian kernel;

[0038] The characteristic Gaussian kernel is calculated as follows:

[0039] ;

[0040] Among them, K(u) is the characteristic Gaussian kernel, u is the standard deviation of the set Gaussian kernel, and e is a natural constant;

[0041] Step S1325: Obtain the density value corresponding to each circle center coordinate respectively;

[0042] Step S1326: Obtain a circle center coordinate density threshold, take the circle center coordinates with a density value greater than the circle center coordinate density threshold as the first type of circle center coordinates, take the circle center coordinates with a density value less than or equal to the circle center coordinate density threshold as the second type of circle center coordinates, and mark the foundation pit area consisting of foundation pits corresponding to multiple first type circle center coordinates as a high-density foundation pit area;

[0043] Step S133: acquiring a depth value corresponding to each foundation pit in the high-density foundation pit area to obtain a plurality of foundation pit depth values, and marking the foundation pit depth value with the largest value as the peak foundation pit depth;

[0044] Step S134: Calculate the average of multiple foundation pit depth values ​​to obtain an average foundation pit depth;

[0045] Step S135: when the a-th mining operation stage is in operation, respectively obtain the foundation pit vibration amplitude value corresponding to each foundation pit in the high-density foundation pit area, obtain multiple vibration amplitude values, and calculate the average of them to obtain the average amplitude of the foundation pit area;

[0046] Step S136: Calculate the regional real-time operation coefficient by using the data of the mining stage division, the peak foundation pit depth, the average foundation pit depth and the average amplitude of the foundation pit area;

[0047] The regional real-time operation coefficient is calculated, and the specific formula is as follows:

[0048] ;

[0049] Among them, Ssx is the regional real-time operation coefficient, Fjs is the peak foundation pit depth, Pjs is the average foundation pit depth, and Kzf is the average amplitude of the foundation pit area.

[0050] Furthermore, the step S2 further includes the following specific steps:

[0051] Step S21: Acquire regional mining operation data, and acquire high-density foundation pit area and mining stage division data according to the regional mining operation data;

[0052] Step S22: when the mining operation is in the ath mining operation stage, obtaining the regional water content coefficient and the characteristic foundation pit quantity area ratio corresponding to the high-density foundation pit area;

[0053] Step S23: Obtaining the regional hydrological risk coefficient;

[0054] Step S24: define the regional hydrological risk coefficient and the ratio of the number and area of ​​characteristic foundation pits as regional hydrological analysis data.

[0055] Furthermore, the step S22 further includes the following specific steps:

[0056] Step S221: randomly selecting a foundation pit in the high-density foundation pit area as a water-containing monitoring sample foundation pit;

[0057] Step S222: Drilling downward from the bottom of the water-containing monitoring sample foundation pit to a target drilling depth, and if no groundwater is found, marking the water-containing monitoring sample foundation pit as a first type foundation pit;

[0058] Step S223: Drilling downward from the bottom of the water-containing monitoring sample foundation pit to the target drilling depth, if groundwater is explored, the water-containing monitoring sample foundation pit is marked as a second type of foundation pit, the drilling depth when groundwater and soil surge occurs is obtained, and the groundwater depth value corresponding to the water-containing monitoring sample foundation pit is obtained, and the surge pressure value when groundwater surges occurs is obtained, and the groundwater surge pressure value corresponding to the water-containing monitoring sample foundation pit is obtained;

[0059] Step S224: acquiring the number of second-type foundation pits to obtain the number value of second-type foundation pits, acquiring the area value of the high-density foundation pit area to obtain the high-density foundation pit area value, and calculating the ratio of the second-type foundation pit number value to the high-density foundation pit area value to obtain the characteristic foundation pit number area ratio;

[0060] Step S225: performing water content monitoring on each foundation pit in the high-density foundation pit area to obtain a plurality of second-type foundation pits, obtaining a groundwater depth value corresponding to each second-type foundation pit to obtain a plurality of groundwater depth values, and averaging the plurality of groundwater depth values ​​to obtain an average sudden inrush depth value of the foundation pit;

[0061] Step S226: acquiring the groundwater surge pressure value corresponding to each second type foundation pit respectively to obtain a plurality of groundwater surge pressure values, and averaging the plurality of groundwater surge pressure values ​​to obtain an average foundation pit surge pressure value;

[0062] Step S227: Calculate the area ratio of the characteristic foundation pit quantity, the average sudden gushing depth value of the foundation pit, and the average sudden gushing pressure value of the foundation pit to obtain the regional water content coefficient corresponding to the a-th mining operation stage;

[0063] The regional moisture coefficient is calculated using the following formula:

[0064] ;

[0065] Among them, Qhsa is the regional moisture coefficient corresponding to the ath mining operation stage, Tsz is the average sudden surge depth of the foundation pit, Tyz is the average sudden surge pressure of the foundation pit, and Tmb is the ratio of the number and area of ​​characteristic foundation pits.

[0066] Furthermore, the step S23 further includes the following specific steps:

[0067] Step S231: when the ath mining operation stage is the first mining operation stage, the regional water content coefficient corresponding to the ath mining operation stage is defined as the regional hydrological risk coefficient;

[0068] Step S232: when the ath mining operation stage is not the first mining operation stage, the regional water content coefficient corresponding to the a-1th mining operation stage is obtained, and the regional water content coefficient corresponding to the ath mining operation stage and the regional water content coefficient corresponding to the a-1th mining operation stage are calculated to obtain the regional hydrological risk coefficient;

[0069] The regional hydrological risk coefficient is calculated using the following formula:

[0070] ;

[0071] Among them, Qfx is the regional hydrological risk coefficient, Qhsa is the regional moisture coefficient corresponding to the a-th mining operation stage, and Qhs (a-1) is the regional moisture coefficient corresponding to the a-1-th mining operation stage.

[0072] Furthermore, the step S3 further includes the following specific steps:

[0073] Step S31: Acquire regional mining operation data, and acquire a high-density foundation pit area according to the regional mining operation data;

[0074] Step S32: obtaining regional hydrological analysis data, and obtaining characteristic foundation pit quantity and area ratio according to the regional hydrological analysis data;

[0075] Step S33: randomly selecting a foundation pit in the high-density foundation pit area as a geological monitoring sample foundation pit;

[0076] Step S34: Perform geological monitoring on the geological monitoring sample foundation pit to obtain the foundation pit geological analysis coefficient corresponding to the geological monitoring sample foundation pit;

[0077] Step S35: In the high-density foundation pit area, the foundation pit geological analysis coefficient corresponding to each foundation pit is obtained to obtain multiple foundation pit geological analysis coefficients, and the obtained foundation pit geological analysis coefficients are averaged to obtain the regional geological risk coefficient.

[0078] Furthermore, the step S34 further includes the following specific steps:

[0079] Step S341: select a rock mass of unit mass in the geological monitoring sample foundation pit as the sample monitoring rock mass, place the sample monitoring rock mass in the water body, make the sample monitoring rock mass in a completely submerged state, wait for the submerged state to reach the characteristic submerged time, take out the sample monitoring rock mass to obtain the wet sample monitoring rock mass, obtain the mass value of the wet sample monitoring rock mass, obtain the volume value corresponding to the wet sample monitoring rock mass by the water displacement method, and obtain the volume value of the sample monitoring rock mass;

[0080] Step S342: placing the sample monitoring rock mass in a characteristic temperature environment, and obtaining the mass value of the sample monitoring rock mass in real time until the mass value of the sample monitoring rock mass maintains a constant weight, obtaining a dry sample monitoring rock mass, and obtaining the mass of the dry sample monitoring rock mass to obtain the mass value of the dry sample rock mass;

[0081] Step S343: The dry sample rock mass value, the sample monitoring rock mass volume value and the wet sample rock mass value are calculated to obtain the rock mass porosity corresponding to the sample monitoring rock mass;

[0082] The rock mass porosity corresponding to the sample monitoring rock mass is calculated, and the specific formula is as follows:

[0083] ;

[0084] Among them, Kxd is the rock mass porosity corresponding to the sample monitoring rock mass, Ms is the wet sample rock mass value, Mg is the dry sample rock mass value, Vy is the sample monitoring rock mass volume value, and ρ is the water density value;

[0085] Step S344: randomly selecting multiple monitoring rock masses in the geological monitoring sample foundation pit, respectively acquiring the rock mass porosity corresponding to each monitoring rock mass, obtaining multiple rock mass porosities, and averaging the obtained multiple rock mass porosities to obtain the average rock mass porosity of the foundation pit;

[0086] Step S345: Acquire the air humidity of the geological monitoring sample foundation pit to obtain the foundation pit air humidity value, and acquire the depth value of the geological monitoring sample foundation pit to obtain the foundation pit depth value;

[0087] Step S346: Calculate the average rock porosity of the foundation pit, the foundation pit depth value, and the foundation pit air humidity value to obtain the foundation pit geological analysis coefficient corresponding to the geological monitoring sample foundation pit;

[0088] The geological analysis coefficient of the foundation pit is calculated, and the specific formula is configured as follows:

[0089] ;

[0090] Among them, Jkx is the geological analysis coefficient of the foundation pit, Sdz is the foundation pit depth value, Wdz is the foundation pit air humidity value, and Kxp is the average rock porosity of the foundation pit.

[0091] Furthermore, the step S4 further includes the following specific steps:

[0092] Step S41: Acquire regional mining operation data, and acquire regional mining risk coefficient according to the regional mining operation data;

[0093] Step S42: obtaining regional hydrological analysis data, and obtaining a regional hydrological risk coefficient according to the regional hydrological analysis data;

[0094] Step S43: Obtaining regional geological risk coefficient;

[0095] Step S44: Calculate the regional mining risk coefficient, the regional hydrological risk coefficient and the regional geological risk coefficient to obtain an overall operation risk coefficient;

[0096] The overall coefficient of operation risk is calculated as follows:

[0097] ;

[0098] Among them, Fxz is the overall coefficient of operation risk, Qcz is the regional mining risk coefficient, Qfx is the regional hydrological risk coefficient, and Qdz is the regional geological risk coefficient;

[0099] Step S45: obtaining an overall operation risk coefficient threshold, comparing the overall operation risk coefficient with the overall operation risk coefficient threshold, and conducting an operation risk assessment on the target mining operation area according to the numerical comparison result;

[0100] The step S45 further includes the following specific steps:

[0101] Step S451: respectively obtaining a regional mining risk coefficient threshold, a regional hydrological risk coefficient threshold, and a regional geological risk coefficient threshold;

[0102] Step S452: Calculate the regional mining risk coefficient threshold, the regional hydrological risk coefficient threshold, and the regional geological risk coefficient threshold to obtain an overall operation risk coefficient threshold;

[0103] Step S453: when the overall coefficient of operation risk is greater than or equal to the overall coefficient threshold of operation risk, the target mining operation area is evaluated as a high water inrush and sand burst risk area;

[0104] Step S454: When the overall operation risk coefficient is less than the overall operation risk coefficient threshold, the target mining operation area is evaluated as a low water inrush and sand burst risk area.

[0105] In summary, due to the adoption of the above technical solution, the beneficial effects of the present invention are:

[0106] 1. The present invention conducts safety assessment on water inrush and sand burst risks by obtaining regional historical operation coefficients and regional real-time operation coefficients, and can combine historical data with real-time data, thereby improving the accuracy of water inrush and sand burst accident assessment in the target mining operation area;

[0107] 2. The present invention can effectively improve the regional targeting of the exploration process and save exploration resources by marking the high-density foundation pit area and focusing on the hydrological exploration and geological exploration in the high-density foundation pit area. BRIEF DESCRIPTION OF THE DRAWINGS

[0108] In order to facilitate understanding by those skilled in the art, the present invention is further described below with reference to the accompanying drawings.

[0109] Figure 1 It is a diagram of the implementation steps of the present invention;

[0110] Figure 2 is the overall system block diagram of the present invention;

[0111] Figure 3 It is a plan view of the target operation area in the present invention. DETAILED DESCRIPTION

[0112] The technical solution of the present invention will be clearly and completely described below in conjunction with the embodiments. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention. Embodiment 1

[0113] See also Figure 1 The present invention provides a technical solution: a method for analyzing factors affecting water inrush and sand burst in high-pressure loose aquifers, comprising the following specific steps:

[0114] Step S1: respectively obtaining the mining stage division data, the regional historical operation coefficient, the regional real-time operation coefficient and the high-density foundation pit area, and obtaining the regional mining risk coefficient by analyzing the regional historical operation coefficient and the regional real-time operation coefficient, and defining the mining stage division data, the regional real-time operation coefficient and the high-density foundation pit area as the regional mining operation data;

[0115] Step S11: Acquire regional mining operation engineering data, divide the regional mining operation engineering into b mining operation stages in chronological order according to the regional mining operation engineering data, and name them as the first mining operation stage to the bth mining operation stage, respectively, to obtain mining stage division data;

[0116] Step S12: when the mining operation is in the ath mining operation stage, historical operation data of the target mining operation area is analyzed to obtain the regional historical operation coefficient;

[0117] The step S12 further includes the following specific steps:

[0118] Step S121: obtaining the rock stratum mining depth corresponding to the target mining operation area before the start of the ath mining operation stage to obtain the first rock stratum mining depth, obtaining the rock stratum target mining depth corresponding to the target mining operation area in the ath mining operation stage to obtain the second rock stratum mining depth, calculating the difference between the second rock stratum mining depth and the first rock stratum mining depth, and taking the absolute value of the obtained difference to obtain the stage rock stratum mining depth;

[0119] Step S122: obtaining the rock stratum mining area corresponding to the target mining operation area before the start of the ath mining operation stage to obtain the first rock stratum mining area, obtaining the rock stratum target mining area corresponding to the target mining operation area in the ath mining operation stage to obtain the second rock stratum mining area, calculating the difference between the first rock stratum mining area and the second rock stratum mining area, and taking the absolute value of the obtained difference to obtain the stage rock stratum mining area;

[0120] Step S123: obtaining the number of mining foundation pits corresponding to the target mining operation area before the start of the ath mining operation stage, obtaining the first foundation pit number value, obtaining the target mining foundation pit number corresponding to the target mining operation area in the ath mining operation stage, obtaining the second foundation pit number;

[0121] Step S124: Calculate the ratio of the first foundation pit quantity value to the first rock stratum mining area to obtain the first foundation pit quantity area ratio, calculate the ratio of the second foundation pit quantity value to the second rock stratum mining area to obtain the second foundation pit quantity area ratio, calculate the difference between the first foundation pit quantity area ratio and the second foundation pit area quantity ratio, and take the absolute value of the obtained difference to obtain the stage foundation pit area quantity ratio;

[0122] Step S125: Calculate the regional historical operation coefficient corresponding to the a-th mining operation stage by using the ratio of the stage rock stratum mining depth, the stage rock stratum mining area and the stage foundation pit area;

[0123] The regional historical operation coefficient corresponding to the a-th mining operation stage is calculated as follows:

[0124] ;

[0125] Among them, Ygz is the regional historical operation coefficient, Jsd is the stage rock mining depth, Jmj is the stage rock mining area, and Mjb is the stage foundation pit area quantity ratio;

[0126] Step S13: when the mining operation is in the ath mining operation stage, real-time operation data analysis is performed on the target mining operation area to obtain the regional real-time operation coefficient;

[0127] The step S13 further includes the following specific steps:

[0128] Step S131: when performing mining operations on the target operation area, a plan view of the target operation area is obtained in real time, and characteristic circles are used to fill the foundation pit areas in the plan view of the target operation area to obtain a plurality of characteristic circles;

[0129] Step S132: Acquire the high-density foundation pit area;

[0130] The step S132 further includes the following specific steps:

[0131] Step S1321: Take the lower left corner vertex of the target operation area plan as the coordinate origin, take the lower frame of the target operation area plan as the coordinate X axis, and draw a straight line perpendicular to the coordinate X axis through the coordinate origin as the coordinate Y axis to obtain the coordinate system of the lower mining operation area;

[0132] Step S1322: Mark each characteristic circle in the lower mining operation area coordinate system, and obtain the center coordinates of each characteristic circle in the lower mining operation area coordinate system to obtain a plurality of center coordinates;

[0133] Step S1323: randomly selecting a coordinate from a plurality of circle center coordinates as a sample circle center coordinate;

[0134] Step S1324: Calculate the density value corresponding to the sample circle center coordinates through the sample circle center coordinates and multiple circle center coordinates;

[0135] Calculate the density value corresponding to the sample center coordinates. The specific formula is as follows:

[0136] ;

[0137] Among them, Mz is the density value corresponding to the coordinates of the sample center, (x i ,y i ) is the sample center coordinate, K() is the characteristic Gaussian kernel, n is the number value corresponding to the characteristic circle in the coordinate system of the down-sampling operation area, x is the horizontal coordinate of the arbitrary circle center coordinate, y is the vertical coordinate of the arbitrary circle center coordinate, and h is the bandwidth parameter set in the characteristic Gaussian kernel;

[0138] Step S1325: Obtain the density value corresponding to each circle center coordinate respectively;

[0139] Step S1326: Obtain a circle center coordinate density threshold, take the circle center coordinates with a density value greater than the circle center coordinate density threshold as the first type of circle center coordinates, take the circle center coordinates with a density value less than or equal to the circle center coordinate density threshold as the second type of circle center coordinates, and mark the foundation pit area consisting of foundation pits corresponding to multiple first type circle center coordinates as a high-density foundation pit area;

[0140] Step S133: acquiring a depth value corresponding to each foundation pit in the high-density foundation pit area to obtain a plurality of foundation pit depth values, and marking the foundation pit depth value with the largest value as the peak foundation pit depth;

[0141] Step S134: Calculate the average of multiple foundation pit depth values ​​to obtain an average foundation pit depth;

[0142] Step S135: when the a-th mining operation stage is in operation, respectively obtain the foundation pit vibration amplitude value corresponding to each foundation pit in the high-density foundation pit area, obtain multiple vibration amplitude values, and calculate the average of them to obtain the average amplitude of the foundation pit area;

[0143] Step S136: Calculate the regional real-time operation coefficient by using the data of the mining stage division, the peak foundation pit depth, the average foundation pit depth and the average amplitude of the foundation pit area;

[0144] The regional real-time operation coefficient is calculated, and the specific formula is as follows:

[0145] ;

[0146] Among them, Ssx is the regional real-time operation coefficient, Fjs is the peak foundation pit depth, Pjs is the average foundation pit depth, and Kzf is the average amplitude of the foundation pit area;

[0147] Step S14: summing the regional historical operation coefficient and the regional real-time operation coefficient to obtain the regional mining risk coefficient;

[0148] Step S15: defining the regional mining risk coefficient and the high-density foundation pit area as regional mining operation data;

[0149] Step S2: performing a regional water content analysis on the target operation area according to the regional down-sampling operation data to obtain regional hydrological analysis data;

[0150] Step S21: Acquire regional mining operation data, and acquire high-density foundation pit area and mining stage division data according to the regional mining operation data;

[0151] Step S22: when the mining operation is in the ath mining operation stage, obtaining the regional water content coefficient corresponding to the high-density foundation pit area;

[0152] The step S22 further includes the following specific steps:

[0153] Step S221: randomly selecting a foundation pit in the high-density foundation pit area as a water-containing monitoring sample foundation pit;

[0154] Step S222: Drilling downward from the bottom of the water-containing monitoring sample foundation pit to a target drilling depth, and if no groundwater is found, marking the water-containing monitoring sample foundation pit as a first type foundation pit;

[0155] Step S223: Drilling downward from the bottom of the water-containing monitoring sample foundation pit to the target drilling depth, if groundwater is explored, the water-containing monitoring sample foundation pit is marked as a second type of foundation pit, the drilling depth when groundwater and soil surge occurs is obtained, and the groundwater depth value corresponding to the water-containing monitoring sample foundation pit is obtained, and the surge pressure value when groundwater surges occurs is obtained, and the groundwater surge pressure value corresponding to the water-containing monitoring sample foundation pit is obtained;

[0156] Step S224: acquiring the number of second-type foundation pits to obtain the number value of second-type foundation pits, acquiring the area value of the high-density foundation pit area to obtain the high-density foundation pit area value, and calculating the ratio of the second-type foundation pit number value to the high-density foundation pit area value to obtain the characteristic foundation pit number area ratio;

[0157] Step S225: performing water content monitoring on each foundation pit in the high-density foundation pit area to obtain a plurality of second-type foundation pits, obtaining a groundwater depth value corresponding to each second-type foundation pit to obtain a plurality of groundwater depth values, and averaging the plurality of groundwater depth values ​​to obtain an average sudden inrush depth value of the foundation pit;

[0158] Step S226: acquiring the groundwater surge pressure value corresponding to each second type foundation pit respectively to obtain a plurality of groundwater surge pressure values, and averaging the plurality of groundwater surge pressure values ​​to obtain an average foundation pit surge pressure value;

[0159] Step S227: Calculate the area ratio of the characteristic foundation pit quantity, the average sudden gushing depth value of the foundation pit, and the average sudden gushing pressure value of the foundation pit to obtain the regional water content coefficient corresponding to the a-th mining operation stage;

[0160] The regional moisture coefficient is calculated using the following formula:

[0161] ;

[0162] Among them, Qhsa is the regional water content coefficient corresponding to the ath mining operation stage, Tsz is the average sudden surge depth value of the foundation pit, Tyz is the average sudden surge pressure value of the foundation pit, and Tmb is the ratio of the number and area of ​​characteristic foundation pits;

[0163] Step S23: Obtaining the regional hydrological risk coefficient;

[0164] The step S23 further includes the following specific steps:

[0165] Step S231: when the ath mining operation stage is the first mining operation stage, the regional water content coefficient corresponding to the ath mining operation stage is defined as the regional hydrological risk coefficient;

[0166] Step S232: when the ath mining operation stage is not the first mining operation stage, the regional water content coefficient corresponding to the a-1th mining operation stage is obtained, and the regional water content coefficient corresponding to the ath mining operation stage and the regional water content coefficient corresponding to the a-1th mining operation stage are calculated to obtain the regional hydrological risk coefficient;

[0167] The regional hydrological risk coefficient is calculated using the following formula:

[0168] ;

[0169] Among them, Qfx is the regional hydrological risk coefficient, Qhsa is the regional water content coefficient corresponding to the a-th mining operation stage, and Qhs (a-1) is the regional water content coefficient corresponding to the a-1-th mining operation stage;

[0170] Step S24: defining the regional hydrological risk coefficient and the characteristic foundation pit quantity-area ratio as regional hydrological analysis data;

[0171] Step S3: Perform regional geological analysis on the target operation area according to the regional mining operation data to obtain the regional geological risk coefficient;

[0172] Step S31: Acquire regional mining operation data, and acquire a high-density foundation pit area according to the regional mining operation data;

[0173] Step S32: obtaining regional hydrological analysis data, and obtaining characteristic foundation pit quantity and area ratio according to the regional hydrological analysis data;

[0174] Step S33: randomly selecting a foundation pit in the high-density foundation pit area as a geological monitoring sample foundation pit;

[0175] Step S34: Perform geological monitoring on the geological monitoring sample foundation pit to obtain the foundation pit geological analysis coefficient corresponding to the geological monitoring sample foundation pit;

[0176] The step S34 further includes the following specific steps:

[0177] Step S341: select a rock mass of unit mass in the geological monitoring sample foundation pit as the sample monitoring rock mass, place the sample monitoring rock mass in the water body, make the sample monitoring rock mass in a completely submerged state, wait for the submerged state to reach the characteristic submerged time, take out the sample monitoring rock mass to obtain the wet sample monitoring rock mass, obtain the mass value of the wet sample monitoring rock mass, obtain the volume value corresponding to the wet sample monitoring rock mass by the water displacement method, and obtain the volume value of the sample monitoring rock mass;

[0178] Step S342: placing the sample monitoring rock mass in a characteristic temperature environment, and obtaining the mass value of the sample monitoring rock mass in real time until the mass value of the sample monitoring rock mass maintains a constant weight, obtaining a dry sample monitoring rock mass, and obtaining the mass of the dry sample monitoring rock mass to obtain the mass value of the dry sample rock mass;

[0179] Step S343: The dry sample rock mass value, the sample monitoring rock mass volume value and the wet sample rock mass value are calculated to obtain the rock mass porosity corresponding to the sample monitoring rock mass;

[0180] The rock mass porosity corresponding to the sample monitoring rock mass is calculated, and the specific formula is as follows:

[0181] ;

[0182] Among them, Kxd is the rock mass porosity corresponding to the sample monitoring rock mass, Ms is the wet sample rock mass value, Mg is the dry sample rock mass value, Vy is the sample monitoring rock mass volume value, and ρ is the water density value;

[0183] Step S344: randomly selecting multiple monitoring rock masses in the geological monitoring sample foundation pit, respectively acquiring the rock mass porosity corresponding to each monitoring rock mass, obtaining multiple rock mass porosities, and averaging the obtained multiple rock mass porosities to obtain the average rock mass porosity of the foundation pit;

[0184] Step S345: Acquire the air humidity of the geological monitoring sample foundation pit to obtain the foundation pit air humidity value, and acquire the depth value of the geological monitoring sample foundation pit to obtain the foundation pit depth value;

[0185] Step S346: Calculate the average rock porosity of the foundation pit, the foundation pit depth value, and the foundation pit air humidity value to obtain the foundation pit geological analysis coefficient corresponding to the geological monitoring sample foundation pit;

[0186] The geological analysis coefficient of the foundation pit is calculated, and the specific formula is configured as follows:

[0187] ;

[0188] Among them, Jkx is the geological analysis coefficient of the foundation pit, Sdz is the foundation pit depth value, Wdz is the foundation pit air humidity value, and Kxp is the average rock porosity of the foundation pit;

[0189] Step S35: in the high-density foundation pit area, the foundation pit geological analysis coefficient corresponding to each foundation pit is obtained to obtain a plurality of foundation pit geological analysis coefficients, and the obtained foundation pit geological analysis coefficients are averaged to obtain the regional geological risk coefficient;

[0190] Step S4: conducting an operation risk assessment on the target mining operation area according to the regional mining operation data, regional hydrological analysis data and regional geological risk coefficient;

[0191] Step S41: Acquire regional mining operation data, and acquire regional mining risk coefficient according to the regional mining operation data;

[0192] Step S42: obtaining regional hydrological analysis data, and obtaining a regional hydrological risk coefficient according to the regional hydrological analysis data;

[0193] Step S43: Obtaining regional geological risk coefficient;

[0194] Step S44: Calculate the regional mining risk coefficient, the regional hydrological risk coefficient and the regional geological risk coefficient to obtain an overall operation risk coefficient;

[0195] The overall coefficient of operation risk is calculated as follows:

[0196] ;

[0197] Among them, Fxz is the overall coefficient of operation risk, Qcz is the regional mining risk coefficient, Qfx is the regional hydrological risk coefficient, and Qdz is the regional geological risk coefficient;

[0198] Step S45: obtaining an overall operation risk coefficient threshold, comparing the overall operation risk coefficient with the overall operation risk coefficient threshold, and conducting an operation risk assessment on the target mining operation area according to the numerical comparison result;

[0199] The step S45 further includes the following specific steps:

[0200] Step S451: respectively obtaining a regional mining risk coefficient threshold, a regional hydrological risk coefficient threshold, and a regional geological risk coefficient threshold;

[0201] Step S452: Calculate the regional mining risk coefficient threshold, the regional hydrological risk coefficient threshold, and the regional geological risk coefficient threshold to obtain an overall operation risk coefficient threshold;

[0202] The overall coefficient threshold of the operation risk is calculated, and the specific formula is as follows:

[0203] ;

[0204] Among them, Fxzy is the overall coefficient threshold of operation risk, Qczy is the regional mining risk coefficient threshold, Qfxy is the regional hydrological risk coefficient threshold, and Qdzy is the regional geological risk coefficient threshold;

[0205] Step S453: when the overall coefficient of operation risk is greater than or equal to the overall coefficient threshold of operation risk, the target mining operation area is evaluated as a high water inrush and sand burst risk area;

[0206] Step S454: When the overall operation risk coefficient is less than the overall operation risk coefficient threshold, the target mining operation area is evaluated as a low water inrush and sand burst risk area.

[0207] In this application, if corresponding calculation formulas appear, the above calculation formulas are all dimensionless and take their numerical calculations. The weight coefficients, proportional coefficients and other coefficients in the formulas are set to a result value obtained by quantifying each parameter. The size of the weight coefficient and the proportional coefficient can be determined as long as it does not affect the proportional relationship between the parameter and the result value. Embodiment 2

[0208] See also Figure 2 , based on another concept of the same invention, a system for analyzing factors affecting water inrush and sand burst in mining under high water pressure loose aquifers is proposed, which is applied to a method for analyzing factors affecting water inrush and sand burst in mining under high water pressure loose aquifers. The influencing factor analysis system includes a mining data module, a geological data module, a hydrological data module, a risk assessment module and a server. The mining data module, the geological data module, the hydrological data module and the risk assessment module are respectively connected to the server, and the server controls the mining data module, the geological data module, the hydrological data module and the risk assessment module respectively;

[0209] The mining data module obtains mining stage division data, regional historical operation coefficient, regional real-time operation coefficient and high-density foundation pit area respectively, and obtains regional mining risk coefficient by analyzing regional historical operation coefficient and regional real-time operation coefficient, and defines mining stage division data, regional real-time operation coefficient and high-density foundation pit area as regional mining operation data;

[0210] Acquire regional mining operation engineering data, divide the regional mining operation engineering into b mining operation stages in chronological order according to the regional mining operation engineering data, and name them as the first mining operation stage to the bth mining operation stage, respectively, to obtain mining stage division data;

[0211] It should be noted here that:

[0212] The b designed here is the quantity value corresponding to the down-mining operation stage, and b is an integer greater than 0;

[0213] When the mining operation is in the ath mining operation stage, historical operation data of the target mining operation area is analyzed to obtain the regional historical operation coefficient;

[0214] The details are as follows:

[0215] Obtain the rock stratum mining depth corresponding to the target mining operation area before the start of the ath mining operation stage to obtain the first rock stratum mining depth, obtain the rock stratum target mining depth corresponding to the target mining operation area in the ath mining operation stage to obtain the second rock stratum mining depth, calculate the difference between the second rock stratum mining depth and the first rock stratum mining depth, and take the absolute value of the obtained difference to obtain the stage rock stratum mining depth;

[0216] It should be noted here that:

[0217] The target mining operation area involved here is the area where the analysis of factors affecting water inrush and sand burst in mining under high water pressure loose aquifers is conducted in this application;

[0218] It should be noted here that:

[0219] The rock mining depth involved here is the maximum depth of rock mining;

[0220] Obtain the rock stratum mining area corresponding to the target mining operation area before the start of the ath mining operation stage to obtain the first rock stratum mining area, obtain the rock stratum target mining area corresponding to the target mining operation area in the ath mining operation stage to obtain the second rock stratum mining area, calculate the difference between the first rock stratum mining area and the second rock stratum mining area, and take the absolute value of the obtained difference to obtain the stage rock stratum mining area;

[0221] Obtain the number of mining foundation pits corresponding to the target mining operation area before the start of the ath mining operation stage, and obtain the first foundation pit number value; obtain the target mining foundation pit number corresponding to the target mining operation area in the ath mining operation stage, and obtain the second foundation pit number;

[0222] The foundation pit involved here is specifically a pit structure with a large depth and relatively small area formed for excavating earth and stone during underground engineering construction;

[0223] Calculate the ratio of the first foundation pit quantity value to the first rock stratum mining area to obtain the first foundation pit quantity area ratio, calculate the ratio of the second foundation pit quantity value to the second rock stratum mining area to obtain the second foundation pit quantity area ratio, calculate the difference between the first foundation pit quantity area ratio and the second foundation pit area quantity ratio, and take the absolute value of the obtained difference to obtain the stage foundation pit area quantity ratio;

[0224] The regional historical operation coefficient corresponding to the a-th mining operation stage is obtained by calculating the ratio of the stage rock stratum mining depth, stage rock stratum mining area and stage foundation pit area;

[0225] The regional historical operation coefficient corresponding to the a-th mining operation stage is calculated as follows:

[0226] ;

[0227] Among them, Ygz is the regional historical operation coefficient, Jsd is the mining depth of the stage rock formation, Jmj is the mining area of ​​the stage rock formation, and Mjb is the number ratio of the stage foundation pit area;

[0228] When the mining operation is in the ath mining operation stage, real-time operation data analysis is performed on the target mining operation area to obtain the regional real-time operation coefficient;

[0229] The details are as follows:

[0230] When the target operation area is mined, a plan view of the target operation area is obtained in real time, and characteristic circles are used to fill the foundation pit areas in the plan view of the target operation area to obtain multiple characteristic circles;

[0231] See also Figure 3 , take the lower left corner vertex of the target operation area plane map as the coordinate origin, take the lower border of the target operation area plane map as the coordinate X-axis, and draw a straight line perpendicular to the coordinate X-axis through the coordinate origin as the coordinate Y-axis to obtain the coordinate system of the lower mining operation area;

[0232] Marking each characteristic circle in the coordinate system of the lower mining operation area, and acquiring the center coordinates of each characteristic circle in the coordinate system of the lower mining operation area to obtain a plurality of center coordinates;

[0233] It should be noted here that:

[0234] In this application, the center of each characteristic circle is a set of two-dimensional data points in the coordinate system of the down-mining operation area, and all exist in the form of rectangular coordinates (x, y);

[0235] Randomly select a coordinate from multiple circle center coordinates as the sample circle center coordinate;

[0236] The density value corresponding to the sample circle center coordinates is calculated by using the sample circle center coordinates and multiple circle center coordinates;

[0237] Calculate the density value corresponding to the sample center coordinates. The specific formula is as follows:

[0238] ;

[0239] Among them, Mz is the density value corresponding to the coordinates of the sample center, (x i ,y i ) is the sample center coordinate, K() is the characteristic Gaussian kernel, n is the number value corresponding to the characteristic circle in the coordinate system of the down-sampling operation area, x is the horizontal coordinate of the arbitrary circle center coordinate, y is the vertical coordinate of the arbitrary circle center coordinate, and h is the bandwidth parameter set in the characteristic Gaussian kernel;

[0240] A coordinate is randomly selected from multiple circle center coordinates as the sample circle center coordinate, and the sample circle center coordinate is (x i ,y i );

[0241] The coordinates of the sample circle center are (x i ,y i ) Calculate the density value corresponding to the sample center coordinates;

[0242] Calculate the density value corresponding to the sample center coordinates. The specific formula is as follows:

[0243] ;

[0244] Among them, Mz is the density value corresponding to the coordinates of the sample center, (x i ,y i ) is the coordinate of the sample circle center, K() is the characteristic Gaussian kernel, n is the number value corresponding to the characteristic circle in the coordinate system of the down-mining operation area, x is the horizontal coordinate of the circle center, and y is the vertical coordinate of the circle center;

[0245] It should be noted here that:

[0246] In this application, the characteristic Gaussian kernel is calculated as follows:

[0247] ;

[0248] Among them, K(u) is the characteristic Gaussian kernel, u is the standard deviation of the set Gaussian kernel, and e is a natural constant;

[0249] Repeat the calculation process of the density value corresponding to the sample circle center coordinates, and obtain the density value corresponding to each circle center coordinate respectively;

[0250] Obtain a circle center coordinate density threshold, take the circle center coordinates with a density value greater than the circle center coordinate density threshold as the first type of circle center coordinates, take the circle center coordinates with a density value less than or equal to the circle center coordinate density threshold as the second type of circle center coordinates, and mark the foundation pit area composed of foundation pits corresponding to multiple first type circle center coordinates as a high-density foundation pit area;

[0251] It should be noted here that:

[0252] The circle center coordinate density threshold involved here is the minimum value of the circle center coordinate density of the foundation pit corresponding to the high-density foundation pit area. The circle center coordinate density threshold here can be specifically set according to actual conditions.

[0253] The depth value corresponding to each foundation pit in the high-density foundation pit area is obtained to obtain multiple foundation pit depth values, and the foundation pit depth value with the largest value is marked as the peak foundation pit depth;

[0254] The average of multiple foundation pit depth values ​​is calculated to obtain the average depth of the foundation pit;

[0255] When the a-th mining operation stage is in operation, the foundation pit vibration amplitude values ​​corresponding to each foundation pit in the high-density foundation pit area are obtained respectively, and multiple vibration amplitude values ​​are obtained, and the average value is calculated to obtain the average amplitude of the foundation pit area;

[0256] The regional real-time operation coefficient is obtained by calculating the data of the mining stage division, the peak foundation pit depth, the average foundation pit depth and the average amplitude of the foundation pit area;

[0257] The regional real-time operation coefficient is calculated, and the specific formula is as follows:

[0258] ;

[0259] Among them, Ssx is the regional real-time operation coefficient, Fjs is the peak foundation pit depth, Pjs is the average foundation pit depth, and Kzf is the average amplitude of the foundation pit area;

[0260] The regional historical operation coefficient and the regional real-time operation coefficient are summed to obtain the regional mining risk coefficient;

[0261] The regional mining risk factor and high-density foundation pit area are defined as regional mining operation data;

[0262] The data acquisition module acquires the regional mining operation data and transmits it to the geological data module, hydrological data module, and risk assessment module;

[0263] The hydrological data module performs regional water content analysis on the target operation area based on the regional down-sampling operation data to obtain regional hydrological analysis data;

[0264] Obtain regional mining operation data, and obtain high-density foundation pit area and mining stage division data based on the regional mining operation data;

[0265] When the mining operation is in the ath mining operation stage, the regional water content coefficient corresponding to the high-density foundation pit area is obtained;

[0266] The details are as follows:

[0267] A foundation pit is randomly selected in the high-density foundation pit area as a water-containing monitoring sample foundation pit;

[0268] Drill from the bottom of the water-bearing monitoring sample foundation pit down to the target drilling depth. If no groundwater is found, the water-bearing monitoring sample foundation pit is marked as a first-type foundation pit;

[0269] Drill downward from the bottom of the water-containing monitoring sample foundation pit to the target drilling depth. If groundwater is explored, the water-containing monitoring sample foundation pit is marked as the second type of foundation pit. The drilling depth when groundwater and soil suddenly surge is obtained to obtain the groundwater depth value corresponding to the water-containing monitoring sample foundation pit. The surge pressure value when groundwater suddenly surges is obtained to obtain the groundwater surge pressure value corresponding to the water-containing monitoring sample foundation pit.

[0270] It should be noted here that:

[0271] In this application, a sudden surge is the first appearance of groundwater after drilling;

[0272] This involves safe drilling where the target drilling depth is groundwater, i.e. there is no risk of water inrush below this depth;

[0273] The number of second type foundation pits is obtained to obtain the number value of second type foundation pits, the area value of high-density foundation pit areas is obtained to obtain the area value of high-density foundation pits, and the ratio of the number value of second type foundation pits to the area value of high-density foundation pits is calculated to obtain the number-area ratio of characteristic foundation pits;

[0274] Perform water content monitoring on each foundation pit in the high-density foundation pit area to obtain multiple second-type foundation pits, obtain groundwater depth values ​​corresponding to each second-type foundation pit to obtain multiple groundwater depth values, and average the multiple groundwater depth values ​​to obtain an average sudden surge depth value of the foundation pit;

[0275] The groundwater surge pressure values ​​corresponding to each second type foundation pit are obtained respectively to obtain a plurality of groundwater surge pressure values, and the plurality of groundwater surge pressure values ​​are averaged to obtain an average foundation pit surge pressure value;

[0276] The regional water content coefficient corresponding to the a-th mining operation stage is obtained by calculating the characteristic foundation pit number area ratio, the average foundation pit sudden gushing depth value and the average foundation pit sudden gushing pressure value;

[0277] The regional moisture coefficient is calculated using the following formula:

[0278] ;

[0279] Among them, Qhsa is the regional water content coefficient corresponding to the ath mining operation stage, Tsz is the average sudden surge depth value of the foundation pit, Tyz is the average sudden surge pressure value of the foundation pit, and Tmb is the ratio of the number and area of ​​characteristic foundation pits;

[0280] When the ath mining operation stage is the first mining operation stage, the regional water content coefficient corresponding to the ath mining operation stage is defined as the regional hydrological risk coefficient;

[0281] When the ath mining operation stage is not the first mining operation stage, the regional water content coefficient corresponding to the a-1th mining operation stage is obtained, and the regional water content coefficient corresponding to the ath mining operation stage and the regional water content coefficient corresponding to the a-1th mining operation stage are calculated to obtain the regional hydrological risk coefficient;

[0282] The regional hydrological risk coefficient is calculated using the following formula:

[0283] ;

[0284] Among them, Qfx is the regional hydrological risk coefficient, Qhsa is the regional water content coefficient corresponding to the a-th mining operation stage, and Qhs (a-1) is the regional water content coefficient corresponding to the a-1-th mining operation stage;

[0285] The regional hydrological risk coefficient and the ratio of the number and area of ​​characteristic foundation pits are defined as regional hydrological analysis data;

[0286] The hydrological data module acquires regional hydrological analysis data and transmits it to the risk assessment module and the geological data module;

[0287] The geological data module conducts regional geological analysis on the target operation area based on the regional mining operation data to obtain the regional geological risk coefficient;

[0288] Obtain regional mining operation data, and obtain high-density foundation pit areas based on the regional mining operation data;

[0289] Obtain regional hydrological analysis data, and obtain the ratio of the number and area of ​​characteristic foundation pits based on the regional hydrological analysis data;

[0290] A foundation pit is randomly selected in the high-density foundation pit area as a geological monitoring sample pit;

[0291] Conduct geological monitoring on the geological monitoring sample foundation pit to obtain the foundation pit geological analysis coefficient corresponding to the geological monitoring sample foundation pit;

[0292] The details are as follows:

[0293] A rock mass of unit mass is selected from a geological monitoring sample foundation pit as a sample monitoring rock mass, and the sample monitoring rock mass is placed in a water body so that the sample monitoring rock mass is in a completely submerged state. When the sample monitoring rock mass is submerged to a characteristic submerged time, the sample monitoring rock mass is taken out to obtain a wet sample monitoring rock mass, and the mass of the wet sample monitoring rock mass is obtained to obtain a wet sample rock mass mass value. The volume value corresponding to the wet sample monitoring rock mass is obtained by a water displacement method to obtain a sample monitoring rock mass volume value;

[0294] It should be noted here that:

[0295] In this application, the unit mass is specifically 5 kg here, and the characteristic immersion time is specifically 30 min;

[0296] The sample monitoring rock mass is placed in a characteristic temperature environment, and the mass value of the sample monitoring rock mass is obtained in real time until the mass value of the sample monitoring rock mass maintains a constant weight, thereby obtaining a dry sample monitoring rock mass, and the mass of the dry sample monitoring rock mass is obtained to obtain the mass value of the dry sample rock mass;

[0297] It should be noted here that:

[0298] In this application, the characteristic temperature environment is specifically 101 degrees Celsius;

[0299] The rock mass porosity corresponding to the sample monitoring rock mass is obtained by calculating the dry sample rock mass mass value, the sample monitoring rock mass volume value and the wet sample rock mass value;

[0300] The rock mass porosity corresponding to the sample monitoring rock mass is calculated, and the specific formula is as follows:

[0301] ;

[0302] Among them, Kxd is the rock mass porosity corresponding to the sample monitoring rock mass, Ms is the wet sample rock mass value, Mg is the dry sample rock mass value, Vy is the sample monitoring rock mass volume value, and ρ is the water density value;

[0303] In the geological monitoring sample foundation pit, multiple monitoring rock masses are randomly selected, and the rock mass porosity corresponding to each monitoring rock mass is obtained respectively to obtain multiple rock mass porosities, and the average of the obtained multiple rock mass porosities is calculated to obtain the average rock mass porosity of the foundation pit;

[0304] The air humidity of the geological monitoring sample foundation pit is obtained to obtain the foundation pit air humidity value, and the depth value of the geological monitoring sample foundation pit is obtained to obtain the foundation pit depth value;

[0305] It should be noted here that:

[0306] In this application, the depth of the foundation pit and the distance between the bottom of the foundation pit and the ground are involved;

[0307] The average rock mass porosity, depth and air humidity of the foundation pit are calculated to obtain the foundation pit geological analysis coefficient corresponding to the geological monitoring sample foundation pit;

[0308] The geological analysis coefficient of the foundation pit is calculated, and the specific formula is configured as follows:

[0309] ;

[0310] Among them, Jkx is the geological analysis coefficient of the foundation pit, Sdz is the foundation pit depth value, Wdz is the foundation pit air humidity value, and Kxp is the average rock porosity of the foundation pit;

[0311] In the high-density foundation pit area, the foundation pit geological analysis coefficient corresponding to each foundation pit is obtained to obtain multiple foundation pit geological analysis coefficients, and the obtained foundation pit geological analysis coefficients are averaged to obtain the regional geological risk coefficient;

[0312] The geological data module acquires the regional geological risk coefficient and transmits it to the risk assessment module;

[0313] The risk assessment module conducts an operational risk assessment on the target mining operation area based on the regional mining operation data, regional hydrological analysis data and regional geological risk factor;

[0314] Obtain regional mining operation data, and obtain regional mining risk coefficient based on the regional mining operation data;

[0315] Obtain regional hydrological analysis data, and obtain regional hydrological risk coefficients based on the regional hydrological analysis data;

[0316] Obtain regional geological risk factors;

[0317] The regional mining risk coefficient, regional hydrological risk coefficient and regional geological risk coefficient are calculated to obtain the overall operation risk coefficient;

[0318] The overall coefficient of operation risk is calculated as follows:

[0319] ;

[0320] Among them, Fxz is the overall coefficient of operation risk, Qcz is the regional mining risk coefficient, Qfx is the regional hydrological risk coefficient, and Qdz is the regional geological risk coefficient;

[0321] Obtaining the threshold value of the overall coefficient of operation risk, comparing the overall coefficient of operation risk with the threshold value of the overall coefficient of operation risk, and conducting an operation risk assessment on the target mining operation area based on the numerical comparison result;

[0322] The details are as follows:

[0323] Obtain the regional mining risk factor threshold, regional hydrological risk factor threshold and regional geological risk factor threshold respectively;

[0324] It should be noted here that:

[0325] In the present application, the regional mining risk coefficient threshold involved here is the minimum regional mining risk coefficient corresponding to the area where the risk of water inrush and sand burst is known to exist, the regional hydrological risk coefficient threshold involved here is the regional hydrological risk coefficient corresponding to the area where the risk of water inrush and sand burst is known to exist, and the regional geological risk coefficient threshold involved here is the minimum regional geological risk coefficient corresponding to the area where the risk of water inrush and sand burst is known to exist. In reality, the mining risk coefficient threshold, the regional hydrological risk coefficient threshold and the regional geological risk coefficient threshold can be specifically set according to the actual situation;

[0326] The regional mining risk coefficient threshold, the regional hydrological risk coefficient threshold and the regional geological risk coefficient threshold are calculated to obtain the overall operation risk coefficient threshold;

[0327] The overall coefficient threshold of the operation risk is calculated, and the specific formula is as follows:

[0328] ;

[0329] Among them, Fxzy is the overall coefficient threshold of operation risk, Qczy is the regional mining risk coefficient threshold, Qfxy is the regional hydrological risk coefficient threshold, and Qdzy is the regional geological risk coefficient threshold;

[0330] When the overall coefficient of operation risk is greater than or equal to the overall coefficient threshold of operation risk, the target mining operation area is assessed as a high water inrush and sand burst risk area;

[0331] When the overall operation risk coefficient is less than the overall operation risk coefficient threshold, the target mining operation area is assessed as a low water inrush and sand burst risk area.

[0332] The preferred embodiments of the present invention disclosed above are only used to help explain the present invention. The preferred embodiments do not describe all the details in detail, nor do they limit the invention to only specific implementation methods. Obviously, many modifications and changes can be made according to the content of this specification. This specification selects and specifically describes these embodiments in order to better explain the principles and practical applications of the present invention, so that those skilled in the art can understand and use the present invention well. The present invention is limited only by the claims and their full scope and equivalents.

Claims

1. The method for analyzing factors affecting water inrush and sand burst in high-pressure loose aquifers is characterized by: include: Step S1: Acquire the mining stage division data, regional historical operation coefficient, regional real-time operation coefficient and high-density foundation pit area, analyze the regional historical operation coefficient and regional real-time operation coefficient to obtain the regional mining risk coefficient, and define the regional mining risk coefficient and the high-density foundation pit area as regional mining operation data; Step S2: Analyze the regional mining operation data to divide the foundation pits in the high-density foundation pit area into the first type of foundation pits and the second type of foundation pits, perform hydrological analysis on the second type of foundation pits, obtain characteristic foundation pit quantity area ratio, average foundation pit sudden gushing depth value and average foundation pit sudden gushing pressure value, analyze the characteristic foundation pit quantity area ratio, average foundation pit sudden gushing depth value and average foundation pit sudden gushing pressure value to obtain the regional hydrological risk coefficient, and define the characteristic foundation pit quantity area ratio and the regional hydrological risk coefficient as regional hydrological analysis data; Step S3: Analyze the regional mining operation data and regional hydrological analysis data to obtain the average rock porosity, pit depth value and pit air humidity value corresponding to the high-density pit area, and calculate the regional geological risk coefficient by using the average rock porosity, pit depth value and pit air humidity value; Step S4: Calculate the regional mining risk coefficient, the regional hydrological risk coefficient and the regional geological risk coefficient to obtain the overall operation risk coefficient, and perform operation risk assessment by comparing the overall operation risk coefficient with the threshold.

2. The method for analyzing factors affecting water inrush and sand burst in mining under high water pressure loose aquifers according to claim 1 is characterized in that: The step S1 further includes the following specific steps: Step S11: Acquire regional mining operation engineering data, divide the regional mining operation engineering into b mining operation stages in chronological order according to the regional mining operation engineering data, and name them as the first mining operation stage to the bth mining operation stage, respectively, to obtain mining stage division data; Step S12: Analyze the historical operation data of the target mining operation area to obtain the historical operation coefficient of the area; Step S13: Analyze the real-time operation data of the target mining operation area to obtain the real-time operation coefficient of the area; Step S14: summing the regional historical operation coefficient and the regional real-time operation coefficient to obtain the regional mining risk coefficient; Step S15: defining the regional mining risk coefficient and the high-density foundation pit area as regional mining operation data.

3. The method for analyzing factors affecting water inrush and sand burst in mining under high water pressure loose aquifers according to claim 2 is characterized in that: The step S12 further includes the following specific steps: step S121: Obtain the rock stratum mining depth corresponding to the target mining operation area before the start of the ath mining operation stage to obtain the first rock stratum mining depth, obtain the rock stratum target mining depth corresponding to the target mining operation area in the ath mining operation stage to obtain the second rock stratum mining depth, calculate the difference between the second rock stratum mining depth and the first rock stratum mining depth, and take the absolute value of the obtained difference to obtain the stage rock stratum mining depth; Step S122: obtaining the rock stratum mining area corresponding to the target mining operation area before the start of the ath mining operation stage to obtain the first rock stratum mining area, obtaining the rock stratum target mining area corresponding to the target mining operation area in the ath mining operation stage to obtain the second rock stratum mining area, calculating the difference between the first rock stratum mining area and the second rock stratum mining area, and taking the absolute value of the obtained difference to obtain the stage rock stratum mining area; Step S123: obtaining the number of mining foundation pits corresponding to the target mining operation area before the start of the ath mining operation stage, obtaining the first foundation pit number value, obtaining the target mining foundation pit number corresponding to the target mining operation area in the ath mining operation stage, obtaining the second foundation pit number; Step S124: Calculate the ratio of the first foundation pit quantity value to the first rock stratum mining area to obtain the first foundation pit quantity area ratio, calculate the ratio of the second foundation pit quantity value to the second rock stratum mining area to obtain the second foundation pit quantity area ratio, calculate the difference between the first foundation pit quantity area ratio and the second foundation pit area quantity ratio, and take the absolute value of the obtained difference to obtain the stage foundation pit area quantity ratio; Step S125: Calculate the regional historical operation coefficient corresponding to the a-th mining operation stage by using the ratio of the stage rock stratum mining depth, the stage rock stratum mining area and the stage foundation pit area; The regional historical operation coefficient corresponding to the a-th mining operation stage is calculated as follows: ； Among them, Ygz is the regional historical operation coefficient, Jsd is the stage rock stratum mining depth, Jmj is the stage rock stratum mining area, and Mjb is the stage foundation pit area quantity ratio.

4. The method for analyzing factors affecting water inrush and sand burst in mining under high water pressure loose aquifers according to claim 2 is characterized in that: The step S13 further includes the following specific steps: Step S131: when performing mining operations on the target operation area, a plan view of the target operation area is obtained in real time, and characteristic circles are used to fill the foundation pit areas in the plan view of the target operation area to obtain a plurality of characteristic circles; Step S132: Acquire the high-density foundation pit area; The step S132 further includes the following specific steps: Step S1321: Take the lower left corner vertex of the target operation area plan as the coordinate origin, take the lower frame of the target operation area plan as the coordinate X axis, and draw a straight line perpendicular to the coordinate X axis through the coordinate origin as the coordinate Y axis to obtain the coordinate system of the lower mining operation area; Step S1322: Mark each characteristic circle in the lower mining operation area coordinate system, and obtain the center coordinates of each characteristic circle in the lower mining operation area coordinate system to obtain a plurality of center coordinates; Step S1323: randomly selecting a coordinate from a plurality of circle center coordinates as a sample circle center coordinate; Step S1324: Calculate the density value corresponding to the sample circle center coordinates through the sample circle center coordinates and multiple circle center coordinates; Calculate the density value corresponding to the sample center coordinates. The specific formula is as follows: ； Among them, Mz is the density value corresponding to the coordinates of the sample center, (x i ,y i ) is the sample center coordinate, K() is the characteristic Gaussian kernel, n is the number value corresponding to the characteristic circle in the coordinate system of the down-sampling operation area, x is the horizontal coordinate of the arbitrary circle center coordinate, y is the vertical coordinate of the arbitrary circle center coordinate, and h is the bandwidth parameter set in the characteristic Gaussian kernel; The characteristic Gaussian kernel is calculated as follows: ； Among them, K(u) is the characteristic Gaussian kernel, u is the standard deviation of the set Gaussian kernel, and e is a natural constant; Step S1325: Obtain the density value corresponding to each circle center coordinate respectively; Step S1326: Obtain a circle center coordinate density threshold, take the circle center coordinates with a density value greater than the circle center coordinate density threshold as the first type of circle center coordinates, take the circle center coordinates with a density value less than or equal to the circle center coordinate density threshold as the second type of circle center coordinates, and mark the foundation pit area composed of the foundation pits corresponding to multiple first type of circle center coordinates as a high-density foundation pit area; Step S133: acquiring a depth value corresponding to each foundation pit in the high-density foundation pit area to obtain a plurality of foundation pit depth values, and marking the foundation pit depth value with the largest value as the peak foundation pit depth; Step S134: Calculate the average of multiple foundation pit depth values ​​to obtain an average foundation pit depth; Step S135: when the a-th mining operation stage is in operation, respectively obtain the foundation pit vibration amplitude value corresponding to each foundation pit in the high-density foundation pit area, obtain multiple vibration amplitude values, and calculate the average of them to obtain the average amplitude of the foundation pit area; Step S136: Calculate the peak foundation pit depth, the average foundation pit depth, and the average amplitude of the foundation pit area to obtain the regional real-time operation coefficient; The regional real-time operation coefficient is calculated, and the specific formula is as follows: ； Among them, Ssx is the regional real-time operation coefficient, Fjs is the peak foundation pit depth, Pjs is the average foundation pit depth, and Kzf is the average amplitude of the foundation pit area.

5. The method for analyzing factors affecting water inrush and sand burst in mining under high water pressure loose aquifers according to claim 1 is characterized in that: The step S2 further includes the following specific steps: Step S21: Acquire regional mining operation data, and acquire high-density foundation pit area and mining stage division data according to the regional mining operation data; Step S22: when the mining operation is in the ath mining operation stage, obtaining the regional water content coefficient and the characteristic foundation pit quantity area ratio corresponding to the high-density foundation pit area; Step S23: Obtaining the regional hydrological risk coefficient; Step S24: define the regional hydrological risk coefficient and the ratio of the number and area of ​​characteristic foundation pits as regional hydrological analysis data.

6. The method for analyzing factors affecting water inrush and sand burst in mining under high water pressure loose aquifers according to claim 5 is characterized in that: The step S22 further includes the following specific steps: Step S221: randomly selecting a foundation pit in the high-density foundation pit area as a water-containing monitoring sample foundation pit; Step S222: Drilling downward from the bottom of the water-containing monitoring sample foundation pit to a target drilling depth, and if no groundwater is found, marking the water-containing monitoring sample foundation pit as a first type foundation pit; Step S223: Drilling downward from the bottom of the water-containing monitoring sample foundation pit to the target drilling depth, if groundwater is explored, the water-containing monitoring sample foundation pit is marked as a second type of foundation pit, the drilling depth when groundwater and soil surge occurs is obtained, and the groundwater depth value corresponding to the water-containing monitoring sample foundation pit is obtained, and the surge pressure value when groundwater surges occurs is obtained, and the groundwater surge pressure value corresponding to the water-containing monitoring sample foundation pit is obtained; Step S224: acquiring the number of second-type foundation pits to obtain the number value of second-type foundation pits, acquiring the area value of the high-density foundation pit area to obtain the high-density foundation pit area value, and calculating the ratio of the second-type foundation pit number value to the high-density foundation pit area value to obtain the characteristic foundation pit number area ratio; Step S225: performing water content monitoring on each foundation pit in the high-density foundation pit area to obtain a plurality of second-type foundation pits, obtaining a groundwater depth value corresponding to each second-type foundation pit to obtain a plurality of groundwater depth values, and averaging the plurality of groundwater depth values ​​to obtain an average sudden inrush depth value of the foundation pit; Step S226: acquiring the groundwater surge pressure value corresponding to each second type foundation pit respectively to obtain a plurality of groundwater surge pressure values, and averaging the plurality of groundwater surge pressure values ​​to obtain an average foundation pit surge pressure value; Step S227: Calculate the area ratio of the characteristic foundation pit quantity, the average sudden gushing depth value of the foundation pit, and the average sudden gushing pressure value of the foundation pit to obtain the regional water content coefficient corresponding to the a-th mining operation stage; The regional moisture coefficient is calculated using the following formula: ； Among them, Qhsa is the regional moisture coefficient corresponding to the ath mining operation stage, Tsz is the average sudden surge depth of the foundation pit, Tyz is the average sudden surge pressure of the foundation pit, and Tmb is the ratio of the number and area of ​​characteristic foundation pits.

7. The method for analyzing factors affecting water inrush and sand burst in mining under high water pressure loose aquifers according to claim 5 is characterized in that: The step S23 further includes the following specific steps: Step S231: when the ath mining operation stage is the first mining operation stage, the regional water content coefficient corresponding to the ath mining operation stage is defined as the regional hydrological risk coefficient; Step S232: when the ath mining operation stage is not the first mining operation stage, the regional water content coefficient corresponding to the a-1th mining operation stage is obtained, and the regional water content coefficient corresponding to the ath mining operation stage and the regional water content coefficient corresponding to the a-1th mining operation stage are calculated to obtain the regional hydrological risk coefficient; The regional hydrological risk coefficient is calculated using the following formula: ； Among them, Qfx is the regional hydrological risk coefficient, Qhsa is the regional moisture coefficient corresponding to the a-th mining operation stage, and Qhs (a-1) is the regional moisture coefficient corresponding to the a-1-th mining operation stage.

8. The method for analyzing factors affecting water inrush and sand burst in mining under high water pressure loose aquifers according to claim 1, characterized in that: The step S3 further includes the following specific steps: Step S31: Acquire regional mining operation data, and acquire a high-density foundation pit area according to the regional mining operation data; Step S32: obtaining regional hydrological analysis data, and obtaining characteristic foundation pit quantity and area ratio according to the regional hydrological analysis data; Step S33: randomly selecting a foundation pit in the high-density foundation pit area as a geological monitoring sample foundation pit; Step S34: Perform geological monitoring on the geological monitoring sample foundation pit to obtain the foundation pit geological analysis coefficient corresponding to the geological monitoring sample foundation pit; Step S35: In the high-density foundation pit area, the foundation pit geological analysis coefficient corresponding to each foundation pit is obtained to obtain multiple foundation pit geological analysis coefficients, and the obtained foundation pit geological analysis coefficients are averaged to obtain the regional geological risk coefficient.

9. The method for analyzing factors affecting water inrush and sand burst in mining under high water pressure loose aquifers according to claim 8, characterized in that: The step S34 further includes the following specific steps: Step S341: select a rock mass of unit mass in the geological monitoring sample foundation pit as the sample monitoring rock mass, place the sample monitoring rock mass in the water body, make the sample monitoring rock mass in a completely submerged state, wait for the submerged state to reach the characteristic submerged time, take out the sample monitoring rock mass to obtain the wet sample monitoring rock mass, obtain the mass value of the wet sample monitoring rock mass, obtain the volume value corresponding to the wet sample monitoring rock mass by the water displacement method, and obtain the volume value of the sample monitoring rock mass; Step S342: placing the sample monitoring rock mass in a characteristic temperature environment, and obtaining the mass value of the sample monitoring rock mass in real time until the mass value of the sample monitoring rock mass maintains a constant weight, obtaining a dry sample monitoring rock mass, and obtaining the mass of the dry sample monitoring rock mass to obtain the mass value of the dry sample rock mass; Step S343: The dry sample rock mass value, the sample monitoring rock mass volume value and the wet sample rock mass value are calculated to obtain the rock mass porosity corresponding to the sample monitoring rock mass; The rock mass porosity corresponding to the sample monitoring rock mass is calculated, and the specific formula is as follows: ； Among them, Kxd is the rock mass porosity corresponding to the sample monitoring rock mass, Ms is the wet sample rock mass value, Mg is the dry sample rock mass value, Vy is the sample monitoring rock mass volume value, and ρ is the water density value; Step S344: randomly selecting multiple monitoring rock masses in the geological monitoring sample foundation pit, respectively acquiring the rock mass porosity corresponding to each monitoring rock mass, obtaining multiple rock mass porosities, and averaging the obtained multiple rock mass porosities to obtain the average rock mass porosity of the foundation pit; Step S345: Acquire the air humidity of the geological monitoring sample foundation pit to obtain the foundation pit air humidity value, and acquire the depth value of the geological monitoring sample foundation pit to obtain the foundation pit depth value; Step S346: Calculate the average rock porosity, the depth and the air humidity of the foundation pit to obtain the regional geological risk coefficient corresponding to the geological monitoring sample foundation pit; The geological analysis coefficient of the foundation pit is calculated, and the specific formula is configured as follows: ； Among them, Jkx is the geological analysis coefficient of the foundation pit, Sdz is the foundation pit depth value, Wdz is the foundation pit air humidity value, and Kxp is the average rock porosity of the foundation pit.

10. The method for analyzing factors affecting water inrush and sand burst in mining under high water pressure loose aquifers according to claim 1, characterized in that: The step S4 further includes the following specific steps: Step S41: Acquire regional mining operation data, and acquire regional mining risk coefficient according to the regional mining operation data; Step S42: obtaining regional hydrological analysis data, and obtaining a regional hydrological risk coefficient according to the regional hydrological analysis data; Step S43: Obtaining regional geological risk coefficient; Step S44: Calculate the regional mining risk coefficient, the regional hydrological risk coefficient and the regional geological risk coefficient to obtain an overall operation risk coefficient; The overall coefficient of operation risk is calculated as follows: ； Among them, Fxz is the overall coefficient of operation risk, Qcz is the regional mining risk coefficient, Qfx is the regional hydrological risk coefficient, and Qdz is the regional geological risk coefficient; Step S45: obtaining an overall operation risk coefficient threshold, comparing the overall operation risk coefficient with the overall operation risk coefficient threshold, and conducting an operation risk assessment on the target mining operation area according to the numerical comparison result; The step S45 further includes the following specific steps: Step S451: respectively obtaining a regional mining risk coefficient threshold, a regional hydrological risk coefficient threshold, and a regional geological risk coefficient threshold; Step S452: Calculate the regional mining risk coefficient threshold, the regional hydrological risk coefficient threshold, and the regional geological risk coefficient threshold to obtain an overall operation risk coefficient threshold; Step S453: when the overall coefficient of operation risk is greater than or equal to the overall coefficient threshold of operation risk, the target mining operation area is evaluated as a high water inrush and sand burst risk area; Step S454: When the overall operation risk coefficient is less than the overall operation risk coefficient threshold, the target mining operation area is evaluated as a low water inrush and sand burst risk area.

Images