Sanding dolomite formation water-rich conditions of the surge prediction method and system

Abstract

The application discloses a method and system for predicting water inrush of sandified dolomite stratum under water-rich conditions, comprising: obtaining first detection data and second detection data; extracting first cross-section data and second cross-section data; the cross-section data is position data of a cross-section with discontinuity; comparing the first cross-section data and the second cross-section data, and taking a region with a difference greater than a preset value between the first cross-section data and the second cross-section data as a target drilling region; obtaining third cross-section data; forming cross-section fitting data of the target cross-section; and predicting construction water inrush risk. The method and system for predicting water inrush of sandified dolomite stratum under water-rich conditions effectively improve the advanced detection capability of sandified dolomite tunnel construction, improve construction safety and reduce construction risk by data fusion on the detection data in front of a tunnel face and accurate data acquisition and completion in combination with the advanced drilling technology.

Description

TECHNICAL FIELD

[0001] The present application relates to intelligent detection technology, in particular to a water-rich condition under the sanding dolomite stratum of the sudden gushing prediction method and system. BACKGROUND

[0002] The characteristics of dolomite sanding generally show the generation and development of fractures at the joint surfaces in the rock mass. Dolomite sanding first occurs from each structural surface in the rock mass, because these structural surfaces allow groundwater to enter the rock mass and further erode and expand from the structural surfaces to other places in the rock mass, eventually causing the entire rock mass to be eroded, i.e., forming sanding dolomite. In this process, both physical and chemical effects exist, complement each other, promote each other, and intensify each other. The geological characteristics of dolomite determine that the dolomite sanding process is uniformly carried out, i.e., the entire rock mass is uniformly decomposed.

[0003] Based on the above mechanism, the sanding dolomite tunnel is more prone to sudden gushing disasters during excavation, affecting construction safety. Although the advanced geological prediction for sanding dolomite tunnels has been greatly developed, there is still a lack of an all-digital accurate prediction method. SUMMARY

[0004] In order to at least overcome the above-mentioned deficiencies in the prior art, the purpose of the present application is to provide a sudden gushing prediction method and system for water-rich sanding dolomite stratum.

[0005] In a first aspect, the present application provides a sudden gushing prediction method for water-rich sanding dolomite stratum, comprising:

[0006] Obtain the detection data of the target section in front of the working face as first detection data through TSP, and obtain the detection data of the target section as second detection data through a geological radar;

[0007] Extract the section data of the target section from the first detection data as first section data, and extract the section data of the target section from the second detection data as second section data; the section data is the position data of the discontinuous section;

[0008] Compare the first section data and the second section data, and regard the area with a difference greater than a preset value between the first section data and the second section data as a target drilling area;

[0009] Drill ahead at the position corresponding to the target drilling area from the working face, and obtain the section data in the ahead drilling as third section data;

[0010] fitting the first section data and the second section data according to the third section data to form section fitting data of the target section;

[0011] obtaining section fitting data corresponding to multiple target sections, and predicting a construction inrush risk according to the multiple section fitting data.

[0012] In the implementation of the embodiments of the present application, a technical solution for comprehensively judging inrush risk based on multiple detection methods is provided. It is understood that the section in the target section refers to the scanning section during detection, and the section in the section data refers to the discontinuous surface existing in the surrounding rock in front. The echo signals detected by the TSP and the geological radar respectively can show the broken condition of the surrounding rock in front of the working face. In order to accurately utilize these data to judge the inrush disaster risk, the embodiments of the present application adopt a data fusion manner to compare the first section data and the second section data. Since the detection methods are different, the first section data and the second section data will have differences, so the embodiments of the present application also adopt actual borehole data to correct the first section data and the second section data.

[0013] Specifically, the regions with large differences between the first section data and the second section data are selected, and advanced boreholes are drilled at these positions. The rock core taken out of the borehole can accurately express the broken condition of the surrounding rock in front, and the first section data and the second section data can be data-fused to form the section fitting data of the target section based on this data. It is understood that the data fitting and fusion manner can adopt the fusion technology in the prior art, such as data interpolation, or curve fitting. The person skilled in the art can select the fitting manner according to the needs, and the embodiments of the present application do not make much repetition. The section fitting data of multiple target sections can roughly know the surrounding rock condition in front of the working face, and the inrush disaster risk of the surrounding rock in front can be analyzed according to the data. The embodiments of the present application fuse the detection data in front of the tunnel working face, and combine the advanced drilling technology to accurately obtain and complete the data, which effectively improves the advanced detection capability of the sanding dolomite tunnel construction, improves the construction safety, and reduces the construction risk.

[0014] In a possible implementation, fitting the first section data and the second section data according to the third section data to form section fitting data of the target section includes:

[0015] determining an expansion range, and expanding the expansion range with the target borehole region as the center to form a fitting range;

[0016] In the fitting range, a point where the first cross section data and the second cross section data coincide is selected as first intermediate cross section data; in a non-fitting range of the target cross section, a point where the first cross section data and the second cross section data coincide is selected as second intermediate cross section data;

[0017] A first cross section development curve is drawn in the fitting range with a cross section position in the third cross section data as a starting point and the first intermediate cross section data as an intermediate fitting point;

[0018] A second cross section development curve is drawn in the fitting range with an intersection of the first cross section development curve at a boundary of the fitting range as a starting point and the second intermediate cross section data as an intermediate fitting point;

[0019] The corresponding first cross section development curve and the second cross section development curve are connected as cross section fitting data of the target cross section.

[0020] In a possible implementation, drawing the first cross section development curve in the fitting range with the cross section position in the third cross section data as the starting point and the first intermediate cross section data as the intermediate fitting point includes:

[0021] When the first cross section development curve is drawn starting from the starting point, a corresponding intermediate fitting point is searched according to a cross section trend at the starting point of the first cross section data and the second cross section data, and a fitting curve is established between the starting point and the intermediate fitting point;

[0022] When the first cross section development curve is drawn starting from the intermediate fitting point, a corresponding next intermediate fitting point is searched according to a cross section trend at the starting point of the first cross section data and the second cross section data, and a fitting curve is established between the starting point and the intermediate fitting point.

[0023] In a possible implementation, the cross section fitting data corresponding to a plurality of target cross sections is acquired, and the construction gushing risk is predicted according to the plurality of cross section fitting data includes:

[0024] The detection data of the target cross section is acquired as third detection data by a transient electromagnetic method;

[0025] Material composition in the cross section fitting data is determined as second judgment reference data according to the third detection data, and cross section distribution data in the cross section fitting data is determined as first judgment reference data;

[0026] The construction gushing risk is predicted according to the first judgment reference data and the second judgment reference data.

[0027] In a possible implementation, judging the material composition in the section in the section fitting data according to the third detection data as second judgment reference data comprises:

[0028] Establishing a spatial coordinate matching relationship between the third detection data and the first judgment reference data;

[0029] Obtaining the apparent resistivity in the third detection data matched by the first judgment reference data;

[0030] Judging the second judgment reference data corresponding to the first judgment reference data according to the apparent resistivity.

[0031] In a second aspect, the embodiments of the present application provide a water-rich condition under the sanding dolomite formation of the water inrush prediction system, comprising:

[0032] The acquisition unit is configured to obtain the detection data of the target section in front of the working face through TSP as first detection data, and obtain the detection data of the target section through the geological radar as second detection data;

[0033] The extraction unit is configured to extract the section data of the target section from the first detection data as first section data, and extract the section data of the target section from the second detection data as second section data; the section data is the position data of the section with discontinuity;

[0034] The comparison unit is configured to compare the first section data and the second section data, and regard the region with a difference greater than a preset value between the first section data and the second section data as a target drilling region;

[0035] The drilling unit is configured to perform advanced drilling on the position corresponding to the target drilling region from the working face, and obtain the section data in the advanced drilling as third section data;

[0036] The fitting unit is configured to fit the first section data and the second section data according to the third section data, and form the section fitting data of the target section;

[0037] The judgment unit is configured to obtain the section fitting data corresponding to a plurality of target sections, and predict the construction water inrush risk according to a plurality of the section fitting data.

[0038] In a possible implementation, the fitting unit is further configured to:

[0039] Determine an expansion range, and expand the expansion range with the target drilling region as the center to form a fitting range;

[0040] In the fitting range, a point where the first cross section data and the second cross section data coincide is selected as first intermediate cross section data; in a non-fitting range of the target cross section, a point where the first cross section data and the second cross section data coincide is selected as second intermediate cross section data;

[0041] A first cross section development curve is drawn in the fitting range with the cross section position in the third cross section data as a starting point and the first intermediate cross section data as an intermediate fitting point;

[0042] A second cross section development curve is drawn in the fitting range with an intersection of the first cross section development curve at a boundary of the fitting range as a starting point and the second intermediate cross section data as an intermediate fitting point;

[0043] The corresponding first cross section development curve and the second cross section development curve are connected as cross section fitting data of the target cross section.

[0044] In a possible implementation, the fitting unit is further configured to:

[0045] When the first cross section development curve is drawn starting from the starting point, a corresponding intermediate fitting point is searched according to a cross section trend at the starting point of the first cross section data and the second cross section data, and a fitting curve is established between the starting point and the intermediate fitting point;

[0046] When the first cross section development curve is drawn starting from the intermediate fitting point, a corresponding next intermediate fitting point is searched according to a cross section trend at the starting point of the first cross section data and the second cross section data, and a fitting curve is established between the starting point and the intermediate fitting point.

[0047] In a possible implementation, the judging unit is further configured to:

[0048] The detection data of the target cross section is acquired by the transient electromagnetic method as third detection data through the acquisition unit;

[0049] The material composition in the cross section fitting data is judged as second judgment reference data according to the third detection data, and the cross section distribution data in the cross section fitting data is judged as first judgment reference data;

[0050] The construction gushing risk is predicted according to the first judgment reference data and the second judgment reference data.

[0051] In a possible implementation, the judging unit is further configured to:

[0052] A spatial coordinate matching relationship between the third detection data and the first judgment reference data is established;

[0053] acquire apparent resistivity in third detection data matched with the first judging reference data;

[0054] determine second judging reference data corresponding to the first judging reference data according to the apparent resistivity.

[0055] Compared with the prior art, the present application has the following advantages and beneficial effects:

[0056] The water-rich condition sanding dolomite stratum gushing prediction method and system can effectively improve the advanced detection capability for sanding dolomite tunnel construction, improve construction safety, and reduce construction risks. BRIEF DESCRIPTION OF DRAWINGS

[0057] The accompanying drawings, which are included to provide a further understanding of the embodiments of the application and are incorporated in and constitute a part of this application, illustrate embodiments of the application and together with the description serve to explain the principles of the application. In the drawings:

[0058] Figure 1 The figure is a schematic diagram of the method steps. DETAILED DESCRIPTION

[0059] To make the objectives, technical solutions, and advantages of the embodiments of the present application clearer, the following will be combined with the accompanying drawings to make a clear and complete description of the technical solutions in the embodiments of the present application. It should be understood that the accompanying drawings are only used for the purpose of illustration and description, and are not used to limit the protection scope of the present application. In addition, it should be understood that the schematic drawings are not drawn according to the actual proportions. The flowcharts in the present application show the operations implemented according to some embodiments of the present application. It should be understood that the operations of the flowcharts can not be implemented in sequence, and the steps without logical context relationship can be reversed in sequence or implemented simultaneously. In addition, one or more other operations can be added to the flowcharts or removed from the flowcharts under the guidance of the content of the present application.

[0060] In addition, the described embodiments are only some of the embodiments of the present application, not all the embodiments. The components of the embodiments of the present application described and shown in the accompanying drawings can be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present application provided in the accompanying drawings is not intended to limit the scope of the claimed present application, but only represents selected embodiments of the present application. Based on the embodiments of the present application, all other embodiments obtained by those skilled in the art without creative work fall within the scope of the present application.

[0061] Please refer to Figure 1A flowchart of a water-rich sand-dolomite stratum gushing prediction method provided by an embodiment of the present application is shown in FIG. 1. The water-rich sand-dolomite stratum gushing prediction method can include the following steps S1-S6.

[0062] S1: Obtain, by TSP, detection data of a target section in front of a working face as first detection data, and obtain, by a geological radar, detection data of the target section as second detection data;

[0063] S2: Extract, from the first detection data, section data of the target section as first section data, and extract, from the second detection data, section data of the target section as second section data; the section data is position data of a discontinuous section;

[0064] S3: Compare the first section data and the second section data, and regard a region in which the first section data and the second section data differ by more than a preset value as a target drilling region;

[0065] S4: Drilling ahead at a position corresponding to the target drilling region from the working face, and obtaining section data in the ahead drilling as third section data;

[0066] S5: Fitting the first section data and the second section data according to the third section data to form section fitting data of the target section;

[0067] S6: Obtaining section fitting data corresponding to a plurality of target sections, and predicting a construction gushing risk according to the plurality of section fitting data.

[0068] When the embodiments of the present application are implemented, a technical solution for comprehensively judging a gushing risk based on a plurality of detection methods is provided. The section data of a target section is detected by TSP and a geological radar. It should be understood that the section in the target section refers to a scanning section during detection, and the section in the section data refers to a discontinuous section in front of surrounding rock. The echo signals detected by TSP and the geological radar can show the broken condition of the surrounding rock in front of the working face. In order to accurately use the data to judge the gushing disaster risk, the embodiments of the present application compare the first section data and the second section data by using a data fusion method. Since the detection methods are different, the first section data and the second section data will differ. Therefore, the embodiments of the present application also correct the first section data and the second section data by using actual drilling data.

[0069] Specifically, a region with a large difference between the first cross-section data and the second cross-section data is selected, and a pre-probing hole is drilled at the position. The core taken from the pre-probing hole can accurately express the broken condition of the surrounding rock in front, and the first cross-section data and the second cross-section data can be fused to form cross-section fitting data of the target cross-section based on the data. It should be understood that the data fusion method can adopt the fusion technology in the prior art, such as data interpolation, or curve fitting. A person skilled in the art can select the fitting method according to the needs, and the present application embodiment does not make much repetition. The surrounding rock conditions in front of the working face can be roughly known through the cross-section fitting data of multiple target cross-sections, and the risk of sudden inrush disaster of the surrounding rock in front can be analyzed according to the data. The present application embodiment fuses the detection data in front of the tunnel working face, and combines the advanced drilling technology to obtain and complete accurate data, thereby effectively improving the advanced detection capability of the sanding dolomite tunnel construction, improving the construction safety, and reducing the construction risk.

[0070] In a possible implementation, fitting the first cross-section data and the second cross-section data according to the third cross-section data to form cross-section fitting data of the target cross-section comprises:

[0071] determining an expansion range and expanding the expansion range with the target pre-probing hole region as the center to form a fitting range;

[0072] in the fitting range, selecting a point where the first cross-section data and the second cross-section data coincide as first intermediate cross-section data; and in the non-fitting range of the target cross-section, selecting a point where the first cross-section data and the second cross-section data coincide as second intermediate cross-section data;

[0073] drawing a first cross-section development curve in the fitting range with the cross-section position in the third cross-section data as a starting point and the first intermediate cross-section data as an intermediate fitting point;

[0074] drawing a second cross-section development curve in the fitting range with an intersection of the first cross-section development curve at the boundary of the fitting range as a starting point and the second intermediate cross-section data as an intermediate fitting point;

[0075] connecting the corresponding first cross-section development curve and the second cross-section development curve as the cross-section fitting data of the target cross-section.

[0076] The embodiment of the application provides a specific fitting scheme. The inventor finds that the crushing process of a certain part is affected by the crushing of other surrounding rocks within a certain range in the process of crushing of the surrounding rock joints, and the crushing of other surrounding rocks within the range is also affected. Therefore, the range is determined as an expansion range, that is, it is determined that the cross section detected by the advanced exploration hole affects the development of the surrounding rock cracks in the expansion range, and the first cross section data and the second cross section data can be corrected in the range according to the accurate third cross section data. Similarly, for the points not in the fitting range, the influence of the third cross section data is generally not considered, and the second intermediate cross section data can be directly selected. For the cross section of the whole tunnel working face, since a certain number of advanced exploration holes can be drilled, the fitting range can cover most of the target cross section to ensure the accuracy of the final fitting result.

[0077] In the data fitting, the cross section position in the third cross section data is taken as the starting point of calibration for curve fitting. The third cross section data contains a plurality of cross section positions, and curve fitting is performed from each cross section position as the starting point to form the first cross section development curve in front of the working face. Similarly, the intersection point of the end of the first cross section development curve in the fitting range, that is, the boundary of the fitting range, can be taken as the starting point to draw the second cross section development curve in the non-fitting range. It should be understood that the change of the starting point position in the curve fitting in the fitting range will affect the distribution of the whole curve, so the influence of the cross section position in the third cross section data on the crushing of other surrounding rocks in the fitting range can be represented in this way.

[0078] In a possible implementation, the first cross section development curve is drawn in the fitting range with the cross section position in the third cross section data as the starting point and the first intermediate cross section data as the intermediate fitting point, including:

[0079] When the first cross section development curve is drawn starting from the starting point, a corresponding intermediate fitting point is searched according to the cross section trend at the starting point of the first cross section data and the second cross section data, and a fitting curve is established between the starting point and the intermediate fitting point;

[0080] When the first cross section development curve is drawn starting from the intermediate fitting point, a corresponding next intermediate fitting point is searched according to the cross section trend at the starting point of the first cross section data and the second cross section data, and a fitting curve is established between the starting point and the intermediate fitting point.

[0081] In a possible implementation, the cross section fitting data corresponding to a plurality of target cross sections is obtained, and the construction gushing risk is predicted according to a plurality of the cross section fitting data, including:

[0082] The detection data of the target cross section is obtained by the transient electromagnetic method as third detection data;

[0083] judging the material composition in the section in the section fitting data according to the third detection data as second judgment reference data, and judging the section distribution data in the section fitting data as first judgment reference data;

[0084] predicting the construction gushing risk according to the first judgment reference data and the second judgment reference data.

[0085] In the implementation of the embodiments of the present application, since the material filled in the fracture of the broken surrounding rock directly affects the gushing risk in tunnel construction, the transient electromagnetic method is introduced to detect the target section. The transient electromagnetic method can identify the water enrichment condition in the region, so as to judge the water content of the fracture filling material in each local region, and further judge the construction gushing risk.

[0086] For example, the judgment of the construction gushing risk can be performed by the following specific scheme:

[0087] According to the second judgment reference data, the target section is divided into a plurality of regions, and the water content of each region and the adjacent region is different;

[0088] According to the first judgment reference data and the water content in the region, the total water content in each region is calculated, and if the calculation result exceeds a preset value, it is judged that there is a construction gushing risk when the region is excavated.

[0089] In a possible implementation, judging the material composition in the section in the section fitting data according to the third detection data as second judgment reference data comprises:

[0090] Establishing a spatial coordinate matching relationship between the third detection data and the first judgment reference data;

[0091] Obtaining the apparent resistivity in the third detection data matched by the first judgment reference data;

[0092] Judging the second judgment reference data corresponding to the first judgment reference data according to the apparent resistivity.

[0093] Based on the same inventive concept, a gushing prediction system for sanding dolomite strata under water enrichment conditions is also provided, comprising:

[0094] The acquisition unit is configured to acquire detection data of a target section in front of a working face through TSP as first detection data, and acquire detection data of the target section through a geological radar as second detection data;

[0095] extracting, by an extraction unit, section data of a target section from the first detection data as first section data and extracting, by the extraction unit, section data of the target section from the second detection data as second section data, the section data being position data of a section where discontinuity exists;

[0096] comparing, by a comparison unit, the first section data and the second section data, and regarding a region where a difference between the first section data and the second section data is greater than a preset value as a target drilling region;

[0097] drilling, by a drilling unit, a pilot hole in advance from a position corresponding to the target drilling region on a working face, and obtaining section data in the pilot hole as third section data;

[0098] fitting, by a fitting unit, the first section data and the second section data according to the third section data to form section fitting data of the target section;

[0099] obtaining, by a judging unit, section fitting data corresponding to a plurality of target sections, and predicting a construction inrush risk according to a plurality of the section fitting data.

[0100] In a possible implementation, the fitting unit is further configured to:

[0101] determine an extension range, and extend the extension range to form a fitting range with the target drilling region as a center;

[0102] select, in the fitting range, a point where the first section data and the second section data coincide as first intermediate section data, and select, in a non-fitting range of the target section, a point where the first section data and the second section data coincide as second intermediate section data;

[0103] draw, with a section position in the third section data as a starting point and with the first intermediate section data as an intermediate fitting point, a first section development curve in the fitting range;

[0104] draw, with an intersection point of the first section development curve at a boundary of the fitting range as a starting point and with the second intermediate section data as an intermediate fitting point, a second section development curve in the fitting range;

[0105] connect the corresponding first section development curve and the second section development curve as the section fitting data of the target section.

[0106] In a possible implementation, the fitting unit is further configured to:

[0107] When the first profile development curve is drawn starting from the starting point, a corresponding intermediate fitting point is searched according to the profile trend at the starting point of the first profile data and the second profile data, and a fitting curve is established between the starting point and the intermediate fitting point;

[0108] When the first profile development curve is drawn starting from the intermediate fitting point, a corresponding next intermediate fitting point is searched according to the profile trend at the starting point of the first profile data and the second profile data, and a fitting curve is established between the starting point and the intermediate fitting point.

[0109] In a possible implementation, the judging unit is further configured to:

[0110] The detection data of the target profile is acquired by the transient electromagnetic method as third detection data through the acquisition unit;

[0111] The material composition in the profile fitting data is judged as second judging reference data according to the third detection data, and the profile distribution data in the profile fitting data is taken as first judging reference data;

[0112] The construction gushing risk is predicted according to the first judging reference data and the second judging reference data.

[0113] In a possible implementation, the judging unit is further configured to:

[0114] The spatial coordinate matching relationship between the third detection data and the first judging reference data is established;

[0115] The apparent resistivity in the third detection data matched by the first judging reference data is acquired;

[0116] The second judging reference data corresponding to the first judging reference data is judged according to the apparent resistivity.

[0117] Those skilled in the art can realize that the units and algorithm steps of each example described in combination with the embodiments disclosed herein can be realized by electronic hardware, computer software or a combination of both. In order to clearly illustrate the interchangeability of hardware and software, the components and steps of each example have been described in the above description in general terms. Whether the functions are realized in hardware or software depends on the specific application and design constraints of the technical solution. The skilled person can use different methods to realize the described functions for each specific application, but such implementation should not be considered beyond the scope of the present application.

[0118] In several embodiments provided in the present application, it should be understood that the disclosed apparatus and method can be implemented by other manners. For example, the apparatus embodiments described above are merely illustrative, for example, the division of the units is merely a logical function division, and actual implementation can have another division manner, for example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the coupling or direct coupling or communication connection between the units shown or discussed can be indirect coupling or communication connection through some interfaces, devices or units, and can also be electrical, mechanical or other form of connection.

[0119] The units described as separate components can or can not be physically separate, and it is obvious to those skilled in the art that the units and algorithm steps of each example described in combination with the embodiments disclosed herein can be realized by electronic hardware, computer software or a combination of both. In order to clearly illustrate the interchangeability of hardware and software, the components and steps of each example have been described in the above description in general terms. Whether the functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. The skilled person can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the present application.

[0120] In addition, each functional unit in each embodiment of the present application can be integrated into one processing unit, or each unit can exist physically, or two or more units can be integrated into one unit. The integrated unit can be realized in the form of hardware or software functional unit.

[0121] The integrated unit, if realized in the form of software functional unit and sold or used as an independent product, can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application essentially or the part that contributes to the prior art, or the whole or part of the technical solutions can be embodied in the form of a software product, which is stored in a storage medium and includes a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a grid device, etc.) to execute all or part of the steps of the methods described in each embodiment of the present application. The foregoing storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various program code storage media.

[0122] The above detailed description of the specific embodiments of the present application has been given to understand the purpose, technical solutions and beneficial effects of the present application. It should be understood that the above description is only a specific embodiment of the present application and is not used to limit the protection scope of the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method for predicting kick in a watered dolomitic formation, characterized in that, The method comprises: obtaining, by TSP, detection data of a target section in front of a working face as first detection data, and obtaining, by a geological radar, detection data of the target section as second detection data; extracting, from the first detection data, section data of the target section as first section data, and extracting, from the second detection data, section data of the target section as second section data; the section data is position data of a section where a discontinuity exists; comparing the first section data and the second section data, and regarding a region where a difference between the first section data and the second section data is greater than a preset value as a target drilling region; drilling ahead at a position corresponding to the target drilling region from the working face, and obtaining section data in the ahead drilling as third section data; fitting the first section data and the second section data according to the third section data to form section fitting data of the target section; obtaining section fitting data corresponding to a plurality of target sections, and predicting a construction gushing risk according to the plurality of section fitting data.

2. The method of kick prediction for watered sanding dolomite formations according to claim 1, characterized in that, The fitting of the first section data and the second section data according to the third section data to form the section fitting data of the target section comprises: determining an extension range, and extending the extension range from the center of the target drilling region to form a fitting range; in the fitting range, selecting a point where the first section data and the second section data coincide as first intermediate section data; and in a non-fitting range of the target section, selecting a point where the first section data and the second section data coincide as second intermediate section data; drawing a first section development curve in the fitting range with the section position in the third section data as a starting point and the first intermediate section data as an intermediate fitting point; drawing a second section development curve in the fitting range with an intersection of the first section development curve at a boundary of the fitting range as a starting point and the second intermediate section data as an intermediate fitting point; connecting the corresponding first section development curve and the second section development curve as the section fitting data of the target section.

3. The method of kick prediction for watered sanding dolomite formations according to claim 2, characterized in that, The drawing of the first section development curve in the fitting range with the section position in the third section data as a starting point and the first intermediate section data as an intermediate fitting point comprises: when drawing the first section development curve from the starting point, searching for a corresponding intermediate fitting point according to a section trend at the starting point of the first section data and the second section data, and establishing a fitting curve between the starting point and the intermediate fitting point; when drawing the first section development curve from the intermediate fitting point, searching for a corresponding next intermediate fitting point according to a section trend at the starting point of the first section data and the second section data, and establishing a fitting curve between the starting point and the intermediate fitting point.

4. The method of kick prediction of a watered-out dolomitic formation according to claim 2, characterized in that, The obtaining of the section fitting data corresponding to a plurality of target sections and the predicting of a construction gushing risk according to the plurality of section fitting data comprises: obtaining, by transient electromagnetic method, detection data of the target section as third detection data; judging material composition in the cross section fitting data according to the third detection data as second judgment reference data, and judging cross section distribution data in the cross section fitting data as first judgment reference data; predicting construction gushing risk according to the first judgment reference data and the second judgment reference data.

5. The method of kick prediction of a watered-out dolomitic formation according to claim 4, characterized in that, judging material composition in the cross section fitting data according to the third detection data as second judgment reference data includes: establishing spatial coordinate matching relationship between the third detection data and the first judgment reference data; obtaining apparent resistivity in the third detection data matched by the first judgment reference data; judging second judgment reference data corresponding to the first judgment reference data according to the apparent resistivity.

6. A kick prediction system for watered-out dolomitic formations using the method of any one of claims 1 to 5, characterized in that, includes: an acquisition unit configured to acquire detection data of a target cross section in front of a working face through TSP as first detection data, and acquire detection data of the target cross section through ground penetrating radar as second detection data; an extraction unit configured to extract cross section data of the target cross section from the first detection data as first cross section data, and extract cross section data of the target cross section from the second detection data as second cross section data; the cross section data is position data of a cross section with discontinuity; a comparison unit configured to compare the first cross section data and the second cross section data, and regard a region with a difference greater than a preset value between the first cross section data and the second cross section data as a target drilling region; a drilling unit configured to perform advanced drilling at a position corresponding to the target drilling region from the working face, and acquire cross section data in the advanced drilling as third cross section data; a fitting unit configured to fit the first cross section data and the second cross section data according to the third cross section data, and form cross section fitting data of the target cross section; a judgment unit configured to acquire cross section fitting data corresponding to a plurality of target cross sections, and predict construction gushing risk according to the plurality of cross section fitting data.

7. The kick prediction system of a watered-out dolomitic formation according to claim 6, characterized in that, the fitting unit is further configured to: determine an expansion range, and expand the expansion range with the target drilling region as a center to form a fitting range; select a point where the first cross section data and the second cross section data coincide as first intermediate cross section data in the fitting range; and select a point where the first cross section data and the second cross section data coincide as second intermediate cross section data in a non-fitting range of the target cross section; draw a first cross section development curve in the fitting range with a cross section position in the third cross section data as a starting point and with the first intermediate cross section data as an intermediate fitting point; draw a second cross section development curve in the fitting range with an intersection of the first cross section development curve at a boundary of the fitting range as a starting point and with the second intermediate cross section data as an intermediate fitting point; connect the corresponding first cross section development curve and the second cross section development curve as cross section fitting data of the target cross section.

8. The kick prediction system of a watered-out dolomitic formation sanding according to claim 7, characterized in that, the fitting unit is further configured to: When the first profile development curve is drawn starting from the starting point, a corresponding intermediate fitting point is searched according to the profile trend at the starting point of the first profile data and the second profile data, and a fitting curve is established between the starting point and the intermediate fitting point; When the first profile development curve is drawn starting from the intermediate fitting point, a corresponding next intermediate fitting point is searched according to the profile trend at the starting point of the first profile data and the second profile data, and a fitting curve is established between the starting point and the intermediate fitting point.

9. The kick prediction system of a watered-out dolomitic formation sanding in accordance with claim 7, wherein, The judging unit is further configured to: acquire, by the transient electromagnetic method, the detection data of the target profile as third detection data; judge, according to the third detection data, the material composition in the profile fitting data as second judgment reference data, and judge the profile distribution data in the profile fitting data as first judgment reference data; predict the construction gushing risk according to the first judgment reference data and the second judgment reference data.

10. The kick prediction system of a watered-out dolomitic formation sanding according to claim 9, characterized in that, The judging unit is further configured to: establish a spatial coordinate matching relationship between the third detection data and the first judgment reference data; acquire the apparent resistivity in the third detection data matched by the first judgment reference data; judge the second judgment reference data corresponding to the first judgment reference data according to the apparent resistivity.

Images