Grouting fissure space positioning method, device and equipment and storage medium

Abstract

This invention provides a method, apparatus, device, and storage medium for spatial location of grouting fractures, relating to the field of geotechnical engineering technology. The method includes: acquiring three-component magnetic anomaly data generated by a fracture diffuser at at least three different observation heights; performing phase transformation processing on each component of the magnetic anomaly data to obtain equivalent vertical magnetization anomaly data corresponding to each observation height; extracting zero-intersection point location information at each observation height based on each equivalent vertical magnetization anomaly data; and calculating the spatial geometric parameters of the fracture diffuser based on the change in each zero-intersection point location information, the spatial geometric parameters including at least the burial depth. This invention enables non-destructive and quantitative location of grouting fractures, which can be used for evaluating grouting effects and optimizing grouting design in engineering projects.

Description

Technical Field

[0001] This invention relates to the field of geotechnical engineering technology, and in particular to a method, apparatus, equipment and storage medium for spatial positioning of grouting fractures. Background Technology

[0002] In rock engineering, grouting is a key means to improve the mechanical properties and permeability of rock mass. The accurate evaluation of grouting effect depends on the effective detection of the spatial location and geometric morphology of fracture diffusers.

[0003] Existing non-destructive testing methods, such as seismic wave methods, ground penetrating radar, and cross-hole tomography, generally suffer from insufficient positioning accuracy, limited detection range, and complex or costly on-site implementation in grouting location, making it difficult to meet the needs of engineering sites for rapid, accurate, and low-cost evaluation of grouting effects. Magnetic methods, as a mature geophysical exploration technique, have advantages such as simple equipment, low cost, and intuitive interpretation. However, traditional magnetic methods are mainly used for detecting large geological bodies or metallic targets, and there is a lack of dedicated magnetic detection and inversion methods for thin-layered grouting fractures, especially under conditions of oblique geomagnetic fields, where the interpretation of magnetic anomalies is quite difficult.

[0004] Therefore, there is an urgent need to design a spatial positioning scheme for grouting fractures that can achieve non-destructive and quantitative inversion of the spatial location and geometry of thin-layer grouting fractures. Summary of the Invention

[0005] This invention provides a method, apparatus, equipment, and storage medium for spatial positioning of grouting fractures, which solves the problems of difficulty in quantitatively locating the spatial morphology of grouting fractures and difficulty in interpreting magnetic anomalies under oblique magnetic field conditions in the prior art.

[0006] This invention provides a method for spatially locating grouting fractures. The method includes: acquiring three-component magnetic anomaly data generated by a fracture diffuser at at least three different observation heights; performing phase transformation processing on each three-component magnetic anomaly data to obtain equivalent vertical magnetization anomaly data corresponding to each observation height; extracting zero-intersection point location information at each observation height based on each equivalent vertical magnetization anomaly data; and calculating the spatial geometric parameters of the fracture diffuser based on the change in each zero-intersection point location information, wherein the spatial geometric parameters include at least the burial depth.

[0007] According to the present invention, a method for spatial positioning of grouting fissures is provided, wherein the fissure diffuser is formed by solidifying a magnetic grouting material, and the magnetic grouting material is a mixture of cement-based grouting material and fine-grained magnetite powder.

[0008] According to the present invention, a spatial positioning method for grouting fractures includes at least three different observation heights: a reference observation surface located on the ground surface, a first observation surface located at a first height, and a second observation surface located at a second height; wherein the second height is greater than the first height, and the projections of the reference observation surface, the first observation surface, and the second observation surface on the horizontal plane coincide with each other.

[0009] According to the present invention, a spatial positioning method for grouting fractures is provided, which extracts the zero-intersection position information at each observation height based on each equivalent vertical magnetization magnetic anomaly data, including: performing symmetry axis analysis on each equivalent vertical magnetization magnetic anomaly data to determine the symmetry axis of the total magnetic field anomaly, wherein the direction of the symmetry axis of the total magnetic field anomaly is a preset survey line direction; selecting a one-dimensional magnetic survey profile along the symmetry axis of the total magnetic field anomaly; and extracting the zero-intersection position information at each observation height on the one-dimensional magnetic survey profile.

[0010] According to the spatial positioning method for grouting fractures provided by the present invention, the zero-intersection point position information at each observation height is extracted, including: determining the zero-point position of the horizontal component of the equivalent vertical magnetization magnetic anomaly data changing from positive to negative or from negative to positive along the preset survey line direction at each observation height, as the horizontal zero-intersection point; and determining the two symmetrical zero-point positions of the vertical component of the equivalent vertical magnetization magnetic anomaly data along the preset survey line direction at each observation height, as a pair of vertical zero-intersection points.

[0011] According to the present invention, a spatial positioning method for grouting fractures is provided, which calculates the spatial geometric parameters of the fracture diffuser body by inverting the change in the position information of each zero intersection point. This includes: calculating the burial depth, dip angle and horizontal projection scale of the fracture diffuser body by analytical relationship based on the rate of change of the position of each horizontal zero intersection point and the rate of change of the distance between each vertical zero intersection point with the observation height.

[0012] According to the spatial location method for grouting fractures provided by the present invention, phase transformation processing is performed on each three-component magnetic anomaly data to obtain equivalent vertical magnetization magnetic anomaly data corresponding to each observation height. The method includes: converting each three-component magnetic anomaly data into magnetic anomaly data under vertical magnetization conditions through frequency domain transformation to obtain equivalent vertical magnetization magnetic anomaly data; wherein, the frequency domain transformation adopts Fourier transform, and the same geomagnetic field parameters are used to construct the phase transformation operator for each observation height.

[0013] The present invention also provides a grouting fracture spatial positioning device, which includes: a data acquisition module for acquiring three-component magnetic anomaly data generated by a fracture diffuser at at least three different observation heights, wherein the fracture diffuser is formed by grouting in the rock mass; a phase transformation module for performing phase transformation processing on each three-component magnetic anomaly data to obtain equivalent vertical magnetization anomaly data corresponding to each observation height; a zero-intersection extraction module for extracting zero-intersection position information at each observation height based on each equivalent vertical magnetization anomaly data; and a joint inversion module for inverting and calculating the spatial geometric parameters of the fracture diffuser based on the change in the zero-intersection position information, wherein the spatial geometric parameters include at least the burial depth.

[0014] 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 computer program, it implements the grouting fracture spatial positioning method described above.

[0015] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the grouting fracture spatial positioning method described in any of the above claims.

[0016] The present invention also provides a computer program product, including a computer program, which, when executed by a processor, implements the grouting fracture spatial positioning method described in any of the above claims.

[0017] The spatial location method, apparatus, equipment, and storage medium for grouting fractures provided by this invention acquire three-component magnetic anomaly data at at least three different observation heights, expanding the magnetic anomaly information from a single observation height into a multi-height observation system containing depth-sensitive information; through phase transformation processing, complex magnetic anomalies under oblique magnetization conditions are converted into equivalent vertical magnetization magnetic anomalies, eliminating the interference of geomagnetic field direction on inversion accuracy; by extracting the zero-intersection point location information at each observation height, the spatial distribution characteristics of the magnetic anomalies are transformed into quantifiable geometric feature points; by inverting and calculating spatial geometric parameters based on the changes in the zero-intersection point location information at each observation height, the difference information between multiple heights is used to achieve a joint quantitative solution for burial depth, dip angle, and horizontal projection scale, thereby realizing non-destructive and quantitative location of grouting fractures. Attached Figure Description

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

[0019] Figure 1 This is a schematic flowchart of a grouting fracture spatial positioning method provided in an embodiment of the present invention.

[0020] Figure 2 This is a schematic diagram of a multi-height magnetic anomaly observation setup provided in an embodiment of the present invention.

[0021] Figure 3 This is a schematic diagram of the structure of a grouting crack spatial positioning device provided in an embodiment of the present invention.

[0022] Figure 4 This is a schematic diagram of the physical structure of an electronic device provided in an embodiment of the present invention. Detailed Implementation

[0023] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention.

[0024] The following is combined with Figure 1 This invention describes a method for spatially locating grouting fractures. For consistency, the entity performing this method will be uniformly referred to as the "system," and will not be described further thereafter.

[0025] Figure 1 This is a flowchart illustrating the grouting fracture spatial positioning method provided in an embodiment of the present invention, as shown below. Figure 1 As shown, the method includes the following: S101. Acquire three-component magnetic anomaly data generated by the fracture diffuser from at least three different observation heights.

[0026] It should be noted that this invention acquires vector magnetic anomaly data containing sensitive information about the depth of grouting fractures through multi-height magnetic anomaly data acquisition, providing a basic data source for subsequent inversion calculations and breaking through the limitation of single observation height on the positioning accuracy of deeply buried grouting bodies.

[0027] In this embodiment of the invention, the fissure diffuser is formed by solidifying a magnetic grouting material, which is a mixture of cement-based grouting material and fine-grained magnetite powder.

[0028] It should be noted that cement-based grouting materials are the conventional grouting base materials used in rock grouting projects. They do not possess identifiable magnetic response characteristics and are classified as magnetically inert materials. Fine-grained magnetite powder, on the other hand, is a powder of iron oxide mineral, which is the material with the highest magnetic susceptibility and strongest remanence among natural ferromagnetic minerals. It has the characteristics of strong acid resistance, fine particle size, and stable performance. As a harmless inert ferromagnetic material, it will not chemically react with cement hydration products and can effectively fill the gaps between cement hydration agglomerates. Its mixing ratio can be flexibly adjusted according to the magnetization intensity required by the target grouting project.

[0029] In some embodiments, after the magnetic grouting material is injected into the fissures of the target rock mass, it diffuses within the fissures under grouting pressure. After the grouting material solidifies, a thin layer of magnetic material with stable magnetization properties is formed within the fissures, i.e., a fissure diffuser. This fissure diffuser can generate a recognizable magnetic anomaly signal in the Earth's magnetic field, providing a signal source for magnetic detection.

[0030] For example, in practical engineering applications, the proportion of magnetite powder can be determined based on the development of fissures in the target rock mass and the diffusion range of the grouting design. The magnetite powder is then uniformly mixed with cement-based grouting material and mixing water to prepare magnetic grouting material. The magnetic grouting material is then injected into the target fissures of the target rock mass under high pressure through grouting boreholes. After the grouting material has completed hydration and solidification, a stable fissure diffuser is formed, and subsequent magnetic anomaly data acquisition can be carried out.

[0031] Specifically, in the preparation process of magnetic grouting materials, the water-cement ratio of the mixing water can be adjusted according to the conventional design requirements of rock mass grouting engineering. The change in the water-cement ratio will not have a significant impact on the magnetization characteristics of the magnetic grouting material after curing. The magnetic susceptibility of the grout body after curing increases linearly with the increase of the amount of magnetite powder added. By adjusting the amount of magnetite powder added, the cured fracture diffuser can generate a magnetic anomaly signal that meets the identification accuracy requirements of magnetic measuring equipment.

[0032] In this embodiment of the invention, at least three different observation heights include: a reference observation surface located on the ground surface, a first observation surface located at a first height, and a second observation surface located at a second height.

[0033] The second height is greater than the first height, and the projections of the reference observation surface, the first observation surface, and the second observation surface on the horizontal plane coincide with each other.

[0034] It should be noted that the reference observation surface is the observation surface that is closest to the fracture diffuser in vertical distance, usually the surface plane of the target grouting area, serving as the reference surface for multi-height observations; the first height and the second height are the vertical distances between the observation surface and the reference observation surface. The selection of the observation height can be determined according to the size of the target grouting area and the expected burial depth of the grouting fracture, ensuring that there are significant differences in depth-sensitive information between the observation data at different heights, and avoiding the inability of the magnetic anomaly change to be effectively identified by the magnetic measurement equipment due to the small interval between observation heights.

[0035] For example, the first height can be a preset ratio of the expected maximum burial depth of the grouting fracture, and the second height can be twice the first height. By setting the height intervals in equal arithmetic or equal proportions, the magnetic anomaly data at different observation heights can have a quantifiable gradient of change, which is convenient for extracting the change in the zero intersection position.

[0036] Specifically, the arrangement of multi-altitude magnetic anomaly observations can be as follows: Figure 2 As shown.

[0037] Figure 2 This is a schematic diagram of a multi-height magnetic anomaly observation layout provided by an embodiment of the present invention. Above the target grouting area 210 after grouting is completed, along a direction perpendicular to the ground surface, a reference observation surface 220, a first observation surface 221, and a second observation surface 222 are arranged sequentially from bottom to top on the ground surface. All three observation surfaces are planes parallel to the ground surface, and the projections of the three observation surfaces on the horizontal plane completely overlap, ensuring that the magnetic measurement points at different observation heights have a one-to-one correspondence in horizontal position. Magnetic sensors for collecting magnetic anomaly data are arranged on each observation surface.

[0038] Among them, the low-position magnetic sensor 230 is deployed on the first observation surface, and the high-position magnetic sensor 231 is deployed on the second observation surface. Magnetic anomaly data at corresponding positions are collected by magnetic sensors at different heights.

[0039] For example, the distance between the first observation surface 221 and the second observation surface 222 is H1; the distance between the first observation surface 221 and the reference observation surface 220 is H2, and H2 is less than H1; the distance between the second observation surface 221 and the reference observation surface 220 is H1+H2.

[0040] Thus, by setting at least three observation heights with overlapping projections, this invention enables simultaneous observation of magnetic anomalies at multiple heights in the same grouting fracture area. It expands the two-dimensional magnetic anomaly information at a single observation height into three-dimensional magnetic response features that include the depth dimension, introducing core depth-sensitive information for subsequent inversion calculations and effectively improving the reliability of locating deeply buried grouting fracture bodies.

[0041] In some embodiments, the acquisition of three-component magnetic anomaly data can be accomplished by a three-component magnetometer, which is a magnetic measurement device that can simultaneously acquire three spatial components of the magnetic anomaly signal. The acquired three-component magnetic anomaly data are vector data of the horizontal x-component, horizontal y-component, and vertical z-component of the magnetic anomaly signal under the geomagnetic field environment.

[0042] For example, in the actual data acquisition process, magnetic measurement points are first set up in a regular grid pattern on each observation surface of the target grouting area. The grid spacing can be adjusted according to the expected diffusion scale of the grouting cracks. Then, the three-component magnetometer is used to collect the original three-component magnetic anomaly data of each measurement point at the corresponding observation height in a preset measurement point order. During the acquisition process, the attitude of the three-component magnetometer is kept fixed to avoid measurement errors caused by changes in the sensor attitude.

[0043] Specifically, before collecting the three-component magnetic anomaly data, it is necessary to first collect and smooth the background magnetic field of the target grouting area to eliminate the influence of environmental electromagnetic interference on the magnetic anomaly data. The background magnetic field collection should be carried out before grouting construction, at the same measuring point location and the same observation height in the target grouting area, to collect the geomagnetic field background data without the influence of the grouting fracture diffuser. After the magnetic anomaly data collection is completed after grouting, the corresponding background magnetic field data is subtracted from the measured original magnetic anomaly data to obtain the three-component magnetic anomaly data generated by the fracture diffuser alone.

[0044] Meanwhile, high-frequency electronic noise during the acquisition process can be eliminated through smoothing by moving average, thereby improving the signal-to-noise ratio of magnetic anomaly data.

[0045] S102. Perform phase transformation processing on each three-component magnetic anomaly data to obtain the equivalent vertical magnetization magnetic anomaly data corresponding to each observation height.

[0046] It should be noted that this invention addresses the objective situation in actual engineering where the geomagnetic field is an oblique magnetic field. Through phase transformation processing, it converts complex magnetic anomaly data under oblique magnetization conditions into equivalent vertical magnetization magnetic anomaly data, eliminates the influence of the geomagnetic field tilt angle on the distribution pattern of magnetic anomalies, satisfies the assumptions of the magnetization direction in the subsequent inversion model, and reduces the difficulty of interpreting magnetic anomalies under oblique magnetic field conditions.

[0047] In some embodiments, each three-component magnetic anomaly data can be converted into magnetic anomaly data under vertical magnetization conditions through frequency domain transformation to obtain equivalent vertical magnetization magnetic anomaly data.

[0048] It should be noted that phase transformation processing is the pole reduction processing in the field of magnetic exploration, that is, the reduction to magnetic pole processing, hereinafter referred to as pole reduction processing. This processing can convert magnetic anomalies under any magnetization direction into magnetic anomalies under vertical magnetization and vertical observation conditions, so that the distribution pattern of magnetic anomalies corresponds to the horizontal projection range of the magnetic body, and eliminates the magnetic anomaly offset and distortion problems caused by oblique magnetization.

[0049] For example, the phase transformation processing is performed in the Fourier frequency domain. First, the three-component magnetic anomaly data in the spatial domain is converted to the frequency domain by Fourier transform. Then, the phase of the frequency domain data is adjusted by the phase transformation operator. Finally, the processed frequency domain data is converted back to the spatial domain by inverse Fourier transform to obtain the equivalent perpendicular magnetization magnetic anomaly data.

[0050] The frequency domain transformation uses Fourier transform, and the same geomagnetic field parameters are used to construct the phase transformation operator at each observation height.

[0051] In this embodiment of the invention, the geomagnetic field parameters include at least the geomagnetic field dip angle and the geomagnetic field declination corresponding to the target grouting area. These parameters can be obtained through a geomagnetic reference field model at the location of the target area or through on-site geomagnetic measurements.

[0052] In addition, for all observation heights of the same target grouting area, the same set of geomagnetic field parameters is used to construct the phase conversion operator to ensure that the equivalent vertical magnetization magnetic anomaly data after conversion at different observation heights have parameter consistency and comparability, and to avoid systematic errors introduced due to differences in geomagnetic field parameters.

[0053] For example, in the actual phase transformation process, the geomagnetic field dip angle and magnetic declination of the target grouting area are first obtained, and a phase transformation operator suitable for three-component magnetic anomaly data is constructed based on these parameters. Then, for the three-component magnetic anomaly data at each observation height, independent frequency domain phase transformation processing is performed. That is, the same phase transformation operator is used to perform pole-shaving processing on the three-component magnetic anomaly data of the reference observation surface, the first observation surface, and the second observation surface to obtain the equivalent vertical magnetization three-component magnetic anomaly data corresponding to each observation height.

[0054] Specifically, during the phase transformation process, the horizontal x-component, horizontal y-component, and vertical z-component of the three-component magnetic anomaly data are subjected to corresponding phase transformation processing to ensure that the three transformed components all satisfy the magnetic field distribution characteristics of vertical magnetization. At the same time, during the frequency domain processing, low-pass filtering can be used to suppress the noise amplification problem in the low-frequency band and improve the stability of the phase transformation results, which is especially suitable for magnetic anomaly data processing in low geomagnetic latitude regions.

[0055] Thus, by performing independent phase transformation processing on the three components of magnetic anomaly data at each observation height, this invention effectively eliminates the distortion effect of the oblique geomagnetic environment on the distribution pattern of magnetic anomalies. This allows the converted equivalent vertical magnetization magnetic anomaly data to directly correspond to the spatial distribution characteristics of the fracture diffuser, significantly reducing the difficulty of subsequent magnetic anomaly inversion while ensuring the consistency of data at different observation heights.

[0056] S103. Based on each equivalent vertical magnetization anomaly data, extract the zero intersection point location information at each observation height.

[0057] Optionally, a symmetry axis analysis can be performed on each equivalent vertical magnetization anomaly data to determine the symmetry axis of the total magnetic field anomaly.

[0058] In this embodiment of the invention, the axis of symmetry of the total magnetic field anomaly corresponds to the axis of symmetry of the direction or tendency of the crack diffuser after grouting. Under vertical magnetization conditions, the distribution pattern of the total magnetic field anomaly is symmetrically distributed with the axis as the center. The magnetic measurement profile selected along the axis can completely reflect the variation law of the magnetic anomaly with the horizontal position and is the optimal profile direction for extracting the zero intersection position.

[0059] The direction of the symmetry axis of the total magnetic field anomaly is the direction of the preset survey line.

[0060] Optionally, the preset measurement line direction can be initially set based on the layout location of the grouting borehole and the expected direction of the fracture, and then corrected by symmetry analysis of the measured total magnetic field anomaly data, and finally the direction of the symmetry axis of the total magnetic field anomaly is determined.

[0061] For example, in the actual symmetry axis analysis process, the expected direction and tendency of the grouting cracks are first determined based on the design data of the grouting project, and the initial direction of the preset survey line is initially set; then, the total magnetic field anomaly data at each observation height is analyzed in two-dimensional plane to find the symmetry center and symmetry axis of the magnetic anomaly distribution, the axis is determined as the symmetry axis of the total magnetic field anomaly, and the direction of the axis is corrected to the final preset survey line direction.

[0062] Specifically, for the equivalent vertical magnetization magnetic anomaly data at each observation height, symmetry axis analysis was performed to verify the consistency of the symmetry axis of the total magnetic field anomaly at different observation heights. If there is a slight deviation in the symmetry axis at different observation heights, it was unified based on the symmetry axis of the reference observation surface to ensure that the one-dimensional magnetic profiles selected at different observation heights are located in the same vertical plane, thus ensuring that the zero intersection point position information has spatial correspondence.

[0063] Furthermore, a one-dimensional magnetic profile is selected along the axis of symmetry of the total magnetic field anomaly.

[0064] For example, after determining the axis of symmetry of the total magnetic field anomaly, a survey line covering the entire distribution range of the magnetic anomaly is selected along the direction of the axis of symmetry. Both ends of the survey line need to extend to the position where the magnetic anomaly signal returns to the level of the background magnetic field, ensuring that the survey line can completely include the rising segment, peak segment and falling segment of the magnetic anomaly. Then, along the survey line, the equivalent vertical magnetization magnetic anomaly component data corresponding to each measuring point is extracted to form a one-dimensional magnetic profile dataset.

[0065] Specifically, for the equivalent vertical magnetization magnetic anomaly data at each observation height, a one-dimensional magnetic profile is extracted along the same preset survey line direction to obtain a one-dimensional magnetic profile dataset corresponding to the reference observation surface, the first observation surface, and the second observation surface, respectively. Each profile dataset contains the correspondence between the horizontal x-component and the vertical z-component magnetic anomaly data and the horizontal position.

[0066] Furthermore, on the one-dimensional magnetic profile, the zero-intersection position information at each observation height is extracted.

[0067] Among them, the zero intersection point refers to the position point where the value of the magnetic anomaly component changes from positive to negative or from negative to positive, that is, the position point where the value of the magnetic anomaly component is 0. Under vertical magnetization conditions, the position of the zero intersection point of the magnetic anomaly component of the thin magnetic plate has a fixed analytical relationship with the burial depth, tilt angle and horizontal distribution scale of the magnetic plate.

[0068] In one alternative implementation, the zero point position of the horizontal component in the equivalent vertical magnetization magnetic anomaly data at each observation height, where it changes from positive to negative or from negative to positive along a preset survey line direction, can be determined as the horizontal zero intersection point.

[0069] In this embodiment of the invention, the horizontal component is the horizontal x component along the preset survey line direction in the equivalent vertical magnetization magnetic anomaly data. On the one-dimensional magnetic survey profile, the zero intersection position of this component is linearly shifted with the change of observation height, and its shift is directly related to the burial depth of the fracture diffuser.

[0070] For example, when extracting the horizontal zero intersection point, first extract the magnetic anomaly values ​​of the horizontal x component and the corresponding sequence data of the horizontal position from the one-dimensional magnetic profile dataset; then traverse the sequence data to find the interval where the magnetic anomaly value changes from positive to negative or from negative to positive, and determine the horizontal position where the magnetic anomaly value is exactly 0 through linear interpolation calculation. This position is the horizontal zero intersection point at the observation height.

[0071] Specifically, for the one-dimensional magnetic profile data of the reference observation plane, the first observation plane, and the second observation plane, the horizontal position coordinates of the horizontal zero intersection point corresponding to each observation height are calculated to form a multi-height horizontal zero intersection point position sequence, which is used for subsequent analysis of the rate of change of the horizontal zero intersection point position at different observation heights.

[0072] In another alternative implementation, the positions of two symmetrical zero points of the vertical component of the equivalent vertical magnetization anomaly data along the preset survey line direction can be determined at each observation height, as a pair of vertical zero intersection points.

[0073] In this embodiment of the invention, the vertical component is the vertical z-component in the equivalent vertical magnetization magnetic anomaly data. On the one-dimensional magnetic profile, this component will have two symmetrically distributed zero intersection points. The distance between the two zero intersection points changes with the observation height, and the amount of change is directly related to the inclination angle of the fracture diffuser and the horizontal projection scale.

[0074] For example, when extracting vertical zero intersections, first extract the magnetic anomaly values ​​of the vertical z component and the corresponding horizontal position sequence data from the one-dimensional magnetic profile dataset; then traverse the sequence data to find the two intervals where the magnetic anomaly values ​​change from positive to negative and from negative to positive, respectively. Through linear interpolation, determine the two horizontal positions where the magnetic anomaly values ​​are 0. These two positions are a pair of vertical zero intersections at the observation height.

[0075] Specifically, for the one-dimensional magnetic profile data of the reference observation plane, the first observation plane, and the second observation plane, the horizontal position coordinates of a pair of vertical zero intersection points corresponding to each observation height are calculated, and the horizontal distance between the two vertical zero intersection points is calculated to form a multi-height vertical zero intersection point distance sequence, which is used for subsequent analysis of the rate of change of the vertical zero intersection point distance with the observation height.

[0076] Thus, by selecting a one-dimensional magnetic profile along the symmetry axis of the total magnetic field anomaly, this invention extracts the position information of the horizontal and vertical zero intersection points at each observation height, transforming the complex magnetic anomaly field distribution into quantifiable geometric feature points. This establishes a direct correlation between the magnetic anomaly features and the geometric parameters of the crack diffuser, significantly simplifying the complexity of subsequent inversion calculations while ensuring the analytical accuracy of the inversion calculations.

[0077] S104. Based on the change in the position information of each zero intersection point, invert and calculate the spatial geometric parameters of the fracture diffuser.

[0078] In this embodiment of the invention, before the inversion calculation, it is necessary to construct an equivalent thin plate model of the grouting fracture. The rock mass fracture after grouting and filling with magnetic grouting material is equivalent to a thin magnetic plate model.

[0079] Among them, the thin-layer magnetic plate model ignores the influence of the crack thickness on the grouting diffusion scale, and the crack diffuser is equivalent to a thin plate geometric model with burial depth, inclination angle and horizontal distribution scale. The magnetization intensity of the model is consistent with the magnetization intensity of the crack diffuser after curing, and the distribution direction of the model matches the direction and tendency of the crack after grouting.

[0080] In this embodiment of the invention, the change in the zero-intersection position information includes the rate of change of the horizontal zero-intersection position at different observation heights and the rate of change of the vertical zero-intersection spacing at different observation heights. Based on the magnetic anomaly forward modeling theory of thin-layer magnetic plate model, the above-mentioned changes have a one-to-one analytical relationship with the burial depth, tilt angle, and horizontal projection scale of the fracture diffuser. Through this analytical relationship, the quantitative inversion of geometric parameters can be achieved.

[0081] In this embodiment of the invention, the spatial geometric parameters include at least the burial depth of the fracture diffuser.

[0082] Among them, the burial depth is the center burial depth of the fracture diffuser, that is, the vertical distance between the center position of the fracture diffuser and the ground surface reference observation surface, reflecting the spatial vertical position of the grouting fracture.

[0083] Optionally, the spatial geometric parameters may also include the inclination angle and horizontal projection scale of the fracture diffuser.

[0084] The inclination angle is the angle between the plane containing the fracture diffuser and the horizontal plane, reflecting the spatial tilt of the grouting fracture; the horizontal projection scale is the horizontal projection length of the fracture diffuser along the preset survey line direction, reflecting the diffusion range of the grouting fracture.

[0085] In some embodiments, the burial depth, dip angle, and horizontal projection scale of the fracture diffuser can be calculated analytically based on the rate of change of the position of each horizontal zero intersection point and the rate of change of the spacing between each vertical zero intersection point with the observation height.

[0086] In some embodiments, the inversion calculation adopts a step-by-step solution method. First, the burial depth parameters are inverted based on the rate of change of the horizontal zero intersection position. Then, the obtained burial depth parameters are substituted into the analytical relationship, and the dip angle parameters and horizontal projection scale parameters are inverted in sequence by combining the rate of change of the vertical zero intersection distance.

[0087] For example, in the actual inversion calculation process, firstly, based on the horizontal zero-intersection positions of the reference observation surface, the first observation surface, and the second observation surface, the offset of the horizontal zero-intersection position between different observation heights is calculated, thereby obtaining the rate of change of the horizontal zero-intersection position with the observation height; then, based on this rate of change, the burial depth parameter of the fracture diffuser is calculated through the analytical relationship of the thin-layer magnetic plate model; after substituting the burial depth parameter into the analytical relationship, combined with the rate of change of the vertical zero-intersection spacing at different observation heights, the inclination parameter of the fracture diffuser is calculated; finally, the burial depth parameter and the inclination parameter are substituted into the analytical relationship to calculate the horizontal projection scale parameter of the fracture diffuser.

[0088] Specifically, after completing the inversion calculation of the burial depth, inclination angle, and horizontal projection scale parameters of the fracture diffuser, the three-dimensional spatial morphology of the grouting fracture can be reconstructed by combining the thin-layer magnetic plate model, and the spatial location, burial depth range, distribution direction and diffusion scale of the grouting fracture can be determined to form a complete spatial positioning result of the grouting fracture.

[0089] Thus, this invention performs step-by-step inversion calculations by varying the position of zero intersection points at multiple heights, making full use of the depth-sensitive information introduced by multi-height observations, breaking through the accuracy limitations of single-observation height inversions, and achieving quantitative solutions for the burial depth, dip angle, and horizontal distribution scale of grouting fractures.

[0090] In the spatial location method for grouting fractures provided by this invention, by acquiring three-component magnetic anomaly data at at least three different observation heights, the magnetic anomaly information at a single observation height is expanded into a multi-height observation system containing depth-sensitive information; through phase transformation processing, the complex magnetic anomaly under oblique magnetization conditions is converted into an equivalent vertical magnetization magnetic anomaly, eliminating the interference of geomagnetic field direction on inversion accuracy; by extracting the zero-intersection point location information at each observation height, the spatial distribution characteristics of the magnetic anomaly are transformed into quantifiable geometric feature points; by inverting and calculating spatial geometric parameters based on the changes in the zero-intersection point location information at each observation height, the joint quantitative solution of burial depth, dip angle, and horizontal projection scale is achieved using the difference information between multiple heights, thereby realizing non-destructive and quantitative location of grouting fractures.

[0091] The grouting fracture spatial positioning device provided by the present invention is described below. The grouting fracture spatial positioning device described below and the grouting fracture spatial positioning method described above can be referred to in correspondence.

[0092] Figure 3 This is a structural diagram of a grouting fracture spatial positioning device provided in an embodiment of the present invention. The grouting fracture spatial positioning device includes: a data acquisition module 301, a phase conversion module 302, a zero intersection point extraction module 303, and a joint inversion module 304.

[0093] The system comprises: a data acquisition module 301, used to acquire three-component magnetic anomaly data generated by fracture diffusers at at least three different observation heights, wherein the fracture diffusers are formed by grouting in the rock mass; a phase transformation module 302, used to perform phase transformation processing on each three-component magnetic anomaly data to obtain equivalent vertical magnetization anomaly data corresponding to each observation height; a zero-intersection extraction module 303, used to extract the zero-intersection position information at each observation height based on each equivalent vertical magnetization anomaly data; and a joint inversion module 304, used to invert and calculate the spatial geometric parameters of the fracture diffusers based on the change in the position information of each zero-intersection, wherein the spatial geometric parameters include at least the burial depth.

[0094] In some embodiments, the zero-intersection extraction module 303 is specifically used for: performing symmetry axis analysis on each equivalent vertical magnetization magnetic anomaly data to determine the symmetry axis of the total magnetic field anomaly, the direction of which is a preset survey line direction; selecting a one-dimensional magnetic survey profile along the symmetry axis of the total magnetic field anomaly; and extracting the zero-intersection position information at each observation height on the one-dimensional magnetic survey profile.

[0095] In some embodiments, the zero-intersection extraction module 303 is specifically used to: determine the zero-point position of the horizontal component in the equivalent vertical magnetization magnetic anomaly data at each observation height, where the horizontal component changes from positive to negative or from negative to positive along the preset survey line direction, as a horizontal zero-intersection point; and determine the two symmetrical zero-point positions of the vertical component in the equivalent vertical magnetization magnetic anomaly data at each observation height, along the preset survey line direction, as a pair of vertical zero-intersection points.

[0096] In some embodiments, the joint inversion module 304 is specifically used to: calculate the burial depth, dip angle and horizontal projection scale of the fracture diffuser by analytical relationships based on the rate of change of the position of each horizontal zero intersection point and the rate of change of the spacing between each vertical zero intersection point with the observation height.

[0097] In some embodiments, the phase conversion module 302 is specifically used to: convert each three-component magnetic anomaly data into magnetic anomaly data under vertical magnetization conditions through frequency domain transformation to obtain equivalent vertical magnetization magnetic anomaly data; wherein, the frequency domain transformation adopts Fourier transform, and the same geomagnetic field parameters are used to construct the phase conversion operator at each observation height.

[0098] In the grouting fracture spatial positioning device provided by this invention, by acquiring three-component magnetic anomaly data at at least three different observation heights, the magnetic anomaly information at a single observation height is expanded into a multi-height observation system containing depth-sensitive information; by using phase transformation processing, the complex magnetic anomaly under oblique magnetization conditions is converted into an equivalent vertical magnetization magnetic anomaly, eliminating the interference of the geomagnetic field direction on the inversion accuracy; by extracting the zero-intersection point position information at each observation height, the spatial distribution characteristics of the magnetic anomaly are transformed into quantifiable geometric feature points; by inverting and calculating the spatial geometric parameters based on the changes in the zero-intersection point position information at each observation height, the joint quantitative solution of burial depth, dip angle, and horizontal projection scale is realized by utilizing the difference information between multiple heights, thereby achieving non-destructive and quantitative positioning of the grouting fracture spatial location.

[0099] Figure 4 An example is a schematic diagram of the physical structure of an electronic device, such as... Figure 4As shown, the electronic device may include: a processor 410, a communications interface 420, a memory 430, and a communications bus 440, wherein the processor 410, the communications interface 420, and the memory 430 communicate with each other through the communications bus 440.

[0100] The processor 410 can call logic instructions in the memory 430 to execute a grouting fracture spatial localization method. This method includes: parsing the device configuration file of the original network to obtain the topology and end-to-end forwarding behavior information of the original network; acquiring three-component magnetic anomaly data generated by the fracture diffuser at at least three different observation heights; performing phase transformation processing on each three-component magnetic anomaly data to obtain equivalent vertical magnetization magnetic anomaly data corresponding to each observation height; extracting zero-intersection point location information at each observation height based on each equivalent vertical magnetization magnetic anomaly data; and calculating the spatial geometric parameters of the fracture diffuser based on the change in each zero-intersection point location information, wherein the spatial geometric parameters include at least the burial depth.

[0101] Furthermore, the logical instructions in the aforementioned memory 430 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, essentially, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0102] On the other hand, the present invention also provides a computer program product, which includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the grouting fracture spatial positioning method provided by the above methods. The method includes: acquiring three-component magnetic anomaly data generated by the fracture diffuser at at least three different observation heights; performing phase transformation processing on each three-component magnetic anomaly data to obtain equivalent vertical magnetization magnetic anomaly data corresponding to each observation height; extracting zero-intersection point position information at each observation height based on each equivalent vertical magnetization magnetic anomaly data; and inverting and calculating the spatial geometric parameters of the fracture diffuser according to the change in each zero-intersection point position information, wherein the spatial geometric parameters include at least the burial depth.

[0103] In another aspect, the present invention also provides a non-transitory computer-readable storage medium storing a computer program thereon. When executed by a processor, the computer program performs the grouting fracture spatial positioning method provided by the above methods. The method includes: acquiring three-component magnetic anomaly data generated by the fracture diffuser at at least three different observation heights; performing phase transformation processing on each three-component magnetic anomaly data to obtain equivalent vertical magnetization magnetic anomaly data corresponding to each observation height; extracting zero-intersection point position information at each observation height based on each equivalent vertical magnetization magnetic anomaly data; and inverting and calculating the spatial geometric parameters of the fracture diffuser based on the change in each zero-intersection point position information, wherein the spatial geometric parameters include at least the burial depth.

[0104] The system embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0105] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0106] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for spatially locating grouting fractures, characterized in that, The method includes: Acquire three-component magnetic anomaly data generated by the fracture diffuser, collected at at least three different observation heights; Each of the three-component magnetic anomaly data is subjected to phase transformation processing to obtain the equivalent vertical magnetization magnetic anomaly data corresponding to each observation height; Based on each of the equivalent vertical magnetization magnetic anomaly data, the zero intersection point location information at each observation height is extracted; Based on the change in the location information of each of the zero intersection points, the spatial geometric parameters of the fracture diffuser are calculated by inversion, and the spatial geometric parameters include at least the burial depth.

2. The method according to claim 1, characterized in that, The fissure diffuser is formed by solidifying a magnetic grouting material, which is a mixture of cement-based grouting material and fine-grained magnetite powder.

3. The method according to claim 1, characterized in that, The at least three different observation altitudes include: a reference observation surface located on the ground surface, a first observation surface located at a first altitude, and a second observation surface located at a second altitude; Wherein, the second height is greater than the first height, and the projections of the reference observation surface, the first observation surface, and the second observation surface on the horizontal plane coincide with each other.

4. The method according to claim 1, characterized in that, The step of extracting the zero-intersection location information at each observation height based on each equivalent vertical magnetization anomaly data includes: A symmetry axis analysis is performed on each of the equivalent vertical magnetization magnetic anomaly data to determine the total magnetic field anomaly symmetry axis, the direction of which is the preset survey line direction; A one-dimensional magnetic profile is selected along the axis of symmetry of the total magnetic field anomaly. On the one-dimensional magnetic profile, the zero-intersection position information at each observation height is extracted.

5. The method according to claim 4, characterized in that, The extraction of the zero-intersection position information at each observation height includes: At each observation height, the zero point position of the horizontal component in the equivalent vertical magnetization anomaly data, which changes from positive to negative or from negative to positive along the preset survey line direction, is determined as the horizontal zero intersection point. At each observation height, determine the positions of two symmetrical zero points of the vertical component in the equivalent vertical magnetization anomaly data along the preset survey line direction, and use them as a pair of vertical zero intersection points.

6. The method according to claim 5, characterized in that, The step of inverting and calculating the spatial geometric parameters of the fracture diffuser based on the change in the position information of each of the zero intersection points includes: Based on the rate of change of the position of each of the horizontal zero intersection points and the rate of change of the distance between each of the vertical zero intersection points with the observation height, the burial depth, inclination angle and horizontal projection scale of the fracture diffuser are calculated through analytical relationships.

7. The method according to claim 1, characterized in that, The step of performing phase transformation processing on each of the three-component magnetic anomaly data to obtain the equivalent vertical magnetization magnetic anomaly data corresponding to each observation height includes: Each of the three-component magnetic anomaly data is converted into magnetic anomaly data under vertical magnetization conditions through frequency domain transformation to obtain the equivalent vertical magnetization magnetic anomaly data. The frequency domain transformation employs Fourier transform, and the same geomagnetic field parameters are used to construct the phase transformation operator at each observation altitude.

8. A grouting fracture spatial positioning device, characterized in that, The device includes: The data acquisition module is used to acquire three-component magnetic anomaly data generated by fracture diffusers from at least three different observation heights, wherein the fracture diffusers are formed by grouting in the rock mass; The phase transformation module is used to perform phase transformation processing on each of the three-component magnetic anomaly data to obtain the equivalent vertical magnetization magnetic anomaly data corresponding to each observation height. The zero-intersection extraction module is used to extract the zero-intersection position information at each observation height based on each of the equivalent vertical magnetization magnetic anomaly data. The joint inversion module is used to invert and calculate the spatial geometric parameters of the fracture diffuser based on the change in the position information of each of the zero intersection points, wherein the spatial geometric parameters include at least the burial depth.

9. An electronic device, characterized in that, The method includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, it implements the method as described in any one of claims 1 to 6.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the method as described in any one of claims 1 to 6.

Images