CN118464128B - Water pressure vibration multi-field coupling water channel detection method, terminal and storage medium - Google Patents

Abstract

The application is suitable for the technical field of coal mine water control and provides a water channel detection method, a terminal and a storage medium for hydraulic shock multi-field coupling, wherein the method comprises the following steps: acquiring microseismic data, supporting pressure of a working surface and water inflow of the working surface; determining microseismic response characteristics based on the microseismic data; determining a migration distribution area of the pressure based on the microseismic response characteristics and the bearing pressure of the working surface; determining a position area of the water passing channel based on the water inflow of the working surface, the pressure migration distribution area and the microseismic response characteristics; the position area of the water channel is used as the main content of the early warning signal. The application can monitor the cracking condition of the surrounding rock of the top and bottom plates in real time in stoping, analyze the water passing channel of the bottom plate of the coal bed, and carry out the range delineation of the water passing channel so as to carry out more accurate and effective early warning.

Description

Water pressure and seismic multi-field coupling water channel detection method, terminal and storage medium

Technical Field

The present application belongs to the technical field of coal mine water prevention and control, and in particular, relates to a water pressure shock multi-field coupled water channel detection method, terminal and storage medium.

Background Art

Coal resources are the main and irreplaceable energy source in my country for a long period of time in the future. Mine water damage is one of the main common disasters in coal mines. There are five types of mine water damage, including surface water damage, pore water damage, fissure water damage, karst water damage and old empty (kiln) water damage. Except for old empty water, the other four types are closely related to the hydrogeological conditions of coal mines.

Most of the major water inrush disasters in my country's coal mines are floor water, especially in pressure mining areas, where the hydrogeological conditions are relatively complex. With the increase in mining depth in recent years, the threat of high-pressure water inrush from Ordovician or Cambrian limestone has become increasingly serious. Even after implementing geophysical exploration, drilling, regional governance, drainage and pressure reduction measures, the risk of floor water inrush during mining is still very high.

Therefore, it is necessary to monitor the fracture of the roof and floor surrounding rocks in real time during mining and analyze the water passage of the coal seam floor.

Summary of the invention

In order to overcome the problems existing in the related technology, the embodiment of the present application provides a water pressure and seismic multi-field coupled water channel detection method, terminal and storage medium, which can monitor the fracture of the roof and floor surrounding rocks in real time during mining, analyze the water channel of the coal seam floor, delineate the range of the water channel, and provide more accurate and effective early warning.

This application is implemented through the following technical solutions:

In a first aspect, the present application provides a water pressure-vibration multi-field coupling water passage detection method, comprising:

Obtain microseismic data, bearing pressure of the working face and water inflow of the working face;

Determine microseismic response characteristics based on microseismic data;

Determine the pressure migration distribution area based on the microseismic response characteristics and the support pressure of the working face;

The location area of the water passage is determined based on the water inflow of the working face, the distribution area of the pressure migration and the microseismic response characteristics; the location area of the water passage serves as the main content of the early warning signal.

In a possible implementation manner of the first aspect, the microseismic response characteristics include the occurrence location, frequency, and energy of the microseismic event;

Based on microseismic data, determine the microseismic response characteristics, including:

Perform quality control processing on microseismic data to obtain microseismic attribute data; microseismic attribute data includes the frequency and energy of microseismic events;

Acquisition of mining engineering data;

Construct a three-dimensional geological model based on microseismic attribute data and mining engineering data;

Based on the 3D geological model, the location of microseismic events is identified.

In a possible implementation of the first aspect, the microseismic events include top plate events and bottom plate events;

Based on the 3D geological model, the location of microseismic events is identified, including:

Based on the three-dimensional geological model, the layers of the top plate events are divided according to the upper three zones and lithology, and the layers of the bottom plate events are divided according to the development of the aquifer;

Based on the interval of the roof event, determine the location of the roof event;

Based on the layer segment of the bottom plate event, the occurrence location of the bottom plate event is determined.

In a possible implementation of the first aspect, determining the pressure migration distribution area based on the microseismic response characteristics and the support pressure of the working surface includes:

Based on the microseismic response characteristics, set the pressure determination index;

When the pressure judgment index is met, a pressure curve chart is drawn based on the support pressure of the working surface;

Based on the pressure determination index and the pressure curve, the pressure situation of the working face is verified according to the mining direction, and the pressure migration distribution area is determined.

In a possible implementation manner of the first aspect, the position area of the water passage includes a plane position area of the water passage and a vertical position area of the water passage;

Based on the water inflow of the working face, the distribution area of the pressure migration and the microseismic response characteristics, the location area of the water channel is determined, including:

Based on the microseismic response characteristics, determine the plane cloud map of the warning layer and the warning layer with the density of deep events and the vertical cloud map of the warning layer and the warning layer with the density of deep events;

Based on the distribution area of pressure migration and the plane cloud map of the early warning layer and the early warning layer and the deep event density, the plane location area of the water channel is determined;

After determining the planar location area of the water passage, the vertical location area of the water passage is determined based on the water inflow of the working face and the vertical cloud map of the early warning layer and the deep event density of the early warning layer.

In a possible implementation of the first aspect, determining the plane location area of the water passage based on the pressure migration distribution area and the early warning layer and the plane cloud map of the early warning layer and the deep event density includes:

Based on the pressure migration distribution area and the density of the plane cloud map of the early warning layer and the early warning layer with the density of the deep event, the initial plane position area of the water channel is determined;

Based on the maximum density range of the early warning layer and the plane cloud map of the early warning layer with deep event density, the plane position area of the water passage is determined from the initial plane position area.

In a possible implementation manner of the first aspect, after determining the initial plane position area of the water passage, the water pressure-vibration multi-field coupling water passage detection method further includes:

When the total frequency increase of microseismic events obtained from the plane cloud map of the early warning layer and the deep event density above the early warning layer exceeds the second preset ratio, the occurrence of microseismic events in the early warning layer continues for a preset time, and the total increase rate of the rear response area and the advance response area in the microseismic monitoring data reaches a preset value, it is determined that the time after this moment is the special attention time for microseismic events;

Accordingly, based on the maximum density range of the plane cloud map of the early warning layer and the early warning layer with deep event density, the plane location area of the water passage is determined from the initial plane location area, including:

During the special attention time, based on the maximum density range of the plane cloud map of the early warning layer and the early warning layer with deep event density, the plane location area of the water passage is determined from the initial plane location area.

In a possible implementation of the first aspect, based on the water inflow of the working face and the warning layer and the vertical cloud map of the warning layer with deep event density, determining the vertical position area of the water passage includes:

In the plane position area of the water channel, based on the maximum density of the vertical cloud map of the early warning layer and the deep event density of the early warning layer, the development position of the water channel at the depth of the coal seam floor is determined;

When the water flow rate of the working face reaches the lowest point of the water flow change curve at a location where the water flow channel is developed deep in the coal seam floor, the vertical position area entering the water flow channel at this time is determined.

In a second aspect, an embodiment of the present application provides a terminal including a memory and a processor, wherein the memory stores a computer program that can be run on the processor, and when the processor executes the computer program, the water pressure-vibration multi-field coupling water channel detection method as in the first aspect is implemented.

In a third aspect, an embodiment of the present application provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, it implements the water pressure-vibration multi-field coupling water channel detection method as in the first aspect.

Compared with the related art, the embodiments of the present application have the following beneficial effects:

The embodiment of the present application can monitor the fracture of the roof and floor surrounding rocks in real time during mining by acquiring microseismic data, water level and volume changes, and the law of mine pressure manifestation, and can accurately analyze the situation of the water channel on the coal seam floor, and accurately predict and determine the location of the water channel, and can identify the possible location area of the water channel in advance, so as to conduct risk assessment and early warning in time. This is of great significance for preventing water disasters and ensuring engineering safety.

The beneficial effects of the second to third aspects mentioned above refer to the beneficial effects of the water pressure-vibration multi-field coupling water passage detection method of the first aspect, and will not be repeated here.

It is to be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the present specification.

BRIEF DESCRIPTION OF THE DRAWINGS

In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the embodiments or related technical descriptions will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

FIG1 is a schematic flow chart of a water pressure-vibration multi-field coupling water passage detection method provided in an embodiment of the present application;

FIG2 is a schematic diagram of the distribution of detectors and monitoring substations provided in an embodiment of the present application;

FIG3 is a schematic diagram of the frequency increase of microseismic events in the R1 layer segment and the change of the sum of daily energy of the top plate provided in an embodiment of the present application;

FIG4 is a schematic diagram of a pressure migration distribution area provided by an embodiment of the present application;

FIG5 is a diagram showing the total frequency of floor events according to an embodiment of the present application;

FIG6 is a plane cloud diagram of the early warning layer and the early warning layer with deep event density provided by an embodiment of the present application;

FIG. 7 is a vertical cloud diagram of the early warning layer and the early warning layer with deep event density provided by an embodiment of the present application;

FIG8 is a schematic diagram of a planar position area of a water passage provided by an embodiment of the present application;

FIG9 is a water volume change curve diagram of water inflow of a working face provided by an embodiment of the present application;

FIG10 is a model of a vertical position area of a water channel provided in one embodiment of the present application;

FIG11 is a schematic structural diagram of a water pressure-vibration multi-field coupling water passage detection device provided in one embodiment of the present application;

FIG. 12 is a schematic diagram of the structure of a terminal provided in an embodiment of the present application.

DETAILED DESCRIPTION

In the following description, specific details such as specific system structures, technologies, etc. are provided for the purpose of illustration rather than limitation, so as to provide a thorough understanding of the embodiments of the present application. However, it should be clear to those skilled in the art that the present application may also be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, devices, circuits, and methods are omitted to prevent unnecessary details from obstructing the description of the present application.

It should be understood that when used in the present specification and the appended claims, the term "comprising" indicates the presence of described features, wholes, steps, operations, elements and/or components, but does not exclude the presence or addition of one or more other features, wholes, steps, operations, elements, components and/or combinations thereof.

It should also be understood that the term “and/or” used in the specification and appended claims refers to any and all possible combinations of one or more of the associated listed items, and includes these combinations.

As used in the specification and appended claims of this application, the term "if" can be interpreted as "when" or "uponce" or "in response to determining" or "in response to detecting", depending on the context. Similarly, the phrase "if it is determined" or "if [described condition or event] is detected" can be interpreted as meaning "uponce it is determined" or "in response to determining" or "uponce [described condition or event] is detected" or "in response to detecting [described condition or event]", depending on the context.

In addition, in the description of the present application specification and the appended claims, the terms "first", "second", "third", etc. are only used to distinguish the descriptions and cannot be understood as indicating or implying relative importance.

References to "one embodiment" or "some embodiments" etc. described in the specification of this application mean that one or more embodiments of the present application include specific features, structures or characteristics described in conjunction with the embodiment. Therefore, the statements "in one embodiment", "in some embodiments", "in some other embodiments", "in some other embodiments", etc. that appear in different places in this specification do not necessarily refer to the same embodiment, but mean "one or more but not all embodiments", unless otherwise specifically emphasized in other ways. The terms "including", "comprising", "having" and their variations all mean "including but not limited to", unless otherwise specifically emphasized in other ways.

In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present application will be clearly and completely described below in conjunction with the accompanying drawings and specific implementation methods. Obviously, the described embodiments are only 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 making creative work are within the scope of protection of the present invention.

Introduction The present invention involves professional terms:

The lower four zones: refers to the division of the rock strata between the overlying confined aquifer and the bottom plate of the mined coal seam into four component zones, including the mine pressure damage zone, the newly added damage zone, the original damage zone and the original high-conductivity zone.

The coal seam roof includes the pseudo roof, direct roof and old roof (basic roof).

Water channels: joints and fissures in rock strata, mining fissures, dissolution fissures in soluble rocks, collapse columns and fault structures, etc.

FIG1 is a schematic diagram of the structure of a water pressure and seismic multi-field coupled water passage detection system provided by an embodiment of the present application. Referring to FIG1 , the water pressure and seismic multi-field coupled water passage detection method is described in detail as follows:

Step 101, obtaining microseismic data, the support pressure of the working face and the water inflow of the working face.

Exemplarily, a microseismic monitoring system is first arranged in the working face. The microseismic monitoring system may include detectors, monitoring substations, etc. As shown in FIG2 , multiple detectors are arranged in the return air lane and machine lane of the working face, with a spacing of 80 to 100 m. Two acquisition substations are generally set.

After the monitoring system collects the original microseismic data, it is uploaded to the ground server via optical fiber, and then pre-processed in professional processing software to perform microseismic data quality control.

In order to accurately obtain the supporting pressure and water inflow of the working face, stress monitoring equipment can be used to monitor the stress state of the coal rock mass in real time. Appropriate water inflow observation methods can be selected according to specific circumstances, such as volumetric method, water pump displacement method, buoy method, water tank water level method, weir measurement method and flow meter method.

Preprocessing can include noise reduction and signal enhancement of the original microseismic data to obtain microseismic data. Among them, noise reduction can be achieved through wavelet denoising functions, such as Daubechies, Symlet, Coiflet and other functions; signal enhancement can be achieved through signal enhancement models, such as unsupervised learning modules, generative adversarial networks and other models. Preprocessing can remove noise and interference, highlight the characteristics of microseismic signals, simplify data representation, and improve computational efficiency.

Step 102: determining microseismic response characteristics based on the microseismic data.

Among them, the microseismic response characteristics include the location, frequency and energy of microseismic events.

Exemplarily, step 102 includes:

First, the microseismic data is quality-controlled to obtain microseismic attribute data, which includes the frequency and energy of microseismic events.

Acquire mining engineering data. Mining engineering data includes mining engineering plan, spatial coordinates of measurement guide points, and histograms of adjacent boreholes on the working face.

Then, a 3D geological model is constructed based on the microseismic attribute data and mining engineering data. The initial 3D geological model can be constructed based on the mining engineering plan, the spatial coordinates of the measurement guide points, the column chart of the adjacent boreholes of the working face, etc., and then the spatial position information in the microseismic attribute data is imported into the initial 3D geological model to form a complete 3D geological model for analyzing microseismic events.

Finally, based on the 3D geological model, the location of microseismic events is identified, which include top plate events and bottom plate events.

Exemplarily, based on the three-dimensional geological model, identifying the location of a microseismic event includes:

Based on the three-dimensional geological model, the layers of the top plate events are divided according to the upper three zones and lithology, and the layers of the bottom plate events are divided according to the development of the aquifer. Take the three-layer thin limestone aquifer and the thick aquifer as examples, see Table 1.

Table 1 Intervals of microseismic events

Based on the layer segments of the top plate events, the occurrence location of the top plate events is determined. Based on the layer segments of the bottom plate events, the occurrence location of the bottom plate events is determined.

In this embodiment, the microseismic response characteristics such as the location, frequency and energy of the microseismic data are extracted through the monitored microseismic data, the geological conditions are analyzed in real time, and the possible bottom plate water inrush risk is predicted.

Step 103, based on the microseismic response characteristics and the support pressure of the working face, determine the pressure migration distribution area.

Periodic pressure refers to the phenomenon that when the coal mining face continues to advance and the span of the old roof reaches a certain length, the old roof in the roof will break and collapse along the coal wall or even in the coal wall under the action of its own weight and the load of the overlying rock strata. This kind of periodic breaking or collapse of the old roof is called periodic pressure of the old roof.

Pressure migration refers to the change or transfer of periodic pressure in time and space during the advancement of the working face. This migration may be related to many factors, such as coal seam inclination, geological structure, mining method, etc. By monitoring the spatial distribution of microseismic events, it is possible to determine the occurrence of periodic pressure in different areas. If microseismic events occur in a certain area, it may indicate that periodic pressure is about to occur in that area.

Therefore, this embodiment provides a pressure migration law, and determines the initial pressure migration distribution area and the final pressure migration distribution area (the final pressure migration distribution area is referred to as the pressure migration distribution area) according to the pressure migration law. The initial pressure migration distribution area is the range enclosed by the three ellipses a, b and c moving along the mining direction as shown in FIG4, and the final pressure migration distribution area is determined, and the spatial distribution of microseismic events along the mining direction as shown in FIG4 is concentrated in the range enclosed by the last ellipse c.

Exemplarily, step 103 includes:

First, based on the microseismic response characteristics, the pressure determination index is set.

The pressure determination index is that the frequency of microseismic events in the roof near the coal seam and the sum of the microseismic event energy in a single day both increase by more than a first preset ratio. The first preset ratio can be set to 20%.

When the pressure judgment index is met, a pressure curve chart is drawn based on the support pressure of the working surface.

Next, based on the pressure judgment index and the pressure curve, the pressure situation of the working face is verified according to the mining direction, and the pressure migration distribution area is determined.

For example, the microseismic events and the pressure curve diagram of the roof near the coal seam are combined, such as the microseismic events (small dots) and the pressure curve diagram (square dots) in Figure 4. The location with a higher density of microseismic events and the overlapping part of the pressure curve diagram are the initial pressure migration distribution area, and the risk factor in the working face increases when pressure comes. Based on the initial pressure migration distribution area, the final pressure migration position is determined along the mining direction as the final pressure migration distribution area.

The pressure migration law includes that when there is a deep water channel (primary fracture damage area) in the working face, the pressure will show migration characteristics. As the mining progresses, the pressure will migrate to the area of the deep water channel, from a to b, to c, and will no longer migrate after c. Therefore, this embodiment determines the pressure migration distribution area c according to the microseismic response characteristics, the pressure curve diagram and the above-mentioned pressure migration law, that is, determines the range of the planar position of the water channel.

Step 104, based on the water inflow of the working face, the pressure migration distribution area and the microseismic response characteristics, determine the location area of the water passage.

Among them, the formation of water channels is closely related to groundwater dynamic conditions, geological structure and lithology. Therefore, the location area of the water channel is used as the main content of the early warning signal.

Exemplarily, the position area of the water passage includes a planar position area of the water passage and a vertical position area of the water passage.

Next, the plane position area of the water passage is mainly determined first, and then the vertical position area of the water passage is determined by extending downward from the plane position area of the water passage.

Step 104 includes:

Step 1041, based on the microseismic response characteristics, determine the plane cloud map of the warning layer and the warning layer with the deep event density and the vertical cloud map of the warning layer and the warning layer with the deep event density.

According to the location of the microseismic monitoring early warning layer, which is a thin layer of limestone with high strength located deeper than the floor mining damage zone and above the floor Ordovician or Cambrian limestone, for example, the second section F2 of the floor aquifer can be determined as the early warning layer of the working face. If regional governance is carried out on the working face, its regional governance layer can also be used as an early warning layer for microseismic monitoring.

According to the four-zone theory under the floor, the floor mining damage zone is the area where the mining pressure has a significant destructive effect on the floor, that is, the shallow floor events below the warning layer are within the normal damage range of mining. The warning layer and its deeper events indicate the participation of water movement in the deep aquifer or the existence of structures such as faults.

Therefore, the warning layer can be selected according to the actual situation. After the warning layer is selected, when the microseismic event develops to the warning layer, an early warning will be issued, such as the interface warning light of the microseismic monitoring system changing color or flashing.

After the warning layer is determined, the warning layer is analyzed with respect to deep microseismic events, and a plane cloud map of the warning layer and the warning layer with respect to the density of deep events and a vertical cloud map of the warning layer and the warning layer with respect to the density of deep events are drawn. For example, if F2 is the warning layer, a plane cloud map and a vertical cloud map of the event density of microseismic events of F2, F3 and F4 are drawn. The plane cloud map is shown in FIG6 , and the vertical cloud map is shown in FIG7 .

Step 1042, based on the pressure migration distribution area and the early warning layer and the early warning layer and the plane cloud map of the deep event density, determine the plane location area of the water passage.

Exemplarily, step 1042 includes:

Based on the pressure migration distribution area and the density of the plane cloud map of the early warning layer and the early warning layer with deep event density, the initial plane position area of the water channel is determined. The initial plane position area of the water channel can generally be selected as the pressure migration distribution area, because the position with a higher density of microseismic events in the pressure migration distribution area is almost the same as the position with a higher density of the plane cloud map of the early warning layer and the early warning layer with deep event density. The initial plane position area of the water channel can also be based on the pressure migration distribution area and limited by setting the density value of the plane cloud map. The limited initial plane position area of the water channel is generally smaller than the pressure migration distribution area.

When the total frequency increase of microseismic events obtained from the plane cloud map of the early warning layer and the deep event density of the early warning layer exceeds the second preset ratio, the occurrence of microseismic events in the early warning layer continues for a preset time, and the total increase rate of the rear response area and the advance response area in the microseismic monitoring data reaches a preset value, it is determined that the time after this moment is the special attention time for microseismic events. This moment is the moment that meets the following conditions: the total frequency increase of microseismic events exceeds the second preset ratio, the occurrence of microseismic events in the early warning layer continues for a preset time, and the total increase rate of the rear response area and the advance response area in the microseismic monitoring data reaches a preset value. The second preset ratio can be set to 70%.

Among them, the rear response area refers to the area behind the monitoring point where the shock wave propagates and affects the monitoring point after the microseismic event occurs. The "rear" here is relative to the mining position of the working face. The mining position is the dotted line position shown in Figure 4. The direction along the mining position is the front, and vice versa, the side of the goaf of the mining is the rear. The advanced response area refers to the distribution range of the microseismic event in the area affected in front of the mining position.

During the special attention time, based on the maximum density range of the plane cloud map of the early warning layer and the early warning layer with deep event density, the plane location area of the water passage is determined from the initial plane location area. Since there may be two or more density concentration areas in the same initial plane location area, the maximum density range of the plane cloud map can be further selected as the plane location area of the water passage by comparing the size of the kernel density, as shown in the range circled by the irregular curve in Figure 8.

Step 1043, after determining the planar position area of the water passage, determine the vertical position area of the water passage based on the water inflow of the working face and the vertical cloud map of the early warning layer and the deep event density of the early warning layer.

Exemplarily, step 1043 includes:

In the plane position area of the water channel, based on the maximum density of the early warning layer and the vertical cloud map of the deep event density of the early warning layer, the development position of the water channel deep in the coal seam floor is determined.

When the water flow from the working face reaches the lowest point of the water flow variation curve at a location where the water flow channel is developed deep in the coal seam floor, the vertical position area entering the water flow channel at this time is determined.

In one embodiment, in order to make the scheme of the present invention clearer, the working face of Ji 17-33200 of Pingmei No. 10 Mine is used as an example for explanation. Ji 17-33200 is the number of the working face of Ji group of Pingmei No. 10 Mine.

Basic conditions of the Ji17-33200 working face of Pingmei No. 10 Mine: The working face coal seam thickness is 0.8 to 2.8m, and the average working face coal seam thickness is 2.4m. The coal seam strike is 219° to 295°, the coal seam dip is 1.5° to 7°, the average coal seam dip is 3.2°, the working face elevation is -886 to -96m, the ground elevation is +150 to +240m, the dip length is 131.2m, and the strike length is 908.6m. The overlying Ji15-33200 goaf, the interlayer distance between the Ji group and the overlying Wu group coal seams is about 180m, and the Wu group has not been mined. The direct water-filled aquifer of the floor is the upper limestone of the Taiyuan Group of the Carboniferous System, and the indirect aquifer is the lower limestone of the Taiyuan Group and the Cambrian limestone, which is under pressure mining.

Ten geophones are arranged in the return air lane and machine lane of the working face, with a spacing of 80-100m. 1#-20# are the serial numbers of the geophones, and there are two acquisition substations, as shown in Figure 2. The layer segments of microseismic events are shown in Table 2.

Table 2 Intervals of microseismic events

Among them, each distance refers to the distance between the roof or floor and the coal seam.

When the pressure judgment index is met, a pressure curve is drawn based on the support pressure of the working face. The frequency increase of microseismic events in the R1 layer and the change in the sum of the daily energy of the roof are used, as shown in the bar chart in Figure 3. The pressure on the working face is determined, and microseisms generally appear 1-2 days in advance. And according to the distribution of microseismic events, there are obvious migration characteristics, as shown in the line chart in Figure 3 and the line (square points) in Figure 4.

According to the migration characteristics, the initial pressure migration distribution area is determined to be a, b and c in Figure 4. The spatial distribution of microseismic events along the mining direction is concentrated in c, so the pressure migration distribution area is determined to be c, and the range of the plane position of the water channel is determined.

Then, based on the microseismic response characteristics, the plane cloud map of the warning layer and the warning layer with deep event density and the vertical cloud map of the warning layer and the warning layer with deep event density are determined.

In this embodiment, the microseismic warning layer is F2, as shown in Table 2. The total frequency change of the bottom plate events is shown in Figure 5. The total frequency of the bottom plate events on July 26 increased by 1242% compared with the previous one. After July 27, warning layer events began to appear in the advanced impact range, and then appeared continuously. The warning layer refers to F3 and F4 in depth. The plane cloud map of the warning layer and the warning layer in depth event density is shown in Figure 6. The vertical cloud map of the warning layer and the warning layer in depth event density is shown in Figure 7. The layer segment identification code in Figure 7 points to the lower boundary of this layer segment.

The plane position area of the water channel is the range encircled by the irregular curve in Figure 8. After the plane position area of the water channel, the vertical position area of the water channel is determined. As shown in Figure 9, the water volume change curve of the water inflow of the working face shows that the water volume of the limestone decreased from September 2 to October 10, and the corresponding mining line entered the abnormal area of the microseismic middle section. Under the influence of the advanced support pressure of the working face, the cracks in the abnormal area were compressed as a whole, and the cracks in the local channel were closed, resulting in a decrease in water volume. After October 10, in the abnormal area of the microseismic middle section, the area of the advanced support pressure-affected area was smaller than the stress release area and gradually decreased. The cracks were reopened, the local channels were formed again, and the water volume increased. Combined with the above situation, the superposition display and comprehensive analysis were performed to locate and explain the possible water channel. The plane position of the deep water channel is 491 to 596 meters in front of the cut-away view and 40 to 62 meters vertically above the coal seam floor. The cracks here are the most developed, and the floor water gushes out from here and then spreads in the upper section of the limestone.

When the water inflow of the working face reaches the lowest point of the V-shaped water volume change curve, the water volume no longer decreases, indicating that the dynamic water volume reduction area and the dynamic water volume increase area are balanced. As mining progresses, the area affected by the advance support pressure is smaller than the stress release area and gradually decreases. The cracks are reopened, the local channels are unblocked again, and the water volume increases again. The change in water inflow of the working face presents a typical V-shaped water volume change curve, and the stress change has undergone a process from being dominated by compressive stress, gradually balanced, and then changed to being dominated by tensile stress. Therefore, when the water inflow of the working face reaches the lowest point of the water volume change curve, the vertical position area entering the water channel at this time is determined. The model of the vertical position area of the water channel is shown in Figure 10.

The water pressure and seismic multi-field coupling water channel detection method provided in the embodiment of the present application can monitor the fracture of the roof and floor surrounding rock in real time during mining by acquiring microseismic data, water level and water volume changes, and the law of mine pressure manifestation, and can accurately analyze it in all directions, use microseismic indicators to judge the pressure, and at the same time, according to the discovery that the pressure manifestation has a spatial migration law and migrates to the deep channel damage area, the pressure migration area is obtained, and the situation of the coal seam floor water channel is analyzed by combining the microseismic early warning layer and its deep event density cloud map. The plane and vertical position of the channel are comprehensively analyzed by multi-information of water volume, mine pressure, and microseismic, so as to achieve accurate prediction and determination of the position of the water channel, and can identify the possible location area of the water channel in advance, so as to conduct risk assessment and early warning in time. This is of great significance for preventing water disasters and ensuring engineering safety.

It should be understood that the order of execution of the above steps does not necessarily mean the order in which they are executed. The order in which each process is executed should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application.

Referring to Figure 11, this embodiment provides a water pressure and seismic multi-field coupling water passage detection device, which implements the water pressure and seismic multi-field coupling water passage detection method as in the above embodiment. The water pressure and seismic multi-field coupling water passage detection device includes a data acquisition module 201, a feature determination module 202, a pressure migration distribution area determination module 203 and a water passage determination module 204.

The data acquisition module 201 is used to acquire microseismic data, the support pressure of the working face and the water inflow of the working face.

The feature determination module 202 is used to determine the microseismic response feature based on the microseismic data.

The pressure migration distribution area determination module 203 is used to determine the pressure migration distribution area based on the microseismic response characteristics and the support pressure of the working face.

The water passage determination module 204 is used to determine the location area of the water passage based on the water inflow of the working face, the pressure migration distribution area and the microseismic response characteristics; the location area of the water passage serves as the main content of the early warning signal.

Exemplarily, the microseismic response characteristics include the occurrence location, frequency and energy of the microseismic event;

The feature determination module 202 is specifically used for:

Perform quality control processing on microseismic data to obtain microseismic attribute data; microseismic attribute data includes the frequency and energy of microseismic events; obtain mining engineering data; construct a three-dimensional geological model based on the microseismic attribute data and mining engineering data; identify the location of microseismic events based on the three-dimensional geological model.

Exemplarily, microseismic events include top plate events and bottom plate events;

In the feature determination module 202, based on the three-dimensional geological model, the location of the microseismic event is identified, including:

Based on the three-dimensional geological model, the layers of roof events are divided according to the upper three zones and lithology, and the layers of floor events are divided according to the development of aquifers; based on the layers of roof events, the occurrence location of roof events is determined; based on the layers of floor events, the occurrence location of floor events is determined.

Exemplarily, the pressure migration distribution area determination module 203 is specifically used to:

Based on the microseismic response characteristics, the pressure determination index is set; when the pressure determination index is met, a pressure curve chart is drawn based on the support pressure of the working face; based on the pressure determination index and the pressure curve chart, the pressure situation of the working face is verified according to the mining direction, and the pressure migration distribution area is determined.

Exemplarily, the position area of the water passage includes a planar position area of the water passage and a vertical position area of the water passage.

The water passage determination module 204 is specifically used for:

Based on the microseismic response characteristics, the plane cloud map of the warning layer and the warning layer with deep event density and the vertical cloud map of the warning layer and the warning layer with deep event density are determined; based on the pressure migration distribution area and the plane cloud map of the warning layer and the warning layer with deep event density, the plane position area of the water passage is determined; after determining the plane position area of the water passage, the vertical position area of the water passage is determined based on the water inflow of the working face and the vertical cloud map of the warning layer and the warning layer with deep event density.

Exemplarily, in the water channel determination module 204, based on the pressure migration distribution area and the early warning layer and the early warning layer and the plane cloud map of the deep event density, the plane location area of the water channel is determined, including:

Based on the pressure migration distribution area and the density of the plane cloud map of the early warning layer and the early warning layer with deep event density, the initial plane position area of the water passage is determined; based on the maximum density range of the plane cloud map of the early warning layer and the early warning layer with deep event density, the plane position area of the water passage is determined from the initial plane position area.

Exemplarily, after determining the initial plane position area of the water passage, the water passage determination module 204 is further used to:

When the total frequency increase of microseismic events obtained from the plane cloud map of the early warning layer and the deep event density above the early warning layer exceeds the second preset ratio, the occurrence of microseismic events in the early warning layer continues for a preset time, and the total increase rate of the rear response area and the advance response area in the microseismic monitoring data reaches a preset value, it is determined that the time after this moment is the special attention time for microseismic events;

Accordingly, based on the maximum density range of the plane cloud map of the early warning layer and the early warning layer with deep event density, the plane location area of the water passage is determined from the initial plane location area, including:

During the special attention time, based on the maximum density range of the plane cloud map of the early warning layer and the early warning layer with deep event density, the plane location area of the water passage is determined from the initial plane location area.

Exemplarily, based on the water inflow of the working face and the warning layer and the vertical cloud map of the warning layer and the deep event density, the vertical position area of the water passage is determined, including:

In the plane position area of the water channel, based on the maximum density of the vertical cloud map of the early warning layer and the deep event density of the early warning layer, the development position of the water channel at the depth of the coal seam floor is determined;

When the water flow rate at the working face reaches the lowest point of the water flow change curve at a location where the water flow channel is developed deep in the coal seam floor, the vertical position area entering the water flow channel at this time is determined.

It should be noted that the information interaction, execution process and other contents between the above-mentioned devices are based on the same concept as the method embodiment of the present application. Their specific functions and technical effects can be found in the method embodiment part and will not be repeated here.

The technicians in the relevant field can clearly understand that for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example for illustration. In practical applications, the above-mentioned function allocation can be completed by different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiment can be integrated in a processing unit, or each unit can exist physically separately, or two or more units can be integrated in one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of this application. The specific working process of the units and modules in the above-mentioned system can refer to the corresponding process in the aforementioned method embodiment, which will not be repeated here.

An embodiment of the present application also provides a terminal, see Figure 12, the terminal 400 may include: at least one processor 410 and a memory 420, the memory 420 stores a computer program 421 that can be executed on at least one processor 410, and when the processor 410 executes the computer program 421, it implements the steps in any of the above-mentioned method embodiments, such as steps 101 to 104 in the embodiment shown in Figure 1, and modules 201 to 204 in the embodiment shown in Figure 11.

Exemplarily, the computer program 421 may be divided into one or more modules/units, one or more modules/units are stored in the memory 420, and are executed by the processor 410 to complete the present application. One or more modules/units may be a series of computer program segments capable of completing specific functions, and the program segments are used to describe the execution process of the computer program 421 in the terminal 400.

Those skilled in the art will understand that FIG12 is merely an example of a terminal and does not constitute a limitation on the terminal. The terminal may include more or fewer components than shown in the figure, or a combination of certain components, or different components, such as input and output devices, network access devices, buses, etc.

The processor 410 may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc.

The memory 420 may be an internal storage unit of the terminal or an external storage device of the terminal, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc. The memory 420 is used to store computer programs and other programs and data required by the terminal. The memory 420 may also be used to temporarily store data that has been output or is to be output.

The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, the bus in the drawings of this application is not limited to only one bus or one type of bus.

The water pressure-vibration multi-field coupled water channel detection method provided in the embodiments of the present application can be applied to terminal devices such as computers, tablet computers, laptops, netbooks, personal digital assistants (PDAs), etc. The embodiments of the present application do not impose any restrictions on the specific types of terminals.

An embodiment of the present application also provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, it can implement the steps in each embodiment of the above-mentioned water pressure shock multi-field coupling water channel detection method.

An embodiment of the present application provides a computer program product. When the computer program product runs on a mobile terminal, the mobile terminal implements the steps in each embodiment of the water-passing channel detection method with multi-field coupling of water pressure and shock when executing the computer program product.

If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the present application implements all or part of the processes in the above-mentioned embodiment method, which can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by the processor, the steps of the above-mentioned method embodiments can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium may at least include: any entity or device that can carry the computer program code to the camera device/terminal device, recording medium, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electrical carrier signal, telecommunication signal and software distribution medium. For example, a USB flash drive, a mobile hard disk, a disk or an optical disk.

In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described or recorded in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this application.

In the embodiments provided in the present application, it should be understood that the disclosed devices/network equipment and methods can be implemented in other ways. For example, the device/network equipment embodiments described above are only schematic. For example, the division of modules or units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

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

Claims (9)

1. The water channel detection method for hydraulic shock multi-field coupling is characterized by comprising the following steps of:

acquiring microseismic data, supporting pressure of a working surface and water inflow of the working surface;

determining microseismic response characteristics based on the microseismic data;

Determining a zone of migration in pressure based on the microseismic response feature and a bearing pressure of the working surface;

Determining a position area of a water passing channel based on the water inflow of the working surface, the pressure migration distribution area and the microseismic response characteristics; the position area of the water passing channel is used as the main content of the early warning signal; the position area of the water channel comprises a plane position area of the water channel and a vertical position area of the water channel;

the determining a location area of the water channel based on the water inflow of the working surface, the pressure migration distribution area and the microseismic response characteristics comprises:

Based on the microseismic response characteristics, determining a planar cloud image of the early warning layer and the early warning layer with the deep event density and a vertical cloud image of the early warning layer and the early warning layer with the deep event density; the early warning layer is positioned at the position of a thin layer of limestone with certain strength above the bottom plate Ore or the chilly limestone with deep mining damage zone;

determining a plane position area of the water channel based on the incoming pressure migration distribution area, the early warning layer and a plane cloud image of the early warning layer with deep event density;

after the plane position area of the water passing channel is determined, the vertical position area of the water passing channel is determined based on the water inflow of the working face and the vertical cloud image of the early warning layer and the early warning layer with deep event density.

2. The method for detecting a water channel coupled by a hydraulic shock multi-field according to claim 1, wherein the microseismic response characteristics comprise occurrence position, frequency and energy of a microseismic event;

the determining the microseismic response feature based on the microseismic data comprises:

Performing quality control processing on the microseismic data to obtain microseismic attribute data; the microseismic attribute data comprises the frequency and energy of microseismic events;

acquiring mining engineering data;

Constructing a three-dimensional geological model based on the microseismic attribute data and the mining engineering data;

and identifying the occurrence position of the microseismic event based on the three-dimensional geological model.

3. The water shock multi-field coupled water channel detection method of claim 2, wherein the microseismic event comprises a roof event and a floor event;

the identifying, based on the three-dimensional geological model, the occurrence location of the microseismic event includes:

Dividing the intervals of the roof events according to the upper three zones and lithology based on the three-dimensional geological model, and dividing the intervals of the floor events according to the development condition of the aquifer;

determining an occurrence location of the roof event based on the interval of the roof event;

Determining an occurrence position of the floor event based on the interval of the floor event.

4. The method of water pressure shock multiple field coupled water channel detection of claim 1, wherein the determining a zone of pressure migration distribution based on the microseismic response characteristics and the bearing pressure of the working surface comprises:

Setting a pressure judgment index based on the microseismic response characteristics;

drawing a graph based on the supporting pressure of the working surface when Fu Gelai presses the judging index;

and verifying the pressure condition of the working surface according to the stoping direction based on the pressure judgment index and the pressure curve diagram, and determining a pressure migration distribution area.

5. The method for detecting a water channel through water vibration multi-field coupling according to claim 2, wherein the determining the planar position area of the water channel based on the incoming pressure migration distribution area and the planar cloud patterns of the early warning layer and the early warning layer with deep event density comprises:

Based on the density of the incoming pressure migration distribution area and the planar cloud patterns of the early warning layer and the early warning layer with the deep event density, determining an initial planar position area of the water channel;

And determining the plane position area of the water channel from the initial plane position area based on the early warning layer and the maximum density range of the plane cloud image of the early warning layer with the deep event density.

6. The water-vibration multi-field coupled water channel detection method as claimed in claim 5, wherein after said determining the initial planar position area of the water channel, the water-vibration multi-field coupled water channel detection method further comprises:

When the total frequency amplification of the microseismic event exceeds a second preset proportion, the occurrence of the microseismic event of the early warning layer continues for a preset time and the total increase rate of the rear response area and the advanced response area range in the microseismic monitoring data reaches a preset value, determining the special attention time of the microseismic event after the moment;

correspondingly, the determining the plane position area of the water channel from the initial plane position area based on the early warning layer and the maximum density range of the plane cloud image of the early warning layer with the deep event density comprises the following steps:

And determining the plane position area of the water channel from the initial plane position area based on the maximum density range of the plane cloud patterns of the early warning layer and the early warning layer in the deep event density in the special attention time.

7. The method for detecting a water channel through water vibration multi-field coupling according to claim 1, wherein determining a vertical position area of the water channel based on the water inflow of the working surface and vertical cloud patterns of the early warning layer and the early warning layer with deep event density comprises:

in the plane position area of the water channel, judging the development position of the water channel in the deep coal seam bottom plate based on the position of the early warning layer and the vertical cloud image with the highest density of the early warning layer in the deep event density;

And determining a vertical position area entering the water passing channel at the moment when the water inflow of the working surface reaches the lowest point of the water inflow change curve at the deep development position of the water passing channel on the coal bed bottom plate.

8. A terminal comprising a memory and a processor, the memory having stored thereon a computer program operable on the processor, wherein the processor when executing the computer program implements the water channel detection method of hydraulic shock multi-field coupling as claimed in any one of claims 1 to 7.

9. A computer readable storage medium storing a computer program, wherein the computer program when executed by a processor implements the water channel detection method of hydraulic shock multi-field coupling according to any one of claims 1 to 7.