Method for detecting direction of deep hole based on slurry diffusion theory

Abstract

The invention discloses a method for detecting the direction of a deep hole based on a slurry diffusion theory, and the method comprises the steps: arranging detection holes based on the geological data of a detection place; the mix proportion of grout is determined, grouting is conducted on the bottom of the detection hole, and grouting parameters including the grout injection rate, the grouting time, the grouting pressure and the grout density are recorded; on the basis of a theoretical formula of slurry diffusion, according to the recorded grouting parameter information, calculating a slurry injection rate time history curve, and according to the two curves, calculating the distance between the probe hole and the cavity; based on the multiple sets of grouting parameters, the distance from each detection hole to the cavity and the radius of the cavity are calculated; and determining the space position of the hole by combining the distance from each detection hole to the hole and the radius of the hole based on a space positioning principle. According to the method, the position and the size of the cavity can be estimated in the three-dimensional space model by comparing and analyzing the grout injection rate time history curve and the monitoring curve in the grouting process in combination with the specific geological condition of the site.

Description

Technical Field

[0001] The present invention relates to the field of underground detection, in particular to the field of cavity position detection, and specifically to a method for detecting the deep cavity position based on slurry diffusion theory. Background Art

[0002] In engineering, underground cavities are a common geological defect that can lead to structural instability, foundation leakage, construction difficulties and other problems. Therefore, determining the location of cavities is of great significance to engineering safety.

[0003] At present, the commonly used methods for cavity detection mainly include ultrasonic detection, remote sensing cavitation, electrical exploration, etc. These methods have their own advantages, but they also have obvious shortcomings. For example, ultrasonic tomography has high requirements for operating technology and environmental conditions, remote sensing technology is easily affected by surface cover, and electrical difference cavitation will be disturbed by soil moisture. In addition, when the cavity volume is small, especially when the cavity is located deep, the detection accuracy of the above methods will be significantly reduced, and it may even be impossible to identify. Summary of the invention

[0004] In order to overcome the above-mentioned deficiencies of the prior art, the present invention provides a method for detecting the position of deep cavities based on slurry diffusion theory. By utilizing drilling grouting technology, combining slurry diffusion theory and spatial positioning technology, a new method for detecting deep cavities is provided.

[0005] According to one aspect of the present invention, a method for detecting the position of deep cavities based on slurry diffusion theory is provided, comprising:

[0006] Arrange the detection holes based on the geological data of the detection site;

[0007] Determine the slurry mix ratio, grout the bottom of the detection hole, and record the grouting parameters including slurry injection rate, grouting time, grouting pressure and slurry density;

[0008] Based on the theoretical formula of slurry diffusion and the recorded grouting parameter information, the slurry injection rate time curve is calculated to determine whether the monitoring curve and the time curve coincide with each other. If they coincide with each other, the grouting parameters are deleted, the detection hole position is re-determined, and the above steps are repeated. If they do not coincide with each other, the grouting parameters are retained, and the positions of other detection holes are determined based on the position of the first detection hole. The above steps are repeated to obtain multiple sets of grouting parameters.

[0009] Based on several groups of grouting parameters, the distance from each detection hole to the cavity and the cavity radius are calculated respectively;

[0010] Based on the principle of spatial positioning, the spatial position of the cavity is determined by combining the distance from each detection hole to the cavity and the cavity radius.

[0011] As a further technical solution, based on the geological data of the detection site, the detection holes are arranged, including: obtaining the geological data of the detection site, evaluating the possibility, development degree and potential orientation of the local cavity, and digging a detection hole in the vicinity of the cavity.

[0012] As a further technical solution, the mix ratio of the slurry is determined, including: determining the mix ratio of the slurry based on the development degree of the detected ground fissures, when the opening of the nearby fissures is large, using thick slurry, and when the opening of the nearby fissures is small, using thin slurry.

[0013] As a further technical solution, the theoretical formula for the slurry diffusion is as follows:

[0014] (1);

[0015] In the formula, is the grouting time; is the number of monitoring points; is the slurry diffusion radius; is the monitoring value of slurry injection rate; For the The monitoring value of the slurry injection rate at the moment; V is the total volume of the slurry diffused in the space within a certain period of time; t is the moment of the grouting monitoring value; The grouting monitoring value a moment of time; The grouting monitoring value - 1 moment of time;

[0016] (2);

[0017] According to Darcy's law, we can get:

[0018] (3);

[0019] In the formula, is the permeability coefficient of slurry in the formation; is the slurry hydraulic gradient; is the spherical diffusion area of ​​the slurry;

[0020] Integrating formula (3), we can get:

[0021] (4);

[0022] In the formula, It is the pressure head of the slurry in the grouting hole; is the pressure head of groundwater; For the Calculated value of injection rate of moment-to-moment time history curve; is the radius of the grouting hole;

[0023] Substituting formula (2) into formula (4), we can get:

[0024] (5).

[0025] As a further technical solution, judging whether the monitoring curve and the time course curve coincide with each other includes: comparing the slurry injection rate time course curve with the detection curve, when the two are basically consistent with each other over time, and then there is a significant difference (such as Figure 2 (a) shows that the slurry flows through the cavity. When the two are always basically consistent (such as Figure 2 (b) indicates that the slurry has not flowed through the cavity and the detection hole position needs to be re-determined.

[0026] As a further technical solution, when the slurry flows through the cavity, it includes: Figure 3 The following stages are shown:

[0027] Phase I ): The slurry has not yet flowed into the cavity, and the time course curve is consistent with the monitoring curve;

[0028] Phase II : The slurry flows into the cavity, and the time-history curve and the monitoring curve show obvious differences, and the former gradually becomes smaller than the latter. The difference between the two continues to expand with time. In formulas (1) and (2), let , calculate the distance from the detection hole to the hole ;

[0029] Phase 3 : The cavities are filled with slurry. Since the slurry continues to flow along the cracks at this time, the theoretical formula of slurry diffusion is applicable again, so the calculation curve and the monitoring curve are distributed in parallel.

[0030] As a further technical solution, the cavity radius The calculation method is as follows:

[0031] (4);

[0032] (5);

[0033] In the formula, is the volume of the cavity, and the cavity is assumed to be spherical; is the calculated value of slurry injection rate; is the monitoring value of slurry injection rate; is the cavity radius; is the injection rate of the time history curve a moment of time; is the injection rate of the time history curve - 1 moment of time;

[0034] Calculate the cavity radius for each detection hole , and then take the average ;

[0035] when When ;

[0036] They are respectively the hole radius corresponding to the first detection hole, the hole radius corresponding to the second detection hole, the hole radius corresponding to the third detection hole, and the hole radius corresponding to the fourth detection hole.

[0037] As a further technical solution, the method for calculating the hole position is as follows:

[0038] Based on the principle of spatial positioning, the coordinates of the center point of the cavity obey the following set of equations:

[0039] (6);

[0040] In the formula, The spatial coordinates of the bottom of the first detection hole; is the spatial coordinate of the bottom of the second detection hole; The spatial coordinates of the bottom of the third detection hole; is the spatial coordinate of the bottom of the fourth detection hole; is the distance between the bottom of the first detection hole and the cavity surface; is the distance between the bottom of the second detection hole and the cavity surface; is the distance between the bottom of the third detection hole and the cavity surface; is the distance between the bottom of the fourth detection hole and the cavity surface; is the coordinate position of the center point of the hole; is the average value of the cavity radius of multiple detection holes;

[0041] The objective function of the hole center coordinates is as follows:

[0042] (7);

[0043] By finding the gradient of formula (9), we can get:

[0044] (8);

[0045] ;

[0046] The gradient optimization method is used to solve the coordinates of the center point of the hole, and we can get:

[0047] (9);

[0048] In the formula, is the number of iterations; is the step length.

[0049] According to one aspect of the present invention specification, there is provided an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method for detecting the position of deep cavities based on slurry diffusion theory as described in any one of claims 1 to 7 is implemented.

[0050] According to one aspect of the present specification, there is provided a non-transitory computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed by a processor, the method for detecting the position of deep cavities based on slurry diffusion theory as described in any one of claims 1 to 7 is implemented.

[0051] Compared with the prior art, the beneficial effect of the present invention is that the present invention can estimate the location and size of the cavity in the three-dimensional space model by comparing and analyzing the slurry injection rate time curve and the monitoring curve during the grouting process in combination with the specific geological conditions on site, thereby providing convenience for engineers to make on-site decisions. Since drilling grouting is a part of most engineering construction, this method can detect cavities without additional detection equipment. At the same time, the present invention can identify small-volume cavities that are difficult to identify with traditional detection technology, providing a new idea for deep cavity detection methods. BRIEF DESCRIPTION OF THE DRAWINGS

[0052] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, a brief introduction is given below to the drawings used in the embodiments or the description of the prior art. Obviously, the drawings described below are 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.

[0053] Figure 1 A schematic flow chart of a method for detecting the position of deep cavities based on slurry diffusion theory provided in an embodiment of the present invention.

[0054] Figure 2 A slurry injection rate time history curve diagram of a method for detecting the position of deep cavities based on slurry diffusion theory provided in an embodiment of the present invention.

[0055] Figure 3 A diagram of the stage division of a slurry injection rate time history curve of a method for detecting the orientation of deep cavities based on slurry diffusion theory provided in an embodiment of the present invention.

[0056] Figure 4 A diagram showing the phase division of the injection rate time history curve of detection hole No. 1 of a method for detecting the orientation of deep cavities based on slurry diffusion theory provided in an embodiment of the present invention.

[0057] Figure 5 A diagram showing the phase division of the injection rate time history curve of detection hole No. 2 of a method for detecting the orientation of deep cavities based on slurry diffusion theory provided in an embodiment of the present invention.

[0058] Figure 6 A diagram showing the phase division of the injection rate time history curve of detection hole No. 3 of a method for detecting the orientation of deep cavities based on slurry diffusion theory provided in an embodiment of the present invention.

[0059] Figure 7 A diagram showing the phase division of the injection rate time history curve of detection hole No. 4 of a method for detecting the orientation of deep cavities based on slurry diffusion theory provided in an embodiment of the present invention.

[0060] Figure 8 A distribution diagram of detection holes and cavity positions of a method for detecting deep cavity positions based on slurry diffusion theory provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0061] The terms "including" and "having" and any variations thereof in the specification and claims of the present invention and the above-mentioned drawings are intended to cover non-exclusive inclusions, for example, a process, method, system, product or apparatus comprising a series of steps or units is not necessarily limited to the steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to these processes, methods, products or apparatuses.

[0062] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention. In addition, the technical features in the various embodiments or single embodiments provided by the present invention are arbitrarily combined with each other to form a new technical solution. This combination is not restricted by the sequence of steps and / or the structural composition mode, but must be based on the ability of ordinary technicians in this field to achieve. When the combination of technical solutions is contradictory or cannot be achieved, it should be considered that this combination of technical solutions does not exist and is not within the scope of protection required by the present invention.

[0063] The embodiment of the present invention provides a method for detecting the position of deep cavities based on slurry diffusion theory, comprising:

[0064] Step 1: Arrange the detection holes based on the geological data of the detection site;

[0065] Step 2, determine the mix ratio of the slurry, grout the bottom of the detection hole, and record the grouting parameters including slurry injection rate, grouting time, grouting pressure and slurry density;

[0066] Step 3, based on the theoretical formula of slurry diffusion and the recorded grouting parameter information, calculate the slurry injection rate time curve, determine whether the monitoring curve and the time curve overlap, if they overlap, delete the grouting parameters, re-determine the detection hole position, and repeat the above steps, if they do not overlap, retain the grouting parameters, determine the positions of other detection holes based on the position of the first detection hole, repeat the above steps, and obtain multiple sets of grouting parameters;

[0067] Step 4, based on several sets of grouting parameters, respectively calculating the distance from each detection hole to the cavity and the cavity radius;

[0068] Step 5: Based on the spatial positioning principle, the spatial position of the cavity is determined in combination with the distance from each detection hole to the cavity and the cavity radius.

[0069] The method for arranging the detection holes in step 1 is specifically as follows: firstly, the geological data of a certain area is obtained to preliminarily estimate the possibility, development degree and potential orientation of the local voids; and a detection hole is dug in the vicinity of the possible voids.

[0070] Determining the mix ratio of the slurry in step 2 includes: determining the mix ratio of the slurry based on the development degree of the detected ground fissures, selecting a thick slurry when the opening of the nearby fissures is large, and selecting a thin slurry when the opening of the nearby fissures is small.

[0071] The derivation process of the formula for the slurry injection rate time history curve in step 3 is as follows:

[0072] The total volume of slurry diffused in space within a certain period of time is:

[0073] (1);

[0074] In the formula, is the grouting time; is the number of monitoring points; is the slurry diffusion radius; is the monitoring value of slurry injection rate; For the The monitoring value of the slurry injection rate at the moment; V is the total volume of the slurry diffused in the space within a certain period of time; t is the moment of the grouting monitoring value; The grouting monitoring value a moment of time; The grouting monitoring value - 1 moment of time;

[0075] (2);

[0076] According to Darcy's law, we can get:

[0077] (3);

[0078] In the formula, is the permeability coefficient of slurry in the formation; is the slurry hydraulic gradient; is the spherical diffusion area of ​​the slurry;

[0079] Integrating formula (3), we can get:

[0080] (4);

[0081] In the formula, It is the pressure head of the slurry in the grouting hole; is the pressure head of groundwater; For the Calculated value of injection rate of moment-to-moment time history curve; is the radius of the grouting hole;

[0082] Substituting formula (2) into formula (4), we can get:

[0083] (5);

[0084] In the formula, For the Calculated value of injection rate of moment-to-moment time history curve; is the permeability coefficient of slurry in the formation; It is the pressure head of the slurry in the grouting hole; is the pressure head of groundwater; For the Monitoring value of slurry injection rate at every moment; is the radius of the grouting hole; is the injection rate of the time history curve a moment of time; is the injection rate of the time history curve -1 moment of time.

[0085] In step 3, it is determined whether the monitoring curve and the time course curve coincide with each other, including: comparing the slurry injection rate time course curve with the detection curve, when the two are basically consistent with each other over time, and then there is a significant difference (such as Figure 2 (a) shows that the slurry flows through the cavity. When the two are always basically consistent (such as Figure 2 (b) indicates that the slurry has not flowed through the cavity and the detection hole position needs to be re-determined.

[0086] When the slurry flows through the cavity, it includes Figure 3 The following stages are shown:

[0087] Phase I ): The slurry has not yet flowed into the cavity, and the time course curve is basically consistent with the monitoring curve;

[0088] Phase II : The slurry flows into the cavity. Since the slurry does not flow completely in the cracks at this time, the theoretical formula of slurry diffusion is no longer applicable. Therefore, there will be obvious differences between the time-history curve and the monitoring curve, and the former is gradually smaller than the latter. The difference between the two continues to expand over time. In formulas (1) and (2), let , calculate the distance from the detection hole to the hole ;

[0089] Phase 3 : The cavity is filled with slurry. Since the slurry continues to flow along the cracks at this time, the theoretical formula of slurry diffusion is applicable again, so the time-history curve and the monitoring curve are distributed in parallel.

[0090] in, It is the dividing point between the two curves from the same injection rate to the different injection rates. It is the dividing point between the two curves from injection rate change to injection rate stability.

[0091] In step 4, the cavity radius The calculation method is as follows:

[0092] (6);

[0093] (7);

[0094] In the formula, is the volume of the cavity, and the cavity is assumed to be spherical; is the calculated value of slurry injection rate; is the monitoring value of slurry injection rate; is the cavity radius; is the injection rate of the time history curve a moment of time; is the injection rate of the time history curve - 1 moment of time;

[0095] Calculate the cavity radius for each detection hole , and then take the average ;

[0096] when When ;

[0097] They are respectively the hole radius corresponding to the first detection hole, the hole radius corresponding to the second detection hole, the hole radius corresponding to the third detection hole, and the hole radius corresponding to the fourth detection hole.

[0098] The calculation method for the hole position in step 6 is as follows:

[0099] Based on the principle of spatial positioning, the coordinates of the center point of the cavity obey the following set of equations:

[0100] (8);

[0101] In the formula, The spatial coordinates of the bottom of the first detection hole; is the spatial coordinate of the bottom of the second detection hole; The spatial coordinates of the bottom of the third detection hole; is the spatial coordinate of the bottom of the fourth detection hole; is the distance between the bottom of the first detection hole and the cavity surface; is the distance between the bottom of the second detection hole and the cavity surface; is the distance between the bottom of the third detection hole and the cavity surface; is the distance between the bottom of the fourth detection hole and the cavity surface; is the coordinate position of the center point of the hole; is the average value of the cavity radius of multiple detection holes;

[0102] The objective function of the hole center coordinates is as follows:

[0103] (9);

[0104] By finding the gradient of formula (9), we can get:

[0105] (10);

[0106] ;

[0107] The gradient optimization method is used to solve the coordinates of the center point of the hole, and we can get:

[0108] (11);

[0109] In the formula, is the number of iterations; is the step length.

[0110] Through the above detection method, the distance between the first detection hole and the hole and the hole radius are calculated. Based on the first detection hole, the positions of the other three detection holes are determined, and then steps S1-S3 are repeated, and the distance between the detection hole and the hole and the hole radius are calculated respectively. Finally, based on the spatial positioning principle, the spatial position of the hole is determined.

[0111] Experimental example:

[0112] According to the geological and hydrological data of a certain area, the orientation of the hole is roughly determined. In the vicinity of the hole, the due east is the x-axis, the due north is the y-axis, and the vertical direction is the z-axis. The coordinates of the center of the ground of the No. 1 detection hole are set as the origin. The depth of the No. 1 detection hole is 20m, so the coordinates of its bottom hole are ,like Figure 8 As shown. Grouting was performed in the No. 1 detection hole, and the grouting rate, grouting time, grouting pressure, and grouting density were recorded. The grouting rate time course curve was calculated by equations (1) and (2) as shown in Figure 4 As shown in the figure, after 35 minutes, there is a significant difference between the calculated injection rate and the monitored value, and the gap continues to expand, which indicates that the No. 1 detection hole is connected to the cavity. The distance between the bottom of the No. 1 detection hole and the cavity surface is calculated by equations (4) and (5): .

[0113] In the vicinity of detection hole No. 1, the spatial coordinates of the bottom of detection holes No. 2, No. 3, and No. 4 are determined in the same way: , , ,like Figure 5 Repeat the above process and calculate the slurry injection rate time history curve through equations (1) and (2) as shown in Figure 5 , Figure 6 and Figure 7 The calculated distances between the bottom of detection holes 2, 3, and 4 and the cavity surface are: , and .

[0114] The cavity radius is calculated by equations (6) and (7): , the corresponding average value is: The hole orientation is calculated by equations (9), (10) and (11): ,like Figure 8 shown.

[0115] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the method for detecting the position of deep cavities based on the slurry diffusion theory is implemented, as shown below:

[0116] Arrange the detection holes based on the geological data of the detection site;

[0117] Determine the slurry mix ratio, grout the bottom of the detection hole, and record the grouting parameters including slurry injection rate, grouting time, grouting pressure and slurry density;

[0118] Based on the theoretical formula of slurry diffusion and the recorded grouting parameter information, the slurry injection rate time curve is calculated to determine whether the monitoring curve and the time curve coincide with each other. If they coincide with each other, the grouting parameters are deleted, the detection hole position is re-determined, and the above steps are repeated. If they do not coincide with each other, the grouting parameters are retained, and the positions of other detection holes are determined based on the position of the first detection hole. The above steps are repeated to obtain multiple sets of grouting parameters.

[0119] Based on several groups of grouting parameters, the distance from each detection hole to the cavity and the cavity radius are calculated respectively;

[0120] Based on the principle of spatial positioning, the spatial position of the cavity is determined by combining the distance from each detection hole to the cavity and the cavity radius.

[0121] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the method for detecting the position of deep cavities based on the slurry diffusion theory is implemented as follows:

[0122] Arrange the detection holes based on the geological data of the detection site;

[0123] Determine the slurry mix ratio, grout the bottom of the detection hole, and record the grouting parameters including slurry injection rate, grouting time, grouting pressure and slurry density;

[0124] Based on the theoretical formula of slurry diffusion and the recorded grouting parameter information, the slurry injection rate time curve is calculated to determine whether the monitoring curve and the time curve coincide with each other. If they coincide with each other, the grouting parameters are deleted, the detection hole position is re-determined, and the above steps are repeated. If they do not coincide with each other, the grouting parameters are retained, and the positions of other detection holes are determined based on the position of the first detection hole. The above steps are repeated to obtain multiple sets of grouting parameters.

[0125] Based on several groups of grouting parameters, the distance from each detection hole to the cavity and the cavity radius are calculated respectively;

[0126] Based on the principle of spatial positioning, the spatial position of the cavity is determined by combining the distance from each detection hole to the cavity and the cavity radius.

[0127] In summary, the present invention can estimate the location and size of the cavity in a three-dimensional space model by comparing and analyzing the slurry injection rate time curve and the monitoring curve during the grouting process in combination with the specific geological conditions on site, thereby providing convenience for engineers to make on-site decisions. Since drilling grouting is a part of most engineering construction, this method can detect cavities without additional detection equipment. At the same time, the present invention can identify small-volume cavities that are difficult to identify with traditional detection technology, providing a new idea for deep cavity detection methods.

[0128] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. 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 replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the technical solutions of the embodiments of the present invention.

Claims

1. A method for detecting the position of deep cavities based on slurry diffusion theory, characterized in that: include: Arrange the detection holes based on the geological data of the detection site; Determine the slurry mix ratio, grout the bottom of the detection hole, and record the grouting parameters including slurry injection rate, grouting time, grouting pressure and slurry density; Based on the theoretical formula of slurry diffusion and the recorded grouting parameter information, the slurry injection rate time curve is calculated to determine whether the monitoring curve and the time curve coincide with each other. If they coincide with each other, the grouting parameters are deleted, the detection hole position is re-determined, and the above steps are repeated. If they do not coincide with each other, the grouting parameters are retained, and the positions of other detection holes are determined based on the position of the first detection hole. The above steps are repeated to obtain multiple sets of grouting parameters. Based on several groups of grouting parameters, the distance from each detection hole to the cavity and the cavity radius are calculated respectively; Based on the principle of spatial positioning, the spatial position of the cavity is determined by combining the distance from each detection hole to the cavity and the cavity radius.

2. According to claim 1, a method for detecting the position of deep cavities based on slurry diffusion theory is characterized in that: Arrange the detection hole based on the geological data of the detection site, including: obtain the geological data of the detection site, evaluate the possibility, development degree and potential location of the local cavity, and dig a detection hole in the vicinity of the cavity.

3. The method for detecting the position of deep cavities based on slurry diffusion theory according to claim 1 is characterized in that: Determine the mix ratio of the slurry, including: determine the mix ratio of the slurry based on the development degree of the detected ground fissures, when the opening of the nearby fissures is large, use thick slurry, and when the opening of the nearby fissures is small, use thin slurry.

4. The method for detecting the position of deep cavities based on slurry diffusion theory according to claim 1 is characterized in that: The theoretical formula for the slurry diffusion is as follows: （1）； In the formula, is the grouting time; is the number of monitoring points; is the slurry diffusion radius; is the monitoring value of slurry injection rate; For the The monitoring value of the slurry injection rate at the moment; V is the total volume of the slurry diffused in the space within a certain period of time; t is the moment of the grouting monitoring value; The grouting monitoring value a moment of time; The grouting monitoring value - 1 moment of time; （2）； According to Darcy's law, we can get: （3）； In the formula, is the permeability coefficient of slurry in the formation; is the slurry hydraulic gradient; is the spherical diffusion area of ​​the slurry; Integrating formula (3), we can get: （4）； In the formula, It is the pressure head of the slurry in the grouting hole; is the pressure head of groundwater; For the Calculated value of injection rate of moment-to-moment time history curve; is the radius of the grouting hole; Substituting formula (2) into formula (4), we can get: （5）。 5. The method for detecting the position of deep cavities based on slurry diffusion theory according to claim 1 is characterized in that: Determine whether the monitoring curve and the time curve coincide, including: comparing the slurry injection rate time curve with the detection curve. When the two are basically consistent with time and then show obvious differences, it indicates that the slurry flows through the hole. When the two are always basically consistent, it indicates that the slurry does not flow through the hole and the detection hole position needs to be re-determined.

6. A method for detecting the position of deep cavities based on slurry diffusion theory according to claim 5, characterized in that: When the slurry flows through the cavity, it includes the following stages: Phase I ): The slurry has not yet flowed into the cavity, and the time course curve is consistent with the monitoring curve; Phase II : The slurry flows into the cavity, and the time-history curve and the monitoring curve show obvious differences, and the former gradually becomes smaller than the latter. The difference between the two continues to expand with time. In formulas (1) and (2), let , calculate the distance from the detection hole to the hole ; Phase 3 : The cavities are filled with slurry. Since the slurry continues to flow along the cracks at this time, the theoretical formula of slurry diffusion is applicable again, so the calculation curve and the monitoring curve are distributed in parallel.

7. The method for detecting the position of deep cavities based on slurry diffusion theory according to claim 1, characterized in that: Hole radius The calculation method is as follows: （6）； （7）； In the formula, is the volume of the cavity, and the cavity is assumed to be spherical; is the calculated value of slurry injection rate; is the monitoring value of slurry injection rate; is the cavity radius; is the injection rate of the time history curve a moment of time; is the injection rate of the time history curve - 1 moment of time; Calculate the cavity radius for each detection hole , and then take the average ; when When ; They are respectively the hole radius corresponding to the first detection hole, the hole radius corresponding to the second detection hole, the hole radius corresponding to the third detection hole, and the hole radius corresponding to the fourth detection hole.

8. The method for detecting the position of deep cavities based on slurry diffusion theory according to claim 1, characterized in that: The calculation method of the hole position is as follows: Based on the principle of spatial positioning, the coordinates of the center point of the cavity obey the following set of equations: （8）； In the formula, The spatial coordinates of the bottom of the first detection hole; is the spatial coordinate of the bottom of the second detection hole; The spatial coordinates of the bottom of the third detection hole; is the spatial coordinate of the bottom of the fourth detection hole; is the distance between the bottom of the first detection hole and the cavity surface; is the distance between the bottom of the second detection hole and the cavity surface; is the distance between the bottom of the third detection hole and the cavity surface; is the distance between the bottom of the fourth detection hole and the cavity surface; is the coordinate position of the center point of the hole; is the average value of the cavity radius of multiple detection holes; The objective function of the hole center coordinates is as follows: （9）； By finding the gradient of formula (9), we can get: （10）； ； The gradient optimization method is used to solve the coordinates of the center point of the hole, and we can get: （11）； In the formula, is the number of iterations; is the step length.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the method for detecting the position of deep cavities based on slurry diffusion theory as described in any one of claims 1 to 7 is implemented.

10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for detecting the position of deep cavities based on slurry diffusion theory as claimed in any one of claims 1 to 7 is implemented.

Images