CN114777633A - Monitoring and analyzing method for closed mining area stage deformation - Google Patents
Monitoring and analyzing method for closed mining area stage deformation Download PDFInfo
- Publication number
- CN114777633A CN114777633A CN202210337255.1A CN202210337255A CN114777633A CN 114777633 A CN114777633 A CN 114777633A CN 202210337255 A CN202210337255 A CN 202210337255A CN 114777633 A CN114777633 A CN 114777633A
- Authority
- CN
- China
- Prior art keywords
- deformation
- closed
- mining area
- mine
- stage
- 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.)
- Granted
Links
- 238000005065 mining Methods 0.000 title claims abstract description 119
- 238000000034 method Methods 0.000 title claims abstract description 82
- 238000012544 monitoring process Methods 0.000 title claims abstract description 32
- 230000007423 decrease Effects 0.000 claims abstract description 34
- 238000012545 processing Methods 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 238000005305 interferometry Methods 0.000 claims description 9
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 230000002123 temporal effect Effects 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 abstract description 13
- 238000000691 measurement method Methods 0.000 abstract 1
- 239000011159 matrix material Substances 0.000 description 5
- 238000005562 fading Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009435 building construction Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000003911 water pollution Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/16—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring the deformation in a solid, e.g. by resistance strain gauge
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
Abstract
The invention provides a method for monitoring and analyzing stage deformation of a closed mining area, which comprises the following steps: acquiring a plurality of synthetic aperture radar images of a closed mining area within a preset time period; performing interference processing on the plurality of synthetic aperture radar images based on a small baseline set interference measurement method, and determining the surface deformation amount and deformation rate of the closed mining area in the preset time period; determining a surface movement phase of the closed mine area based on the surface deformation amount and the deformation rate, wherein the surface movement phase comprises a start phase, an active phase and a decline phase; determining stability of the shut-down mine area based on the surface movement phase. The method can refine the earth surface moving process of the closed mining area, and improves the reliability of the stability evaluation of the closed mining area earth surface and the scientificity of the earth surface moving process analysis.
Description
Technical Field
The invention relates to the technical field of radar remote sensing and image processing, in particular to a method for monitoring and analyzing stage deformation of a closed mining area.
Background
With the continuous development of the technology and the adjustment of the energy structure, coal mining gradually goes from original distributed mining to refined and centralized mining, and a large number of mining areas with laggard energy production face to be closed or abandoned. This shift also causes many problems such as water pollution and environmental pollution caused by closing the mine, and thus the subsequent treatment of the closed mine becomes an urgent problem to be solved.
In the prior art, deformation monitoring is carried out on a mining area generally by selecting the mining area which is being mined, then evaluation and prediction are carried out on mining risks by combining mining propulsion work, the mining safety and the mining orderliness of the mining area are ensured, deformation monitoring is rarely carried out on a closed mining area, actually, the stability of the ground surface of the closed mining area can be detected by carrying out deformation monitoring on the closed mining area, a lot of valuable information can be provided, and the method is crucial to the subsequent hazard evaluation of the closed mining area and the preparation of a processing scheme.
Based on the above problems, the technical scheme for monitoring the closed mining area is provided in the prior art, but in the analysis process of the closed mining area, the stability of the ground surface is evaluated through the integral deformation rate, and the influence of different ground surface movement stages of the closed mining area on the stability of the closed mining area is not considered, so that the technical problem that the evaluation of the stability of the ground surface of the closed mining area is not fine enough is caused.
Disclosure of Invention
In view of this, it is necessary to provide a method for monitoring and analyzing the deformation of the closed mining area in a stage manner, so as to solve the technical problem in the prior art that the evaluation of the stability of the ground surface of the closed mining area is not fine enough.
The invention provides a closed mining area stage deformation monitoring and analyzing method, which comprises the following steps:
acquiring a plurality of synthetic aperture radar images of a closed mining area within a preset time period;
performing interference processing on the plurality of synthetic aperture radar images based on a small baseline set interferometry method, and determining the surface deformation amount and deformation rate of the closed mining area in the preset time period;
determining a surface movement phase of the closed mine area based on the surface deformation amount and the deformation rate, wherein the surface movement phase comprises a start phase, an active phase and a decline phase;
and determining the stability of the closed mining area according to the surface movement stage.
In some possible implementations, the determining the surface movement phase of the closed mine area based on the surface deformation amount and the deformation rate includes:
when the deformation amount of the earth surface of the closed mining area is greater than or equal to the initial deformation amount and the deformation rate is less than the threshold deformation rate, the earth surface movement process of the closed mining area is the starting stage;
when the deformation rate is greater than or equal to a threshold deformation rate, the process of closing the earth surface of the mining area is the active stage;
and when the deformation rate is smaller than a threshold deformation rate, or the surface deformation amount is smaller than or equal to a threshold deformation amount in a threshold time period, the surface movement process of the closed mining area is in the decline stage.
In some possible implementations, the interference processing on the plurality of synthetic aperture radar images based on the small baseline set interferometry method to determine the amount and rate of surface deformation of the closed mine area within the preset time period includes:
dividing the plurality of synthetic aperture radar images into a plurality of small baseline sets based on a spatial baseline threshold and a temporal baseline threshold, wherein the small baseline sets comprise a plurality of synthetic aperture radar image pairs;
determining a deformation time sequence of each synthetic aperture radar image pair in the plurality of synthetic aperture radar image pairs based on a least square method;
and obtaining the surface deformation quantity and the deformation rate of the closed mining area in the preset time period based on a singular value decomposition method and the deformation time sequence.
In some possible implementations, the determining the stability of the closed mine area based on the surface deformation amount, the deformation rate and the preset staging rule includes:
judging whether the ground surface moving stage for closing the mining area is the decline stage or not based on the ground surface deformation quantity and the deformation rate;
if the stage of closing the ground surface movement of the mining area is the decline stage, the mining area is closed stably; and if the ground surface movement stage of the closed mining area is the starting stage or the active stage, the closed mining area is unstable.
In some possible implementations, the closed-mine-area stage deformation monitoring and analyzing method further includes:
determining a decline start time of the decline stage for the surface movement process of the closed mining area based on the surface deformation amount and the deformation rate;
determining a mine closing time for the closed mine based on the decline start time.
In some possible implementations, the method for monitoring and analyzing deformation in a closed mining area stage further includes:
obtaining the closing time of the mine area;
determining the decline starting time of the decline stage as the earth surface moving process of the closed mining area based on the earth surface deformation quantity, the deformation rate and the preset stage division rule;
determining whether the mine shut-down time is reasonable based on the decline start time.
In some possible implementations, the closed-mine-area stage deformation monitoring and analyzing method further includes:
obtaining various types of surface ground objects of the closed mining area;
constructing corresponding relations between potential damage levels of the various types of surface features and various stages of the surface movement process;
determining a target stage of the earth surface movement process based on the earth surface deformation quantity and the deformation rate;
and determining the potential damage level of the multiple types of surface ground objects based on the target stage and the corresponding relation.
In some possible implementations, after the interferometric processing the plurality of synthetic aperture radar images based on the small baseline set interferometry method determines the amount and rate of deformation of the surface of the closed mine area within the preset time period, the method further includes:
obtaining a reference surface deformation quantity of the closed mining area based on a global satellite navigation system leveling method;
converting the surface deformation into a vertical surface deformation;
and determining whether the precision of the surface deformation quantity meets the precision requirement or not according to the reference surface deformation quantity and the vertical surface deformation quantity.
In some possible implementations, the preset time period covers a mine shut-off time of the mine area.
The beneficial effects of adopting the above embodiment are: according to the monitoring and analyzing method for closed mining area staged deformation, the ground surface moving stage for closing the mining area is determined based on the ground surface deformation quantity and the deformation rate, and the ground surface moving process for closing the mining area is divided into the starting stage, the active stage and the decline stage based on the preset staging rule, so that the ground surface moving process for closing the mining area can be refined. And the stability of closing the mining area can be determined by dividing the regular ground surface moving stage based on the ground surface deformation quantity, the deformation rate and the preset stage, so that the reliability of evaluation on the ground surface stability of the closed mining area and the scientificity of ground surface moving process analysis can be improved.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
FIG. 1 is a schematic flow chart of an embodiment of a closed-seam staged deformation monitoring and analyzing method according to the present invention;
FIG. 2 is a schematic flow chart of another embodiment of S103 of FIG. 1 according to the present invention;
FIG. 3 is a schematic flow chart of another embodiment of S102 of FIG. 1;
FIG. 4 is a schematic flow chart of one embodiment of S104 of FIG. 1;
FIG. 5 is a schematic flow chart diagram illustrating one embodiment of determining a closing time for a mine area provided by the present invention;
FIG. 6 is a schematic flow chart illustrating an embodiment of the present invention for determining whether the closing time of a mine is reasonable;
FIG. 7 is a schematic flow chart illustrating one embodiment of the present invention for determining a level of surface damage closing a mine area;
fig. 8 is a schematic flow chart of an embodiment of verifying the accuracy of the surface deformation quantity provided by the present invention.
Detailed Description
The technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. It should be apparent that the described embodiments are only some embodiments of the present invention, and not all embodiments. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
It should be understood that the schematic drawings are not drawn to scale. The flowcharts used in this invention illustrate operations performed in accordance with some embodiments of the present invention. It should be understood that the operations of the flow diagrams may be performed out of order, and that steps without logical context may be reversed in order or performed concurrently. One skilled in the art, under the direction of this summary, may add one or more other operations to, or remove one or more operations from, the flowchart.
Some of the block diagrams shown in the figures are functional entities and do not necessarily correspond to physically or logically separate entities. These functional entities may be implemented in the form of software, or in one or more hardware modules or integrated circuits, or in different networks and/or processor systems and/or microcontroller systems.
Reference herein to "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment can be included in at least one embodiment of the invention. The appearances of the phrase in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. It is explicitly and implicitly understood by one skilled in the art that the embodiments described herein may be combined with other embodiments.
The embodiment of the invention provides a method for monitoring and analyzing the stage deformation of a closed mining area, which is respectively explained below.
Fig. 1 is a schematic flow chart of an embodiment of a method for monitoring and analyzing closed-mine-area periodic deformation according to the present invention, and as shown in fig. 1, the method for monitoring and analyzing closed-mine-area periodic deformation includes:
s101, obtaining a plurality of Synthetic Aperture Radar (SAR) images for closing a mining area within a preset time period;
s102, carrying out interference processing on a plurality of Synthetic Aperture Radar images based on a small baseline set interferometry (SBAS-InSAR) method, and determining the surface deformation amount and deformation rate of a closed mining area in a preset time period;
s103, determining an earth surface movement stage for closing the mining area based on the earth surface deformation quantity and the deformation rate, wherein the earth surface movement stage comprises a starting stage, an active stage and a decline stage;
and S104, determining the stability of closing the mining area according to the ground surface moving stage. Compared with the prior art, the monitoring and analyzing method for the closed mining area stage deformation provided by the embodiment of the invention has the advantages that the ground surface moving stage for closing the mining area is determined based on the ground surface deformation quantity and the deformation rate, and the ground surface moving process comprises the starting stage, the active stage and the decline stage, so that the ground surface moving process for closing the mining area can be refined. And the stability of the closed mining area is determined through the ground surface movable stage, so that the reliability of ground surface stability evaluation of the closed mining area and the scientificity of ground surface movement process analysis can be improved.
It should be understood that: the multiple synthetic aperture radar images are the same orbit data, the preset time period covers the closing time of the mining area, and the synthetic aperture radar images completely cover the closing mining area.
By setting the closing time of the mining area covered and closed by the preset time period, the synthetic aperture radar image completely covers and closes the mining area, and the reliability of analyzing the stability of the mining area can be ensured.
It should be noted that: in order to reduce the amount of computation and improve the efficiency of determining the amount of surface deformation and the deformation speed, in some embodiments of the present invention, after step S101, a plurality of synthetic aperture radar images may be further clipped, so as to reduce the data amount of the synthetic aperture radar images, thereby improving the efficiency of determining the amount of surface deformation and the deformation speed.
It should also be noted that: the synthetic aperture radar image in the embodiment of the invention is an image subjected to VV polarization.
In some embodiments of the present invention, as shown in fig. 2, step S103 comprises:
s201, when the deformation amount of the ground surface of the closed mining area is larger than or equal to the initial deformation amount and the deformation rate is smaller than a threshold deformation rate, closing the ground surface moving process of the mining area as a starting stage;
s202, when the deformation rate is larger than or equal to the threshold deformation rate, closing the earth surface movement process of the mining area to be an active stage;
and S203, when the deformation rate is smaller than the threshold deformation rate or the surface deformation amount is smaller than or equal to the threshold deformation amount in the threshold time period, closing the surface movement process of the mining area to be a decline stage.
In a specific embodiment of the invention, the initial amount of deformation is 10mm, the threshold amount of deformation is 30mm, and the threshold time period is six months.
It should be understood that: the initial deformation amount, the threshold deformation rate, the threshold time period and the threshold deformation amount can be adjusted according to specific mining areas, and are not described in detail herein.
In some embodiments of the present invention, as shown in fig. 3, step S102 includes:
s301, dividing the multiple synthetic aperture radar images into multiple small baseline sets based on a space baseline threshold value and a time baseline threshold value, wherein the small wiring sets comprise multiple synthetic aperture radar image pairs;
s302, determining a deformation time sequence of each synthetic aperture radar image pair in a plurality of synthetic aperture radar image pairs based on a least square method;
s303, obtaining the surface deformation quantity and the deformation rate of the closed mining area in a preset time period based on the singular value decomposition method and the deformation time sequence.
It should be understood that: the spatial baseline threshold and the temporal baseline threshold may be adjusted according to actual situations, and are not specifically limited herein.
The basic principle of the interferometry method based on the small baseline set is as follows:
let tAAnd tB(tAPrior to tB) Two synthetic aperture radar images at the moment generate a Kth interference image with the interference phase of
(k-1, 2, … M). Neglecting the influence of atmospheric and terrain phases in the azimuth-distance pixel coordinate system (x, r), the interference phase can be expressed as:
where λ is the radar wavelength, d (t)AX, r) and d (t)BX, r) are each tAAnd tBThe amount of deformation in the radar realization direction at the moment.
Let the master and slave images of each synthetic aperture radar image pair subjected to interference processing be sequentially expressed as: (m)1，s1)，…，(mk，sk)，…，(mM，sM) And satisfies:
mk＞sk k＝1，2，...M
all differential interference phases can then be composed as the following observation equation:
the above formula can be expressed in matrix form as:
Aφ＝δφ
in the formula, A is a coefficient matrix which is an M multiplied by N matrix, N +1 is the total number of the synthetic aperture radar images, and delta phi is each differential interference phase value; phi is the phase value of each synthetic aperture radar image.
The most ideal situation is that all data are grouped into one group, when M is more than N, the rank of the matrix A is N, and the estimated value of phi can be solved by adopting a least square method:
in the formula (I), the compound is shown in the specification,
is an estimate of phi.
And when the rank of the matrix A is less than N, obtaining the surface deformation quantity and the deformation rate by using a singular value decomposition method.
In some embodiments of the present invention, as shown in fig. 4, step S104 comprises:
s401, judging whether the ground surface movement stage of the closed mining area is a decline stage or not based on the ground surface deformation quantity and the deformation rate;
s402, if the stage of closing the ground surface movement of the mining area is a decline stage, closing the stability of the mining area; and if the ground surface movement stage of the closed mining area is a starting stage or an active stage, the closed mining area is unstable.
The method and the device determine whether the closed mining area is stable in real time based on the ground surface moving process where the closed mining area is located, can improve the fine evaluation of the stability of the closed mining area, and improve the reliability of the stability evaluation of the closed mining area.
Further, in some scenarios, it is desirable to obtain a mine closing time for closing the mine, and the mine closing time should be approximately the same as the fading start time of the fading phase, so in some embodiments of the present invention, as shown in fig. 5, the method for monitoring and analyzing deformation in the closed mine phase further includes:
s501, determining the fading start time of closing the ground surface movement process of the mining area to be a fading stage based on the ground surface deformation quantity and the deformation rate;
and S502, determining mine area closing time for closing the mine area based on the decline starting time.
The mining area closing time for closing the mining area is determined based on the decline starting time, so that an obtaining method can be provided for the mining area closing time, and the mining area closing time can be conveniently obtained.
Further, in other application scenarios, it is necessary to evaluate whether the closing time of the mine area is reasonable, and since the closing time of the mine area should be approximately the same as the decay start time of the decay phase, in some embodiments of the present invention, as shown in fig. 6, the method for monitoring and analyzing the deformation in the mine area closing phase further includes:
s601, obtaining mine area closing time for closing the mine area;
s602, determining the regression starting time of closing the surface movement process of the mining area to be a regression stage based on the surface deformation quantity and the deformation rate;
and S603, determining whether the closing time of the mining area is reasonable or not based on the decline starting time.
It should be understood that: step S603 specifically includes: and judging whether the difference between the decline start time and the mining area closing time is greater than a threshold difference, if the difference between the decline start time and the mining area closing time is greater than the threshold difference, the mining area closing time is unreasonable, and if the difference between the decline start time and the mining area closing time is less than or equal to the threshold difference, the mining area closing time is reasonable.
In order to improve the refinement of the evaluation of the surface damage of the closed mine, in some embodiments of the invention, as shown in fig. 7, the closed mine stage deformation monitoring and analyzing method further includes:
s701, obtaining various types of surface ground objects for closing the mining area;
s702, constructing corresponding relations between potential damage levels of multiple types of surface features and all stages of the surface movement process;
s703, determining a target stage of the earth surface moving process based on the earth surface deformation quantity and the deformation rate;
and S704, determining potential damage levels of the multi-class ground features based on the target stages and the corresponding relations.
In the embodiment of the present invention, if the mining area is closed and includes three types of surface features, i.e. roads, buildings and arable land, the corresponding relationship of step S702 is as shown in table 1:
TABLE 1 correspondence of potential damage levels to various stages of the earth's surface movement process
| Initial stage | Active phase | Decline stage | |
| Road | Ⅰ | Ⅱ | Ⅰ |
| Building construction | Ⅰ | Ⅱ | Ⅰ |
| Cultivation of land | Ⅱ | Ⅲ | Ⅰ |
Through the corresponding relation, the potential damage levels of various earth surface ground objects in all stages of the current earth surface moving process can be obtained, the potential damage levels can be referred by workers, fine evaluation on the potential damage of the closed mining area is realized, and the safety and the reliability of the closed mining area are improved.
In order to avoid the technical problem of inaccurate evaluation of the reliability of the closed mining area when the surface deformation amount and the deformation rate determined based on the small baseline set interferometry are inaccurate, in some embodiments of the present invention, as shown in fig. 8, after step S102, the method further includes:
s801, obtaining a reference surface deformation quantity for closing a mining area based on a Global Navigation Satellite System (GNSS) leveling method;
s802, converting the surface deformation into a vertical surface deformation;
and S803, determining whether the precision of the surface deformation quantity meets the precision requirement according to the reference surface deformation quantity and the vertical surface deformation quantity.
It should be understood that: step S103 is performed when the accuracy of the surface deformation quantity satisfies the accuracy requirement.
According to the embodiment of the invention, the accuracy of the surface deformation is verified through the reference surface deformation obtained based on the GNSS leveling method, so that the reliability of the surface deformation can be improved, and the accuracy of reliability evaluation on closed mining areas is further improved.
In an embodiment of the present invention, the conversion formula in step S802 is:
H＝D/cosθ
in the formula, H is a vertical surface deformation amount; d is the deformation of the earth surface; cos θ is the radar incidence angle of the SAR image.
It should be understood that: the surface deformation is the visual line deformation.
Those skilled in the art will appreciate that all or part of the flow of the method implementing the above embodiments may be implemented by instructing relevant hardware (such as a processor, a controller, etc.) by a computer program, and the computer program may be stored in a computer readable storage medium. The computer readable storage medium is a magnetic disk, an optical disk, a read-only memory or a random access memory.
The method for monitoring and analyzing the closed mine section staged deformation provided by the invention is described in detail, a specific example is applied in the method for explaining the principle and the implementation mode of the invention, and the description of the embodiment is only used for helping to understand the method and the core idea of the invention; meanwhile, for those skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed, and in summary, the content of the present specification should not be construed as limiting the present invention.
Claims (9)
1. A closed mining area stage deformation monitoring and analyzing method is characterized by comprising the following steps:
acquiring a plurality of synthetic aperture radar images of a closed mining area within a preset time period;
performing interference processing on the plurality of synthetic aperture radar images based on a small baseline set interferometry method, and determining the surface deformation amount and deformation rate of the closed mining area in the preset time period;
determining a surface movement phase of the closed mine area based on the surface deformation amount and the deformation rate, wherein the surface movement phase comprises a start phase, an active phase and a decline phase;
determining stability of the shut-down mine area based on the surface movement phase.
2. The closed mine phased deformation monitoring and analysis method of claim 1, wherein the determining the ground movement phase of the closed mine based on the ground deformation amount and the deformation rate comprises:
when the deformation amount of the earth surface of the closed mining area is greater than or equal to the initial deformation amount and the deformation rate is less than the threshold deformation rate, the earth surface movement process of the closed mining area is the starting stage;
when the deformation rate is larger than or equal to a threshold deformation rate, the process of closing the earth surface of the mining area is the active stage;
and when the deformation rate is smaller than a threshold deformation rate, or the surface deformation amount is smaller than or equal to a threshold deformation amount in a threshold time period, the surface movement process of the closed mining area is in the decline stage.
3. The closed mine staged deformation monitoring and analyzing method according to claim 1, wherein the interference processing of the plurality of synthetic aperture radar images based on the small baseline set interferometry method to determine the amount and rate of deformation of the ground surface of the closed mine within the preset time period comprises:
dividing the plurality of synthetic aperture radar images into a plurality of small baseline sets based on a spatial baseline threshold and a temporal baseline threshold, wherein the small baseline sets comprise a plurality of synthetic aperture radar image pairs;
determining a deformation time sequence of each synthetic aperture radar image pair in the plurality of synthetic aperture radar image pairs based on a least square method;
and obtaining the surface deformation quantity and the deformation rate of the closed mining area in the preset time period based on a singular value decomposition method and the deformation time sequence.
4. The closed mine staged deformation monitoring and analysis method according to claim 1, wherein the determining the stability of the closed mine according to the surface movement stage comprises:
judging whether the ground surface movement stage of the closed mining area is the decline stage;
if the earth surface movement stage of the closed mining area is the decline stage, the closed mining area is stable; and if the ground surface movement stage of the closed mining area is the starting stage or the active stage, the closed mining area is unstable.
5. The closed-mine staged deformation monitoring and analysis method according to claim 1, further comprising:
determining a decline start time of the decline stage as a course of surface movement of the shut-down mine area based on the amount of surface deformation and the rate of deformation;
determining a mine area closure time for the closed mine area based on the decline start time.
6. The closed-mine staged deformation monitoring and analysis method according to claim 1, further comprising:
obtaining the closing time of the mine area;
determining a decline start time of the decline stage as a course of surface movement of the shut-down mine area based on the amount of surface deformation and the rate of deformation;
determining whether the mine shut down time is reasonable based on the decline start time.
7. The closed mine phased deformation monitoring and analysis method according to claim 1, further comprising:
obtaining various types of surface ground objects of the closed mining area;
constructing corresponding relations between potential damage levels of the various types of surface ground objects and stages of the surface moving process;
determining a target stage of the earth surface movement process based on the earth surface deformation quantity and the deformation rate;
and determining the potential damage level of the multiple types of surface ground objects based on the target stage and the corresponding relation.
8. The closed mine staged deformation monitoring and analyzing method according to claim 1, wherein after the interferometric processing of the plurality of synthetic aperture radar images based on the small baseline set interferometry method determines the amount and rate of deformation of the ground surface of the closed mine within the preset time period, further comprising:
obtaining a reference surface deformation quantity of the closed mining area based on a global satellite navigation system leveling method;
converting the surface deformation amount into a vertical surface deformation amount;
and determining whether the precision of the surface deformation quantity meets the precision requirement or not according to the reference surface deformation quantity and the vertical surface deformation quantity.
9. The closed mine phased deformation monitoring and analysis method of claim 1, wherein the predetermined time period covers a mine closing time of the mine area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210337255.1A CN114777633B (en) | 2022-04-02 | Method for monitoring and analyzing staged deformation of closed mining area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210337255.1A CN114777633B (en) | 2022-04-02 | Method for monitoring and analyzing staged deformation of closed mining area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114777633A true CN114777633A (en) | 2022-07-22 |
| CN114777633B CN114777633B (en) | 2024-05-14 |
Family
ID=
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103728604A (en) * | 2013-11-19 | 2014-04-16 | 中国国土资源航空物探遥感中心 | Broadband synthetic aperture radar sub-band interferometric data processing method |
| CN105444730A (en) * | 2015-11-12 | 2016-03-30 | 中国矿业大学 | Time-space characteristic and cross-border mining identification method for multi-source data monitoring mining area deformation |
| CN106066478A (en) * | 2016-05-27 | 2016-11-02 | 中国矿业大学 | Merge pixel skew to follow the tracks of and the mining area surface deformation calculation method of short baseline set |
| CN106226767A (en) * | 2016-07-12 | 2016-12-14 | 中南大学 | Mining area three-D sequential deformation monitoring method based on single radar imagery geometry SAR image |
| CN108010280A (en) * | 2017-12-26 | 2018-05-08 | 成都理工大学 | A kind of sudden Loess Landslide method for early warning and its application |
| CN108416150A (en) * | 2018-03-15 | 2018-08-17 | 河北工业大学 | A kind of bevel edge slope early alarming and forecasting method with rheology and acceleration deformation behaviour |
| CN109738892A (en) * | 2019-01-24 | 2019-05-10 | 中南大学 | A kind of mining area surface high-spatial and temporal resolution three-dimensional deformation estimation method |
| CN110888130A (en) * | 2019-10-30 | 2020-03-17 | 华东师范大学 | Coal mine area ground surface deformation monitoring method based on lifting rail time sequence InSAR |
| CN111650579A (en) * | 2020-06-12 | 2020-09-11 | 中南大学 | InSAR mining area three-dimensional deformation estimation method and device for rock migration parameter adaptive acquisition and medium |
| CN112359807A (en) * | 2020-09-27 | 2021-02-12 | 中建东设岩土工程有限公司 | Device and method for monitoring influence of deep foundation pit excavation on vertical deformation of subway |
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103728604A (en) * | 2013-11-19 | 2014-04-16 | 中国国土资源航空物探遥感中心 | Broadband synthetic aperture radar sub-band interferometric data processing method |
| CN105444730A (en) * | 2015-11-12 | 2016-03-30 | 中国矿业大学 | Time-space characteristic and cross-border mining identification method for multi-source data monitoring mining area deformation |
| CN106066478A (en) * | 2016-05-27 | 2016-11-02 | 中国矿业大学 | Merge pixel skew to follow the tracks of and the mining area surface deformation calculation method of short baseline set |
| CN106226767A (en) * | 2016-07-12 | 2016-12-14 | 中南大学 | Mining area three-D sequential deformation monitoring method based on single radar imagery geometry SAR image |
| CN108010280A (en) * | 2017-12-26 | 2018-05-08 | 成都理工大学 | A kind of sudden Loess Landslide method for early warning and its application |
| CN108416150A (en) * | 2018-03-15 | 2018-08-17 | 河北工业大学 | A kind of bevel edge slope early alarming and forecasting method with rheology and acceleration deformation behaviour |
| CN109738892A (en) * | 2019-01-24 | 2019-05-10 | 中南大学 | A kind of mining area surface high-spatial and temporal resolution three-dimensional deformation estimation method |
| CN110888130A (en) * | 2019-10-30 | 2020-03-17 | 华东师范大学 | Coal mine area ground surface deformation monitoring method based on lifting rail time sequence InSAR |
| CN111650579A (en) * | 2020-06-12 | 2020-09-11 | 中南大学 | InSAR mining area three-dimensional deformation estimation method and device for rock migration parameter adaptive acquisition and medium |
| CN112359807A (en) * | 2020-09-27 | 2021-02-12 | 中建东设岩土工程有限公司 | Device and method for monitoring influence of deep foundation pit excavation on vertical deformation of subway |
Non-Patent Citations (1)
| Title |
|---|
| 刘金龙;郭华东;宋瑞;ZBIGNIEW PERSKI;刘广;赵峰;范景辉;闫世勇;阮智星;: "多模式雷达在矿区沉降监测中的应用研究", 遥感技术与应用, no. 04, 15 August 2012 (2012-08-15) * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Brêda et al. | Assimilation of satellite altimetry data for effective river bathymetry | |
| Agram et al. | A noise model for InSAR time series | |
| Setianto et al. | Comparison of kriging and inverse distance weighted (IDW) interpolation methods in lineament extraction and analysis | |
| Choi et al. | Geologic inheritance and earthquake rupture processes: The 1905 M≥ 8 Tsetserleg‐Bulnay strike‐slip earthquake sequence, Mongolia | |
| Fan et al. | A new model for three-dimensional deformation extraction with single-track InSAR based on mining subsidence characteristics | |
| Hill et al. | Combination of geodetic observations and models for glacial isostatic adjustment fields in Fennoscandia | |
| Yang et al. | Deriving time-series three-dimensional displacements of mining areas from a single-geometry InSAR dataset | |
| CN104459696A (en) | SAR interference baseline precise estimating method based on flat-earth phase | |
| CN104035138A (en) | Accurate and rapid calculation method for global and local ocean disturbing gravity | |
| CN109978279B (en) | Ocean surface temperature area prediction method | |
| Chen et al. | A generalization of inverse distance weighting method via kernel regression and its application to surface modeling | |
| Zhou et al. | Accuracy analysis of dam deformation monitoring and correction of refraction with robotic total station | |
| Xue et al. | A closed-form of Newton method for solving over-determined pseudo-distance equations | |
| Zhang et al. | Long-term ground multi-level deformation fusion and analysis based on a combination of deformation prior fusion model and OTD-InSAR for longwall mining activity | |
| Nowicki et al. | Assimilation of the satellite SST data in the 3D CEMBS model | |
| Zheng et al. | Joint probability integral method and TCPInSAR for monitoring mining time-series deformation | |
| Revel et al. | A framework for estimating global‐scale river discharge by assimilating satellite altimetry | |
| Wang et al. | Analysis and prediction of regional land subsidence with InSAR technology and machine learning algorithm | |
| Mason et al. | Interpolation of an intertidal digital elevation model from heighted shorelines: a case study in the western Wash | |
| Igarashi et al. | Maximum tsunami height prediction using pressure gauge data by a Gaussian process at Owase in the Kii Peninsula, Japan | |
| Nguyen et al. | Spatial interpolation of meteorologic variables in Vietnam using the Kriging method | |
| US20200083679A1 (en) | System and method for analyzing survivability of an infrastructure link | |
| Porkoláb et al. | Active tectonics of the Circum-Pannonian region in the light of updated GNSS network data | |
| CN114777633B (en) | Method for monitoring and analyzing staged deformation of closed mining area | |
| CN114777633A (en) | Monitoring and analyzing method for closed mining area stage deformation |
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 | ||
| GR01 | Patent grant |