CN112097733A - Surface deformation monitoring method combining InSAR technology and geographic detector - Google Patents
Surface deformation monitoring method combining InSAR technology and geographic detector Download PDFInfo
- Publication number
- CN112097733A CN112097733A CN202010733973.1A CN202010733973A CN112097733A CN 112097733 A CN112097733 A CN 112097733A CN 202010733973 A CN202010733973 A CN 202010733973A CN 112097733 A CN112097733 A CN 112097733A
- Authority
- CN
- China
- Prior art keywords
- surface deformation
- factor
- driving
- driving force
- factors
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Geophysics And Detection Of Objects (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
The invention discloses a surface deformation monitoring method combining an InSAR technology and a geographic detector, which mainly comprises the following steps: firstly, acquiring a surface deformation rate by utilizing an InSAR technology; secondly, constructing grid surface deformation driving factors and random point attribute information; and finally, detecting driving force by utilizing single-factor and multi-factor interaction of the geographic detector. The invention converts the qualitative analysis of the original earth surface deformation driving force into the quantitative analysis, and provides powerful data support and a new idea for disaster prevention and reduction.
Description
Technical Field
The invention relates to the field of surface deformation monitoring of InSAR technology, in particular to a surface deformation monitoring method combining the InSAR technology and a geographic detector.
Background
Ground subsidence caused by factors such as underground water mining and tectonic movement has become an environmental geological disaster commonly developed in cities. The non-uniform settlement can cause damage to infrastructure and large-scale engineering in different degrees, and especially threatens the stability of a high-speed railway with extremely strict requirements on the stability of a roadbed structure, the smoothness of continuous bridges and rails, so that the monitoring of the deformation of the earth surface is very important. The traditional earth surface deformation monitoring modes (GPS, level and the like) have long period and low monitoring density, and meanwhile, the analysis of the earth surface deformation mechanism is lack of quantitative analysis. Therefore, a large-range and high-precision ground surface deformation monitoring means is needed to effectively monitor ground surface settlement, a quantitative driving force analysis method is needed to analyze the mechanism of the ground surface settlement, and the method has great significance for urban disaster prevention and reduction, risk prevention and the like.
Synthetic Aperture Radar interferometry (InSAR) is a remote sensing technology capable of obtaining surface centimeter or even millimeter-scale deformation, and has the advantages of all weather, coverage area, high precision, low cost and the like. The InSAR technology is widely applied to deformation monitoring of surface deformation monitoring, volcano, landslide, underground mining and the like, and has good application foundation and application prospect.
A Geo-Detector (Geo-Detector) is capable of quantitatively detecting spatial dissimilarity to reveal statistical methods of its driving force. The core idea of the method is as follows: if a certain factor has a significant spatial impact on the occurrence of a certain phenomenon, then the factor should have similarity to the spatial distribution of the phenomenon. The method can not only detect the driving force of a single factor, but also judge the strength of the interaction of multiple factors.
The invention utilizes InSAR technology to monitor the deformation of the earth surface and puts forward the problem of applying a geographic detector to quantitatively analyze the mechanism of the deformation of the earth surface. Therefore, the earth surface deformation monitoring method combining the InSAR technology and the geographic detector is provided.
In the prior art, the following documents are mainly relevant to the present application:
document 1 discloses an invention of a method and a system for splicing time series data of ground settlement monitored by adaptive multi-source InSAR, applied to chenbei et al, university of capital, application No.: 201911018550.5. the method combines external data to unify spatial reference, reduces system errors, constructs an adaptive time sequence model, avoids introduction of common fitting model errors, improves accuracy of spliced data, has strong applicability to large-scale splicing of ground settlement time sequence data of multi-source InSAR monitoring results, and obviously improves accuracy of the spliced data.
Document 2 invention, entitled "a method for detecting surface subsidence" applied to geodetic sea-land geotechnical engineering development ltd, zhangbaohua, of Tianjin, application No.: 201511034902.8. the method utilizes a total station, a computer and the like to measure the ground settlement, avoids personnel launching operation and improves the working efficiency.
Document 3, wangquan quails of geological science research institute of Shandong province, an invention of a method for constructing a ground subsidence GPS automatic monitoring station, application number: 201911057119.1. according to the invention, high-precision observation can be carried out on ground settlement by building one GPS reference station and six GPS monitoring stations.
Document 5 invention of the present application, entitled "a driving force evaluation method for mountain torrent disaster distribution in super large area", applied by Tianjin university school park, application number: 201911125521.9. according to the method, various factors influencing landslide are firstly constructed, then the factors are detected by a geographic detector, and the driving factors of the space-time distribution pattern of the mountain torrent disaster are quantitatively analyzed, so that the influence of each driving factor is more reliable, and the evaluation result is more real.
Document 6, a geographic detector final discretization method based on simulated annealing, applied by the institute of geoscience and resource research, xu Dong, of the Chinese academy of sciences, application No.: 201810678269.3. the method comprises the steps of firstly initially calculating the geographic detector q, then optimizing the breakpoint vector by using a simulated annealing algorithm, and obtaining the optimal breakpoint vector and the maximum geographic detector q statistic corresponding to the optimal breakpoint vector after optimization.
The above patent applications do not relate to quantitative analysis of combination of surface deformation and driving force, so the invention provides a surface deformation monitoring method combining InSAR technology and a geographic detector, realizes single-factor and multi-factor interactive quantitative analysis of a surface deformation mechanism, and develops qualitative analysis into quantitative analysis.
Disclosure of Invention
The purpose of the invention is as follows: the invention aims to solve the quantification of the large-range ground surface deformation mechanism explanation, and provides a ground surface deformation monitoring method combining an InSAR technology and a geographic detector, so as to improve the ground surface deformation monitoring range, quantitatively analyze the single-factor and multi-factor driving force of ground surface deformation, and provide scientific basis for urban disaster prevention and reduction and ground surface settlement risk prevention.
The technical scheme is as follows: the invention relates to a surface deformation monitoring method combining InSAR technology and a geographic detector, which sequentially comprises the following steps:
(1) and (2) monitoring the surface deformation of the research area in a large range by utilizing InSAR technologies (including time series InSAR technologies such as PS-InSAR and SBAS-InSAR) to obtain the surface deformation rate and the time series deformation quantity of the research area.
(2) And constructing the space-time distribution of the surface deformation driving factors, including natural factors (geological faults, elevations, gradients, underground water mining, precipitation, air temperature and the like) and social factors (population density, building density, distance from a railway, distance from a road, land utilization types and the like), and preprocessing.
(3) And dividing random points of the fishing net production, taking point density and processing efficiency into consideration when selecting the random points, selecting a proper number of points, extracting driving factors and values of sedimentation rate, and exporting the values into a table.
(4) The driving force analysis is carried out on the earth surface deformation rate by utilizing a single-factor and multi-factor interaction method of the geographic detector, the earth surface deformation rate is used as a dependent variable, the constructed driving factors are used as independent variables, and the data are imported into geographic detector software for analysis and processing to obtain a single-factor driving force and a multi-factor interaction driving force, so that earth surface deformation mechanism parameters are revealed quantitatively, and the influence factors of the earth surface deformation of the research area are explained powerfully.
The invention has the beneficial effects that: the surface deformation monitoring method combining the InSAR technology and the geographic detector disclosed by the invention can monitor the surface subsidence in a large range and high precision through the InSAR technology, can effectively reduce the economic investment of surface subsidence monitoring, can quantitatively research the driving force of the surface subsidence by analyzing the driving force through the geographic detector, and provides theoretical and technical support for preventing the surface subsidence disaster.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Fig. 2 is a plot of the amount of surface subsidence monitored by InSAR technology.
Fig. 3 is a grid and random point plot.
Detailed Description
An earth surface deformation monitoring method combining InSAR technology and a geographic detector comprises the following implementation steps:
the method comprises the following steps: and acquiring the surface deformation rate of the research area by using an SBAS-InSAR technology.
(1) Data preprocessing: the 32 scene Sentinel-1A image was imported into ENVI and cropped according to the study area.
(2) And (3) generating a connection graph: the time baseline is 200 days, the space baseline is set to be 45% of the critical baseline, and the image with the optimal Doppler centroid frequency at 2015, 9 months and 22 days is selected as the super main image.
(3) Differential interference pattern generation: according to the resolution and the incidence angle of the azimuth direction and the distance direction of the image, the vision is calculated to be 1: and 4, selecting a Delaunay MCF as a phase unwrapping method, and screening the interferogram by using a baseline tool to remove the interferogram with poor quality.
(4) Track refining and re-leveling: selecting a certain number of control points with high coherence in a research area according to certain priori knowledge, estimating residual phases in the initial unwrapping phases according to phase information of the control points, and removing residual terrain phases.
(5) And (3) SBAS inversion: in the first inversion, a linear model is used for estimating deformation rate and residual topography, secondary unwrapping is carried out, and the input interference pattern is optimized so as to be convenient for next processing. In the second inversion, low-pass filtering and high-pass filtering are used for removing the atmospheric phase, and the final deformation rate and the deformation quantity are obtained.
(6) And (3) geocoding: and converting the result of the SBAS inversion into a geographic coordinate system by a coherence coefficient of 0.85.
(7) And (4) visualizing the result: and (5) outputting the vector file subjected to geocoding by drawing with ArcGIS software.
Step two: and constructing a ground surface deformation driving factor spatial distribution map.
The factors influencing the ground surface settlement are various, factors such as underground water, geological faults, land coverage types, built-up areas, roads, air temperature and rainfall are researched to be the reasons of ground surface deformation according to the existing documents, and the geological faults, the land coverage types, the built-up areas, the roads, the air temperature and the rainfall are finally constructed to be used as a ground surface settlement driving factor spatial distribution map in consideration of data sources.
Step three: and extracting the driving factors and the surface deformation attributes of the random points by using the fishing net.
And dividing the fishing net, generating random points in the fishing net, wherein the grid size is according to the size of a research area, and the city can be divided into grids of 500m-2000m, and provinces can be divided into kilometer-level grids. The number of points should be considered in combination with accuracy and computational complexity.
In this embodiment, a grid of 1 × 1km is used, and 30 random points are generated in each grid, resulting in 34 grids (fig. 2) and 1020 random points. When the driving factor is numerical, the reclassification is performed, and generally, a reclassification method in ArcGIS is used. And extracting driving factors and sedimentation rate attribute values by using the grids, eliminating points without attributes, finally generating 812 effective points, and transferring the attributes of the points to Excel.
Step four: and (4) carrying out driving force analysis on the deformation rate of the earth surface by using a geographic detector.
The driving force of the single factor of the earth surface deformation rate is detected by the differentiation and factor detection module of the geographic detector (table 1), and the driving force of different factor combinations of the earth surface deformation rate is detected by the interaction detection module (table 2).
TABLE 1 Single factor detection results
TABLE 2 factor combination drive analysis
Claims (5)
- The earth surface deformation monitoring method combining the InSAR technology and the geographic detector is characterized by comprising the following steps: the method comprises the steps of obtaining an SAR image of a research area, constructing a ground surface deformation driving factor spatial distribution map of the research area, dividing fishing net production random points, and analyzing driving force of ground surface deformation rate by using a single-factor and multi-factor interaction method of a geographic detector.
- 2. The method of claim 1, wherein the InSAR technique is used to obtain the earth deformation rate of the research area by the steps of connectivity map generation, interferometry, deplature, phase unwrapping, inversion, and geocoding.
- 3. The method for constructing the spatial distribution map of the surface deformation driving factors of the research area according to claim 1, wherein the spatial distribution map of the surface deformation driving factors of the research area comprises natural factors (geological faults, elevation, gradient, underground water mining, precipitation, air temperature and the like) and social factors (population density, building density, distance to railway, distance to road, land utilization type and the like).
- 4. The divided fishing net production random point according to claim 1, wherein the selection of the random point should be considered for point density and processing efficiency, and an appropriate number of points are selected to extract the properties of the deformation rate and the driving factor according to claims 2 and 3 and are exported as a table.
- 5. The method for driving force analysis of earth surface deformation rate by using single-factor and multi-factor interaction method of geographic detector as claimed in claim 1 is characterized in that according to claims 2, 3 and 4, earth surface deformation rate is used as dependent variable, constructed driving factor is used as independent variable, and the deformation rate and the attribute data of the driving factor are introduced into geographic detector software for analysis and processing to obtain single-factor driving force and multi-factor interaction driving force, so that earth surface deformation mechanism parameters are revealed quantitatively, and influence factors of earth surface deformation in a research area are explained powerfully.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010733973.1A CN112097733A (en) | 2020-07-28 | 2020-07-28 | Surface deformation monitoring method combining InSAR technology and geographic detector |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010733973.1A CN112097733A (en) | 2020-07-28 | 2020-07-28 | Surface deformation monitoring method combining InSAR technology and geographic detector |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112097733A true CN112097733A (en) | 2020-12-18 |
Family
ID=73750508
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010733973.1A Pending CN112097733A (en) | 2020-07-28 | 2020-07-28 | Surface deformation monitoring method combining InSAR technology and geographic detector |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112097733A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113251914A (en) * | 2021-04-06 | 2021-08-13 | 兰州交通大学 | Surface deformation prediction method combining InSAR technology and long-term memory neural network |
| CN113281742A (en) * | 2021-06-02 | 2021-08-20 | 西南交通大学 | SAR landslide early warning method based on landslide deformation information and meteorological data |
| CN113343563A (en) * | 2021-05-27 | 2021-09-03 | 中交第二公路勘察设计研究院有限公司 | Landslide susceptibility evaluation method based on automatic sample selection and surface deformation rate |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101770027A (en) * | 2010-02-05 | 2010-07-07 | 河海大学 | Ground surface three-dimensional deformation monitoring method based on InSAR and GPS data fusion |
| CN102938095A (en) * | 2012-11-23 | 2013-02-20 | 湖南科技大学 | Mining subsidence monitoring and early warning method based on multi-source data |
| CN108801211A (en) * | 2018-04-08 | 2018-11-13 | 山东农业大学 | A kind of sinking land in coalmining areas INSAR image method for quickly identifying |
| CN110109105A (en) * | 2019-04-23 | 2019-08-09 | 中国电力科学研究院有限公司 | A method of the InSAR technical monitoring Ground Deformation based on timing |
| CN110888130A (en) * | 2019-10-30 | 2020-03-17 | 华东师范大学 | Coal mine area ground surface deformation monitoring method based on lifting rail time sequence InSAR |
-
2020
- 2020-07-28 CN CN202010733973.1A patent/CN112097733A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101770027A (en) * | 2010-02-05 | 2010-07-07 | 河海大学 | Ground surface three-dimensional deformation monitoring method based on InSAR and GPS data fusion |
| CN102938095A (en) * | 2012-11-23 | 2013-02-20 | 湖南科技大学 | Mining subsidence monitoring and early warning method based on multi-source data |
| CN108801211A (en) * | 2018-04-08 | 2018-11-13 | 山东农业大学 | A kind of sinking land in coalmining areas INSAR image method for quickly identifying |
| CN110109105A (en) * | 2019-04-23 | 2019-08-09 | 中国电力科学研究院有限公司 | A method of the InSAR technical monitoring Ground Deformation based on timing |
| CN110888130A (en) * | 2019-10-30 | 2020-03-17 | 华东师范大学 | Coal mine area ground surface deformation monitoring method based on lifting rail time sequence InSAR |
Non-Patent Citations (1)
| Title |
|---|
| 石鹏远等: "应用地理探测器改进地面沉降危险性评估模型的研究", 《中国地质灾害与防治学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113251914A (en) * | 2021-04-06 | 2021-08-13 | 兰州交通大学 | Surface deformation prediction method combining InSAR technology and long-term memory neural network |
| CN113343563A (en) * | 2021-05-27 | 2021-09-03 | 中交第二公路勘察设计研究院有限公司 | Landslide susceptibility evaluation method based on automatic sample selection and surface deformation rate |
| CN113343563B (en) * | 2021-05-27 | 2022-05-03 | 中交第二公路勘察设计研究院有限公司 | Landslide susceptibility evaluation method based on automatic sample selection and surface deformation rate |
| CN113281742A (en) * | 2021-06-02 | 2021-08-20 | 西南交通大学 | SAR landslide early warning method based on landslide deformation information and meteorological data |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111142119A (en) | Mine geological disaster dynamic identification and monitoring method based on multi-source remote sensing data | |
| CN112097733A (en) | Surface deformation monitoring method combining InSAR technology and geographic detector | |
| Chen et al. | Investigating land subsidence and its causes along Beijing high-speed railway using multi-platform InSAR and a maximum entropy model | |
| Xue et al. | Novel SfM-DLT method for metro tunnel 3D reconstruction and Visualization | |
| CN104133996A (en) | Ground settlement risk grade evaluation method based on cloud model and data field | |
| CN107657618B (en) | Automatic extraction method of regional scale erosion gully based on remote sensing image and topographic data | |
| CN105627938A (en) | Pavement asphalt thickness detection method based on vehicle-mounted laser scanning spot cloud | |
| CN112882032B (en) | Method and device for dynamically monitoring geological disaster SAR in key area of gas pipeline | |
| CN114966689A (en) | Coastal city time sequence InSAR settlement monitoring analysis method, device, equipment and medium | |
| CN113251914A (en) | Surface deformation prediction method combining InSAR technology and long-term memory neural network | |
| CN103235845A (en) | Water and soil erosion monitoring method for highway construction | |
| Yang et al. | A PSI targets characterization approach to interpreting surface displacement signals: A case study of the Shanghai metro tunnels | |
| Hu et al. | Methods for monitoring fast and large gradient subsidence in coal mining areas using SAR images: A review | |
| Wang et al. | Analysis and prediction of regional land subsidence with InSAR technology and machine learning algorithm | |
| Guo et al. | Utilization of 3D laser scanning for stability evaluation and deformation monitoring of landslides | |
| CN114674277A (en) | Method for monitoring surface subsidence of deep mining and evaluating land reclamation suitability by combining whole-field survey lines | |
| Yang et al. | A new technical pathway for extracting high accuracy surface deformation information in coal mining areas using UAV LiDAR data: An example from the Yushen mining area in western China | |
| Janowski et al. | Assessment of cliff stability after the demolition of the engineering facilities | |
| Xiao et al. | A review of remote sensing applications in urban planning and management in China | |
| An et al. | Ground subsidence monitoring in based on UAV-LiDAR technology: a case study of a mine in the Ordos, China | |
| Gu et al. | Ground point extraction using self-adaptive-grid and point to surface comparison | |
| Xianghang et al. | Integrated space-air-ground early detection technologies and applicationfor potential landslide of transmission line corridor | |
| CN115287984B (en) | Road surface flatness evaluation method and system based on spaceborne InSAR | |
| CN116403116B (en) | D-InSAR earth surface deformation space-time characteristic fine sensing method for comprehensively planning multiple observation indexes | |
| Lato et al. | Predicting the Future by Mapping the Past: Revolutionary Innovations in Lidar Change Detection Analysis are Enabling Regional Scale Mapping and Identification of Threats From Geohazards |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20201218 |
|
| WD01 | Invention patent application deemed withdrawn after publication |