Rock stability assessment methods, media and equipment based on three-dimensional geological modeling

Abstract

The present invention discloses a rock formation stability assessment method, medium and equipment based on three-dimensional geological modeling, and relates to the technical field of rock formation assessment. The rock formation stability assessment method based on three-dimensional geological modeling mainly includes: obtaining a three-dimensional geological model and obtaining rock formation parameters to be assessed according to data of the area to be assessed, thereby using finite element simulation software to establish a numerical model; according to engineering requirements and data of the area to be assessed, using the numerical model and an improved stability analysis formula to perform numerical simulation to obtain a comprehensive analysis result; comparing the comprehensive analysis result with actual observation data, adjusting the parameters of the three-dimensional geological model and the numerical model, and using finite element simulation software to perform numerical simulation according to the adjusted model to obtain the rock formation stability result of the area to be assessed. Implementing the rock formation stability assessment method, medium and equipment based on three-dimensional geological modeling provided by the present invention can improve the accuracy of rock formation stability assessment.

Description

Technical Field

[0001] The present invention relates to the technical field of rock formation evaluation, and more specifically, to a rock formation stability evaluation method, medium and equipment based on three-dimensional geological modeling. Background Art

[0002] The 3D geological modeling-based assessment technology is relatively mature and is widely used in the fields of engineering geology and geotechnical engineering. The existing technology relies on high-precision geological and physical data and performs 3D geological modeling and numerical simulation analysis through modeling software.

[0003] The above technology can describe the underground structure and rock formation characteristics in detail, simulate the mechanical behavior of rock formations under different stress conditions, and provide a scientific evaluation basis for rock formation stability. However, the existing technical solutions cannot reflect the changes in rock formation characteristics and environmental conditions over time, which may lead to inaccurate stability assessment results, lack of dynamics and temporal and spatial adaptability. Because the assessment results cannot fully reflect the individual differences of rock formations, it is easy to ignore the personalized characteristics of rock formations. In the relevant existing technologies, all parameters are usually regarded as definite values, and the uncertainty and randomness of the parameters are not considered. The uncertainty analysis is insufficient, resulting in inaccurate assessment results. Similarly, due to its poor accuracy, it is difficult to adapt to complex geological conditions and cannot provide a comprehensive risk assessment. Therefore, it is urgent to propose a rock formation stability assessment method, system and electronic equipment based on three-dimensional geological modeling. Summary of the invention

[0004] The purpose of the present invention is to provide a rock formation stability assessment method, medium and equipment based on three-dimensional geological modeling, which can improve the accuracy of rock formation stability assessment.

[0005] The present invention provides a rock formation stability assessment method based on three-dimensional geological modeling, comprising the following steps: S1: acquiring data of a region to be assessed, and obtaining a three-dimensional geological model according to the data of the region to be assessed; S2: acquiring parameters of the rock formation to be assessed according to the three-dimensional geological model, and establishing a numerical model according to the parameters of the rock formation to be assessed using finite element simulation software; S3: acquiring engineering requirements, and performing numerical simulation using a numerical model and an improved stability analysis formula according to the engineering requirements and the data of the region to be assessed, to obtain a comprehensive analysis result; S4: acquiring actual observation data, comparing the comprehensive analysis result with the actual observation data, and adjusting the parameters of the three-dimensional geological model and the numerical model according to the comparison result to obtain an adjusted model; S5: performing numerical simulation using finite element simulation software according to the adjusted model to obtain a rock formation stability result of the region to be assessed.

[0006] Furthermore, the data of the area to be evaluated in the above-mentioned rock formation stability assessment method based on three-dimensional geological modeling include geological exploration data, remote sensing data, environmental data and historical data; the geological exploration data include drilling data, geological maps, seismic reflection data and core samples of the area to be evaluated; the remote sensing data include satellite images and LIDAR data of the area to be evaluated; the environmental data include groundwater level data and meteorological data of the area to be evaluated; the historical data include earthquake record data and landslide record data of the area to be evaluated; the parameters of the rock formation to be evaluated include physical parameters, mechanical parameters, geological structure parameters, hydrogeological parameters and environmental impact parameters; the physical parameters include density, porosity and water content of the rock formation to be evaluated; the mechanical parameters include uniaxial compressive strength, triaxial compressive strength, shear strength, elastic modulus and Poisson's ratio of the rock formation to be evaluated; the geological structure parameters include fault parameters and joint parameters; the hydrogeological parameters include permeability parameters and groundwater pressure parameters; the environmental impact parameters include temperature and rainfall.

[0007] Furthermore, step S1 of the above-mentioned rock formation stability assessment method based on three-dimensional geological modeling specifically includes: S11: acquiring data of the area to be assessed, preprocessing the data of the area to be assessed, and obtaining preprocessed data of the area to be assessed; S12: performing geological construction using modeling software based on the preprocessed data of the area to be assessed, and obtaining a three-dimensional geological model.

[0008] Furthermore, step S2 of the above-mentioned rock formation stability assessment method based on three-dimensional geological modeling specifically includes: S21: obtaining rock formation parameters to be assessed based on the three-dimensional geological model; S22: importing the three-dimensional geological model into finite element simulation software, using structured grids, setting grid division parameters, running grid division tools, and generating a grid of a numerical model; S23: obtaining a numerical model based on the rock formation parameters to be assessed and the grid of the numerical model.

[0009] Furthermore, step S3 of the above-mentioned rock formation stability assessment method based on three-dimensional geological modeling specifically includes: S31: obtaining engineering requirements, and obtaining simulation time range, simulation space range, simulation characteristics and simulation parameters according to the engineering requirements and the data of the area to be assessed; S32: performing numerical simulation using a numerical model according to the simulation time range, simulation space range, simulation characteristics and simulation parameters to obtain simulation process data, and the simulation process data includes stress, strain and displacement during the simulation process; S33: calculating the comprehensive stability score of each time and space position according to the simulation process data using an improved stability analysis formula to obtain a comprehensive analysis result.

[0010] Furthermore, step S33 of the above-mentioned rock formation stability assessment method based on three-dimensional geological modeling specifically includes: according to the simulation process data, using the improved stability analysis formula, calculating the comprehensive stability score of each time and space position, and obtaining a comprehensive analysis result, such as the formula:

[0011]

[0012]

[0013]

[0014]

[0015] in, F ( t , x , y , z ) is the rock formation stability score, which changes with time t and spatial position ( x , y , z )change; w i ( t , x , y , z ) is the i The dynamic weight factor of each feature reflects the change of the feature over time and space; x i ' is the i The standardized feature values; α i and β i is a constant used to adjust the eigenvalue; It is j Dynamic weight factors, including temporal and spatial variables; d j ' is the j The normalized dynamic eigenvalues; θ k It is k The weighting factors of the uncertainty parameters; u k ' is the k The standardized uncertainty parameter value; , , and is a constant used to adjust the dynamic eigenvalues ​​and uncertainty parameters; m , n and p are the number of features, dynamic weight factors, and uncertainty parameters, respectively; w i ( t ) is the time factor, which is used to consider seasonal changes and long-term trends. ai and b i is a constant, ω is the frequency, φ i is the phase angle; w i ( x , y , z ) is a spatial factor, which is used to consider the differences in geographical locations and geological conditions. x 0, y 0, z 0) is the center position, σ i is the scale parameter, N is the number of simulations, f (.) indicates a simulation function. u k It is k uncertainty parameter values, is a random variable.

[0016] Furthermore, step S4 of the above-mentioned rock formation stability assessment method based on three-dimensional geological modeling specifically includes: S41: obtaining actual observation data, comparing the comprehensive analysis results with the actual observation data, and obtaining model accuracy assessment results; S42: adjusting the parameters and structure of the three-dimensional geological model according to the model accuracy assessment results to obtain a new three-dimensional geological model; adjusting the physical and mechanical parameters of the numerical model based on the model accuracy assessment results and the new three-dimensional geological model to obtain adjusted parameters and structural data; S43: performing numerical simulation using finite element simulation software based on the adjusted parameters and structural data and the numerical model to obtain new simulation results, comparing the new simulation results with the actual observation data again, and verifying and adjusting the three-dimensional geological model and the numerical model according to the results of the re-comparison to obtain the adjusted model.

[0017] Furthermore, step S5 of the above-mentioned rock formation stability assessment method based on three-dimensional geological modeling specifically includes: S51: according to the adjusted model, setting the final simulation conditions, confirming the boundary conditions and initial conditions, and synchronously setting the gravity, ground stress, and groundwater pressure loading conditions; using finite element simulation software, performing numerical simulation to obtain simulation results; the simulation results include stress distribution diagrams, strain distribution diagrams, and displacement distribution diagrams; S52: according to the simulation results, obtaining the rock formation stability results of the area to be assessed.

[0018] The present invention also provides a computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor, the steps of the rock formation stability assessment method based on three-dimensional geological modeling are implemented.

[0019] The present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of the above-mentioned rock formation stability assessment method based on three-dimensional geological modeling are implemented.

[0020] The implementation of the rock formation stability assessment method, medium and equipment based on three-dimensional geological modeling provided by the present invention has the following beneficial effects:

[0021] The present invention obtains data of the area to be evaluated and constructs a three-dimensional geological model, obtains parameters of the rock formation to be evaluated based on the model and establishes a numerical model based on the parameters, and uses an improved stability analysis formula to perform stability analysis on the area to be evaluated. The improved stability analysis formula increases the personalized characteristics of the rock formation by introducing special weight factors for the rock formation part, and simultaneously introduces dynamic weight factors to increase dynamics and time-space adaptability, while taking time and space variables into consideration and adding uncertainty parameters, thereby deriving a comprehensive analysis result. The method makes targeted limitations on the assessment of rock formation stability, and different assessment results can be obtained for rock formations in different locations in the same area, thereby greatly improving the accuracy of the model results. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] The present invention will be further described below with reference to the accompanying drawings and embodiments, in which:

[0023] Figure 1 It is a flow chart of a rock formation stability assessment method based on three-dimensional geological modeling provided by the present invention;

[0024] Figure 2 A schematic diagram of the execution steps of the rock formation stability assessment method based on three-dimensional geological modeling provided by the present invention;

[0025] Figure 3 A schematic diagram of a rock formation stability assessment method based on three-dimensional geological modeling provided by the present invention;

[0026] Figure 4 A schematic diagram of an electronic device provided by the present invention;

[0027] Among them, 41 is an electronic device, 411 is a memory, 412 is a communication module, 413 is a bus, and 414 is a processor. DETAILED DESCRIPTION

[0028] In order to have a clearer understanding of the technical features, purposes and effects of the present invention, specific embodiments of the present invention are now described in detail with reference to the accompanying drawings.

[0029] Figure 1 A schematic diagram of a rock formation stability assessment method based on three-dimensional geological modeling in this embodiment is shown. In this embodiment, the rock formation stability assessment method based on three-dimensional geological modeling includes the following steps:

[0030] S1: Acquire data of the area to be evaluated, and obtain a three-dimensional geological model based on the data of the area to be evaluated;

[0031] Specifically, the data of the area to be evaluated include geological exploration data, remote sensing data, environmental data and historical data; the geological exploration data include drilling data, geological maps, seismic reflection data and core samples of the area to be evaluated; the remote sensing data include satellite images and LIDAR data of the area to be evaluated; the environmental data include groundwater level data and meteorological data of the area to be evaluated; the historical data include earthquake record data and landslide record data of the area to be evaluated;

[0032] Specifically, step S1 includes:

[0033] S11: Acquire the data of the area to be evaluated, preprocess the data of the area to be evaluated, and obtain the preprocessed data of the area to be evaluated;

[0034] Specifically, the above preprocessing includes format conversion, denoising, interpolation and standardization. Format conversion is used to convert data from different sources into a unified format. Denoising uses filtering technology to remove noise and outliers in the data. Interpolation is used to interpolate missing data to ensure data integrity. Standardization is used to standardize data to make them in the same dimension and range.

[0035] S12: Based on the pre-processed data of the area to be evaluated, the geological structure is carried out using modeling software to obtain a three-dimensional geological model;

[0036] Specifically, the above geological structure includes stratigraphic structure and fault structure, and the three-dimensional geological model includes a detailed model, and the detailed model is used to reflect the complexity of the underground structure;

[0037] S2: According to the three-dimensional geological model, the parameters of the rock formation to be evaluated are obtained, and according to the parameters of the rock formation to be evaluated, a numerical model is established using finite element simulation software;

[0038] Specifically, the above-mentioned rock formation parameters to be evaluated include physical parameters, mechanical parameters, geological structure parameters, hydrogeological parameters and environmental impact parameters; physical parameters include density, porosity and water content of the rock formation to be evaluated; mechanical parameters include uniaxial compressive strength, triaxial compressive strength, shear strength, elastic modulus and Poisson's ratio of the rock formation to be evaluated; geological structure parameters include fault parameters and joint parameters; hydrogeological parameters include permeability parameters and groundwater pressure parameters; environmental impact parameters include temperature and rainfall;

[0039] Specifically, step S2 includes: S21: obtaining the parameters of the rock formation to be evaluated according to the three-dimensional geological model; S22: importing the three-dimensional geological model into the finite element simulation software, using the structured grid, setting the grid division parameters, running the grid division tool, and generating the grid of the numerical model; S23: obtaining the numerical model according to the parameters of the rock formation to be evaluated and the grid of the numerical model;

[0040] Specifically, the above step S2 can also be implemented by the following steps: importing a three-dimensional geological model and converting it into ANSYS format, creating an analysis project in ANSYS; using a structured grid, setting grid division parameters, running a grid division tool, and generating a grid for a numerical model; in ANSYS, assigning the acquired physical parameters and mechanical parameters to the rock formation to be evaluated; in ANSYS, assigning structural parameters based on the acquired geological structure information; in ANSYS, assigning hydrogeological parameters to the model; in ANSYS, assigning environmental impact parameters to the model; setting boundary conditions and initial conditions; boundary conditions, setting the boundary conditions of the model and determining the specific location and type of the boundary conditions; defining initial conditions in ANSYS, setting initial stress states and initial groundwater level initial conditions; preliminary verification, running preliminary simulations, checking the accuracy of the numerical model, and based on the preliminary results, adjusting the model parameters and grid division to obtain a numerical model;

[0041] S3: Obtain engineering requirements, and perform numerical simulation using numerical models and improved stability analysis formulas based on engineering requirements and data of the area to be evaluated to obtain comprehensive analysis results;

[0042] Specifically, step S3 includes:

[0043] S31: Obtain engineering requirements, and obtain simulation time range, simulation space range, simulation characteristics and simulation parameters according to the engineering requirements and the data of the area to be evaluated;

[0044] S32: performing numerical simulation using a numerical model according to a simulation time range, a simulation space range, simulation characteristics and simulation parameters to obtain simulation process data, where the simulation process data includes stress, strain and displacement during the simulation process;

[0045] S33: According to the simulation process data, the comprehensive stability score of each time and space position is calculated using the improved stability analysis formula to obtain a comprehensive analysis result;

[0046] Specifically, step S33 includes: according to the simulation process data, using the improved stability analysis formula, calculating the comprehensive stability score of each time and space position, and obtaining a comprehensive analysis result, such as the formula:

[0047]

[0048]

[0049]

[0050]

[0051] in, F ( t , x , y , z ) is the rock formation stability score, which changes with time t and spatial position ( x , y , z )change; w i ( t , x , y , z ) is the i The dynamic weight factor of each feature reflects the change of the feature over time and space; x i ' is the i The standardized feature values; α i and β i is a constant used to adjust the eigenvalue; It is j Dynamic weight factors, including temporal and spatial variables; d j ' is the j The normalized dynamic eigenvalues; θ k It is k The weighting factors of the uncertainty parameters; u k ' is the k The standardized uncertainty parameter value; , , and is a constant used to adjust the dynamic eigenvalues ​​and uncertainty parameters; m , n and p are the number of features, dynamic weight factors, and uncertainty parameters, respectively; w i ( t ) is the time factor, which is used to consider seasonal changes and long-term trends. a i and b i is a constant, ω is the frequency,φ i is the phase angle; w i ( x , y , z ) is a spatial factor, which is used to consider the differences in geographical locations and geological conditions. x 0, y 0, z 0) is the center position, σ i is the scale parameter, N is the number of simulations, f (.) indicates a simulation function. u k It is k uncertainty parameter values, is a random variable;

[0052] S4: obtaining actual observation data, comparing the comprehensive analysis results with the actual observation data, and adjusting the parameters of the three-dimensional geological model and the numerical model according to the comparison results to obtain an adjusted model;

[0053] Specifically, step S4 includes:

[0054] S41: Obtain actual observation data, compare the comprehensive analysis results with the actual observation data, and obtain the model accuracy evaluation results;

[0055] Specifically, step S41 includes: collecting actual observation data and historical data, comparing the comprehensive analysis results and numerical simulation results with the actual observation data, calculating the error of the simulation results, and evaluating the accuracy of the model; specifically, there are two comparisons, namely, graphic comparison and numerical comparison; wherein the graphic comparison compares the stress, strain and displacement distribution diagrams of the simulation results with the actual measurement data diagram; for the stress concentration area, compare whether the stress concentration area in the actual observation data is consistent with the simulation result; for the strain distribution, compare whether the strain distribution in the actual observation is consistent with the simulation result; for the displacement difference, compare whether the displacement in the actual observation is consistent with the simulation result; the numerical comparison is to calculate the numerical difference between the simulation result and the actual data; the model accuracy evaluation result is obtained by combining the graphic comparison and the numerical comparison; specifically, the evaluation error is evaluated using the root mean square error, such as the formula:

[0056]

[0057] In the formula, is the actual observed value, is the analog value, is the sample size;

[0058] S42: adjusting the parameters and structure of the three-dimensional geological model according to the model accuracy evaluation result to obtain a new three-dimensional geological model; adjusting the physical and mechanical parameters of the numerical model according to the model accuracy evaluation result and the new three-dimensional geological model to obtain adjusted parameters and structural data;

[0059] Specifically, step S42 includes: making adjustments based on the results obtained from the above verification, identifying the deviation between the numerical simulation results and the actual observation data, adjusting the parameters and structure of the three-dimensional geological model based on the verification results, and adjusting the physical and mechanical parameters of the numerical model based on the verification results; wherein adjusting the parameters and structure of the three-dimensional geological model includes: correcting the geological structure, based on the actual observation data, correcting the faults, joints and other geological structures in the model; adjusting the stratum thickness, based on the actual measurement data, adjusting the stratum thickness in the model; optimizing the grid division, based on the verification results, optimizing the grid division; adjusting the physical and mechanical parameters of the numerical model according to the fir tree verification results, specifically including: physical parameter adjustment, based on the actual measurement data, adjusting the density, porosity and water content parameters; mechanical parameter adjustment, based on the actual observation data, adjusting the uniaxial compressive strength, triaxial uniaxial compressive strength, shear strength, elastic modulus and Poisson's ratio parameters; hydrogeological parameter adjustment, based on the actual groundwater level and permeability data, adjusting the permeability coefficient and groundwater pressure; finally, obtaining the adjusted parameter and structure data;

[0060] S43: performing numerical simulation using finite element simulation software according to the adjusted parameters, structural data and numerical model to obtain new simulation results, comparing the new simulation results with the actual observation data again, and verifying and adjusting the three-dimensional geological model and the numerical model according to the results of the re-comparison to obtain an adjusted model;

[0061] S5: Based on the adjusted model, numerical simulation is performed using finite element simulation software to obtain the rock stability results of the area to be evaluated;

[0062] Specifically, step S5 includes:

[0063] S51: According to the adjusted model, set the final simulation conditions, confirm the boundary conditions and initial conditions, and simultaneously set the gravity, ground stress, and groundwater pressure loading conditions; use finite element simulation software to perform numerical simulation and obtain simulation results; the simulation results include stress distribution diagrams, strain distribution diagrams, and displacement distribution diagrams;

[0064] S52: Obtaining the rock stability result of the area to be evaluated according to the simulation result;

[0065] Specifically, the rock stability results of the above-mentioned area to be evaluated include stress concentration areas, areas with large deformation and potential unstable areas.

[0066] In some embodiments, the above-mentioned rock formation stability assessment method based on three-dimensional geological modeling can also be implemented in the following ways;

[0067] Figure 2 A schematic diagram of the execution steps of the rock formation stability assessment method based on three-dimensional geological modeling provided in an embodiment of the present application. In this embodiment, the rock formation stability assessment method based on three-dimensional geological modeling can be executed by an electronic device, and the electronic device can be a server (such as a cloud server, a server cluster, etc.) or a terminal (such as a personal computer, a laptop, etc.), which is not limited here; in this embodiment of the present application, the rock formation stability assessment method based on three-dimensional geological modeling includes step S100, step S200, step S300, step S400 and step S500;

[0068] Step S100, obtaining data of the area to be evaluated and constructing a three-dimensional geological model; in order to construct a three-dimensional geological model and obtain benchmark data, the electronic device can execute step S100; in this embodiment, the data of the area to be evaluated includes geological exploration data, remote sensing data, environmental data and historical data; the geological exploration data includes drilling data, geological maps, seismic reflection data and core samples of the area to be evaluated; the remote sensing data includes satellite images and LIDAR data of the area to be evaluated; the environmental data includes groundwater level data and meteorological data of the area to be evaluated; the historical data includes seismic record data and landslide record data of the area to be evaluated; wherein the drilling data is used to provide vertical profile information of the rock formation, including lithology, stratigraphic position and thickness data; the geological map is used to show the distribution and relationship of surface geological units; the seismic reflection data is used to identify underground structures, specifically faults and interlayer interfaces; the core samples are used to obtain detailed rock physics and mechanics nature; satellite images are used to provide large-scale surface information; LIDAR data are used to provide high-precision surface height data, targeting areas with complex terrain; groundwater level data specifically include: groundwater level and its changes have an impact on rock stability; meteorological data specifically include rainfall and temperature, which directly affect the mechanical properties of rock formations and groundwater level changes; earthquake record data specifically record the frequency and intensity of historical earthquake activities; landslide record data specifically record the occurrence of historical landslide events and their impacts; among them, it should be explained that drilling data: core samples are obtained through drilling projects to record the depth, thickness and lithology of rock formations; geological maps and seismic reflection data: obtained from geological surveys and seismic exploration projects; remote sensing data and LIDAR data: obtained from satellite remote sensing and aerial surveys; environmental data: obtained from meteorological and hydrological monitoring stations; historical geological activity data: obtained from earthquake and landslide record databases;

[0069] Step S200, obtaining the parameters of the rock formation to be evaluated based on the model, and establishing a numerical model based on the parameters; in order to realize numerical simulation, the electronic device may execute step S200;

[0070] Step S200 specifically includes: preprocessing, which includes format conversion, denoising, interpolation and standardization processing; wherein format conversion: converting data from different sources into a unified format; denoising: using filtering technology to remove noise and outliers in the data; interpolation: interpolating missing data to ensure data integrity; standardization: standardizing the data to make them in the same dimension and range;

[0071] Constructing a three-dimensional geological model, inputting the pre-processed data into a modeling software, and constructing a three-dimensional geological model based on the input data, wherein the construction includes strata, faults and other geological structures, and refining the preliminary model, wherein the refined model is used to reflect the complexity of the underground structure;

[0072] Among them, by obtaining multi-source data of the area to be evaluated, building a high-precision three-dimensional geological model, and obtaining detailed rock formation parameters, a solid data foundation can be provided for rock formation stability assessment; this step not only ensures the scientificity and accuracy of the assessment, but also lays the foundation for subsequent numerical simulation and stability analysis;

[0073] It should be noted that the parameters of the rock formation to be evaluated include physical parameters, mechanical parameters, geological structure parameters, hydrogeological parameters and environmental impact parameters; physical parameters include density, porosity and water content of the rock formation to be evaluated; mechanical parameters include uniaxial compressive strength, triaxial compressive strength, shear strength, elastic modulus and Poisson's ratio of the rock formation to be evaluated; geological structure parameters include fault parameters and joint parameters; hydrogeological parameters include permeability parameters and groundwater pressure parameters; environmental impact parameters include temperature and rainfall;

[0074] It should be noted that ANSYS is used for numerical simulation, a three-dimensional geological model is imported and converted into ANSYS format, and an analysis project is created in ANSYS; a structured grid is used, grid division parameters are set, and a grid division tool is run to generate a grid for the numerical model; in ANSYS, the physical parameters and mechanical parameters are assigned to the rock layer to be evaluated; in ANSYS, structural parameters are assigned based on the obtained geological structure information; in ANSYS, hydrogeological parameters are assigned to the model; in ANSYS, environmental impact parameters are assigned to the model; boundary conditions and initial conditions are set; boundary conditions, the boundary conditions of the model are set and the specific location and type of boundary conditions are determined; initial conditions are defined in ANSYS, and initial stress states and initial groundwater levels are set; preliminary verification, preliminary simulation is run to check the accuracy of the numerical model, and based on the preliminary results, model parameters and grid division are adjusted;

[0075] Step S300, performing stability analysis on the area to be evaluated, wherein the stability analysis adopts an improved stability analysis formula and outputs a comprehensive analysis result; in order to construct the improved stability analysis formula and perform stability analysis processing, the electronic device may execute step S300;

[0076] Among them, the improved stability analysis formula is constructed by introducing special weight factors of the rock formation and dynamic weight factors simultaneously. The dynamic weight factors include time and space variables and uncertainty parameters are added to the analysis formula. Specifically,

[0077]

[0078] in, F ( t , x , y , z ) is the rock formation stability score, which changes with time t and spatial position ( x , y , z )change; w i ( t , x , y , z ) is the i Dynamic weight factors of features; x i ' is the i The standardized feature values; α i and β i is a constant used to adjust the eigenvalue; It is j Dynamic weight factors, including temporal and spatial variables; d j ' is the j The normalized dynamic eigenvalues; θ k It is k The weighting factors of the uncertainty parameters; u k ' is the k The standardized uncertainty parameter value; , , and is a constant used to adjust the dynamic eigenvalues ​​and uncertainty parameters; m , n and p are the number of features, dynamic weight factors, and uncertainty parameters, respectively;

[0079] It should be noted that the score is determined based on the threshold, and the specific threshold standards are set as follows: low risk: F < 1.0; medium risk: 1.0 ≤ F < 2.0; high risk: F ≥ 2.0;

[0080] Among them, the improved stability analysis formula includes: introducing dynamic weight factor and Reflects the changes of characteristics over time and space. The time factor is used to consider seasonal changes and long-term trends. Specifically:

[0081]

[0082] In the formula, a i and b i is a constant, ω is the frequency, φ i is the phase angle;

[0083] Spatial factors are used to account for differences in geographical locations and geological conditions, specifically:

[0084]

[0085] In the formula, a i is a constant, ( x 0, y 0, z 0) is the center position, σ i is the scale parameter;

[0086] Introducing uncertainty parameters θ k and u k , Monte Carlo simulation method is used to evaluate the uncertainty of model results, specifically:

[0087]

[0088] In the formula, is the number of simulations, is a random variable;

[0089] The stability analysis of the area to be assessed includes setting the time range and space range of the simulation based on the geological environment to be assessed and the engineering requirements, and simultaneously confirming the characteristics and parameters of the simulation; running the simulation, inputting all parameters and conditions in ANSYS, running the numerical simulation, and simultaneously recording the specific data generated during the simulation, where the data includes stress, strain and displacement during the simulation; using the improved stability analysis formula, calculating and outputting the comprehensive stability score for each time and space position;

[0090] Exemplarily, a set of experimental data is substituted for verification;

[0091] Physical parameters: Density ( ρ ):2.5g / cm 3 , porosity ( φ ): 15%, moisture content ( w ):10%;

[0092] Mechanical parameters: Uniaxial compressive strength ( σ c ):50MPa, shear strength ( τ ):20MPa, elastic modulus ( E c ):30GPa, Poisson's ratio ( ν ):0.25;

[0093] Geological structural parameters: Fault influence ( F ): There is a fault, the fault dip angle is 45°, and the joint spacing is 1.5m;

[0094] Hydrogeological parameters: Permeability coefficient ( k ):1×10 -7 m / s , groundwater pressure ( u ):05MPa;

[0095] Environmental impact parameters: Temperature ( T ):20℃, rainfall( P ):100mm / year;

[0096] Dynamic weight factor: time factor ( t ): consider seasonal changes, spatial factors (x, y, z): consider geographical location changes;

[0097] Normalize the above values ​​to: ρ ′=0.8, φ ′=0.7, w ′=0.6, σ c ′=0.9, τ ′=0.5, E′=0.85, ν ′=0.4, F ′=0.75, k ′=0.65, u ′=0.55, T ′=0.7, P ′=0.6;

[0098] Dynamic weight factor normalization, specifically: w ρ ( t, x, y, z )=0.2, w φ ( t, x, y, z )=0.15, w w ( t, x, y, z )=0.1, w σc ( t, x, y, z )=0.25; w τ ( t, x, y, z )=0.05, w E ( t, x, y, z )=0.2, w ν ( t, x, y, z )=0.05, w F ( t, x, y, z )=0.2; w k ( t, x, y, z )=0.1, w u ( t, x, y, z )=0.1, w T ( t, x, y, z )=0.15, w P ( t, x, y, z )=0.1;

[0099] Other parameter settings are: α i =1, β i =2, θ k =0.1, α uk =1, β uk =2;

[0100] Substitute the above data and use the improved stability analysis formula in step S300 to perform calculations, specifically:

[0101]

[0102] Calculate the contribution of each term, specifically:

[0103]

[0104] Further calculations yield:

[0105]

[0106] Finally, we conclude:

[0107]

[0108] Based on the scoring results, it was concluded that the stability score of the rock formation was 1.425, which was in the medium risk range and required enhanced monitoring and management to ensure that the rock formation remained stable during future engineering processes;

[0109] Step S400, verifying the analysis result, and making adjustments based on the verification result, wherein the adjustments include adjusting the three-dimensional geological model and the numerical model parameters; in order to verify and confirm the analysis result, the electronic device may execute step S400;

[0110] Step S400 specifically includes: verifying the results, including collecting actual observation data and historical data, comparing the numerical simulation results with the actual data, calculating the error of the simulation results, and evaluating the accuracy of the model; comparing the numerical simulation results with the actual observation data to determine the accuracy of the model, a total of two comparisons, namely: graphic comparison, comparing the stress, strain and displacement distribution diagram of the simulation results with the actual measurement data diagram; numerical comparison, calculating the numerical difference between the simulation results and the actual data; evaluating the error, using the root mean square error for evaluation, specifically:

[0111]

[0112] In the formula, is the actual observed value, is the analog value, is the sample size;

[0113] Adjustments are made based on the results of the verification, the deviations between the numerical simulation results and the actual observation data are identified, the parameters and structure of the geological model are adjusted based on the verification results, and the physical and mechanical parameters of the numerical model are adjusted based on the verification results; for stress concentration areas, whether the stress concentration areas in the actual observation data are consistent with the simulation results is determined through comparative analysis and error evaluation; for strain distribution, whether the strain distribution in the actual observation is consistent with the simulation results is determined; for displacement differences, whether the displacement in the actual observation is consistent with the simulation results is determined; adjust the parameters and structure of the three-dimensional geological model, correct the geological structure, and correct the faults, joints and other geological structures in the model based on the actual observation data; adjust the stratum thickness, and adjust the stratum thickness in the model based on the actual measurement data. Layer thickness; Optimize mesh division, based on the verification results, optimize mesh division; According to the verification results, adjust the physical and mechanical parameters of the numerical model, including: physical parameter adjustment, based on the actual measurement data, adjust the density, porosity and water content parameters; mechanical parameter adjustment, based on the actual observation data, adjust the uniaxial compressive strength, triaxial uniaxial compressive strength, shear strength, elastic modulus and Poisson's ratio parameters; hydrogeological parameter adjustment, based on the actual data of groundwater level and permeability, adjust the permeability coefficient and groundwater pressure; Update the adjusted parameters and structure to the numerical model, re-run the simulation, re-run the numerical simulation in ANSYS, compare the new results, compare the new simulation results with the actual observation data again, and verify the adjustment effect;

[0114] Step S500, confirming the rock formation stability result of the area to be evaluated based on the adjusted model; in order to confirm the stability result of the area to be evaluated, the electronic device may execute step S500;

[0115] Step S500 specifically includes: setting the final simulation conditions, confirming the boundary conditions and initial condition settings, and synchronously setting the gravity, ground stress, and groundwater pressure loading conditions; running the final numerical simulation in ANSYS, and recording the data generated by the process; analyzing the simulation results, wherein the results include stress distribution diagrams, strain distribution diagrams, and displacement distribution diagrams; based on the analysis of the simulation results, identifying stress concentration areas, areas with large deformation, and potential unstable areas.

[0116] This embodiment obtains data of the area to be evaluated, constructs a three-dimensional geological model, obtains parameters of the rock formation to be evaluated based on the model, and establishes a numerical model based on the parameters to perform stability analysis on the area to be evaluated. The stability analysis adopts an improved stability analysis formula. The improved stability analysis formula introduces special weight factors of the rock formation part and simultaneously introduces dynamic weight factors. The dynamic weight factors include time and space variables. Uncertainty parameters are added to the analysis formula, and comprehensive analysis results are output. The analysis results are verified, and adjustments are made based on the verification results. The stability results of the rock formations in the area to be evaluated are confirmed based on the adjusted model. This method is based on special weight factors of the rock formation part, combined with time and space variables, and introduces uncertainty parameters to improve the accuracy of the model results. It is not a simple increase or decrease of formula calculation items, but a targeted limitation is made for the rock formation stability assessment, even for rock formations in the same region but different locations. There are different assessment results.

[0117] See also Figure 3 , Figure 3 A schematic diagram of a rock formation stability assessment system based on three-dimensional geological modeling is provided for an embodiment of the present application; in this embodiment, the rock formation stability assessment system based on three-dimensional geological modeling may include: a modeling unit, used to obtain data of the area to be assessed and construct a three-dimensional geological model; a numerical simulation unit, used to obtain parameters of the rock formation to be assessed based on the model, and to establish a numerical model based on the parameters; an analysis unit, used to perform stability analysis on the area to be assessed, wherein the stability analysis adopts an improved stability analysis formula and outputs a comprehensive analysis result; a verification unit, used to verify the analysis result and make adjustments based on the verification result, wherein the adjustment includes adjusting the three-dimensional geological model and the numerical model parameters; a confirmation unit, used to confirm the rock formation stability result of the area to be assessed based on the adjusted model.

[0118] See also Figure 4 , Figure 4 A schematic diagram of an electronic device 41 provided for an embodiment of the present application, in this embodiment, the electronic device 41 can be a terminal, such as a tablet computer, a personal computer, etc.; the electronic device 41 can also be a server, such as a cloud server, a server cluster, etc., which is not limited here; exemplarily, the electronic device 41 may include: a communication module 412 connected to the outside world through a network, one or more processors 414 for executing program instructions, a bus 413 and different forms of memory 411, such as a disk, ROM or RAM, or any combination thereof; the memory 411, the communication module 412, and the processor 414 can be connected through the bus 413; exemplarily, programs are stored in the memory 411; the processor 414 can call and run these programs from the memory 411, so that the rock stability assessment method based on three-dimensional geological modeling can be implemented by running the program.

[0119] In the embodiments provided in the present application, it should be understood that the disclosed methods and systems can be implemented in other ways. The system embodiments described above are only schematic. For example, the division of units is only a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some communication interfaces, devices or units, which can be electrical, mechanical or other forms. In addition, the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the present embodiment. In the present application, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply that there is any such actual relationship or order between these entities or operations. The above description is only an embodiment of the present application and is not intended to limit the scope of protection of the present application. For those skilled in the art, the present application may have various changes and variations. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present application should be included in the scope of protection of the present application. The embodiments of the present invention are described above in conjunction with the accompanying drawings, but the present invention is not limited to the above-mentioned specific embodiments. The above-mentioned specific embodiments are merely illustrative and not restrictive. Under the enlightenment of the present invention, those skilled in the art can also make many forms without departing from the scope of protection of the purpose of the present invention and the claims, which all belong to the protection of the present invention.

Claims

1. A rock formation stability assessment method based on three-dimensional geological modeling, characterized in that: The following steps are involved: S1: Acquire data of the area to be evaluated, and obtain a three-dimensional geological model according to the data of the area to be evaluated; S2: obtaining the parameters of the rock formation to be evaluated according to the three-dimensional geological model, and establishing a numerical model using finite element simulation software according to the parameters of the rock formation to be evaluated; S3: obtaining engineering requirements, and performing numerical simulation using the numerical model and the improved stability analysis formula according to the engineering requirements and the data of the area to be evaluated, to obtain a comprehensive analysis result; S4: acquiring actual observation data, comparing the comprehensive analysis result with the actual observation data, and adjusting the parameters of the three-dimensional geological model and the numerical model according to the comparison result to obtain an adjusted model; S5: According to the adjusted model, numerical simulation is performed using finite element simulation software to obtain the rock stability results of the area to be evaluated; Step S3 specifically includes: S31: Acquire engineering requirements, and obtain simulation time range, simulation space range, simulation characteristics and simulation parameters according to the engineering requirements and the data of the area to be evaluated; S32: performing numerical simulation using the numerical model according to the simulation time range, simulation space range, simulation characteristics and simulation parameters to obtain simulation process data, wherein the simulation process data includes stress, strain and displacement during the simulation process; S33: Calculate the comprehensive stability score of each time and space position according to the simulation process data using an improved stability analysis formula to obtain a comprehensive analysis result; Step S33 specifically includes: according to the simulation process data, using the improved stability analysis formula, calculating the comprehensive stability score of each time and space position, and obtaining a comprehensive analysis result, such as the formula: ， ， ， ， in, F ( t , x , y , z ) is the rock formation stability score, which changes with time t and spatial position ( x , y , z )change; w i ( t , x , y , z ) is the i The dynamic weight factor of each feature reflects the change of the feature over time and space; x i ' is the i The standardized feature values; α i and β i is a constant used to adjust the eigenvalue; It is j Dynamic weight factors, including temporal and spatial variables; d j ' is the j The normalized dynamic eigenvalues; θ k It is k The weighting factors of the uncertainty parameters; u k ' is the k The standardized uncertainty parameter value; , , and is a constant used to adjust the dynamic eigenvalues ​​and uncertainty parameters; m , n and p are the number of features, dynamic weight factors, and uncertainty parameters, respectively; w i ( t ) is the time factor, which is used to consider seasonal changes and long-term trends. a i and b i is a constant, ω is the frequency, φ i is the phase angle; w i ( x , y , z ) is a spatial factor, which is used to consider the differences in geographical locations and geological conditions. x 0, y 0, z 0) is the center position, σ i is the scale parameter, N is the number of simulations, f (.) indicates a simulation function. u k It is k uncertainty parameter values, is a random variable.

2. The rock formation stability assessment method based on three-dimensional geological modeling according to claim 1 is characterized in that: The data of the area to be evaluated include geological exploration data, remote sensing data, environmental data and historical data; the geological exploration data include drilling data, geological maps, seismic reflection data and core samples of the area to be evaluated; the remote sensing data include satellite images and LIDAR data of the area to be evaluated; the environmental data include groundwater level data and meteorological data of the area to be evaluated; the historical data include earthquake record data and landslide record data of the area to be evaluated; the rock formation parameters to be evaluated include physical parameters, mechanical parameters, geological structure parameters, hydrogeological parameters and environmental impact parameters; the physical parameters include density, porosity and water content of the rock formation to be evaluated; the mechanical parameters include uniaxial compressive strength, triaxial compressive strength, shear strength, elastic modulus and Poisson's ratio of the rock formation to be evaluated; the geological structure parameters include fault parameters and joint parameters; the hydrogeological parameters include permeability parameters and groundwater pressure parameters; the environmental impact parameters include temperature and rainfall.

3. The rock formation stability assessment method based on three-dimensional geological modeling according to claim 1, characterized in that: Step S1 specifically includes: S11: Acquire the data of the area to be evaluated, and preprocess the data of the area to be evaluated to obtain the preprocessed data of the area to be evaluated; S12: Based on the pre-processed data of the area to be evaluated, using modeling software, geological structure is performed to obtain a three-dimensional geological model.

4. The rock formation stability assessment method based on three-dimensional geological modeling according to claim 1, characterized in that: Step S2 specifically includes: S21: acquiring rock formation parameters to be evaluated according to the three-dimensional geological model; S22: importing the three-dimensional geological model into finite element simulation software, using structured grids, setting grid division parameters, running grid division tools, and generating grids for the numerical model; S23: Obtaining a numerical model according to the rock formation parameters to be evaluated and the grid of the numerical model.

5. The rock formation stability assessment method based on three-dimensional geological modeling according to claim 1, characterized in that: Step S4 specifically includes: S41: Acquire actual observation data, compare the comprehensive analysis result with the actual observation data, and obtain a model accuracy evaluation result; S42: adjusting the parameters and structure of the three-dimensional geological model according to the model accuracy evaluation result to obtain a new three-dimensional geological model; adjusting the physical and mechanical parameters of the numerical model according to the model accuracy evaluation result and the new three-dimensional geological model to obtain adjusted parameters and structural data; S43: Perform numerical simulation using finite element simulation software based on the adjusted parameters and structural data and the numerical model to obtain new simulation results, compare the new simulation results with the actual observation data again, and verify and adjust the three-dimensional geological model and the numerical model based on the results of the re-comparison to obtain an adjusted model.

6. The rock formation stability assessment method based on three-dimensional geological modeling according to claim 1, characterized in that: Step S5 specifically includes: S51: according to the adjusted model, setting the final simulation conditions, confirming the boundary conditions and initial conditions, and synchronously setting the gravity, ground stress, and groundwater pressure loading conditions; using finite element simulation software, performing numerical simulation to obtain simulation results; the simulation results include stress distribution diagrams, strain distribution diagrams, and displacement distribution diagrams; S52: According to the simulation results, the rock stability results of the area to be evaluated are obtained.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the rock formation stability assessment method based on three-dimensional geological modeling as described in any one of claims 1 to 6 are implemented.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the steps of the rock formation stability assessment method based on three-dimensional geological modeling as described in any one of claims 1-6 are implemented.

Images