Multi-source data fused refined treatment decision-making method for complex stratum disaster source

Abstract

The invention belongs to the technical field of tunnel construction geological disaster prevention and control, and discloses a multi-source data fused refined treatment decision-making method for a complex stratum disaster source, which comprises the following steps: collecting and fusing multi-source geological data, and constructing a three-dimensional geological model; generating a disaster source risk dynamic assessment and treatment scheme; based on a fluid-solid coupling similarity theory, verifying the preliminary treatment scheme by adopting a physical model test, and determining an optimal treatment scheme; the optimal treatment scheme is executed, and the treatment process is dynamically regulated and controlled; after treatment, the treatment effect is evaluated through posterior data, and the effect data is fed back to the three-dimensional geologic model and the knowledge base, so that the dynamic updating of the model and the self-learning of the decision-making system are realized. By the adoption of the treatment decision method, the problems that a traditional method depends on experience, information is one-sided, and treatment is extensive are solved, advanced accurate forecasting and refined and personalized treatment of complex stratum disaster sources are achieved, and the safety and efficiency of tunnel construction are remarkably improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of tunnel construction geological disaster prevention, and particularly relates to a fine treatment decision method for complex stratum disaster sources based on multi-source data fusion. BACKGROUND

[0002] In the process of tunnel construction, disaster sources (such as karst caves, fractures, water-rich zones, etc.) in complex strata (such as karst strata, water-rich broken strata, etc.) are prone to cause water inrush, mud inrush, surrounding rock collapse and other accidents, which seriously threaten construction safety and engineering progress. Traditional disaster source detection methods mostly rely on single geophysical prospecting technology or geological drilling, and have problems of low detection accuracy and limited coverage. At the same time, the treatment scheme is often based on experience, lacks comprehensive support of multi-source data, and is difficult to achieve fine treatment, resulting in poor treatment effect or resource waste.

[0003] In the prior art, although some methods try to combine multiple geophysical prospecting technologies for disaster source detection, they do not form a systematic multi-source data fusion mechanism, and cannot effectively integrate geophysical prospecting data, geological data and construction monitoring data. In addition, the combination degree of model test and numerical simulation is insufficient, and it is difficult to accurately reflect the dynamic response law of disaster sources in the construction process, resulting in a lack of scientific basis for treatment decision.

[0004] Therefore, an integrated method for fine detection and treatment decision of disaster sources based on multi-source data fusion is urgently needed. SUMMARY

[0005] The purpose of the present application is to provide a fine treatment decision method for complex stratum disaster sources based on multi-source data fusion, which realizes fine management and control of complex stratum disaster sources through multi-source data acquisition, fusion modeling, model test verification, risk assessment and treatment decision.

[0006] To achieve the above purpose, the present application provides a fine treatment decision method for complex stratum disaster sources based on multi-source data fusion, comprising the following steps: Step S1, a plurality of geophysical prospecting methods and geological drilling are combined to collect multi-source geological data of the construction area; Step S2, the obtained multi-source geological data is fused at data level, feature level and decision level to construct a three-dimensional geological model reflecting the stratum structure and spatial distribution of disaster sources; Step S3, based on the three-dimensional geological model, the response of surrounding rock and disaster sources under construction disturbance is analyzed by numerical simulation to dynamically evaluate the disaster risk level; combined with the knowledge base and case base, a preliminary treatment scheme is generated; Step S4, based on the fluid-solid coupling similarity theory, a physical model test is used to verify the preliminary treatment scheme, and the treatment scheme is optimized according to the test results to determine the optimal treatment scheme; Step S5, during the construction process, the optimal treatment scheme is executed, and the real-time parameters of the construction machinery and the field monitoring data are used to dynamically control the treatment process; Step S6, after treatment, the treatment effect is evaluated through post-data, and the effect data is fed back to the three-dimensional geological model and knowledge base, realizing dynamic updating of the model and self-learning of the decision system.

[0007] Preferably, in step S1, multi-scale and multi-dimensional data acquisition of complex stratum disaster sources is realized through space-air-ground-hole integrated survey technology, and the specific process is as follows: Step S11, in the early stage of the project, water area transient electromagnetic method or sea area multi-channel shallow seismic method is used to scan and detect the project area, and the distribution of karst development section is macroscopically mastered, and the project area is preliminarily divided into different sections of weak, medium and strong karst development; Step S12, according to the division result, drill holes are arranged in different sections and corresponding geophysical method combinations are selected, specifically including: (1) Karst strong / medium development area: elastic wave cross-hole CT as the core method; (2) Weak karst development area: supplemented by tube wave detection method and multi-frequency borehole sonar; Step S13, during the multi-source data acquisition process, control geological drill holes are arranged at all geophysical points, core sampling, logging are carried out, and rock mineral composition analysis and physical and mechanical tests are carried out to obtain the real and in-situ parameters of rock-soil mass, providing calibration points for geophysical interpretation.

[0008] Preferably, in step S2, the obtained multi-source heterogeneous data is fused and processed to build a three-dimensional geological model, and the specific process is as follows: Step S21, data level fusion; Align and calibrate the original data or preliminary processing results of different geophysical methods in a unified three-dimensional spatial coordinate system; Step S22, feature level fusion; Firstly, the feature information of the disaster source is extracted from various data, including: extracting low-speed abnormal area and wave velocity value from elastic wave CT data; extracting lithology change surface, RQD index and actual revealed karst cave position from drill hole data; extracting three-dimensional point cloud of karst cave wall from borehole sonar data; Then, the obtained feature information is combined into a comprehensive feature vector for jointly describing the same geological unit or disaster source; Step S23, decision level fusion and three-dimensional geological model construction; Step S231, basic geological model construction: based on drill hole data and geological profile, through rock-soil mass automatic / semi-automatic modeling tool, create terrain surface, stratum unit, weathering surface basic geological model; Step S232, disaster source model construction: the wave velocity three-dimensional numerical simulation analysis model obtained by elastic wave CT inversion is imported into the geological three-dimensional survey design system; the system automatically circumscribes the development range of the karst cave in the three-dimensional space by identifying the wave velocity abnormal threshold; the range is used as a constraint to construct a three-dimensional model of the disaster source by using discrete point modeling or a karst cave professional modeling tool; Step S233, model integration: the tunnel design model is integrated with the basic geological model and the disaster source model to construct a three-dimensional geological model, and the spatial relative position relationship between the tunnel and the disaster source is intuitively displayed.

[0009] Preferably, in step S3, the dynamic influence of the tunnel construction on the surrounding rock and the disaster source is simulated and predicted in advance by using the three-dimensional geological model constructed by numerical simulation technology, and the specific process is as follows: Step S311, simulation model establishment; Firstly, the three-dimensional geological model constructed is directly imported into the numerical simulation software; different geological bodies in the model are automatically assigned with corresponding rock and soil mechanical parameters obtained through experiments; Then, the whole process of shield tunneling is simulated, including: the disturbance of cutter head rotation cutting to the soil in front, the friction between the shield shell and the surrounding rock, the formation of the shield tail gap, and the pressure and filling effect of synchronous grouting; Finally, the fluid-structure coupling analysis process is defined in the model to simulate the interaction between groundwater and rock and soil in the water-rich karst stratum; Step S312, after calculation, the response information of the surrounding rock and the disaster source in the whole construction process is obtained, and multi-physical field response analysis is carried out, specifically including: (1) stress field evolution: analysis of the plastic zone range of the surrounding rock, development law, and stress concentration and release at key positions; (2) displacement field evolution: prediction of ground subsidence tank, tunnel convergence deformation, and advanced displacement of soil in front of the working face; (3) seepage field evolution: simulation of the seepage path and pressure distribution change of confined water in the karst cave-surrounding rock-tunnel system, and prediction of potential water inrush channels; (4) comprehensive identification: through the analysis of the coupling effect of the physical field, it is identified whether the water-resisting rock mass will be unstable, seepage unstable or overall damaged under the construction disturbance.

[0010] Preferably, in step S3, based on the output results of numerical simulation, the disaster risk is quantitatively and hierarchically dynamically evaluated, and the specific process is as follows: The evaluation indexes of safety factor, plastic zone penetration, displacement threshold, and seepage velocity / pressure gradient are established; the improved risk matrix method is used to combine the probability and consequences of disaster occurrence; wherein, the calculation method of risk value is: ; in, It represents the probability of disaster occurrence, based on numerical simulation results and determined through a comprehensive assessment of multiple indicators; This represents the plastic zone penetration rate, i.e., the number of plastic zone elements / the total number of elements; when hour, Otherwise, the value is taken proportionally. This indicates the proportion of monitoring points whose displacement exceeds the critical displacement threshold; This indicates the proportion of seepage velocity or pressure gradient exceeding a critical value; , and Let be the weight of each indicator, and ; ; in, Indicates the severity of the disaster consequences; Indicates the severity of the disaster's consequences for personnel; Indicates the severity of the consequences of equipment damage; Indicates the severity of the consequences of disasters during the construction period; Indicates the severity of the consequences of an environmental disaster; ; in, Indicates the risk value; Set decision threshold The risk levels are classified as follows: when At that time, the risk level was low, and normal tunneling was underway; when At that time, the risk level was medium, and monitoring was strengthened and preparations were made for treatment. when If the risk level is high, take immediate action and implement the treatment plan. when At that time, the risk level was extremely high, and the design scheme was optimized. Based on the characteristics and risk assessment level of the disaster source, one or more preliminary treatment plans are matched and recommended from the knowledge base; after generating multiple preliminary treatment plans, a multi-objective decision-making method is used for comparison and selection; the comprehensive scoring function is as follows: ; in, Indicates the first The overall score of each treatment plan; Indicates the first The weight of each evaluation indicator; Indicates the first The treatment plan was in the first a normalized score on one evaluation index; denotes the total number of evaluation indexes; Finally, the comprehensive score is selected The highest scheme is selected as the optimal preliminary treatment scheme. Then, the recommended treatment scheme and its parameters are re-substituted into the numerical model for digital twin simulation to predict the effect after treatment; through the cycle of calculation-comparison-optimization, the optimal preliminary treatment scheme is selected from multiple alternative schemes.

[0011] Preferably, in step S4, the preliminary treatment scheme is first verified based on the fluid-structure coupling similarity theory using physical model tests, and the specific process is as follows: Step S411, based on the fluid-structure coupling similarity theory, the target mechanical parameters of the ideal similar material of the model test are calculated; Step S412, a CBCS type fluid-structure coupling similar material is used, which is composed of aggregate, binder, adjusting agent and mixing agent; Through proportioning test, the final proportioning is determined, so that the physical and mechanical parameters of the prepared similar material strictly meet the target mechanical parameter values of the ideal similar material, ensuring the similarity of the model and the prototype in mechanics and seepage behavior; Then, a physical model test system is constructed to simulate the process of shield tunneling and treatment in water-rich karst strata, which specifically includes: Step S421, a three-dimensional visual model box is assembled using modular steel structure, and multiple layers of visual observation windows are arranged on the side; a shield machine excavation entrance is reserved, and the boundary distance from the tunnel diameter is greater than 3 times to eliminate the boundary effect; Step S422, a composite EPB shield tunneling simulation system is established, which has the functions of cutter head rotation tunneling, hydraulic jacking, soil pressure warehouse pressure monitoring and double auger unearthing; the tunneling parameters are set and automatically controlled according to the similarity ratio, and are recorded in real time; Step S423, based on the real cave shape determined by the three-dimensional geological model, a disaster source model similar in geometry to the prototype is printed using 3D printing technology; Step S424, a water pressure intelligent regulation and control and multi-element information monitoring system is constructed, a water pressure loading device is connected with the disaster source model to apply and maintain a constant water pressure similar to the prototype water pressure; a micro soil pressure cell, a seepage pressure gauge, a waterproof strain brick and a grating multi-point displacement meter are embedded inside the model body to monitor the evolution of stress, seepage pressure, strain and displacement in real time; a three-dimensional laser scanner is used to scan the model surface at high frequency and non-contact to construct the settlement cloud picture of the entire surface during the construction process; all tunneling parameters of the shield machine are recorded, and the unearthing state and the process of sudden gushing water are recorded automatically; Finally, the implementation, verification and optimization of the treatment scheme, specifically including: Step S431, test implementation: First, the generated preliminary treatment scheme is implemented in a model test; then, the shield machine is started and tunneling is performed according to the set parameters until the disaster source influence area is passed through; Step S432, scheme verification: Tunneling is performed under the premise of applying the treatment scheme; by comparing the monitoring data before and after treatment, the reliability of the scheme is verified; Step S433, scheme optimization and optimal scheme determination: Parameter inversion: if the treatment effect does not fully meet the expectations, the deficiencies of the treatment scheme are analyzed based on the monitored failure data in the test; Iterative optimization: based on the analysis results, the parameters of the treatment scheme are adjusted, and the test verification is performed again.

[0012] Preferably, in step S5, an intelligent construction execution system is built, and the specific process is as follows: Step S511, scheme parameter digitalization issuance: The key process parameters of the optimal treatment scheme are digitized and directly issued to the corresponding intelligent construction equipment control system; Step S512, intelligent equipment precise execution: Intelligent grouting: after receiving the instruction, the grouting equipment automatically grouts according to the preset pressure-flow curve, and records the cumulative grouting amount in real time, ensuring that the slurry is accurately filled to the designed range; Shield intelligent tunneling: when the shield machine passes through the disaster source influence area, the control system preferentially adopts the recommended parameters for automatic tunneling mode, reducing the uncertainty of human operation; Mechanized support installation: using automatic equipment such as a robot, the support components are installed according to the designed position in the three-dimensional model.

[0013] Preferably, a digital sensory system is built to comprehensively and real-time perceive the construction environment and structure state, real-time collect shield machine parameters and grouting system data, and build an excavation face stability index and a grouting fullness index; Based on the real-time parameters of the shield machine, an excavation face stability index is built As shown below: ; Wherein, represents the measured soil bin pressure; represents the theoretically calculated static soil pressure; represents the measured cutterhead torque; represents the set torque; represents the actual soil volume; represents the theoretical excavation soil volume; , and These are weighting coefficients, reflecting the degree of influence of each parameter on stability; The regulatory decision is as follows: when At that time, the condition was stable, and tunneling continued according to the original parameters; when If the excavation face is at risk of instability, the system will automatically prompt you to increase the pressure in the soil chamber. Or reduce the tunneling speed; when At this time, mud cake formation or over-excavation may occur, and the system will automatically prompt you to reduce torque. Or adjust the dosage of soil amendment; Based on real-time grouting data, a grouting fullness index is constructed. As shown below: ; in, This indicates the cumulative actual grouting volume; This represents the theoretical grouting volume calculated based on the shield tail gap; This represents the integral of the actual grouting pressure over time. This represents the integral of the set grouting pressure over time; Decision-making rules: like If the grouting is complete, it can be continued normally. like Then dynamic control will be activated: the system will adjust according to... Automatically increases grouting pressure ;in, The gain coefficient is until... Meets the standards; Step S522, Geological and structural state perception; Internal sensing: In the rock mass surrounding and ahead of the tunnel, deploy or utilize advanced geological prediction systems to acquire in real time: microseismic / acoustic emission signals to monitor the development of rock mass fractures; deep displacement to monitor stratum slippage; and pore water pressure changes to monitor seepage field disturbances. Surface sensing: Inside the tunnel, convergence meters and total stations are used to automatically monitor tunnel convergence deformation; on the surface, surveying robots, GB-InSAR, or 3D laser scanners are used to continuously monitor changes in the surface settlement field.

[0014] Preferably, the real-time sensed data is compared with the prediction model and control thresholds to achieve dynamic optimization and adjustment of construction parameters, specifically including: Step S531, Data Fusion and Status Diagnosis: The acquired multi-element real-time data is fused and analyzed in a unified data platform, and is compared with the safety threshold predicted by numerical simulation in real time to diagnose the current construction state; In step S532, the system generates a control instruction or issues a warning to an operator based on the diagnosis result, and gives a control suggestion: Shield parameter control: If the diagnosis finds that the front rock mass is weak and deforms greatly, the system automatically or prompts the operator to appropriately reduce the tunneling speed and increase the soil chamber pressure to stabilize the excavation face; Treatment parameter strengthening: If the diagnosis finds that the deformation of the grouting reinforcement area is still beyond expectation, the system instructs the intelligent grouting system to dynamically increase the grouting pressure or perform compensatory grouting in a specific area; Support dynamic strengthening: If the initial support stress or deformation is close to the warning value, the system prompts to strengthen the support parameters in the next cycle; Step S533, risk warning and emergency intervention: A multi-level warning mechanism is established, and when the risk level reaches orange or red after multi-element data fusion, the system not only adjusts the parameters, but also automatically triggers an audible and light alarm, and sets an emergency shutdown interlock to force manual inspection and decision-making.

[0015] Preferably, in step S6, the effect evaluation and model updating, the specific process is as follows: Step S61, objectively and quantitatively evaluate the effect of the executed treatment scheme through the post-verification data collected after construction; Step S611, post-verification data collection; (1) Direct exposure verification includes: drilling core verification and borehole television / optical imaging; (2) Non-destructive testing verification includes: cross-hole CT re-measurement and seismic wave method / surface wave method; (3) Long-term stability monitoring: continue to collect tunnel convergence, ground settlement, and support structure stress monitoring data for a period of time to confirm whether the deformation and stress have stabilized and are far below the control value; Step S612, establish an effect evaluation index system based on the post-verification data, specifically including: (1) Filling degree index: based on coring and borehole television, evaluate the filling percentage of slurry to cavities / cracks; (2) Strength improvement rate: the ratio of the rock mass strength / wave velocity value after treatment to the value before treatment; (3) Integrity index: area reduction rate or wave velocity improvement rate of abnormal areas in the CT image after treatment; (4) Stability index: deformation convergence rate and final stable value after treatment; According to the above indexes, the comprehensive effect score of treatment is given to give quantitative conclusions of excellent, qualified, and unqualified; wherein, the treatment effect quantitative scoring formula is: ; wherein, represents the comprehensive effect score; and respectively represent the elastic wave velocities before and after treatment; and respectively represent the areas of low-velocity abnormal zones in the CT images before and after treatment; represents the final displacement value after treatment; represents the allowable displacement value; , and respectively represent the weights of intensity, integrity, and stability; Step S62, dynamic updating of the three-dimensional geological model and the numerical analysis model; Step S621, updating of the geological attribute model; Updating of the disaster source model: updating the attributes of the disaster source model after treatment from untreated to grout-filled or reinforced, and correcting the geometric boundaries and physical and mechanical parameters according to the verification data; Updating of the rock mass parameter partition: according to the strength improvement rate and the wave velocity improvement rate, the rock mass quality grade and the mechanical parameters of the reinforced area are partitioned and corrected in the three-dimensional geological model; Step S622, correction of the numerical analysis model; The updated three-dimensional geological model is reimported into the numerical simulation software, and the actual deformation and stress data monitored after treatment are taken as the target to perform parameter inversion on the constitutive model parameters in the numerical model, so that the results of numerical simulation are consistent with the actual situation; Step S63, self-learning and evolution of the knowledge base and the case base; Successful treatment experience is converted into system storage data as historical experience data; Step S631, creating a new successful treatment case, the data structure of which includes: (1) problem description: geological conditions before treatment, disaster source characteristics, risk assessment level, etc.; (2) treatment scheme: all details of the final adopted optimal treatment scheme; (3) effect data: all collected post-treatment data and effect evaluation conclusions; (4) context information: project name, stratum type, hydrological condition; Step S632, optimization and generation of knowledge rules; Verification and strengthening of existing rules: if the current successful case highly matches a recommended rule in the knowledge base, the confidence of the rule will be improved; New rule generation: if a new type of optimization treatment scheme is adopted this time and excellent results are obtained, a new treatment rule is automatically or semi-automatically generated; Failure lesson record: if the treatment effect is unqualified, the case is also recorded in the warehouse, and the reason is analyzed, which is used to issue a warning when similar situations are encountered in the future to avoid repeating the same mistakes.

[0016] Therefore, the fine treatment decision method of the fusion multi-source data of the complex stratum disaster source has the beneficial effects as follows: (1) Through multi-source geophysical data fusion and three-dimensional geological modeling, the precise positioning of the disaster source is realized, which can not only determine the position, but also accurately depict the three-dimensional spatial form, scale, filling state and spatial relationship with the tunnel, overcoming the one-sidedness of the information of a single geophysical method.

[0017] (2) Using numerical simulation based on the real geological model, the response of the surrounding rock and the disaster source under different construction schemes can be dynamically and quantitatively predicted before construction. This makes the risk assessment of water inrush, collapse and other risks no longer stay in qualitative experience judgment, but based on quantitative indexes such as safety factor, displacement field and plastic zone range, realizing the calculable and predictable risk.

[0018] (3) The scheme generation combines the prediction results of numerical simulation and expert knowledge base, and is verified and optimized through physical model test. This ensures that the treatment scheme is a repeatedly verified, personalized optimal solution, which ensures safety while improving economy.

[0019] The technical solutions of the present application will be further described in detail below with the help of the accompanying drawings and examples. BRIEF DESCRIPTION OF DRAWINGS

[0020] Figure 1 is the flow chart of the fine treatment decision method of the fusion multi-source data of the complex stratum disaster source of the present application; Figure 2 is the flow chart of the multi-source data acquisition of the present application; Figure 3 is the flow chart of the fusion processing of the multi-source data of the present application; Figure 4 is the three-dimensional visualization model box schematic diagram of the present application; Figure 5 is the cloud image of the ground settlement caused by shield construction. DETAILED DESCRIPTION

[0021] The technical solutions of the present application will be further described in detail below with the help of the accompanying drawings and examples.

[0022] As Figure 1As shown, the application is a fine treatment decision method of complex stratum disaster source fusion multi-source data, which comprises the following steps: Step S1, multi-source data acquisition: comprehensive multi-geophysical methods and geological drilling are used to collect multi-source geological data of the construction area.

[0023] Through the "space-air-ground-hole" integrated survey technology, multi-scale and multi-dimensional data acquisition of complex stratum, especially marine karst disaster source, is realized, such as Figure 2 As shown.

[0024] Step S11, in the early stage of the project, the water transient electromagnetic method or the marine multi-channel shallow seismic method is used to scan and detect the project area on a large scale, and the distribution of the karst development section is macroscopically mastered. The project area is preliminarily divided into different sections such as "karst micro-development, medium development, and strong development".

[0025] Step S12, according to the division result, drill holes are arranged in different sections and corresponding geophysical method combinations are selected, which specifically include: (1) Karst strong / medium development area: elastic wave cross-hole CT as the core method.

[0026] A PVC pipe with an outer diameter of ≥91mm is lowered into the drill hole. After pulling out the casing, one drill hole is used as a transmitting hole and the other as a receiving hole. The elastic wave is excited by using an electric spark source at a preset point distance (1.0m), and the signal is received by using a string detector in the receiving hole. The drill hole spacing is controlled within 25 meters to ensure signal penetration and data quality.

[0027] The elastic wave cross-hole CT method uses the good conductivity of seawater to promote the energy release of the electric spark source, overcomes the shielding effect of electromagnetic method in the sea area, and can accurately find out the planar distribution of inter-hole karst caves and karst fissures.

[0028] (2) Karst micro-development area: supplemented by tube wave detection method and multi-frequency borehole sonar.

[0029] The tube wave detection method uses a super-magnetic vibration source to generate Stoneley waves in the drill hole. By analyzing the reflection waves of the wave impedance difference interface around the hole, the un-revealed karst caves within a radius of 0.6-1.0 meters around the drill hole are detected. The multi-frequency borehole sonar lowers the sonar transducer to the karst cave section and performs 360° horizontal rotation scanning. By emitting sound wave pulses and receiving echo signals, the three-dimensional spatial form of the single drill hole revealed karst cave is directly measured and fitted.

[0030] Step S13, during the multi-source data acquisition process, control geological drill holes are arranged at all geophysical points, rock cores are taken, recorded, and rock mineral composition analysis and physical and mechanical tests are carried out to obtain the real and in-situ parameters of the rock-soil body, and provide "calibration points" for geophysical interpretation.

[0031] Step S2, multi-source data fusion: data-level, feature-level or decision-level fusion is performed on the acquired multi-source geological data to construct a three-dimensional geological model reflecting the stratum structure and spatial distribution of disaster sources, as shown in FIG. 2. Figure 3

[0032] Step S21, data-level fusion.

[0033] Raw data or preliminary processing results of different geophysical methods (such as elastic wave CT and tube wave) are aligned and calibrated in a unified three-dimensional spatial coordinate system.

[0034] Step S22, feature-level fusion.

[0035] Firstly, feature information of disaster sources is extracted from various data, including: low-velocity anomaly zone and wave velocity value are extracted from elastic wave CT data; lithological change surface, RQD index and actually exposed karst cave position are extracted from drilling data; and karst cave cavity wall three-dimensional point cloud is extracted from drilling sonar data.

[0036] Then, the obtained feature information is combined into a comprehensive feature vector for jointly describing a same geological unit or disaster source.

[0037] Step S23, decision-level fusion and three-dimensional geological model construction.

[0038] Step S231, basic geological model construction: based on drilling data and geological profile, a terrain surface, a stratum unit, a weathering surface and other basic geological models are created through an automatic / semi-automatic modeling tool.

[0039] Step S232, three-dimensional model construction of disaster sources (karst caves): the wave velocity three-dimensional numerical simulation model obtained by elastic wave CT inversion is imported into a geological (rock-soil) three-dimensional survey design system (such as GeoStation). The system automatically delineates the development range of karst caves in three-dimensional space by identifying the wave velocity anomaly threshold. With this range as a constraint, a discrete point building or a karst cave professional modeling tool is used to construct a three-dimensional model of karst caves with irregular shape and accurate spatial position, and the model is divided into full filling, half filling and no filling types.

[0040] Step S233, model integration: the tunnel design model (such as BIM model) is integrated with the basic geological model and the disaster source model to construct a three-dimensional geological model, which intuitively shows the spatial relative position relationship between the tunnel and the disaster sources, realizes comprehensive analysis from macro to micro, from two-dimensional profile to three-dimensional space, from single data interpretation to mutual verification of multi-source information, and provides a reliable and detailed “transparent geological” base map for subsequent risk assessment and treatment decision.

[0041] ​Step S3, disaster source risk dynamic assessment and treatment scheme generation: based on the three-dimensional geological model, the response of the surrounding rock and the disaster source under construction disturbance is analyzed through numerical simulation, and the disaster risk level is dynamically assessed; combined with the knowledge base and the case base, the preliminary treatment scheme is generated.

[0042] Step S31, construction disturbance response analysis based on numerical simulation.

[0043] Using the high-precision three-dimensional geological model constructed, the dynamic influence of tunnel construction (such as shield tunneling) on surrounding rock and disaster source is simulated and predicted in advance through numerical simulation technology.

[0044] Step S311, simulation model establishment.

[0045] Firstly, the three-dimensional geological model (including accurate karst cave morphology, stratum interface, fault, etc.) constructed in step S2 is directly imported into numerical simulation software (such as ANSYS, FLAC3D, ABAQUS, etc.). Different geological bodies in the model are automatically assigned with corresponding rock and soil mechanical parameters (such as elastic modulus, Poisson's ratio, cohesion, internal friction angle, permeability coefficient, etc.) obtained through experiments in step S2.

[0046] Then, the whole process of shield tunneling is accurately simulated, including but not limited to: disturbance of cutter head rotation cutting to the front soil, friction between shield shell and surrounding rock, formation of shield tail gap, pressure and filling effect of synchronous grouting.

[0047] Finally, the fluid-structure coupling analysis step is defined in the model to simulate the interaction between groundwater and rock-soil in water-rich karst stratum, especially the seepage path and pressure change of confined cave water under construction disturbance.

[0048] Step S312, multi-physical field response analysis.

[0049] Through calculation, the key response information of surrounding rock and disaster source during the whole construction process is obtained, including: (1) Stress field evolution: analysis of the plastic zone range of surrounding rock, development law, and stress concentration and release of key parts (such as tunnel vault, side wall, working face, rock column between karst cave and tunnel).

[0050] (2) Displacement field evolution: prediction of ground subsidence tank, tunnel convergence deformation, and advanced displacement of soil in front of the working face.

[0051] (3) Seepage field evolution: simulation of the seepage path and pressure distribution change of confined water in the cave-surrounding rock-tunnel system, prediction of potential water inrush channel.

[0052] (4) Comprehensive identification: Through the analysis of the coupling effect of the above physical fields, it is identified whether the water-resisting rock mass will be unstable, seepage unstable (pipe flow) or overall damage under construction disturbance.

[0053] Step S32, dynamic assessment of disaster risk.

[0054] Based on the output results of numerical simulation, the disaster risk is quantitatively and hierarchically dynamically assessed.

[0055] Step S321, establishment of risk evaluation index system.

[0056] Establish the evaluation index with safety factor, plastic zone penetration, displacement threshold, seepage velocity / pressure gradient as the core.

[0057] Step S322, dynamic risk level division.

[0058] The improved risk matrix method is used to combine the probability of disaster occurrence (based on the ease of instability in numerical simulation) and the consequences (based on the scale of disaster and the impact on the project).

[0059] The calculation method of risk value is: ; Wherein, represents the probability of disaster occurrence, which is based on the numerical simulation results and determined by multi-index comprehensive identification; represents the plastic zone penetration rate, i.e. the number of plastic zone units / total number of units; when , ; otherwise, the value is taken in proportion; represents the proportion of displacement of monitoring points exceeding the critical displacement threshold; represents the proportion of seepage velocity or pressure gradient exceeding the critical value; , and are the weights of each index, and , the weights are determined by analytic hierarchy process (AHP) or historical data inversion.

[0060] ; Wherein, represents the severity of disaster consequences; represents the severity of personnel disaster consequences; represents the severity of equipment disaster consequences; represents the severity of construction period disaster consequences; represents the severity of environmental disaster consequences; each index is assigned a grade according to historical cases and expert experience (such as 1-5 points).

[0061] ; wherein, represents the risk value.

[0062] Setting decision threshold , risk level classification is carried out as follows: When , the risk level is low risk, and normal excavation is carried out. When , the risk level is medium risk, and monitoring is strengthened and treatment is prepared. When , the risk level is high risk, and immediate treatment measures are taken and the treatment scheme is implemented. When , the risk level is extremely high risk, and the design scheme is optimized.

[0063] This evaluation is "dynamic", which means that as the position of the shield tunneling face changes, the disaster sources in different sections in front will be re-evaluated in real time (or in stages).

[0064] Step S33, intelligently generating a preliminary treatment scheme.

[0065] The evaluation results are combined with expert knowledge to intelligently output a targeted treatment scheme.

[0066] Step S331, knowledge base and case base construction.

[0067] Knowledge base: store structured treatment method rules, for example: IF (risk level = "high" and disaster type = "high-pressure water-rich fracture zone") THEN (recommended method = "advanced curtain grouting").

[0068] Case base: store a large number of historical engineering cases, including disaster geological conditions, successful treatment schemes (such as grouting pressure, slurry ratio, support parameters) and treatment effects.

[0069] Step S332, scheme generation and comparison.

[0070] The system matches and recommends one or more preliminary treatment schemes from the knowledge base according to the characteristics (type, size, water pressure, and spatial relationship with the tunnel) of the current disaster source and the risk assessment level. The scheme not only indicates the method (such as "full-face advanced deep-hole grouting"), but also gives the range of key process parameters.

[0071] After generating multiple preliminary treatment schemes, a multi-objective decision method is used for comparison; among them, the comprehensive scoring function is as follows: ; wherein, represents the comprehensive score of the th treatment scheme; Indicates the first The weight of each evaluation indicator; Indicates the first The treatment plan was in the first Normalized scores on each evaluation indicator; This indicates the total number of evaluation indicators.

[0072] Ultimately, the overall score is selected. The highest-level plan is considered the optimal initial treatment plan.

[0073] Then, the recommended treatment plan and its parameters are re-introduced into the numerical model for "digital twin" simulation to predict the post-treatment effects (such as whether the plastic zone is reduced, whether the displacement is controlled, and whether the safety factor is met). Through a cycle of "calculation-comparison-optimization," the optimal preliminary treatment plan is selected from multiple alternatives.

[0074] Step S4: Verification and optimization of treatment plan: Based on the fluid-structure interaction similarity theory, physical model tests are used to verify the preliminary treatment plan, and the treatment plan is optimized based on the test results to determine the optimal treatment plan.

[0075] Step S41, Fluid-structure Interaction Similarity Theory and Material Development: The actual engineering prototype (Prototype, P) is scaled down to a laboratory model (Model, M) according to scientific principles.

[0076] Step S411: Based on the fluid-structure interaction similarity theory, calculate the target mechanical parameters (density, compressive strength, elastic modulus, permeability coefficient, etc.) of the ideal similar material required for the model test.

[0077] Step S412, Development of similar materials: The CBCS-type fluid-structure interaction similar material is used, and its composition includes: Aggregates: standard sand, calcium carbonate, iron powder; Cementitious agents: white cement, chlorinated paraffin; Modifier: Silicone oil; Mixing agent: water.

[0078] Through rigorous proportioning tests (such as uniaxial compression, Brazilian splitting, triaxial compression, and permeation tests), the final proportions are determined so that the physical and mechanical parameters (density, compressive / tensile strength, cohesion, internal friction angle, elastic modulus, and permeability coefficient) of the prepared similar materials strictly meet the target values ​​calculated in step S41, thereby ensuring a high degree of similarity between the model and the prototype in terms of mechanics and permeation behavior.

[0079] Step S42: Construct a physical model test system to simulate the shield tunneling and treatment process in water-rich karst strata.

[0080] Step S421, three-dimensional visual model box.

[0081] The size is 3m (length) x 2m (width) x 2m (height), and multiple visual observation windows are arranged on the side surface, as shown in Figure 4 .

[0082] The shield tunneling machine excavation entrance is reserved, and the boundary distance from the tunnel hole diameter is greater than 3 times to eliminate the boundary effect.

[0083] Step S422, composite EPB shield tunneling simulation system: The system has the functions of cutter head rotation tunneling, hydraulic jacking, soil pressure bin pressure monitoring, and double auger soil removal.

[0084] Core features: Its tunneling parameters (thrust, torque, speed, and tunneling speed) can be accurately set and automatically controlled according to the similarity ratio, and can be recorded in real time.

[0085] Step S423, real form disaster source simulation: Based on the real cave form determined in step S2, a 3D printing technology is used to print a disaster source model similar in geometry to the prototype.

[0086] This method solves the problem that traditional pre-embedding methods cannot form complex and irregularly shaped caves and are difficult to remove.

[0087] Step S424, water pressure intelligent regulation and multi-element information monitoring system: Water pressure regulation: Connect the water pressure loading device with the disaster source model to accurately apply and maintain a constant water pressure similar to the prototype water pressure.

[0088] Multi-element data monitoring: Internal physical field: Micro soil pressure cells, osmotic pressure gauges, waterproof strain bricks, and grating multi-point displacement meters are buried at key positions in the model body (such as the tunnel periphery, the front of the tunnel face, and the rock pillar between the cave and the tunnel) to monitor the evolution of stress, osmotic pressure, strain, and displacement in real time.

[0089] Ground subsidence: A three-dimensional laser scanner is used to scan the model surface at high frequency and non-contact, and a subsidence cloud chart of the entire surface during construction is constructed, as shown in Figure 5 .

[0090] Construction parameters: Automatically record all tunneling parameters of the shield tunneling machine.

[0091] Macro phenomena: Record the state of soil removal (such as changes in water content) and the process of sudden gushing water.

[0092] Step S43, implementation, verification, and optimization of treatment scheme.

[0093] Step S431, test implementation: First, according to the preliminary treatment scheme, the model test is accurately implemented.

[0094] Then, the shield machine is started and the tunneling is carried out according to the set parameters until the disaster source influence area is passed.

[0095] Step S432, scheme verification: Under the premise of applying the treatment scheme, the tunneling is carried out. By comparing the monitoring data before and after the treatment, the reliability of the scheme is verified.

[0096] Verification standard: After treatment, whether the displacement of surrounding rock and ground subsidence is effectively controlled within the safety threshold.

[0097] After treatment, whether the plastic zone of surrounding rock is not damaged.

[0098] After treatment, whether the osmotic pressure is not sharply increased and no gushing water disaster occurs.

[0099] After treatment, whether the tunneling parameters of the shield machine are stable.

[0100] Step S433, scheme optimization and optimal scheme determination: Parameter inversion: If the treatment effect does not completely meet the expectation (for example, there is still a large deformation), the shortcomings of the treatment scheme (such as insufficient grouting pressure, insufficient reinforcement range) are analyzed according to the "failure" data monitored in the test.

[0101] Iterative optimization: Based on the analysis results, the parameters of the treatment scheme are adjusted (such as increasing the grouting pressure, optimizing the slurry ratio, increasing the support stiffness), and the test is re-performed.

[0102] Through the iterative process of "test-evaluation-optimization-retest", an optimal treatment scheme and its complete process parameters are finally determined, which are proved to be safe, economic and efficient in the model test, greatly improving the reliability and success rate of the treatment scheme, and providing the most direct technical support for actual construction.

[0103] Step S5, intelligent construction and dynamic control: In the construction process, the optimal treatment scheme is implemented, and the real-time parameters of the construction machinery and the on-site monitoring data are used to dynamically control the treatment process.

[0104] Step S51, intelligent construction execution system.

[0105] The optimal treatment scheme determined in step S4 is accurately implemented in the construction site through digital and intelligent construction equipment and systems.

[0106] Step S511, digitalization and issuance of scheme parameters: The key process parameters of the optimal treatment scheme (such as grouting pressure and flow, shield tunneling parameters, support structure installation position and timing, etc.) are digitized and directly issued to the corresponding intelligent construction equipment control system.

[0107] Step S512, precise execution of intelligent equipment: Intelligent grouting: After receiving the instruction, the grouting equipment automatically grouts according to the preset pressure-flow curve and records the cumulative grouting amount in real time, ensuring accurate filling of the slurry to the design range.

[0108] Shield intelligent tunneling: When the shield machine passes through the disaster source affected area, the control system preferentially adopts the recommended parameters from steps S3 / S4 for automatic tunneling mode, reducing human operation uncertainty.

[0109] Mechanized support installation: Using robotic arms and other automated equipment, the designed positions in the three-dimensional model are quickly and accurately installed with support members such as steel frames and anchor rods.

[0110] Step S52, real-time sensing of multi-element information throughout the process.

[0111] A "digital sensory" system is constructed to comprehensively and real-time sense the construction environment and structure state.

[0112] Step S521, sensing of construction machinery parameters: Real-time collection of shield machine parameters and grouting system data, construction of excavation face stability index and grouting fullness index.

[0113] Based on real-time parameters of the shield machine, the excavation face stability index is constructed As follows: ; Where, represents the measured soil bin pressure; represents the theoretically calculated static soil pressure; represents the measured cutterhead torque; represents the set torque; represents the actual soil volume; represents the theoretical excavation soil volume; , and are weight coefficients reflecting the influence degree of each parameter on stability, which can be obtained through machine learning of historical data.

[0114] Control decision, as follows: When , the state is stable, and the original parameters are used for tunneling; When When the excavation face is at risk of instability, the system automatically prompts to increase the soil chamber pressure or reduce the excavation speed; When , it may occur that the mud cake is broken or over-excavation occurs, the system automatically prompts to reduce the torque or adjust the amount of spoil modifier.

[0115] Based on real-time grouting data, construct grouting fullness index , as follows: ; wherein, represents the cumulative actual grouting amount; represents the theoretical grouting amount calculated according to the shield tail gap; represents the integral of the actual grouting pressure over time (reflecting the effect of pressure duration); represents the integral of the set grouting pressure over time.

[0116] Decision rule: If , it is considered that the grouting is full, and it can continue normally; If , start dynamic control: the system automatically increases the grouting pressure by , wherein is the gain coefficient, until meets the standard.

[0117] Step S522, geological and structural state perception (“digital nerve”): Internal perception: In the rock mass around and in front of the tunnel, the advanced geological prediction system is laid or used to obtain in real time: microseismic / acoustic emission signals to monitor the development of rock mass fracture; deep displacement (through fixed inclinometer or optical fiber) to monitor stratum slip; pore water pressure change to monitor seepage field disturbance.

[0118] Surface perception: In the tunnel, convergence meters and total stations are used to automatically monitor tunnel convergence deformation; on the ground, measurement robots, GB-InSAR or three-dimensional laser scanners are used to continuously monitor changes in the ground subsidence field.

[0119] Step S53, dynamic control based on data fusion.

[0120] Compare the real-time perceived data with the prediction model and control threshold to achieve dynamic optimization and adjustment of construction parameters.

[0121] Step S531, data fusion and state diagnosis: ​The multi-element real-time data obtained in step S52 is fused and analyzed in a unified data platform. Through real-time comparison with the safety threshold predicted by numerical simulation in step S3, the current construction state is diagnosed.

[0122] Step S532, intelligent feedback and parameter regulation: Based on the above diagnosis results, the system generates regulation instructions or issues warnings to the operators, and gives regulation suggestions.

[0123] Shield parameter regulation: If the diagnosis finds that the front rock mass is weak and deforms greatly, the system can automatically or prompt the operator to appropriately reduce the tunneling speed and increase the soil chamber pressure to stabilize the excavation face.

[0124] Treatment parameter strengthening: If the diagnosis finds that the deformation of the grouting reinforcement area is still beyond expectation, the system can instruct the intelligent grouting system to dynamically increase the grouting pressure or perform compensatory grouting in a specific area.

[0125] Support dynamic strengthening: If the initial support stress or deformation is close to the warning value, the system will prompt to strengthen the support parameters in the next cycle, such as thickening the sprayed concrete and increasing the steel frame spacing.

[0126] Step S533, risk warning and emergency intervention: A multi-level warning mechanism is established. When the risk level reaches "orange" or "red" after multi-element data fusion, the system not only adjusts the parameters, but also automatically triggers sound and light alarms, and can set up an emergency shutdown interlock (such as automatically stopping the shield tunneling) to force manual inspection and decision-making to prevent accidents from expanding.

[0127] Step S54, real-time mapping of the digital twin platform.

[0128] The entire process of step S52 and step S53 is visually mapped in real time on the digital twin platform. The data collected on site drives the three-dimensional model in real time, enabling managers to intuitively see the real-time state and interaction of the shield machine, monitoring points, and disaster sources in the virtual space, providing the most intuitive decision support for dynamic regulation.

[0129] Based on the above process, the construction process is no longer a static execution of a predetermined plan, but a "self-adaptive system" that can perceive the environment, diagnose the state, predict risks, and optimize itself, ultimately achieving fine and intelligent dynamic treatment of complex stratum disaster sources.

[0130] Step S6, effect evaluation and model updating: After treatment, the treatment effect is evaluated through post-data, and the effect data is fed back to the three-dimensional geological model and knowledge base to realize dynamic updating of the model and self-learning of the decision-making system.

[0131] Step S61, multi-dimensional post-data collection and effect quantitative evaluation.

[0132] Through the post-construction collected "posterior data", the treatment scheme effect performed in step S5 is objectively and quantitatively evaluated.

[0133] Step S611, posterior data collection.

[0134] (1) Direct exposure verification includes: drilling core verification and borehole television / optical imaging.

[0135] Drilling core verification: verification drilling is performed at the key positions in the grouting reinforcement area or after passing through the disaster source. By analyzing the extracted core, the distribution and filling condition of the slurry vein are observed, and the wave velocity value and uniaxial compressive strength of the core are tested and compared with the data before treatment.

[0136] Borehole television / optical imaging: borehole camera equipment is lowered into the verification hole to visually inspect the integrity of the hole wall, the filling condition of the karst cave and the shape of the slurry solid.

[0137] (2) Non-destructive testing verification includes: cross-hole CT re-measurement and seismic wave method / surface wave method.

[0138] Cross-hole CT re-measurement: elastic wave CT or resistivity CT test is performed again in the treatment area to generate wave velocity / resistivity cloud chart after treatment. By comparing the cloud charts before and after treatment, it is directly shown whether the low-velocity anomaly area / high-conductivity anomaly area disappears, shrinks or weakens in strength.

[0139] Seismic wave method / surface wave method: such methods are used to evaluate the overall uniformity and strength improvement of the surrounding rock after treatment in a large range on the surface above the tunnel or in the hole.

[0140] (3) Long-term stability monitoring: continue to collect monitoring data such as tunnel convergence, ground subsidence and support structure stress for a period of time to confirm whether the deformation and stress have tended to be stable and far below the control value.

[0141] Step S612, effect quantitative evaluation.

[0142] An effect evaluation index system based on posterior data is established, which specifically includes: (1) Filling degree index: based on core and borehole television, the filling percentage of slurry to cavity / crack is evaluated.

[0143] (2) Strength improvement rate: the ratio of the rock mass strength / wave velocity value after treatment to the value before treatment.

[0144] (3) Integrity index: area reduction rate or wave velocity improvement rate of abnormal area in CT image after treatment.

[0145] (4) Stability index: deformation convergence rate and final stable value after treatment.

[0146] According to the above indexes, the comprehensive effect score of treatment is given a quantitative conclusion of "excellent, qualified, unqualified". Among them, the quantitative scoring formula of treatment effect is: ; Among them, represents the comprehensive effect score; and respectively represent the elastic wave velocity before and after treatment; and respectively represent the area of low-speed abnormal area in the CT image before and after treatment; represents the final displacement value after treatment; represents the allowable displacement value; , and respectively represent the weight of strength, integrity and stability.

[0147] Step S62, dynamic updating of three-dimensional geological model and numerical analysis model.

[0148] The posterior data is taken as "true data" to correct the previous model, so as to approach the reality infinitely.

[0149] Step S621, geological property model updating.

[0150] Disaster source model updating: update the disaster source model attribute after treatment from "untreated" to "slurry filling" or "reinforcement", and correct its geometric boundary and physical and mechanical parameters according to the verification data (such as CT retest results).

[0151] Rock mass parameter partition updating: according to the strength improvement rate and wave velocity improvement rate, the rock mass quality grade and mechanical parameters of the reinforced area are partitioned and corrected in the three-dimensional geological model.

[0152] Step S622, numerical analysis model correction.

[0153] The updated three-dimensional geological model is imported into the numerical simulation software again, and the constitutive model parameters (such as the mechanical behavior of grouting improvement body) in the numerical model are fine-tuned (parameter inversion) with the actual deformation and stress data monitored after treatment as the target, so that the results of numerical simulation are highly consistent with the actual situation.

[0154] After calibration, the reliability and accuracy of the model in predicting the risk of future construction stage (such as the next interval) will be significantly improved.

[0155] Step S63, self-learning and evolution of knowledge base and case base.

[0156] Each successful treatment experience is converted into system storage data as historical experience data.

[0157] Step S631, case structuring and storage.

[0158] A new successful treatment case is created, and its data structure includes: (1) Problem description: geological conditions before treatment, characteristics of disaster sources, risk assessment level, etc.

[0159] (2) Treatment scheme: all details of the optimal treatment scheme finally adopted (construction method, equipment, materials, all key process parameters).

[0160] (3) Effect data: all post-treatment data collected in step S61 and effect evaluation conclusion.

[0161] (4) Context information: project name, stratum type, hydrological conditions, etc.

[0162] Step S632, knowledge rule optimization and generation.

[0163] Verify and strengthen existing rules: if the current successful case matches a recommended rule in the knowledge base, the confidence of the rule will be improved.

[0164] Generate new rules: if an innovative or highly optimized treatment scheme is adopted and excellent results are achieved, the system can automatically or semi-automatically generate a new treatment rule.

[0165] Failure lesson record: if the treatment effect is "unqualified", the case will also be recorded in the database, and the reasons will be analyzed to issue a warning in the future to avoid repeating the same mistake.

[0166] Among them, the confidence (CF) is dynamically updated as follows: ; Among them, and respectively represent the updated and updated confidence; represents the normalized value (mapped to 0-1) of the effect score of this case ; represents the learning rate ( ), which controls the speed of confidence update.

[0167] Based on this, successful cases will improve the confidence of the rules used, and failure cases will reduce the confidence, thereby realizing the self-learning and evolution of the system.

[0168] Therefore, the application adopts the above-mentioned fine treatment decision method of complex stratum disaster source fusion multi-source data, realizes the fine management and control of complex stratum disaster source through multi-source data acquisition, fusion modeling, model test verification, risk assessment and treatment decision.

[0169] It should be pointed out finally that the above embodiments are only used to illustrate the technical solutions of the present application but not to limit it, and although the present application has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that the technical solutions of the present application can still be modified or replaced equivalently, and these modifications or equivalent replacements should not make the modified technical solutions deviate from the spirit and scope of the technical solutions of the present application.

Claims

1. A refined treatment decision-making method based on the fusion of multi-source data on hazard sources in complex geological formations, characterized in that, Includes the following steps: Step S1: Integrate multiple geophysical exploration methods and geological drilling to collect multi-source geological data of the construction area; Step S2: Perform data-level, feature-level, and decision-level fusion on the acquired multi-source geological data to construct a three-dimensional geological model that reflects the stratigraphic structure and spatial distribution of disaster sources; Step S3: Based on the three-dimensional geological model, the response of the surrounding rock and the disaster source under construction disturbance is analyzed through numerical simulation, and the disaster risk level is dynamically assessed. By combining the knowledge base and case library, a preliminary treatment plan is generated; Step S4: Based on the fluid-structure interaction similarity theory, a physical model test is used to verify the preliminary treatment plan, and the treatment plan is optimized based on the test results to determine the optimal treatment plan. Step S5: During the construction process, implement the optimal treatment plan and use the real-time parameters of the construction machinery and on-site monitoring data to dynamically control the treatment process. Step S6: After treatment, the treatment effect is evaluated through posterior data, and the effect data is fed back to the three-dimensional geological model and knowledge base to realize the dynamic updating of the model and the self-learning of the decision-making system.

2. The refined treatment decision-making method for complex geological disaster sources by fusing multi-source data as described in claim 1, characterized in that, In step S1, multi-scale and multi-dimensional data acquisition of complex geological hazard sources is achieved through integrated air-space-ground-hole exploration technology. The specific process is as follows: Step S11: In the early stage of the project, the transient electromagnetic method in the water or the multichannel shallow seismic method in the sea is used to conduct a range scan and exploration of the project site to macroscopically understand the distribution of karst development sections and preliminarily divide the project area into different sections of karst micro-development, moderate development and strong development. Step S12: Based on the division results, drill holes are strategically deployed in different sections, and appropriate combinations of geophysical exploration methods are selected, specifically including: (1) Karst areas with strong and moderate development: Elastic wave cross-hole CT as the core method; (2) Karst micro-development zone: supplemented by tube wave detection method and multi-frequency borehole sonar; Step S13: During the multi-source data acquisition process, control geological boreholes are set up at all geophysical exploration points to collect and record rock cores, and to conduct rock and mineral composition analysis and physical and mechanical tests to obtain the real and in-situ parameters of the rock and soil mass, providing calibration points for geophysical interpretation.

3. The refined treatment decision-making method for complex geological disaster sources by fusing multi-source data as described in claim 1, characterized in that, In step S2, the acquired multi-source heterogeneous data are fused to construct a three-dimensional geological model. The specific process is as follows: Step S21: Data-level fusion; The raw data or preliminary processing results of different geophysical exploration methods are aligned and calibrated in a unified three-dimensional spatial coordinate system. Step S22: Feature-level fusion; First, characteristic information of disaster sources is extracted from various types of data, specifically including: extracting low-velocity anomaly zones and their wave velocity values ​​from elastic wave CT data; extracting lithological change surfaces, RQD indices, and the actual locations of exposed karst caves from borehole data; and extracting three-dimensional point clouds of karst cave walls from borehole sonar data. Then, the acquired feature information is combined into a comprehensive feature vector to jointly describe the same geological unit or disaster source; Step S23: Decision-level fusion and 3D geological model construction; Step S231, Basic geological model construction: Based on borehole data and geological profile maps, create basic geological models of topographic surfaces, stratigraphic units, and weathering surfaces using automatic or semi-automatic rock and soil modeling tools; Step S232, Disaster Source Model Construction: Import the three-dimensional numerical simulation analysis model of wave velocity obtained by elastic wave CT inversion into the geological three-dimensional exploration and design system; the system automatically delineates the development range of karst caves in three-dimensional space by identifying wave velocity anomaly thresholds; using this range as a constraint, construct the three-dimensional model of the disaster source using discrete point volume building or karst cave professional modeling tools; Step S233, Model Integration: Integrate the tunnel design model with the basic geological model and the disaster source model to construct a three-dimensional geological model, which intuitively shows the spatial relative position of the tunnel and the disaster source.

4. The refined treatment decision-making method for complex geological disaster sources by fusing multi-source data according to claim 1, characterized in that, In step S3, using the constructed three-dimensional geological model, numerical simulation technology is employed to simulate and predict the dynamic impact of tunnel construction on the surrounding rock and hazard sources. The specific process is as follows: Step S311: Simulation model establishment; First, the constructed three-dimensional geological model is directly imported into the numerical simulation software; different geological bodies in the model are automatically assigned the corresponding rock and soil mechanical parameters obtained through experiments; Then, the entire shield tunneling process was simulated, including: the disturbance of the soil in front by the rotating cutterhead, the friction between the shield shell and the surrounding rock, the formation of the shield tail void, and the pressure and filling effect of synchronous grouting. Finally, a fluid-structure interaction analysis process is defined in the model to simulate the interaction between groundwater and soil in water-rich karst strata. Step S312: After calculation, obtain the response information of the surrounding rock and disaster source throughout the entire construction process, and perform multi-physics response analysis, specifically including: (1) Stress field evolution: Analyze the range and development law of the plastic zone of the surrounding rock, as well as the stress concentration and release in key parts; (2) Displacement field evolution: predicting surface settlement troughs, tunnel convergence deformation, and advance displacement of soil in front of the tunnel face; (3) Evolution of seepage field: Simulate the seepage path and pressure distribution changes of confined water in the cave-surrounding rock-tunnel system to predict potential water inrush channels; (4) Comprehensive identification: By analyzing the coupling effect of the physical field, it can be identified whether the water-retaining rock mass will become unstable, seepage unstable or completely destroyed under construction disturbance.

5. The refined treatment decision-making method for complex geological disaster sources by fusing multi-source data according to claim 4, characterized in that, In step S3, based on the output results of the numerical simulation, a quantitative and hierarchical dynamic assessment of disaster risk is performed. The specific process is as follows: Evaluation indicators were established with safety factor, plastic zone continuity, displacement mutation threshold, seepage velocity, and pressure gradient as the core; an improved risk matrix method was adopted to combine the probability and consequences of disaster occurrence; the risk value was calculated as follows: ； in, It represents the probability of disaster occurrence, based on numerical simulation results and determined through a comprehensive assessment of multiple indicators; This represents the plastic zone penetration rate, i.e., the number of plastic zone elements / the total number of elements; when hour, Otherwise, the value will be taken proportionally. This indicates the proportion of monitoring points whose displacement exceeds the critical displacement threshold; This indicates the proportion of seepage velocity or pressure gradient exceeding a critical value; , and Let be the weight of each indicator, and ; ； in, Indicates the severity of the disaster consequences; Indicates the severity of the disaster's consequences for personnel; Indicates the severity of the consequences of equipment damage; Indicates the severity of the consequences of disasters during the construction period; Indicates the severity of the consequences of an environmental disaster; ； in, Indicates the risk value; Set decision threshold The risk levels are classified as follows: when At that time, the risk level was low, and normal tunneling was underway; when At that time, the risk level was medium, and monitoring was strengthened and preparations were made for treatment. when If the risk level is high, take immediate action and implement the treatment plan. when At that time, the risk level was extremely high, and the design scheme was optimized. Based on the characteristics and risk assessment level of the disaster source, one or more preliminary treatment plans are matched and recommended from the knowledge base; after generating multiple preliminary treatment plans, a multi-objective decision-making method is used for comparison and selection; the comprehensive scoring function is as follows: ； in, Indicates the first The overall score of each treatment plan; Indicates the first The weight of each evaluation indicator; Indicates the first The treatment plan was in the first Normalized scores on each evaluation indicator; Indicates the total number of evaluation indicators; Ultimately, the overall score is selected. The highest-level plan is considered the optimal initial treatment plan. Then, the recommended treatment plan and its parameters are substituted back into the numerical model for digital twin simulation to predict the effect after treatment; through a cycle of calculation-comparison-optimization, the optimal preliminary treatment plan is selected from multiple alternative plans.

6. The refined treatment decision-making method for complex geological disaster sources by fusing multi-source data according to claim 1, characterized in that, In step S4, the preliminary treatment plan is first verified by physical model experiments based on the fluid-structure interaction similarity theory. The specific process is as follows: Step S411: Based on the fluid-structure interaction similarity theory, calculate the target mechanical parameters of the ideal similar material in the model test; Step S412: Use CBCS type fluid-structure interaction similar material, the composition of which includes: aggregate, binder, modifier and mixing agent; Through proportioning experiments, the final proportions are determined so that the physical and mechanical parameters of the prepared similar materials strictly meet the target mechanical parameter values ​​of the ideal similar materials, ensuring the similarity between the model and the prototype in terms of mechanics and percolation behavior; Then, a physical model test system was constructed to simulate the shield tunneling and treatment process in water-rich karst strata, specifically including: Step S421: Assemble a 3D visualization model box using modular steel structure, with multiple visualization observation windows on the side; reserve an excavation entrance for the tunnel boring machine and ensure that its boundary distance is greater than 3 times the tunnel diameter to eliminate boundary effects; Step S422: Establish a composite EPB shield tunneling simulation system, which has the functions of cutterhead rotation tunneling, hydraulic jacking, earth pressure chamber pressure monitoring, and double auger soil removal; its tunneling parameters are set and automatically controlled according to the similarity ratio, and are collected and recorded in real time. Step S423: Based on the actual karst cave morphology determined by the three-dimensional geological model, use 3D printing technology to print a disaster source model that is geometrically similar to the prototype; Step S424: Construct a water pressure intelligent control and multi-dimensional information monitoring system. Connect the water pressure loading device to the disaster source model to apply and maintain a constant water pressure similar to the prototype water pressure. Embed miniature earth pressure cells, piezometers, waterproof strain bricks, and grating multi-point displacement gauges in the model body to monitor the evolution of stress, seepage pressure, strain, and displacement in real time. Use a 3D laser scanner to perform high-frequency, non-contact scanning on the model surface to construct a settlement cloud map of the entire surface during construction. Automatically record all tunneling parameters of the tunnel boring machine, record the soil excavation status, and record the process of sudden water inrush. Finally, the implementation, verification, and optimization of the treatment plan specifically include: Step S431, Test Implementation: First, the generated preliminary treatment plan is implemented in a model test; then, the tunnel boring machine is started and tunneled according to the set parameters until it passes through the disaster source's affected area. Step S432, Solution Verification: Under the premise of implementing the treatment plan, tunneling is carried out; the reliability of the plan is verified by comparing the monitoring data before and after the treatment. Step S433, Scheme Optimization and Determination of Optimal Scheme: Parameter inversion: If the treatment effect does not fully meet expectations, the shortcomings of the treatment plan can be analyzed based on the failure data monitored in the experiment. Iterative optimization: Based on the back analysis results, adjust the parameters of the treatment plan and conduct new experiments to verify it.

7. The refined treatment decision-making method for complex geological disaster sources by fusing multi-source data according to claim 1, characterized in that, In step S5, an intelligent construction execution system is established, and the specific process is as follows: Step S511: Digital distribution of solution parameters: The key process parameters of the optimal treatment plan are digitized and directly sent to the corresponding intelligent construction equipment control system. Step S512, Intelligent equipment executes precisely: Intelligent grouting: After receiving the instruction, the grouting equipment automatically performs grouting according to the preset pressure-flow curve and records the cumulative grouting volume in real time to ensure that the grout is accurately filled into the design range; Intelligent tunneling: When the tunnel boring machine passes through the disaster-affected area, the control system prioritizes the use of recommended parameters for automatic tunneling mode, reducing the uncertainty of human operation; Mechanized support installation: Using robotic automated equipment, support components are installed according to the design positions in the three-dimensional model.

8. The refined treatment decision-making method for complex geological disaster sources by fusing multi-source data according to claim 7, characterized in that, A digital sensory system is constructed to perceive the construction environment and structural status in an all-round and real-time manner, collect shield machine parameters and grouting system data in real time, and construct the excavation face stability index and grouting fullness index. Based on real-time parameters of the tunnel boring machine, an excavation face stability index is constructed. As shown below: ； in, This indicates the measured pressure in the earth chamber; This represents the theoretically calculated earth pressure at rest. This indicates the measured torque of the cutter head; Indicates the set torque; This indicates the actual amount of soil excavated; This indicates the theoretical volume of soil to be excavated. , and These are weighting coefficients, reflecting the degree of influence of each parameter on stability; The regulatory decision is as follows: when At that time, the condition was stable, and tunneling continued according to the original parameters; when If the excavation face is at risk of instability, the system will automatically prompt you to increase the pressure in the soil chamber. Or reduce the tunneling speed; when At this time, mud cake formation or over-excavation may occur, and the system will automatically prompt you to reduce torque. Or adjust the dosage of soil amendment; Based on real-time grouting data, a grouting fullness index is constructed. As shown below: ； in, This indicates the cumulative actual grouting volume; This represents the theoretical grouting volume calculated based on the shield tail gap; This represents the integral of the actual grouting pressure over time. This represents the integral of the set grouting pressure over time; Decision-making rules: like If the grouting is complete, it can be continued normally. like Then dynamic control will be activated: the system will adjust according to... Automatically increases grouting pressure ;in, The gain coefficient is until... Meets the standards; Step S522, Geological and structural state perception; Internal sensing: In the rock mass surrounding and ahead of the tunnel, deploy or utilize advanced geological prediction systems to acquire in real time: microseismic / acoustic emission signals to monitor the development of rock mass fractures; deep displacement to monitor stratum slippage; and pore water pressure changes to monitor seepage field disturbances. Surface sensing: Inside the tunnel, convergence meters and total stations are used to automatically monitor tunnel convergence deformation; on the surface, surveying robots, GB-InSAR, or 3D laser scanners are used to continuously monitor changes in the surface settlement field.

9. The refined treatment decision-making method for complex geological disaster sources by fusing multi-source data according to claim 8, characterized in that, By comparing real-time sensed data with prediction models and control thresholds, dynamic optimization and adjustment of construction parameters can be achieved, specifically including: Step S531, Data Fusion and Status Diagnosis: The acquired multi-dimensional real-time data is integrated and analyzed in a unified data platform; the current construction status is diagnosed by comparing it in real time with the safety thresholds predicted by numerical simulation. Step S532: Based on the diagnostic results, the system generates control commands or issues warnings to operators and provides control suggestions: Shield tunneling parameter control: If the diagnosis finds that the rock mass ahead is weak and deformed, the system will automatically or prompt the operator to appropriately reduce the tunneling speed and increase the soil chamber pressure to stabilize the excavation face; Enhanced treatment parameters: If the diagnosis finds that the deformation in the grouting reinforcement area still exceeds expectations, the system instructs the intelligent grouting system to dynamically increase the grouting pressure or perform compensatory grouting in a specific area; Dynamic reinforcement of support: If the initial support stress or deformation is detected to reach the warning value, the system will prompt to strengthen the support parameters in the next cycle; Step S533, Risk Warning and Emergency Intervention: A multi-level early warning mechanism is established. When the risk level is determined to be orange or red after the fusion of multiple data, the system will not only adjust the parameters, but also automatically trigger an audible and visual alarm and set up an emergency shutdown interlock to force manual inspection and decision-making.

10. The refined treatment decision-making method for complex geological disaster sources by fusing multi-source data according to claim 1, characterized in that, In step S6, the effect evaluation and model update are carried out, and the specific process is as follows: Step S61: Objectively and quantitatively evaluate the effectiveness of the treatment plan by using the post-construction data collected. Step S611, posterior data collection; (1) Direct exposure verification includes: core drilling verification and in-hole television / optical imaging; (2) Non-destructive testing verification includes: cross-hole CT re-measurement and seismic wave / surface wave method; (3) Long-term stability monitoring: Continue to collect data on tunnel convergence, surface settlement and support structure stress for a period of time to confirm whether the deformation and stress have stabilized and are less than the control value; Step S612: Establish an effect evaluation index system based on posterior data, specifically including: (1) Filling percentage index: Based on core sampling and in-hole television, evaluate the percentage of voids / cracks filled by the grout; (2) Strength improvement rate: The ratio of rock mass strength / wave velocity value after treatment to the value before treatment; (3) Integrity index: the area reduction rate or wave velocity increase rate of the abnormal area in the CT image after treatment; (4) Stability indicators: deformation convergence rate and final stable value after treatment; Based on the above indicators, a quantitative conclusion of excellent, qualified, or unqualified is given for the overall treatment effect score; the quantitative scoring formula for treatment effect is as follows: ； in, Indicates the overall effect score; and These represent the elastic wave velocities before and after treatment, respectively. and These represent the areas of the low-velocity abnormality region in CT images before and after treatment, respectively. This represents the final displacement value after treatment; Indicates the allowable displacement value; , and These represent the weights for strength, integrity, and stability, respectively. Step S62: Dynamic updating of the three-dimensional geological model and numerical analysis model; Step S621: Update the geological attribute model; Disaster source model update: Update the attributes of the treated disaster source model from untreated to grout-filled or reinforced, and correct the geometric boundaries and physical and mechanical parameters based on the verification data; Rock mass parameter zoning update: Based on the strength improvement rate and wave velocity improvement rate, the rock mass quality grade and mechanical parameters of the reinforced area are zoning correction in the three-dimensional geological model; Step S622: Numerical analysis model correction; The updated 3D geological model was re-imported into the numerical simulation software. Using the actual deformation and stress data monitored after treatment as the target, the constitutive model parameters in the numerical model were inverted to make the numerical simulation results match the actual situation. Step S63: Self-learning and evolution of the knowledge base and case library; Successful treatment experiences are transformed into system-stored data as historical experience data. Step S631: Create a new successful treatment case, whose data structure includes: (1) Problem description: geological conditions, disaster source characteristics, and risk assessment level before treatment; (2) Treatment plan: All details of the optimal treatment plan to be adopted; (3) Results data: All posterior data collected and results evaluation conclusions; (4) Contextual information: Project name, geological type, hydrological conditions; Step S632: Knowledge rule optimization and generation; Validate and strengthen existing rules: If the current successful case highly matches a recommendation rule in the knowledge base, the confidence of that rule will be increased; Generate new rule: If a new optimized treatment plan is adopted and excellent results are achieved, the system will automatically or semi-automatically generate a new treatment rule. Failure Lessons Record: If the treatment is unsatisfactory, the case will be recorded in the database and the reasons will be analyzed, so as to issue a warning when similar situations are encountered in the future.

Images