Insar-based measurement method for angle of critical deformation of coal mining area, and system

Abstract

The present invention relates to the technical field of prevention and control of geological disasters in mountainous coal mining areas. Disclosed are an InSAR-based measurement method for the angle of critical deformation of a coal mining area, and a system. The method comprises: acquiring a mining-induced surface deformation field image by means of small baseline subset-interferometric synthetic aperture radar; in the mining-induced surface deformation field image, separately drawing profile lines along a coal seam strike and a coal seam dip, and selecting target deformation points on the profile lines as points to be monitored; according to the incident angle, the flight direction and the azimuth angle of a satellite, obtaining vertical and horizontal deformation values; using a probability integral method to calculate parameters such as inclination, curvature and horizontal deformation; and drawing a curve graph on the basis of the vertical deformation of a target region and the above parameters, so as to obtain a corresponding angle of critical deformation. The InSAR-based measurement method for angle of critical deformation of a coal mining area provided by the present invention has the advantages of improved temporal resolution and reduced phase noises and the like, can achieve high precision and high efficiency, and can generate clear and visible Excel tables for surface deformation data of mining areas and corresponding profile data.

Description

A coal mine movement angle measurement method and system based on InSAR Technical Field

[0001] The present invention relates to the technical field of coal mine geological disaster prevention and control in mountainous areas, and in particular to a coal mine movement angle measurement method and system based on InSAR. Background Art

[0002] Ground subsidence in coal mining areas is one of the most common and typical types of geological hazards. It results from continuous surface deformation caused by underground mining, and carries the risk of secondary hazards such as landslides, collapses, and mudslides. The large-scale development of coal resources has made an indelible contribution to the development of my country's national economy and provided a strong foundation for the rapid economic and social development of the country. Data from the National Bureau of Statistics shows that in 2022, my country produced 4.5 billion tons of raw coal, a year-on-year increase of 9.0%, and imported 290 million tons of coal, a year-on-year decrease of 9.2%.

[0003] However, with the widespread exploitation of coal resources, geological disasters have become frequent in my country's coal mining areas. This damage can cause deformation and displacement of the rock mass, ultimately leading to surface subsidence above the mined-out areas. The key to resolving these issues lies in determining mining subsidence parameters. These parameters are related to factors such as the displacement basin, the amount of subsidence, the strike and dip of the coal seam, and the characteristics of the overlying rock strata. By taking two sections on the plane of the displacement basin after coal seam mining, which align with the strike and dip of the coal seam, the positional relationship between the displacement basin and the two sections can be visualized, and the displacement angle can be calculated based on this relationship.

[0004] Currently, commonly used methods for obtaining surface movement parameters include numerical simulation, similar materials, and rock formation movement observation stations. We have invented an approach using InSAR technology to obtain these parameters. Common methods for observing surface subsidence in the study area include leveling, Global Positioning System (GPS), and total station measurements. However, these methods struggle to accurately capture the complete characteristics of spatial surface subsidence and suffer from limitations such as a small detection range, high cost, and significant labor consumption. InSAR technology is widely used to obtain large-scale, high-precision deformation data. One study used time-series D-InSAR to monitor subsidence caused by underground coal mining activities in mountainous areas, obtaining parameters such as the impact angle, pull angle, and starting distance at different mining stages. Other researchers have developed a new method based on time-series synthetic aperture radar imagery to monitor mining-induced surface subsidence, thereby determining the average subsidence rate of the mine over many years. However, D-InSAR technology is susceptible to decoherence and atmospheric delay, making it incapable of providing continuous time-series deformation data.

[0005] Summary of the Invention

[0006] In view of the above-mentioned problems, the present invention is proposed.

[0007] Therefore, the technical problem solved by the present invention is that the existing traditional coal mining area movement angle measurement method has the problems of low accuracy and low efficiency, making it difficult to monitor and predict the surface movement and deformation of the coal mining area in real time and accurately, resulting in the inability to timely prevent and deal with geological disasters caused by surface movement.

[0008] To solve the above technical problems, the present invention provides the following technical solutions: a method for measuring the movement angle in a coal mine area based on InSAR, comprising:

[0009] Differential interferometry short baseline set time series analysis technology is used to obtain mining-induced surface deformation field images; in the mining-induced surface deformation field images, profile lines are drawn along the strike and dip of the coal seam, and n target deformation points are selected on the profile lines as monitoring points; based on the satellite's incident angle, satellite flight direction, and azimuth, the monitoring point data are two-dimensionally decomposed, ignoring the north-south deformation contribution, to obtain vertical and horizontal deformation values; combined with the two-dimensional decomposition vertical and horizontal deformation results, the inclination, curvature, and horizontal deformation parameters are calculated according to the probability integral method; a curve graph is drawn based on the vertical deformation of the target area and the above parameters, and the vertical deformation within 10 mm is used as the research range. Combined with the strike and dip profiles, the corresponding movement angle value is obtained.

[0010] As a preferred solution of the InSAR-based coal mining area movement angle measurement method described in the present invention, the acquisition of mining-induced surface deformation field images includes image correction, interference processing, coherent phase analysis and data interpretation; the image correction includes atmospheric correction and radiation correction, the interference processing includes calculating interference images at different time points, the coherent phase analysis includes obtaining the coherent phase history of each pixel and detecting surface displacement signals, and the data interpretation includes obtaining the uplift and subsidence phenomena of the target area based on the data; the selection of n target deformation points as monitoring points includes taking the length of the coal seam strike and dip profile lines passing through all goafs, and the profile line lengths are equal, so that the number of monitoring points taken is equal.

[0011] As a preferred solution of the InSAR-based coal mine area movement angle measurement method of the present invention, the two-dimensional decomposition is represented as follows:

[0012] Among them, θ represents the incident angle of the satellite image, represents the angle between the satellite flight direction and the north direction, α-3 / 2π represents the angle between the projection of the satellite line of sight on the ground and the north direction, T1 and T2 represent the acquisition time of the two satellite data in the same ascending orbit, represents the radar line of sight deformation from T1 to T2, d U Represents the deformation in the vertical direction, d Erepresents the deformation in the horizontal east-west direction, and the deformation in the line of sight of the descending radar is T1′ and T2′ represent the data acquisition time, θ′ represents the incident angle of the satellite descending image, and α and α′ represent the angles between the satellite flight direction and the north direction, respectively.

[0013] As a preferred solution of the InSAR-based coal mine movement angle measurement method of the present invention, the probability integral method is expressed as:

[0014] W0=mq cosα

[0015] Among them, W0 is the maximum surface subsidence value, r z is the main influencing radius, λ is the horizontal integrated distance, m is the mining thickness of the ore body, q is the settlement coefficient, α is the inclination of the ore layer; A is a constant, H is the mining depth, tanβ is the tangent value of the main influencing angle, i represents the inclination, k represents the curvature, U represents the horizontal displacement value, B z represents the horizontal movement coefficient, and ε represents the horizontal deformation.

[0016] As a preferred solution of the InSAR-based coal mining area movement angle measurement method described in the present invention, the movement angle is calculated under full mining conditions, and a set of critical deformation values ​​are used to determine the dangerous movement boundary using the current standard, i = ± 3 mm / m, ε = ± 2 mm / m, k = ± 0.2*10 -3 / m to generate an intersection, take the outermost point of all critical intersections within 10mm, and project it vertically to the ground point and the line connecting the boundary point of the goaf and the horizontal line, and the angle between them is the moving angle.

[0017] As a preferred solution of the InSAR-based coal mining area movement angle measurement method described in the present invention, the movement angle includes, on the basis of the defined movement angle point, if the vertical settlement curve exceeds 10 mm at the outermost side of the profile, selecting the outermost point of the intersection of the critical deformation curve, and the angle between the intersection line of the line passing through the point and the outermost intersection point of the target coal seam and the horizontal line is the movement angle.

[0018] As a preferred solution of the InSAR-based coal mine movement angle measurement method described in the present invention, the movement angle also includes when there is a discontinuity in the goaf through which the profile line passes, based on the time period of full mining, corresponding to the fluctuation trend of the coal seam and settlement curve, selecting the intersection of the critical deformation value and vertically projecting it on the ground point, and the angle between the intersection of the line passing through the point and the outermost intersection of the target coal seam and the horizontal line is the movement angle.

[0019] Another object of the present invention is to provide an InSAR-based coal mine area movement angle measurement system, which can solve the problem of accurately and efficiently monitoring and analyzing surface movement and deformation in coal mine areas under complex geological conditions by constructing a coal mine area movement angle measurement system, thereby improving the ability to prevent and deal with geological disasters.

[0020] To solve the above technical problems, the present invention provides the following technical solutions: a coal mine area movement angle measurement system based on InSAR, comprising: a data acquisition module, a deformation field analysis module, a deformation parameter calculation module and a movement angle value extraction module; the data acquisition module is used to acquire InSAR satellite data, perform atmospheric and radiation correction, calculate interference images at different time points, analyze the coherent phase history of each pixel, and detect surface displacement signals; the deformation field analysis module is used to analyze the processed data using differential interferometry short baseline set timing analysis technology, extract surface deformation information, draw profile lines on the coal seam strike and dip, and select monitoring points; the deformation parameter calculation module is used to calculate the inclination, curvature and horizontal deformation parameters based on the two-dimensional decomposition results combined with the probability integral method, taking into account the satellite incident angle, flight direction and azimuth; the movement angle value extraction module is used to draw a curve graph and determine the movement angle value based on the vertical deformation and other calculated parameters, and present the data in the form of a chart.

[0021] A computer device includes a memory and a processor, wherein the memory stores a computer program, and is characterized in that when the processor executes the computer program, the steps of the above-mentioned coal mine area movement angle measurement method based on InSAR are implemented.

[0022] A computer-readable storage medium having a computer program stored thereon, characterized in that when the computer program is executed by a processor, the steps of the above-mentioned coal mine area movement angle measurement method based on InSAR are implemented.

[0023] Beneficial effects of the present invention: Compared with the traditional D-InSAR technology, the coal mine area movement angle measurement method provided by the present invention adopts the SBAS-InSAR technology that inherits the advantages of the traditional InSAR method, and at the same time has the advantages of improved time resolution and reduced phase noise, etc., and solves the problems of D-InSAR technology that are prone to decoherence, susceptibility to atmospheric delay, and inability to provide continuous time series deformation data.

[0024] Compared with traditional numerical simulation, similar materials, rock movement method observation station and other technologies, the SBAS-InSAR technology adopted in this method has the advantages of high precision and high efficiency, which solves a series of problems such as the limitations of traditional methods such as small detection range, high cost and high manpower consumption.

[0025] Compared to using SBAS-InSAR alone, this method combines SBAS-InSAR with GIS. SBAS-InSAR extracts surface deformation, while GIS classifies and displays the surface deformation values ​​extracted by SBAS-InSAR, extracting profile data. This results in an Excel spreadsheet showing clearly visible surface deformation data for the mining area and corresponding profile data. This solves a range of issues, including data display and effective aggregation.

[0026] Compared with the single 1 probability integration method, this method effectively integrates SBAS-InSAR technology and the probability integration method. The detection accuracy of SBAS-InSAR technology can reach the millimeter level. Combined with the advantages of fast convergence and high accuracy of the probability integration method, the obtained motion parameter results are more reliable. BRIEF DESCRIPTION OF THE DRAWINGS

[0027] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0028] FIG1 is an overall flow chart of a method for measuring movement angle in a coal mine area based on InSAR, provided by one embodiment of the present invention.

[0029] FIG2 is a comparison diagram of load balancing rates of multiple algorithms under the same task request in a method for measuring movement angle in a coal mine area based on InSAR provided by a second embodiment of the present invention.

[0030] FIG3 is a diagram showing the running time of a coal mine area movement angle measurement method based on InSAR at different repetition times provided by the second embodiment of the present invention.

[0031] FIG4 is a schematic diagram of strike curvature and settlement curves of a coal mine movement angle measurement method based on InSAR provided by the fourth embodiment of the present invention.

[0032] FIG5 is a schematic diagram of strike movement angle calculation of a coal mine movement angle measurement method based on InSAR provided by the fourth embodiment of the present invention.

[0033] FIG6 is a schematic diagram of a dip curvature and settlement curve of a coal mine movement angle measurement method based on InSAR provided by the fourth embodiment of the present invention.

[0034] FIG7 is a schematic diagram of a coal mine area movement angle measurement method based on InSAR provided by the fourth embodiment of the present invention. DETAILED DESCRIPTION

[0035] To make the above-mentioned objects, features, and advantages of the present invention more clearly understood, the following detailed description of the specific embodiments of the present invention is given in conjunction with the accompanying drawings. It is obvious that the described embodiments are only part of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by ordinary persons in this field without creative work should fall within the scope of protection of the present invention.

[0036] In the following description, many specific details are set forth to facilitate a full understanding of the present invention. However, the present invention may also be implemented in other ways different from those described herein. Those skilled in the art may make similar generalizations without violating the connotation of the present invention. Therefore, the present invention is not limited to the specific embodiments disclosed below.

[0037] Example 1

[0038] 1 , which is an embodiment of the present invention, provides a method for measuring movement angle in a coal mine area based on InSAR, including:

[0039] Differential interferometry short baseline set time series analysis technology is used to obtain mining-induced surface deformation field images.

[0040] In the mining-induced surface deformation field image, section lines are drawn along the strike and dip of the coal seam, and n target deformation points are selected on the section lines as monitoring points.

[0041] According to the satellite's incident angle, flight direction, and azimuth, the monitoring point data is subjected to a two-dimensional decomposition ignoring the north-south deformation contribution to obtain the vertical and horizontal deformation values.

[0042] Combining the two-dimensional decomposition of vertical and horizontal deformation results, the tilt, curvature and horizontal deformation parameters are calculated using the probability integration method.

[0043] A curve graph is drawn based on the vertical deformation of the target area and the above parameters. The vertical deformation within 10 mm is taken as the research range. Combined with the strike and dip profiles, the corresponding movement angle value is obtained.

[0044] Acquiring mining-induced surface deformation field images includes image correction, interferometric processing, coherent phase analysis, and data interpretation.

[0045] Image correction includes atmospheric correction and radiation correction, interference processing includes calculating interference images at different time points, coherent phase analysis includes obtaining the coherent phase history of each pixel and detecting surface displacement signals, and data interpretation includes obtaining the uplift and subsidence phenomena of the target area based on the data.

[0046] Select n target deformation points as monitoring points, including the coal seam strike and dip profile lines, and take the length of all goafs through which the profile lines are equal, so that the number of monitoring points taken is equal.

[0047] The two-dimensional decomposition is expressed as,

[0048] Among them, θ represents the incident angle of the satellite image, represents the angle between the satellite flight direction and the north direction, and α-3 / 2π represents the angle between the projection of the satellite line of sight on the ground and the north direction. T1 and T2 represent the acquisition time of the two satellite data in the same ascending orbit. represents the radar line of sight deformation from T1 to T2, d U Represents the deformation in the vertical direction, d E Represents the deformation in the horizontal east-west direction. The deformation in the line of sight of the descending radar is T1′ and T2′ represent the data acquisition time, θ′ represents the incident angle of the satellite descending image, and α and α′ represent the angles between the satellite flight direction and the north direction, respectively.

[0049] The probability integral method is expressed as,

[0050] W0=mq cosα

[0051] Among them, W0 is the maximum surface subsidence value, r z is the main influencing radius, λ is the horizontal integrated distance, m is the mining thickness of the ore body, q is the settlement coefficient, α is the inclination of the ore layer; A is a constant, H is the mining depth, and tanβ is the tangent value of the main influencing angle. i represents the inclination, k represents the curvature, U represents the horizontal displacement value, and B z represents the horizontal movement coefficient, and ε represents the horizontal deformation.

[0052] The movement angle is calculated under full mining conditions, and the critical deformation value and the current standard use a set of critical deformation values ​​to determine the dangerous movement boundary, i = ± 3mm / m, ε = ± 2mm / m, k = ± 0.2*10 -3 / m to generate an intersection, take the outermost point of all critical intersections within 10mm, and project it vertically to the ground point and the line connecting the boundary point of the goaf and the horizontal line, and the angle between them is the moving angle.

[0053] The movement angle includes taking points based on the defined movement angle. If the vertical settlement curve exceeds 10 mm at the outermost side of the section, the outermost point of the intersection of the critical deformation curve is selected, and the angle between the line passing through this point and the outermost intersection point of the target coal seam and the intersection line of the horizontal line is the movement angle.

[0054] The movement angle also includes when there is a discontinuity in the goaf passed by the profile line. According to the time period of full mining, the intersection of the critical deformation value is selected and vertically projected onto the ground point corresponding to the fluctuation trend of the coal seam and settlement curve. The angle between the intersection of the line passing through this point and the outermost intersection of the target coal seam and the horizontal line is the movement angle.

[0055] Example 2

[0056] 3 , which is an embodiment of the present invention, provides a coal mine movement angle measurement system based on InSAR, including:

[0057] Data acquisition module, deformation field analysis module, deformation parameter calculation module and movement angle value extraction module.

[0058] The data acquisition module is used to acquire InSAR satellite data, perform atmospheric and radiation correction, calculate interferometric images at different time points, analyze the coherent phase history of each pixel, and detect surface displacement signals.

[0059] The deformation field analysis module is used to analyze the processed data using differential interferometry short baseline set time series analysis technology, extract surface deformation information, draw profile lines on the coal seam strike and dip, and select monitoring points.

[0060] The deformation parameter calculation module is used to calculate the tilt, curvature and horizontal deformation parameters based on the two-dimensional decomposition results combined with the probability integration method, taking into account the satellite's incident angle, flight direction and azimuth.

[0061] The movement angle value extraction module is used to draw a curve graph and determine the movement angle value based on the vertical deformation and other calculated parameters, and present the data in the form of a chart.

[0062] Example 3

[0063] An embodiment of the present invention is different from the previous two embodiments in that:

[0064] If the functions are implemented in the form of 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 the present invention, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0065] The logic and / or steps represented in the flowcharts or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing the logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (e.g., a computer-based system, a system including a processor, or other system that can fetch and execute instructions from an instruction execution system, apparatus, or device). For purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transport a program for use by, or in conjunction with, an instruction execution system, apparatus, or device.

[0066] More specific examples (a non-exhaustive list) of computer-readable media include the following: an electrical connection with one or more wires (electronic devices), a portable computer disk cartridge (magnetic devices), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), a fiber optic device, and a portable compact disc read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium, followed by editing, deciphering, or processing in another suitable manner as necessary, and then stored in a computer memory.

[0067] It should be understood that various parts of the present invention can be implemented using hardware, software, firmware, or a combination thereof. In the above-described embodiments, multiple steps or methods can be implemented using software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented using hardware, as in another embodiment, any one of the following technologies known in the art or a combination thereof can be used: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application-specific integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0068] Example 4

[0069] 4 to 7 , an embodiment of the present invention provides a method for measuring movement angles in coal mines based on InSAR. To verify the beneficial effects of the present invention, scientific demonstration is conducted through economic benefit calculations and simulation experiments.

[0070] The present invention overcomes the shortcoming that mining subsidence parameters can only be obtained through ground monitoring stations or actual measurements. The method based on the fusion of SBAS-InSAR technology and the probability integration method can obtain mining subsidence parameters more quickly and accurately.

[0071] As shown in Figures 4 to 7, the present invention proposes a single coal seam movement angle measurement method based on SBAS-InSAR, which includes 128 monitoring points deployed along the strike and dip profiles of the coal seam, two-dimensional decomposition of the data at each monitoring point, calculation of three critical deformation curves, and finally, combining two corresponding geological profiles drawn by Rhino6, combined with the settlement curve, three critical deformation curves and geological profile map, the movement angle of the corresponding coal mining area in the mountainous area can be obtained according to the mining subsidence data: under the conditions of single coal seam mining, the strike movement angles are 48.96° and 54.07°; the dip movement angles are 67.66° and 61.56° respectively.

[0072] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit the present invention. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical solutions of the present invention, which should all be included in the scope of the claims of the present invention.

Claims

1. A method for measuring the movement angle of a coal mine area based on InSAR, characterized in that: include: The image of the mining-induced surface deformation field is obtained by using the short baseline set time series analysis technology of differential interferometry; In the mining-induced surface deformation field image, section lines are drawn along the strike and dip of the coal seam, and n target deformation points are selected as monitoring points on the section lines; According to the satellite's incident angle, satellite's flight direction, and azimuth, the monitoring point data is subjected to a two-dimensional decomposition ignoring the north-south deformation contribution to obtain the vertical and horizontal deformation values; Combining the vertical and horizontal deformation results of two-dimensional decomposition, the tilt, curvature and horizontal deformation parameters are calculated according to the probability integral method; A curve graph is drawn based on the vertical deformation of the target area and the above parameters. The vertical deformation within 10 mm is taken as the research range. The corresponding movement angle value is obtained by combining the strike and dip profiles.

2. The coal mine area movement angle measurement method based on InSAR according to claim 1, characterized in that: The acquisition of mining-induced surface deformation field images includes image correction, interference processing, coherent phase analysis and data interpretation; The corrected image includes atmospheric correction and radiation correction, the interferometric processing includes calculating interferometric images at different time points, the coherent phase analysis includes obtaining the coherent phase history of each pixel and detecting the surface displacement signal, and the data interpretation includes obtaining the uplift and subsidence phenomena of the target area according to the data; The said selecting n target deformation points as monitoring points includes the length of the coal seam strike and dip profile lines passing through all goafs, and the profile lines are equal in length, so that the number of monitoring points taken is equal.

3. The coal mine area movement angle measurement method based on InSAR as claimed in claim 2, characterized in that: The two-dimensional decomposition is expressed as, Among them, θ represents the incident angle of the satellite image, represents the angle between the satellite flight direction and the north direction, α-3 / 2π represents the angle between the projection of the satellite line of sight on the ground and the north direction, T1 and T2 represent the acquisition time of the satellite data of the same ascending orbit, represents the radar line of sight deformation from T1 to T2, d U represents the deformation in the vertical direction, d E represents the deformation in the horizontal east-west direction, and the deformation in the line of sight of the descending radar is T1′ and T2′ represent the data acquisition time, θ′ represents the incident angle of the satellite descending image, and α and α′ represent the angles between the satellite flight direction and the north direction, respectively.

4. The coal mine area movement angle measurement method based on InSAR as claimed in claim 3, characterized in that: The probability integral method is expressed as, W0=mq cosα Among them, W0 is the maximum surface subsidence value, rz is the main influence radius, λ is the horizontal integrated distance, m is the mining thickness of the ore body, q is the settlement coefficient, α is the inclination of the ore layer; A is a constant, H is the mining depth, tanβ is the tangent value of the main influence angle, i represents the inclination, k represents the curvature, U represents the horizontal displacement value, B z represents the horizontal movement coefficient, and ε represents the horizontal deformation.

5. The coal mine area movement angle measurement method based on InSAR as claimed in claim 4, characterized in that: The movement angle is calculated under full mining conditions. The critical deformation value and the current standard use a set of critical deformation values ​​to determine the dangerous movement boundary, i = ± 3mm / m, ε = ± 2mm / m, k = ± 0.2*10 -3 / m generates an intersection, and the outermost point of all critical intersections within 10mm is taken. The angle between the vertical projection to the ground point and the line connecting the boundary point of the goaf and the horizontal line is the moving angle.

6. The coal mine area movement angle measurement method based on InSAR as claimed in claim 5, characterized in that: The movement angle includes, based on the defined movement angle point, if the vertical settlement curve exceeds 10 mm at the outermost side of the profile, selecting the outermost point of the intersection of the critical deformation curve, and the angle between the line passing through this point and the outermost intersection point of the target coal seam and the intersection line of the horizontal line is the movement angle.

7. The coal mine area movement angle measurement method based on InSAR as claimed in claim 6, characterized in that: The movement angle also includes when there is a discontinuity in the goaf through which the profile line passes, based on the time period of full mining, corresponding to the fluctuation trend of the coal seam and settlement curve, selecting the intersection of the critical deformation value and vertically projecting it on the ground point, and the angle between the line passing through the point and the outermost intersection of the target coal seam and the horizontal line is the movement angle.

8. A system using the coal mine area movement angle measurement method based on InSAR as claimed in any one of claims 1 to 7, characterized in that: include: Data acquisition module, deformation field analysis module, deformation parameter calculation module and movement angle value extraction module; The data acquisition module is used to acquire InSAR satellite data, perform atmospheric and radiation correction, calculate interferometric images at different time points, analyze the coherent phase history of each pixel, and detect surface displacement signals; The deformation field analysis module is used to analyze the processed data using differential interferometry short baseline set time series analysis technology, extract surface deformation information, draw profile lines on the coal seam strike and dip, and select monitoring points; The deformation parameter calculation module is used to calculate the tilt, curvature and horizontal deformation parameters based on the two-dimensional decomposition result combined with the probability integration method, taking into account the satellite incident angle, flight direction and azimuth; The movement angle value extraction module is used to draw a curve graph and determine the movement angle value according to the vertical deformation and other calculated parameters, and present the data in the form of a chart.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the coal mine movement angle measurement method based on InSAR described in any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the coal mine movement angle measurement method based on InSAR described in any one of claims 1 to 7 are implemented.

Images