Method and system for predicting underground water level of karst mountainous area

Abstract

The invention discloses a karst mountainous area underground water level prediction method and system, and relates to the technical field of underground water level prediction.An underground water level prediction model of a karst mountainous area is constructed, underground excavation simulation is carried out, the prediction efficiency of the prediction model is analyzed, when the prediction efficiency is poor, subsequent optimization is carried out, and when the prediction effect is good, the prediction efficiency is improved. According to the simulation result, the excavation position where the underground water level changes greatly in the excavation process is determined, possible dangerous phenomena are predicted, then a prevention scheme is made, and reference is provided for protection of follow-up actual excavation. The underground water level change can be rapidly and accurately obtained in subsequent actual construction, protection is carried out in time, accidents such as water inrush and mud inrush are reduced, the life safety of constructors and equipment safety are guaranteed, the phenomenon that construction is interrupted or extra treatment measures need to be taken is reduced, and therefore the project progress is guaranteed, and the project cost is reduced.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of groundwater level prediction, in particular to a karst mountainous area groundwater level prediction method and system. BACKGROUND

[0002] Karst mountainous areas have developed underground caves and fissures, and are rich in groundwater with complex runoff. Excavation in karst mountainous areas may damage the original groundwater runoff channels and aquifers, causing groundwater to flow into the construction area and triggering gushing and water inrush accidents. In addition, changes in groundwater level may have long-term adverse effects on the structure, such as causing structure leakage and corrosion. Predicting the groundwater level during excavation can help take effective maintenance measures during the engineering operation phase to ensure the long-term safe use of the project.

[0003] Prior art such as the invention patent with publication number CN119322987A discloses a groundwater level prediction method and system based on spatio-temporal attention gate neural network, which includes: constructing a historical groundwater level spatio-temporal sequence dataset; performing principal component analysis on the spatio-temporal sequence dataset to obtain a reduced dimension spatio-temporal sequence dataset; for the reduced dimension spatio-temporal sequence dataset, a plurality of spatio-temporal feature matrices are constructed by a sliding window method; a trained groundwater level prediction model is obtained by training the groundwater level prediction model using the plurality of spatio-temporal feature matrices; the spatio-temporal sequence data of the historical groundwater level at the to-be-measured location is obtained, the spatio-temporal feature matrix of the groundwater level is constructed, and the trained groundwater level prediction model is used to predict the groundwater level at the to-be-measured location; the present application realizes fast and accurate prediction of the groundwater level, reduces the modeling cost, and greatly improves the generalization ability, prediction accuracy and interpretability of the groundwater level prediction model.

[0004] Prior art such as the invention patent with publication number CN114881323A discloses a deep neural network-based groundwater level prediction and updating method for foundation pit dewatering area, which is used for predicting and updating the 2-D spatial distribution of groundwater level in the foundation pit dewatering area. The method includes: collecting groundwater level observation data of each observation point, constructing a deep neural network model, and establishing a functional relationship through a deep neural network algorithm to predict and update the groundwater level of the target area, providing information for the development of foundation pit construction work. Based on the observation data, a variety of loss functions are fused to solve the problems of single spatial position of prediction points, low accuracy of observation data, large data result dispersion, poor real-time updating effect, and inability to timely and fully reflect the actual situation of water level changes in the past neural network model for groundwater level prediction, which leads to low construction efficiency and inability to obtain accurate data.

[0005] The construction and optimization of the underground water level prediction model disclosed in detail in the above scheme, but the karst mountainous area has a complex geological structure, and the prediction model involves the fusion of multiple data during prediction, thereby affecting the prediction efficiency of the underground water level, but the speed of the underground water level prediction affects the excavation engineering, however, the above scheme does not simulate test the prediction efficiency of the prediction model after the construction of the prediction model is completed, the prediction speed of the prediction model cannot be mastered, the underground water level change cannot be quickly and accurately obtained in subsequent actual construction, thereby the protection cannot be timely performed, the water inrush, mud inrush and other accidents cannot be reduced, and the life safety and equipment safety of construction personnel cannot be ensured, in addition, if the underground water level prediction is not timely, the construction may be interrupted or additional treatment measures need to be taken due to the water level change, thereby the engineering progress is delayed and the engineering cost is increased. SUMMARY

[0006] In view of the above technical deficiencies, the purpose of the present application is to provide a karst mountainous area underground water level prediction method and system.

[0007] To solve the above technical problems, the present application adopts the following technical scheme: in a first aspect, the present application provides a karst mountainous area underground water level prediction method, S1, a water level prediction model is constructed, a three-dimensional model of a specified karst mountainous area and an underground excavation plan of the specified karst mountainous area are obtained, underground excavation simulation is performed on the specified karst mountainous area, and water level prediction is performed on the underground excavation simulation by using the water level prediction model; S2, based on the water level prediction of the underground excavation simulation, the water level prediction efficiency of the water level prediction model is analyzed, when the water level prediction efficiency is poor, S4 is executed, and when the water level prediction efficiency is good, S3 is executed; S3, the water level prediction of the underground excavation simulation is performed by using the water level prediction model, the preventive measures for the underground excavation of the specified karst mountainous area are formulated, and then S4 is executed; S4, when the water level prediction efficiency is poor, the feedback is fed back to the visual display terminal; the preventive measures for the underground excavation of the specified karst mountainous area are fed back to the visual display terminal, and a prevention prompt is performed.

[0008] In a second aspect, the present application provides a karst mountainous area underground water level prediction system, comprising: a model construction and prediction module for constructing a water level prediction model, obtaining a three-dimensional model of a specified karst mountainous area and an underground excavation plan of the specified karst mountainous area, performing underground excavation simulation on the specified karst mountainous area, and performing water level prediction on the underground excavation simulation by using the water level prediction model; a model prediction evaluation module for analyzing the water level prediction efficiency of the water level prediction model based on the water level prediction of the underground excavation simulation, executing a prediction execution module when the water level prediction efficiency is poor, and executing a prediction result processing module when the water level prediction efficiency is good; The prediction result processing module is used to predict the water level of the underground excavation simulation using the water level prediction model, formulate preventive measures for underground excavation in designated karst mountainous areas, and then execute the prediction execution module. The prediction execution module is used to provide feedback to the visualization display terminal when the water level prediction efficiency is poor; it also sends the specified preventive measures for underground excavation in karst mountainous areas to the visualization display terminal and provides preventive prompts.

[0009] The beneficial effects of this invention are as follows: This application provides a method and system for predicting groundwater levels in karst mountainous areas. A groundwater level prediction model for karst mountainous areas is constructed, and underground excavation simulation is performed. The prediction efficiency of the model is analyzed. When the prediction efficiency is poor, subsequent optimization is carried out. When the prediction effect is good, based on the simulation results, the excavation locations with significant groundwater level changes during the excavation process are identified, and potential dangerous phenomena are predicted. Then, prevention plans are formulated to provide a reference for protection during subsequent actual excavation. By evaluating the prediction efficiency of the model, this application ensures that groundwater level changes can be quickly and accurately obtained during subsequent actual construction, enabling timely protection, reducing the occurrence of accidents such as water inrush and mudslide, protecting the lives and equipment of construction personnel, and reducing the occurrence of construction interruptions or the need for additional treatment measures. This ensures project progress and reduces project costs. Attached Figure Description

[0010] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0011] Figure 1 This is a schematic diagram of the implementation steps of the method of the present invention.

[0012] Figure 2 This is a schematic diagram of the system structure connection of the present invention. Detailed Implementation

[0013] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention. Example 1:

[0014] See Figure 1As shown, a groundwater level prediction method for karst mountain areas includes the following steps: S1, constructing a water level prediction model, simultaneously obtaining a three-dimensional model of a specified karst mountain area and an underground excavation plan for the specified karst mountain area, simulating underground excavation in the specified karst mountain area, and using the water level prediction model to predict the water level in the underground excavation simulation. In a specific embodiment, the specific process of constructing the water level prediction model is as follows: S11, through geological exploration, obtain geological data and hydrogeological data of the designated karst mountain area, construct a conceptual model of the designated karst mountain area, and use numerical simulation software to divide the designated karst mountain area into grids and establish a three-dimensional groundwater flow numerical model. Numerical simulation software such as Visual MODFLOW.

[0015] In the above, geological exploration includes drilling, geophysical exploration and other methods, geological data includes stratigraphic distribution, rock type, karst development characteristics, etc., and hydrogeological data includes water flow path and flow velocity changes, etc.

[0016] A conceptual model is an abstract and simplified qualitative description of a groundwater system in karst mountainous areas. It is mainly used to help understand and analyze the operating mechanism of the groundwater system and to provide a basic framework for subsequent model building. Specifically, it includes boundary determination, aquifer characteristic description, water flow process analysis, generalization of hydrogeological conditions, and identification of major influencing factors.

[0017] Boundary Determination: Clearly define the scope of the groundwater system in karst mountainous areas, including its horizontal and vertical boundaries. For example, determine the boundaries of the groundwater recharge zone, runoff zone, and discharge zone, as well as the top and bottom interfaces of the aquifer, based on geological structures, topography, and other factors.

[0018] Aquifer characteristics description: This involves classifying and describing the characteristics of aquifers in karst mountainous areas, including their lithology, permeability, and porosity. For example, it distinguishes between karst-developed carbonate aquifers and relatively impermeable non-carbonate strata, and describes the unique water storage and conduction space characteristics of carbonate aquifers, such as dissolution pores, fissures, and caves formed by karst processes.

[0019] Water flow process analysis: This involves analyzing the flow process of groundwater in karst mountainous areas, including recharge, runoff, and discharge. For example, it considers the recharge of groundwater by atmospheric precipitation through surface infiltration and karst conduits, the runoff of groundwater along a certain hydraulic gradient in the aquifer, and the discharge through springs, rivers, and other pathways.

[0020] Hydrogeological condition generalization: This involves reasonably simplifying and generalizing complex actual hydrogeological conditions. For example, the topographic relief of karst mountain areas can be generalized into a relatively simple topographic model, and the heterogeneity of aquifers can be averaged to a certain extent to facilitate subsequent calculations.

[0021] Identification of key influencing factors: Determine the main factors affecting groundwater level changes in karst mountain areas, such as meteorological factors (including precipitation and evaporation), geological factors (including stratigraphic lithology and geological structure), and human activities (including groundwater extraction and engineering construction). Clarify how these factors interact with the groundwater system to influence the dynamic changes in groundwater levels.

[0022] The specific process of the conceptual model has been publicly disclosed and can be found on the Internet, so it will not be elaborated here.

[0023] S12. Based on the three-dimensional groundwater flow numerical model, long-term observation data of a specified karst mountain area are obtained from the data center. Correlation analysis is performed on the long-term observation data of the specified karst mountain area. Statistical methods are used to establish a statistical model for groundwater level prediction. The parameters of the model are determined through regression analysis. The statistical methods mentioned above include multiple linear regression and time series analysis. Statistical models for groundwater level prediction, for example, use precipitation and river flow as independent variables and groundwater level as the dependent variable. The specific process is publicly available online and will not be elaborated upon here.

[0024] S13. Optimize and correct the model using unused data from long-term observation data to obtain a water level prediction model.

[0025] In the above, the long-term observation data is divided into the model building dataset and the model optimization dataset. The model is built using the data in the model building dataset, and the unused data in the long-term observation data is the model optimization dataset. The specific model optimization process has been disclosed in existing technology and will not be elaborated here.

[0026] The long-term observation data consists of data collected over a long period of time from multiple groundwater level monitoring points set up in designated karst mountain areas using various monitoring equipment. This data includes underground geological structure, groundwater level, groundwater pH, groundwater conductivity, groundwater dissolved oxygen, surface water level, surface water flow rate, and surface water temperature. The monitoring equipment includes water level gauges, multi-parameter water quality analyzers, and ground-penetrating radar.

[0027] In another specific embodiment, the process of simulating underground excavation in a designated karst mountain area is as follows: obtain each excavation location, the excavation parameters of each excavation location, and the excavation depths for which groundwater prediction is required at each excavation location from the underground excavation plan, and use each excavation depth for which groundwater prediction is required as the prediction depth. The excavation parameters include the size and depth of the excavation.

[0028] In the three-dimensional model of the specified karst mountain area, excavation simulation is carried out sequentially according to each excavation location and its excavation parameters. When the underground excavation simulation is carried out, when the excavation depth of each excavation location reaches the predicted depth, the groundwater level at each excavation location is predicted by the water level prediction model. At the same time, the start time and end time of the prediction of the groundwater level at each excavation location and the predicted depth are obtained from the computer background records and used as the initial time and end time of the groundwater level prediction at each predicted depth in each excavation location. The actual predicted depth of the predicted groundwater level at each excavation location is obtained from the three-dimensional model of the specified karst mountain area. The actual predicted depth of the predicted groundwater level at each excavation location is subtracted from the actual predicted depth of the predicted groundwater level at each excavation location to obtain the excavation depth of each excavation location at the predicted groundwater level. The groundwater level at each excavation location when the excavation depth reaches each predicted depth, the excavation depth at each excavation location when the groundwater level at each predicted depth is predicted, and the initial and end times of the predicted groundwater level at each predicted depth in each excavation location are used as the water level prediction simulation data.

[0029] S2. Based on underground excavation simulation, predict the water level and analyze the water level prediction efficiency of the water level prediction model. If the water level prediction efficiency is poor, execute S4; if the water level prediction efficiency is good, execute S3. In a specific embodiment, the water level prediction efficiency of the water level prediction model is specifically analyzed as follows: using the groundwater level at each excavation location when the excavation depth reaches each prediction depth in the water level prediction simulation data, the prediction efficiency influence coefficient of each excavation location at each prediction depth is analyzed, and the prediction duration of the groundwater level at each prediction depth in each excavation location is obtained based on the initial time and end time of the groundwater level prediction at each prediction depth in each excavation location. The prediction efficiency influence coefficient of each excavation location at each prediction depth, the prediction time of groundwater level at each prediction depth in each excavation location, and the excavation depth of each excavation location at each prediction depth when predicting groundwater level are input into the water level prediction efficiency formula, and the water level prediction efficiency analysis results of the water level prediction model are output. The results of the water level prediction efficiency analysis include YES and NO. When the result is YES, it indicates that the water level prediction efficiency is good, and vice versa.

[0030] The analysis process for the prediction efficiency influence coefficient of each excavation location at each prediction depth, as described above, is as follows: The groundwater level at each excavation location when the excavation depth reaches the predicted depth is denoted as . , where r represents the number of each excavation location, y represents the number of each predicted depth, and both r and y are positive integers; Using analytical formulas: The prediction efficiency influence coefficient of the r-th excavation location at the y-th prediction depth is obtained. In the formula, and Let represent the groundwater level at the r-th excavation location when the excavation depth reaches the (y-1)-th predicted depth and the groundwater level at the r-th predicted depth, respectively. This represents the groundwater level at the r-th excavation location before excavation. This indicates the threshold value set for the difference in groundwater level changes. Let Y represent the prediction efficiency requirement factor for the r-th excavation location, and let R represent the number of prediction depths and the number of excavation locations, respectively. in, .

[0031] It should be noted that groundwater levels at each excavation location were collected using water level gauges before excavation. The threshold for the difference in groundwater level variation is a permissible value for groundwater level variation, jointly determined by multiple professionals based on the characteristics and capabilities of karst mountain areas. It can be 5 or 8, and no specific numerical limit is imposed here.

[0032] Preferably, the water level prediction efficiency formula is: In the formula, This indicates the results of the water level prediction efficiency analysis of the water level prediction model. This represents the excavation depth at the r-th excavation location when the groundwater level is predicted at the y-th predicted depth. This represents the depth difference between the predicted depths at the r-th excavation location, where T is the set prediction duration threshold. Indicates the prediction efficiency threshold. This represents the prediction duration of the groundwater level at the y-th predicted depth in the r-th excavation location.

[0033] It should be noted that the depth difference between the predicted depths at each excavation location is the average of the differences between adjacent predicted depths at each excavation location.

[0034] The prediction duration threshold is the maximum permitted prediction duration set by multiple professionals based on testing needs. It can be 1 minute or 3 minutes, and no specific numerical limit is set here.

[0035] The prediction efficiency threshold is a benchmark value for evaluating whether the prediction efficiency is good. When it is greater than or equal to the prediction efficiency threshold, it indicates that the prediction efficiency is good, and vice versa. It is set by multiple professionals according to the test requirements. It can be 2 or 6. There is no specific numerical limit here.

[0036] It should be noted that, Figure 1 and Figure 2 YES and NO represent water level prediction efficiency that is good and water level prediction efficiency that is poor, respectively.

[0037] S3. Use the water level prediction model to predict the water level in the underground excavation simulation, formulate preventive measures for underground excavation in designated karst mountainous areas, and then execute S4. In a specific embodiment, the process of formulating preventive measures for underground excavation in designated karst mountain areas is as follows: S31, using water level prediction simulation data, screening out each target excavation location, and obtaining geological data, groundwater level parameters, underground environmental data and excavation parameters of each target excavation location as construction correlation data for each excavation location, obtaining construction records corresponding to each historical excavation in each karst mountain area from the data center, and predicting the set of dangerous phenomena for each target excavation location; Preferably, the target excavation locations are selected by means of the following process: if the absolute value of the difference between the groundwater level at at least one predicted depth and the groundwater level at an adjacent predicted depth in a certain excavation location is greater than the threshold value of the difference in groundwater level changes, then the excavation location is taken as the target excavation location, thereby selecting the target excavation locations.

[0038] It should be noted that the geological data includes geological structures and topography. The groundwater level parameters are the groundwater levels at each predicted depth. The groundwater environmental data includes groundwater pH, groundwater conductivity, and dissolved oxygen.

[0039] Preferably, the specific process for predicting the set of hazardous phenomena at each target excavation location is as follows: geological data, groundwater level parameters, underground environmental data, and excavation parameters for each historical excavation location are obtained from the construction records corresponding to each historical excavation in each karst mountain area, which serve as construction correlation data for each historical excavation location. At the same time, the combination of hazardous phenomena, the hazard level of each hazardous phenomenon, prevention plan, and prevention level are obtained from the construction records corresponding to each historical excavation in each karst mountain area. Calculate the similarity between the construction association data of each historical excavation location in each karst mountain area and the construction association data of each target excavation location. Use the historical excavation locations of each historical excavation location in each karst mountain area with a similarity greater than a preset similarity threshold as reference excavation locations for each target excavation location. It should be noted that the similarity threshold and the threshold for groundwater level change difference are set similarly, and will not be elaborated here. The similarity between the construction association data of each historical excavation location corresponding to each historical excavation in each karst mountain area and the construction association data of each target excavation location is 1 minus the absolute value of the difference between the construction association data of each historical excavation location corresponding to each historical excavation in each karst mountain area and the construction association data of each target excavation location, and the ratio of the construction association data of each target excavation location.

[0040] Extract the hazardous phenomenon combinations corresponding to each reference excavation location for each target excavation location, integrate the hazardous phenomenon combinations corresponding to each reference excavation location for each target excavation location, and use the integrated result as the hazardous phenomenon set for each target excavation location.

[0041] S32. Based on the set of dangerous phenomena at each target excavation location, select a set of target prediction schemes for each target excavation location as preventive measures.

[0042] Preferably, the specific process of S32 is as follows: S321, based on the combination of hazardous phenomena corresponding to each reference excavation location, the hazard level of each hazardous phenomenon, the prevention plan and the prevention level, obtain the hazardous phenomena and the prevention level corresponding to each prevention plan in each target excavation location; at the same time, obtain the hazard level of each hazardous phenomenon corresponding to each prevention plan in each target excavation location, and cluster the hazard levels, and use the clustering result as the target hazard level of each hazardous phenomenon, thereby obtaining the target hazard level and the prevention level of each hazardous phenomenon corresponding to each prevention plan in each target excavation location. It should be noted that dangerous phenomena include sudden water and mud inrush, ground subsidence, and instability of surrounding rock. Prevention measures include strengthening the support structure, setting up a drainage system, and grouting to seal groundwater.

[0043] Hazard Level Assessment: An assessment index system is established, specifically by selecting key indicators related to hazardous phenomena, such as karst ratio, uniaxial compressive strength of rock, groundwater pressure, and overburden thickness. Professionals assign corresponding weights based on their impact on the degree of hazard. For example: karst ratio weight 0.3, uniaxial compressive strength of rock weight 0.25, groundwater pressure weight 0.3, and overburden thickness weight 0.15. Karst ratio, uniaxial compressive strength of rock, groundwater pressure, and overburden thickness can be obtained using ground-penetrating radar, pressure testing machines, hydraulic gauges, and drilling equipment, respectively.

[0044] Taking karst rate as an example, karst rate is divided into four levels: less than 5%, 5%-15%, 15%-30%, and greater than 30%, corresponding to danger scores of 100, 80, 60, and 30, respectively. The danger score corresponding to the karst rate is obtained based on the karst rate, and then multiplied by the karst rate weight of 0.3 to obtain the danger score of the karst rate. Following the process of obtaining the danger score of the karst rate, the danger scores of the uniaxial compressive strength of the rock, the groundwater pressure, and the overburden thickness are obtained separately, and then accumulated to obtain the total danger score.

[0045] Assume there are four hazard levels: Level 1 has a total hazard score range of 0-30, Level 2 has a total hazard score range of 31-60, Level 3 has a total hazard score range of 61-80, and Level 4 has a total hazard score range of 81-100. By comparing the total hazard score with the total hazard score range of each level, the hazard level is obtained.

[0046] The process of obtaining prevention levels is similar to that of obtaining hazard levels, and will not be elaborated upon here.

[0047] S322. Compare the hazardous phenomena prevented by each prevention scheme at each target excavation location with the hazardous phenomena in the hazardous phenomenon set at each target excavation location. Obtain the number of hazardous phenomena prevented by each prevention scheme at each target excavation location that are the same as the hazardous phenomena in the hazardous phenomenon set. Take this number as the number of effective hazardous phenomena prevented by each prevention scheme at each target excavation location. Take each hazardous phenomenon that is the same as the hazardous phenomena in the hazardous phenomenon set as each effective hazardous phenomenon. Extract the target hazard level and prevention level of each effective hazardous phenomenon prevented by each prevention scheme at each target excavation location. Calculate the priority value of each prevention scheme at each target excavation location. Select the prevention scheme with the highest priority value as the first target prevention scheme for each target excavation location. Preferably, the target hazard level and prevention level of each effective hazard phenomenon prevented by each prevention scheme at each target excavation location are calculated by weighted average. The calculation result is used as the target hazard level and prevention level of each effective hazard phenomenon prevented by each prevention scheme at each target excavation location. The number of effective hazard phenomena prevented by each prevention scheme at each target excavation location, the target hazard level of the effective hazard phenomenon, and the prevention level of the effective hazard phenomenon are respectively denoted as... , and Where f represents the number of each target excavation location, and p represents the number of each prevention plan, both f and p are positive integers. The priority value is calculated using the following formula: In the formula, This represents the priority value of the p-th prevention plan among the f-th target excavation locations. , These represent the weighting factors for the target hazard level and the prevention level, respectively.

[0048] in, , In the formula, This represents the maximum difference between the target hazard levels of each effective hazard phenomenon prevented by the p-th prevention scheme at the f-th target excavation location. This represents the minimum target hazard level among all the effective hazard phenomena prevented by the p-th prevention scheme at the f-th target excavation location. This represents the maximum difference between the prevention levels of each effective hazard phenomenon corresponding to the p-th prevention scheme at the f-th target excavation location. F represents the minimum prevention level among the prevention levels of each effective hazard phenomenon corresponding to the p-th prevention scheme in the f-th target excavation location, where F and P represent the number of target excavation locations and the number of prevention schemes, respectively.

[0049] S323. Remove the effective hazardous phenomena corresponding to the target prevention scheme from the hazardous phenomenon set of each target excavation location to obtain the second hazardous phenomenon set of each target excavation location. Then, compare the remaining hazardous phenomena corresponding to each prevention scheme in each target excavation location with the second hazardous phenomenon set of each target excavation location, and obtain the second target prevention scheme of each target excavation location according to the steps of S322. S324. Repeat S322-S323 until all hazardous phenomena in the hazardous phenomenon set of each target excavation location are removed. Stop repeating and obtain the target prevention plan for each target excavation location. The set of target prevention plans is used as the target prediction plan set.

[0050] This application's embodiments provide a reference for subsequent excavation construction protection by formulating preventive measures, which is conducive to preparing protective equipment in advance, reducing the negative impact of groundwater level changes, improving construction response speed, and ensuring construction safety.

[0051] S4. When the water level prediction efficiency is poor, feedback is sent to the visualization display terminal; the preventive measures for underground excavation in karst mountainous areas are sent to the visualization display terminal, and preventive prompts are given. Example 2:

[0052] See Figure 2 As shown, a groundwater level prediction system for karst mountain areas includes: The model building and prediction module is used to build a water level prediction model, obtain a 3D model of a specified karst mountain area and an underground excavation plan for the specified karst mountain area, simulate underground excavation in the specified karst mountain area, and use the water level prediction model to predict the water level of the underground excavation simulation. The model prediction and evaluation module is used for water level prediction based on underground excavation simulation. It analyzes the water level prediction efficiency of the water level prediction model. When the water level prediction efficiency is poor, the prediction execution module is executed. When the water level prediction efficiency is good, the prediction result processing module is executed. The prediction result processing module is used to predict the water level of the underground excavation simulation using the water level prediction model, formulate preventive measures for underground excavation in designated karst mountainous areas, and then execute the prediction execution module. The prediction execution module is used to provide feedback to the visualization display terminal when the water level prediction efficiency is poor; it also sends the specified preventive measures for underground excavation in karst mountainous areas to the visualization display terminal and provides preventive prompts.

[0053] The above description is merely an example and illustration of the concept of the present invention. Those skilled in the art can make various modifications or additions to the specific embodiments described or use similar methods to replace them, as long as they do not deviate from the concept of the invention or exceed the scope defined in this specification, they should all fall within the protection scope of the present invention.

Claims

1. A karst mountain groundwater level prediction method, characterized in that, The method comprises the following steps: S1, constructing a water level prediction model, obtaining a three-dimensional model of a specified karst mountain area and an underground excavation plan of the specified karst mountain area, simulating underground excavation of the specified karst mountain area, and predicting the water level of the underground excavation simulation by using the water level prediction model; S2, based on the water level prediction of the underground excavation simulation, analyzing the water level prediction efficiency of the water level prediction model, when the water level prediction efficiency is poor, executing S4, when the water level prediction efficiency is good, executing S3; S3, predicting the water level of the underground excavation simulation by using the water level prediction model, formulating preventive measures for the underground excavation of the specified karst mountain area, and then executing S4; S4, when the water level prediction efficiency is poor, feeding back to a visual display terminal; feeding the preventive measures for the underground excavation of the specified karst mountain area to the visual display terminal, and giving a prevention prompt. 2.The karst mountain groundwater level prediction method according to claim 1, characterized in that, The specific process of constructing the water level prediction model is: S11, obtaining geological data and hydrogeological data of the specified karst mountain area through geological exploration, constructing a conceptual model of the specified karst mountain area, dividing the specified karst mountain area into grids by using numerical simulation software, and establishing a three-dimensional numerical model of underground water flow; S12, based on the three-dimensional numerical model of underground water flow, obtaining long-term observation data of the specified karst mountain area from a data center, performing correlation analysis on the long-term observation data of the specified karst mountain area, establishing a statistical model for predicting the underground water level by using a statistical method, and determining the parameters of the model through regression analysis; S13, optimizing and correcting the model by using the data not used in the long-term observation data, and obtaining the water level prediction model. 3.The karst mountain groundwater level prediction method according to claim 1, characterized in that, The specific process of simulating the underground excavation of the specified karst mountain area is: obtaining each excavation position, excavation parameters of each excavation position, and each excavation depth needing to predict the underground water in each excavation position from the underground excavation plan, and taking each excavation depth needing to predict the underground water as each prediction depth; in the three-dimensional model of the specified karst mountain area, excavating in sequence according to each excavation position and its excavation parameters, when the excavation depth of each excavation position reaches each prediction depth during the underground excavation simulation, predicting the underground water level of each excavation position when the excavation depth of each excavation position reaches each prediction depth by using the water level prediction model, and obtaining the start prediction time and the end prediction time of the water level prediction model predicting the underground water level of each prediction depth in each excavation position from the computer background record as the initial time and the end time of the underground water level prediction of each prediction depth in each excavation position; and obtaining the actual prediction depth of the underground water level prediction end time of each prediction depth in each excavation position from the three-dimensional model of the specified karst mountain area, subtracting each prediction depth from the actual prediction depth of the underground water level prediction end time of each prediction depth in each excavation position, and obtaining the excavation depth of each prediction depth in each excavation position during the underground water level prediction; taking the underground water level of each excavation position when the excavation depth reaches each prediction depth, the excavation depth of each prediction depth in each excavation position during the underground water level prediction, and the initial time and the end time of the underground water level prediction of each prediction depth in each excavation position as the water level prediction simulation data.

4. The karst mountain groundwater level prediction method of claim 3, wherein, The specific analysis process of analyzing the water level prediction efficiency of the water level prediction model is: The water level prediction simulation data is used to simulate the underground water level of each excavation position when the excavation depth reaches each predicted depth, analyze the prediction efficiency influence coefficient of each excavation position at each predicted depth, and obtain the prediction time length of the underground water level of each excavation position at each predicted depth according to the initial time and the end time of the underground water level prediction of each excavation position at each predicted depth; The prediction efficiency influence coefficient of each excavation position at each predicted depth, the prediction time length of the underground water level of each excavation position at each predicted depth, and the excavation depth of each excavation position at each predicted depth are input into the water level prediction efficiency formula, and the water level prediction efficiency analysis result of the water level prediction model is output; The water level prediction efficiency analysis result includes YES and NO. When the water level prediction efficiency analysis result is YES, it indicates that the water level prediction efficiency is good, otherwise it indicates that the water level prediction efficiency is poor.

5. The karst mountain groundwater level prediction method of claim 4, wherein, The analysis process of the prediction efficiency influence coefficient of each excavation position at each predicted depth is as follows: The underground water level when the excavation depth of each excavation position reaches each predicted depth is denoted as wherein r represents the number of each excavation position, y represents the number of each predicted depth, and r and y are both positive integers. The predicted efficiency influence coefficient of the rth excavation position at the yth predicted depth is obtained by using an analysis formula: , wherein, and respectively represent the underground water level when the excavation depth of the rth excavation position reaches the (y-1)th predicted depth and the underground water level when the excavation depth reaches the 1st predicted depth, represents the underground water level when the rth excavation position is not excavated, represents a threshold value of the set underground water level change difference, represents a predicted efficiency demand factor of the rth excavation position, Y and R respectively represent the number of predicted depths and the number of excavation positions. wherein .

6. The karst mountain groundwater level prediction method of claim 5, wherein, The water level prediction efficiency formula is: , wherein, represents the water level prediction efficiency analysis result of the water level prediction model, represents the excavation depth of the rth excavation position when the water level is predicted at the yth prediction depth, represents the depth difference between the prediction depths in the rth excavation position, and T is a set prediction time threshold, represents a prediction efficiency threshold, represents the prediction time of the water level at the yth prediction depth in the rth excavation position.

7. The karst mountain groundwater level prediction method of claim 3, wherein, The specific process of formulating the preventive measures for underground excavation in the designated karst mountainous area is as follows: S31, using water level prediction simulation data, screening out each target excavation position, and obtaining geological data, underground water level parameters, underground environment data and excavation parameters of each target excavation position as construction correlation data of each excavation position, obtaining construction records corresponding to each historical excavation of each karst mountainous area from the data center, and predicting a dangerous phenomenon set of each target excavation position; S32, based on the dangerous phenomenon set of each target excavation position, screening out a target prediction scheme set of each target excavation position as a preventive measure.

8. The karst mountain groundwater level prediction method of claim 7, wherein, The specific process of predicting the dangerous phenomenon set of each target excavation position is as follows: Obtain geological data, underground water level parameters, underground environment data and excavation parameters of each historical excavation position from the construction records corresponding to each historical excavation of each karst mountainous area as the construction correlation data of each historical excavation position, and at the same time obtain the dangerous phenomenon combination of each historical excavation position, the dangerous level of each dangerous phenomenon, the preventive scheme and the preventive level from the construction records corresponding to each historical excavation of each karst mountainous area; Calculate the similarity of the construction correlation data of each historical excavation position corresponding to each historical excavation of each karst mountainous area and the construction correlation data of each target excavation position, and take each historical excavation position corresponding to each historical excavation of each karst mountainous area with a similarity greater than a preset similarity threshold as each reference excavation position of each target excavation position; Extract the dangerous phenomenon combination of each target excavation position corresponding to each reference excavation position, integrate the dangerous phenomenon combination of each target excavation position corresponding to each reference excavation position, and the integration result is taken as the dangerous phenomenon set of each target excavation position.

9. The karst mountain groundwater level prediction method of claim 8, wherein, The specific process of S32 is as follows: S321, based on each target excavation position corresponding to each reference excavation position, the combination of dangerous phenomena, the dangerous level of each dangerous phenomenon, the prevention scheme and the prevention level, obtain each prevention scheme in each target excavation position corresponding to each dangerous phenomenon and the prevention level; At the same time, obtain each dangerous phenomenon corresponding to each prevention scheme in each target excavation position, and cluster each dangerous level, and the clustering result is the target dangerous level of each dangerous phenomenon, so as to obtain the target dangerous level of each dangerous phenomenon corresponding to each prevention scheme in each target excavation position and the prevention level of each dangerous phenomenon; S322, compare each dangerous phenomenon corresponding to each prevention scheme in each target excavation position with each dangerous phenomenon in the dangerous phenomenon set of each target excavation position, obtain the number of dangerous phenomena corresponding to each prevention scheme in each target excavation position which is the same as the dangerous phenomenon set, as the number of effective dangerous phenomena corresponding to each prevention scheme in each target excavation position, and the dangerous phenomenon which is the same as the dangerous phenomenon set is regarded as each effective dangerous phenomenon, and the target dangerous level of each effective dangerous phenomenon corresponding to each prevention scheme in each target excavation position and the prevention level of each dangerous phenomenon are extracted, the priority value of each prevention scheme in each target excavation position is calculated, and the prevention scheme with the maximum priority value is selected as the first target prevention scheme of each target excavation position; S323, remove the effective dangerous phenomenon corresponding to the target prevention scheme in the dangerous phenomenon set of each target excavation position, obtain the second dangerous phenomenon set of each target excavation position, then compare each dangerous phenomenon corresponding to each remaining prevention scheme in each target excavation position with the second dangerous phenomenon set of each target excavation position, and obtain the second target prevention scheme of each target excavation position according to the steps of S322; S324, repeat S322-S323 until the dangerous phenomena in the dangerous phenomenon set of each target excavation position are removed, stop repeating, and obtain each target prevention scheme of each target excavation position, and the set of each target prevention scheme is regarded as the target prediction scheme set.

10. A karst mountain groundwater level prediction system for performing the karst mountain groundwater level prediction method according to any one of claims 1 to 9, characterized by, It comprises: The model construction and prediction module is used for constructing the water level prediction model, obtaining the three-dimensional model of the specified karst mountainous area and the underground excavation plan of the specified karst mountainous area, simulating the underground excavation of the specified karst mountainous area, and predicting the water level of the underground excavation simulation by using the water level prediction model; The model prediction evaluation module is used for analyzing the water level prediction efficiency of the water level prediction model based on the water level prediction of the underground excavation simulation, executing the prediction execution module when the water level prediction efficiency is poor, and executing the prediction result processing module when the water level prediction efficiency is good; The prediction result processing module is used for predicting the water level of the underground excavation simulation by using the water level prediction model, formulating the prevention measures for the underground excavation of the specified karst mountainous area, and then executing the prediction execution module; The prediction execution module is used for feeding back to the visual display terminal when the water level prediction efficiency is poor; the prevention measures for the underground excavation of the specified karst mountainous area are fed into the visual display terminal, and the prevention prompt is performed.

Images