Method and device for determining regional horizontal stress based on point-surface in-situ data fusion

Abstract

The application relates to a regional horizontal stress determination method and device based on point-surface in-situ data fusion. The specific scheme is as follows: N measuring points and the coordinate values of the N measuring points in a target region of a coal mine working face are determined according to a preset rule; the maximum horizontal principal stress value and the minimum horizontal principal stress value of each of the N measuring points are obtained; a measuring point is determined as a center measuring point, and at least one adjacent measuring point of the center measuring point is determined; the maximum horizontal principal stress migration direction of the center measuring point and the minimum horizontal principal stress migration direction of the center measuring point are respectively determined; the maximum horizontal principal stress change rate and the minimum horizontal principal stress change rate between any two adjacent measuring points are respectively determined; and the maximum horizontal principal stress vector distribution data and the minimum horizontal principal stress vector distribution data of the target region are respectively determined. The application can determine the continuous horizontal stress distribution of a coal mining region.

Description

Technical Field

[0001] This application relates to the field of geostress measurement technology in underground coal mine engineering, and in particular to a method and apparatus for determining regional horizontal stress based on point-surface in-situ data fusion. Background Technology

[0002] In related technologies, during underground coal mining, the regional mining stress generated by roadway excavation and face mining significantly impacts roadway stability. In particular, abnormal concentrations of horizontal stress can easily lead to severe dynamic disasters such as rockbursts and coal and gas outbursts. Accurately determining the direction and absolute value of horizontal stress distribution in the mining area through in-situ testing is crucial; however, there is currently no reliable in-situ measurement method for the absolute value of regional horizontal stress. Water-induced fracturing and stress relief methods, as relatively reliable in-situ measurement methods for localized in-situ stress, are widely used in production. Typically, the stress results from one or two local measuring points are selected to guide the design and production work of the entire mining area or face. However, in actual engineering, the physical properties of coal and rock masses exhibit significant variations, and are also affected by coal seam undulations and geological structures. Localized stress testing at select locations is insufficient to meet engineering needs. Seismic wave CT technology can detect large areas of the working face. Its principle is to reflect the relative vertical stress accumulation distribution in the area by observing the changes in the propagation velocity of seismic waves in the coal and rock mass, but it cannot obtain the absolute value of regional horizontal stress. Summary of the Invention

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

[0004] According to a first aspect of the embodiments of this application, a method for determining regional horizontal stress based on point-surface in-situ data fusion is provided, wherein the stress includes a maximum horizontal principal stress and a minimum horizontal principal stress, the method comprising:

[0005] According to preset rules, N measuring points and their respective coordinate values ​​are determined within the target area of ​​the coal mine working face; wherein, N is an integer greater than or equal to 9;

[0006] Obtain the maximum and minimum horizontal principal stress values ​​for each of the N measuring points;

[0007] For each measuring point, the measuring point is determined as the central measuring point, and at least one adjacent measuring point of the central measuring point is determined;

[0008] Based on the coordinates of the central measuring point, the maximum horizontal principal stress of the central measuring point, and the minimum horizontal principal stress of the central measuring point, as well as the coordinates of each of the adjacent measuring points, the maximum horizontal principal stress of each of the adjacent measuring points, and the minimum horizontal principal stress of each of the adjacent measuring points, the migration direction of the maximum horizontal principal stress of the central measuring point and the migration direction of the minimum horizontal principal stress of the central measuring point are determined respectively.

[0009] Based on the maximum and minimum horizontal principal stress values ​​of each of the N measuring points, the maximum and minimum horizontal principal stress change rates between any two adjacent measuring points are determined respectively.

[0010] Based on the migration direction of the maximum and minimum horizontal principal stresses at each measuring point, as well as the rate of change of the maximum and minimum horizontal principal stresses between any two adjacent measuring points, the maximum and minimum horizontal principal stress vector distribution data of the target area are determined respectively.

[0011] According to one embodiment of this application, determining the maximum horizontal principal stress vector distribution data and the minimum horizontal principal stress vector distribution data of the target area based on the maximum horizontal principal stress migration direction and the minimum horizontal principal stress migration direction of each measuring point, as well as the rate of change of the maximum horizontal principal stress and the rate of change of the minimum horizontal principal stress between any two adjacent measuring points, includes:

[0012] Connect the head and tail of the maximum horizontal principal stress migration direction of each of the N measuring points to obtain the maximum horizontal principal stress migration direction of the target area.

[0013] Based on the migration direction of the maximum horizontal principal stress in the target area and the rate of change of the maximum horizontal principal stress between any two adjacent measuring points, the maximum horizontal principal stress vector distribution data of the target area is determined.

[0014] Connect the beginning and end of the minimum horizontal principal stress migration direction of each of the N measuring points to obtain the minimum horizontal principal stress migration direction of the target area.

[0015] Based on the migration direction of the minimum horizontal principal stress in the target area and the rate of change of the minimum horizontal principal stress between any two adjacent measuring points, the minimum horizontal principal stress vector distribution data of the target area are determined.

[0016] According to one embodiment of this application, determining the migration direction of the maximum horizontal principal stress and the migration direction of the minimum horizontal principal stress of the central measuring point based on the coordinate values ​​of the central measuring point, the maximum horizontal principal stress value of the central measuring point, and the minimum horizontal principal stress value of the central measuring point, as well as the coordinate values ​​of each of the adjacent measuring points, the maximum horizontal principal stress value of each of the adjacent measuring points, and the minimum horizontal principal stress value of each of the adjacent measuring points, respectively, includes:

[0017] Based on the coordinates of the central measuring point, the maximum horizontal principal stress of the central measuring point, the coordinates of each of the adjacent measuring points, and the maximum horizontal principal stress of each of the adjacent measuring points, the directional derivative is calculated to obtain the migration direction of the maximum horizontal principal stress between the central measuring point and each of the adjacent measuring points.

[0018] The maximum value among the maximum horizontal principal stress migration directions between the central measuring point and at least one adjacent measuring point around it is determined as the maximum horizontal principal stress migration direction of the central measuring point.

[0019] Based on the coordinates of the central measuring point, the maximum horizontal principal stress of the central measuring point, the coordinates of each of the adjacent measuring points, and the maximum horizontal principal stress of each of the adjacent measuring points, the directional derivative is calculated to obtain the migration direction of the maximum horizontal principal stress between the central measuring point and each of the adjacent measuring points.

[0020] The maximum value among the minimum horizontal principal stress migration directions between the central measuring point and the at least one adjacent measuring point is determined as the minimum horizontal principal stress migration direction of the central measuring point.

[0021] According to one embodiment of this application, after determining the migration direction of the maximum horizontal principal stress and the migration direction of the minimum horizontal principal stress of the central measuring point based on the coordinate value of the central measuring point, the maximum horizontal principal stress value of the central measuring point, and the minimum horizontal principal stress value of the central measuring point, as well as the coordinate value of each of the adjacent measuring points, the maximum horizontal principal stress value of each of the adjacent measuring points, and the minimum horizontal principal stress value of each of the adjacent measuring points, the method further includes:

[0022] For each central measuring point, the maximum rate of change of the horizontal principal stress at the central measuring point in the direction of maximum horizontal principal stress migration and the minimum rate of change of the horizontal principal stress at the central measuring point in the direction of minimum horizontal principal stress migration are determined respectively.

[0023] According to a preset interpolation factor k, based on the maximum horizontal principal stress change rate of the N measuring points in their respective maximum horizontal principal stress migration directions and the minimum horizontal principal stress change rate of the N measuring points in their respective minimum horizontal principal stress migration directions, interpolation processing is performed on the maximum and minimum horizontal principal stress values ​​of the N measuring points respectively to obtain the maximum and minimum horizontal principal stress values ​​of kN measuring points respectively; where k is an integer greater than 0;

[0024] The maximum horizontal principal stress migration direction and the minimum horizontal principal stress migration direction of each of the N measuring points are evenly divided to obtain the maximum horizontal principal stress migration direction and the minimum horizontal principal stress migration direction of each of the kN measuring points after the even division.

[0025] According to one embodiment of this application, after determining the migration direction of the maximum horizontal principal stress and the migration direction of the minimum horizontal principal stress of the central measuring point based on the coordinate value of the central measuring point, the maximum horizontal principal stress value of the central measuring point, and the minimum horizontal principal stress value of the central measuring point, as well as the coordinate value of each of the adjacent measuring points, the maximum horizontal principal stress value of each of the adjacent measuring points, and the minimum horizontal principal stress value of each of the adjacent measuring points, the method further includes:

[0026] The maximum horizontal principal stress variation rate at the central measuring point along the x-direction and the maximum horizontal principal stress variation rate at the central measuring point along the y-direction are determined; wherein, the maximum horizontal principal stress variation rate at the central measuring point along the x-direction and the maximum horizontal principal stress variation rate at the central measuring point along the y-direction are both obtained by decomposing the maximum horizontal principal stress variation rate between the central measuring point and the adjacent measuring point in the direction of migration of the maximum horizontal principal stress along the x and y directions, respectively.

[0027] The rate of change of the maximum horizontal principal stress at the central measuring point along the x-direction is integrated to obtain a first integral result; the rate of change of the maximum horizontal principal stress at the central measuring point along the y-direction is integrated to obtain a second integral result.

[0028] Based on the first integration result and the second integration result, the expression for the maximum horizontal stress migration continuity function between the central measuring point and its adjacent measuring points in the direction of maximum horizontal stress migration is determined.

[0029] Determine the minimum rate of change of the horizontal principal stress at the central measuring point along the x-direction and the minimum rate of change of the horizontal principal stress at the central measuring point along the y-direction;

[0030] The rate of change of the minimum horizontal principal stress at the central measuring point along the x-direction is integrated to obtain the third integral result; the rate of change of the minimum horizontal principal stress at the central measuring point along the y-direction is integrated to obtain the fourth integral result.

[0031] Based on the third and fourth integral results, the expression for the minimum horizontal stress migration continuity function between the central measuring point and its adjacent measuring points in the direction of minimum horizontal stress migration is determined.

[0032] According to one embodiment of this application, the determination of the maximum horizontal principal stress distribution vector data and the minimum horizontal principal stress distribution data of the target area based on the maximum horizontal principal stress migration direction and the minimum horizontal principal stress migration direction of each measuring point, as well as the maximum horizontal principal stress change rate and the minimum horizontal principal stress change rate between any two adjacent measuring points, includes:

[0033] Based on the minimum horizontal stress migration continuity function expression between each measuring point and its two adjacent measuring points in the minimum horizontal stress migration direction, the first function value and the second function value of each measuring point are determined respectively.

[0034] In response to the fact that the first function value and the second function value are the same, and the continuous function expression of minimum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of minimum horizontal stress migration satisfies the function continuity condition, the continuous function expressions of minimum horizontal stress migration corresponding to the first function value and the second function value are merged to obtain the merged continuous function expression of minimum horizontal stress migration.

[0035] In response to the difference between the first function value and the second function value, and / or the continuous function expression of minimum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of minimum horizontal stress migration does not satisfy the function continuity condition, the continuous function expression of minimum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of minimum horizontal stress migration is determined as the first piecewise function;

[0036] The merged minimum horizontal stress migration continuous function expression and the first piecewise function are determined as the minimum horizontal principal stress vector distribution function of the target region;

[0037] Based on the maximum horizontal stress migration continuity function expression between each measuring point and its two adjacent measuring points in the direction of maximum horizontal stress migration, the third and fourth function values ​​of each measuring point are determined respectively;

[0038] In response to the fact that the third function value is the same as the fourth function value, and the maximum horizontal stress migration continuous function expression between the measuring point and its two adjacent measuring points in the maximum horizontal stress migration direction satisfies the function continuity condition, the maximum horizontal stress migration continuous function expressions corresponding to the third function value and the fourth function value are merged to obtain the merged maximum horizontal stress migration continuous function expression.

[0039] In response to the fact that the third function value is different from the fourth function value, and / or the continuous function expression of the maximum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of maximum horizontal stress migration does not satisfy the function continuity condition, the continuous function expression of the maximum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of maximum horizontal stress migration is determined as the second piecewise function.

[0040] The combined maximum horizontal stress migration continuous function expression and the second piecewise function are determined as the maximum horizontal principal stress vector distribution function of the target region.

[0041] According to one embodiment of this application, determining the rate of change of the maximum and minimum horizontal principal stress between any two adjacent measuring points based on the maximum and minimum horizontal principal stress values ​​of the N measuring points includes:

[0042] The maximum horizontal principal stress value at the central measuring point is subtracted from the maximum horizontal principal stress value at each of the adjacent measuring points to obtain a first difference; based on the first difference, the rate of change of the maximum horizontal principal stress at the central measuring point and at each of the adjacent measuring points is determined.

[0043] The minimum horizontal principal stress value of the central measuring point is subtracted from the minimum horizontal principal stress value of each of the adjacent measuring points to obtain a second difference; based on the second difference, the rate of change of the minimum horizontal principal stress of the central measuring point and each of the adjacent measuring points is determined.

[0044] According to one embodiment of this application, determining N measuring points within a target area of ​​the coal mine working face according to preset rules includes:

[0045] The target area is evenly divided into multiple blocks, and the coordinate regions of each of the multiple blocks are determined.

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

[0047] Based on the coordinate regions of each of the N blocks, the coordinates of the center point of each of the N blocks are determined respectively;

[0048] Based on the coordinates of the center point of each of the N blocks, N measurement points within the target area are determined.

[0049] According to a second aspect of the embodiments of this application, a device for determining regional horizontal stress based on point-surface in-situ data fusion is provided, the device comprising:

[0050] The first determining module is used to determine N measuring points and their respective coordinate values ​​within a target area of ​​a coal mine working face according to preset rules; wherein, N is an integer greater than or equal to 9;

[0051] The acquisition module is used to acquire the maximum and minimum horizontal principal stress values ​​of each of the N measuring points;

[0052] The second determining module is used to determine each measuring point as a central measuring point and to determine at least one adjacent measuring point of the central measuring point.

[0053] The third determining module is used to determine the migration direction of the maximum horizontal principal stress of the central measuring point and the migration direction of the minimum horizontal principal stress of the central measuring point based on the coordinate value of the central measuring point, the maximum horizontal principal stress value of the central measuring point and the minimum horizontal principal stress value of the central measuring point, as well as the coordinate value of each of the adjacent measuring points, the maximum horizontal principal stress value of each of the adjacent measuring points and the minimum horizontal principal stress value of each of the adjacent measuring points.

[0054] The fourth determining module is used to determine the rate of change of the maximum and minimum horizontal principal stress between any two adjacent measuring points based on the maximum and minimum horizontal principal stress values ​​of each of the N measuring points.

[0055] The fifth determining module is used to determine the maximum horizontal principal stress vector distribution data and the minimum horizontal principal stress vector distribution data of the target area based on the maximum horizontal principal stress migration direction and the minimum horizontal principal stress migration direction of each measuring point, as well as the maximum horizontal principal stress change rate and the minimum horizontal principal stress change rate between any two adjacent measuring points.

[0056] According to one embodiment of this application, the fifth determining module includes:

[0057] The first connection submodule is used to connect the head and tail of the maximum horizontal principal stress migration direction of each of the N measuring points to obtain the maximum horizontal principal stress migration direction of the target area.

[0058] The first determining submodule is used to determine the maximum horizontal principal stress vector distribution data of the target area based on the migration direction of the maximum horizontal principal stress in the target area and the rate of change of the maximum horizontal principal stress between any two adjacent measuring points.

[0059] The second connection submodule is used to connect the beginning and end of the minimum horizontal principal stress migration direction of each of the N measuring points to obtain the minimum horizontal principal stress migration direction of the target area.

[0060] The second determining submodule is used to determine the minimum horizontal principal stress vector distribution data of the target area based on the minimum horizontal principal stress migration direction of the target area and the minimum horizontal principal stress change rate between any two adjacent measuring points.

[0061] The technical solution provided by the embodiments of this application brings at least the following beneficial effects: by fusing the in-situ measured data of seismic wave CT detection and local point stress testing in the target area, the maximum and minimum horizontal principal stress vector distribution data in the target area of ​​the coal mine working face are determined, realizing the quantitative assessment of regional horizontal stress, intuitively providing reliable data support for mine production design, thereby significantly improving the efficiency of mine design and production, greatly reducing the risk of disaster occurrence, and reducing the workload of disaster prevention and mitigation.

[0062] It should be understood that the above general description and the following detailed description are exemplary and explanatory only, and do not limit this application. Attached Figure Description

[0063] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this application and, together with the description, serve to explain the principles of this application, and do not constitute an undue limitation of this application.

[0064] Figure 1 This is a flowchart of a method for determining regional horizontal stress based on point-surface in-situ data fusion in an embodiment of this application;

[0065] Figure 2 This is a structural block diagram of a regional horizontal stress determination device based on point-surface in-situ data fusion according to an embodiment of this application;

[0066] Figure 3 This refers to the wave velocity distribution cloud map and wave velocity gradient information of the seismic wave CT detection in the embodiments of this application;

[0067] Figure 4 This is a schematic diagram showing the location of blocks and the arrangement of measuring points in an embodiment of this application.

[0068] Figure 5 This is a schematic diagram of the stress migration direction at the target measuring point in an embodiment of this application;

[0069] Figure 6 This is a schematic diagram of the distribution of the regional horizontal stress migration direction in an embodiment of this application.

[0070] Figure Labels

[0071] 4- Wave velocity gradient anomaly region; 5- Working surface coordinate system; 6- Wave velocity gradient boundary line; 7- Block; 8- Measuring point; 9- Central measuring point; 10- Stress direction of measuring point; 11- Horizontal stress migration direction. Detailed Implementation

[0072] To enable those skilled in the art to better understand the technical solutions of this application, the technical solutions in the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings.

[0073] It should be noted that the terms "first," "second," etc., used in the specification, claims, and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of this application described herein can be implemented in orders other than those illustrated or described herein. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with this application. Rather, they are merely examples of apparatuses and methods consistent with some aspects of this application as detailed in the appended claims.

[0074] It should be noted that in underground coal mining engineering, the regional mining stress generated by roadway excavation and working face mining has a significant impact on roadway stability. In particular, abnormal concentration of horizontal stress can easily lead to serious dynamic disasters such as coal mine rock bursts and coal and gas outbursts. Accurately determining the distribution direction and absolute value of horizontal stress in the mining area through in-situ testing is crucial, but currently there is no reliable in-situ measurement method for the absolute value of regional horizontal stress. Water pressure fracturing and stress relief are widely used as relatively reliable in-situ measurement methods for localized in-situ stress in production lines. Generally, the stress results from one or two local measuring points are used to guide the design and production work of the entire mining area or working face. However, in actual engineering, the physical properties of coal and rock masses vary significantly, and are also affected by coal seam undulations and geological structures. Localized stress testing at individual locations is insufficient to meet engineering needs. Seismic wave CT technology can achieve large-scale area detection at the working face. Its principle is to reflect the relative vertical stress accumulation distribution in the area by the variation in the propagation speed of seismic waves in the coal and rock mass, but it cannot obtain the absolute value of regional horizontal stress.

[0075] To address the aforementioned issues, this application proposes a method and apparatus for determining regional horizontal stress based on point-surface in-situ data fusion. This method can determine the maximum and minimum horizontal principal stress vector distribution data within the target area of ​​a coal mine working face by fusing in-situ measured data from seismic CT detection and local point stress testing in the target area. This enables quantitative assessment of regional horizontal stress, providing reliable data support for mine production design, thereby significantly improving mine design and production efficiency, greatly reducing the risk of disasters, and reducing the workload of disaster prevention and mitigation.

[0076] Figure 1 This is a flowchart of a method for determining regional horizontal stress based on point-surface in-situ data fusion in an embodiment of this application.

[0077] like Figure 1 As shown, the method for determining regional horizontal stress based on point-surface in-situ data fusion includes:

[0078] Step 101: Determine N measuring points and their respective coordinate values ​​within the target area of ​​the coal mine working face according to preset rules.

[0079] In some embodiments of this application, N is an integer greater than or equal to 9.

[0080] In some embodiments of this application, step 101 includes:

[0081] Step a1: Divide the target area evenly into multiple blocks and determine the coordinate regions of each block.

[0082] It should be noted that the block can cover all areas except for the abnormal area as much as possible.

[0083] Step a2: Compare the coordinate regions of each of the multiple blocks with the preset region. In response to the coordinate region falling into the preset coordinate region, remove the block corresponding to that coordinate region to obtain N blocks.

[0084] In this embodiment of the application, the preset coordinate region is a wave velocity anomaly region.

[0085] It should be noted that, as Figure 3 , Figure 4As shown, when selecting the location of measuring point block 7, it is necessary to avoid wave velocity anomaly area 4, microseismic event clusters, and areas with well-developed local geological structures in the wave velocity and wave velocity gradient cloud map, thereby enhancing the reliability and representativeness of the regional test results. It can be understood that the four types of regions in the wave velocity and wave velocity gradient distribution cloud map obtained from the original seismic wave CT detection are defined with wave velocity thresholds of v1, v2, v3, and v4, respectively, and the region between two adjacent wave velocity gradient boundaries 6 is considered as one type of region. Optionally, the coordinates (x, y) of measuring point 8 can be coordinates on the working surface coordinate system 5.

[0086] As one possible implementation example, the coordinate regions of each of the multiple blocks are compared with a preset region. If a coordinate region falls within the preset coordinate region, it indicates that the block does not meet the requirements and is removed. After filtering multiple blocks using the above method, the remaining blocks are the aforementioned N blocks.

[0087] Optionally, the aforementioned blocks can be squares with a preset side length, and the preset value should be as small as possible to obtain as many blocks as possible.

[0088] Step a3: Based on the coordinate regions of each of the N blocks, determine the coordinates of the center point of each of the N blocks.

[0089] Step a4: Based on the coordinates of the center points of the N blocks, determine the N measurement points within the target area.

[0090] As an example of possible implementation, based on the coordinate regions of each of the N blocks, the coordinates of the center point of each of the N blocks are determined, and the points corresponding to the coordinates of the center point of each of the N blocks are determined as N measurement points in the target area, and each measurement point represents the block to which it belongs.

[0091] Step 102: Obtain the maximum and minimum horizontal principal stress values ​​for each of the N measuring points.

[0092] Optionally, boreholes can be drilled at selected measuring points, and then in-situ geostress tests can be performed to obtain the maximum and minimum horizontal principal stress values ​​for each of the N measuring points.

[0093] Step 103: For each measuring point, determine the measuring point as the central measuring point, and determine at least one adjacent measuring point of the central measuring point.

[0094] Step 104: Based on the coordinates of the central measuring point, the maximum horizontal principal stress of the central measuring point, and the minimum horizontal principal stress of the central measuring point, as well as the coordinates of each adjacent measuring point, the maximum horizontal principal stress of each adjacent measuring point, and the minimum horizontal principal stress of each adjacent measuring point, determine the migration direction of the maximum horizontal principal stress of the central measuring point and the migration direction of the minimum horizontal principal stress of the central measuring point, respectively.

[0095] In some embodiments of this application, step 104 further includes:

[0096] Step b1: Based on the coordinates of the central measuring point, the maximum horizontal principal stress of the central measuring point, the coordinates of each adjacent measuring point, and the maximum horizontal principal stress of each adjacent measuring point, determine the migration direction of the maximum horizontal principal stress between the central measuring point and each adjacent measuring point.

[0097] In some embodiments of this application, step b1 specifically includes: calculating the directional derivative based on the coordinates of the center measuring point, the maximum horizontal principal stress of the center measuring point, the coordinates of each adjacent measuring point, and the maximum horizontal principal stress of each adjacent measuring point, to obtain the stress vector change rate between the center measuring point and each adjacent measuring point.

[0098] Step b2: The migration direction corresponding to the maximum value of the stress vector change rate between the central measuring point and at least one adjacent measuring point is determined as the migration direction of the maximum horizontal principal stress of the central measuring point.

[0099] As an example of a possible implementation, the coordinates of the central measuring point are (i, j), the coordinates of a measuring point adjacent to the central block are (x, y), and the directional derivative of the maximum horizontal principal stress migration at the central measuring point is σ. H (i, j) is calculated using the following formula:

[0100]

[0101] Where d is the length of a single block; x and y are the coordinates of adjacent blocks, respectively.

[0102] Understandably, since the central measuring point is used to represent the block it belongs to, d is the length of a single block.

[0103] like Figure 5 As shown, the coordinates of the central measuring point 9 are (i, j). There are 8 adjacent measuring points. The coordinate values ​​of the 8 adjacent measuring points are assigned to (x, y), and the derivative of the stress at the central measuring point with respect to each adjacent measuring point is calculated along the direction 10. The maximum value is then determined to obtain the result.

[0104] Step b3: Based on the coordinates of the central measuring point, the minimum horizontal principal stress of the central measuring point, the coordinates of each adjacent measuring point, and the minimum horizontal principal stress of each adjacent measuring point, determine the direction of minimum horizontal principal stress migration between the central measuring point and each adjacent measuring point.

[0105] In some embodiments of this application, step b3 specifically includes: calculating the directional derivative based on the coordinates of the center measuring point, the minimum horizontal principal stress of the center measuring point, the coordinates of each adjacent measuring point, and the minimum horizontal principal stress of each adjacent measuring point, to obtain the migration direction of the minimum horizontal principal stress between the center measuring point and each adjacent measuring point.

[0106] Step b4: Determine the maximum value among the minimum horizontal principal stress migration directions between the center measuring point and at least one adjacent measuring point as the minimum horizontal principal stress migration direction of the center measuring point.

[0107] As an example of a possible implementation, the coordinates of the central measuring point are (i, j), the coordinates of a measuring point adjacent to the central block are (x, y), and the direction of migration of the minimum horizontal principal stress at the central measuring point is... It is calculated using the following formula:

[0108]

[0109] Where d is the length of a single block; x and y are the coordinates of adjacent blocks, respectively.

[0110] In some embodiments of this application, after step 104, the method further includes:

[0111] Step c1: For each central measuring point, determine the maximum rate of change of the horizontal principal stress in the direction of maximum horizontal principal stress migration and the minimum rate of change of the horizontal principal stress in the direction of minimum horizontal principal stress migration.

[0112] Step c2: According to the preset interpolation factor k, based on the maximum horizontal principal stress change rate of N measuring points in their respective maximum horizontal principal stress migration direction and the minimum horizontal principal stress change rate of N measuring points in their respective minimum horizontal principal stress migration direction, interpolate the maximum horizontal principal stress value and the minimum horizontal principal stress value of each of the N measuring points to obtain the maximum horizontal principal stress value and the minimum horizontal principal stress value of each of the kN measuring points.

[0113] In this embodiment of the application, k is an integer greater than 0.

[0114] As an example of a possible implementation, the above interpolation method could be a numerical interpolation method such as Kriging interpolation.

[0115] Step c3: Divide the maximum horizontal principal stress migration direction and the minimum horizontal principal stress migration direction of each of the N measuring points into equal parts to obtain the maximum horizontal principal stress migration direction and the minimum horizontal principal stress migration direction of each of the kN measuring points after equal parts processing.

[0116] As a possible implementation example, the maximum and minimum horizontal principal stress migration directions at N measuring points are equally divided using the equal-division method, thus obtaining the maximum and minimum horizontal principal stress migration directions at kN measuring points. Taking the maximum horizontal principal stress as an example, the equal-division formula is as follows:

[0117] Step 105: Based on the maximum and minimum horizontal principal stress values ​​of each of the N measuring points, determine the rate of change of the maximum and minimum horizontal principal stress between any two adjacent measuring points.

[0118] In some embodiments of this application, step 105 includes:

[0119] The maximum horizontal principal stress value at the center measuring point is subtracted from the maximum horizontal principal stress value at each adjacent measuring point to obtain the first difference. Based on the first difference, the rate of change of the maximum horizontal principal stress at the center measuring point and at each adjacent measuring point is determined. The minimum horizontal principal stress value at the center measuring point is subtracted from the minimum horizontal principal stress value at each adjacent measuring point to obtain the second difference. Based on the second difference, the rate of change of the minimum horizontal principal stress at the center measuring point and at each adjacent measuring point is determined.

[0120] In some embodiments of this application, after step 105, the method further includes:

[0121] Step d1: Determine the maximum rate of change of the horizontal principal stress at the center measuring point along the x-direction and the maximum rate of change of the horizontal principal stress at the center measuring point along the y-direction.

[0122] In this embodiment, the maximum horizontal principal stress change rate of the center measuring point along the x-direction and the maximum horizontal principal stress change rate of the center measuring point along the y-direction are obtained by decomposing the maximum horizontal principal stress change rate between the center measuring point and the adjacent measuring point in the direction of maximum horizontal principal stress migration along the x and y directions, respectively.

[0123] For example, to solve for the rate of change of the maximum horizontal stress along the stress migration direction of the block represented by each measuring point. H And the rate of change of stress σ' in the direction of stress migration H Decomposed into the rate of change σ' along the x-direction Hx and the rate of change along the y-axis σ' Hy Let's take block (i, j) as an example: Assume that the stress migration direction of block (i, j) is the same as that of block (i-1, j-1).

[0124] The expression for the rate of change of stress σ'H along the stress migration direction in block (i, j) is:

[0125]

[0126] Block (i, j) is decomposed into the stress change rate σ' along the x-direction. Hx The expression is:

[0127]

[0128] Block (i, j) is decomposed into the stress change rate σ' along the y-direction. Hy The expression is:

[0129]

[0130] Where d is the side length of the aforementioned block.

[0131] Step d2: Integrate the rate of change of the maximum horizontal principal stress at the center measuring point along the x-direction to obtain the first integral result; Integrate the rate of change of the maximum horizontal principal stress at the center measuring point along the y-direction to obtain the second integral result.

[0132] Step d3: Based on the first and second integral results, determine the expression for the continuous function of maximum horizontal stress migration between the central measuring point and its adjacent measuring points in the direction of maximum horizontal stress migration.

[0133] Step d4: Determine the minimum rate of change of the horizontal principal stress at the center measuring point along the x-direction and the minimum rate of change of the horizontal principal stress at the center measuring point along the y-direction.

[0134] Step d5: Integrate the rate of change of the minimum horizontal principal stress at the center measuring point along the x-direction to obtain the third integral result; Integrate the rate of change of the minimum horizontal principal stress at the center measuring point along the y-direction to obtain the fourth integral result.

[0135] Step d6: Based on the third and fourth integral results, determine the expression for the minimum horizontal stress migration continuity function between the central measuring point and its adjacent measuring points in the direction of minimum horizontal stress migration.

[0136] Understandably, compared to directly integrating the maximum and minimum horizontal principal stress change rates, decomposing them along the x-axis and y-axis respectively reduces computational complexity, improves efficiency, and facilitates unmanned operation and automated data processing. Using integration allows for the extraction of continuous stress migration data between two adjacent measuring points based on the stress migration data at each point.

[0137] Step 106: Based on the migration direction of the maximum and minimum horizontal principal stresses at each measuring point, as well as the rate of change of the maximum and minimum horizontal principal stresses between any two adjacent measuring points, determine the maximum and minimum horizontal principal stress vector distribution data of the target area.

[0138] In some embodiments of this application, such as Figure 6 As shown, the arrow direction indication step 106 includes:

[0139] Step e1: Connect the head and tail of the maximum horizontal principal stress migration direction 11 of each of the N measuring points 8 to obtain the maximum horizontal principal stress migration direction of the target area.

[0140] Step e2: Based on the migration direction of the maximum horizontal principal stress in the target area and the rate of change of the maximum horizontal principal stress between any two adjacent measuring points 8, determine the vector distribution data of the maximum horizontal principal stress in the target area.

[0141] Step e3: Connect the beginning and end of the minimum horizontal principal stress migration directions of the N measuring points 8 to obtain the minimum horizontal principal stress migration direction of the target area.

[0142] Step e4: Based on the migration direction of the minimum horizontal principal stress in the target area and the rate of change of the minimum horizontal principal stress between any two adjacent measuring points 8, determine the vector distribution data of the minimum horizontal principal stress in the target area.

[0143] In some embodiments of this application, step 106 includes:

[0144] Step f1: Based on the minimum horizontal stress migration continuity function expression between each measuring point and its two adjacent measuring points in the minimum horizontal stress migration direction, determine the first function value and the second function value of each measuring point.

[0145] As an example of a possible implementation, besides the measuring points at both ends of the aforementioned minimum horizontal stress migration direction, each measuring point has two adjacent measuring points along its minimum horizontal stress migration direction; that is, the connected measuring points of the i-th measuring point are the (i-1)-th and (i+1)-th measuring points. A first minimum horizontal stress migration continuity function expression is used to represent the minimum horizontal stress distribution between the i-th and (i-1)-th measuring points. Substituting the coordinates of the i-th measuring point into this expression yields the first function value for the i-th measuring point. A second minimum horizontal stress migration continuity function expression is used to represent the minimum horizontal stress distribution between the i-th and (i+1)-th measuring points. Substituting the coordinates of the i-th measuring point into this expression yields the second function value for the i-th measuring point.

[0146] Step f2: In response to the fact that the first function value and the second function value are the same, and the continuous function expression of minimum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of minimum horizontal stress migration satisfies the function continuity condition, the continuous function expressions of minimum horizontal stress migration corresponding to the first function value and the second function value are merged to obtain the merged continuous function expression of minimum horizontal stress migration.

[0147] Step f3: In response to the difference between the first function value and the second function value, and / or the continuous function expression of minimum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of minimum horizontal stress migration does not satisfy the function continuity condition, the continuous function expression of minimum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of minimum horizontal stress migration is determined as the first piecewise function.

[0148] Step f4: The merged minimum horizontal stress migration continuous function expression and the first piecewise function are determined as the minimum horizontal principal stress vector distribution function of the target region.

[0149] For example, consider the pairwise connection of three adjacent blocks i-1, i, and i+1. The stress function f(x, y) between blocks i-1 and i (i-1≤x<i, i-1≤y<i) and the stress function g(x, y) between blocks i and i+1 (i≤x<i+1, i≤y<i+1) have the same value f(x0, y0) and g(x0, y0) at the connection point (x0, y0) (i.e., the measurement point corresponding to block i). If both f(x0, y0) and g(x0, y0) are the same and satisfy the function continuity condition, they can be directly merged into a new function. If the function values ​​f(x0, y0) and g(x0, y0) at the connection point (x0, y0) are different or do not satisfy the function continuity condition, then the connection point is considered discontinuous, and f(x0, y0) and g(x0, y0) are treated as piecewise functions.

[0150] Step f5: Based on the maximum horizontal stress migration continuity function expression between each measuring point and its two adjacent measuring points in the maximum horizontal stress migration direction, determine the third function value and the fourth function value of each measuring point.

[0151] Step f6: In response to the fact that the third function value and the fourth function value are the same, and the continuous function expression of the maximum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of maximum horizontal stress migration satisfies the function continuity condition, the continuous function expressions of the maximum horizontal stress migration corresponding to the third function value and the fourth function value are merged to obtain the merged continuous function expression of the maximum horizontal stress migration.

[0152] Step f7: In response to the difference between the third function value and the fourth function value, and / or the continuous function expression of the maximum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of maximum horizontal stress migration does not satisfy the function continuity condition, the continuous function expression of the maximum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of maximum horizontal stress migration is determined as the second piecewise function.

[0153] Step f8: The combined expression of the maximum horizontal stress migration continuous function and the second piecewise function are determined as the maximum horizontal principal stress vector distribution function of the target region.

[0154] It is understandable that the method for determining the maximum horizontal principal stress vector distribution function is the same as that for the minimum horizontal principal stress vector distribution function, and will not be elaborated here.

[0155] According to the regional horizontal stress determination method based on point-surface in-situ data fusion in this application embodiment, N measuring points and their respective coordinate values ​​are determined within the target area of ​​the coal mine working face according to preset rules; the maximum and minimum horizontal principal stress values ​​of each of the N measuring points are obtained; for each measuring point, the measuring point is determined as the center measuring point, and at least one adjacent measuring point of the center measuring point is determined; based on the coordinate value of the center measuring point, the maximum and minimum horizontal principal stress values ​​of the center measuring point, and the coordinate value, maximum and minimum horizontal principal stress values ​​of each adjacent measuring point, the method determines the regional horizontal stress. The maximum and minimum horizontal principal stress migration directions at the central measuring point are determined based on the principal stress values. The rate of change of the maximum and minimum horizontal principal stress between any two adjacent measuring points is also determined based on the maximum and minimum horizontal principal stress migration directions at each measuring point, as well as the rate of change of the maximum and minimum horizontal principal stress between any two adjacent measuring points. Furthermore, the maximum and minimum horizontal principal stress vector distribution data for the target area are determined based on the migration directions of the maximum and minimum horizontal principal stresses at each measuring point, and the rate of change of the maximum and minimum horizontal principal stresses between any two adjacent measuring points. By fusing seismic CT detection data and in-situ measured data from local point stress tests in the target area, the understanding of horizontal stress over a large area of ​​the mining area or working face is significantly improved. This enables quantitative assessment of regional horizontal stress, providing reliable data support for mine production design, thereby significantly improving mine design and production efficiency, greatly reducing the risk of disasters, and reducing the workload of disaster prevention and mitigation. Simultaneously, it provides a methodological basis for continuous online monitoring of horizontal stress in the working face area.

[0156] Figure 2 This is a flowchart of a regional horizontal stress determination device based on point-surface in-situ data fusion in an embodiment of this application.

[0157] like Figure 2 As shown, the regional horizontal stress determination device based on point-surface in-situ data fusion includes:

[0158] The first determining module 201 is used to determine N measuring points and their respective coordinate values ​​within the target area of ​​the coal mine working face according to preset rules; where N is an integer greater than or equal to 9.

[0159] The acquisition module 202 is used to acquire the maximum and minimum horizontal principal stress values ​​of each of the N measuring points;

[0160] The second determining module 203 is used to determine each measuring point as the central measuring point and to determine at least one adjacent measuring point of the central measuring point.

[0161] The third determining module 204 is used to determine the migration direction of the maximum horizontal principal stress of the center measuring point and the migration direction of the minimum horizontal principal stress of the center measuring point based on the coordinate value of the center measuring point, the maximum horizontal principal stress value of the center measuring point and the minimum horizontal principal stress value of the center measuring point, as well as the coordinate value of each adjacent measuring point, the maximum horizontal principal stress value of each adjacent measuring point and the minimum horizontal principal stress value of each adjacent measuring point.

[0162] The fourth determining module 205 is used to determine the rate of change of the maximum and minimum horizontal principal stress between any two adjacent measuring points based on the maximum and minimum horizontal principal stress values ​​of each of the N measuring points.

[0163] The fifth determining module 206 is used to determine the maximum horizontal principal stress vector distribution data and the minimum horizontal principal stress vector distribution data of the target area based on the maximum horizontal principal stress migration direction and the minimum horizontal principal stress migration direction of each measuring point, as well as the maximum horizontal principal stress change rate and the minimum horizontal principal stress change rate between any two adjacent measuring points.

[0164] According to one embodiment of this application, the fifth determining module includes:

[0165] The first connection submodule is used to connect the head and tail of the maximum horizontal principal stress migration direction of N measurement points to obtain the maximum horizontal principal stress migration direction of the target area.

[0166] The first determination submodule is used to determine the maximum horizontal principal stress vector distribution data of the target area based on the migration direction of the maximum horizontal principal stress in the target area and the rate of change of the maximum horizontal principal stress between any two adjacent measuring points.

[0167] The second connection submodule is used to connect the head and tail of the minimum horizontal principal stress migration direction of each of the N measuring points to obtain the minimum horizontal principal stress migration direction of the target area.

[0168] The second determination submodule is used to determine the minimum horizontal principal stress vector distribution data of the target area based on the minimum horizontal principal stress migration direction of the target area and the minimum horizontal principal stress change rate between any two adjacent measuring points.

[0169] According to the embodiment of this application, the regional horizontal stress determination device based on point-surface in-situ data fusion determines N measuring points and their respective coordinate values ​​within a target area of ​​a coal mine working face according to preset rules; obtains the maximum and minimum horizontal principal stress values ​​for each of the N measuring points; for each measuring point, determines the measuring point as the center measuring point, and determines at least one adjacent measuring point of the center measuring point; based on the coordinate values ​​of the center measuring point, the maximum and minimum horizontal principal stress values ​​of the center measuring point, and the coordinate values, maximum and minimum horizontal principal stress values ​​of each adjacent measuring point, the device determines the maximum and minimum horizontal principal stress values ​​of each adjacent measuring point. The maximum and minimum horizontal principal stress migration directions at the central measuring point are determined based on the principal stress values. The rate of change of the maximum and minimum horizontal principal stress between any two adjacent measuring points is also determined based on the maximum and minimum horizontal principal stress migration directions at each measuring point, as well as the rate of change of the maximum and minimum horizontal principal stress between any two adjacent measuring points. Furthermore, the maximum and minimum horizontal principal stress vector distribution data for the target area are determined based on the migration directions of the maximum and minimum horizontal principal stresses at each measuring point, and the rate of change of the maximum and minimum horizontal principal stresses between any two adjacent measuring points. By fusing seismic CT detection data and in-situ measured data from local point stress tests in the target area, the understanding of horizontal stress over a large area of ​​the mining area or working face is significantly improved. This enables quantitative assessment of regional horizontal stress, providing reliable data support for mine production design, thereby significantly improving mine design and production efficiency, greatly reducing the risk of disasters, and reducing the workload of disaster prevention and mitigation. Simultaneously, it provides a methodological basis for continuous online monitoring of horizontal stress in the working face area.

[0170] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of that feature. In the description of this invention, "a plurality of" means at least two, such as two, three, etc., unless otherwise explicitly specified.

[0171] In this invention, unless otherwise explicitly specified and limited, "above" or "below" the second feature can mean that the first feature is in direct contact with the second feature, or that the first feature is in indirect contact with the second feature through an intermediate medium. Furthermore, "above," "over," and "on top" of the second feature can mean that the first feature is directly above or diagonally above the second feature, or simply that the first feature is at a higher horizontal level than the second feature. "Below," "below," and "under" the second feature can mean that the first feature is directly below or diagonally below the second feature, or simply that the first feature is at a lower horizontal level than the second feature.

[0172] In this invention, the terms "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., refer to a specific feature, structure, material, or characteristic described in connection with that embodiment or example, which is included in at least one embodiment or example of the invention. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.

[0173] Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present invention. Those skilled in the art can make changes, modifications, substitutions and variations to the above embodiments within the scope of the present invention.

Claims

1. A method for determining regional horizontal stress based on point-surface in-situ data fusion, characterized in that, The stress includes the maximum horizontal principal stress and the minimum horizontal principal stress, and the method includes: According to preset rules, N measuring points and their respective coordinate values ​​are determined within the target area of ​​the coal mine working face; wherein, N is an integer greater than or equal to 9; Obtain the maximum and minimum horizontal principal stress values ​​for each of the N measuring points; For each measuring point, the measuring point is determined as the central measuring point, and at least one adjacent measuring point of the central measuring point is determined; Based on the coordinates of the central measuring point, the maximum horizontal principal stress of the central measuring point, and the minimum horizontal principal stress of the central measuring point, as well as the coordinates of each of the adjacent measuring points, the maximum horizontal principal stress of each of the adjacent measuring points, and the minimum horizontal principal stress of each of the adjacent measuring points, the migration direction of the maximum horizontal principal stress of the central measuring point and the migration direction of the minimum horizontal principal stress of the central measuring point are determined respectively. Based on the maximum and minimum horizontal principal stress values ​​of each of the N measuring points, the maximum and minimum horizontal principal stress change rates between any two adjacent measuring points are determined respectively. Based on the migration direction of the maximum and minimum horizontal principal stresses at each measuring point, as well as the rate of change of the maximum and minimum horizontal principal stresses between any two adjacent measuring points, the maximum and minimum horizontal principal stress vector distribution data of the target area are determined respectively. Specifically, determining the maximum and minimum horizontal principal stress vector distribution data of the target region based on the migration direction of the maximum and minimum horizontal principal stresses at each measuring point, and the rate of change of the maximum and minimum horizontal principal stresses between any two adjacent measuring points, includes: Based on the minimum horizontal stress migration continuity function expression between each measuring point and its two adjacent measuring points in the minimum horizontal stress migration direction, the first function value and the second function value of each measuring point are determined respectively. In response to the fact that the first function value and the second function value are the same, and the continuous function expression of minimum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of minimum horizontal stress migration satisfies the function continuity condition, the continuous function expressions of minimum horizontal stress migration corresponding to the first function value and the second function value are merged to obtain the merged continuous function expression of minimum horizontal stress migration. In response to the difference between the first function value and the second function value, and / or the continuous function expression of minimum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of minimum horizontal stress migration does not satisfy the function continuity condition, the continuous function expression of minimum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of minimum horizontal stress migration is determined as the first piecewise function; The merged minimum horizontal stress migration continuous function expression and the first piecewise function are determined as the minimum horizontal principal stress vector distribution function of the target region.

2. The method according to claim 1, characterized in that, The determination of the maximum and minimum horizontal principal stress vector distribution data of the target region based on the migration direction of the maximum and minimum horizontal principal stresses at each measuring point, and the rate of change of the maximum and minimum horizontal principal stresses between any two adjacent measuring points, includes: Connect the head and tail of the maximum horizontal principal stress migration direction of each of the N measuring points to obtain the maximum horizontal principal stress migration direction of the target area. Based on the migration direction of the maximum horizontal principal stress in the target area and the rate of change of the maximum horizontal principal stress between any two adjacent measuring points, the maximum horizontal principal stress vector distribution data of the target area is determined. Connect the beginning and end of the minimum horizontal principal stress migration direction of each of the N measuring points to obtain the minimum horizontal principal stress migration direction of the target area. Based on the migration direction of the minimum horizontal principal stress in the target area and the rate of change of the minimum horizontal principal stress between any two adjacent measuring points, the minimum horizontal principal stress vector distribution data of the target area are determined.

3. The method according to claim 2, characterized in that, The determination of the migration direction of the maximum horizontal principal stress and the migration direction of the minimum horizontal principal stress at the central measuring point, based on the coordinates of the central measuring point, the maximum horizontal principal stress of the central measuring point, and the minimum horizontal principal stress of the central measuring point, as well as the coordinates of each adjacent measuring point, the maximum horizontal principal stress of each adjacent measuring point, and the minimum horizontal principal stress of each adjacent measuring point, respectively, includes: Based on the coordinates of the central measuring point, the maximum horizontal principal stress of the central measuring point, the coordinates of each of the adjacent measuring points, and the maximum horizontal principal stress of each of the adjacent measuring points, the directional derivative is calculated to obtain the migration direction of the maximum horizontal principal stress between the central measuring point and each of the adjacent measuring points. The maximum value among the maximum horizontal principal stress migration directions between the central measuring point and at least one adjacent measuring point around it is determined as the maximum horizontal principal stress migration direction of the central measuring point. Based on the coordinates of the central measuring point, the maximum horizontal principal stress of the central measuring point, the coordinates of each of the adjacent measuring points, and the maximum horizontal principal stress of each of the adjacent measuring points, the directional derivative is calculated to obtain the migration direction of the maximum horizontal principal stress between the central measuring point and each of the adjacent measuring points. The maximum value among the minimum horizontal principal stress migration directions between the central measuring point and the at least one adjacent measuring point is determined as the minimum horizontal principal stress migration direction of the central measuring point.

4. The method according to claim 3, characterized in that, After determining the migration direction of the maximum horizontal principal stress and the migration direction of the minimum horizontal principal stress of the central measuring point based on the coordinate values ​​of the central measuring point, the maximum horizontal principal stress value of the central measuring point, and the minimum horizontal principal stress value of the central measuring point, as well as the coordinate values ​​of each of the adjacent measuring points, the maximum horizontal principal stress value of each of the adjacent measuring points, and the minimum horizontal principal stress value of each of the adjacent measuring points, the process further includes: For each central measuring point, the maximum rate of change of the horizontal principal stress at the central measuring point in the direction of maximum horizontal principal stress migration and the minimum rate of change of the horizontal principal stress at the central measuring point in the direction of minimum horizontal principal stress migration are determined respectively. According to a preset interpolation factor k, based on the maximum horizontal principal stress change rate of the N measuring points in their respective maximum horizontal principal stress migration directions and the minimum horizontal principal stress change rate of the N measuring points in their respective minimum horizontal principal stress migration directions, interpolation processing is performed on the maximum and minimum horizontal principal stress values ​​of the N measuring points respectively to obtain the maximum and minimum horizontal principal stress values ​​of kN measuring points respectively; where k is an integer greater than 0; The maximum horizontal principal stress migration direction and the minimum horizontal principal stress migration direction of each of the N measuring points are evenly divided to obtain the maximum horizontal principal stress migration direction and the minimum horizontal principal stress migration direction of each of the kN measuring points after the even division.

5. The method according to claim 3, characterized in that, After determining the migration direction of the maximum horizontal principal stress and the migration direction of the minimum horizontal principal stress of the central measuring point based on the coordinate values ​​of the central measuring point, the maximum horizontal principal stress value of the central measuring point, and the minimum horizontal principal stress value of the central measuring point, as well as the coordinate values ​​of each of the adjacent measuring points, the maximum horizontal principal stress value of each of the adjacent measuring points, and the minimum horizontal principal stress value of each of the adjacent measuring points, the process further includes: The maximum horizontal principal stress variation rate at the central measuring point along the x-direction and the maximum horizontal principal stress variation rate at the central measuring point along the y-direction are determined; wherein, the maximum horizontal principal stress variation rate at the central measuring point along the x-direction and the maximum horizontal principal stress variation rate at the central measuring point along the y-direction are both obtained by decomposing the maximum horizontal principal stress variation rate between the central measuring point and the adjacent measuring point in the direction of migration of the maximum horizontal principal stress along the x and y directions, respectively. The rate of change of the maximum horizontal principal stress at the central measuring point along the x-direction is integrated to obtain a first integral result; the rate of change of the maximum horizontal principal stress at the central measuring point along the y-direction is integrated to obtain a second integral result. Based on the first integration result and the second integration result, the expression for the maximum horizontal stress migration continuity function between the central measuring point and its adjacent measuring points in the direction of maximum horizontal stress migration is determined. Determine the minimum rate of change of the horizontal principal stress at the central measuring point along the x-direction and the minimum rate of change of the horizontal principal stress at the central measuring point along the y-direction; The rate of change of the minimum horizontal principal stress at the central measuring point along the x-direction is integrated to obtain the third integral result; the rate of change of the minimum horizontal principal stress at the central measuring point along the y-direction is integrated to obtain the fourth integral result. Based on the third and fourth integral results, the expression for the minimum horizontal stress migration continuity function between the central measuring point and its adjacent measuring points in the direction of minimum horizontal stress migration is determined.

6. The method according to claim 5, characterized in that, Based on the migration direction of the maximum and minimum horizontal principal stresses at each measuring point, and the rate of change of the maximum and minimum horizontal principal stresses between any two adjacent measuring points, the maximum and minimum horizontal principal stress vector distribution data of the target region are determined, including: Based on the maximum horizontal stress migration continuity function expression between each measuring point and its two adjacent measuring points in the direction of maximum horizontal stress migration, the third and fourth function values ​​of each measuring point are determined respectively; In response to the fact that the third function value is the same as the fourth function value, and the maximum horizontal stress migration continuous function expression between the measuring point and its two adjacent measuring points in the maximum horizontal stress migration direction satisfies the function continuity condition, the maximum horizontal stress migration continuous function expressions corresponding to the third function value and the fourth function value are merged to obtain the merged maximum horizontal stress migration continuous function expression. In response to the fact that the third function value is different from the fourth function value, and / or the continuous function expression of the maximum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of maximum horizontal stress migration does not satisfy the function continuity condition, the continuous function expression of the maximum horizontal stress migration between the measuring point and its two adjacent measuring points in the direction of maximum horizontal stress migration is determined as the second piecewise function. The combined maximum horizontal stress migration continuous function expression and the second piecewise function are determined as the maximum horizontal principal stress vector distribution function of the target region.

7. The method according to claim 1, characterized in that, The determination of the rate of change of the maximum and minimum horizontal principal stresses between any two adjacent measuring points, based on the maximum and minimum horizontal principal stresses of each of the N measuring points, includes: The maximum horizontal principal stress value at the central measuring point is subtracted from the maximum horizontal principal stress value at each of the adjacent measuring points to obtain a first difference; based on the first difference, the rate of change of the maximum horizontal principal stress at the central measuring point and at each of the adjacent measuring points is determined. The minimum horizontal principal stress value of the central measuring point is subtracted from the minimum horizontal principal stress value of each of the adjacent measuring points to obtain a second difference; based on the second difference, the rate of change of the minimum horizontal principal stress of the central measuring point and each of the adjacent measuring points is determined.

8. The method according to claim 1, characterized in that, The determination of N measuring points within the target area of ​​the coal mine working face according to preset rules includes: The target area is evenly divided into multiple blocks, and the coordinate regions of each of the multiple blocks are determined. The coordinate regions of each of the multiple blocks are compared with a preset region. In response to a coordinate region falling into the preset coordinate region, the block corresponding to that coordinate region is removed to obtain the N blocks; wherein, the preset coordinate region is a wave velocity anomaly region. Based on the coordinate regions of each of the N blocks, the coordinates of the center point of each of the N blocks are determined respectively; Based on the coordinates of the center point of each of the N blocks, N measurement points within the target area are determined.

9. A device for determining regional horizontal stress based on point-surface in-situ data fusion, characterized in that, The apparatus, applicable to the method of any one of claims 1-8, comprises: The first determining module is used to determine N measuring points and their respective coordinate values ​​within a target area of ​​a coal mine working face according to preset rules; wherein, N is an integer greater than or equal to 9; The acquisition module is used to acquire the maximum and minimum horizontal principal stress values ​​of each of the N measuring points; The second determining module is used to determine each measuring point as a central measuring point and to determine at least one adjacent measuring point of the central measuring point. The third determining module is used to determine the migration direction of the maximum horizontal principal stress of the central measuring point and the migration direction of the minimum horizontal principal stress of the central measuring point based on the coordinate value of the central measuring point, the maximum horizontal principal stress value of the central measuring point and the minimum horizontal principal stress value of the central measuring point, as well as the coordinate value of each of the adjacent measuring points, the maximum horizontal principal stress value of each of the adjacent measuring points and the minimum horizontal principal stress value of each of the adjacent measuring points. The fourth determining module is used to determine the rate of change of the maximum and minimum horizontal principal stress between any two adjacent measuring points based on the maximum and minimum horizontal principal stress values ​​of each of the N measuring points. The fifth determining module is used to determine the maximum horizontal principal stress vector distribution data and the minimum horizontal principal stress vector distribution data of the target area based on the maximum horizontal principal stress migration direction and the minimum horizontal principal stress migration direction of each measuring point, as well as the maximum horizontal principal stress change rate and the minimum horizontal principal stress change rate between any two adjacent measuring points.

10. The apparatus according to claim 9, characterized in that, The fifth determining module includes: The first connection submodule is used to connect the head and tail of the maximum horizontal principal stress migration direction of each of the N measuring points to obtain the maximum horizontal principal stress migration direction of the target area. The first determining submodule is used to determine the maximum horizontal principal stress vector distribution data of the target area based on the migration direction of the maximum horizontal principal stress in the target area and the rate of change of the maximum horizontal principal stress between any two adjacent measuring points. The second connection submodule is used to connect the beginning and end of the minimum horizontal principal stress migration direction of each of the N measuring points to obtain the minimum horizontal principal stress migration direction of the target area. The second determining submodule is used to determine the minimum horizontal principal stress vector distribution data of the target area based on the minimum horizontal principal stress migration direction of the target area and the minimum horizontal principal stress change rate between any two adjacent measuring points.

Images