Geotechnical engineering safety monitoring method based on muon detection technology

Abstract

The present invention relates to the technical field of geotechnical engineering safety monitoring, and provides a geotechnical engineering safety monitoring method based on muon detection technology. The method comprises: arranging a plurality of muon detection devices around or inside a geological body to be observed to form an array, and measuring in real time the muon observation flux of multi-angle cosmic rays passing through said geological body; converting the muon observation flux into the density length, and obtaining the three-dimensional density length distribution of said geological body; and analyzing the stress, strain, and deformation. In the present invention, muon particles are used to perform geological detection of geotechnical engineering, and the muons are particles having extremely high penetration capability, and can penetrate through the thick layer of the earth, so that detailed information about the interior of the geotechnical engineering is provided. In addition, the present invention has the advantages of being free of vibration and electromagnetic noise interference, and can implement accurate and rapid dynamic monitoring of the geotechnical engineering. The present invention has a wide application prospect in the fields of geological exploration, mineral exploration, underground resource development, and the like.

Description

A geotechnical engineering safety monitoring method based on muon detection technology

[0001] This application claims priority to the Chinese patent application filed with the China Patent Office on October 27, 2023, with application number 202311408536.2 and invention name “A method for geotechnical engineering safety monitoring based on muon detection technology”, the entire contents of which are incorporated by reference into this application. Technical Field

[0002] The present invention belongs to the technical field of geotechnical engineering safety monitoring, and in particular relates to a geotechnical engineering safety monitoring method based on muon detection technology. Background Art

[0003] In geotechnical engineering, various safety issues, such as landslides and collapses, can arise due to factors such as geological conditions, engineering design, and construction operations. These issues can severely impact the safety of personnel, equipment, and the surrounding environment. Therefore, safety monitoring of geotechnical engineering is crucial and one of the most pressing issues in current engineering construction.

[0004] Because geotechnical engineering projects often present hidden geotechnical problems and dynamic changes during construction and operation, such as hidden underground fault fracture zones, karst caves, and other unfavorable geological bodies, changes in physical properties caused by changes in moisture content during rainfall can significantly reduce the strength and stability of the geotechnical body, making it very likely to trigger geological disasters. To ensure the quality and safety of geotechnical engineering projects, it is necessary to monitor the response of the geotechnical body structure to changes in the underground geological environment. Scientific and reasonable monitoring methods can not only improve construction quality and efficiency, but also have a very important practical significance for ensuring project safety.

[0005] Currently, there are three main types of geotechnical engineering safety monitoring methods: direct geometric methods, indirect geometric methods, and physical methods. These methods primarily focus on directly monitoring geotechnical stress, strain, and displacement parameters. However, these methods have the following significant drawbacks: First, these methods rely on single-point monitoring and fail to comprehensively cover all possible monitoring points, leading to the neglect of some key monitoring points and an inability to effectively monitor overall geotechnical changes. Second, these methods require significant human and material resources, which increases project costs. Furthermore, human factors can lead to inaccurate monitoring results. For example, monitoring personnel may fail to promptly identify and address potential safety hazards, or may misreport monitoring data due to fatigue or negligence. Furthermore, due to the unique characteristics of geotechnical engineering, some monitoring methods fail to accurately reflect actual project conditions, thus affecting the accuracy of monitoring results. For example, in some unique geological conditions or complex construction environments, existing monitoring equipment and technologies cannot provide sufficient information for effective safety monitoring. Finally, these methods fail to promptly identify and address potential safety issues, increasing project risks. This is a significant drawback when real-time feedback and adjustments to construction plans are required.

[0006] In recent years, large-scale, complex numerical calculations have been widely used to predict deformation in geotechnical engineering, to some extent addressing the shortcomings of existing single-point monitoring technologies. However, this technology still has limitations. First, the accuracy of the finite element model directly affects the accuracy of the simulation results. If the model is not established accurately or is simplified, the prediction results will deviate from the actual situation, thereby affecting the accuracy of the decision. Second, finite element models typically contain many parameters, such as material properties, geometry, boundary conditions, etc. The selection and setting of these parameters can greatly affect the simulation results. Therefore, the selection and setting of parameters requires professional knowledge and experience, otherwise it will lead to inaccurate simulation results. Finally, the establishment of finite element models is often based on theories and assumptions. When the internal structure of the geotechnical mass is not clearly understood, it cannot fully reflect the various complex situations in the actual construction process, resulting in a certain deviation between the simulation results and the actual situation.

[0007] In summary, the existing geotechnical engineering monitoring technology has at least the following shortcomings: (1) The current monitoring technology has a small monitoring range and it is difficult to obtain the overall geological characteristics; (2) Some current monitoring equipment has weak noise and interference resistance capabilities, and is difficult to construct in the field; (3) When using finite element analysis, geophysical parameters are not fully combined, and the prediction accuracy is not high.

[0008] Summary of the Invention

[0009] In view of the defects of the existing technology, the present invention provides a geotechnical engineering safety monitoring method based on muon detection technology, which can effectively solve the various shortcomings of the current traditional monitoring technology methods.

[0010] The technical solution adopted in the present invention is as follows:

[0011] The present invention provides a geotechnical engineering safety monitoring method based on muon detection technology, comprising the following steps:

[0012] Step 1: Multiple muon detection devices are deployed around or inside the observed geological body to form an array, and the muon observation flux of cosmic rays passing through the observed geological body is detected in real time from multiple angles.

[0013] Step 2: Establish a conversion formula between the density length of the observed geological body and the muon observation flux, and convert the muon observation flux of the observed geological body obtained in step 1 into the density length;

[0014] Step 3, based on the density length of the observed geological body obtained in step 2, a geophysical method is used to perform inversion to obtain a three-dimensional density length distribution of the observed geological body;

[0015] Step 4: performing a drilling test on the observed geological body to obtain a true density length of a portion of the geological body, and using the true density length of the portion of the geological body to correct the three-dimensional density length distribution of the observed geological body obtained in step 3 to obtain a corrected three-dimensional density length distribution of the observed geological body;

[0016] Step 5, establishing a three-dimensional geological model of the observed geological body, importing the corrected three-dimensional density length distribution of the observed geological body into the three-dimensional geological model, and obtaining a three-dimensional geological model of the observed geological body having a three-dimensional density length distribution;

[0017] Step 6: Analyze the internal stress, strain, and deformation of the three-dimensional geological model of the observed geological body with a three-dimensional density length distribution using the finite element method. If the internal stress, strain, and deformation of the three-dimensional geological model of the observed geological body exceed a preset threshold, a safety warning is issued.

[0018] Preferably, in step 1, the muon observation flux is the muon decay flux, and its detection method is:

[0019] The muon flux detected by the muon detection device before the muon passes through the geological body in a certain direction is: The muon flux after passing through the geological body is Then the muon attenuation flux in this direction is defined as Where θ is the zenith angle, is the azimuth.

[0020] Preferably, in step 2, the conversion formula between the density length of the observed geological body and the muon observation flux is established as follows:

[0021] Through numerical simulation and experimental data calibration, the conversion formula between the density length of different rocks and the observed muon flux is obtained as follows:

[0022] Where: a, b, c are fitting parameters; a, b, c are calibrated according to the measured data; Λ is the density length.

[0023] Preferably, in step 3, the geophysical method includes an alternating direction multiplier method and a conjugate gradient algorithm.

[0024] Preferably, in step 5, establishing a three-dimensional geological model of the observed geological body includes:

[0025] Step 5.1: Obtain the surface elevation image of the observed geological body:

[0026] Download the world surface elevation image map, import the geospatial position data of the left and right sections of the observed geological body into the world surface elevation image map, select the observed geological body area, download the selected observed geological body area, and obtain the surface elevation image map of the observed geological body;

[0027] Step 5.2, conversion of surface elevation image:

[0028] Importing the obtained surface elevation image of the observed geological body into map drawing software for projection conversion to obtain a surface elevation data file of the observed geological body;

[0029] Step 5.3, geophysical information processing of the observed geological body section:

[0030] Obtaining a two-dimensional geological profile along the left and right sections of the observed geological body; wherein the obtained two-dimensional geological profile includes both section information and geophysical information; importing the two-dimensional geological profile along the left and right sections of the observed geological body into a CAD program; using the CAD program to separate the section information and the geophysical information, and respectively deriving section profiles of the left and right sections of the observed geological body and geophysical density maps of the left and right sections of the observed geological body;

[0031] Step 5.4, CATIA 3D geological modeling:

[0032] Import the surface elevation data file of the observed geological body obtained in step 5.2 into the reverse engineering surface module of CATIA to generate a point cloud, and then transfer the point cloud into the quick surface module to generate the geological body surface surface, that is, the left and right cross-section entities of the geological body;

[0033] Map the cross-sectional views and geophysical density maps corresponding to the left and right sections of the geological body obtained in step 5.3 to the left and right section entities of the geological body to obtain the left and right section model views of the geological body;

[0034] Different topographic maps, including lithology maps, density maps, and fault maps, are extracted from the left and right section model maps of the geological body. The extracted left and right section entities of the geological body are spliced ​​with the topographic maps to form a three-dimensional geological model, completing the three-dimensional modeling of the observed geological body.

[0035] Preferably, step 6 specifically includes:

[0036] Step 6.1, import the three-dimensional geological model of the observed geological body into the three-dimensional finite element analysis software, and import the three-dimensional density length distribution of the observed geological body into the material library;

[0037] Step 6.2: Use 3D finite element analysis software to mesh the 3D geological model, add boundary conditions, and solve to view the stress, strain, and deformation inside the geological body.

[0038] Step 6.3: Based on the three-dimensional density length distribution of the observed geological body obtained by muon monitoring inversion, finite element analysis is continuously performed to obtain the latest internal stress, strain and deformation of the observed geological body;

[0039] Step 6.4, by setting the safety warning threshold, when the stress, strain and deformation data inside the observed geological body exceed the safety warning threshold, it is determined that there is a safety risk and an early warning message is issued.

[0040] The geotechnical engineering safety monitoring method based on muon detection technology provided by the present invention has the following advantages:

[0041] This technology uses muon particles for geotechnical geological surveys. Muons are highly penetrating particles that can penetrate thick layers of the Earth, providing detailed information about the interior of geotechnical structures. Muons are also immune to vibration and electromagnetic noise, enabling accurate and rapid dynamic monitoring of geotechnical structures. This technology has broad application prospects in geological exploration, mineral prospecting, and underground resource development.

[0042] Figures in the specification

[0043] FIG1 is a schematic flow chart of a geotechnical engineering safety monitoring method based on muon detection technology provided by the present invention;

[0044] FIG2 is a schematic diagram of a multi-angle array muon detection arrangement provided by the present invention;

[0045] FIG3 is a diagram showing simulation and experimental data of the conversion relationship between the standard rock density length and the muon observation flux provided by the present invention;

[0046] FIG4 is a flow chart of three-dimensional geological modeling provided by the present invention;

[0047] FIG5 is a flow chart of monitoring and early warning provided by the present invention. DETAILED DESCRIPTION

[0048] In order to make the technical problems, technical solutions and beneficial effects solved by the present invention more clearly understood, the present invention is further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0049] The present invention provides a geotechnical engineering safety monitoring method based on muon detection technology, with reference to FIG1 . Specifically, the method comprises the following steps:

[0050] Step 1: Use multiple muon detection devices to surround or inside the observed geological body to form an array, and detect in real time the muon observation flux of multi-angle cosmic rays passing through the observed geological body.

[0051] Specifically, as shown in Figure 2, a schematic diagram of the multi-angle array muon detection arrangement is provided. The specific calculation formula for the muon flux observed by a single muon detector is as follows:

[0052] Where: θ is the zenith angle, is the azimuth.

[0053] For the muon detection device The number of muons detected in the direction.

[0054] S is the detection area, Ω is the solid angle, and t is the measurement time. The solid angle (Ω, unit: sr) is the projected area of ​​any object on the unit sphere with the observation point as the center of the unit sphere. Here, it refers to the projected area of ​​the pixel element in a receiving unit onto the unit sphere with the center of the receiving unit as the center, which is used to constrain the spatial angular range of the muon ray.

[0055] is the zenith angle θ and azimuth angle of the muon detector The detected muon flux is the number of muons that arrive at a certain position on the sea level from a certain angle (including zenith angle and azimuth angle) with a certain energy within a unit area dS, time dt, and solid angle dΩ.

[0056] Since the natural muon flux is affected by environmental factors such as temperature, air pressure, and solar activity, the observed muon flux is unstable. When the muon flux is converted into density length, the density length result may produce a large error. Therefore, in the present invention, the muon flux is actually the muon attenuation flux. The muon flux mentioned in the following description refers to the muon attenuation flux. Specifically, the present invention directly measures the attenuation flux of muons before and after passing through the geological body through the muon detection device to avoid the influence of environmental fluctuations on the data results. Assume that the flux of muons before passing through the geological body in a certain direction is The flux after passing through the geological body is Then the muon attenuation flux in this direction is defined as

[0057] Step 2: Establish a conversion formula between the density length of the observed geological body and the muon observation flux, and convert the muon observation flux of the observed geological body obtained in step 1 into the density length;

[0058] Specifically, muons will experience a certain amount of energy loss when passing through geological bodies. The energy loss equation is:

[0059] Where: E is the muon energy, x is its position along the path, ρ(x) is the density of the rock, α and β are the rock ionization and radiation energy losses, respectively.

[0060] By changing the formula, we can get:

[0061] Where: Λ is the density length.

[0062] As shown in Figure 3, the conversion relationship between the density length of a standard rock and the observed muon flux is simulated and experimental data. The present invention obtains the empirical formula for the conversion between the density length of different rocks and the observed muon flux through numerical simulation and experimental data:

[0063] Among them, a, b, and c are fitting parameters; a, b, and c are calibrated according to the measured data.

[0064] Therefore, the present invention considers factors such as altitude, longitude and latitude, solar activity, and lithology, and calibrates the conversion formula between the density length of the observed geological body and the muon observation flux through actual on-site observation data to obtain the final conversion formula.

[0065] The advantage of converting the observed muon flux of the observed geological body into density length through the empirical formula is that density length is a very important concept in the field of muon detection. Density length, also known as opacity, is the integral of density along the ray path. Linear inversion based on opacity data can quickly reflect changes in the relative density of the detected target.

[0066] Step 3: Based on the density length of the observed geological body obtained in step 2, a geophysical method is used to perform inversion to obtain a three-dimensional density length distribution of the observed geological body.

[0067] Specifically, after converting the measured muon flux into density length, geophysical methods (including but not limited to the alternating direction multiplier method (ADMM), conjugate gradient and other algorithms) are used for inversion to obtain the three-dimensional density length distribution of the observed geological body.

[0068] Furthermore, from the perspective of the algorithm, for the model parameter m (the vector composed of the density values ​​of each voxel in the model), the density length d = F(m) corresponding to the model is forward calculated, hoping that d can be as close as possible to the observed density length d obs , that is, we hope that the bi-norm of the difference between the two can be as small as possible. The objective function is as follows:

[0069] Solving this problem can employ, but is not limited to, the extended alternating direction method of multipliers (eADMM). eADMM is a computational framework for solving distributed convex optimization problems. Through a decomposition and coordination process, it decomposes a large global problem into multiple smaller, more easily solved local subproblems. The solutions to these subproblems are then coordinated to yield a solution to the larger global problem. Using the alternating direction method of multipliers, the three-dimensional density-length distribution of the observed geological volume can be inverted.

[0070] Step 4: Conduct a drilling test on the observed geological body to obtain the true density length of part of the geological body, and use the true density length of part of the geological body to correct the three-dimensional density length distribution of the observed geological body obtained in step 3 to obtain a corrected three-dimensional density length distribution of the observed geological body.

[0071] Furthermore, the correction method is specifically as follows:

[0072] Step 4.1: Conduct a partial drilling test on the observed geological body to obtain the true density length of the partial geological body.

[0073] In step 4.2, the three-dimensional density length distribution of the observed geological body obtained by inversion in step 3 is subjected to multiple comparison correction with the true density length of some geological bodies obtained by drilling tests. The correction method can adopt the Bonferroni method, the False Discovery Rate (FDR) method, the Random-field Theory (RFT) method, etc.

[0074] Step 5: Establish a three-dimensional geological model of the observed geological body, import the corrected three-dimensional density length distribution of the observed geological body into the three-dimensional geological model, and obtain a three-dimensional geological model of the observed geological body with a three-dimensional density length distribution.

[0075] Furthermore, as shown in FIG4 , which is a flow chart of 3D geological modeling, the 3D geological modeling method is as follows:

[0076] Step 5.1, download the world surface elevation image map, import the kml file of the geospatial location data of the left and right sections of the observed geological body into the world surface elevation image map, select the observed geological body area, download the selected observed geological body area, and obtain the surface elevation image map of the observed geological body.

[0077] Step 5.2, conversion of surface elevation image:

[0078] The downloaded surface elevation image of the observed geological body is imported into the map drawing software Globalmapper for projection conversion, and the required surface elevation data file of the observed geological body in xyz format is exported.

[0079] Step 5.3, geophysical information processing of the observed geological body section:

[0080] Obtain a two-dimensional geological profile along the left and right sections of the observed geological body; wherein the obtained two-dimensional geological profile includes both section information and geophysical information; import the two-dimensional geological profile along the left and right sections of the observed geological body into CAD; use CAD to separate the section information and the geophysical information, and respectively derive the section profiles of the left and right sections of the observed geological body and the geophysical density maps of the left and right sections of the observed geological body.

[0081] Save the cross-sectional views of the left and right sections of the observed geological body and the geophysical density maps of the left and right sections of the observed geological body as dxf format files; among them, there are two cross-sectional views, corresponding to the left section and the right section of the geological body respectively; there are two geophysical density maps, corresponding to the left section and the right section of the geological body respectively.

[0082] Step 5.4, CATIA 3D geological modeling:

[0083] Import the surface elevation data file of the observed geological body obtained in step 5.2 into the reverse engineering surface module of CATIA to generate a point cloud, and then transfer the point cloud to the quick surface module to generate the surface surface of the geological body, that is, the left and right cross-section entities of the geological body.

[0084] Map the cross-sectional views and geophysical density maps corresponding to the left and right sections of the geological body obtained in step 5.3 to the left and right section entities of the geological body to obtain the left and right section model views of the geological body.

[0085] Different topographic maps, including lithology maps, density maps, and fault maps, are extracted from the left and right section model maps of the geological body. The extracted left and right section entities of the geological body are spliced ​​with the topographic maps to form a three-dimensional geological model, completing the three-dimensional modeling of the observed geological body.

[0086] Step 6: Analyze the internal stress, strain, and deformation of the three-dimensional geological model of the observed geological body with a three-dimensional density length distribution using the finite element method. If the internal stress, strain, and deformation of the three-dimensional geological model of the observed geological body exceed a preset threshold, a safety warning is issued.

[0087] Furthermore, as shown in FIG5 , which is a monitoring and early warning flow chart, in step 6, the early warning method is:

[0088] Step 6.1: Import the three-dimensional geological model of the geological body obtained in step 5.4 into the three-dimensional finite element analysis software in *.dxf format, and import the three-dimensional density length distribution of the observed geological body into the material library.

[0089] Step 6.2: Use 3D finite element analysis software to mesh the 3D geological model, add boundary conditions, and solve to view the stress, strain, and deformation inside the geological body.

[0090] In step 6.3, based on the three-dimensional density length distribution of the observed geological body obtained by muon monitoring inversion, finite element analysis is continuously performed to obtain the latest internal stress, strain and deformation of the observed geological body.

[0091] Step 6.4, by setting the safety warning threshold, when the stress, strain and deformation data inside the observed geological body exceed the safety warning threshold, it is determined that there is a safety risk and an early warning message is issued.

[0092] The geotechnical engineering safety monitoring method based on muon detection technology provided by the present invention has the following advantages:

[0093] This technology uses muon particles for geotechnical geological surveys. Muons are highly penetrating particles that can penetrate thick layers of the Earth, providing detailed information about the interior of geotechnical structures. Muons are also immune to vibration and electromagnetic noise, enabling accurate and rapid dynamic monitoring of geotechnical structures. This technology has broad application prospects in geological exploration, mineral prospecting, and underground resource development.

[0094] The above is only a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principles of the present invention. These improvements and modifications should also be considered within the scope of protection of the present invention.

Claims

1. A geotechnical engineering safety monitoring method based on muon detection technology, characterized in that: The following steps are involved: Step 1: multiple muon detection devices are arranged around or inside the observed geological body to form an array, and the muon observation flux of multi-angle cosmic rays passing through the observed geological body is detected in real time; Step 2, establishing a conversion formula between the density length of the observed geological body and the muon observation flux, and converting the muon observation flux of the observed geological body obtained in step 1 into the density length; Step 3, based on the density length of the observed geological body obtained in step 2, a geophysical method is used to perform inversion to obtain a three-dimensional density length distribution of the observed geological body; Step 4, performing a drilling test on the observed geological body to obtain a true density length of a portion of the geological body, and using the true density length of the portion of the geological body to correct the three-dimensional density length distribution of the observed geological body obtained in step 3 to obtain a corrected three-dimensional density length distribution of the observed geological body; Step 5, establishing a three-dimensional geological model of the observed geological body, importing the corrected three-dimensional density length distribution of the observed geological body into the three-dimensional geological model, and obtaining a three-dimensional geological model of the observed geological body with a three-dimensional density length distribution; Step 6: Analyze the internal stress, strain and deformation of the three-dimensional geological model of the observed geological body with three-dimensional density length distribution by finite element method. If the internal stress, strain and deformation of the three-dimensional geological model of the observed geological body exceed the preset threshold, issue a safety warning.

2. A geotechnical engineering safety monitoring method based on muon detection technology according to claim 1, characterized in that: In step 1, the muon observation flux is the muon decay flux, and its detection method is: The muon flux detected by the muon detection device before the muon passes through the geological body in a certain direction is The muon flux after passing through the geological body is Then the muon attenuation flux in this direction is defined as Where θ is the zenith angle, is the azimuth.

3. A geotechnical engineering safety monitoring method based on muon detection technology according to claim 2, characterized in that: In step 2, the conversion formula between the density length of the observed geological body and the muon observation flux is established as follows: Through numerical simulation and experimental data calibration, the conversion formula between the density length of different rocks and the observed muon flux is obtained as follows: Where: a, b, c are fitting parameters; a, b, c are calibrated according to measured data; Λ is the density length.

4. The method for geotechnical engineering safety monitoring based on muon detection technology according to claim 1 is characterized in that: In step 3, the geophysical method includes an alternating direction multiplier method and a conjugate gradient algorithm.

5. The method for geotechnical engineering safety monitoring based on muon detection technology according to claim 1 is characterized in that: In step 5, establishing a three-dimensional geological model of the observed geological body includes: Step 5.1, obtain the surface elevation image of the observed geological body: Download the world surface elevation image map, import the geospatial position data of the left and right sections of the observed geological body into the world surface elevation image map, select the observed geological body area, download the selected observed geological body area, and obtain the surface elevation image map of the observed geological body; Step 5.2, conversion of surface elevation image map: Importing the obtained surface elevation image of the observed geological body into the map drawing software for projection conversion to obtain the surface elevation data file of the observed geological body; Step 5.3, geophysical information processing of the observed geological body section: Obtain a two-dimensional geological profile along the left and right sections of the observed geological body; wherein the obtained two-dimensional geological profile includes both section information and geophysical information; import the two-dimensional geological profile along the left and right sections of the observed geological body into CAD; separate the section information and the geophysical information by using CAD, and respectively derive the section diagram of the left and right sections of the observed geological body and the geophysical density map of the left and right sections of the observed geological body; Step 5.4, CATIA 3D geological modeling: Import the surface elevation data file of the observed geological body obtained in step 5.2 into the reverse engineering surface module of CATIA to generate a point cloud, and then transfer the point cloud to the rapid surface module to generate the surface surface of the geological body, which is the left and right cross-section entities of the geological body; Map the cross-sectional diagram and geophysical density map corresponding to the left and right sections of the geological body obtained in step 5.3 to the left and right section entities of the geological body to obtain the left and right section model diagrams of the geological body; Different topographic maps, including lithology map, density map and fault map, are extracted from the left and right section model maps of the geological body. The extracted left and right section entities of the geological body are spliced ​​with the topographic map to form a three-dimensional geological model, thus completing the three-dimensional modeling of the observed geological body.

6. The method for geotechnical engineering safety monitoring based on muon detection technology according to claim 1 is characterized in that: Step 6 specifically includes: Step 6.1, import the three-dimensional geological model of the observed geological body into the three-dimensional finite element analysis software, and import the three-dimensional density length distribution of the observed geological body into the material library; Step 6.2, the 3D finite element analysis software is used to mesh the 3D geological model, add boundary conditions, and solve to view the internal stress, strain, and deformation of the geological body; Step 6.3, based on the three-dimensional density length distribution of the observed geological body obtained by muon monitoring inversion, finite element analysis is continuously performed to obtain the latest internal stress, strain and deformation of the observed geological body; Step 6.4, by setting the safety warning threshold, when the internal stress, strain and deformation data of the observed geological body exceed the safety warning threshold, it is determined that there is a safety risk and a warning message is issued.

Images