Fracturing horizon detection method, device and system for thick and hard roof area of mine stratum

Abstract

The invention provides a fracturing horizon detection method, device and system for a mine stratum thick and hard roof area, and the method comprises the steps: carrying out the micro-motion signal collection of a current mine stratum area based on a plurality of detectors disposed in the current mine stratum area, obtaining time sequence signal data containing background noise of the current mine stratum area; performing Rayleigh surface wave extraction on the time sequence signal data to obtain a phase velocity frequency dispersion curve corresponding to the time sequence signal data; performing inversion according to the phase velocity dispersion curve to obtain underground shear wave velocity structure information of the current mine stratum area; and according to the underground transverse wave velocity structure information, fracturing horizon detection is carried out on the current mine stratum area, and a fracturing horizon detection result of the current mine stratum area is obtained. According to the method, the distribution change information of the rock stratum is comprehensively and accurately detected, so that the judgment precision of the fracturing horizon is improved.

Description

Technical Field

[0001] The present invention relates to the field of mine detection technology, and in particular to a method, device and system for detecting fracturing layers in a thick hard roof area of ​​a mine. Background Art

[0002] If the hard and difficult-to-collapse roof is not broken and collapsed in time, it will cause a large area of ​​suspended roof, which will cause the continuous accumulation of stress and energy, resulting in strong mining pressure on the working face, roof fall, spalling and hydraulic support damage. Therefore, accurate identification of fracturing layers is the key factor in determining the effect of fracturing pressure relief and is the most important part of fracturing design.

[0003] Among the related technologies, three-dimensional seismic exploration has high accuracy, but has a long cycle and high cost. At the engineering fracturing site, methods such as borehole core logging, column charts, borehole wall observation, and geophysical logging are usually used to determine the fracturing target layer. However, these detection methods can only provide local stratigraphic structure data and cannot fully detect the distribution change information of the rock formation. The detection range is small, which affects the accuracy of the regional rock fracturing layer determination. Summary of the invention

[0004] The present invention provides a method, device and system for detecting the fracturing layer in a thick hard roof area of ​​a mine stratum, so as to solve the problem that the prior art can only provide local stratum structure data, cannot comprehensively detect the regional stratum distribution change information, and cannot accurately perform fracturing layer detection, thereby achieving comprehensive and accurate detection of stratum distribution change information, and further improving the accuracy of fracturing layer determination.

[0005] The present invention provides a method for detecting a fracturing layer in a thick hard roof area of ​​a mine, comprising: Based on multiple geophones deployed in the current mine formation area, micro-motion signals are collected in the current mine formation area to obtain time series signal data containing background noise in the current mine formation area; Performing Rayleigh surface wave extraction on the time series signal data to obtain a phase velocity dispersion curve corresponding to the time series signal data; Inverting the phase velocity dispersion curve to obtain underground shear wave velocity structure information of the current mine formation area; According to the underground shear wave velocity structure information, the fracturing layer detection is performed on the current mine formation area to obtain the fracturing layer detection result of the current mine formation area.

[0006] According to a method for detecting a fracturing layer in a thick hard roof area of ​​a mine stratum provided by the present invention, a fracturing layer detection is performed on the current mine stratum area according to the underground shear wave velocity structure information to obtain a fracturing layer detection result of the current mine stratum area, including: Constructing a shear wave velocity profile of the current mine formation area according to the underground shear wave velocity structure information; Analyzing the shear wave velocity profile to obtain stratigraphic structure information of the current mine stratigraphic area; According to the stratum structure information, the fracturing layer detection is performed on the stratum area of ​​the current mine to obtain the fracturing layer detection result of the stratum area of ​​the current mine.

[0007] According to a method for detecting a fracturing layer in a thick hard roof area of ​​a mine stratum provided by the present invention, a shear wave velocity profile of the current mine stratum area is constructed according to the underground shear wave velocity structure information, comprising: Calibrate the shear wave velocity in the underground shear wave velocity structure information according to the well logging data of the current mine formation area to obtain the calibrated underground shear wave velocity structure information; A shear wave velocity profile of the current mine stratum area is constructed based on the calibrated underground shear wave velocity structure information.

[0008] According to a method for detecting a fracturing layer in a thick hard roof area of ​​a mine stratum provided by the present invention, the shear wave velocity profile of the current mine stratum area is constructed according to the calibrated underground shear wave velocity structure information, comprising: According to the phase velocity dispersion curve, obtaining a variation curve between the shear wave velocity in the calibrated underground shear wave velocity structure information and the depth of the current mine formation area; Preprocessing the variation curve to obtain the shear wave velocity profile; Wherein, the preprocessing includes data interpolation processing and / or smooth difference processing.

[0009] According to a method for detecting fracturing layers in a thick hard roof area of ​​a mine stratum provided by the present invention, the multiple geophones deployed in the current mine stratum area are formed by rolling deployment of multiple geophones deployed in the previous mine stratum area; Each detector deployed in the current mine formation area is coupled to the ground surface of the current mine formation area.

[0010] According to a method for detecting a fracturing layer in a thick hard roof area of ​​a mine stratum provided by the present invention, the method collects micro-vibration signals of the current mine stratum area based on multiple detectors deployed in the current mine stratum area to obtain time series signal data containing background noise of the current mine stratum area, including: Dividing the multiple geophones deployed in the current mine formation area into multiple groups according to the positions and numbers of the geophones deployed in the current mine formation area; In the current mine formation area, obtaining the measurement points corresponding to each of the groups; Based on the detectors in each group, micro-motion signals are collected at the measuring points corresponding to each group in the current mine formation area to obtain the time series signal data.

[0011] The present invention also provides a fracturing layer detection device for a thick hard roof area of ​​a mine, comprising: An acquisition module, for acquiring micro-motion signals of the current mine formation area based on a plurality of geophones deployed in the current mine formation area, and obtaining time series signal data of the current mine formation area including background noise; An extraction module, used for performing Rayleigh surface wave extraction on the time series signal data to obtain a phase velocity dispersion curve corresponding to the time series signal data; An inversion module, used to perform inversion according to the phase velocity dispersion curve to obtain underground shear wave velocity structure information of the current mine formation area; The detection module is used to perform fracturing layer detection on the current mine stratum area according to the underground shear wave velocity structure information to obtain the fracturing layer detection result of the current mine stratum area.

[0012] The present invention also provides a fracturing layer detection system for a thick hard roof area in a mine stratum, comprising a fracturing layer detection device for a thick hard roof area in a mine stratum as described in any one of the above items and a plurality of detectors; The fracturing layer detection device is communicatively connected with each of the detectors.

[0013] 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, wherein when the processor executes the program, a method for detecting fracturing layers in a thick hard roof area of ​​a mine stratum as described in any one of the above-mentioned methods is implemented.

[0014] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements a method for detecting fracturing layers in a thick hard roof area of ​​a mine stratum as described in any one of the above.

[0015] The present invention also provides a computer program product, comprising a computer program, which, when executed by a processor, implements the method for detecting fracturing layers in a thick hard roof area of ​​a mine stratum as described in any one of the above.

[0016] The method, device and system for detecting the fracturing layer in the thick hard roof area of ​​the mine stratum provided by the present invention deploy multiple detectors in the current mine stratum area to comprehensively collect the micro-motion signals in the current mine stratum area and convert them into time series signal data, thereby extracting Rayleigh surface wave information reflecting the stratum characteristics, generating a phase velocity dispersion curve, and using the phase velocity dispersion curve for inversion, thereby comprehensively and accurately revealing the underground shear wave velocity structure, and can comprehensively and accurately detect the distribution change information of the rock stratum, effectively overcoming the limitation of only being able to obtain local stratum data, and based on the obtained underground shear wave velocity structure information, can more accurately detect the fracturing layer, which not only expands the detection range, but also significantly improves the detection accuracy, thereby improving the accuracy of fracturing layer determination. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. 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 creative work.

[0018] Figure 1 It is a schematic flow chart of the method for detecting fracturing layers in a thick hard roof area of ​​a mine stratum provided by the present invention.

[0019] Figure 2 It is a distribution diagram of the stratum changes in the mine stratum area provided by the present invention.

[0020] Figure 3 It is a distribution schematic diagram of the phase velocity dispersion curve provided by the present invention.

[0021] Figure 4 It is a schematic diagram of the distribution of the variation curve between the shear wave velocity and the depth provided by the present invention.

[0022] Figure 5 It is a schematic diagram of the distribution of the shear wave velocity profile provided by the present invention.

[0023] Figure 6 It is a distribution diagram of multiple detectors provided by the present invention.

[0024] Figure 7 It is a structural schematic diagram of a fracturing layer detection device for a thick hard roof area of ​​a mine stratum provided by the present invention.

[0025] Figure 8 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0026] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings 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 of 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.

[0027] In some areas, the coal seams with hard roof account for about 30%, covering more than 50% of the mining area. If the hard and difficult-to-collapse roof is not broken or collapsed in time, it will cause a large area of ​​suspended roof, causing continuous accumulation of stress and energy, resulting in strong mine pressure on the working face, roof falls, slabs and damage to hydraulic supports. Hydraulic fracturing is a green and safe pressure relief and weakening technology. In recent years, coal mine rock fracturing technology has developed rapidly. Accurate identification of key fracturing layers is a key factor in determining the pressure relief effect of fracturing and is the most important link in fracturing design.

[0028] The logging technology in the oil industry is difficult to directly replicate downhole for multi-parameter measurement of formations. Downhole logging currently mainly obtains rock parameters such as natural gamma and resistivity, and lacks parameters of rock mechanics and strength related to hydraulic fracturing.

[0029] Therefore, currently, the main methods used to identify the target layer for fracturing are core logging, columnar graphs, borehole wall observation, geophysical logging, etc., which can only provide local stratigraphic structure and strength data, and cannot fully detect the information on rock distribution changes. There is a lack of comprehensive testing and analysis of information such as ground stress, physical and mechanical properties of rock formations, rock formation structure, rock fracturing ability, and rock distribution changes. In addition, the determination of fracturing layers still relies on experience to a certain extent, and lacks quantitative design methods, resulting in low detection accuracy and small detection range, which affects the accuracy of fracturing layer determination.

[0030] For example, core drilling is a commonly used exploration method in solid mineral geological exploration. Its main purpose is to use drilling equipment to obtain physical samples for core analysis or experimental physical mechanics and chemical experimental research to provide local formation strength data. Complementing the drilling method is ground comprehensive logging, which can reduce the workload of drilling and coring and reduce exploration costs. However, logging can only indirectly and conditionally reflect one aspect of the geological characteristics of the rock formation, and cannot fully understand the underground geological features and top plate geological information. Therefore, in general, core drilling and comprehensive logging require drilling. Although the single well survey accuracy is high, its construction cost is high, the exploration cycle is long, and the coverage area is limited. It can only provide local formation strength data, and cannot fully detect the distribution change information of the rock formation. The detection accuracy is low and the detection range is small, which affects the accuracy of fracturing layer determination.

[0031] For another example, the selection of regional fracturing layers is mainly based on the mine column chart. The constraint range of a single borehole is insufficient, and the control range is above 300-500m. Due to its age and the fact that it was provided by a third party, the column chart data is questionable, and the control range does not match the accuracy required for fracturing construction. Even for two boreholes 300m apart, the rock properties of the top and bottom plates are quite different, and it is difficult to draw a comparison map of the coal and rock layers, which has a great impact on the construction of the geological modeling system. The test data of a single borehole can only provide local formation strength data, and cannot fully detect the distribution change information of the rock formation. The detection accuracy is low and the detection range is small, which affects the accuracy of the fracturing layer determination.

[0032] In summary, the fracturing layer detection methods used in related technologies are mostly local testing methods. The test data of a single borehole can only provide local formation strength data, and cannot fully detect the distribution change information of the rock formation. The detection accuracy is low and the detection range is small. It is impossible to quickly detect and model the rock formation in the area to be fractured to achieve accurate fracturing layer determination, thus affecting the accuracy of fracturing layer determination.

[0033] In response to the above problems, the present embodiment provides a method for detecting fracturing layers in a thick hard roof area of ​​a mine stratum. The method provided by the method performs detection imaging through multiple detectors deployed in the current mine stratum area to collect micro-motion signals in the current mine stratum area, and extracts the dispersion curve of the surface wave (such as Rayleigh surface wave) from the micro-motion signal based on the theory of stationary random processes. By inverting the dispersion curve, more comprehensive velocity structure information of the underground medium is obtained, and based on the velocity structure information, regional strata detection is quickly achieved, so as to finely and comprehensively characterize the rock stratum distribution change information and stratum structure information in the current mine stratum area. This not only avoids large-scale drilling coring and logging, and greatly reduces the detection cost, but also can quickly achieve regional detection and three-dimensional velocity modeling, thereby accurately determining the key fracturing layers.

[0034] Figure 1 Schematic diagram of the process of the method for detecting the fracturing layer in the thick hard roof area of ​​the mine stratum provided by the present invention. Figure 1 As shown, the method includes step 110 , step 120 , step 130 and step 140 .

[0035] Step 110, based on multiple detectors deployed in the current mine formation area, collect micro-motion signals in the current mine formation area to obtain time series signal data containing background noise in the current mine formation area.

[0036] Figure 2 : is a distribution diagram of the stratum changes in the mine stratum area provided by the present invention. Among them, the current mine stratum area is the mine stratum area where the fracturing layer detection is required in the current cycle. Figure 2As shown, the current mine stratum area may be an area where the stratum changes drastically, or an area where the reservoir physical properties are poor, etc., and this embodiment does not specifically limit this.

[0037] Each cycle may be divided according to a preset time period, such as 30 minutes or 1 hour as one cycle, etc., which is not specifically limited in this embodiment.

[0038] There are multiple detectors deployed in the current mine formation area. Each detector is fully coupled with the ground in the current mine formation area. The tail vertebra is inserted as solidly, stably, correctly, deeply and straightly as possible. Each detector is distributed at different measuring points in the current mine formation area, and is used to synchronously collect micro-motion signals at multiple measuring points to eliminate signal inconsistency problems caused by time differences.

[0039] Here, the specific number of detectors can be configured according to actual needs, such as 50. These detectors are arranged according to a preset geometric layout to ensure that the target detection area can be fully covered and the micro-motion signals from deep in the formation can be effectively received. Among them, the channel spacing between the detectors and the geometric layout arrangement can be configured according to actual needs. For example, the channel spacing is set to 5 meters, and the geometric layout arrangement is that the detectors for detection at each measuring point form a straight-line array, which helps to enhance the signal reception capability and reduce the signal interference caused by the multipath effect, thereby simplifying the data processing process and improving the data quality.

[0040] Optionally, when conducting fracturing layer detection in the current mine formation area, multiple detectors deployed in the current mine formation area can be used to collect micro-motion signals in the current mine formation area to obtain the micro-motion signals of each measuring point in the current mine formation area at different sampling times, thereby forming time series signal data based on the micro-motion signals of each measuring point at different sampling times.

[0041] Here, the method of forming the timing signal data includes directly constructing the timing signal data based on the micro-motion signal of each measuring point at different sampling times; or constructing the timing signal data after filtering, denoising and other pre-processing of the micro-motion signal of each measuring point at different sampling times. This embodiment does not make specific limitations on this.

[0042] Step 120: extract Rayleigh surface waves from the time series signal data to obtain a phase velocity dispersion curve corresponding to the time series signal data.

[0043] Figure 3 It is a distribution diagram of the phase velocity dispersion curve provided by the present invention. The horizontal coordinate of the phase velocity dispersion curve is the shear wave velocity value (or Vel), and the vertical coordinate is the frequency (or Freq).

[0044] Optionally, after acquiring the time series signal data, the time series signal data may be processed in the transform domain to extract the dispersion characteristics of the Rayleigh surface wave, and then obtain the phase velocity dispersion curve corresponding to the time series signal data. The specific distribution form is as follows: Figure 3 shown.

[0045] The so-called transform domain processing can be implemented based on algorithms such as spatial autocorrelation method, frequency-wavenumber method, etc., and can be adaptively selected according to factors such as geological conditions, exploration requirements and data quality in the current mine formation area. This embodiment does not make specific limitations on this.

[0046] Among them, the spatial autocorrelation method extracts the dispersion characteristics of Rayleigh surface waves by analyzing the similarity (i.e., autocorrelation) of micromotion signals between different receiving points; specifically, it can be based on time series signal data, calculate the autocorrelation coefficient of signals between different measuring points, and draw an autocorrelation diagram of the autocorrelation coefficient and the distance between measuring points. Finally, the phase velocity dispersion curve of Rayleigh surface waves is extracted from the autocorrelation diagram through fitting or inversion methods.

[0047] The frequency-wavenumber method extracts the dispersion characteristics of Rayleigh surface waves by transforming the micromotion signal from the time-space domain to the frequency-wavenumber domain; specifically, it can be a transformation method such as two-dimensional Fourier transform on the time series signal data to transform it from the time-space domain to the frequency-wavenumber domain. In the frequency-wavenumber domain, the phase velocity dispersion curve is extracted by identifying the energy distribution characteristics of the Rayleigh surface wave.

[0048] Step 130, performing inversion according to the phase velocity dispersion curve to obtain underground shear wave velocity structure information of the current mine formation area.

[0049] Optionally, after obtaining the phase velocity dispersion curve, the phase velocity dispersion curve can be inverted to obtain underground shear wave velocity structure information of the current mine formation area. The underground shear wave velocity structure information is used to characterize the shear wave velocity values ​​at different depths and locations underground.

[0050] The implementation method of inversion here includes: establishing an initial geological model based on known geological information of the mine stratigraphic area, such as drilling data, geological maps, etc. The model should include basic parameters such as the thickness, lithology, and density of the stratigraphic area. According to the characteristics of the dispersion curve and the complexity of the geological model, a suitable inversion algorithm is selected, such as the damped least squares method of singular value decomposition, genetic algorithm, simulated annealing algorithm, etc., and forward simulation is performed using the initial geological model and the inversion algorithm to calculate the phase velocity dispersion inference curve, and the phase velocity dispersion inference curve is compared with the phase velocity dispersion curve obtained by detection to adjust the parameters of the initial geological model so that the phase velocity dispersion inference curve gradually approaches the phase velocity dispersion curve obtained by detection, thereby obtaining the final geological model, and then extracting the underground shear wave velocity structure information of the current mine stratigraphic area from the final address model.

[0051] Step 140: perform fracturing layer detection on the current mine formation area according to the underground shear wave velocity structure information to obtain a fracturing layer detection result of the current mine formation area.

[0052] Optionally, after obtaining the underground shear wave velocity structure information, the underground shear wave velocity structure information can be analyzed to obtain the stratigraphic structure information of the current mine stratigraphic area, thereby optimizing the detection of thick, hard and complete strata as the key layer for fracturing (or fracturing layer) in the current mine stratigraphic area based on the stratigraphic structure information, thereby obtaining the fracturing layer detection result of the current mine stratigraphic area. The fracturing layer detection result here includes identification information indicating whether there is a key layer for fracturing in the current mine stratigraphic area, and the location information of the key layer for fracturing.

[0053] The method provided in this embodiment deploys multiple detectors in the current mine formation area to comprehensively collect micro-vibration signals in the current mine formation area and convert them into time series signal data, thereby extracting Rayleigh surface wave information reflecting the formation characteristics, generating a phase velocity dispersion curve, and using the phase velocity dispersion curve for inversion, thereby comprehensively and accurately revealing the underground shear wave velocity structure, which can comprehensively and accurately detect the distribution change information of the rock formation, effectively overcome the limitation of only being able to obtain local formation data, and based on the obtained underground shear wave velocity structure information, it can more accurately detect the fracturing layer, which not only expands the detection range, but also significantly improves the detection accuracy, thereby improving the accuracy of fracturing layer determination.

[0054] In some embodiments, the step 140 of performing fracturing layer detection on the current mine formation area according to the underground shear wave velocity structure information to obtain the fracturing layer detection result of the current mine formation area further includes: Constructing a shear wave velocity profile of the current mine formation area according to the underground shear wave velocity structure information; Analyzing the shear wave velocity profile to obtain stratigraphic structure information of the current mine stratigraphic area; According to the stratum structure information, the fracturing layer detection is performed on the stratum area of ​​the current mine to obtain the fracturing layer detection result of the stratum area of ​​the current mine.

[0055] Optionally, a shear wave velocity profile of the current mine stratigraphic area is constructed directly based on the underground shear wave velocity structure information; or the underground shear wave velocity structure information is pre-processed, such as calibration, outlier verification, etc., and then a shear wave velocity profile of the current mine stratigraphic area is constructed based on the processed underground shear wave velocity structure information. This embodiment does not make any specific limitations on this.

[0056] Subsequently, according to the velocity change characteristics on the shear wave velocity profile, the underground strata are divided into different velocity layers, and on the basis of the velocity layer division, the interfaces between the strata are further identified to obtain the distribution range, thickness and morphology of the strata, as well as the lithology, physical properties, strength, stability and other stratum structural information. Based on the stratum structural information, the thick, hard and complete strata are optimized for detection in the current mine stratum area as the key layer (or fracturing layer) for fracturing, thereby obtaining the fracturing layer detection results in the current mine stratum area.

[0057] The method provided in this embodiment can quickly realize the velocity structure detection of the regional strata of the mine roof through the inversion of the dispersion curve, so as to help accurately determine the stratigraphic structure information of the current mine stratigraphic area, thereby accurately determining the thick, hard and complete roof strata above the coal seam based on the stratigraphic structure information of the current mine stratigraphic area, and formulating a fracturing plan, avoiding large-scale drilling coring and logging, and greatly reducing the detection cost.

[0058] In some embodiments, constructing the shear wave velocity profile of the current mine formation area according to the underground shear wave velocity structure information includes: Calibrate the shear wave velocity in the underground shear wave velocity structure information according to the well logging data of the current mine formation area to obtain the calibrated underground shear wave velocity structure information; A shear wave velocity profile of the current mine formation area is constructed based on the calibrated underground shear wave velocity structure information.

[0059] Among them, the magnitude and change of shear wave velocity can reflect the continuity, integrity, compressibility, etc. of the formation. Although the apparent S wave velocity (or shear wave velocity) obtained by the phase velocity dispersion curve is not equal to the actual apparent S wave velocity value, it has a physical parameter with velocity dimension, and its change trend can reflect the actual S wave velocity change, so geological diagnosis can be performed through the profile of the apparent S wave velocity. Therefore, in order to improve the accuracy of fracturing layer detection, it is also necessary to calibrate the apparent S wave velocity obtained by the phase velocity dispersion curve in combination with the logging data of the current mine formation area to eliminate the error or deviation of the system measurement, thereby obtaining the absolute apparent S wave velocity of the formation. The calibrated underground shear wave velocity data is projected onto the profile frame, and the shear wave velocity change curves at different depths are drawn on the profile according to the depth and location information of the data, thereby obtaining the shear wave velocity profile.

[0060] In some embodiments, constructing the shear wave velocity profile of the current mine formation area according to the calibrated underground shear wave velocity structure information includes: According to the phase velocity dispersion curve, obtaining a variation curve between the shear wave velocity in the calibrated underground shear wave velocity structure information and the depth of the current mine formation area; Preprocessing the variation curve to obtain the shear wave velocity profile; Wherein, the preprocessing includes data interpolation processing and / or smooth difference processing.

[0061] Figure 4 It is a schematic diagram of the distribution of the variation curve between the shear wave velocity and the depth provided by the present invention.

[0062] Figure 5 It is a schematic diagram of the distribution of the shear wave velocity profile provided by the present invention.

[0063] First, select the key frequency point representing the underground shear wave velocity from the phase velocity dispersion curve, or determine a frequency range according to the analysis needs. Then, read the corresponding shear wave velocity value at this frequency point or frequency range, and match these velocities with their corresponding depth information in the underground strata. Subsequently, using these matched shear wave velocity and depth data, a curve of shear wave velocity changing with depth is generated through drawing or data processing software, as shown in the following figure. Figure 4 As shown in the figure, the variation curve between the shear wave velocity and the depth data intuitively shows the law of the underground shear wave velocity changing with depth in the current mine formation area.

[0064] Subsequently, data interpolation processing is performed on the acquired variation curve between shear wave velocity and depth to fill the blank areas caused by missing data or uneven measurement intervals to ensure the continuity of the curve, and / or smooth difference processing is applied to reduce the small fluctuations and noise on the curve through averaging or filtering technology to make the curve smoother and easier to analyze, thereby obtaining a clear and continuous shear wave velocity profile, as shown in the following example. Figure 5 As shown. Based on the shear wave velocity profile, the variation law of underground shear wave velocity with depth can be accurately reflected.

[0065] The method provided in this embodiment obtains a curve of calibrated underground shear wave velocity changing with the depth of the current mine formation area based on the phase velocity dispersion curve, and performs preprocessing steps including data interpolation and / or smooth difference on the curve, thereby ensuring the integrity and smoothness of the data, thereby drawing a clear and accurate shear wave velocity profile, and further realizing the detection of fracturing layers in the mine formation area more conveniently and accurately.

[0066] In some embodiments, the multiple geophones deployed in the current mine formation area are formed by rolling deployment of multiple geophones deployed in the previous mine formation area; The multiple geophones deployed in the current mine formation area are arranged in a straight line array at preset intervals; Each detector deployed in the current mine formation area is coupled to the ground surface of the current mine formation area.

[0067] Optionally, when the number of detectors is valid, the detection of fracturing layers in different mine formation areas can be achieved through rolling detection. Taking the last detection as an example, the deployment of multiple detectors in the current mine formation area can be achieved by using a preset track spacing and a single measuring point with a preset number of detectors to form a straight-line array, so that after the last detection is completed, the multiple detectors are rolled back as a whole by a preset distance, and each detector is kept coupled with the ground in the current mine formation area, thereby achieving the deployment of multiple detectors in the current mine formation area. The preset track spacing, preset number, and preset rolling distance here can all be adaptively set according to actual needs.

[0068] It should be noted that before the detector collects data, it is necessary to number the detectors in ascending order, and record the detector numbers and their positions in the array accurately and in detail. When collecting data, it is necessary to ensure that the detectors are fully coupled to the ground, and the tail vertebrae are inserted as firmly, stably, correctly, deeply, and straightly as possible.

[0069] The following describes the deployment of multiple detectors using a specific example.

[0070] Figure 6FIG. 1 is a schematic diagram of the distribution of multiple detectors provided by the present invention. Figure 6 As shown in the figure, 50 geophones are used to collect data simultaneously, the track spacing between each geophone is 5 meters, and a single measuring point uses 31 geophones to form a straight array for data collection, so the observation radius is 75 meters. Accordingly, a single arrangement can measure 20 measuring points. After the collection of one arrangement of data (that is, the micro-seismic signal of a mine formation area) is completed, the interval distance of 20 geophones is rolled back to continue the next micro-seismic signal collection of the mine formation area.

[0071] The method provided in this embodiment, through a rolling deployment strategy, multiple detectors deployed in the current mine formation area can continuously and efficiently cover and detect a wider underground area, significantly improving the flexibility and efficiency of data acquisition. At the same time, ensuring that each detector is tightly coupled with the ground in the current mine formation area helps to accurately capture the weak signals of crustal micromotion, improve the accuracy and reliability of the data, and thus improve the accuracy and convenience of fracturing layer detection.

[0072] In some embodiments, the micro-motion signal acquisition of the current mine formation area is performed based on multiple geophones deployed in the current mine formation area to obtain time series signal data containing background noise of the current mine formation area, including: Dividing the multiple geophones deployed in the current mine formation area into multiple groups according to the positions and numbers of the geophones deployed in the current mine formation area; In the current mine formation area, obtaining the measurement points corresponding to each of the groups; Based on the detectors in each group, micro-motion signals are collected at the measuring points corresponding to each group in the current mine formation area to obtain the time series signal data.

[0073] Optionally, before the geophones collect data, the geophones are numbered in ascending order, and the geophone numbers and their positions in the array are accurately and detailedly recorded. Therefore, when collecting data, the positions and numbers of the geophones deployed in the current mine formation area can be obtained based on this, and the multiple geophones deployed in the current mine formation area can be divided into multiple groups according to the positions and numbers.

[0074] Subsequently, in the current mine formation area, the measuring points corresponding to each group are associated and acquired, and multiple micro-motion signal acquisitions are performed at the measuring points corresponding to each group in the current mine formation area using the detectors in each group, thereby constructing time series signal data.

[0075] For example, Figure 6 As shown, the 15 detectors arranged in the front form a group, and they jointly perform micro-motion signal detection at measuring point 1.

[0076] The method provided in this embodiment accurately divides the detectors in the current mine formation area into multiple groups, and determines the corresponding measuring points based on the positions of the detectors in each group, so as to conduct targeted micro-motion signal collection. It can efficiently and accurately obtain time series signal data reflecting the formation characteristics, and significantly improve the accuracy and efficiency of fracturing layer detection.

[0077] The following is a description of the fracturing layer detection device for the thick hard roof area of ​​the mine strata provided by the present invention. The fracturing layer detection device for the thick hard roof area of ​​the mine strata described below and the fracturing layer detection method for the thick hard roof area of ​​the mine strata described above can be referenced to each other.

[0078] Figure 7 Schematic diagram of the structure of the fracturing layer detection device for the thick hard roof area of ​​the mine stratum provided by the present invention; Figure 7 As shown, the device includes a collection module 710 , an extraction module 720 , an inversion module 730 and a detection module 740 .

[0079] The acquisition module 710 is used to acquire micro-motion signals in the current mine formation area based on multiple detectors deployed in the current mine formation area, and obtain time series signal data containing background noise in the current mine formation area; The extraction module 720 is used to extract Rayleigh surface waves from the time series signal data to obtain a phase velocity dispersion curve corresponding to the time series signal data; The inversion module 730 is used to perform inversion according to the phase velocity dispersion curve to obtain underground shear wave velocity structure information of the current mine formation area; The detection module 740 is used to perform fracturing layer detection on the current mine formation area according to the underground shear wave velocity structure information, and obtain the fracturing layer detection result of the current mine formation area.

[0080] The device provided in this embodiment deploys multiple detectors in the current mine formation area to comprehensively collect micro-motion signals in the current mine formation area and convert them into time series signal data, thereby extracting Rayleigh surface wave information reflecting the formation characteristics, generating a phase velocity dispersion curve, and using the phase velocity dispersion curve for inversion, thereby comprehensively and accurately revealing the underground shear wave velocity structure, and can comprehensively and accurately detect the distribution change information of the rock formation, effectively overcoming the limitation of only being able to obtain local formation data, and based on the obtained underground shear wave velocity structure information, it can more accurately detect the fracturing layer, which not only expands the detection range, but also significantly improves the detection accuracy, thereby improving the accuracy of fracturing layer determination.

[0081] In some embodiments, the present embodiment also provides a fracturing layer detection system for a thick hard roof area in a mine stratum, characterized in that it comprises a fracturing layer detection device for a thick hard roof area in a mine stratum and multiple detectors provided in the above embodiments; the fracturing layer detection device is communicatively connected with each of the detectors, and the fracturing layer detection device is used to interact with the detectors to implement the fracturing layer detection method for a thick hard roof area in a mine stratum provided in the above embodiments. It will not be repeated here, and please refer to the contents described in the above embodiments for details.

[0082] The system provided in this embodiment deploys multiple detectors in the current mine formation area to comprehensively collect micro-vibration signals in the current mine formation area and convert them into time series signal data, thereby extracting Rayleigh surface wave information reflecting the formation characteristics, generating a phase velocity dispersion curve, and using the phase velocity dispersion curve for inversion, thereby comprehensively and accurately revealing the underground shear wave velocity structure, and can comprehensively and accurately detect the distribution change information of the rock formation, effectively overcoming the limitation of only being able to obtain local formation data, and based on the obtained underground shear wave velocity structure information, it can more accurately detect the fracturing layer, which not only expands the detection range, but also significantly improves the detection accuracy, thereby improving the accuracy of fracturing layer determination.

[0083] Figure 8 An example of a physical structure diagram of an electronic device is shown in FIG. Figure 8 As shown, the electronic device may include: a processor 810, a communication interface 820, a memory 830 and a communication bus 840, wherein the processor 810, the communication interface 820 and the memory 830 communicate with each other through the communication bus 840. The processor 810 may call the logic instructions in the memory 830 to execute a method for detecting a fracturing layer in a thick hard roof area of ​​a mine stratum, the method comprising: based on a plurality of detectors deployed in the current mine stratum area, collecting micro-seismic signals in the current mine stratum area to obtain time series signal data containing background noise in the current mine stratum area; extracting Rayleigh surface waves from the time series signal data to obtain a phase velocity dispersion curve corresponding to the time series signal data; inverting according to the phase velocity dispersion curve to obtain underground shear wave velocity structure information of the current mine stratum area; and detecting a fracturing layer in the current mine stratum area according to the underground shear wave velocity structure information to obtain a fracturing layer detection result in the current mine stratum area.

[0084] In addition, the logic instructions in the above-mentioned memory 830 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0085] On the other hand, the present invention also provides a computer program product, which includes a computer program. The computer program can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the method for detecting the fracturing layer in the thick hard roof area of ​​the mine stratum provided by the above-mentioned methods. The method includes: based on multiple detectors deployed in the current mine stratum area, collecting micro-motion signals in the current mine stratum area to obtain time series signal data containing background noise in the current mine stratum area; extracting Rayleigh surface waves from the time series signal data to obtain a phase velocity dispersion curve corresponding to the time series signal data; inverting according to the phase velocity dispersion curve to obtain underground shear wave velocity structure information of the current mine stratum area; based on the underground shear wave velocity structure information, performing fracturing layer detection on the current mine stratum area to obtain a fracturing layer detection result of the current mine stratum area.

[0086] On the other hand, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, is implemented to execute the method for detecting a fracturing layer in a thick hard roof area of ​​a mine stratum provided by the above-mentioned methods, the method comprising: based on a plurality of detectors deployed in the current mine stratum area, collecting micro-motion signals in the current mine stratum area to obtain time series signal data containing background noise in the current mine stratum area; extracting Rayleigh surface waves from the time series signal data to obtain a phase velocity dispersion curve corresponding to the time series signal data; inverting according to the phase velocity dispersion curve to obtain underground shear wave velocity structure information of the current mine stratum area; and performing fracturing layer detection on the current mine stratum area based on the underground shear wave velocity structure information to obtain a fracturing layer detection result of the current mine stratum area.

[0087] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0088] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0089] 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 make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for detecting fracturing layers in a thick hard roof area of ​​a mine, characterized in that: include: Based on multiple geophones deployed in the current mine formation area, micro-motion signals are collected in the current mine formation area to obtain time series signal data containing background noise in the current mine formation area; Performing Rayleigh surface wave extraction on the time series signal data to obtain a phase velocity dispersion curve corresponding to the time series signal data; Inverting the phase velocity dispersion curve to obtain underground shear wave velocity structure information of the current mine formation area; According to the underground shear wave velocity structure information, the fracturing layer detection is performed on the current mine formation area to obtain the fracturing layer detection result of the current mine formation area.

2. The method for detecting the fracturing layer in the thick hard roof area of ​​a mine according to claim 1, characterized in that: According to the underground shear wave velocity structure information, the fracturing layer detection is performed on the current mine formation area to obtain the fracturing layer detection result of the current mine formation area, including: Constructing a shear wave velocity profile of the current mine formation area according to the underground shear wave velocity structure information; Analyzing the shear wave velocity profile to obtain stratigraphic structure information of the current mine stratigraphic area; According to the stratum structure information, the fracturing layer detection is performed on the stratum area of ​​the current mine to obtain the fracturing layer detection result of the stratum area of ​​the current mine.

3. The method for detecting the fracturing layer in the thick hard roof area of ​​a mine according to claim 2, characterized in that: According to the underground shear wave velocity structure information, a shear wave velocity profile of the current mine formation area is constructed, including: Calibrate the shear wave velocity in the underground shear wave velocity structure information according to the well logging data of the current mine formation area to obtain the calibrated underground shear wave velocity structure information; A shear wave velocity profile of the current mine stratum area is constructed based on the calibrated underground shear wave velocity structure information.

4. The method for detecting the fracturing layer in the thick hard roof area of ​​a mine according to claim 3, characterized in that: The step of constructing a shear wave velocity profile of the current mine formation area according to the calibrated underground shear wave velocity structure information includes: According to the phase velocity dispersion curve, obtaining a variation curve between the shear wave velocity in the calibrated underground shear wave velocity structure information and the depth of the current mine formation area; Preprocessing the variation curve to obtain the shear wave velocity profile; Wherein, the preprocessing includes data interpolation processing and / or smooth difference processing.

5. The method for detecting the fracturing layer in the thick hard roof area of ​​a mine according to any one of claims 1 to 4, characterized in that: The multiple geophones deployed in the current mine formation area are formed by rolling deployment of multiple geophones deployed in the previous mine formation area; Each detector deployed in the current mine formation area is coupled to the ground surface of the current mine formation area.

6. The method for detecting fracturing layers in a thick hard roof area of ​​a mine according to any one of claims 1 to 4, characterized in that: The method of collecting micro-motion signals of the current mine stratum region based on multiple geophones deployed in the current mine stratum region to obtain time series signal data containing background noise of the current mine stratum region includes: Dividing the multiple geophones deployed in the current mine formation area into multiple groups according to the positions and numbers of the geophones deployed in the current mine formation area; In the current mine formation area, obtaining the measurement points corresponding to each of the groups; Based on the detectors in each group, micro-motion signals are collected at the measuring points corresponding to each group in the current mine formation area to obtain the time series signal data.

7. A fracturing layer detection device for a thick hard roof area of ​​a mine, characterized in that: include: An acquisition module, for acquiring micro-motion signals of the current mine stratum region based on a plurality of geophones deployed in the current mine stratum region, and obtaining time series signal data of the current mine stratum region including background noise; An extraction module, used for performing Rayleigh surface wave extraction on the time series signal data to obtain a phase velocity dispersion curve corresponding to the time series signal data; An inversion module, used to perform inversion according to the phase velocity dispersion curve to obtain underground shear wave velocity structure information of the current mine formation area; The detection module is used to perform fracturing layer detection on the current mine stratum area according to the underground shear wave velocity structure information to obtain the fracturing layer detection result of the current mine stratum area.

8. A fracturing layer detection system for thick hard roof areas in mine strata, characterized in that: It comprises a fracturing layer detection device for a thick hard roof area of ​​a mine stratum as claimed in claim 7 and a plurality of geophones; The fracturing layer detection device is communicatively connected with each of the detectors.

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 fracturing layer in the thick hard roof area of ​​the mine stratum as described in any one of claims 1 to 6 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 fracturing layers in a thick hard roof area of ​​a mine stratum as described in any one of claims 1 to 6 is implemented.

Images