Geological fracture exploration and detection method and system for geological engineering

Abstract

The present invention relates to the field of gravity exploration technology, and more specifically to a method and system for geological fracture exploration and detection for geological engineering. The method comprises: obtaining gravity data corresponding to a point to be detected and gravity data corresponding to neighboring points in each preset direction; determining an initial fracture rate corresponding to each preset direction and its opposite direction; determining a reference usable index corresponding to each preset direction; determining a target fracture rate corresponding to the point to be detected; if the target fracture rate corresponding to the point to be detected is greater than a preset fracture threshold, determining the target fracture direction based on the distribution of the initial fracture rates corresponding to all preset directions; determining a final fracture rate corresponding to the point to be detected based on gravity changes in the target fracture direction, and judging whether the point to be detected is a fracture source location based on the final fracture rate. The present invention achieves geological fracture exploration and detection through gravity exploration, and improves the accuracy of geological fracture exploration and detection.

Description

Technical Field

[0001] The present invention relates to the field of gravity exploration technology, and in particular to a geological fracture exploration and detection method and system for geological engineering. Background Art

[0002] Geological fractures are structural defects in rocks or underground formations caused by external stresses, resulting in cracks or fractures. Their detection plays a crucial role in mineral resource development, groundwater management, oil and gas extraction, and earthquake early warning. It also provides a scientific basis for engineering construction and environmental protection. Therefore, the exploration and detection of geological fractures for geological engineering is crucial.

[0003] Generally, gravity exploration is used to explore geological fractures. Specifically, when the gravity value of the exploration point is abnormal compared with the gravity values ​​of other exploration points around it, the exploration point is considered to be a fracture location point. Among them, the exploration point is also called the detection point. When the difference between the gravity value of the exploration point and the gravity values ​​of other exploration points around it is large, it is often considered that the gravity value of the exploration point is abnormal compared with the gravity values ​​of other exploration points around it. The fracture source location point can be the source location point where the fracture occurs, which is a type of fracture location point. However, when judging whether the exploration point is a fracture source location point based on whether the gravity value of the exploration point is abnormal, the following technical problems often occur:

[0004] Since underground cracks are usually accompanied by changes in the physical properties of rock formations, especially density differences, if there are areas of low or high density underground in a certain area, such as cracks, voids, and faults, it will often lead to abnormal changes in gravity values. Therefore, when the gravity value of the exploration point is abnormal compared with the gravity values ​​of other surrounding exploration points, the exploration point is not necessarily the crack source location point, it may also be a void location point or a fault location point. At this time, if the exploration point is directly determined to be the crack source location point, it may lead to misjudgment of the crack source location point, resulting in poor accuracy in geological crack exploration and detection. Summary of the Invention

[0005] In order to solve the technical problem of poor accuracy in geological fracture exploration and detection, the present invention proposes a geological fracture exploration and detection method and system for geological engineering.

[0006] In a first aspect, the present invention provides a geological fracture exploration and detection method for geological engineering, the method comprising:

[0007] Obtain the gravity data corresponding to the point to be detected and the gravity data corresponding to each of its neighboring points in a preset direction;

[0008] Based on the difference between the gravity data corresponding to the point to be detected and the gravity data corresponding to its neighboring points in each preset direction, the initial crack rate corresponding to each preset direction is determined. Similarly, the initial crack rate corresponding to the opposite direction of each preset direction is determined.

[0009] Determining a reference usable index corresponding to each preset direction according to a difference between an initial crack rate corresponding to each preset direction and an initial crack rate corresponding to the opposite direction;

[0010] Determine the target crack rate corresponding to the point to be tested based on the initial crack rates corresponding to all preset directions and the reference available indicators;

[0011] If the target crack rate corresponding to the point to be detected is greater than the preset crack threshold, the target crack direction is determined based on the distribution of the initial crack rates corresponding to all preset directions;

[0012] According to the gravity change in the target crack direction, the final crack rate corresponding to the detection point is determined, and based on the final crack rate, it is judged whether the detection point is the crack source location point.

[0013] In combination with the first aspect above, in one possible implementation, determining the initial crack rate corresponding to each preset direction based on the difference between the gravity data corresponding to the to-be-detected point and the gravity data corresponding to its neighboring points in each preset direction includes:

[0014] Determine any preset direction as a marked direction, and determine the absolute value of the difference between the gravity data corresponding to the point to be detected and the gravity data corresponding to the neighboring points in the marked direction as the target gravity difference corresponding to the marked direction;

[0015] An initial crack rate corresponding to the marked direction is determined according to a target gravity difference corresponding to the marked direction, wherein the target gravity difference is positively correlated with the initial crack rate.

[0016] In conjunction with the first aspect above, in one possible implementation, determining the reference usable indicator corresponding to each preset direction based on the difference between the initial crack rate corresponding to each preset direction and the initial crack rate corresponding to the opposite direction thereof includes:

[0017] Determine any preset direction as a marking direction, and determine the opposite direction of the marking direction as a reference direction;

[0018] Determine the absolute value of the difference between the initial crack rate corresponding to the marked direction and the initial crack rate corresponding to the reference direction as the possible change difference of the crack corresponding to the marked direction;

[0019] According to the possible crack change difference corresponding to the marked direction, a reference usable index corresponding to the marked direction is determined, wherein the possible crack change difference is negatively correlated with the reference usable index.

[0020] In combination with the first aspect above, in a possible implementation, determining the target crack rate corresponding to the point to be detected based on the initial crack rates corresponding to all preset directions and the reference available index includes:

[0021] Determine the accumulated values ​​of the reference available indicators corresponding to all preset directions as the available representative indicator;

[0022] Determine the target available weight corresponding to each preset direction by the ratio of the reference available index corresponding to each preset direction to the available representative index;

[0023] The product of the target available weight corresponding to each preset direction and the initial crack rate is determined as the target sub-crack factor corresponding to each preset direction;

[0024] The target crack rate corresponding to the point to be detected is determined according to the target sub-crack factors corresponding to all preset directions, wherein the target sub-crack factor is positively correlated with the target crack rate.

[0025] In conjunction with the first aspect above, in a possible implementation, determining the target crack direction according to the distribution of the initial crack rates corresponding to all preset directions includes:

[0026] According to the distribution of initial crack rates corresponding to all preset directions, a crack rate variation function and its error sum function are constructed;

[0027] The crack rate variation function is solved by the least squares method, and the solution with the minimum error and function is taken as the final solution to obtain the target crack direction.

[0028] In combination with the first aspect above, in a possible implementation, constructing a crack rate variation function and its error sum function based on the distribution of initial crack rates corresponding to all preset directions includes:

[0029] The value range of is [0°, 180°]. When it belongs to [0°, 90°], the crack rate variation function is:

[0030] ;

[0031] ;

[0032] ;in, is the crack rate change function value when the independent variable is a, which is equal to the fitted value of the initial crack rate corresponding to the preset direction with the corresponding angle a; a is the angle corresponding to different preset directions; k is the slope of the fitted line contained in the crack rate change function; 、 and is the intercept of different fitting straight lines included in the crack rate variation function;

[0033] when When it belongs to (90°, 180°], the crack rate variation function is:

[0034] ;

[0035] ;

[0036] ;

[0037] The error sum function of the crack rate variation function is:

[0038] ; Wherein, Q is the error and function value of the crack rate variation function; N is the number of preset directions; i is the serial number of the preset direction; is the reference available indicator corresponding to the i-th preset direction; is the crack rate variation function value when the independent variable takes the angle corresponding to the i-th preset direction; is the initial crack rate corresponding to the i-th preset direction.

[0039] In conjunction with the first aspect above, in a possible implementation, determining the final crack rate corresponding to the to-be-detected point based on the gravity change in the target crack direction includes:

[0040] Obtaining gravity data corresponding to a preset number of preset position points in the target crack direction of the point to be detected and the neighboring points of these preset position points in each preset direction;

[0041] Determine the absolute value of the difference between the gravity data corresponding to each preset position point and the gravity data corresponding to its neighboring points in each preset direction as the temporary gravity difference of each preset position point in each preset direction;

[0042] Determine the cumulative value of the temporary gravity difference of each preset position point in all preset directions as the reference gravity difference corresponding to each preset position point;

[0043] The final crack rate corresponding to the point to be detected is determined according to the reference gravity differences corresponding to a preset number of preset position points in the target crack direction.

[0044] In combination with the first aspect above, in a possible implementation, determining the final crack rate corresponding to the to-be-detected point based on the reference gravity differences corresponding to a preset number of preset position points in the target crack direction includes:

[0045] Determine the distance between the point to be detected and each preset position point as the target distance corresponding to each preset position point;

[0046] The final crack rate corresponding to the point to be detected is determined based on the target crack rate corresponding to the point to be detected, and the reference gravity difference and target distance corresponding to a preset number of preset position points in the target crack direction, wherein the target crack rate and the reference gravity difference are both positively correlated with the final crack rate, and the target distance is negatively correlated with the final crack rate.

[0047] In conjunction with the first aspect above, in a possible implementation, judging whether the to-be-detected point is a crack source location point based on the final crack rate includes:

[0048] If the final crack rate corresponding to the to-be-detected point is greater than a preset abnormal threshold, the to-be-detected point is determined to be a crack source location point.

[0049] In a second aspect, the present invention provides a geological fracture exploration and detection system for geological engineering, comprising a processor and a memory, wherein the processor is configured to process instructions stored in the memory to implement a geological fracture exploration and detection method for geological engineering. Specifically, the system comprises:

[0050] A gravity data acquisition module is used to acquire the gravity data corresponding to the point to be detected and the gravity data corresponding to the neighboring points in each preset direction;

[0051] An initial crack rate determination module is used to determine the initial crack rate corresponding to each preset direction based on the difference between the gravity data corresponding to the point to be detected and the gravity data corresponding to its neighboring points in each preset direction. Similarly, the initial crack rate corresponding to the opposite direction of each preset direction is determined;

[0052] A reference available index determination module is used to determine a reference available index corresponding to each preset direction based on a difference between an initial crack rate corresponding to each preset direction and an initial crack rate corresponding to the opposite direction;

[0053] A target crack rate determination module is used to determine the target crack rate corresponding to the point to be detected based on the initial crack rates corresponding to all preset directions and reference available indicators;

[0054] A target crack direction determination module is used to determine the target crack direction based on the distribution of initial crack rates corresponding to all preset directions if the target crack rate corresponding to the point to be detected is greater than a preset crack threshold;

[0055] The determination and judgment module is used to determine the final crack rate corresponding to the point to be detected according to the gravity change in the target crack direction, and based on the final crack rate, determine whether the point to be detected is the crack source location point.

[0056] In a third aspect, a server is provided, comprising a memory and a processor. The memory is configured to store executable program code, and the processor is configured to call and execute the executable program code from the memory, so that the device executes the method of the first aspect or any possible implementation of the first aspect.

[0057] In a fourth aspect, a computer program product is provided, comprising: a computer program code, which, when executed on a computer, enables the computer to execute the method in the first aspect or any possible implementation of the first aspect.

[0058] In a fifth aspect, a computer-readable storage medium is provided, which stores a computer program code. When the computer program code runs on a computer, the computer executes the method in the above-mentioned first aspect or any possible implementation of the first aspect.

[0059] The present invention has the following beneficial effects:

[0060] The geological crack exploration and detection method for geological engineering of the present invention realizes geological crack exploration and detection through gravity exploration, solves the technical problem of poor accuracy of geological crack exploration and detection, and improves the accuracy of geological crack exploration and detection. Specifically, when performing geological crack exploration and detection, the present invention not only considers the difference between the gravity data corresponding to the point to be detected and the gravity data corresponding to the neighboring points in each preset direction, but also comprehensively considers multiple features related to the law of geological crack change, such as reference available indicators, target crack rate and gravity change in the target crack direction, thereby quantifying the final crack rate corresponding to the point to be detected, and based on the final crack rate, it can relatively objectively judge whether the point to be detected is a crack source location point, which can improve the accuracy of geological crack exploration and detection to a certain extent. Among them, the crack source location point can be the source location point where the crack is generated, which is a type of crack location point. BRIEF DESCRIPTION OF THE DRAWINGS

[0061] In order to more clearly illustrate the technical solutions and advantages of the embodiments of the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the prior art descriptions. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0062] Figure 1 This is a flow chart of the geological fracture exploration and detection method for geological engineering of the present invention;

[0063] Figure 2 Schematic diagram of the distribution of the point to be detected and multiple neighboring points of the present invention;

[0064] Figure 3 This is a schematic diagram of the structure of the geological fracture exploration and detection system for geological engineering of the present invention;

[0065] Figure 4 A schematic structural diagram of a computer device according to the present invention;

[0066] Among them, the reference numerals include: 201, a first solid dot; 202, a second solid dot; 203, a third solid dot; 204, a fourth solid dot. DETAILED DESCRIPTION

[0067] To further illustrate the technical means and effects employed by the present invention to achieve its intended objectives, the following, in conjunction with the accompanying drawings and preferred embodiments, describes in detail the specific implementations, structures, features, and effects of the technical solutions proposed by the present invention. In the following description, references to "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. Furthermore, specific features, structures, or characteristics of one or more embodiments may be combined in any suitable manner.

[0068] Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.

[0069] refer to Figure 1 , showing the process of some embodiments of the geological fracture exploration and detection method for geological engineering of the present invention. The geological fracture exploration and detection method for geological engineering includes the following steps:

[0070] Step S1: Obtain gravity data corresponding to the point to be detected and gravity data corresponding to neighboring points in each preset direction.

[0071] The point to be detected may be a location where crack detection is to be performed. The preset direction may be a pre-set direction, and the corresponding angle may be in the range of [0°, 180°]. The number of preset directions may be pre-set, for example, the number of preset directions may be 19. It should be noted that the more preset directions are set, the higher the accuracy of subsequent crack detection. If the angle corresponding to the preset direction is in the range of [0°, 180°], and the number of preset directions is 19, then the 19 preset directions may be, in order: 0° direction, 10° direction, 20° direction, 30° direction, 40° direction, 50° direction, 60° direction, 70° direction, 80° direction, 90° direction, 100° direction, 110° direction, 120° direction, 130° direction, 140° direction, 150° direction, 160° direction, 170° direction, and 180° direction. The neighborhood point may be a location point pre-set around the point to be detected, and the distances between different neighborhood points and the point to be detected may be the same. For example, the distance between the neighboring point and the point to be detected may be 10 meters. The gravity data may be a gravity value collected by a gravity measuring instrument. The gravity measuring instrument is mainly used to collect gravity, and may be a gravimeter.

[0072] As an example, Figure 2 As shown, the first solid dot 201 can represent the point to be detected, the second solid dot 202 can represent the neighboring point in the 0° direction of the point to be detected, the third solid dot 203 can represent the neighboring point in the 90° direction of the point to be detected, and the fourth solid dot 204 can represent the neighboring point in the 180° direction of the point to be detected; the distance between the first solid dot 201 and the second solid dot 202 can be equal to the distance between the first solid dot 201 and the third solid dot 203; the distance between the first solid dot 201 and the third solid dot 203 can be equal to the distance between the first solid dot 201 and the third solid dot 203. 03 can be equal to the distance between the first solid dot 201 and the fourth solid dot 204; the gravity data corresponding to the point to be detected, the gravity data corresponding to the neighboring points in the 0° direction of the point to be detected, the gravity data corresponding to the neighboring points in the 90° direction of the point to be detected, and the gravity data corresponding to the neighboring points in the 180° direction of the point to be detected can be collected in sequence by installing gravimeters at the position points represented by the first solid dot 201, the second solid dot 202, the third solid dot 203 and the fourth solid dot 204.

[0073] Step S2, based on the difference between the gravity data corresponding to the point to be detected and the gravity data corresponding to its neighboring points in each preset direction, determine the initial crack rate corresponding to each preset direction. Similarly, determine the initial crack rate corresponding to the opposite direction of each preset direction.

[0074] The opposite direction of the preset direction may be a direction completely opposite to the preset direction on the same straight line. For example, the opposite direction of the 0° direction may be the 180° direction. The opposite direction of the 10° direction may be the 190° direction.

[0075] As an example, this step may include the following steps:

[0076] In the first step, any preset direction is determined as the marked direction, and the absolute value of the difference between the gravity data corresponding to the above-mentioned point to be detected and the gravity data corresponding to the neighboring point on the above-mentioned marked direction is determined as the target gravity difference corresponding to the above-mentioned marked direction.

[0077] The second step is to determine the initial crack rate corresponding to the above-marked direction based on the target gravity difference corresponding to the above-marked direction.

[0078] Among them, the target gravity difference can be positively correlated with the initial crack rate.

[0079] It should be noted that because geological fractures are often accompanied by the alternation of different materials, such as the alternating distribution of softer rock layers, harder rock layers, and different minerals, the density differences of these materials often lead to changes in the gravity field at different locations. Therefore, the gravity at the fracture location often differs from the gravity at the surrounding locations. Therefore, the greater the initial fracture rate corresponding to the marked direction, the greater the relative gravity difference between the detection point and the neighboring points in the marked direction. This often indicates that the gravity value at the detection point is more likely to have a larger gravity anomaly than other surrounding gravity exploration points, and the greater the possibility that the detection point is a fracture location.

[0080] For example, the formula for determining the initial crack rate corresponding to the preset direction can be:

[0081] ;in, is the initial crack rate corresponding to the i-th preset direction. i is the serial number of the preset direction. It is the absolute value function. is the gravity data corresponding to the neighborhood point in the i-th preset direction. G is the gravity data corresponding to the point to be detected. is the target gravity difference corresponding to the i-th preset direction.

[0082] In the third step, similarly, determine the initial crack rate corresponding to the opposite direction of each preset direction.

[0083] It should be noted that the method for obtaining the initial crack rate corresponding to the reverse direction can be the same as the method for obtaining the initial crack rate corresponding to the preset direction, and will not be repeated here.

[0084] Step S3: determining a reference usable index corresponding to each preset direction according to a difference between an initial crack rate corresponding to each preset direction and an initial crack rate corresponding to the opposite direction thereof.

[0085] As an example, this step may include the following steps:

[0086] In the first step, any preset direction is determined as the marking direction, and the opposite direction of the marking direction is determined as the reference direction.

[0087] In the second step, the absolute value of the difference between the initial crack rate corresponding to the above-mentioned marked direction and the initial crack rate corresponding to the above-mentioned reference direction is determined as the possible change difference of the crack corresponding to the above-mentioned marked direction.

[0088] The third step is to determine the reference available indicators corresponding to the above-mentioned marked directions based on the possible changes in the cracks corresponding to the above-mentioned marked directions.

[0089] Among them, the possible change difference of cracks can be negatively correlated with the reference available indicators.

[0090] For example, the formula for determining the reference available indicator corresponding to the preset direction may be:

[0091] ;in, is the reference availability index corresponding to the i-th preset direction. i is the sequence number of the preset direction. is the normalization function. It is the absolute value function. is the initial crack rate corresponding to the i-th preset direction. is the initial crack rate corresponding to the opposite direction of the i-th preset direction. is the possible change difference of the crack corresponding to the i-th preset direction. It is a pre-set factor greater than 0, mainly used to prevent the denominator from being 0, such as, It can be 0.001.

[0092] It should be noted that cracks usually present a banded structure, wherein the banded structure is also called a linear structure. If the point to be detected is the crack source position point, there is often a crack line that runs through the point to be detected, and the gravity anomaly at the point farther away from the point to be detected on the crack line is relatively smaller, and the gravity anomaly at the point closer to the point to be detected is relatively larger. In other words, the gravity anomalies at the symmetrical positions on both sides of the point to be detected on the crack line are often similar. The crack source position point can be the source position point where the crack is generated, which is a type of crack position point. The neighboring points in the preset direction and the neighboring points in the opposite direction of the preset direction can often represent the symmetrical positions on both sides of the point to be detected. When The smaller it is, the more similar the gravity anomalies at the symmetrical positions on both sides of the detection point on the straight line in the i-th preset direction are. It often means that the straight line in the i-th preset direction is more consistent with the distribution law of the crack line, and the straight line in the i-th preset direction is more likely to be a crack line passing through the detection point, and the detection point is more likely to be the crack source location point.

[0093] Step S4: determining the target crack rate corresponding to the point to be detected based on the initial crack rates corresponding to all preset directions and the reference available indicators.

[0094] As an example, this step may include the following steps:

[0095] In the first step, the accumulated values ​​of the reference available indicators corresponding to all preset directions are determined as the available representative indicators.

[0096] In the second step, the ratio of the reference available index corresponding to each preset direction to the above available representative index is determined as the target available weight corresponding to each preset direction.

[0097] In the third step, the product of the target available weight corresponding to each preset direction and the initial crack rate is determined as the target sub-crack factor corresponding to each preset direction.

[0098] The fourth step is to determine the target crack rate corresponding to the point to be detected based on the target sub-crack factors corresponding to all preset directions.

[0099] Among them, the target sub-crack factor can be positively correlated with the target crack rate.

[0100] For example, the formula for determining the target crack rate corresponding to the point to be detected can be:

[0101] ; Where y is the target crack rate corresponding to the point to be detected. is a normalization function. N is the number of preset directions. i is the index of the preset direction. is the reference available index corresponding to the i-th preset direction. B is the available representative index, which is the cumulative value of the reference available indexes corresponding to all preset directions. is the initial crack rate corresponding to the i-th preset direction. is the target available weight corresponding to the i-th preset direction. is the target sub-crack factor corresponding to the i-th preset direction.

[0102] It should be noted that when The smaller it is, the more similar the gravity anomalies at the symmetrical positions on both sides of the detection point on the straight line where the i-th preset direction is located are. This often means that the straight line where the i-th preset direction is located is more likely to be a crack line passing through the detection point, and the detection point is more likely to be the crack source location point. This often means that the position points on the straight line where the i-th preset direction are located are more valuable for analysis, for example, the neighborhood points on the i-th preset direction are more valuable for analysis. Can be used as When The larger the value, the greater the relative gravity difference between the target point and the neighboring points in the i-th preset direction. This indicates that the gravity value at the target point is more likely to have a larger gravity anomaly than other surrounding gravity exploration points, and the greater the likelihood that the target point is a fracture location. Therefore, the larger the value of y, the more likely the gravity value at the target point is to have a larger gravity anomaly than other surrounding gravity exploration points, and the more likely the target point is a fracture source location.

[0103] Secondly, the greater the target crack rate corresponding to the test point, the more likely it is that the gravity value at the test point has a significant gravity anomaly compared to other surrounding gravity exploration points. Therefore, if the target crack rate corresponding to the test point is greater than the preset crack threshold, it often indicates that the gravity value at the test point has a significant gravity anomaly compared to other surrounding gravity exploration points, which often indicates that the test point can be preliminarily judged to be the location of the crack source. If the target crack rate corresponding to the test point is less than or equal to the preset crack threshold, it often indicates that the gravity value at the test point has no significant gravity anomaly compared to other surrounding gravity exploration points, which often indicates that the test point can be preliminarily judged to be a location where no cracks have occurred.

[0104] Step S5: If the target crack rate corresponding to the point to be detected is greater than the preset crack threshold, the target crack direction is determined according to the distribution of the initial crack rates corresponding to all preset directions.

[0105] The preset crack threshold may be a preset maximum target crack rate allowed when the gravity anomaly is considered to be absent and not a crack source. For example, the preset crack threshold may be 0.7.

[0106] It should be noted that for multiple locations at the same distance from the crack source, the closer they are to the crack line, the greater the gravity anomaly. In other words, if the crack line is the line from the 0° direction to the 90° direction, the gravity anomaly tends to increase from the 0° direction to the direction corresponding to the crack line. In this case, a linear fit can be performed on the gravity anomaly from the 0° direction to the direction corresponding to the crack line to obtain a fitted line. The gravity anomaly tends to decrease from the direction corresponding to the crack line to the direction 90° behind the crack line, and another fitted line can be obtained. The gravity anomaly tends to increase from the direction 90° behind the crack line to the direction 180° behind the crack line, and another fitted line can be obtained from the gravity anomaly from the direction 90° behind the crack line to the direction 180° behind the crack line. In this case, a linear fit can be performed on the gravity anomaly from the direction 90° behind the crack line to the direction 180° behind the crack line to obtain another fitted line, resulting in a total of three fitted lines.

[0107] Similarly, if the crack line is a straight line from the straight line in the 90° direction to the straight line in the 180° direction, then the gravity anomaly in the 90° from the 0° direction to the direction corresponding to the crack line often decreases. At this time, the gravity anomaly in the 90° from the 0° direction to the direction corresponding to the crack line can often be linearly fitted to obtain a fitting straight line; the gravity anomaly in the 90° from the direction corresponding to the crack line to the direction corresponding to the crack line often increases. At this time, the gravity anomaly in the 90° from the direction corresponding to the crack line to the direction corresponding to the crack line can often be linearly fitted to obtain another fitting straight line; the gravity anomaly in the direction corresponding to the crack line to the 180° direction often decreases. At this time, the gravity anomaly in the direction corresponding to the crack line to the 180° direction can often be linearly fitted to obtain another fitting straight line, and a total of 3 fitting straight lines can be obtained.

[0108] As an example, this step may include the following steps:

[0109] The first step is to construct a crack rate variation function and its error sum function based on the distribution of the initial crack rates corresponding to all preset directions, which can include the following sub-steps:

[0110] The first sub-step, The value range of is [0°, 180°]. When it belongs to [0°, 90°], the crack rate variation function can be:

[0111] ;

[0112] ;

[0113] ;in, is the crack rate change function value when the independent variable is a. It is equal to the fitted value of the initial crack rate corresponding to the preset direction with the corresponding angle a. a is the angle corresponding to different preset directions. k is the slope of the fitted line included in the crack rate change function. 、 and are the intercepts of different fitting straight lines contained in the crack rate variation function.

[0114] In the second sub-step, when When it belongs to (90°, 180°], the crack rate variation function can be:

[0115] ;

[0116] ;

[0117] .

[0118] In the third sub-step, the error sum function of the crack rate variation function can be:

[0119] Where Q is the error and function value of the crack rate variation function. N is the number of preset directions. i is the sequence number of the preset direction. It is the reference available indicator corresponding to the i-th preset direction. It is the crack rate variation function value when the independent variable takes the angle corresponding to the i-th preset direction. is the initial crack rate corresponding to the i-th preset direction.

[0120] In the second step, the crack rate variation function is solved by the least squares method, and the solution with the minimum error and function is taken as the final solution to obtain the target crack direction.

[0121] Among them, the error and function are solved when they are minimum. The corresponding direction is the target crack direction.

[0122] It should be noted that the unknowns in the crack rate variation function can be solved by the least squares method. ,k, 、 and .

[0123] Step S6: determining the final crack rate corresponding to the point to be detected according to the gravity change in the target crack direction, and judging whether the point to be detected is a crack source location point based on the final crack rate.

[0124] As an example, this step may include the following steps:

[0125] The first step is to obtain gravity data corresponding to a preset number of preset position points in the target crack direction of the above-mentioned detection point and the neighborhood points of these preset position points in each preset direction.

[0126] The preset number may be a pre-set number, which may be 10. The preset position points may be position points every 10 meters in the target crack direction. For example, the first preset position point in the target crack direction of the point to be detected may be a neighboring point in the target crack direction of the point to be detected. The neighboring point of the preset position point in the preset direction, that is, the neighboring point in the preset direction of the preset position point, may be a position point 10 meters away from the preset position point in the preset direction.

[0127] It should be noted that the method for acquiring the gravity data corresponding to the preset position point and its neighborhood points may be the same as the method for acquiring the gravity data corresponding to the point to be detected, and will not be described in detail here.

[0128] In the second step, the absolute value of the difference between the gravity data corresponding to each preset position point and the gravity data corresponding to its neighboring points in each preset direction is determined as the temporary gravity difference of each preset position point in each preset direction.

[0129] In the third step, the accumulated value of the temporary gravity difference of each preset position point in all preset directions is determined as the reference gravity difference corresponding to each preset position point.

[0130] The fourth step, based on the reference gravity differences corresponding to a preset number of preset position points in the target crack direction, determines the final crack rate corresponding to the to-be-detected point, which may include the following sub-steps:

[0131] In the first sub-step, the distance between the point to be detected and each preset position point is determined as the target distance corresponding to each preset position point.

[0132] In the second sub-step, a final crack rate corresponding to the point to be detected is determined based on the target crack rate corresponding to the point to be detected, and the reference gravity difference and target distance corresponding to a preset number of preset position points in the target crack direction.

[0133] The target crack rate and the reference gravity difference may both be positively correlated with the final crack rate, and the target distance may be negatively correlated with the final crack rate.

[0134] For example, the formula for determining the final crack rate corresponding to the point to be detected can be:

[0135] ; Where F is the final crack rate corresponding to the point to be detected. is a normalized function. y is the target crack rate corresponding to the point to be detected. n is the preset number. j is the serial number of the preset position point in the target crack direction. is the reference gravity difference corresponding to the j-th preset position point in the target crack direction. is the target distance corresponding to the jth preset position point in the target crack direction.

[0136] It should be noted that when When it is larger, it often indicates that the gravity value at the j-th preset position point in the target fracture direction is more likely to have a larger gravity anomaly compared with other surrounding gravity exploration points, which often indicates that the j-th preset position point in the target fracture direction is more likely to be a fracture location point. Can be used as The closer the jth preset location is to the target point, the greater its reference value. If the target point is a fracture source, the gravity anomaly corresponding to the jth preset location will be greater. A larger value for y indicates a higher probability of a larger gravity anomaly at the target point compared to other surrounding gravity exploration points, indicating a higher probability of the target point being a fracture source. Therefore, a larger value for F indicates a higher probability of the target point being a fracture source.

[0137] In the fifth step, if the final crack rate corresponding to the above-mentioned point to be detected is greater than the preset abnormal threshold, the above-mentioned point to be detected is determined to be the crack source location point.

[0138] The preset anomaly threshold may be a preset maximum final crack rate allowed when a gravity anomaly is considered to exist but not a crack source. For example, the preset anomaly threshold may be 0.6.

[0139] refer to Figure 3 Based on the same inventive concept as the above-mentioned method embodiment, the present invention provides a geological fracture exploration and detection system for geological engineering. The system includes a memory and a processor. The processor is used to process instructions stored in the memory to implement the steps of the geological fracture exploration and detection method for geological engineering, which may specifically include:

[0140] Gravity data acquisition module 301, used to acquire gravity data corresponding to the point to be detected and gravity data corresponding to neighboring points in each preset direction;

[0141] An initial crack rate determination module 302 is configured to determine an initial crack rate corresponding to each preset direction based on the difference between the gravity data corresponding to the point to be detected and the gravity data corresponding to its neighboring points in each preset direction. Similarly, an initial crack rate corresponding to the opposite direction of each preset direction is determined.

[0142] A reference usable index determining module 303 is configured to determine a reference usable index corresponding to each preset direction based on a difference between an initial crack rate corresponding to each preset direction and an initial crack rate corresponding to an opposite direction;

[0143] The target crack rate determination module 304 is used to determine the target crack rate corresponding to the point to be detected based on the initial crack rates corresponding to all preset directions and the reference available indicators;

[0144] The target crack direction determination module 305 is configured to determine the target crack direction based on the distribution of initial crack rates corresponding to all preset directions if the target crack rate corresponding to the point to be detected is greater than a preset crack threshold;

[0145] The determination and judgment module 306 is used to determine the final crack rate corresponding to the point to be detected according to the gravity change in the target crack direction, and to judge whether the point to be detected is a crack source location point based on the final crack rate.

[0146] Figure 4 FIG. 1 is a schematic diagram of the structure of a computer device provided by an embodiment of the present invention. For example, Figure 4 As shown, the computer device 400 includes: a memory 401, a processor 402, and a computer program 403 stored in the memory 401 and running on the processor 402, wherein when the processor 402 executes the computer program 403, the computer device can execute any one of the geological fracture exploration and detection methods for geological engineering introduced above.

[0147] Based on the same inventive concept as the above-described method embodiments, the present invention provides a server comprising a memory and a processor. The memory is configured to store executable program code, and the processor is configured to retrieve and execute the executable program code from the memory, thereby enabling the server to perform any of the above-described geological fracture exploration and detection methods for geological engineering.

[0148] Based on the same inventive concept as the above-mentioned method embodiment, the present invention provides a computer program product, which includes: computer program code, which, when running on a computer, enables the computer to execute any one of the above-mentioned geological fracture exploration and detection methods for geological engineering.

[0149] Based on the same inventive concept as the above-mentioned method embodiment, the present invention provides a computer-readable storage medium, which stores computer program code. When the computer program code is run on a computer, the computer executes any one of the above-mentioned geological fracture exploration and detection methods for geological engineering.

[0150] In summary, when conducting geological fracture exploration and detection, the present invention not only considers the difference between the gravity data corresponding to the point to be detected and the gravity data corresponding to its neighboring points in each preset direction, but also comprehensively considers multiple characteristics related to the change law of geological fractures, such as reference available indicators, target fracture rate and gravity change in the target fracture direction, thereby quantifying the final fracture rate corresponding to the point to be detected. Based on the final fracture rate, it is possible to relatively objectively judge whether the point to be detected is the fracture source location point, which can improve the accuracy of geological fracture exploration and detection to a certain extent.

[0151] The above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. These modifications or replacements do not deviate the essence of the corresponding technical solutions from the scope of the technical solutions of the embodiments of the present invention, and should all be included in the scope of protection of the present invention.

Claims

1. A geological fracture exploration and detection method for geological engineering, characterized in that: The following steps are involved: Obtain the gravity data corresponding to the point to be detected and the gravity data corresponding to each of its neighboring points in a preset direction; Based on the difference between the gravity data corresponding to the point to be detected and the gravity data corresponding to its neighboring points in each preset direction, the initial crack rate corresponding to each preset direction is determined. Similarly, the initial crack rate corresponding to the opposite direction of each preset direction is determined. Determining a reference usable index corresponding to each preset direction according to a difference between an initial crack rate corresponding to each preset direction and an initial crack rate corresponding to the opposite direction; Determine the target crack rate corresponding to the point to be tested based on the initial crack rates corresponding to all preset directions and the reference available indicators; If the target crack rate corresponding to the point to be detected is greater than the preset crack threshold, the target crack direction is determined based on the distribution of the initial crack rates corresponding to all preset directions; According to the gravity change in the target crack direction, the final crack rate corresponding to the detection point is determined, and based on the final crack rate, whether the detection point is the crack source location point is determined; The formula for determining the initial crack rate corresponding to the preset direction is: ;in, is the initial crack rate corresponding to the i-th preset direction; i is the serial number of the preset direction; It is the absolute value function; is the gravity data corresponding to the neighborhood point in the i-th preset direction; G is the gravity data corresponding to the point to be detected; The formula for determining the reference available indicator corresponding to the preset direction is: ;in, is the reference available indicator corresponding to the i-th preset direction; is the normalization function; is the initial crack rate corresponding to the i-th preset direction; is the initial crack rate corresponding to the opposite direction of the i-th preset direction; is a pre-set factor greater than 0; The method of determining the final crack rate corresponding to the point to be detected according to the gravity change in the target crack direction includes: Obtaining gravity data corresponding to a preset number of preset position points in the target crack direction of the point to be detected and the neighboring points of these preset position points in each preset direction; Determine the absolute value of the difference between the gravity data corresponding to each preset position point and the gravity data corresponding to its neighboring points in each preset direction as the temporary gravity difference of each preset position point in each preset direction; Determine the cumulative value of the temporary gravity difference of each preset position point in all preset directions as the reference gravity difference corresponding to each preset position point; The final crack rate corresponding to the point to be detected is determined according to the reference gravity differences corresponding to a preset number of preset position points in the target crack direction.

2. A geological fracture exploration and detection method for geological engineering according to claim 1, characterized in that: The determining of the initial crack rate corresponding to each preset direction based on the difference between the gravity data corresponding to the to-be-detected point and the gravity data corresponding to the neighboring points in each preset direction includes: Determine any preset direction as a marked direction, and determine the absolute value of the difference between the gravity data corresponding to the point to be detected and the gravity data corresponding to the neighboring points in the marked direction as the target gravity difference corresponding to the marked direction; An initial crack rate corresponding to the marked direction is determined according to a target gravity difference corresponding to the marked direction, wherein the target gravity difference is positively correlated with the initial crack rate.

3. A geological fracture exploration and detection method for geological engineering according to claim 1, characterized in that: Determining the reference usable index corresponding to each preset direction according to the difference between the initial crack rate corresponding to each preset direction and the initial crack rate corresponding to the opposite direction thereof includes: Determine any preset direction as a marking direction, and determine the opposite direction of the marking direction as a reference direction; Determine the absolute value of the difference between the initial crack rate corresponding to the marked direction and the initial crack rate corresponding to the reference direction as the possible change difference of the crack corresponding to the marked direction; According to the possible crack change difference corresponding to the marked direction, a reference usable index corresponding to the marked direction is determined, wherein the possible crack change difference is negatively correlated with the reference usable index.

4. A geological fracture exploration and detection method for geological engineering according to claim 1, characterized in that: The step of determining the target crack rate corresponding to the point to be detected based on the initial crack rates corresponding to all preset directions and the reference available indicators includes: Determine the accumulated values ​​of the reference available indicators corresponding to all preset directions as the available representative indicator; Determine the target available weight corresponding to each preset direction by the ratio of the reference available index corresponding to each preset direction to the available representative index; The product of the target available weight corresponding to each preset direction and the initial crack rate is determined as the target sub-crack factor corresponding to each preset direction; The target crack rate corresponding to the point to be detected is determined according to the target sub-crack factors corresponding to all preset directions, wherein the target sub-crack factor is positively correlated with the target crack rate.

5. A geological fracture exploration and detection method for geological engineering according to claim 1, characterized in that: The step of determining the target crack direction based on the distribution of the initial crack rates corresponding to all preset directions includes: According to the distribution of initial crack rates corresponding to all preset directions, a crack rate variation function and its error sum function are constructed; The crack rate variation function is solved by the least squares method, and the solution with the minimum error and function is taken as the final solution to obtain the target crack direction.

6. A geological fracture exploration and detection method for geological engineering according to claim 5, characterized in that: The method of constructing a crack rate variation function and its error sum function based on the distribution of the initial crack rates corresponding to all preset directions includes: The value range of is [0°, 180°]. When it belongs to [0°, 90°], the crack rate variation function is: ； ； ;in, is the crack rate change function value when the independent variable is a, which is equal to the fitted value of the initial crack rate corresponding to the preset direction with the corresponding angle a; a is the angle corresponding to different preset directions; k is the slope of the fitted line contained in the crack rate change function; 、 and is the intercept of different fitting straight lines included in the crack rate variation function; when When it belongs to (90°, 180°], the crack rate variation function is: ； ； ； The error sum function of the crack rate variation function is: ; Wherein, Q is the error and function value of the crack rate variation function; N is the number of preset directions; i is the serial number of the preset direction; is the reference available indicator corresponding to the i-th preset direction; is the crack rate variation function value when the independent variable takes the angle corresponding to the i-th preset direction; is the initial crack rate corresponding to the i-th preset direction.

7. A geological fracture exploration and detection method for geological engineering according to claim 1, characterized in that: The determining of the final crack rate corresponding to the to-be-detected point according to the reference gravity differences corresponding to a preset number of preset position points in the target crack direction includes: Determine the distance between the point to be detected and each preset position point as the target distance corresponding to each preset position point; The final crack rate corresponding to the point to be detected is determined based on the target crack rate corresponding to the point to be detected, and the reference gravity difference and target distance corresponding to a preset number of preset position points in the target crack direction, wherein the target crack rate and the reference gravity difference are both positively correlated with the final crack rate, and the target distance is negatively correlated with the final crack rate.

8. A geological fracture exploration and detection method for geological engineering according to claim 1, characterized in that: The step of judging whether the point to be detected is a crack source location point based on the final crack rate includes: If the final crack rate corresponding to the to-be-detected point is greater than a preset abnormal threshold, the to-be-detected point is determined to be a crack source location point.

9. A geological fracture exploration and detection system for geological engineering, characterized in that: The method comprises a processor and a memory, wherein the processor is used to process instructions stored in the memory to implement a geological fracture exploration and detection method for geological engineering according to any one of claims 1 to 8.

Images