Coal mine goaf overburden rock three zone identification method and device
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- CHINA COAL SCI & ENG ECOLOGICAL ENVIRONMENT TECH CO LTD
- Filing Date
- 2023-05-11
- Publication Date
- 2026-06-09
Smart Images
Figure CN116498308B_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of coal mining technology, and in particular to a method and apparatus for identifying the three zones of overburden in coal mine goaf. Background Technology
[0002] In the remediation of coal mine goaf areas, determining the distribution of underground goaf areas and the extent of overlying strata damage is a crucial exploration task. Among various exploration methods, drilling, as the most direct method of exploration and verification, plays a vital role in determining the distribution and extent of the three zones (bending zone, fault zone, and collapse zone) of the overlying strata in the goaf area. In related technologies, the delineation of these three zones is usually carried out by specialized technicians after drilling is completed, analyzing observation data, combining field records, and previous geological data and standards to comprehensively determine the boundaries of the three zones. As the number of boreholes increases, the workload of technicians will increase significantly, resulting in a simultaneous increase in time and labor costs, and even affecting the project schedule. In addition, the level and experience of technicians can lead to differences in the analysis results among different personnel. Summary of the Invention
[0003] This application aims to at least partially address one of the technical problems in the related art.
[0004] Therefore, the first aspect of this application proposes a method for identifying the three zones of overburden in coal mine goaf, including:
[0005] Obtain drilling data, empirical values of fault zone height, empirical values of caving zone height, borehole depth, bottom plate position constraints, and coal mining thickness from the drilling site;
[0006] Based on the Bayesian change point detection method, the change points of the drilling data are obtained;
[0007] The lower limit depth of the bending zone is obtained based on the borehole depth, the empirical value of the fracture zone height, the empirical value of the collapse zone height, and the change points of the drilling data.
[0008] The lower limit depth of the caving zone is determined based on the borehole depth, the bottom plate position constraints, the coal mining thickness, and the change points of the drilling data.
[0009] The lower limit depth of the fracture zone is determined based on the empirical values of the fracture zone height, the empirical values of the collapse zone height, the lower limit depth of the bending zone, the lower limit depth of the collapse zone, and the change points of the drilling data.
[0010] In some embodiments of this application, the change points of the drilling data include multiple change points of flushing fluid leakage; obtaining the lower limit depth of the bending zone based on the borehole depth, the empirical value of the fracture zone height, the empirical value of the collapse zone height, and the change points of the drilling data includes: determining a first depth range based on the empirical values of the fracture zone height and the collapse zone height; determining multiple depths corresponding to the multiple change points of flushing fluid leakage; and determining the minimum depth among the multiple depths that falls within the first depth range as the lower limit depth of the bending zone.
[0011] In some embodiments of this application, the change points of the drilling data include multiple change points of the drilling speed; determining the lower limit depth of the caving zone based on the borehole depth, the bottom plate position constraints, the coal mining thickness, and the change points of the drilling data includes: determining a second depth range based on the borehole depth and the bottom plate position constraints; and determining multiple depths corresponding to the multiple change points of the drilling speed.
[0012] The smallest depth among the plurality of depths that falls within the second depth range is determined as the depth of the goaf floor; the lower limit depth of the caving zone is obtained based on the depth of the goaf floor and the coal thickness.
[0013] In some embodiments of this application, determining the lower limit depth of the fracture zone based on the empirical value of the fracture zone height, the empirical value of the collapse zone height, the lower limit depth of the curvature zone, the lower limit depth of the collapse zone, and the change point includes: determining a third depth range based on the empirical value of the fracture zone height, the empirical value of the collapse zone height, the lower limit depth of the curvature zone, and the lower limit depth of the collapse zone; when there is a change point within the third depth range, obtaining the error value from each change point within the third depth range to the boundary of the third depth range; and determining the depth corresponding to the change point with the smallest error value as the lower limit depth of the fracture zone.
[0014] In some embodiments of this application, the method further includes: when there is no change point within the third depth range, taking the midpoint of the third depth range as the lower limit depth of the fracture zone.
[0015] In some embodiments of this application, the method for obtaining change points in the drilling data based on Bayesian change point detection includes: filtering the drilling data, identifying change points in the filtered drilling data using the Bayesian change point detection method, and obtaining the change points.
[0016] The second aspect of this application discloses a device for identifying the three zones of overburden in a coal mine goaf, comprising:
[0017] The first acquisition module is used to acquire drilling data, empirical values of fault zone height, empirical values of caving zone height, borehole depth, bottom plate position constraints, and coal mining thickness at the drilling site.
[0018] The second acquisition module is used to acquire the change points of the drilling data based on the Bayesian change point detection method;
[0019] The first determining module is used to obtain the lower limit depth of the bending zone based on the borehole depth, the empirical value of the fracture zone height, the empirical value of the collapse zone height, and the change points of the drilling data;
[0020] The second determining module is used to determine the lower limit depth of the caving zone based on the borehole depth, the bottom plate position constraints, the coal mining thickness, and the change points of the drilling data;
[0021] The third determining module is used to determine the lower limit depth of the fracture zone based on the empirical value of the fracture zone height, the empirical value of the collapse zone height, the lower limit depth of the bending zone, the lower limit depth of the collapse zone, and the change points of the drilling data.
[0022] In some embodiments of this application, the second acquisition module is specifically used to: filter the drilling data, identify change points in the filtered drilling data using the Bayesian change point detection method, and obtain the change points.
[0023] A third aspect of this application provides an electronic device comprising: a processor; and a memory for storing processor-executable instructions; wherein the instructions are executed by the processor to enable the processor to perform the method described in the first aspect above.
[0024] The fourth aspect of this application provides a non-transitory computer-readable storage medium, wherein when instructions in the storage medium are executed by a processor of an electronic device, the electronic device is able to perform the method described in the first aspect above.
[0025] According to the technical solution of this application, by detecting changes in on-site drilling data and combining geological data and empirical formulas to identify the boundaries of the three zones, the accuracy of the three-zone identification can be effectively improved. This method is unaffected by regional or borehole differences and has a wide range of applications. Furthermore, this application can reduce the workload of technical personnel, saving time and manpower costs.
[0026] Additional aspects and advantages of this application will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of this application. Attached Figure Description
[0027] The above and / or additional aspects and advantages of this application will become apparent and readily understood from the following description of the embodiments taken in conjunction with the accompanying drawings, wherein:
[0028] Figure 1 A flowchart illustrating a method for identifying the three zones of overburden in a coal mine goaf, as provided in an embodiment of this application;
[0029] Figure 2 This is a schematic diagram showing the distribution of the three zones of overlying strata in a goaf area;
[0030] Figure 3 This is a flowchart illustrating another method for identifying the three zones of overburden in a coal mine goaf, as provided in an embodiment of this application.
[0031] Figure 4a A flushing fluid leakage curve provided for an embodiment of this application;
[0032] Figure 4b Drilling rate curves provided for embodiments of this application;
[0033] Figure 5a This is a schematic diagram showing the changes in flushing fluid leakage in the embodiments of this application;
[0034] Figure 5b This is a schematic diagram of drilling rate variation points provided in an embodiment of this application;
[0035] Figure 6 This is a schematic diagram of the three-band division results provided in an embodiment of this application;
[0036] Figure 7 A comparison diagram of the three-band identification boundary and the actual boundary provided for the embodiments of this application;
[0037] Figure 8 This is a schematic diagram of a coal mine goaf overburden three-zone identification device provided in an embodiment of this application. Detailed Implementation
[0038] The embodiments of this application are described in detail below. Examples of these embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and intended to explain this application, and should not be construed as limiting this application.
[0039] This application proposes a method and apparatus for identifying the three zones of overburden in a coal mine goaf. Specifically, the method and apparatus for identifying the three zones of overburden in a coal mine goaf according to embodiments of this application are described below with reference to the accompanying drawings.
[0040] Figure 1 This is a flowchart illustrating a method for identifying the three zones of overburden in a coal mine goaf, as provided in an embodiment of this application. Figure 1As shown, the method for identifying the three zones of overburden in the goaf of this coal mine includes the following steps:
[0041] Step 101: Obtain drilling data, empirical values of fault zone height, empirical values of caving zone height, borehole depth, bottom plate position constraints, and coal mining thickness from the drilling site.
[0042] Figure 2 This is a schematic diagram showing the distribution of the three zones of overlying strata in a goaf. (For example...) Figure 2 As shown in the figure, the area includes a bend zone 1, a fault zone 2, and a collapse zone 3. Here, d0 is the surface depth, d1 is the lower limit depth of the bend zone, d2 is the lower limit depth of the fault zone, d3 is the lower limit depth of the collapse zone, and d4 is the depth of the goaf floor.
[0043] Drilling data at the drilling site can be automatically acquired in real time through a field monitoring system. In some embodiments of this application, drilling data may include flushing volume, leakage volume, drilling speed, drilling parameters, etc. The empirical values for fault zone height and caving zone height depend on the characteristics of the overlying strata; the empirical value for fault zone height H′... li and empirical value of caving zone height H′ m The empirical formula can be referenced as follows:
[0044]
[0045]
[0046] Where M is the coal mining thickness and k is the coefficient of fracture expansion of the roof strata.
[0047] Step 102: Based on the Bayesian change point detection method, obtain the change points of the drilling data.
[0048] In one possible implementation, assuming that the data in all three bands conform to a Gaussian distribution, and that the mean and precision are unknown, the corresponding expression for the exponential family conjugate prior probability density function is:
[0049]
[0050] Given initial parameter values α = 0.1, β = 0.01, μ = 0, and λ = 1, based on the characteristics of the conjugate distribution of the exponential family, the expression for parameter iterative update can be obtained as follows:
[0051]
[0052]
[0053]
[0054]
[0055]
[0056]
[0057]
[0058]
[0059] Where n = 1, For the newly observed data point value x t+1 t is the step size of the data point, and r is the running length of the current change point interval.
[0060] Given an exponential family of discrete priors λ gap =50, the risk function H is expressed as:
[0061]
[0062] Based on the above parameters, changes in drilling data can be detected.
[0063] Optionally, in some embodiments of this application, the drilling data can be filtered, such as by median filtering, with the filter window size set to 5. The Bayesian change point detection method is used to identify change points in the filtered drilling data, thereby reducing the interference of noise data and improving the accuracy of change point identification.
[0064] Step 103: Based on the empirical values of borehole depth, fracture zone height, collapse zone height, and changes in drilling data, obtain the lower limit depth of the bending zone.
[0065] In one possible implementation, the depth range of the lower limit depth d1 of the bending zone can be determined based on empirical values of the fault zone height, the caving zone height, and the borehole depth. Then, based on the obtained depth range of the lower limit depth d1 and the detected points of change, the lower limit depth d1 of the bending zone is determined.
[0066] Step 104: Determine the lower limit depth of the caving zone based on the borehole depth, bottom plate position constraints, coal mining thickness, and changes in drilling data.
[0067] It should be noted that the floor position constraint condition refers to the depth of the goaf floor being within a certain distance above the borehole depth. In one possible implementation, the goaf floor depth d4 can be determined based on the borehole depth, the floor position constraint condition, and the point of change. The lower limit depth d3 of the caving zone is then determined based on the goaf floor depth d4 and the coal thickness.
[0068] Step 105: Determine the lower limit depth of the fault zone based on the empirical values of the fault zone height, the empirical values of the caving zone height, the lower limit depth of the bending zone, the lower limit depth of the caving zone, and the change points of the drilling data.
[0069] In some embodiments of this application, the depth range of the lower limit depth d2 of the fracture zone can be determined based on empirical values of the fracture zone height, the collapse zone height, the lower limit depth d1 of the bend zone, and the lower limit depth d3 of the collapse zone. If there are variation points in the drilling data that fall within the depth range of the lower limit depth d2 of the fracture zone, then the minimum depth value corresponding to the variation point falling within that depth range can be determined as the lower limit depth d2 of the fracture zone.
[0070] It should be noted that, since the coal mine goaf overburden three-zone identification method described in this application embodiment only relies on pre-drilling geological data, empirical formulas and data observed during on-site drilling, this method can not only automatically identify the three-zone boundary depth in real time, but also be used for offline identification of the three-zone boundary depth.
[0071] The coal mine goaf overburden three-zone identification method according to the embodiments of this application identifies the three-zone boundaries by detecting changes in on-site drilling data and combining geological data and empirical formulas. This effectively improves the accuracy and efficiency of three-zone identification, is unaffected by regional or borehole differences, and has a wide range of applications. Furthermore, this application reduces the workload of technical personnel, saving time and manpower costs.
[0072] Figure 3 This is a flowchart illustrating another method for identifying the three zones of overburden in a coal mine goaf, as provided in an embodiment of this application. Figure 3 As shown, the method for identifying the three zones of overburden in the goaf of this coal mine includes the following steps:
[0073] Step 301: Obtain drilling data, empirical values of fault zone height, empirical values of caving zone height, borehole depth, bottom plate position constraints, and coal mining thickness from the drilling site.
[0074] Step 302: Based on the Bayesian change point detection method, obtain the change points of the drilling data. These change points include multiple change points of flushing fluid loss and multiple change points of drilling speed.
[0075] Step 303: Determine the first depth range based on the empirical values of the fault zone height and the caving zone height.
[0076] As an example, given a borehole depth d, the lower limit depth d1 of the bend zone is distributed within a certain range H above the borehole depth d. u The calculation method is as follows:
[0077] H u =(H′) li +H′ m )×(1+n)
[0078] Here, is defined as the window index, determined based on previous drilling data. Based on the range H... u The first depth range can be determined as dH u ~d.
[0079] Step 304: Determine the multiple depths corresponding to multiple points of change in flushing fluid leakage.
[0080] Step 305: Determine the minimum depth that falls within the first depth range from the multiple depths as the lower limit depth of the bending zone.
[0081] Step 306: Determine the second depth range based on the borehole depth and the positional constraints of the base plate.
[0082] As an example, suppose the constraint on the location of the bottom plate is that the depth d4 of the bottom plate in the goaf is H above the design hole depth. d Within m. Then the second depth range where the bottom plate depth d4 of the goaf is located is dH. d ~d.
[0083] Step 307: Determine the multiple depths corresponding to the multiple points of change in drilling speed.
[0084] Step 308: Determine the depth of the goaf floor that falls within the second depth range from the multiple depths.
[0085] Step 309: Obtain the lower limit depth of the caving zone based on the depth of the goaf floor and the thickness of the coal seam.
[0086] As an example, if the coal seam thickness is given as M before drilling, then the lower limit depth of the caving zone, d3, is d4-M.
[0087] Step 310: Determine the third depth range based on the empirical values of the fault zone height, the caving zone height, the lower limit depth of the bending zone, and the lower limit depth of the caving zone.
[0088] In some embodiments of this application, the height of the fracture zone is based on the empirical value H′. li Empirical value of landslide zone height H′ m The lower limit depth d1 of the bend zone and the lower limit depth d3 of the collapse zone can be used to determine the upper limit depth min(d′) of the third depth range. 21 ,d′ 22 ) and lower bound depth max(d′ 21 ,d′ 22 ).in,
[0089]
[0090] Step 311: When there are change points within the third depth range, obtain the error value from each change point within the third depth range to the boundary of the third depth range.
[0091] It should be noted that the change points in this step include all change points in the drill-down data. As an example, assume that d... c This represents the depth value at each point of change falling within the third depth range. The formula for calculating the error value from each point of change within the third depth range to the boundary of the third depth range can be found below:
[0092]
[0093] Step 312: Determine the depth corresponding to the point of smallest error change as the lower limit depth of the fault zone.
[0094] Step 313: When there is no change point within the third depth range, the midpoint of the third depth range is taken as the lower limit depth of the fault zone.
[0095] As an example, if there are no points of change within the third depth range, the formula for calculating d2 as the lower limit depth of the fault zone can be found as follows:
[0096] d2=(d′ 21 -d′ 22 ) / 2
[0097] In the embodiments of this application, steps 301-302 can be implemented in any of the embodiments of this application. This application does not make specific limitations on this and will not elaborate further.
[0098] The coal mine goaf overburden three-zone identification method according to the embodiments of this application identifies the three-zone boundaries by detecting changes in flushing volume, leakage, and drilling speed, combined with geological data and empirical formulas. This further improves the accuracy and efficiency of three-zone identification, making the method unaffected by regional or borehole differences and with a wide range of applications. Furthermore, this application can reduce the workload of technical personnel, saving time and manpower costs.
[0099] To better understand the method for identifying the three zones of overburden in coal mine goaf proposed in this application, an example is provided below to describe the method. Taking an exploration well in a certain mining area as an example, the roof of the goaf is composed of sandy mudstone and mudstone, with average uniaxial saturated compressive strengths of 40.77 MPa and 31.66 MPa, respectively, belonging to relatively hard rocks. The stratum dip angle α = 6°, the coefficient of fragmentation of the roof strata is k = 1.35, and the coal seam thickness is M = 2.8 m. The exploration borehole depth is 570 m, and the coal seam floor is H above the designed borehole depth. d =Within 10m (i.e., the constraint condition of the base plate position). Select the leakage rate of the on-site flushing fluid ( Figure 4a ) and drilling speed monitoring data ( Figure 4b ).
[0100] Based on the Bayesian change point detection method, the change points of flushing fluid leakage are obtained. Figure 5a ) and the point of change in drilling speed ( Figure 5b ).like Figure 5a As shown, the flushing fluid loss data yielded four variation points at depths of 494.22m, 511.36m, 515.36m, and 543.5m. The drilling speed data identified six variation points at depths of 269.72m, 505.68m, 510.42m, 515.71m, 544.26m, and 562.87m.
[0101] Select the empirical formulas for fault zones and caving zones from the standard to determine the empirical value H′ of the fault zone height. li and empirical value of caving zone height H′ m :
[0102]
[0103]
[0104] H′ li Taking the average of the upper and lower limits, and setting the window exponent n = 0.8, then:
[0105] H u =(52.24+8.04)×1.8=108.504m
[0106] According to H u The first depth range can be determined to be between 461.496 and 570 m. The point where the minimum depth of the flushing fluid leakage falls within this first depth range is defined as the lower limit depth d1 of the bend. Figure 5a It can be seen that the minimum depth change point within the range of 460 to 570 meters is 494.22 meters. Therefore, the lower limit depth d1 of the curved zone is determined to be 494.22 meters, and the height of the curved zone is 494.22 meters.
[0107] Based on the positional constraints of the bottom plate, it can be seen that the coal seam bottom plate extends upwards from the borehole depth H. d If the depth is within 10m, then the second depth range where the goaf floor depth d4 is located is between 560 and 570m. The point where the drilling speed variation falls within this second depth range and has the smallest depth is determined as the goaf floor depth d4. Figure 5b It can be seen that the minimum depth change point within the range of 560-570m is 562.87m, so the depth d4 of the goaf floor is determined to be 562.87m. Subtracting the coal mining thickness of 2.8m from the depth d4 of the goaf floor, the lower boundary d3 of the caving zone is obtained as 560.07m.
[0108] The lower limit depth of the bending zone d1, the lower limit depth of the caving zone d3, and the empirical value of the fault zone height H′ are derived from the above.li and empirical value of caving zone height H′ m It can be determined that the lower limit depth d2 of the fault zone is located within the third depth range of 551.46–552.03 m. Figure 5a and Figure 5b It can be seen that there are no points of change within this boundary. Therefore, by taking the average of the two boundary values, we can obtain that the lower limit depth d2 of the fault zone is 551.745m, the height of the fault zone is 57.525m, and the height of the collapse zone is 8.325m.
[0109] At this point, the three altitude zones and their corresponding depth boundaries have been determined. (See the results for reference.) Figure 6 For comparison with actual results, please refer to [the provided text]. Figure 7 The identification errors for the three zones were +0.67m, +0.645m, and +0.52m, respectively, while the errors for the heights of the bending zone, fracture zone, and collapse zone were +0.136%, -0.043%, and -1.479%, respectively. Application results demonstrate that the proposed method for identifying the three zones of overlying strata in coal mine goaf areas is accurate and feasible.
[0110] Figure 8 This is a schematic diagram of a coal mine goaf overburden three-zone identification device provided in an embodiment of this application. Figure 8 As shown, the coal mine goaf overburden three-zone identification device includes: a first acquisition module 801, a second acquisition module 802, a first determination module 803, a second determination module 804, and a third determination module 805.
[0111] in,
[0112] The first acquisition module 801 is used to acquire drilling data, empirical values of fault zone height, empirical values of caving zone height, borehole depth, bottom plate position constraints, and coal mining thickness at the drilling site.
[0113] The second acquisition module 802 is used to acquire the change points of the drilling data based on the Bayesian change point detection method.
[0114] The first determining module 803 is used to obtain the lower limit depth of the bending zone based on the empirical values of borehole depth, fracture zone height, collapse zone height, and changes in drilling data.
[0115] The second determining module 804 is used to determine the lower limit depth of the caving zone based on the borehole depth, bottom plate position constraints, coal mining thickness, and changes in drilling data.
[0116] The third determining module 805 is used to determine the lower limit depth of the fault zone based on the empirical values of the fault zone height, the empirical values of the caving zone height, the lower limit depth of the bending zone, the lower limit depth of the caving zone, and the change points of the drilling data.
[0117] In some embodiments of this application, the second acquisition module 802 is specifically used to: filter the drilling data, identify change points in the filtered drilling data using a Bayesian change point detection method, and obtain change points.
[0118] In some embodiments of this application, the change points of the drilling data include multiple change points of flushing fluid leakage; the first determining module 803 is specifically used to: determine a first depth range based on empirical values of fracture zone height and collapse zone height; determine multiple depths corresponding to multiple change points of flushing fluid leakage; and determine the minimum depth among the multiple depths that falls within the first depth range as the lower limit depth of the bend zone.
[0119] In some embodiments of this application, the change points of drilling data include multiple change points of drilling speed; the second determining module 804 is specifically used to: determine a second depth range based on the borehole depth and bottom plate position constraints; determine multiple depths corresponding to multiple change points of drilling speed; determine the smallest depth among the multiple depths that falls within the second depth range as the bottom plate depth of the goaf; and obtain the lower limit depth of the caving zone based on the bottom plate depth of the goaf and the coal thickness.
[0120] In some embodiments of this application, the third determining module 805 is specifically used to: determine a third depth range based on empirical values of fracture zone height, collapse zone height, lower limit depth of bending zone, and lower limit depth of collapse zone; when there are change points within the third depth range, obtain the error value from each change point within the third depth range to the boundary of the third depth range; and determine the depth corresponding to the change point with the smallest error value as the lower limit depth of the fracture zone.
[0121] Regarding the apparatus in the above embodiments, the specific manner in which each module performs its operation has been described in detail in the embodiments related to the method, and will not be elaborated upon here.
[0122] To implement the above embodiments, this application also proposes an electronic device, including: a processor and a memory for storing processor-executable instructions. These instructions are executed by the processor to enable the processor to perform the aforementioned method for identifying the three zones of overburden in coal mine goafs.
[0123] To implement the above embodiments, this application also proposes a non-transitory computer-readable storage medium, which, when the instructions in the storage medium are executed by the processor of an electronic device, enables the electronic device to perform the aforementioned method for identifying the three zones of overburden in coal mine goaf.
[0124] In the description of this specification, the references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., refer to specific features, structures, materials, or characteristics described in connection with that embodiment or example, which are included in at least one embodiment or example of this application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.
[0125] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of that feature. In the description of this application, "multiple" means at least two, such as two, three, etc., unless otherwise explicitly specified.
[0126] Any process or method description in the flowchart or otherwise herein can be understood as representing a module, segment, or portion of code comprising one or more executable instructions for implementing custom logic functions or processes, and the scope of the preferred embodiments of this application includes additional implementations in which functions may be performed not in the order shown or discussed, including substantially simultaneously or in reverse order depending on the functions involved, as should be understood by those skilled in the art to which embodiments of this application pertain.
[0127] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (such as a computer-based system, a processor-included system, or other system that can fetch and execute instructions from, an instruction execution system, apparatus, or device). For the purposes of this specification, "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transmit programs for use by, or in conjunction with, an instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of computer-readable media include: an electrical connection having one or more wires (electronic device), a portable computer disk drive (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and editable read-only memory (EPROM or flash memory), fiber optic devices, and portable optical disc read-only memory (CDROM). Alternatively, the computer-readable medium may be paper or other suitable media on which the program can be printed, since the program can be obtained electronically, for example, by optically scanning the paper or other medium, followed by editing, interpreting, or otherwise processing as necessary, and then stored in a computer memory.
[0128] It should be understood that various parts of this application can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiments, multiple steps or methods can be implemented using software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware as in another embodiment, it can be implemented using any one or a combination of the following techniques known in the art: discrete logic circuits having logic gates for implementing logical functions on data signals, application-specific integrated circuits (ASICs) having suitable combinational logic gates, programmable gate arrays (PGAs), field-programmable gate arrays (FPGAs), etc.
[0129] Those skilled in the art will understand that all or part of the steps of the methods in the above embodiments can be implemented by a program instructing related hardware. The program can be stored in a computer-readable storage medium, and when executed, the program includes one or a combination of the steps of the method embodiments.
[0130] Furthermore, the functional units in the various embodiments of this application can be integrated into a processing module, or each unit can exist physically separately, or two or more units can be integrated into a module. The integrated module can be implemented in hardware or as a software functional module. If the integrated module is implemented as a software functional module and sold or used as an independent product, it can also be stored in a computer-readable storage medium.
[0131] The storage medium mentioned above can be a read-only memory, a disk, or an optical disk, etc. Although embodiments of this application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting this application. Those skilled in the art can make changes, modifications, substitutions, and variations to the above embodiments within the scope of this application.
Claims
1. A method for identifying the three zones of overburden in a coal mine goaf, characterized in that, Includes the following steps: Obtain drilling data, empirical values of fault zone height, empirical values of caving zone height, borehole depth, bottom plate position constraints, and coal mining thickness from the drilling site; Based on the Bayesian change point detection method, the change points of the drilling data are obtained; The lower limit depth of the bending zone is obtained based on the borehole depth, the empirical value of the fracture zone height, the empirical value of the collapse zone height, and the change points of the drilling data. The lower limit depth of the caving zone is determined based on the borehole depth, the bottom plate position constraints, the coal mining thickness, and the change points of the drilling data. The lower limit depth of the fracture zone is determined based on the empirical values of the fracture zone height, the empirical values of the collapse zone height, the lower limit depth of the bending zone, the lower limit depth of the collapse zone, and the change points of the drilling data.
2. The method according to claim 1, characterized in that, The changes in the drilling data include multiple changes in flushing fluid loss; obtaining the lower limit depth of the bending zone based on the borehole depth, the empirical value of the fracture zone height, the empirical value of the collapse zone height, and the changes in the drilling data includes: The first depth range is determined based on the empirical values of the fault zone height and the collapse zone height. Determine multiple depths corresponding to multiple points of change in the amount of flushing fluid leakage; The depth that falls within the first depth range and is the smallest among the multiple depths corresponding to multiple points of change in the amount of flushing fluid leakage is determined as the lower limit depth of the bending zone.
3. The method according to claim 2, characterized in that, The changes in the drilling data include multiple changes in drilling speed; determining the lower limit depth of the caving zone based on the borehole depth, the bottom plate position constraints, the coal mining thickness, and the changes in the drilling data includes: The second depth range is determined based on the borehole depth and the positional constraints of the base plate; Determine multiple depths corresponding to multiple points of change in the drilling rate; The depth that falls within the second depth range and is the smallest among the multiple depths corresponding to the multiple points of change of the drilling speed is determined as the depth of the goaf floor. The lower limit depth of the caving zone is obtained based on the depth of the goaf floor and the thickness of the coal seam.
4. The method according to claim 1, characterized in that, The lower limit depth of the fault zone is determined based on the empirical values of the fault zone height, the empirical values of the caving zone height, the lower limit depth of the curvature zone, the lower limit depth of the caving zone, and the variation points of the drilling data, including: The third depth range is determined based on the empirical values of the fracture zone height, the collapse zone height, the lower limit depth of the bend zone, and the lower limit depth of the collapse zone. When there are points of change in the drilling data within the third depth range, obtain the error value from each point of change within the third depth range to the boundary of the third depth range; The depth corresponding to the point of smallest error change is determined as the lower limit depth of the fracture zone.
5. The method according to claim 4, characterized in that, The method further includes: When there is no point of change in the drilling data within the third depth range, the midpoint of the third depth range is taken as the lower limit depth of the fracture zone.
6. The method according to claim 1, characterized in that, The Bayesian change point detection method acquires the change points of the drilling data, including: The drilling data is filtered, and the Bayesian change point detection method is used to identify change points in the filtered drilling data to obtain the change points.
7. A device for identifying the three zones of overburden in a coal mine goaf, characterized in that, include: The first acquisition module is used to acquire drilling data, empirical values of fault zone height, empirical values of caving zone height, borehole depth, bottom plate position constraints, and coal mining thickness at the drilling site. The second acquisition module is used to acquire the change points of the drilling data based on the Bayesian change point detection method; The first determining module is used to obtain the lower limit depth of the bending zone based on the borehole depth, the empirical value of the fracture zone height, the empirical value of the collapse zone height, and the change points of the drilling data; The second determining module is used to determine the lower limit depth of the caving zone based on the borehole depth, the bottom plate position constraints, the coal mining thickness, and the change points of the drilling data; The third determining module is used to determine the lower limit depth of the fracture zone based on the empirical value of the fracture zone height, the empirical value of the collapse zone height, the lower limit depth of the bending zone, the lower limit depth of the collapse zone, and the change points of the drilling data.
8. The apparatus according to claim 7, characterized in that, The second acquisition module is specifically used for: The drilling data is filtered, and the Bayesian change point detection method is used to identify change points in the filtered drilling data to obtain the change points.
9. An electronic device, characterized in that, include: processor; A memory for storing processor-executable instructions; wherein the instructions are executed by the processor to enable the processor to perform the method of any one of claims 1-6.
10. A non-transitory computer-readable storage medium, characterized in that, When the instructions in the storage medium are executed by the processor of the electronic device, the electronic device is able to perform the method of any one of claims 1-6.