Mining area earth surface accumulated risk automatic monitoring and early warning method based on InSAR technology

Abstract

The invention relates to the technical field of geological engineering monitoring, in particular to a mining area earth surface accumulated risk automatic monitoring and early warning method based on the InSAR technology. Importing a database and downloading a task to obtain original data; determining a synthetic aperture radar data range according to the obtained original data, generating an intensity graph of a super main image, and generating a reference InSAR deformation field; calculating a deformation value of a date corresponding to the intensity graph to generate a deformation field image, outputting a deformation field txt file, and monitoring an InSAR deformation field; according to the method, an early warning mechanism is set, the accumulated risk of the mining area earth surface is evaluated in real time, in this way, the method not only has the advantages in the aspects of large-range, high-precision and non-contact monitoring based on the InSAR technology, but also achieves efficient, precise and continuous monitoring of the deformation of the mining area earth surface by optimizing the data processing flow and the automatic algorithm, and the safety of the mining area earth surface is improved. And a solid technical defense line and an early warning barrier are built for safe mining of a mining area.

Description

Technical Field

[0001] The present invention relates to the field of geological engineering monitoring technology, and in particular to an automatic monitoring and early warning method for surface cumulative risk in a mining area based on InSAR technology. Background Art

[0002] In the process of mining mineral resources, the complexity and variability of geological structures and the continuous disturbance caused by mining activities inevitably lead to a series of geological disasters such as landslides and collapses. These disasters not only seriously erode the ecological environment of the mining area, but also directly threaten the safety of production operations and the quality of life of surrounding residents, posing a potential safety hazard that cannot be ignored. Traditional deformation monitoring methods, such as the global positioning system (GPS) and leveling technology, have shown certain advantages in their respective fields, but have also encountered significant bottlenecks in practical applications.

[0003] Although GPS technology, especially the setting of continuous observation stations, can capture surface deformation data in real time, its high initial investment cost and limited deployment density limit its ability to achieve comprehensive and detailed monitoring in a vast mining area, and it is difficult to meet the monitoring needs of high sensitivity to cumulative risks. Although leveling technology has excellent accuracy in vertical deformation monitoring, it is limited by its complex and time-consuming operation. It is difficult to respond quickly to the rapid changes in the mining environment, and the efficiency bottleneck is obvious. With the rapid development of remote sensing technology, synthetic aperture radar interferometry (InSAR) technology has gradually occupied a place in the field of mining deformation monitoring with its unique advantages of wide coverage, high accuracy and non-contact monitoring. DInSAR technology has initially realized the monitoring of mining deformation through differential interferometry processing, but its single observation limitation makes it difficult to fully capture the dynamic change process under continuous mining in the mining area. Although PSInSAR technology realizes long-term series monitoring by focusing on permanent scatterers, the interference pattern is prone to lose coherence in high deformation gradient areas, resulting in discontinuity or even loss of deformation information, which limits its application scope. SBASInSAR technology effectively improves the continuity of deformation field through the small baseline subset strategy, becoming a powerful tool for surface deformation monitoring in mining areas. However, for a long time, policy restrictions on data acquisition, high costs, and insufficient observation frequency have still made SBASInSAR technology face challenges such as complex processing strategies and low automation in the dynamic monitoring of surface cumulative risks in mining areas.

[0004] With the significant improvement in the performance of the Sentinel-1SAR satellite, its precise orbit determination and ultra-high registration rate have opened up a new path for automated time-series InSAR processing. At the same time, significant progress in SAR processing technology at the algorithm level, such as the development of advanced platforms such as GAMMA and GMTSAR, has further promoted the in-depth application of InSAR technology in the field of deformation monitoring.

[0005] In view of this, the present invention aims to overcome the limitations of the existing technology and proposes an automatic monitoring and early warning method for the cumulative risk of the surface of the mining area based on InSAR technology. This method is not only based on InSAR technology, which has advantages in large-scale, high-precision, and non-contact monitoring, but also achieves efficient, accurate, and continuous monitoring of the surface deformation of the mining area by optimizing the data processing process and automation algorithm, thus building a solid technical defense line and early warning barrier for the safe mining of the mining area. Summary of the invention

[0006] The purpose of the present invention is to provide an automated monitoring and early warning method for surface cumulative risks in mining areas based on InSAR technology. This method is not only based on InSAR technology and has advantages in large-scale, high-precision, and non-contact monitoring, but also optimizes data processing procedures and automation algorithms to achieve efficient, accurate, and continuous monitoring of surface deformation in mining areas, thereby building a solid technical defense line and early warning barrier for safe mining in mining areas.

[0007] To achieve the above purpose, the present invention adopts an automatic monitoring and early warning method for surface cumulative risk in mining areas based on InSAR technology, which comprises the following steps:

[0008] Import the database and download the task to obtain the original data;

[0009] Confirm the synthetic aperture radar data range based on the acquired raw data, generate the intensity map of the super master image, and generate the reference InSAR deformation field;

[0010] Calculate the deformation value of the intensity map corresponding to the date to generate the deformation field image, and output the deformation field txt file to monitor the InSAR deformation field;

[0011] Establish an early warning mechanism to conduct real-time assessment of the accumulated risks on the surface of the mining area.

[0012] Among them, in the step of importing the database and downloading the task to obtain the original data:

[0013] According to the scope of the mining area, download the digital elevation model of the mining area;

[0014] Preset time span and download synthetic aperture radar data according to the time span;

[0015] Import the database into the script and download the task, the downloaded task is the original data; the data in the database includes processing file path, temporary file, regular expression, time, SSL, command line parameter, data processing, subprocess call, XML parsing and date and time.

[0016] Among them, in the step of downloading the task:

[0017] Initialize and set the download file list, cookie storage path, ASF OAuth information, and SSL context;

[0018] Obtain and verify user information;

[0019] Download files and handle errors and redirects;

[0020] Iterate through the list of files and download them.

[0021] Among them, in the step of confirming the synthetic aperture radar data range according to the acquired original data, generating the intensity map of the super main image, and generating the reference InSAR deformation field:

[0022] Carry out radiometric calibration on the downloaded SAR data according to the scope of the mining area and convert the SAR data into different formats;

[0023] Generate a time-space baseline pairing file and an intensity map of the super master image based on the time-space baseline, and obtain an interference map;

[0024] After phase unwrapping of the interference pattern, the unwrapping residual in the deformation field is corrected, and the unwrapping phase image is obtained through the atmospheric stratification effect;

[0025] The meta-information of the corrected unwrapped phase image is written into the tab file to generate the time series deformation field.

[0026] Among them, in the step of generating a time-space baseline pairing file and an intensity map of a super master image according to the time-space baseline, and obtaining an interference map:

[0027] Select a super master image according to the time-space baseline, and generate a time-space baseline pairing file and an intensity map of the super master image;

[0028] Based on the registration interference script, the SLC files were batch processed according to the baseline pairing files, the multi-view ratio was set to 10:2, and the interference map was obtained.

[0029] Among them, in the step of calculating the deformation value of the intensity map corresponding to the date to generate the deformation field image, and outputting the deformation field txt file to monitor the InSAR deformation field:

[0030] Select the paired data that meet the requested spatial and temporal baselines, and update the corresponding baseline map;

[0031] Process the new paired data from the baseline map until each interference product is geocoded;

[0032] The deformation value of the corresponding date of the calculated intensity map is used to generate a deformation field image;

[0033] Output the deformation field txt file of the intensity map corresponding to the date.

[0034] Among them, in the step of processing the new paired data in the baseline map until geocoding each interference product:

[0035] Geocoding processing includes registration and interference, phase unwrapping and correction, and geocoding.

[0036] Among them, in the steps of establishing an early warning mechanism and conducting real-time assessment of the accumulated risks on the surface of the mining area:

[0037] Preset thresholds and a risk assessment model are used to conduct real-time assessment of the cumulative surface risks in the mining area based on multi-dimensional factors.

[0038] Among them, in the step of real-time assessment of the cumulative risk of the surface of the mining area based on multi-dimensional factors:

[0039] Multidimensional factors include deformation rate, deformation range, and rainfall conditions.

[0040] Among them, in the step of real-time assessment of the cumulative risk of the surface of the mining area based on multi-dimensional factors:

[0041] If abnormal deformation is detected and the risk level exceeds the preset threshold, the early warning mechanism will be triggered immediately, and a warning and deformation analysis report will be sent to the management personnel.

[0042] The invention discloses an automatic monitoring and early warning method for accumulated risk on the surface of a mining area based on InSAR technology, which includes importing a database and downloading tasks to obtain raw data; confirming the data range of synthetic aperture radar according to the obtained raw data, generating an intensity map of a super main image, and generating a reference InSAR deformation field; calculating the deformation value of the intensity map corresponding to the date to generate a deformation field image, and outputting a deformation field txt file to monitor the InSAR deformation field; establishing an early warning mechanism to conduct real-time assessment of accumulated risk on the surface of a mining area; and through the above-mentioned method, not only is it based on InSAR technology and has advantages in large-scale, high-precision, and non-contact monitoring, but also by optimizing the data processing flow and the automation algorithm, it realizes efficient, accurate, and continuous monitoring of the surface deformation of the mining area, and builds a solid technical defense line and early warning barrier for safe mining in the mining area. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. 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 creative work.

[0044] Figure 1It is a flow chart of the steps of the method for automatic monitoring and early warning of surface cumulative risk in mining areas based on InSAR technology of the present invention.

[0045] Figure 2 It is a step flow chart of S100 of the present invention.

[0046] Figure 3 It is a step flow chart of S130 of the present invention.

[0047] Figure 4 It is a step flow chart of S200 of the present invention.

[0048] Figure 5 It is a step flow chart of S220 of the present invention.

[0049] Figure 6 It is a step flow chart of S300 of the present invention. DETAILED DESCRIPTION

[0050] Here, exemplary embodiments are described in detail, and examples thereof are shown in the accompanying drawings. When the following description refers to the drawings, unless otherwise indicated, the same numbers in different drawings represent the same or similar elements. The implementations described in the following exemplary embodiments do not represent all implementations consistent with the present application.

[0051] The terms used in this application are for the purpose of describing specific embodiments only and are not intended to limit this application. The singular forms of "a", "said" and "the" used in this application and the appended claims are also intended to include plural forms unless the context clearly indicates other meanings. It should also be understood that the term "and / or" used herein refers to and includes any or all possible combinations of one or more associated listed items.

[0052] It should be understood that although the terms first, second, third, etc. may be used in the present application to describe various information, these information should not be limited to these terms. These terms are only used to distinguish the same type of information from each other. For example, without departing from the scope of the present application, the first information may also be referred to as the second information, and similarly, the second information may also be referred to as the first information. Depending on the context, the word "if" as used herein may be interpreted as "at the time of" or "when" or "in response to determining".

[0053] See also Figure 1 to Figure 6 ,in Figure 1 This is a flow chart of the steps of the automatic monitoring and early warning method of surface cumulative risk in mining areas based on InSAR technology. Figure 2 This is the step flow chart of S100. Figure 3 This is the step flow chart of S130. Figure 4 This is the step flow chart of S200. Figure 5 This is the step flow chart of S220. Figure 6 It is a step flow chart of S300.

[0054] The present invention provides an automatic monitoring and early warning method for surface cumulative risk in mining areas based on InSAR technology, comprising the following steps:

[0055] S100: Import the database and download the task to obtain the original data;

[0056] S200: confirming the synthetic aperture radar data range according to the acquired raw data, generating an intensity map of the super main image, and generating a reference InSAR deformation field;

[0057] S300: Calculate the deformation value of the intensity map corresponding to the date to generate a deformation field image, and output the deformation field txt file to monitor the InSAR deformation field;

[0058] S400: Establish an early warning mechanism to conduct real-time assessment of the accumulated risks on the surface of the mining area.

[0059] In this implementation, firstly, the database is imported and the task is downloaded to obtain the original data; then, the synthetic aperture radar data range is confirmed according to the obtained original data, and the intensity map of the super main image is generated to generate the reference InSAR deformation field; then, the deformation value of the intensity map corresponding to the date is calculated to generate the deformation field image, and the deformation field txt file is output to monitor the InSAR deformation field; finally, an early warning mechanism is established to conduct real-time assessment of the cumulative risk of the surface of the mining area; through the above method, not only based on InSAR technology, it has advantages in large-scale, high-precision, and non-contact monitoring, but also by optimizing the data processing process and the automation algorithm, it realizes efficient, accurate, and continuous monitoring of the surface deformation of the mining area, and builds a solid technical defense line and early warning barrier for the safe mining of the mining area.

[0060] Furthermore, in the step of importing the database and downloading the task to obtain the original data:

[0061] S110: according to the scope of the mining area, downloading the digital elevation model of the mining area;

[0062] S120: Preset a time span, and download synthetic aperture radar data according to the time span;

[0063] S130: Importing a database into the script and downloading a task, the downloaded task is original data; wherein the data in the database includes processing file paths, temporary files, regular expressions, time, SSL, command line parameters, data processing, subprocess calls, XML parsing and date and time.

[0064] Furthermore, in the step of downloading the task:

[0065] S131: Initialize and set the download file list, cookie storage path, ASF OAuth information, and SSL context;

[0066] S132: Obtain and verify user information;

[0067] S133: download files and handle errors and redirection data;

[0068] S134: traverse the file list and download.

[0069] In this embodiment, efficient batch acquisition of raw data is performed; the specific process is: first, according to the scope of the mining area, the latest published digital elevation model (DEM) with a resolution of 30 meters is downloaded for the mining area, and it is used as a reference DEM for subsequent deformation analysis.

[0070] The download URL is: https: / / dataspace.copernicus.eu / .

[0071] Then, a time span is preset, and a Python script is used to batch download synthetic aperture radar data from the Alaska Satellite Facility (ASF) website according to the time span. Various databases are imported into the Python script, and tasks are downloaded, and the downloaded tasks are raw data; the data in the database includes processing file paths, temporary files, regular expressions, time, SSL, command line parameters, data processing, subprocess calls, XML parsing, and date and time.

[0072] The core of the Python script is to handle download tasks, including: (1) initialization, setting the download file list, cookie storage path, ASF OAuth information, and SSL context; it also processes command line parameters and supports reading download links from Metalink or CSV files. (2) Obtain and verify the user's Earthdata login cookie for subsequent download operations. (3) Responsible for using cookies to download files and handle various possible errors and redirection data. (4) Traverse the file list and download. Through the above download tasks, not only the efficiency of data collection is greatly improved, but also the timeliness and consistency of the data are ensured.

[0073] Further, in the step of confirming the synthetic aperture radar data range according to the acquired raw data, generating the intensity map of the super main image, and generating the reference InSAR deformation field:

[0074] S210: performing radiometric calibration on the downloaded synthetic aperture radar data according to the scope of the mining area, and converting the synthetic aperture radar data into a format;

[0075] S220: Generate a time-space baseline pairing file and an intensity map of the super master image according to the time-space baseline, and obtain an interference map;

[0076] S230: After performing phase unwrapping on the interference pattern, correcting the unwrapping residual in the deformation field, and obtaining an unwrapped phase image through the atmospheric stratification effect;

[0077] S240: Write the meta-information of the corrected unwrapped phase image into a tab file to generate a time series deformation field.

[0078] Furthermore, in the step of generating a time-space baseline pairing file and an intensity map of a super master image according to the time-space baseline and obtaining an interference map:

[0079] S221: selecting a super master image according to the time-space baseline, and generating a time-space baseline pairing file and an intensity map of the super master image;

[0080] S222: Based on the registration interference script, batch process the SLC files according to the baseline pairing files, set the multi-view ratio to 10:2, and obtain the interference map.

[0081] In this embodiment, a reference InSAR deformation field is accurately generated, wherein the specific steps are: preprocessing the raw data, generating an interference pattern, unwrapping and correcting the phase, and producing a time-series InSAR result.

[0082] The raw data is preprocessed: the scope of SAR data to be processed is determined according to the scope of the mining area. Therefore, it is necessary to manually edit the reference burst file according to the actual situation, and then perform radiometric calibration on the downloaded SAR image data based on the extract_S1_from_zipfiles.sh script, and convert the received SAR signal into the standard format SLC that can be processed by the GAMMA software.

[0083] Generate interferogram: First, select_master.sh is used to select the most suitable super master image according to the time-space baseline and generate the time-space baseline pairing file (bperp_temp.list). Generate the intensity map (*.mli and *.mli.par) of the super master image; then, based on the intensity map (*.mli and *.mli.par) and the digital elevation model (DEM), after reference frame and format conversion, it is used as the basic geographic reference for subsequent deformation analysis to ensure the accuracy and reliability of deformation measurement. Use the registration interferometer script TOPS_Coreg_InSAR.sh to batch process SLC files according to the baseline pairing file (bperp_temp.list), with a multi-view ratio of 10:2. This step will obtain a large number of interferograms (*.diff).

[0084] Unwrapping phase and correction: Since the phase in the interferogram is entangled, that is, the phase value is limited to between -π and π, a phase unwrapping operation is required to recover the true phase change. Use mcf_with_multi_cpx.sh to perform phase unwrapping of the interferogram. This script implements the unwrapping of the complex-valued interferogram using unw_model based on the minimum cost rate method to further optimize the unwrapped complex-valued interferogram. The algorithm unwraps the complex-valued interferogram using a pre-existing unwrapping phase model. For each pixel, the unwrapped phase value is set to a value that is in the interval (-PI, PI) in the model and is consistent with the complex-valued interferogram, that is, re-wrapping the unwrapped interferogram will produce the phase value of the original interferogram, except for a constant phase offset from the reference value that sets the unwrapped phase to the reference position. The prerequisite for this model is that the deviation of the model from the true unwrapped phase is only a small value. For cases where the deviation is greater than PI, the unwrapped phase will be set to an incorrect ambiguity, that is, the model will deviate from the correct unwrapping result by a multiple of 2*PI. Then, the unw_correction_fit command is used to correct the unwrapped residuals in the deformation field again; the mk_atm_2d command is used to correct the atmospheric stratification effect in the phase, and finally a large number of unwrapped phase maps (*.unw) are obtained. In this process, in order to ensure the reliability of the reference InSAR deformation field, the pixels with coherence less than 0.2 are removed.

[0085] Output time series InSAR results: Use the mk_tab command to write all the corrected unwrapped phase map (*.unw) metadata into a tab file, and then use the mb command to process and generate the time series deformation field. Use disp_map to convert the phase signal into a deformation signal. Then, geocoding is performed to obtain the reference InSAR deformation field in the latitude and longitude reference frame. Finally, the processed deformation field data is output in an appropriate format, such as GeoTIFF, txt, NetCDF, etc. In order to carry out dynamic monitoring and abnormal warning in the future, the deformation field data is saved in txt format here.

[0086] Further, in the step of calculating the deformation value of the intensity map corresponding to the date to generate the deformation field image, and outputting the deformation field txt file to monitor the InSAR deformation field:

[0087] S301: Select paired data that meets the requested spatial and temporal baselines, and update the corresponding baseline map;

[0088] S302: Processing the new paired data in the baseline map until geocoding is performed on each interference product; wherein the geocoding process includes registration and interference, unwrapping phase and correction, and geocoding;

[0089] S303: generating a deformation field image using the deformation value of the calculated intensity map corresponding to the date;

[0090] S304: Output the deformation field txt file of the date corresponding to the intensity map.

[0091] In this embodiment, after the reference InSAR deformation field generation is realized, the present invention integrates the second step InSAR deformation field processing shell script and command to realize the automatic processing of SAR data preprocessing, registration, interference, unwrapping, correction, extraction, etc. The following process is proposed to realize the automatic dynamic update of the interference baseline map: after the image is updated, the script update_baseline_master.sh selects the paired data that meets the requested spatial and temporal baseline, and then updates the corresponding baseline map (bperp_temp.list). Ensure that each image SLC forms an interference map with the four previous and later SLCs, and all SLCs will be registered and sampled under the geometric framework of the super master image. Then execute each and only new paired data in the baseline map bperp_temp.list for processing until each interference product is geocoded: registration and interference, unwrapping phase and correction, geocoding. All InSAR products are geocoded on the grid defined in the parameter file so that all satellite geometries are similar. The deformation value corresponding to the date of the calculated updated image is used to generate a deformation field image, which can intuitively display the spatial distribution and change trend of the surface deformation in the mining area. And output the deformation field txt file of the corresponding date of the updated image for use in the next step of dynamic warning.

[0092] Furthermore, in the step of establishing an early warning mechanism and conducting real-time assessment of the cumulative risks on the surface of the mining area:

[0093] A preset threshold is used, and a risk assessment model is adopted to conduct a real-time assessment of the cumulative risk of the surface of the mining area based on multi-dimensional factors; the multi-dimensional factors include deformation rate, deformation range, and rainfall conditions.

[0094] Furthermore, in the step of real-time assessment of the cumulative risk of the surface of the mining area based on multi-dimensional factors:

[0095] If abnormal deformation is detected and the risk level exceeds the preset threshold, the early warning mechanism will be triggered immediately, and a warning and deformation analysis report will be sent to the management personnel.

[0096] In this embodiment, after the automatic processing of the deformation field is completed, the real-time risk warning and decision support stage is entered. The risk assessment model is adopted to comprehensively consider multi-dimensional factors such as deformation rate, deformation range, rainfall conditions, etc., and the cumulative risk of the surface of the mining area is evaluated in real time. Once the deformation abnormality is detected or the risk level exceeds the preset threshold, the early warning mechanism is immediately triggered, and early warning prompts are issued to managers through various methods such as graphical interface, email, and text messages, and a detailed deformation analysis report is attached to provide urgent data support for safe mining in the mining area.

[0097] The present invention has the following advantages:

[0098] 1. Large-scale and high-efficiency monitoring:

[0099] Compared with traditional ground survey and GPS technology, InSAR technology has a wider monitoring coverage, which can achieve comprehensive coverage of mining areas and surrounding areas. This large-scale monitoring capability helps to capture more comprehensive deformation information and provide more comprehensive data support for mining safety assessment. In addition, through Python scripts and optimized data processing processes, efficient batch acquisition and rapid processing of Sentinel-1SAR data are achieved. It not only simplifies the cumbersome manual operations in traditional data processing, but also significantly improves data processing efficiency, making the acquisition of surface deformation information in mining areas more timely and accurate.

[0100] (II) Reliability and scalability:

[0101] The present invention adopts phase correction and optimization processing in the process of generating the reference InSAR deformation field. At the same time, by setting a strict coherence threshold, low-quality data points are eliminated, further enhancing the stability and reliability of the deformation field. In addition, the technical solution of the present invention has high flexibility and scalability. By adjusting the parameters and configurations in the Python script, it can easily adapt to the deformation monitoring needs of different mining areas and different time spans. In addition, with the continuous development and improvement of InSAR technology, the present invention can further integrate more advanced functions and algorithms to improve monitoring performance and early warning accuracy.

[0102] (III) High degree of automation:

[0103] The present invention realizes the full chain of automated processing from data preprocessing to deformation field output, including the automated dynamic update of the interference baseline map. This highly automated processing mechanism greatly reduces the burden of manual intervention, improves the continuity and real-time nature of monitoring work, and ensures the consistency of monitoring results.

[0104] (IV) Fast early warning response:

[0105] The present invention has a built-in advanced risk assessment model that can comprehensively consider multi-dimensional risk factors and conduct real-time assessment of the cumulative risk of the surface of the mining area. Once abnormal deformation is detected or the risk level exceeds the preset threshold, the early warning mechanism is immediately activated and early warning information is sent to management personnel through multiple channels. This rapid early warning response capability provides timely early warning signals and decision-making support for safe mining in the mining area, which helps to take timely measures to avoid or reduce potential geological disaster risks.

[0106] In summary, the advantages of the present invention are not only reflected in the improvement of technological innovation and data processing efficiency, but also in the comprehensive benefits it brings to safe mining and ecological environment protection in mining areas. Through efficient and accurate deformation monitoring and real-time early warning mechanism, the present invention helps to reduce the probability and degree of geological disasters in mining areas, and ensure the production safety of mining areas and the safety of life and property of surrounding residents. At the same time, it also provides strong technical support and guarantee for the sustainable development of mining areas.

[0107] Those skilled in the art will readily appreciate other embodiments of the present application after considering the specification and practicing the contents disclosed herein. The present application is intended to cover any variations, uses or adaptations of the present application, which follow the general principles of the present application and include common knowledge or customary technical means in the art that are not disclosed in the present application.

[0108] It should be understood that the present application is not limited to the exact construction that has been described above and shown in the drawings, and that various modifications and changes may be made without departing from the scope thereof.

Claims

1. An automatic monitoring and early warning method for surface cumulative risk in mining areas based on InSAR technology, characterized in that: The steps include: Import the database and download the task to obtain the original data; Confirm the synthetic aperture radar data range based on the acquired raw data, generate the intensity map of the super master image, and generate the reference InSAR deformation field; Calculate the deformation value of the intensity map corresponding to the date to generate the deformation field image, and output the deformation field txt file to monitor the InSAR deformation field; Establish an early warning mechanism to conduct real-time assessment of the accumulated risks on the surface of the mining area.

2. The method for automatic monitoring and early warning of surface cumulative risk in mining areas based on InSAR technology as claimed in claim 1, characterized in that: In the steps of importing the database and downloading the task to obtain the original data: According to the scope of the mining area, download the digital elevation model of the mining area; Preset time span and download synthetic aperture radar data according to the time span; Import the database into the script and download the task, the downloaded task is the original data; the data in the database includes processing file path, temporary file, regular expression, time, SSL, command line parameter, data processing, subprocess call, XML parsing and date and time.

3. The method for automatic monitoring and early warning of surface cumulative risk in mining areas based on InSAR technology as claimed in claim 2, characterized in that: In the steps of downloading tasks: Initialize and set the download file list, cookie storage path, ASF OAuth information, and SSL context; Obtain and verify user information; Download files and handle errors and redirects; Iterate through the list of files and download them.

4. The method for automatic monitoring and early warning of surface cumulative risk in mining areas based on InSAR technology as claimed in claim 1, characterized in that: In the steps of confirming the synthetic aperture radar data range based on the acquired raw data, generating the intensity map of the super master image, and generating the reference InSAR deformation field: Carry out radiometric calibration on the downloaded SAR data according to the scope of the mining area and convert the SAR data into different formats; Generate a time-space baseline pairing file and an intensity map of the super master image based on the time-space baseline, and obtain an interference map; After phase unwrapping of the interference pattern, the unwrapping residual in the deformation field is corrected, and the unwrapping phase image is obtained through the atmospheric stratification effect; The meta-information of the corrected unwrapped phase image is written into the tab file to generate the time series deformation field.

5. The method for automatic monitoring and early warning of surface cumulative risk in mining areas based on InSAR technology as claimed in claim 4, characterized in that: In the step of generating a time-space baseline pairing file and an intensity map of a super master image according to the time-space baseline and obtaining an interference map: Select a super master image according to the time-space baseline, and generate a time-space baseline pairing file and an intensity map of the super master image; Based on the registration interference script, the SLC files were batch processed according to the baseline pairing files, the multi-view ratio was set to 10:2, and the interference map was obtained.

6. The method for automatic monitoring and early warning of surface cumulative risk in mining areas based on InSAR technology as claimed in claim 1, characterized in that: In the steps of calculating the deformation value of the intensity map corresponding to the date to generate the deformation field image, and outputting the deformation field txt file to monitor the InSAR deformation field: Select the paired data that meet the requested spatial and temporal baselines, and update the corresponding baseline map; Process the new paired data from the baseline map until each interference product is geocoded; The deformation value of the corresponding date of the calculated intensity map is used to generate a deformation field image; Output the deformation field txt file of the intensity map corresponding to the date.

7. The method for automatic monitoring and early warning of surface cumulative risk in mining areas based on InSAR technology as claimed in claim 6, characterized in that: In the steps that process the new paired data from the baseline map until each interferometric product is geocoded: Geocoding processing includes registration and interference, phase unwrapping and correction, and geocoding.

8. The method for automatic monitoring and early warning of surface cumulative risk in mining areas based on InSAR technology as claimed in claim 1, characterized in that: In the steps of establishing an early warning mechanism and conducting real-time assessment of the accumulated risks on the surface of the mining area: Preset thresholds and a risk assessment model are used to conduct real-time assessment of the cumulative surface risks in the mining area based on multi-dimensional factors.

9. The method for automatic monitoring and early warning of surface cumulative risk in mining areas based on InSAR technology as claimed in claim 8, characterized in that: In the steps of real-time assessment of the cumulative risk of the surface of the mining area based on multi-dimensional factors: Multidimensional factors include deformation rate, deformation range, and rainfall conditions.

10. The method for automatic monitoring and early warning of surface cumulative risk in mining areas based on InSAR technology as claimed in claim 8, characterized in that: In the steps of real-time assessment of the cumulative risk of the surface of the mining area based on multi-dimensional factors: If abnormal deformation is detected and the risk level exceeds the preset threshold, the early warning mechanism will be triggered immediately, and a warning and deformation analysis report will be sent to the management personnel.

Images