Single-track InSAR three-dimensional deformation solution method and device for subsidence area

Abstract

The present invention relates to the field of deformation monitoring technology, and in particular to a single-track InSAR three-dimensional deformation solution method and device for a subsidence area, wherein the method comprises: calculating a preliminary subsidence field of a target subsidence area, and correcting the preliminary subsidence field based on the surface displacement field of the radar line of sight direction LOS of the single-track synthetic aperture radar interferometry InSAR to obtain a final subsidence field; calculating the horizontal displacement coefficient between the horizontal displacement of each position in the target subsidence area and the surface subsidence gradient based on the final subsidence field and the surface displacement field; calculating the horizontal displacement of the target subsidence area in the target direction using the horizontal displacement coefficient, and combining the final subsidence field with the horizontal displacement in the target direction to obtain a three-dimensional deformation field of the target subsidence area. The present invention can combine a priori models with observation data to obtain a three-dimensional deformation field of a target subsidence area, and can provide more accurate surface deformation data when dealing with complex deformation characteristics of goaf areas, thereby improving the monitoring capability and safety of goaf areas.

Description

Technical Field

[0001] The present invention relates to the technical field of deformation monitoring, and in particular to a method and a device for solving single-track InSAR three-dimensional deformation in a subsidence area. Background Art

[0002] Mining subsidence areas are widely distributed in my country. With the acceleration of social modernization, especially in mining and infrastructure construction, many projects have to be carried out in these areas. Therefore, three-dimensional deformation monitoring and evaluation of mining subsidence areas has become extremely necessary.

[0003] Among the related technologies, InSAR technology is widely used in surface deformation monitoring. It mainly monitors and evaluates the deformation of the collapse area by acquiring single-dimensional deformation data in LOS (line of sight direction).

[0004] However, the single-dimensional deformation data obtained by InSAR technology in related technologies has high limitations and cannot meet the needs of more comprehensive three-dimensional deformation research required when facing the increasingly complex deformation characteristics of goaf areas in actual engineering. How to convert InSAR LOS displacement data into three-dimensional deformation solution of the surface needs to be solved urgently. Summary of the invention

[0005] The present invention provides a single-track InSAR three-dimensional deformation solution method and device for a subsidence area, so as to solve the problems that the single-dimensional deformation data obtained by InSAR technology in the related technology has high limitations and cannot meet the more comprehensive three-dimensional deformation research needs required when facing the increasingly complex deformation characteristics of goaf areas in actual engineering, and how to convert InSAR LOS displacement data into three-dimensional deformation solution of the surface.

[0006] The first aspect of the present invention provides a single-track InSAR three-dimensional deformation solution method for a subsidence area, comprising the following steps: calculating a preliminary subsidence field of a target subsidence area, and correcting the preliminary subsidence field based on a surface displacement field in a radar line of sight direction LOS of a single-track synthetic aperture radar interferometry InSAR to obtain a final subsidence field; calculating a horizontal movement coefficient between a horizontal displacement of each position in the target subsidence area and a surface subsidence gradient based on the final subsidence field and the surface displacement field; calculating the horizontal displacement of a target direction of the target subsidence area using the horizontal movement coefficient, and obtaining a three-dimensional deformation field of the target subsidence area by combining the final subsidence field with the horizontal displacement in the target direction.

[0007] Optionally, in one embodiment of the present invention, before calculating the preliminary settlement field of the target subsidence area, it also includes: deriving the correspondence between the surface settlement gradient and the horizontal displacement of the target subsidence area according to preset prior information; based on the correspondence, deriving the relationship between the surface displacement field and the surface settlement gradient field in combination with SAR imaging features and geometric structure; and calculating the preliminary settlement field according to the relationship.

[0008] Optionally, in one embodiment of the present invention, the preliminary settlement field is corrected based on the surface displacement field in the radar line of sight direction LOS of single-track synthetic aperture radar interferometry InSAR to obtain a final settlement field, including: obtaining random errors and distribution characteristics of the preliminary settlement field; based on the surface displacement field, processing the random errors and distribution characteristics to correct the preliminary settlement field to obtain the final settlement field.

[0009] Optionally, in one embodiment of the present invention, the use of the horizontal movement coefficient to calculate the horizontal displacement of the target direction of the target subsidence area includes: when the actual working surface of the target subsidence area is the first preset surface, using the horizontal movement coefficient to calculate the horizontal displacement in the east-west direction and the horizontal displacement in the north-south direction of the target subsidence area to obtain the horizontal displacement in the target direction; or, obtaining the angle between the target long axis direction of the actual working surface and the due north direction, and calculating the horizontal displacement along the target long axis direction of the actual working surface and the horizontal displacement in the vertical direction according to the horizontal displacement in the west direction and the horizontal displacement in the north-south direction of the actual working surface and the angle to obtain the horizontal displacement in the target direction.

[0010] Optionally, in one embodiment of the present invention, the calculation formula of the horizontal movement coefficient between the horizontal displacement and the sedimentation gradient can be expressed as:

[0011]

[0012] in, is the horizontal shift coefficient, are the coordinate points on the working surface of the subsidence area, is the InSAR LOS displacement field to the surface, is the corrected settlement field of the target area, and are the resolutions of the spatial deformation map in the east-west and north-south directions, is the azimuth of the SAR orbit, is the angle of incidence.

[0013] The second aspect of the present invention provides an evaluation device for multi-form transient voltage stability of an electric power system, comprising: a first calculation module, used to calculate a preliminary settlement field of a target subsidence area, and to correct the preliminary settlement field based on a surface displacement field in a radar line of sight direction LOS of a single-track synthetic aperture radar interferometry InSAR to obtain a final settlement field; a second calculation module, used to calculate a horizontal movement coefficient between a horizontal displacement of each position in the target subsidence area and a surface settlement gradient based on the final settlement field and the surface displacement field; a solution module, used to calculate a horizontal displacement in a target direction of the target subsidence area using the horizontal movement coefficient, and to obtain a three-dimensional deformation field of the target subsidence area by combining the final settlement field and the horizontal displacement in the target direction.

[0014] Optionally, in one embodiment of the present invention, it also includes: a first derivation module, which is used to derive the correspondence between the surface settlement gradient and the horizontal displacement of the target subsidence area according to preset prior information before calculating the preliminary settlement field of the target subsidence area; a second derivation module, which is used to derive the relationship between the surface displacement field and the surface settlement gradient field based on the correspondence relationship in combination with SAR imaging features and geometric structure; and a third calculation module, which is used to calculate the preliminary settlement field according to the relationship.

[0015] Optionally, in one embodiment of the present invention, the first calculation module includes: an acquisition unit for acquiring random errors and distribution characteristics of the preliminary settlement field; a correction unit for processing the random errors and distribution characteristics based on the surface displacement field to correct the preliminary settlement field to obtain a final settlement field.

[0016] Optionally, in one embodiment of the present invention, the solution module includes: a first calculation unit, used to calculate the horizontal displacement in the east-west direction and the north-south direction of the target subsidence area using the horizontal movement coefficient when the actual working surface of the target subsidence area is the first preset surface, so as to obtain the horizontal displacement in the target direction; a second calculation unit, used to obtain the angle between the target long axis direction of the actual working surface and the north direction, and calculate the horizontal displacement along the target long axis direction of the actual working surface and the horizontal displacement in the vertical direction according to the horizontal displacement in the west direction and the horizontal displacement in the north-south direction of the actual working surface and the angle, so as to obtain the horizontal displacement in the target direction.

[0017] Optionally, in one embodiment of the present invention, the calculation formula of the horizontal movement coefficient between the horizontal displacement and the sedimentation gradient can be expressed as:

[0018]

[0019] in, is the horizontal shift coefficient, are the points on the working surface of the subsidence area, is the displacement field of InSARLOS toward the surface, is the corrected settlement field of the target area, and are the resolutions of the spatial deformation map in the east-west and north-south directions, is the azimuth of the SAR orbit, is the angle of incidence.

[0020] A third aspect of the present invention provides an electronic device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the single-track InSAR three-dimensional deformation solution method for a subsidence area as described in the above embodiment.

[0021] A fourth aspect of the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, which, when executed by a processor, implements the above-mentioned method for solving the three-dimensional deformation of a subsidence area with a single track InSAR.

[0022] A fifth aspect of the present invention provides a computer program product, including a computer program, which, when executed, is used to implement the above-mentioned method for solving the three-dimensional deformation of a subsidence area with a single track InSAR.

[0023] The embodiment of the present invention can combine the prior model and the InSAR LOS observation data to realize the single-track InSAR three-dimensional deformation solution of the surface, thereby calculating the three-dimensional deformation field of the target subsidence area. As a result, the limitations of the traditional method are overcome, and more accurate surface deformation data can be provided when dealing with the complex deformation characteristics of the goaf area, improving the monitoring ability and safety of the goaf subsidence area, and further ensuring the safe operation of the infrastructure in the target subsidence area, greatly improving the engineering safety, and also providing the necessary technical support for optimizing resource management and environmental protection. As a result, the single-dimensional deformation data obtained by InSAR technology in the related technology has high limitations and cannot meet the more comprehensive three-dimensional deformation research needs required in actual engineering when facing the increasingly complex deformation characteristics of the goaf, and how to convert the InSAR LOS displacement data into the three-dimensional deformation solution of the surface.

[0024] Additional aspects and advantages of the present invention will be given in part in the following description and in part will be obvious from the following description, or will be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] The above and / or additional aspects and advantages of the present invention will become apparent and easily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

[0026] Figure 1 A flowchart of a method for solving three-dimensional deformation of a subsidence area using a single-track InSAR according to an embodiment of the present invention;

[0027] Figure 2 A schematic diagram of calculation parameters of an embodiment of the present invention;

[0028] Figure 3 A flowchart of a single-track InSAR three-dimensional deformation solution method based on a priori model combined with observation data according to an embodiment of the present invention;

[0029] Figure 4 A preliminary settlement field according to an embodiment of the present invention Schematic diagram of calculation results;

[0030] Figure 5 The matched settlement field data of an embodiment of the present invention is Schematic diagram of the results;

[0031] Figure 6 The corrected settlement field data of an embodiment of the present invention is Schematic diagram of the results;

[0032] Figure 7 is a horizontal shift coefficient of an embodiment of the present invention Schematic diagram of calculation results;

[0033] Figure 8 The east-west horizontal movement of an embodiment of the present invention Schematic diagram of calculation results;

[0034] Fig. 9 The north-south horizontal movement of an embodiment of the present invention Schematic diagram of calculation results;

[0035] Fig.10 It is a schematic diagram of the structure of a single-track InSAR three-dimensional deformation solving device for a subsidence area according to an embodiment of the present invention;

[0036] Fig.11 It is a schematic diagram of the structure of an electronic device provided according to an embodiment of the present invention.

[0037] Reference numerals:

[0038] 10-Single-track InSAR three-dimensional deformation solution device for subsidence area: 100-first calculation module, 200-second calculation module, 300-solution module; 1101-memory, 1102-processor, 1103-communication interface. DETAILED DESCRIPTION

[0039] Embodiments of the present invention are described in detail below, examples of which are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present invention, and should not be construed as limiting the present invention.

[0040] The following describes the evaluation method and device of the multi-form transient voltage stability of the power system of the embodiment of the present invention with reference to the accompanying drawings. In view of the high limitation of the single-dimensional deformation data obtained by the InSAR technology in the related technology mentioned in the above background technology, it cannot meet the more comprehensive three-dimensional deformation research needs required when facing the increasingly complex deformation characteristics of the goaf area in actual engineering, and the problem of how to convert the InSAR LOS displacement data into the three-dimensional deformation solution of the surface, the present invention provides a single-track InSAR three-dimensional deformation solution method for the subsidence area, in which the single-track InSAR three-dimensional deformation solution of the surface can be realized based on the prior model and the observation data of InSAR LOS, so as to calculate the three-dimensional deformation field of the target subsidence area. In this way, the limitations of the traditional method are overcome, and more accurate surface deformation data can be provided when dealing with the complex deformation characteristics of the goaf area, so as to improve the monitoring capability and safety of the goaf subsidence area, and further ensure the safe operation of the infrastructure in the target subsidence area, greatly improving the engineering safety, and also providing the necessary technical support for optimizing resource management and environmental protection. This solves the problem that the single-dimensional deformation data obtained by InSAR technology in related technologies has high limitations and cannot meet the needs of more comprehensive three-dimensional deformation research required when facing increasingly complex deformation characteristics of goaf areas in actual engineering, and how to convert InSAR LOS displacement data into three-dimensional deformation solution of the surface.

[0041] Specifically, Figure 1 The present invention provides a flowchart of a method for solving three-dimensional deformation of a subsidence area using a single-track InSAR.

[0042] like Figure 1 As shown, the single-track InSAR three-dimensional deformation solution method for the subsidence area includes the following steps:

[0043] In step S101, a preliminary settlement field of the target subsidence area is calculated, and the preliminary settlement field is corrected based on the surface displacement field in the radar line of sight direction LOS of single-track synthetic aperture radar interferometry InSAR to obtain a final settlement field.

[0044] It can be understood that the target subsidence area here can be understood as the subsidence area that needs to be solved for the three-dimensional deformation field, wherein the subsidence area here refers to the mining subsidence area.

[0045] It can be understood by professionals in this technical field that InSAR (Interferometric Synthetic Aperture Radar), that is, synthetic aperture radar interferometry technology, is a radar technology used in surveying and remote sensing. InSAR can use synthetic aperture radar to perform coherent processing on the data of two complex-valued images (images with both amplitude and phase) observed in the same area to obtain surface elevation information. For example, it can obtain two coherent single-view complex images with a certain viewing angle difference in the same area through repeated observations by two SAR antennas or one antenna, and obtain surface elevation information from their interferometric phase information, thereby reconstructing a digital elevation model of the surface. Single-track InSAR is synthetic aperture radar interferometry in single-track mode, for example, single-track dual-antenna lateral mode or single-track dual-antenna longitudinal mode.

[0046] The LOS of a single track InSAR can be understood here as the LOS (Line Of Sight) measured by InSAR, that is, the deformation of the radar line of sight. For example, when the radar antenna transmits a microwave signal to the ground and receives its echo, the slight deformation of the surface target will cause the phase of the echo signal to change. By comparing the phase information in SAR images at different times or different viewing angles, the LOS deformation information of the surface target can be extracted.

[0047] In some embodiments, before calculating the final three-dimensional deformation field of the target subsidence area, the settlement field of the target subsidence area may be calculated first. The preliminary settlement field here refers to the settlement field of the target subsidence area directly calculated without adjustment or correction, wherein the settlement field may be understood to refer to the settlement phenomenon and range formed in the target area due to the sinking of the ground or soil caused by some reason such as natural factors or human activities.

[0048] In view of the fact that the preliminary settlement field calculated directly is insufficient in data precision and accuracy and has certain errors, the embodiment of the present invention can also use the surface displacement field in the radar line of sight direction LOS of the single-track synthetic aperture radar interferometry InSAR to make certain corrections to the preliminary settlement field, so as to obtain the corrected settlement field, i.e., the final settlement field, whose data precision and accuracy meet the requirements of calculating the three-dimensional deformation field of the target settlement area.

[0049] Next, the calculation process and correction process of the preliminary settlement field in the embodiment of the present invention are further explained.

[0050] Optionally, in one embodiment of the present invention, before calculating the preliminary settlement field of the target subsidence area, it also includes: deriving the correspondence between the surface settlement gradient and the horizontal displacement of the target subsidence area based on preset prior information; based on the correspondence, combining SAR imaging features and geometric structure to derive the relationship between the surface displacement field and the surface settlement gradient field; and calculating the preliminary settlement field based on the relationship.

[0051] It can be understood that the preset prior information here can be understood as some pre-established prior information used in calculating the preliminary settlement field, such as random medium theory and some normal distribution models.

[0052] In some embodiments, when calculating the preliminary settlement field of the target subsidence area, the present invention is mainly implemented based on certain prior information. According to certain prior information, the present invention can derive the corresponding relationship between the surface settlement gradient and the horizontal displacement of the target subsidence area, and then combine the SAR imaging characteristics and geometric structure to derive the relationship between the surface displacement field and the surface settlement gradient field in the InSAR LOS direction, and finally calculate the preliminary settlement field.

[0053] Taking random medium theory and Knothe normal distribution model as prior information, when calculating the corresponding relationship between surface settlement and horizontal displacement in the target subsidence area, the corresponding relationship between the unit horizontal displacement and the surface settlement gradient can be expressed as, but not limited to, the following:

[0054] ,

[0055] in, is the unit horizontal displacement, B is the unit horizontal shift coefficient, T is the horizontal gradient of unit settlement, represents the unit settlement, The unit settles in x Directional derivative.

[0056] Next, the embodiment of the present invention can deduce the relationship between the surface displacement field and the surface subsidence gradient field in the InSAR LOS direction based on the correspondence between the unit horizontal displacement and the surface subsidence gradient combined with the SAR imaging characteristics and geometric structure, which can be expressed as but not limited to:

[0057] ,

[0058] ,

[0059] ,

[0060] ,

[0061] Among them, Δ E and Δ N are the resolutions of the spatial deformation map in the east-west and north-south directions, α is the azimuth of the SAR orbit, θ is the angle of incidence, is the InSAR LOS displacement field to the surface, is the horizontal shift coefficient, is the calculated current settlement field.

[0062] Finally, the embodiment of the present invention may, but is not limited to, use a related calculation method of random medium theory PIM (random medium theory and probability integral method) to determine the initial settlement field , the calculation method can be expressed as but not limited to:

[0063] ,

[0064] ,

[0065]

[0066]

[0067]

[0068] in, For the equivalent mining height in the middle compaction area, The equivalent height of the boundary void area is is the subsidence coefficient of the middle compaction zone, is the sinking coefficient of the boundary void area, r is the influence radius, θ 0 is the mining impact propagation angle, is the coal seam inclination, and Figure 2 FIG. 1 is a schematic diagram of calculation parameters of an embodiment of the present invention, as shown in FIG. Figure 2 As shown, the parameters , , , , , , , , , , , , , The meaning of Figure 2 .

[0069] After obtaining the preliminary sedimentation field, the embodiment of the present invention can correct the preliminary sedimentation field to obtain a final sedimentation field that can be used. Next, the correction process of the preliminary sedimentation field in the embodiment of the present invention is further explained.

[0070] Optionally, in one embodiment of the present invention, the preliminary settlement field is corrected based on the surface displacement field in the radar line of sight direction LOS of single-track synthetic aperture radar interferometry InSAR to obtain the final settlement field, including: obtaining random errors and distribution characteristics of the preliminary settlement field; based on the surface displacement field, processing random errors and distribution characteristics to correct the preliminary settlement field to obtain the final settlement field.

[0071] It can be understood based on the relevant descriptions of other embodiments that the present invention can, but is not limited to, use the surface displacement field in the InSAR LOS direction to correct the preliminary settlement field to obtain the final settlement field.

[0072] During the actual implementation process, the error between the calculated preliminary settlement field and the actual situation is relatively large. Therefore, when correcting the preliminary settlement field, the present invention can first obtain the random errors and distribution characteristics of the preliminary settlement field, and then use the surface displacement field to process these random errors to achieve the correction of the preliminary settlement field and obtain the final settlement field.

[0073] For example, the correction process of the preliminary settlement field can be expressed as follows, but is not limited to:

[0074] The present invention can firstly The data is thinned to a certain extent, and The data resolution remains consistent and matches The coordinate system of the data. Among them, thinning is also called downsampling or data simplification, which is the process of minimizing the number of data points while ensuring that the shape or main features of the data remain unchanged. The number of data after thinning will be greatly reduced, but it can basically reflect the basic shape characteristics of the original data or curve, saving space and time for subsequent processing.

[0075] Next, the relationship between the displacement in the LOS direction and the surface subsidence gradient can be simplified according to the SAR imaging geometry. The formula can be, but is not limited to, expressed as follows:

[0076] ,

[0077] in, is the surface displacement field in the InSAR LOS direction, is the calculated current settlement field, k is the proportional factor obtained by linking the measured data, is the incident angle of SAR, is the error term.

[0078] Furthermore, the embodiment of the present invention can calculate the objective function according to the calculated current settlement field. The objective function formula can be, but is not limited to, expressed as follows:

[0079] ,

[0080] in, Represents the error between the observed value and the predicted value.

[0081] Finally, the embodiment of the present invention can calculate the gradient according to the objective function and use the gradient descent method to update the settlement field to obtain the final settlement field. The specific calculation method can be but is not limited to the following:

[0082] ,

[0083]

[0084] in, This is the settlement site before the update. For the updated settlement site, The learning rate that controls the update step size.

[0085] It should be noted that the calculation process of the settlement field involves multiple updates (iterations), so the calculation process When the update (iteration) stops, it can represent the final settlement field in the embodiment of the present invention. The distinction is made here only for the convenience of explanation, but it is not a specific limitation of the present invention.

[0086] Step S102, based on the final settlement field and the surface displacement field, calculate the horizontal displacement coefficient between the horizontal displacement of each position in the target subsidence area and the surface settlement gradient. The calculation formula of the horizontal displacement coefficient between the horizontal displacement and the settlement gradient can be, but is not limited to, expressed as:

[0087]

[0088] in, is the horizontal shift coefficient, are the points on the working surface of the subsidence area, is the InSAR LOS displacement field to the surface, is the corrected settlement field of the target area, i.e. the final settlement field. and is the resolution of the spatial deformation map in the east-west and north-south directions, is the azimuth of the SAR orbit, is the angle of incidence.

[0089] In other embodiments, after obtaining the final settlement field, the embodiments of the present invention can calculate the horizontal displacement coefficient between the horizontal displacement of each point in the target settlement area and the surface settlement gradient based on the final settlement field combined with the surface displacement field in the InSAR LOS direction, so as to calculate the three-dimensional deformation field of the target settlement area using the horizontal movement coefficient.

[0090] For example, the present invention can be carried out by modifying the settling field and InSAR LOS displacement field , calculate the horizontal displacement coefficient between each point and the sedimentation gradient , the calculation formula of the horizontal shift coefficient can be expressed as but not limited to:

[0091]

[0092] in, is the horizontal shift coefficient, are the points on the working surface of the subsidence area, is the InSAR LOS displacement field to the surface, is the corrected settlement field of the target area, i.e. the final settlement field. and is the resolution of the spatial deformation map in the east-west and north-south directions, is the azimuth of the SAR orbit, is the angle of incidence.

[0093] In addition, the embodiments of the present invention can also combine certain experience constraints on this basis to regulate The calculation results can be improved in accuracy and rationality. The specific empirical constraints can be selected or adjusted by professional and technical personnel in this field according to actual conditions and are not specifically limited here.

[0094] Step S103, using the horizontal movement coefficient to calculate the horizontal displacement of the target subsidence area in the target direction, combining the final settlement field with the horizontal displacement in the target direction to obtain a three-dimensional deformation field of the target subsidence area.

[0095] As a possible implementation method, after calculating the final settlement field and the horizontal movement coefficient, the embodiment of the present invention can calculate the horizontal displacement of the target and the settlement area in the target direction through the horizontal movement coefficient and obtain the three-dimensional deformation field of the target settlement area in combination with the final settlement field.

[0096] The three-dimensional deformation field can be understood as a field that can describe the deformation of the target subsidence area in three-dimensional space. For example, the three-dimensional deformation field can usually be represented by a series of mathematical models, such as deformation matrix, deformation vector, etc. In the deformation field, each point of the object has a corresponding deformation amount, which can be displacement, rotation, scaling, etc.

[0097] Next, the process of how to use the horizontal movement coefficient to calculate the horizontal displacement of the target subsidence area in the target direction in the embodiment of the present invention is further explained.

[0098] Optionally, in one embodiment of the present invention, the horizontal displacement of the target subsidence area in the target direction is calculated using the horizontal movement coefficient, including: when the actual working surface of the target subsidence area is the first preset surface, the horizontal displacement of the target subsidence area in the east-west direction and the horizontal displacement of the north-south direction are calculated using the horizontal movement coefficient to obtain the horizontal displacement in the target direction; or, the angle between the target long axis direction of the actual working surface and the due north direction is obtained, and the horizontal displacement along the target long axis direction of the actual working surface and the horizontal displacement in the vertical direction are calculated according to the horizontal displacement in the west direction and the horizontal displacement in the north-south direction of the actual working surface as well as the angle to obtain the horizontal displacement in the target direction.

[0099] It can be understood that the first preset surface here refers to a surface pre-established when calculating the horizontal displacement of a target subsidence area in a certain direction, for example, a working surface with a very standard actual distribution in the directions of southeast, northwest, etc.

[0100] During the actual implementation process, when the working surface of the target subsidence area is the first preset surface, that is, the working surface facing the southeast, northwest, and northeast, the present invention can use the horizontal movement coefficient to calculate the horizontal displacement of the target subsidence area in the east-west and north-south directions respectively, and use it as the horizontal displacement in the target direction to determine the three-dimensional deformation field of the target subsidence area.

[0101] Specifically, the embodiment of the present invention can use the calculated Further calculation of horizontal displacement in the east-west and north-south directions and , the calculation formula can be expressed as but not limited to:

[0102]

[0103] in, and Points East-west and north-south horizontal displacements, and Points The east-west and north-south horizontal gradients of subsidence, For point The settlement, For point The horizontal shift coefficient, and are the resolutions of the spatial deformation map in the east-west and north-south directions, respectively.

[0104] Furthermore, considering that the actual working surface of the target subsidence area is not distributed in the east-west-north-south direction, the embodiment of the present invention can assume that the angle between the long axis direction and the north direction is and Calculate the horizontal displacement along the long axis of the working surface and vertical displacement , the specific relationship can be expressed as but not limited to:

[0105]

[0106] in, It is the angle between the major axis and due north.

[0107] The present invention is described in detail with reference to a specific embodiment below.

[0108] Figure 3 FIG. 1 is a flow chart of a single-track InSAR three-dimensional deformation solution method based on a priori model combined with observation data according to an embodiment of the present invention. Figure 3 As shown:

[0109] (1) Data acquisition and study area description:

[0110] Take a coal mine as an example: suppose the coal mine contains two mined working faces, 2205 and 2207, and the monitoring period is from November 202x to May 202x. The study area uses Sentinel-1 satellite images, and obtains image data every 12 days to monitor the surface deformation in the study area. During the monitoring period, both working faces have completed the mining work. At the same time, GNSS measurement points are set on the surface above the 2205 working face to obtain surface displacement data, and GNSS data are used to adjust and verify the model parameters.

[0111] (2) Prior model construction and computational data acquisition:

[0112] The PIM (priori model) method was used to model the study area. Based on the surface displacement information of the GNSS measurement points, a settlement model covering the study area and a certain range of its periphery was obtained through inversion calculation. The model resolution was 5 meters, which was consistent with the spatial resolution of the Sentinel-1 InSAR data, so as to facilitate subsequent data matching and fusion. Figure 4 A preliminary settlement field according to an embodiment of the present invention Schematic diagram of calculation results, preliminary settlement field data calculated in the embodiment of the present invention That is Figure 4 As shown, this preliminary settlement field can serve as the basis for subsequent corrections and analyses.

[0113] (3) Matching of prior model data with InSAR data:

[0114] In order to ensure the spatial correspondence between the prior model data and the InSAR monitoring data, the embodiment of the present invention can, but is not limited to, use the nearest neighbor matching method to The data is matched to the InSAR data coordinate system. Taking each sampling point of the InSAR monitoring data as the benchmark, find Match the points at the same position in the data to obtain the matched settlement field data .like Figure 5 As shown, Figure 5 That is, the matched settlement field data of an embodiment of the present invention Schematic diagram of the results.

[0115] (4) Settlement field correction:

[0116] The matched prior model data is matched by the least squares method After correction, the corrected settlement field is obtained. The correction process is based on InSAR observation data. And the matching settlement area , using the following objective function to minimize the error, the correction algorithm can be expressed as follows but is not limited to:

[0117] ，

[0118] ，

[0119] ，

[0120] in, represents the error between the observed value and the predicted value, This is the settlement site before the update. For the updated settlement site, To control the learning rate parameter of the update step size, It can be determined according to the observation angle of the satellite. It is obtained based on GNSS measured data.

[0121] like Figure 6 As shown, Figure 6 That is, the corrected settlement field data of an embodiment of the present invention Schematic diagram of the results.

[0122] (5) Coefficient between horizontal displacement and sedimentation gradient Calculation:

[0123] After the correction of the settlement field Then, the proportional factor is calculated using the InSAR observation data and the corrected sedimentation field data. , which is used for the subsequent calculation of horizontal displacement. The calculation formula can be, but is not limited to, as follows:

[0124]

[0125] in, is the horizontal shift coefficient, are the points on the working surface of the subsidence area, is the displacement field of InSARLOS toward the surface, is the corrected settlement field of the target area, and are the resolutions of the spatial deformation map in the east-west and north-south directions, is the azimuth of the SAR orbit, is the angle of incidence.

[0126] The outermost circle is little affected by the goaf. Assign a value of 0 and combine it with the empirical constraint specification Calculation results, such as Figure 7 As shown, Figure 7 That is, the horizontal shift coefficient of one embodiment of the present invention Schematic diagram of calculation results.

[0127] (6) Calculation of east-west and north-south horizontal displacements:

[0128] Utilization coefficient and the corrected settlement field , calculate the horizontal displacement of the surface in the east and north directions, the formula can be but is not limited to the following:

[0129]

[0130] in, and Points East-west and north-south horizontal displacements, and Points The east-west and north-south horizontal gradients of subsidence, For point The settlement, For point The horizontal shift coefficient, ΔE and ΔN It is the resolution of the spatial deformation map in the east-west and north-south directions.

[0131] (4) For the outermost data points, the impact of the goaf is small. and The value is 0. The calculation result is as follows Figure 8 and Fig. 9 As shown, Figure 8 The east-west horizontal movement of an embodiment of the present invention Schematic diagram of calculation results. Fig. 9 The north-south horizontal movement of an embodiment of the present invention Schematic diagram of calculation results.

[0132] According to the single-track InSAR three-dimensional deformation solution method for the subsidence area proposed in the embodiment of the present invention, the single-track InSAR three-dimensional deformation solution of the surface can be realized based on the prior model and the InSAR LOS observation data, and the two can be combined to realize the three-dimensional deformation solution of the surface, so as to calculate the three-dimensional deformation field of the target subsidence area. Therefore, the limitations of the traditional method are overcome, and more accurate surface deformation data can be provided when dealing with the complex deformation characteristics of the goaf area, and the monitoring ability and safety of the goaf subsidence area can be improved, and then the safe operation of the infrastructure in the target subsidence area can be guaranteed, which greatly improves the engineering safety, and also provides the necessary technical support for optimizing resource management and environmental protection. Therefore, the single-dimensional deformation data obtained by InSAR technology in the related technology has high limitations and cannot meet the more comprehensive three-dimensional deformation research needs required in actual engineering when facing the increasingly complex deformation characteristics of the goaf area, and how to convert the InSAR LOS displacement data into the three-dimensional deformation solution of the surface.

[0133] Next, a single-track InSAR three-dimensional deformation solving device for a subsidence area according to an embodiment of the present invention will be described with reference to the accompanying drawings.

[0134] Fig.10 It is a structural schematic diagram of a subsidence area single-track InSAR three-dimensional deformation solution device according to an embodiment of the present invention.

[0135] like Fig.10 As shown, the subsidence area single-track InSAR three-dimensional deformation solving device 10 includes: a first calculation module 100, a second calculation module 200 and a solving module 300.

[0136] The first calculation module 100 is used to calculate the preliminary settlement field of the target settlement area, and to correct the preliminary settlement field based on the surface displacement field in the radar line of sight direction LOS of the single-track synthetic aperture radar interferometry InSAR to obtain the final settlement field.

[0137] The second calculation module 200 is used to calculate the horizontal displacement coefficient between the horizontal displacement of each point in the target subsidence area and the settlement gradient based on the final settlement field and the surface displacement field.

[0138] The solution module 300 is used to calculate the horizontal displacement of the target subsidence area in the target direction by using the horizontal movement coefficient, and obtain the three-dimensional deformation field of the target subsidence area by combining the final settlement field and the horizontal displacement in the target direction.

[0139] Optionally, in one embodiment of the present invention, it further includes: a first derivation module, a second derivation module and a third calculation module.

[0140] Among them, the first derivation module is used to derive the corresponding relationship between the surface settlement gradient and the horizontal displacement of the target subsidence area according to preset prior information before calculating the preliminary settlement field of the target subsidence area.

[0141] The second derivation module is used to derive the relationship between the surface displacement field and the surface settlement gradient field based on the corresponding relationship and in combination with the SAR imaging features and the geometric structure.

[0142] The third calculation module is used to calculate the preliminary settlement field according to the relationship.

[0143] Optionally, in one embodiment of the present invention, the first calculation module 100 includes: an acquisition unit and a correction unit.

[0144] Wherein, the acquisition unit is used to obtain the random error and distribution characteristics of the preliminary settlement field;

[0145] The correction unit is used to process random errors and their distribution characteristics based on the surface displacement field to correct the preliminary settlement field and obtain the final settlement field.

[0146] Optionally, in one embodiment of the present invention, the solving module 300 includes: a first calculating unit and a second calculating unit.

[0147] The first calculation unit is used to calculate the horizontal displacement in the east-west direction and the horizontal displacement in the north-south direction of the target subsidence area by using the horizontal movement coefficient when the actual working surface of the target subsidence area is the first preset surface, so as to obtain the horizontal displacement in the target direction;

[0148] The second calculation unit is used to obtain the angle between the target long axis direction of the actual working surface and the due north direction, and calculate the horizontal displacement along the target long axis direction of the actual working surface and the horizontal displacement in the vertical direction according to the horizontal displacement in the west direction and the horizontal displacement in the north-south direction of the actual working surface and the angle to obtain the horizontal displacement in the target direction.

[0149] Optionally, in one embodiment of the present invention, the calculation formula of the horizontal movement coefficient between the horizontal displacement and the sedimentation gradient can be expressed as but not limited to:

[0150]

[0151] in, is the horizontal shift coefficient, are the points on the working surface of the subsidence area, is the InSAR LOS displacement field to the surface, is the corrected settlement field of the target area, and is the resolution of the spatial deformation map in the east-west and north-south directions, is the azimuth of the SAR orbit, is the angle of incidence.

[0152] It should be noted that the above explanation of the embodiment of the method for solving the three-dimensional deformation of a subsidence area with a single track InSAR is also applicable to the device for solving the three-dimensional deformation of a subsidence area with a single track InSAR in this embodiment, and will not be repeated here.

[0153] According to the single-track InSAR three-dimensional deformation solution device for the subsidence area proposed in the embodiment of the present invention, the single-track InSAR three-dimensional deformation solution of the surface can be realized based on the prior model and the observation data of InSAR LOS, and the two can be combined to realize the three-dimensional deformation solution of the surface, so as to calculate the three-dimensional deformation field of the target subsidence area. Therefore, the limitations of the traditional method are overcome, and more accurate surface deformation data can be provided when dealing with the complex deformation characteristics of the goaf area, and the monitoring ability and safety of the goaf subsidence area can be improved, and then the safe operation of the infrastructure in the target subsidence area can be guaranteed, which greatly improves the engineering safety, and also provides the necessary technical support for optimizing resource management and environmental protection. Therefore, the single-dimensional deformation data obtained by InSAR technology in the related technology has high limitations and cannot meet the more comprehensive three-dimensional deformation research needs required in actual engineering when facing the increasingly complex deformation characteristics of the goaf area, and how to convert the InSAR LOS displacement data into the three-dimensional deformation solution of the surface.

[0154] Fig.11 A schematic diagram of the structure of an electronic device provided by an embodiment of the present invention. The electronic device may include:

[0155] A memory 1101 , a processor 1102 , and a computer program stored in the memory 1101 and executable on the processor 1102 .

[0156] When the processor 1102 executes the program, the subsidence area single-track InSAR three-dimensional deformation solution method provided in the above embodiment is implemented.

[0157] Furthermore, the electronic device further comprises:

[0158] The communication interface 1103 is used for communication between the memory 1101 and the processor 1102 .

[0159] The memory 1101 is used to store computer programs that can be executed on the processor 1102 .

[0160] The memory 1101 may include a high-speed RAM memory, and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.

[0161] If the memory 1101, the processor 1102 and the communication interface 1103 are implemented independently, the communication interface 1103, the memory 1101 and the processor 1102 can be connected to each other through a bus and communicate with each other. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Fig.11 Only one thick line is used in the diagram, but this does not mean that there is only one bus or only one type of bus.

[0162] Optionally, in a specific implementation, if the memory 1101, the processor 1102 and the communication interface 1103 are integrated on a chip, the memory 1101, the processor 1102 and the communication interface 1103 can communicate with each other through an internal interface.

[0163] The processor 1102 may be a central processing unit (CPU), or an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present invention.

[0164] An embodiment of the present invention further provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-mentioned method for solving the three-dimensional deformation of a subsidence area using a single track InSAR.

[0165] An embodiment of the present invention further provides a computer program product, including a computer program, which can run computer instructions. When the computer instructions are executed by a processor, the method for solving the single-track InSAR three-dimensional deformation of a subsidence area provided by an embodiment of the present invention is implemented.

[0166] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or N embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradiction.

[0167] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined as "first" and "second" may explicitly or implicitly include at least one of the features. In the description of the present invention, the meaning of "N" is at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0168] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, fragment or portion of code comprising one or N executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present invention includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the embodiments of the present invention belong.

[0169] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute instructions), or in combination with these instruction execution systems, devices or apparatuses. For the purposes of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in combination with these instruction execution systems, devices or apparatuses. More specific examples (non-exhaustive list) of computer-readable media include the following: an electrical connection with one or N wirings (electronic device), a portable computer disk box (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically by optically scanning the paper or other medium and then editing, interpreting or processing in other suitable ways as necessary and then storing it in a computer memory.

[0170] It should be understood that the various parts of the present invention can be implemented by hardware, software, firmware or a combination thereof. In the above embodiment, the N steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. If implemented by hardware, as in another embodiment, it can be implemented by any one or a combination of multiple of the following technologies known in the art: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0171] A person skilled in the art may understand that all or part of the steps in the method for implementing the above-mentioned embodiment may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, which, when executed, includes one or a combination of the steps of the method embodiment.

[0172] In addition, each functional unit in each embodiment of the present invention may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

[0173] The storage medium mentioned above may be a read-only memory, a magnetic disk or an optical disk, etc. Although the embodiments of the present invention have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limiting the present invention. A person of ordinary skill in the art may change, modify, replace and modify the above embodiments within the scope of the present invention.

Claims

1. A method for calculating the three-dimensional deformation of a subsidence area using a single-track InSAR, characterized in that: The following steps are involved: Calculate the preliminary subsidence field of the target subsidence area, and correct the preliminary subsidence field based on the surface displacement field in the radar line of sight direction LOS of single-track synthetic aperture radar interferometry InSAR to obtain the final subsidence field; Based on the final settlement field and the surface displacement field, calculating the horizontal movement coefficient between the horizontal displacement of each position in the target settlement area and the surface settlement gradient; Calculating the horizontal displacement of the target subsidence area in the target direction by using the horizontal movement coefficient, and obtaining the three-dimensional deformation field of the target subsidence area by combining the final settlement field and the horizontal displacement in the target direction; The method of correcting the preliminary settlement field based on the surface displacement field in the radar line of sight direction LOS of the single-track synthetic aperture radar interferometry InSAR to obtain the final settlement field includes: obtaining the random error and distribution characteristics of the preliminary settlement field; based on the surface displacement field, processing the random error and distribution characteristics to correct the preliminary settlement field to obtain the final settlement field; The calculation formula of the horizontal movement coefficient between the horizontal displacement and the sedimentation gradient is: Wherein, B(i,j) is the horizontal movement coefficient, (i,j) is the coordinate point on the working surface of the subsidence area, LOS(i,j) is the InSAR LOS displacement field to the surface, and W A (i, j) is the corrected settlement field of the target area, ΔE and ΔN are the resolutions of the spatial deformation map in the east-west and north-south directions, respectively, α is the azimuth of the SAR orbit, and θ is the incident angle.

2. The method for calculating the three-dimensional deformation of a subsidence area using a single track InSAR according to claim 1, characterized in that: Before calculating the preliminary settlement field of the target subsidence area, it also includes: Derivation of the corresponding relationship between the surface subsidence gradient and the horizontal displacement according to preset prior information; Based on the corresponding relationship, the relationship between the surface displacement field and the surface subsidence gradient field is derived in combination with SAR imaging features and geometric structure; The preliminary settlement field is calculated according to the relationship.

3. The method for calculating the three-dimensional deformation of a subsidence area using a single track InSAR according to claim 1, characterized in that: The calculating the horizontal displacement of the target subsidence area in the target direction by using the horizontal movement coefficient includes: When the actual working surface of the target subsidence area is the first preset surface, the horizontal displacement of the target subsidence area in the east-west direction and the horizontal displacement in the north-south direction are calculated using the horizontal movement coefficient to obtain the horizontal displacement in the target direction; Alternatively, the angle between the target long axis direction of the actual working surface and the due north direction is obtained, and the horizontal displacement along the target long axis direction of the actual working surface and the horizontal displacement in the vertical direction are calculated according to the horizontal displacement in the west direction and the horizontal displacement in the north-south direction of the actual working surface and the angle to obtain the horizontal displacement in the target direction.

4. A single-track InSAR three-dimensional deformation solution device for subsidence areas, characterized in that: include: The first calculation module is used to calculate the preliminary settlement field of the target settlement area, and to correct the preliminary settlement field based on the surface displacement field in the radar line of sight direction LOS of the single-track synthetic aperture radar interferometry InSAR to obtain the final settlement field; A second calculation module is used to calculate the horizontal displacement coefficient between the horizontal displacement of each position in the target subsidence area and the surface subsidence gradient based on the final subsidence field and the surface displacement field; A solution module, used to calculate the horizontal displacement of the target subsidence area in the target direction by using the horizontal movement coefficient, and obtain the three-dimensional deformation field of the target subsidence area by combining the final settlement field and the horizontal displacement in the target direction; The first calculation module includes: an acquisition unit for acquiring the random error and distribution characteristics of the preliminary settlement field; a correction unit for processing the random error and distribution characteristics based on the surface displacement field to correct the preliminary settlement field to obtain the final settlement field; The calculation formula of the horizontal movement coefficient between the horizontal displacement and the sedimentation gradient is: Wherein, B(i,j) is the horizontal movement coefficient, (i,j) is the coordinate point on the working surface of the subsidence area, LOS(i,j) is the InSAR LOS displacement field to the surface, and W A (i, j) is the corrected settlement field of the target area, ΔE and ΔN are the resolutions of the spatial deformation map in the east-west and north-south directions, respectively, α is the azimuth of the SAR orbit, and θ is the incident angle.

5. The device for calculating the three-dimensional deformation of a subsidence area using a single track InSAR according to claim 4, characterized in that: Also includes: A first derivation module is used to derive the corresponding relationship between the surface settlement gradient and the horizontal displacement according to preset prior information before calculating the preliminary settlement field of the target subsidence area; A second derivation module is used to derive the relationship between the surface displacement field and the surface sedimentation gradient field based on the corresponding relationship in combination with SAR imaging features and geometric structure; The third calculation module is used to calculate the preliminary settlement field according to the relationship.

6. An electronic device, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method for solving the three-dimensional deformation of a subsidence area using a single track InSAR as claimed in any one of claims 1 to 3.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that: The computer program is executed by a processor to implement the method for solving the three-dimensional deformation of a subsidence area using a single-track InSAR method as described in any one of claims 1 to 3.

8. A computer program product, comprising a computer program, characterized in that When the computer program is executed, it is used to implement the single-track InSAR three-dimensional deformation solution method for a subsidence area as described in any one of claims 1 to 3.

Images