A method and system for characterizing stress field transfer direction of coal seam region

Abstract

This invention discloses a method and system for characterizing the stress field transfer direction in coal seams, relating to the field of coal seam safety monitoring technology. The method includes: establishing a quantitative coupling relationship between stress and wave velocity in coal and rock samples through laboratory testing; inverting the three-dimensional wave velocity field of the working face based on microseismic monitoring data, and obtaining the corresponding three-dimensional stress field distribution according to the coupling relationship; calculating the transfer vector of the anomaly area by tracking the temporal changes of the centroid of the stress anomaly area, and obtaining the stress gradient change vector by calculating the temporal changes of the stress field spatial gradient; performing consistency analysis on the transfer vector and the stress gradient change vector, and determining the corresponding direction as the dominant transfer direction of the stress field in the coal seam region when the consistency meets preset conditions. This invention, through dual-path dynamic analysis and cross-validation, can characterize the transfer direction of the coal seam stress field at a regional scale, providing a new technical approach for the monitoring and early warning of coal and rock dynamic disasters.

Description

Technical Field

[0001] This invention relates to the field of coal seam safety monitoring technology, and in particular to a method and system for characterizing the direction of stress field transfer in a coal seam region. Background Technology

[0002] During underground coal mining, factors such as mining disturbances, overburden structure evolution, and mechanical vibrations cause continuous redistribution and dynamic transfer of the stress field in the coal seam and its surrounding rock. The direction of stress field transfer directly affects the location and development trend of dynamic disasters such as coal and gas outbursts, rock bursts, and roof instability. Therefore, effectively characterizing the direction of stress field transfer in coal seams is one of the key technical issues in coal mine safety production and disaster early warning.

[0003] Existing coal seam stress monitoring technologies mainly rely on point measurement methods such as borehole stress gauges, strain gauges, or anchor bolt stress monitoring. These methods have high deployment and maintenance costs, limited monitoring range, difficulty in obtaining the spatial distribution characteristics of regional stress fields, and insufficient ability to characterize dynamic evolution processes, making it difficult to meet the analytical needs of regional stress transfer behavior under complex mining conditions.

[0004] In recent years, the development of microseismic monitoring and tomographic imaging technologies has made it possible to invert the three-dimensional wave velocity and stress fields of coal seams and surrounding rocks at a regional scale. However, existing technologies mainly focus on the static inversion and anomaly identification of wave velocity or stress fields. Systematic and reliable characterization methods are still lacking for the temporal transfer characteristics of stress anomaly zones and the dominant transfer direction of the regional stress field they reflect. In particular, how to stably determine the transfer direction of the regional stress field under complex and noisy monitoring conditions remains a problem that urgently needs to be solved by existing technologies. Summary of the Invention

[0005] To address the aforementioned technical problems in existing technologies, embodiments of the present invention provide a method and system for characterizing the stress field transfer direction in coal seams. This invention is not limited to the inversion of wave velocity fields or stress fields themselves, but rather improves the reliability and stability of stress field transfer direction characterization by determining the dominant transfer direction of the regional stress field through dual-path dynamic analysis of the stress field evolution process and based on consistency judgment. The technical solution is as follows: On the one hand, a method for characterizing the stress field transfer direction in a coal seam region is provided. The method includes: acquiring microseismic monitoring data of the working face of the coal seam to be monitored, and obtaining a three-dimensional wave velocity field of the working face based on the microseismic monitoring data; converting the three-dimensional wave velocity field into a corresponding three-dimensional stress field based on a pre-established correspondence between coal and rock stress and wave velocity; analyzing the three-dimensional stress field at different times during mining to determine stress anomaly regions, and obtaining a transfer vector of the stress anomaly regions based on the spatial position changes of the stress anomaly regions at adjacent times; calculating the spatial gradient based on the three-dimensional stress field, and obtaining a stress gradient change vector based on the changes of the spatial gradient at adjacent times; performing a direction consistency analysis on the transfer vector and the stress gradient change vector, and determining the dominant transfer direction of the stress field in the coal seam region based on the consistency analysis results.

[0006] Optionally, the relationship between coal and rock stress and wave velocity is obtained by performing wave velocity tests on coal samples and surrounding rock samples under different stress conditions and fitting the results using a quadratic function model.

[0007] Optionally, the three-dimensional stress field is analyzed at different times during the mining process to determine stress anomaly areas. Specifically, this includes: comparing the stress value of each unit in the three-dimensional stress field with the dynamic average stress value of the region; marking a stress anomaly point when the stress value of a certain unit exceeds a preset proportion of the dynamic average stress value; and combining spatially adjacent stress anomaly points through three-dimensional connectivity analysis to form a continuous stress anomaly area.

[0008] Optionally, the transfer vector of the stress anomaly region is determined based on the position change of the centroid coordinates of the stress anomaly region at adjacent moments during the mining process; the centroid coordinates are obtained by calculating the average value of the coordinates of all discrete points within the corresponding stress anomaly region.

[0009] Optionally, the stress gradient change vector is obtained by calculating the spatial gradient field of the three-dimensional stress field at adjacent time points. x , y , z The changes in the three directions are combined to obtain the result.

[0010] Optionally, the directional consistency analysis is based on the angle relationship between the transfer vector and the stress gradient change vector, including: calculating the cosine value of the angle between the transfer vector and the stress gradient change vector; comparing the cosine value with a preset consistency threshold; when the cosine value is greater than or equal to the preset consistency threshold, determining that the transfer vector and the stress gradient change vector are in the same direction, and taking the direction of the transfer vector and the stress gradient change vector as the dominant transfer direction of the stress field in the coal seam region.

[0011] On the other hand, a system for characterizing the stress field transfer direction in a coal seam region is also provided, used to implement the method for characterizing the stress field transfer direction in a coal seam region as described in any one of claims 1-6; the system includes: a wave velocity field acquisition module, used to acquire microseismic monitoring data of the working face of the coal seam to be monitored, and to obtain the three-dimensional wave velocity field of the working face based on the microseismic monitoring data; a stress field inversion module, used to convert the three-dimensional wave velocity field into a corresponding three-dimensional stress field based on a pre-established correspondence between coal and rock stress and wave velocity; an anomaly area analysis module, used to analyze the three-dimensional stress field at different times during mining, determine the stress anomaly area, and obtain the transfer vector of the stress anomaly area according to the spatial position change of the stress anomaly area at adjacent times; a stress gradient analysis module, used to calculate the spatial gradient based on the three-dimensional stress field, and obtain the stress gradient change vector according to the change of the spatial gradient at adjacent times; and a stress field transfer direction determination module, used to perform directional consistency analysis on the transfer vector and the stress gradient change vector, and determine the dominant transfer direction of the stress field in the coal seam region according to the consistency analysis results.

[0012] On the other hand, an electronic device is also provided, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method provided in the embodiments of the present invention.

[0013] On the other hand, a computer-readable storage medium is also provided, wherein program code is stored in the computer-readable storage medium, and the program code can be called by a processor to execute the method provided in the embodiments of the present invention.

[0014] This invention provides a method and system for characterizing the stress field transfer direction in a coal seam region. Compared with the prior art, this invention has the following advantages: (1) A high-resolution stress field at the regional scale is constructed by microseismic inversion and stress-wave velocity relationship to achieve non-contact and continuous stress field acquisition.

[0015] (2) By using the dual-path independent analysis and cross-validation of the stress anomaly zone transfer trajectory and stress field gradient change trend, the local anomalies and noise interference in the single-path analysis are effectively suppressed, and the robustness and reliability of direction identification are significantly improved.

[0016] (3) Through the consistency judgment mechanism, the dominant transfer direction of stress field in coal seam region can be stably and reliably identified.

[0017] (4) It can provide refined and dynamic technical support for coal mine production regulation, early warning of power disasters and evaluation of gas outburst mitigation effects. Attached Figure Description

[0018] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0019] Figure 1 This is a flowchart of a method for characterizing the stress field transfer direction in a coal seam region, provided in an embodiment of the present invention. Figure 2 This is a schematic diagram of a stress field transfer direction characterization system for coal seam regions provided in an embodiment of the present invention. Detailed Implementation

[0020] The technical solution of the present invention will now be described with reference to the accompanying drawings.

[0021] In embodiments of the present invention, words such as "exemplarily," "for example," etc., are used to indicate that something is an example, illustration, or description. Any embodiment or design described as "exemplary" in the present invention should not be construed as being more preferred or advantageous than other embodiments or designs. Specifically, the use of the word "exemplary" is intended to present the concept in a concrete manner. Furthermore, in embodiments of the present invention, the meaning expressed by "and / or" can be both, or either one.

[0022] To make the technical problems, technical solutions and advantages of the present invention clearer, a detailed description will be given below in conjunction with the accompanying drawings and specific embodiments.

[0023] Figure 1 This is a flowchart illustrating a method for characterizing the stress field transfer direction in a coal seam region according to an embodiment of the present invention. Figure 1 As shown, the method specifically includes the following steps: Step S102: Obtain microseismic monitoring data of the working face of the coal seam to be monitored, and obtain the three-dimensional wave velocity field of the working face based on the microseismic monitoring data.

[0024] Step S104: Based on the pre-established correspondence between coal and rock stress and wave velocity, the three-dimensional wave velocity field is converted into the corresponding three-dimensional stress field.

[0025] Step S106: Analyze the three-dimensional stress field at different times during the mining process, determine the stress anomaly region, and obtain the transfer vector of the stress anomaly region based on the spatial position change of the stress anomaly region at adjacent times.

[0026] Step S108: Calculate the spatial gradient based on the three-dimensional stress field, and obtain the stress gradient change vector according to the change of the spatial gradient at adjacent time points.

[0027] Step S110: Perform a direction consistency analysis on the transfer vector and the stress gradient change vector, and determine the dominant transfer direction of the stress field in the coal seam region based on the consistency analysis results.

[0028] Specifically, step S102 further includes: deploying a microseismic pickup at the working face of the coal seam to be monitored, and acquiring microseismic data of the working face through microseismic monitoring; and performing inversion on the microseismic data based on a tomographic imaging algorithm to obtain a three-dimensional wave velocity field covering the working face of the coal seam to be monitored.

[0029] Specifically, microseismic pickups are deployed in an array along the working face to ensure effective coverage of the coal seam area. Based on a preset time interval (e.g., 7 days), the time-period data collected by the microseismic pickups deployed on the coal seam working face are input into a tomographic imaging algorithm for processing. Wave velocity inversion is performed using the Simultaneous Iterative Reconstruction Technique (SIRT). The iteration terminates when the residual converges to less than 0.01 or reaches the maximum number of iterations (50), thus obtaining the three-dimensional wave velocity field distribution for the corresponding time period.

[0030] Specifically, in step S104, the correspondence between coal and rock stress and wave velocity is obtained by performing wave velocity tests on coal and surrounding rock samples under different stress conditions and fitting the results using a quadratic function model. Specifically, a stress-wave velocity function model is established based on the stress and wave velocity test data of the target test samples. The target test samples include roof rock samples, coal samples, and floor rock samples from the working face of the coal seam to be monitored.

[0031] Specifically, step S104 further includes the following steps: Stress-wave velocity tests were conducted on the target test samples (roof rock sample, coal sample, and floor rock sample) in the laboratory.

[0032] Specifically, roof rock samples, coal samples, and floor rock samples are collected sequentially along the vertical direction of the working face of the coal seam to be monitored. After the samples are processed into standard cylindrical specimens, ultrasonic wave velocity tests are conducted on the target test samples under different loading stress conditions in the laboratory to obtain the corresponding stress-wave velocity test data.

[0033] Specifically, the experimental testing method is as follows: The sample is processed into a cylinder with a diameter of 50 mm and a height of 100 mm; a uniaxial loading device is used to load the sample according to the incremental stress step of 0→5→10→15→20 MPa, and the sample is kept stable for 10 minutes at each step. At the same time, a longitudinal ultrasonic transducer is used to measure the wave velocity along the axial and radial directions to collect stress-wave velocity data; each stress step is repeated three times and the average value is taken to ensure the reliability of the data.

[0034] Based on the above test data, a regression fitting method is used to establish a quantitative correspondence between stress and wave velocity. Preferably, the stress-wave velocity function model is in quadratic function form, and its mathematical expression is:

[0035] σ The stress of the target test sample, V p The wave velocity of the target test sample. m , n and c These are the fitting parameters.

[0036] Based on the established stress-wave velocity function model, the three-dimensional wave velocity field obtained in step S102 is converted into the corresponding three-dimensional stress field distribution.

[0037] Specifically, step S106, which involves analyzing the three-dimensional stress field at different times during the mining process to determine areas of stress anomalies, includes the following steps: The stress value of each element in the three-dimensional stress field is compared with the dynamic average stress value of the region. When the stress value of a certain unit exceeds the preset proportion of the dynamic average stress value, it is marked as a stress anomaly point; By combining spatially adjacent stress anomaly points through three-dimensional connectivity analysis, a continuous stress anomaly region is formed.

[0038] Specifically, the preset ratio is set based on historical monitoring data or experience of the mining area, and is preferably 20%.

[0039] Specifically, each voxel (or mesh element) within the three-dimensional stress field is compared to the average stress value of its region. An anomaly is identified when the stress value exceeds 120% of the average stress. Through three-dimensional connectivity analysis, spatially adjacent anomalies are grouped to form three-dimensional stress anomaly regions, which are then marked. This preset ratio of 120% is based on statistical analysis of working face stress inversion results from 10 outburst coal seam mining operations in a certain mining area, effectively identifying anomaly regions. It can be slightly adjusted based on actual mining experience.

[0040] Specifically, the transfer vector of the stress anomaly region is determined based on the positional change of the centroid coordinates of the stress anomaly region at adjacent moments during the mining process; the centroid coordinates are obtained by calculating the average value of the coordinates of all discrete points within the corresponding stress anomaly region.

[0041] Specifically, based on the identified stress anomaly regions, their centroid coordinates are calculated, and an anomaly transfer vector is constructed to characterize the transfer path of the stress anomaly regions. V .

[0042] Specifically, the minimum inertia ellipsoid fitting method is used to construct a corresponding anomaly region fitting body for each stress anomaly region, and the coordinates of discrete points within the fitting body are extracted. x i , y i , z i Based on the coordinates of discrete points, calculate the centroid coordinates of the fitted volume. x , y , z ):

[0043] in, N This represents the total number of discrete points within the fitted body.

[0044] By comparing the changes in the centroid coordinates of the fitted volume at adjacent time points, an anomaly transfer vector is constructed. V .

[0045] Specifically, the stress gradient change vector in step S108 is obtained by calculating the spatial gradient field of the three-dimensional stress field at adjacent time points. x , y , z The changes in the three directions are combined to obtain the result.

[0046] Specifically, the directional consistency analysis in step S110 is based on the angle relationship between the transfer vector and the stress gradient change vector. Specifically, step S110 also includes the following steps: Step S1101: Calculate the cosine value R of the angle between the transfer vector and the stress gradient change vector. R is calculated as follows:

[0047] in, V Represents the anomaly transition vector. G This represents the stress gradient change vector.

[0048] Step S1102: Compare the cosine value of the included angle with a preset consistency threshold; when the cosine value of the included angle is greater than or equal to the preset consistency threshold, determine that the direction of the transfer vector is consistent with the direction of the stress gradient change vector, and take the direction of the transfer vector and the stress gradient change vector as the dominant transfer direction of the stress field in the coal seam region.

[0049] Specifically, the preset consistency threshold is calibrated based on historical data of the mining area, and is preferably 0.75.

[0050] Specifically, when RWhen the value is greater than 0.75, the direction is considered consistent, and the dominant transfer direction of the coal seam stress field is determined. This threshold of 0.75 is also determined based on data calibration of a longwall mining face in a certain outburst-prone area, and can be slightly adjusted according to actual mining experience.

[0051] To illustrate the feasibility of the coal seam stress field transfer direction characterization method in this embodiment of the invention, a simplified data evolution example is given below: The origin of the coordinate system is the starting point of the ventilation roadway at the working face. O (0, 0, 0), assuming different times during the mining process t 1 and t 2. Fitted centroid coordinates of a certain stress anomaly region V 1. V 2 are respectively: July 1st of a certain year t Moment 1: V 1=( x 1, y 1, z 1) = (12.0m, 8.5m, 3.2m), 10 days later on July 11th. t Time 2: V 2=( x 2, y 2, z 2) = (20.2m, 10.5m, 3.8m). The anomaly transfer vector Δ is calculated from this. V : Δ V=V 2- V 1 = (8.2m, 2.0m, 0.6m) Meanwhile, assuming that the spatial stress gradient change vector of the stress field calculated in this region is Δ G : Δ G= (2.0, 1.6m, 0.6m) Calculating the consistency coefficient between two vectors based on the vector cosine formula. R :

[0052] Consistency coefficient R If the value is greater than the preset threshold of 0.75, it indicates that the direction of the anomaly transfer vector is consistent with the direction of gradient change, and this direction is determined to be the dominant transfer direction of the stress field in the coal seam region.

[0053] The above examples clearly demonstrate the complete process of obtaining the transfer vector from the centroid coordinates of the anomaly zone, calculating the consistency coefficient by combining the gradient vector, and finally determining the dominant transfer direction of the stress field in the coal seam region, thus proving that the method of the present invention can be practically implemented and obtain quantitative results.

[0054] Figure 2 This is a schematic diagram of a stress field transfer direction characterization system for a coal seam region according to an embodiment of the present invention. Figure 2 As shown, the system includes: The wave velocity field acquisition module 10 is used to acquire microseismic monitoring data of the working face of the coal seam to be monitored, and to obtain the three-dimensional wave velocity field of the working face based on the microseismic monitoring data. The stress field inversion module 20 is used to convert the three-dimensional wave velocity field into the corresponding three-dimensional stress field based on the pre-established correspondence between coal and rock stress and wave velocity. The anomaly analysis module 30 is used to analyze the three-dimensional stress field at different times during the mining process, determine the stress anomaly area, and obtain the transfer vector of the stress anomaly area based on the spatial position change of the stress anomaly area at adjacent times. The stress gradient analysis module 40 is used to calculate the spatial gradient based on the three-dimensional stress field and obtain the stress gradient change vector according to the change of the spatial gradient at adjacent time points. The stress field transfer direction determination module 50 is used to perform directional consistency analysis on the transfer vector and the stress gradient change vector, and determine the dominant transfer direction of the stress field in the coal seam region based on the consistency analysis results.

[0055] The present invention also provides an electronic device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method provided in the embodiments of the present invention.

[0056] The present invention also provides a computer-readable storage medium storing program code, which can be called by a processor to execute the method provided in the embodiments of the present invention.

[0057] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

[0058] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working processes of the devices, apparatuses, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0059] In the several embodiments provided by this invention, it should be understood that the disclosed devices, apparatuses, and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative; for instance, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another device, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces; the indirect coupling or communication connection between devices or units may be electrical, mechanical, or other forms.

[0060] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0061] In addition, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0062] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0063] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for characterizing the direction of stress field transfer in a coal seam region, characterized in that, The method includes: Acquire microseismic monitoring data of the working face of the coal seam to be monitored, and obtain the three-dimensional wave velocity field of the working face based on the microseismic monitoring data; Based on the pre-established correspondence between coal and rock stress and wave velocity, the three-dimensional wave velocity field is converted into the corresponding three-dimensional stress field. The three-dimensional stress field is analyzed at different times during the mining process to determine the stress anomaly region, and the transfer vector of the stress anomaly region is obtained based on the spatial position change of the stress anomaly region at adjacent times. The spatial gradient is calculated based on the three-dimensional stress field, and the stress gradient change vector is obtained according to the change of the spatial gradient at adjacent time points. A directional consistency analysis is performed on the transfer vector and the stress gradient change vector. Based on the consistency analysis results, the dominant transfer direction of the stress field in the coal seam region is determined.

2. The method according to claim 1, characterized in that, The relationship between coal and rock stress and wave velocity was obtained by testing the wave velocity of coal and surrounding rock samples under different stress conditions and fitting the results using a quadratic function model.

3. The method according to claim 1, characterized in that, The three-dimensional stress field was analyzed at different times during the mining process to determine areas of stress anomalies, specifically including: The stress value of each element in the three-dimensional stress field is compared with the dynamic average stress value of the region. When the stress value of a certain unit exceeds a preset proportion of the dynamic average stress value, it is marked as a stress anomaly point; By combining spatially adjacent stress anomaly points through three-dimensional connectivity analysis, a continuous stress anomaly region is formed.

4. The method according to claim 1, characterized in that, The transfer vector of the stress anomaly region is determined based on the position change of the centroid coordinates of the stress anomaly region at adjacent moments during the mining process; the centroid coordinates are obtained by calculating the average value of the coordinates of all discrete points within the corresponding stress anomaly region.

5. The method according to claim 1, characterized in that, The stress gradient change vector is obtained by calculating the spatial gradient field of the three-dimensional stress field at adjacent time points. x , y , z The changes in the three directions are combined to obtain the result.

6. The method according to claim 1, characterized in that, The directional consistency analysis is based on the angle relationship between the transfer vector and the stress gradient change vector, and includes: Calculate the cosine of the angle between the transfer vector and the stress gradient change vector; The cosine value of the included angle is compared with a preset consistency threshold. When the cosine value of the included angle is greater than or equal to the preset consistency threshold, it is determined that the transfer vector is consistent with the direction of the stress gradient change vector, and the direction of the transfer vector and the stress gradient change vector is taken as the dominant transfer direction of the stress field in the coal seam region.

7. A system for characterizing the direction of stress field transfer in a coal seam region, characterized in that, A system for implementing the method for characterizing the stress field transfer direction in a coal seam region as described in any one of claims 1-6; the system comprises: The wave velocity field acquisition module is used to acquire microseismic monitoring data of the working face of the coal seam to be monitored, and to invert the three-dimensional wave velocity field of the working face based on the microseismic monitoring data. The stress field inversion module is used to convert the three-dimensional wave velocity field into a corresponding three-dimensional stress field based on the pre-established correspondence between coal and rock stress and wave velocity. The anomaly zone analysis module is used to analyze the three-dimensional stress field at different times during the mining process, determine the stress anomaly zone, and obtain the transfer vector of the stress anomaly zone based on the spatial position change of the stress anomaly zone at adjacent times. The stress gradient analysis module is used to calculate the spatial gradient based on the three-dimensional stress field, and to obtain the stress gradient change vector according to the change of the spatial gradient at adjacent time points. The stress field transfer direction determination module is used to perform directional consistency analysis on the transfer vector and the stress gradient change vector, and determine the dominant transfer direction of the stress field in the coal seam region based on the consistency analysis results.

8. An electronic device, characterized in that, include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements the method as claimed in any one of claims 1-6.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium contains program code that can be invoked by a processor to execute the method as described in any one of claims 1 to 6.

Images