Method and device for determining regional vertical stress based on point-plane in-situ data fusion

Abstract

This application relates to a method and device for determining the regional vertical stress based on point-plane in-situ data fusion. The specific solution is as follows: Determine multiple target layers within the target area of the coal mining face according to preset rules, determine the target detection plane based on the multiple target layers, and respectively determine multiple measurement points and their respective horizontal coordinate values on the target detection plane. Obtain the wave velocity distribution information of the target area, the coal and rock mass density data of each measurement point, the wave velocity data of each measurement point, and the vertical stress data of each measurement point respectively. Based on the corresponding relationships between the coal and rock mass density data of each measurement point, the wave velocity data corresponding to each measurement point, the vertical stress data of each measurement point, and the wave velocity gradient of the target area, respectively determine the vertical stress distribution information of each of the multiple target layers. Based on the vertical stress distribution information of each of the multiple target layers, determine the vertical stress distribution information of the target area. This application can determine the continuous vertical stress distribution in the coal mining area.

Description

Technical Field

[0001] This application relates to the technical field of in-situ stress measurement in underground coal mine engineering areas, and particularly to a method and device for determining regional vertical stress based on point-plane in-situ data fusion. Background Art

[0002] In related technologies, the physical property variation characteristics of coal and rock masses are prominent. At the same time, affected by the undulation of coal seams and geological structures, the local vertical stress measured at individual positions is difficult to meet the engineering requirements. The seismic wave CT technology can realize the detection of a large area in the working face. Its principle is to reflect the relative stress accumulation distribution in the area through the change characteristics of the seismic wave propagation speed in the coal and rock masses, but it cannot obtain continuous numerical values of the regional vertical stress. Summary of the Invention

[0003] Therefore, this application provides a method and device for determining regional vertical stress based on point-plane in-situ data fusion. The technical solution of this application is as follows:

[0004] According to the first aspect of the embodiments of this application, a method for determining regional vertical stress based on point-plane in-situ data fusion is provided. The method includes:

[0005] Determine multiple target layers in the target area of the coal mine working face according to a preset rule, determine multiple target detection surfaces based on the multiple target layers, and respectively determine multiple measurement points and the three-dimensional coordinates of each of the multiple measurement points on each target detection surface;

[0006] Respectively obtain the wave velocity distribution information of the target area, the coal and rock mass density data of each of the multiple measurement points, the wave velocity data of each of the multiple measurement points, and the vertical stress data of each of the multiple measurement points; wherein, the wave velocity distribution information includes a wave velocity distribution nephogram and wave velocity gradient information;

[0007] Based on the corresponding relationships between the coal and rock mass density data of each of the multiple measurement points, the wave velocity data corresponding to each of the multiple measurement points, the vertical stress data of each of the multiple measurement points, and the wave velocity gradient of the target area, respectively determine the vertical stress distribution information of each of the multiple target detection surfaces; wherein, the vertical stress distribution information includes a vertical stress distribution nephogram and vertical stress information;

[0008] Based on the vertical stress distribution information of each of the multiple target detection surfaces, determine the vertical stress distribution information of the target area.

[0009] According to an embodiment of the present application, the vertical stress distribution information of each of the plurality of target detection surfaces is determined respectively based on the corresponding relationship between the coal and rock mass density data of each of the plurality of measurement points, the wave velocity data corresponding to each of the plurality of measurement points, the vertical stress data of each of the plurality of measurement points, and the wave velocity gradient of the target area, including:

[0010] For each target layer, based on the corresponding relationship between the coal and rock mass density data of each of the plurality of measurement points in the target layer, the wave velocity data corresponding to each of the plurality of measurement points, the vertical stress data of each of the plurality of measurement points, and the wave velocity gradient of the target layer, the vertical stress data corresponding to each wave velocity data is determined;

[0011] Based on the vertical stress data corresponding to each wave velocity data, the wave velocity distribution data in the wave velocity distribution cloud map is replaced with the vertical stress distribution data to obtain the vertical stress distribution cloud map of the target layer.

[0012] According to an embodiment of the present application, determining a plurality of target layers in the target area of the coal mining face according to a preset rule and determining a plurality of target detection surfaces based on the plurality of target layers includes:

[0013] The thickness of each target layer is obtained respectively;

[0014] Based on the thickness of each target layer, the middle position of each target layer is determined as a target detection surface respectively;

[0015] The interface between every two adjacent target layers is determined as a target detection surface respectively;

[0016] Based on the thickness of each target layer, the z-axis coordinate of each target detection surface is determined respectively.

[0017] According to an embodiment of the present application, determining a plurality of target layers in the target area of the coal mining face according to a preset rule includes:

[0018] The total thickness of the coal seam in the target area is obtained;

[0019] The preset number of target layers is obtained;

[0020] The total thickness of the coal seam is divided by the preset number of target layers to obtain the coal seam thickness of each target layer.

[0021] According to an embodiment of the present application, the wave velocity distribution information of the target area is obtained by the following method:

[0022] Determine an excitation point at the middle position of each target layer; the x-axis coordinates of each of the excitation points are the same; the y-axis coordinates of each of the excitation points are the same; each of the excitation points is located on the side of the seismic wave excitation end roadway close to the seismic wave receiving end roadway;

[0023] Determine a receiving point at the middle position of each target layer; the x-axis coordinates of each of the receiving points are the same; the y-axis coordinates of each of the receiving points are the same; each of the receiving points is located on the side of the seismic wave receiving end roadway close to the working face;

[0024] Obtain the seismic wave data transmitted by the receiving end geophones corresponding to each target layer respectively;

[0025] Based on the seismic wave data corresponding to each target layer, determine the wave velocity distribution information of the target area.

[0026] According to an embodiment of the present application, the determining the vertical stress distribution information of the target area based on the vertical stress distribution information of each of the multiple coal seams further includes:

[0027] Based on the three-dimensional coordinates of each of the multiple measurement points, divide the measurement points with the same z-axis coordinate value into one group;

[0028] For each group, based on the vertical stress data of each of the multiple measurement points, determine the vertical stress difference between every two adjacent measurement points respectively;

[0029] Based on the vertical stress differences between every two adjacent measurement points, calculate the vertical stress change rate of each measurement point in the z-axis direction respectively;

[0030] Integrate the vertical stress change rate of the measurement points in the z-axis direction to obtain an integration result;

[0031] Based on the integration result of each measurement point, determine the vertical stress continuous function expression between each measurement point and its adjacent measurement point in the z-axis direction;

[0032] Merge the vertical stress continuous function expressions of multiple measurement points in each group to obtain the vertical stress function expression of each group;

[0033] According to the vertical stress function expression of each group, determine the vertical stress distribution information of the target area.

[0034] According to an embodiment of the present application, the respectively determining multiple measurement points and the three-dimensional coordinates of each of the multiple measurement points on each target detection surface includes:

[0035] For each target detection surface, the target is evenly divided into multiple blocks in layers, and the coordinate regions of the multiple blocks are respectively determined;

[0036] The coordinate regions of the multiple blocks are compared with a preset region. In response to the coordinate region falling within the preset coordinate region, the block corresponding to the coordinate region is removed to obtain the multiple blocks; wherein, the preset coordinate region is a wave velocity anomaly region;

[0037] Based on the coordinate regions of the multiple blocks, the central point coordinates of the multiple blocks are respectively determined;

[0038] The points corresponding to the central point coordinates of the multiple blocks are determined as the measuring points of the target detection surface.

[0039] According to the second aspect of the embodiments of the present application, a device for determining regional vertical stress based on point-surface in-situ data fusion is provided. The device includes:

[0040] A first determination module, configured to determine multiple target layers in a target area of a coal mining face according to a preset rule, determine multiple target detection surfaces based on the multiple target layers, and respectively determine multiple measuring points and the three-dimensional coordinates of the multiple measuring points on each target detection surface;

[0041] An acquisition module, configured to respectively acquire the wave velocity distribution information of the target area, the coal and rock mass density data of the multiple measuring points, the wave velocity data of the multiple measuring points, and the vertical stress data of the multiple measuring points; wherein, the wave velocity distribution information includes a wave velocity distribution nephogram and wave velocity gradient information;

[0042] A second determination module, configured to respectively determine the vertical stress distribution information of the multiple target detection surfaces based on the corresponding relationship between the coal and rock mass density data of the multiple measuring points, the wave velocity data corresponding to the multiple measuring points, the vertical stress data of the multiple measuring points, and the wave velocity gradient of the target area; wherein, the vertical stress distribution information includes a vertical stress distribution nephogram and vertical stress information;

[0043] A third determination module, configured to determine the vertical stress distribution information of the target area based on the vertical stress distribution information of the multiple target detection surfaces.

[0044] According to the third aspect of the embodiments of the present application, an electronic device is provided, including: a processor and a memory communicatively connected to the processor;

[0045] The memory stores computer-executable instructions;

[0046] The processor executes the computer-executable instructions stored in the memory to implement the method according to any one of the first aspect.

[0047] According to a fourth aspect of the embodiments of the present application, there is provided a computer-readable storage medium, characterized in that the computer-readable storage medium stores computer-executable instructions, and when the computer-executable instructions are executed by a processor, they are used to implement the method according to any one of the first aspect.

[0048] The technical solutions provided by the embodiments of the present application at least bring the following beneficial effects:

[0049] By determining a plurality of target stratifications in a target area of a coal mining face according to a preset rule, determining a plurality of target detection surfaces based on the plurality of target stratifications, and respectively determining a plurality of measurement points and the three-dimensional coordinates of each of the plurality of measurement points on each target detection surface. Respectively obtain the wave velocity distribution information of the target area, the density data of the coal and rock mass of each of the plurality of measurement points, the wave velocity data of each of the plurality of measurement points, and the vertical stress data of each of the plurality of measurement points. Based on the corresponding relationships between the density data of the coal and rock mass of each of the plurality of measurement points, the wave velocity data corresponding to each of the plurality of measurement points, the vertical stress data of each of the plurality of measurement points, and the wave velocity gradient of the target area, respectively determine the vertical stress distribution information of each of the plurality of target detection surfaces. Based on the vertical stress distribution information of each of the plurality of target detection surfaces, determine the vertical stress distribution information of the target area. Thus, the quantitative evaluation of the regional vertical stress is realized, which intuitively provides reliable data support for the mine production design, further significantly improves the mine design and production efficiency, greatly reduces the risk of disasters, and avoids the generation of ineffective workload.

[0050] It should be understood that the above general description and the following detailed description are only exemplary and explanatory, and cannot limit the present application. Description of the Drawings

[0051] The drawings here are incorporated into the specification and constitute a part of this specification, showing embodiments consistent with the present application, and are used together with the specification to explain the principles of the present application, and do not constitute an improper limitation to the present application.

[0052] Figure 1 It is a flowchart of a method for determining regional vertical stress based on point-plane in-situ data fusion in an embodiment of the present application;

[0053] Figure 2 It is a schematic diagram of the layout method for coal seam seismic wave CT stratification detection proposed in an embodiment of the present application;

[0054] Figure 3 It is a wave velocity distribution cloud map and wave velocity gradient obtained by seismic wave CT detection proposed in an embodiment of the present application;

[0055] Figure 4Schematic vertical section diagram of the layout of vertical stress measurement points proposed in the embodiments of the present application;

[0056] Figure 5 Schematic horizontal section diagram of the layout of vertical stress measurement points proposed in the embodiments of the present application;

[0057] Figure 6 Schematic connection diagram of the vertical stress direction proposed in the embodiments of the present application;

[0058] Figure 7 Structural block diagram of a regional vertical stress determination device based on point-plane in-situ data fusion in the embodiments of the present application;

[0059] Figure 8 Block diagram of an electronic device in the embodiments of the present application.

[0060] Reference Signs

[0061] 4 - Wave velocity gradient anomaly region; 5 - Working face coordinate system; 6 - Wave velocity gradient boundary line; 7 - Block; 8 - Local measurement point; 10 - Vertical stress direction of the measurement point; 11 - First layer; 12 - Second layer; 13 - Third layer; 14 - Receiving side roadway; 15 - Exciting side roadway; 16 - Geophone; 17 - Bolt; 18 - Anchorage section; 19 - Borehole; 20 - Explosive; 21 - Vertical section line; 23 - First target layer; 24 - Second target layer; 25 - Interface. Detailed implementation manners

[0062] In order to enable those of ordinary skill in the art to better understand the technical solutions of the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings.

[0063] It should be noted that the terms "first", "second", etc. in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and do not have to be used to describe a specific order or sequence. It should be understood that such used data can be interchanged under appropriate circumstances so that the embodiments of the present application described here can be implemented in an order different from those illustrated or described here. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present application. On the contrary, they are only examples of devices and methods consistent with some aspects of the present application as detailed in the appended claims.

[0064] It should be noted that in the related art, the physical property variation characteristics of coal and rock masses are prominent. At the same time, affected by the undulation of coal seams and geological structures, the local vertical stress measured at individual positions is difficult to meet the engineering requirements. The seismic wave CT technology can realize the detection of a large area of the working face. Its principle is to reflect the relative stress accumulation distribution in this area through the change characteristics of the seismic wave propagation speed in the coal and rock masses, but it cannot obtain the continuous numerical value of the regional vertical stress.

[0065] Based on the above problems, the present application proposes a method and device for determining the regional vertical stress based on the fusion of point and surface in-situ data. It can realize determining multiple target layers in the target area of the coal mine working face according to a preset rule, determining multiple target detection surfaces based on the multiple target layers, and respectively determining multiple measuring points and the three-dimensional coordinates of each of the multiple measuring points on each target detection surface. Respectively obtain the wave velocity distribution information of the target area, the coal and rock mass density data of each of the multiple measuring points, the wave velocity data of each of the multiple measuring points, and the vertical stress data of each of the multiple measuring points. Based on the corresponding relationship between the coal and rock mass density data of each of the multiple measuring points, the wave velocity data corresponding to each of the multiple measuring points, the vertical stress data of each of the multiple measuring points, and the wave velocity gradient of the target area, respectively determine the vertical stress distribution information of each of the multiple target detection surfaces. Based on the vertical stress distribution information of each of the multiple target detection surfaces, determine the vertical stress distribution information of the target area. Thus, the quantitative evaluation of the regional vertical stress is realized, intuitively providing reliable data support for the mine production design, further significantly improving the mine design and production efficiency, greatly reducing the risk of disasters, and reducing the ineffective workload.

[0066] Figure 1 It is a flowchart of a method for determining the regional vertical stress based on the fusion of point and surface in-situ data in an embodiment of the present application.

[0067] As Figure 1 shown, the method for determining the regional vertical stress based on the fusion of point and surface in-situ data includes:

[0068] Step 101, determine multiple target layers in the target area of the coal mine working face according to a preset rule, determine multiple target detection surfaces based on the multiple target layers, and respectively determine multiple measuring points and the three-dimensional coordinates of each of the multiple measuring points on each target detection surface.

[0069] Among them, in some embodiments of the present application, determining multiple target layers in the target area of the coal mine working face according to a preset rule and determining multiple target detection surfaces based on the multiple target layers includes:

[0070] Step a1, respectively obtain the thickness of each target layer;

[0071] Among them, in some embodiments of the present application, step a1 includes:

[0072] Step a11, obtaining the total thickness of the coal seam in the target area.

[0073] Step a12, obtaining the preset target number of layers.

[0074] It can be understood that the preset target number of layers can be set in advance according to actual needs. The more the preset target number of layers, the higher the data accuracy.

[0075] Step a13, dividing the total thickness of the coal seam by the preset target number of layers to obtain the thickness of each target layer.

[0076] Step a2, based on the thickness of each target layer, respectively determine the middle position of each target layer as a target detection surface.

[0077] Step a3, respectively determine the interface between every two adjacent target layers as a target detection surface.

[0078] For example, as Figure 4 shown, the middle positions of the first target layer 23 and the second target layer 24 are respectively determined as a target detection surface. The interface 25 between the first target layer 23 and the second target layer 24 is determined as a target detection surface.

[0079] It should be noted that, to increase the continuity of the measurement points, the interface between every two adjacent target layers is respectively determined as a target detection surface.

[0080] Step a4, based on the thickness of each target layer, respectively determine the z-axis coordinates of each target detection surface.

[0081] Among them, in some embodiments of the present application, determining multiple measurement points and the three-dimensional coordinates of each measurement point on each target detection surface respectively includes:

[0082] Step b1, for each target detection surface, evenly divide the target detection surface into multiple blocks, and respectively determine the coordinate regions of the multiple blocks.

[0083] It should be noted that the blocks can cover as much as possible all the areas except the abnormal areas.

[0084] Step b2, compare the coordinate regions of the multiple blocks with the preset region, and in response to the coordinate region falling into the preset coordinate region, remove the block corresponding to the coordinate region to obtain multiple blocks.

[0085] Among them, in the embodiments of the present application, the preset coordinate region is the wave velocity abnormal region.

[0086] It should be noted that, asFigure 3 , Figure 5 As shown, when selecting the location of the measuring point block 7, it is necessary to avoid the abnormal velocity area 4, the microseismic event cluster area and the local geological structure development area in the velocity and velocity gradient cloud map, so as to enhance the reliability and representativeness of the regional test results. It can be understood that the four types of areas of the velocity and velocity gradient distribution cloud map obtained by the original seismic wave CT detection are defined as the velocity thresholds in the region as v1, v2, v3, and v4, respectively, and the area between two adjacent velocity gradient boundaries 6 is a type of area. Optionally, the coordinates (x, y) of the local measuring point 8 can be the coordinates on the working surface coordinate system 5.

[0087] It should be noted that the local measuring point 8 is one of the multiple measuring points, and an example is given for describing one of the multiple measuring points.

[0088] As a possible implementation example, the coordinate areas of the multiple blocks are compared with the preset area. If the coordinate area falls within the preset coordinate area, it means that the block does not meet the requirements and is removed. After the multiple blocks are screened according to the above method, the remaining blocks are the above multiple blocks.

[0089] Optionally, the blocks may be squares with a side length of a preset value, and the preset value is as small as possible to obtain as many blocks as possible.

[0090] Step b3, based on the coordinate areas of the multiple blocks, respectively determine the coordinates of the center points of the multiple blocks.

[0091] Step b4, determining the points corresponding to the coordinates of the center points of the multiple blocks as the measuring points of the target detection surface.

[0092] As an example of possible implementation, based on the coordinate areas of the multiple blocks, the coordinates of the center points of the multiple blocks are determined respectively, and the points corresponding to the coordinates of the center points of the multiple blocks are determined as multiple measuring points in the target area, and each measuring point is used to represent the block to which it belongs.

[0093] Step 102, respectively obtain wave velocity distribution information of the target area, coal rock density data of each of the multiple measuring points, wave velocity data of each of the multiple measuring points, and vertical stress data of each of the multiple measuring points.

[0094] In this embodiment of the present application, the wave velocity distribution information includes a wave velocity distribution cloud map and wave velocity gradient information.

[0095] Among them, in the embodiments of the present application, Figure 3 As shown, the wave velocity distribution information includes the wave velocity distribution cloud map and the wave velocity gradient information.

[0096] Optionally, seismic wave CT in-situ detection can be carried out in the target area of the coal mining face to obtain the wave velocity and wave velocity gradient distribution data of the target area, and then the wave velocity distribution information of the target area and the wave velocity values of each measurement point can be obtained.

[0097] Optionally, boreholes can be drilled at the selected measurement point positions, and then in-situ stress tests can be carried out to obtain the vertical stress of each measurement point.

[0098] In some embodiments of the present application, the wave velocity distribution information of the target area is obtained by the following method:

[0099] Step c1, determine an excitation point at the middle position of each target layer respectively.

[0100] Wherein, in the embodiments of the present application, the x-axis coordinates of each excitation point are the same; the y-axis coordinates of each excitation point are the same; each excitation point is located on the side of the seismic wave excitation end roadway close to the seismic wave receiving end roadway.

[0101] Step c2, determine a receiving point at the middle position of each target layer respectively.

[0102] Wherein, in the embodiments of the present application, the x-axis coordinates of each receiving point are the same; the y-axis coordinates of each receiving point are the same; each receiving point is located on the side of the seismic wave receiving end roadway close to the working face.

[0103] Step c3, obtain the seismic wave data sent by the receiving end geophones corresponding to each target layer respectively.

[0104] Step c4, determine the wave velocity distribution information of the target area based on the seismic wave data corresponding to each target layer.

[0105] For example, as Figure 2 shown, (L - 1) boreholes 19 corresponding to L layers are arranged in the excitation side roadway 15 and the receiving side roadway 14. The arrangement method of the first layer 11 is the same as the conventional single coal seam CT detection arrangement method. The CT detection arrangement methods of the second layer 12 and the third layer 13 are: taking the position 1 - 2 m away from the positive rib of the roadway roof as the borehole opening position, and the boreholes 19 are inclined towards the inside of the working face. The principle for determining the borehole arrangement parameters of the second layer 12 and the third layer 13 is: the ends of each borehole 19 are located at half of the thickness position of the corresponding layer; the horizontal ends of the boreholes 19 at the ends of each layer are located on the same vertical line as the horizontal end of the borehole 19 in the first layer 11; the ends of the boreholes 19 avoid the locally fractured area of the rock stratum, and if it cannot be avoided, grouting reinforcement needs to be carried out on the bottom area of the hole.

[0106] The specific arrangement method of the borehole 19 on the excitation side is as follows: The borehole 19 on the excitation side is used for charging and blasting. The explosive 20 is loaded at the bottom of the borehole. The amount of explosive loaded is the same as that of the first layer. The remaining length of the borehole except the charging section is sealed with stemming to prevent the explosive 20 from slipping and deviating from the predetermined position.

[0107] The specific arrangement method of the borehole 19 on the receiving side is as follows: The borehole is used to arrange the bolt 17. A certain length outward from the bottom of the borehole is the anchorage section 18, and the length except the anchorage section 18 is the free section of the bolt 17. A matching locking device is arranged at the orifice of the borehole 19 and a pre-tightening force is applied to make the bolt 18 fully locked and stressed. The geophone 16 is arranged on the bolt 17 exposed at the orifice of the borehole 19 in the roadway.

[0108] Step 103: Based on the corresponding relationships between the coal and rock mass density data of each of the multiple measurement points, the wave velocity data corresponding to each of the multiple measurement points, the vertical stress data of each of the multiple measurement points, and the wave velocity gradient of the target area, respectively determine the vertical stress distribution information of each of the multiple target detection surfaces.

[0109] Among them, in the embodiments of the present application, the vertical stress distribution information includes a vertical stress distribution nephogram and vertical stress information.

[0110] In some embodiments of the present application, step 103 includes:

[0111] Step d1: For each target layer, based on the corresponding relationships between the coal and rock mass density data of each of the multiple measurement points in the target layer, the wave velocity data corresponding to each of the multiple measurement points, the vertical stress data of each of the multiple measurement points, and the wave velocity gradient of the target layer, determine the vertical stress data corresponding to each wave velocity data.

[0112] Step d2: Based on the vertical stress data corresponding to each wave velocity data, replace the wave velocity distribution data in the wave velocity distribution nephogram with the vertical stress distribution data to obtain the vertical stress distribution nephogram of the target layer.

[0113] As a possible example, the coordinates of a certain measurement point are determined as (x, y, z), the rock density of this measurement point is determined as D(x, y, z), and the vertical stress σv(x, y, z) of a certain measurement point can be obtained in real time through in-situ in-situ stress testing. Since there is a correlation between the vertical stress of a certain measurement point, the wave velocity value of this measurement point, and the core density value of this measurement point, the correlation expression can be expressed as:

[0114] σ v (x, y, z) = D(x, y, z) · v(x, y, z)

[0115] Where: σv(x, y, z) is the measured vertical stress value at the measuring point location; D(x, y, z) is the core density value taken at the measuring point location; v(x, y, z) is the wave velocity parameter to which the block belongs, which can be obtained by corresponding to the seismic wave CT detection cloud map and the wave velocity gradient.

[0116] For example, as Figure 3 、 Figure 5 shown, the wave velocity thresholds in four regions of the wave velocity and wave velocity gradient distribution cloud map obtained from the original seismic wave CT detection are defined as v1, v2, v3, and v4 respectively. The region between two adjacent wave velocity gradient boundaries 6 is a type of region. The vertical stress σv(x, y, z) and the corresponding density D(x, y, z) of the local measuring point 8 are known. Based on the conversion relationship between the above σv(x, y, z), D(x, y, z), and v(x, y, z), the corresponding relationship between the vertical stress σv(x, y, z) of each local measuring point 8 and the wave velocity v(x, y, z) corresponding to the position of the block 7 to which it belongs is calculated. By replacing the wave velocity in the wave velocity distribution cloud map of the existing entire region, an isogram or distribution cloud map of the absolute value of the regional vertical stress can be obtained. Optionally, (x, y, z) can be the coordinates on the working face coordinate system 5.

[0117] Step 104, based on the vertical stress distribution information of each of the multiple target detection surfaces, determine the vertical stress distribution information of the target area.

[0118] In some embodiments of the present application, step 104 includes:

[0119] Step e1, based on the three-dimensional coordinates of each of the multiple measuring points, divide the measuring points with the same z-axis coordinate value into one group.

[0120] Step e2, for each group, based on the vertical stress data of each of the multiple measuring points, respectively determine the vertical stress difference between every two adjacent measuring points.

[0121] Step e3, based on the vertical stress difference between every two adjacent measuring points, respectively calculate the vertical stress change rate of each measuring point in the z-axis direction.

[0122] For example, as Figure 5 shown, the vertical stress change rate σ′ v (i, j + 1, m) of the block (i, j + 1, m) is calculated by the following formula:

[0123]

[0124] where d is the side length of the block.

[0125] Step e4: Integrate the vertical stress change rate of the measurement points along the z-axis direction to obtain the integration result.

[0126] Step e5: Based on the integration results of each measurement point, determine the expression of the vertical stress continuous function between each measurement point and its adjacent measurement point in the z-axis direction.

[0127] Step e6: Merge the expressions of the vertical stress continuous functions of multiple measurement points in each group to obtain the expression of the vertical stress function for each group.

[0128] It can be understood that, as Figure 6 shown, since the vertical stress directions 10 of multiple measurement points in each group are consistent and the expressions of their respective vertical stress continuous functions are continuous and on the same straight line, they can be directly merged to obtain the expression of the vertical stress function for each group.

[0129] It should be noted that Figure 6 the schematic diagram of the connection of the vertical stress directions shown can be Figure 4 the schematic cross-sectional view at the vertical profile line 21 in

[0130] Step e7: Based on the expression of the vertical stress function for each group, determine the vertical stress distribution information of the target area.

[0131] As a possible implementation example, the set of the expressions of the vertical stress functions for all groups is the vertical stress distribution information of the target area.

[0132] According to the method for determining the regional vertical stress based on point-plane in-situ data fusion in the embodiments of the present application, multiple target stratifications are determined within the target area of the coal mining face according to preset rules, multiple target detection planes are determined based on the multiple target stratifications, and multiple measurement points and the three-dimensional coordinates of each measurement point are respectively determined on each target detection plane. The wave velocity distribution information of the target area, the coal and rock mass density data of each measurement point, the wave velocity data of each measurement point, and the vertical stress data of each measurement point are respectively obtained. Based on the corresponding relationships between the coal and rock mass density data of each measurement point, the corresponding wave velocity data of each measurement point, the vertical stress data of each measurement point, and the wave velocity gradient of the target area, the vertical stress distribution information of each target detection plane is respectively determined. Based on the vertical stress distribution information of each target detection plane, the vertical stress distribution information of the target area is determined. Thus, the quantitative evaluation of the regional vertical stress is realized, providing reliable data support for the mine production design intuitively, further significantly improving the mine design and production efficiency, greatly reducing the risk of disasters, and avoiding the generation of ineffective workload.

[0133] Figure 7It is a flowchart of a device for determining regional vertical stress based on point-plane in-situ data fusion in an embodiment of this application.

[0134] As Figure 7 shown, the device for determining regional vertical stress based on point-plane in-situ data fusion includes:

[0135] The first determination module 701 is configured to determine a plurality of target layers in a target area of a coal mining face according to a preset rule, determine a plurality of target detection planes based on the plurality of target layers, and respectively determine a plurality of measurement points and the three-dimensional coordinates of each of the plurality of measurement points on each target detection plane;

[0136] The acquisition module 702 is configured to respectively acquire the wave velocity distribution information of the target area, the coal and rock mass density data of each of the plurality of measurement points, the wave velocity data of each of the plurality of measurement points, and the vertical stress data of each of the plurality of measurement points; wherein, the wave velocity distribution information includes a wave velocity distribution nephogram and wave velocity gradient information;

[0137] The second determination module 703 is configured to respectively determine the vertical stress distribution information of each of the plurality of target detection planes based on the corresponding relationship between the coal and rock mass density data of each of the plurality of measurement points, the wave velocity data corresponding to each of the plurality of measurement points, the vertical stress data of each of the plurality of measurement points, and the wave velocity gradient of the target area; wherein, the vertical stress distribution information includes a vertical stress distribution nephogram and vertical stress information;

[0138] The third determination module 704 is configured to determine the vertical stress distribution information of the target area based on the vertical stress distribution information of each of the plurality of target detection planes.

[0139] According to the device for determining regional vertical stress based on point-plane in-situ data fusion in an embodiment of this application, by determining a plurality of target layers in a target area of a coal mining face according to a preset rule, determining a plurality of target detection planes based on the plurality of target layers, and respectively determining a plurality of measurement points and the three-dimensional coordinates of each of the plurality of measurement points on each target detection plane. Respectively acquire the wave velocity distribution information of the target area, the coal and rock mass density data of each of the plurality of measurement points, the wave velocity data of each of the plurality of measurement points, and the vertical stress data of each of the plurality of measurement points. Based on the corresponding relationship between the coal and rock mass density data of each of the plurality of measurement points, the wave velocity data corresponding to each of the plurality of measurement points, the vertical stress data of each of the plurality of measurement points, and the wave velocity gradient of the target area, respectively determine the vertical stress distribution information of each of the plurality of target detection planes. Based on the vertical stress distribution information of each of the plurality of target detection planes, determine the vertical stress distribution information of the target area. Thus, the quantitative evaluation of regional vertical stress is realized, which intuitively provides reliable data support for mine production design, further significantly improves the mine design and production efficiency, greatly reduces the risk of disasters, and avoids generating ineffective workload.

[0140] Figure 8The block diagram of an electronic device in an embodiment of the present application. As Figure 8 shown, the electronic device may include: a transceiver 81, a processor 82, and a memory 83.

[0141] The processor 82 executes the computer-executable instructions stored in the memory, enabling the processor 82 to execute the solutions in the above embodiments. The processor 82 may be a general-purpose processor, including a central processing unit CPU, a network processor (NP), etc.; it may also be a digital signal processor DSP, an application-specific integrated circuit ASIC, a field-programmable gate array FPGA, or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.

[0142] The memory 83 is connected to the processor 82 through a system bus and completes communication therebetween. The memory 83 is used to store computer program instructions.

[0143] The transceiver 81 may be used to obtain the task to be run and the configuration information of the task to be run.

[0144] The system bus may be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus, etc. The system bus may be divided into an address bus, a data bus, a control bus, etc. For the sake of simplicity, only a thick line is used to represent it in the figure, but it does not mean that there is only one bus or one type of bus. The transceiver is used to implement communication between the database access device and other computers (such as clients, read-write libraries, and read-only libraries). The memory may include a random access memory (RAM), and may also include a non-volatile memory.

[0145] The electronic device provided by the embodiment of the present application may be the terminal device in the above embodiment.

[0146] The embodiment of the present application further provides a chip for running instructions. The chip is used to execute the technical solution of the message processing method in the above embodiment.

[0147] The embodiment of the present application further provides a computer-readable storage medium. The computer-readable storage medium stores computer instructions. When the computer instructions are run on a computer, the computer is enabled to execute the technical solution of the message processing method in the above embodiment.

[0148] An embodiment of the present application further provides a computer program product, which includes a computer program stored in a computer-readable storage medium. At least one processor can read the computer program from the computer-readable storage medium, and when the at least one processor executes the computer program, the technical solution of the message processing method in the above embodiment can be implemented.

[0149] After considering the specification and practicing the invention disclosed herein, those skilled in the art will readily conceive of other embodiments of the present application. The present application is intended to cover any variations, uses, or adaptations of the present application, which follow the general principles of the present application and include common general knowledge or conventional technical means in the technical field not disclosed in the present application. The specification and embodiments are only regarded as exemplary, and the true scope and spirit of the present application are pointed out by the following claims.

[0150] It should be understood that the present application is not limited to the exact structures described above and shown in the drawings, and various modifications and changes can be made without departing from its scope. The scope of the present application is only limited by the appended claims.

Claims

1. A method for determining regional vertical stress based on point-surface in-situ data fusion, characterized in that The method includes: Determine a plurality of target layers within a target area of a coal mining face according to a preset rule, determine a plurality of target detection surfaces based on the plurality of target layers, and respectively determine a plurality of measurement points and the three-dimensional coordinates of each of the plurality of measurement points on each target detection surface; Obtain the wave velocity distribution information of the target area, the coal and rock mass density data of each of the plurality of measurement points, the wave velocity data of each of the plurality of measurement points, and the vertical stress data of each of the plurality of measurement points respectively; wherein, the wave velocity distribution information includes a wave velocity distribution nephogram and wave velocity gradient information; Based on the corresponding relationships between the coal and rock mass density data of each of the plurality of measurement points, the wave velocity data corresponding to each of the plurality of measurement points, the vertical stress data of each of the plurality of measurement points, and the wave velocity gradient of the target area, determine the vertical stress distribution information of each of the plurality of target detection surfaces respectively; wherein, the vertical stress distribution information includes a vertical stress distribution nephogram and vertical stress information; Based on the vertical stress distribution information of each of the plurality of target detection surfaces, determine the vertical stress distribution information of the target area.

2. The method according to claim 1, wherein The step of determining the vertical stress distribution information of each of the plurality of target detection surfaces respectively based on the corresponding relationships between the coal and rock mass density data of each of the plurality of measurement points, the wave velocity data corresponding to each of the plurality of measurement points, the vertical stress data of each of the plurality of measurement points, and the wave velocity gradient of the target area includes: For each target layer, based on the corresponding relationships between the coal and rock mass density data of the plurality of measurement points in the target layer, the wave velocity data corresponding to the plurality of measurement points, the vertical stress data of the plurality of measurement points, and the wave velocity gradient of the target layer, determine the vertical stress data corresponding to each wave velocity data; Based on the vertical stress data corresponding to each wave velocity data, replace the wave velocity distribution data in the wave velocity distribution nephogram with vertical stress distribution data to obtain the vertical stress distribution nephogram of the target layer.

3. The method according to claim 1, characterized in that, The step of determining a plurality of target layers within a target area of a coal mining face according to a preset rule and determining a plurality of target detection surfaces based on the plurality of target layers includes: Obtain the thickness of each target layer respectively; Based on the thickness of each target layer, determine the middle position of each target layer as a target detection surface respectively; Determine the interface between every two adjacent target layers as a target detection surface respectively; Based on the thickness of each target layer, determine the z-axis coordinate of each of the target detection surfaces respectively.

4. The method according to claim 3, wherein The step of obtaining the thickness of each target layer respectively includes: Obtain the total thickness of the coal seam in the target area; Obtain the preset number of target layers; Divide the total thickness of the coal seam by the preset number of target layers to obtain the coal seam thickness of each target layer.

5. The method according to claim 1, wherein The wave velocity distribution information of the target area is obtained by the following method: Determine an excitation point at the middle position of each target layer respectively; the x-axis coordinates of each of the excitation points are the same; the y-axis coordinates of each of the excitation points are the same; each of the excitation points is located on the side of the seismic wave excitation end roadway close to the seismic wave receiving end roadway; Determine a receiving point at the middle position of each target layer; the x-axis coordinates of each of the receiving points are the same; the y-axis coordinates of each of the receiving points are the same; each of the receiving points is located on the side of the roadway of the seismic wave receiving end close to the working face; Obtain the seismic wave data transmitted by the receiving end geophones corresponding to each target layer respectively; Based on the seismic wave data corresponding to each target layer, determine the wave velocity distribution information of the target area.

6. The method according to claim 1, wherein The determining the vertical stress distribution information of the target area based on the vertical stress distribution information of each of the multiple target detection surfaces further includes: Based on the three-dimensional coordinates of each of the multiple measurement points, divide the measurement points with the same z-axis coordinate value into a group; For each group, based on the vertical stress data of each of the multiple measurement points, determine the vertical stress difference between every two adjacent measurement points respectively; Based on the vertical stress differences between every two adjacent measurement points, calculate the vertical stress change rate of each measurement point in the z-axis direction respectively; Perform integral processing on the vertical stress change rate of the measurement points in the z-axis direction to obtain an integral result; Based on the integral result of each measurement point, determine the vertical stress continuous function expression between each measurement point and its adjacent measurement point in the z-axis direction; Merge the vertical stress continuous function expressions of the multiple measurement points in each group to obtain the vertical stress function expression of each group; Based on the vertical stress function expression of each group, determine the vertical stress distribution information of the target area.

7. The method according to claim 1, wherein The determining multiple measurement points and the three-dimensional coordinates of each of the multiple measurement points on each target detection surface respectively includes: For each target detection surface, evenly divide the target detection surface into multiple blocks, and determine the coordinate regions of each of the multiple blocks respectively; Compare the coordinate regions of each of the multiple blocks with a preset region, and in response to the coordinate region falling into the preset coordinate region, remove the block corresponding to the coordinate region to obtain the multiple blocks; wherein, the preset coordinate region is a wave velocity abnormal region; Based on the coordinate regions of each of the multiple blocks, determine the center point coordinates of each of the multiple blocks respectively; Determine the points corresponding to the center point coordinates of each of the multiple blocks as the measurement points of the target detection surface.

8. A device for determining regional vertical stress based on point-plane in-situ data fusion, characterized in that, The device includes: A first determination module, configured to determine multiple target layers in a target area of a coal mine working face according to a preset rule, determine multiple target detection surfaces based on the multiple target layers, and determine multiple measurement points and the three-dimensional coordinates of each of the multiple measurement points on each target detection surface respectively; An acquisition module, configured to acquire the wave velocity distribution information of the target area, the coal and rock mass density data of each of the multiple measurement points, the wave velocity data of each of the multiple measurement points, and the vertical stress data of each of the multiple measurement points respectively; wherein, the wave velocity distribution information includes a wave velocity distribution nephogram and wave velocity gradient information; A second determination module, configured to respectively determine the vertical stress distribution information of each of the plurality of target detection surfaces based on the corresponding relationship between the coal and rock mass density data of each of the plurality of measurement points, the wave velocity data corresponding to each of the plurality of measurement points, the vertical stress data of each of the plurality of measurement points, and the wave velocity gradient of the target area; wherein the vertical stress distribution information includes a vertical stress distribution nephogram and vertical stress information. A third determination module, configured to determine the vertical stress distribution information of the target area based on the vertical stress distribution information of each of the plurality of target detection surfaces.

9. An electronic device, characterized in that, Comprising: A processor, and a memory communicatively connected to the processor; The memory stores computer-executable instructions; The processor executes the computer-executable instructions stored in the memory to implement the method according to any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, Computer-executable instructions are stored in the computer-readable storage medium, and when the computer-executable instructions are executed by the processor, they are used to implement the method according to any one of claims 1-7.

Images