CN110657742A - Aquifer deformation signal separation method, aquifer deformation signal separation device, aquifer deformation signal separation equipment and readable storage medium - Google Patents
Aquifer deformation signal separation method, aquifer deformation signal separation device, aquifer deformation signal separation equipment and readable storage medium Download PDFInfo
- Publication number
- CN110657742A CN110657742A CN201910942684.XA CN201910942684A CN110657742A CN 110657742 A CN110657742 A CN 110657742A CN 201910942684 A CN201910942684 A CN 201910942684A CN 110657742 A CN110657742 A CN 110657742A
- Authority
- CN
- China
- Prior art keywords
- deformation
- aquifer
- time sequence
- sequence
- signal separation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the invention discloses a method, a device and equipment for separating aquifer deformation signals and a computer readable storage medium. The method comprises the steps of obtaining a ground surface time sequence deformation sequence of a region to be detected based on an SBAS-InSAR technology, and separating elastic deformation time sequence information and inelastic deformation time sequence information of the ground surface time sequence deformation sequence by adopting a continuous wavelet transformation separation method so as to realize deformation signal separation of an aquifer of the region to be detected. The method adopts the SBAS-InSAR method to obtain the time sequence signal, can reduce the geometric and time decorrelation effect, obtain a wider time sequence time-varying signal, and improve the space-time resolution of monitoring on a large regional scale so as to obtain the higher-precision long-time sequence earth surface settlement deformation information; by adopting the CWT separation method, the elastic and inelastic deformation information of the aquifer can be more accurately separated, the independent inversion of the elastic and inelastic water release coefficients of the aquifer is facilitated, and the establishment of a high-precision underground water-ground settlement prediction model is facilitated.
Description
Technical Field
The embodiment of the invention relates to the field of surface deformation monitoring, in particular to a method, a device and equipment for separating aquifer deformation signals and a computer-readable storage medium.
Background
With the rapid development of urban economy, a series of geological disasters such as serious ground subsidence, surface subsidence and the like caused by large-scale urban construction and excessive underground water extraction occur, and the strengthening of urban surface subsidence monitoring is an effective means for preventing the occurrence of surface subsidence disasters.
In the related art, a signal based on a GPS is generally adopted to monitor the ground surface subsidence state, a plurality of single-point signals are monitored to obtain some discrete data, and then the acquired signal data are analyzed by adopting a wavelet change method. The wavelet transformation method is that a function uses a certain wavelet basis function to do localized mathematical transformation in time and space, the time information of the original function under the wavelet basis can be obtained by the translation of the wavelet basis, and then the frequency information is obtained by scaling the scale of the wavelet basis.
The method not only consumes a great deal of manpower and financial resources, but also is difficult to acquire data, and in addition, the wavelet transformation method has poor effect in analyzing and processing non-stationary signals, can only be used for carrying out fuzzy inversion on the range of the general water release coefficient, and cannot clearly separate the elastic deformation and the non-elastic deformation information of the aquifer caused by elastic compression.
Disclosure of Invention
The invention aims to provide a method, a device, equipment and a computer readable storage medium for separating aquifer deformation signals, which can accurately separate elastic deformation and inelastic deformation information of an aquifer caused by elastic compression.
In order to solve the above technical problems, embodiments of the present invention provide the following technical solutions:
the embodiment of the invention provides a method for separating aquifer deformation signals, which comprises the following steps:
acquiring a ground surface time sequence deformation sequence of a region to be detected based on an SBAS-InSAR technology;
and separating the elastic deformation time sequence information and the inelastic deformation time sequence information of the surface time sequence deformation sequence by adopting a continuous wavelet change separation method so as to realize the separation of deformation signals of the aquifer of the region to be detected.
Optionally, after the elastic deformation time sequence information and the inelastic deformation time sequence information of the earth surface time sequence deformation sequence are separated by using a continuous wavelet transformation separation method, the method further includes:
and inverting the elastic water release coefficient and the inelastic water release coefficient of the aquifer by utilizing a pre-constructed ground subsidence soil layer compression model.
Optionally, the ground settlement soil layer compression model is as follows:
calculating the settlement of the aquifer using the following formula:
ΔbInSAR＝-miΔtA；
wherein A is the thickness of the aquifer,. DELTA.bInSARIs the amount of compression of the aquifer over a time Δ t, SskAnd h is the water release coefficient of the aquifer, wherein the water release coefficient comprises an elastic water release coefficient and an inelastic water release coefficient, and h is the water level value of the aquifer water head.
Optionally, the obtaining of the surface time sequence deformation sequence of the region to be detected based on the SBAS-InSAR technology includes:
acquiring all interferograms of a region to be detected acquired by the SBAS-InSAR technology;
unwrapping the wrapping phase of each interference pattern by using a three-dimensional unwrapping method to obtain phase information of each interference pattern;
and calculating corresponding average deformation speed according to the phase information of each interference pattern, and performing time domain high-pass filtering and space domain low-pass filtering to generate a surface time sequence deformation sequence of the region to be detected.
Optionally, after the obtaining of all interferograms of the area to be measured by the SBAS-InSAR technology, the method further includes:
carrying out self-adaptive filtering on each interference pattern in space, and selecting pixels meeting conditions as initial high coherence points;
calculating the standard deviation of the noise phase of the high coherence point of each interferogram, and discarding the high coherence points which do not meet the deviation condition;
the high coherence points in each interferogram are error corrected using the SRTM3 DEM.
Optionally, the separating the elastic deformation time sequence information and the inelastic deformation time sequence information of the surface time sequence deformation sequence by using a continuous wavelet transformation separation method includes:
separating the earth surface time sequence deformation sequence into a long period and a short period by using a continuous wavelet change separation method;
and calculating the intensity range of each period by utilizing the Fourier change, and calculating the intensity value of each period of the corresponding basic wave according to each intensity range and the basic wave fitting metadata.
Another aspect of an embodiment of the present invention provides an aquifer deformation signal separation apparatus, including:
the system comprises a surface time sequence deformation sequence acquisition module, a time sequence deformation sequence acquisition module and a time sequence deformation sequence acquisition module, wherein the surface time sequence deformation sequence acquisition module is used for acquiring a surface time sequence deformation sequence of a region to be detected based on an SBAS-InSAR technology;
and the signal separation module is used for separating the elastic deformation time sequence information and the inelastic deformation time sequence information of the surface time sequence deformation sequence by adopting a continuous wavelet change separation method so as to realize the deformation signal separation of the aquifer of the region to be detected.
Optionally, the method further includes:
and the coefficient inversion module is used for inverting the elastic water release coefficient and the inelastic water release coefficient of the aquifer by utilizing a pre-constructed ground subsidence soil layer compression model.
An embodiment of the present invention further provides an aquifer deformation signal separation device, which includes a processor, and the processor is configured to implement the steps of the aquifer deformation signal separation method according to any one of the preceding items when executing the computer program stored in the memory.
Finally, an embodiment of the present invention provides a computer readable storage medium, in which an aquifer deformation signal separation program is stored, and the aquifer deformation signal separation program, when executed by a processor, implements the steps of the aquifer deformation signal separation method according to any one of the previous items.
The technical scheme provided by the application has the advantages that the SBAS-InSAR method is adopted to obtain the time sequence signal of the area to be detected, so that the geometric and time decorrelation effect can be reduced, the time sequence time-varying signal of a wider area can be obtained, and the space-time resolution of monitoring on a large area scale is improved, so that the ground surface settlement deformation information of a high-precision long-time sequence can be obtained; the method has the advantages that the CWT method is adopted to separate the elastic time sequence deformation signals and the inelastic time sequence deformation signals of the aquifer, more accurate parameter basis is provided for the current situation of difficult separation caused by delay between the elastic time sequence deformation signals and the inelastic InSAR time sequence deformation signals of the aquifer, the problem of poor effect when the related technology processes the unstable signals is solved, the elastic deformation information and the inelastic deformation information of the aquifer can be separated more accurately, independent inversion of the elastic water release coefficients and the inelastic water release coefficients of the aquifer is facilitated, and then the high-precision underground water-ground settlement prediction model is facilitated to be established.
In addition, the embodiment of the invention also provides a corresponding implementation device, equipment and a computer readable storage medium for the aquifer deformation signal separation method, so that the method has higher practicability, and the device, the equipment and the computer readable storage medium have corresponding advantages.
Drawings
In order to more clearly illustrate the embodiments or technical solutions of the present invention, the drawings used in the description of the embodiments or the prior art will be briefly described 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 that other drawings can be obtained based on these drawings without creative efforts.
Fig. 1 is a schematic flow chart of a method for separating aquifer deformation signals according to an embodiment of the present invention;
FIG. 2 is a diagram of the change of ground subsidence of a monitoring point No. 1 of a test field based on wavelet analysis (sampling interval: 15 days), wherein a very coordinate represents a time interval and an ordinate table represents frequency;
fig. 3 is a ground settlement sequence diagram (sampling interval: 15 days) of a monitoring point No. 1 in an experimental site, where the abscissa represents time interval and the ordinate represents ground settlement;
FIG. 4 is a schematic flow chart of another method for separating aquifer deformation signals according to an embodiment of the invention;
FIG. 5 is a diagram illustrating a parameter inversion result according to an embodiment of the present invention;
fig. 6 is a block diagram of an embodiment of an aquifer deformation signal separation device according to an embodiment of the present invention;
fig. 7 is a block diagram of another embodiment of an aquifer deformation signal separation device according to an embodiment of the invention.
Detailed Description
In order that those skilled in the art will better understand the disclosure, the invention will be described in further detail with reference to the accompanying drawings and specific embodiments. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The terms "first," "second," "third," "fourth," and the like in the description and claims of this application and in the above-described drawings are used for distinguishing between different objects and not for describing a particular order. Furthermore, the terms "comprising" and "having," as well as any variations thereof, are intended to cover non-exclusive inclusions. For example, a process, method, system, article, or apparatus that comprises a list of steps or elements is not limited to only those steps or elements but may include other steps or elements not expressly listed.
Having described the technical solutions of the embodiments of the present invention, various non-limiting embodiments of the present application are described in detail below.
Referring to fig. 1, fig. 1 is a schematic flow chart of a method for separating an aquifer deformation signal according to an embodiment of the present invention, where the embodiment of the present invention may include the following:
s101: and acquiring the surface time sequence deformation sequence of the area to be detected based on the SBAS-InSAR technology.
S102: and separating the elastic deformation time sequence information and the inelastic deformation time sequence information of the earth surface time sequence deformation sequence by adopting a continuous wavelet change separation method so as to realize the separation of deformation signals of the aquifer of the region to be detected.
The SBAS (Satellite-Based Augmentation System) carries a Satellite navigation Augmentation signal transponder through a geostationary orbit Satellite, and broadcasts various correction information such as ephemeris error, Satellite clock error, ionospheric delay and the like to a user, thereby realizing improvement of positioning accuracy of the original Satellite navigation System.
The method comprises the steps that an InSAR (Synthetic Aperture Radar interference) transmits microwaves to a target area by utilizing a Radar, then echoes reflected by the target are received, an SAR complex image pair imaged in the same target area is obtained, if a coherence condition exists between the complex image pairs, an interference diagram can be obtained by conjugate multiplication of the SAR complex image pair, and the path difference of the microwaves in two times of imaging is obtained according to the phase value of the interference diagram, so that the terrain, the landform and the surface micro-change of the target area are calculated.
In S101, acquiring all interferograms of a region to be detected by an SBAS-InSAR technology; filtering all interferograms in space by using band-pass filtering as self-adaptive filtering, and selecting pixels meeting conditions (such as stable pixels) as initial high coherence points after filtering is finished; the noise phase standard deviation of the high coherence points of each interferogram is calculated and high coherence points that do not satisfy the deviation condition, e.g., high coherence points with a noise phase standard deviation greater than 1 radian, typically unstable points, which can be discarded.
Alternatively, the following formula can be used to extract the high coherence point:
in the formula, Ψx,iThe interference phase of the x pixel in the ith interference pattern, N is the number of the interference patterns,estimates of the spatially coherent part, such as surface subsidence information, atmospheric disturbances and orbital errors,for spatial incoherent partial view error, gammaxThe coherence coefficient is the time dimension of the high coherence point, theta is the satellite incidence angle, and u is the incoherent property.
After discarding some of the high coherence points from each interferogram, the high coherence points in each interferogram are error corrected using the SRTM3 DEM (a global 30 meter resolution digital elevation model).
And (3) unwrapping the wrapping phase of each interference pattern by using a three-dimensional unwrapping method to obtain phase information of each scene during image acquisition, estimating the average deformation speed of each interference pattern, performing time domain high-pass filtering and space domain low-pass filtering, separating deformation information from errors such as atmosphere and noise, and finally generating a surface time sequence deformation sequence of the region to be detected.
In S102, time series x (n) { x) of continuous wavelet transform separation method (CWT)n}n＝1Is the time sequence of InSAR deformation and water level change:
W(a,n)＝Ψ(a,n)X(n)；
where Ψ is a circulant matrix of n × n, Ψ is a first row vector of Ψ, a is a scaling function, each row vector is a derivative of a gaussian wavelet function that can be used with degree m through a cycle relative to a previous element, and depending on a dimensionless time parameter η, can be defined as follows:
by setting the value of m, a high spectral resolution suitable for identifying high frequency, low amplitude signals in a time series is provided.
Separating the earth surface time sequence deformation sequence into a long period and a short period by using a continuous wavelet change separation method; and calculating the intensity range of each period by utilizing the Fourier change, sorting the fitting metadata of the basic wave, and calculating the intensity value of each period of the corresponding basic wave according to each intensity range and the fitting metadata of the basic wave.
In order to confirm that the technical scheme provided by the application can realize accurate classification of aquifer elastic and inelastic time sequence deformation signals, please refer to fig. 2-3, and it can be seen from the graphs, the vertical InSAR and water level change time sequence pixel wavelet power spectrum analysis graph has maximum settlement or swelling in each feature, and color and contour represent the intensity of the amplitude or frequency component of the wavelet coefficient. The vertical axis represents the frequency (reciprocal of period), the horizontal axis represents time, and the color represents the intensity of the frequency. In order to make the intensity variation more apparent here is the logarithm of the real intensity. It can be seen that there are 2 high intensity regions at approximately 0.15Hz and 0.20 Hz. Therefore, the acquired InSAR time sequence deformation signal is separated into a long period and a short period by using a CWT continuous wavelet transform method, the long period signal is assumed to be an inelastic deformation signal caused by aquifer compression, and the short period signal is assumed to be an aquifer deformation signal caused by elastic compression, so that the independent inversion of elastic and inelastic water release coefficients is achieved.
In the technical scheme provided by the embodiment of the invention, the SBAS-InSAR method is adopted to obtain the time sequence signal of the area to be detected, so that the geometric and time decorrelation effect can be reduced, a wider time sequence time-varying signal can be obtained, and the space-time resolution of monitoring on a large area scale is improved, so that the ground surface settlement deformation information of a higher-precision long-time sequence can be obtained; the method has the advantages that the CWT method is adopted to separate the elastic time sequence deformation signals and the inelastic time sequence deformation signals of the aquifer, more accurate parameter basis is provided for the current situation of difficult separation caused by delay between the elastic time sequence deformation signals and the inelastic InSAR time sequence deformation signals of the aquifer, the problem of poor effect when the related technology processes the unstable signals is solved, the elastic deformation information and the inelastic deformation information of the aquifer can be separated more accurately, independent inversion of the elastic water release coefficients and the inelastic water release coefficients of the aquifer is facilitated, and then the high-precision underground water-ground settlement prediction model is facilitated to be established.
Based on the above embodiments, please refer to fig. 4, which may further include:
s103: and (3) inverting the elastic water release coefficient and the inelastic water release coefficient of the aquifer by utilizing the pre-constructed ground subsidence soil layer compression model.
The ground subsidence soil layer compression model is a model representing the relationship between the ground subsidence and the water release coefficient of the aquifer, and the subsidence of the aquifer can be calculated by using the following formula:
ΔbInSAR＝-miΔtA；
wherein A is the thickness of the aquifer,. DELTA.bInSARIs the amount of compression of the aquifer over time Δ t, SskIs the water release coefficient of the aquifer, the water release coefficient comprises an elastic water release coefficient and an inelastic water release coefficient, and h is the water level value of the aquifer water head.
Calculating an elastic water release coefficient and a non-elastic water release coefficient of the aquifer according to the elastic deformation time sequence information and the non-elastic deformation time sequence information obtained by separation in the S102; when the ground settlement amount (the value predicted by the underground water-ground settlement prediction model) is known, the ground settlement soil layer compression model is utilized to obtain the elastic water release coefficient and the inelastic water release coefficient of the aquifer through inversion, and the accuracy of the underground water-ground settlement prediction model can be verified according to the elastic water release coefficient and the inelastic water release coefficient of the aquifer which are obtained through comparison of the inverted elastic water release coefficient and the inelastic water release coefficient of the aquifer and the elastic water release coefficient and the inelastic water release coefficient of the aquifer which are obtained through calculation. Or inputting the elastic water release coefficient and the inelastic water release coefficient of the aquifer into a ground subsidence soil layer compression model, calculating the subsidence of the ground, and comparing the calculated subsidence with the subsidence predicted by the underground water-ground subsidence prediction model according to the model to verify the accuracy of the underground water-ground subsidence prediction model.
As shown in fig. 5, the abscissa is year (from 2007 to 2009), the ordinate is a ground settlement value, the curve corresponding to InSAR in the graph is a ground settlement value actually monitored by InSAR, the curve corresponding to Model is a ground settlement simulation value (parameters obtained by the inversion method of the present application are substituted into settlement simulation software), and as can be seen from the graph, the accuracy of both is within ± 5mm, which can meet the requirement of monitoring accuracy, so the technical solution proposed in the present application is feasible.
From the above, the embodiment of the invention accurately inverts the elastic and inelastic water release coefficients of the aquifer, and is beneficial to establishing a high-precision underground water-ground settlement prediction model.
The embodiment of the invention also provides a corresponding implementation device for the aquifer deformation signal separation method, so that the method has higher practicability. In the following, the aquifer deformation signal separation device provided by the embodiment of the invention is introduced, and the aquifer deformation signal separation device described below and the aquifer deformation signal separation method described above may be referred to correspondingly.
Referring to fig. 6, fig. 6 is a block diagram of an aquifer deformation signal separation device according to an embodiment of the present invention, in a specific implementation manner, the aquifer deformation signal separation device may include:
the earth surface time sequence deformation sequence acquisition module 601 is used for acquiring an earth surface time sequence deformation sequence of the region to be detected based on the SBAS-InSAR technology;
and the signal separation module 602 is configured to separate elastic deformation timing sequence information and inelastic deformation timing sequence information of the surface timing deformation sequence by using a continuous wavelet transform separation method, so as to realize deformation signal separation of an aquifer in the region to be detected.
Optionally, in some embodiments of this embodiment, referring to fig. 7, the apparatus may further include a coefficient inversion module 603, for example, configured to invert the elastic water release coefficient and the inelastic water release coefficient of the aquifer by using a pre-constructed ground subsidence soil layer compression model.
In some embodiments of this embodiment, the coefficient inversion module 603 may further include a model construction sub-module, for example, configured to calculate the settling amount of the aquifer by using the following formula:
ΔbInSAR＝-miΔtA；
wherein A is the thickness of the aquifer,. DELTA.bInSARIs the amount of compression of the aquifer over a time Δ t, SskAnd h is the water release coefficient of the aquifer, wherein the water release coefficient comprises an elastic water release coefficient and an inelastic water release coefficient, and h is the water level value of the aquifer water head.
Optionally, the surface time-series deformation sequence acquisition module 601 may further include:
the acquisition submodule is used for acquiring all interferograms of the area to be detected acquired by the SBAS-InSAR technology;
the phase unwrapping sub-module is used for unwrapping the wrapping phase of each interference pattern by using a three-dimensional unwrapping method so as to acquire phase information of each interference pattern;
and the time sequence signal generation submodule is used for calculating corresponding average deformation speed according to the phase information of each interference pattern, performing time domain high-pass filtering and space domain low-pass filtering, and generating a surface time sequence deformation sequence of the region to be detected.
In this embodiment, for example, the method may further include:
the high coherence point extraction submodule is used for carrying out self-adaptive filtering on each interference pattern in space and selecting pixels meeting conditions as initial high coherence points;
the deletion submodule is used for calculating the standard deviation of the noise phase of the high coherent points of each interferogram and discarding the high coherent points which do not meet the deviation condition;
and the error correction submodule is used for carrying out error correction on the high-coherence points in each interferogram by utilizing the SRTM3 DEM.
In some other embodiments, the signal separation module 602 may further include:
the separation submodule is used for separating the ground surface time sequence deformation sequence into a long period and a short period by using a continuous wavelet change separation method;
and the intensity calculation submodule is used for calculating the intensity range of each period by utilizing Fourier change and calculating the intensity value of each period of the corresponding basic wave according to each intensity range and the basic wave fitting metadata.
The functions of the functional modules of the aquifer deformation signal separation device according to the embodiment of the present invention may be specifically implemented according to the method in the above method embodiment, and the specific implementation process may refer to the related description of the above method embodiment, and will not be described herein again.
From the above, the embodiment of the invention accurately separates the information of the elastic deformation and the inelastic deformation of the aquifer caused by the elastic compression.
The embodiment of the invention also provides aquifer deformation signal separation equipment, which specifically comprises:
a memory for storing a computer program;
a processor for executing a computer program to implement the steps of the aquifer deformation signal separation method according to any one of the above embodiments.
The functions of the functional modules of the aquifer deformation signal separation device according to the embodiment of the present invention may be specifically implemented according to the method in the above method embodiment, and the specific implementation process may refer to the related description of the above method embodiment, and will not be described herein again.
From the above, the embodiment of the invention accurately separates the information of the elastic deformation and the inelastic deformation of the aquifer caused by the elastic compression.
An embodiment of the present invention further provides a computer-readable storage medium, in which an aquifer deformation signal separation program is stored, where the aquifer deformation signal separation program is executed by a processor, and the method according to any one of the above embodiments is performed.
The functions of the functional modules of the computer-readable storage medium according to the embodiment of the present invention may be specifically implemented according to the method in the foregoing method embodiment, and the specific implementation process may refer to the related description of the foregoing method embodiment, which is not described herein again.
From the above, the embodiment of the invention accurately separates the information of the elastic deformation and the inelastic deformation of the aquifer caused by the elastic compression.
The embodiments are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same or similar parts among the embodiments are referred to each other. The device disclosed by the embodiment corresponds to the method disclosed by the embodiment, so that the description is simple, and the relevant points can be referred to the method part for description.
Those of skill would further appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the various illustrative components and steps have been described above generally in terms of their functionality in order to clearly illustrate this interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
The steps of a method or algorithm described in connection with the embodiments disclosed herein may be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module may reside in Random Access Memory (RAM), memory, Read Only Memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.
The method, the device, the equipment and the computer readable storage medium for separating the aquifer deformation signal provided by the invention are described in detail above. The principles and embodiments of the present invention are explained herein using specific examples, which are presented only to assist in understanding the method and its core concepts. It should be noted that, for those skilled in the art, it is possible to make various improvements and modifications to the present invention without departing from the principle of the present invention, and those improvements and modifications also fall within the scope of the claims of the present invention.
Claims (10)
1. An aquifer deformation signal separation method is characterized by comprising the following steps:
acquiring a ground surface time sequence deformation sequence of a region to be detected based on an SBAS-InSAR technology;
and separating the elastic deformation time sequence information and the inelastic deformation time sequence information of the surface time sequence deformation sequence by adopting a continuous wavelet change separation method so as to realize the separation of deformation signals of the aquifer of the region to be detected.
2. The aquifer deformation signal separation method according to claim 1, further comprising, after the separating the elastic deformation time-series information and the inelastic deformation time-series information of the surface time-series deformation sequence by using the continuous wavelet transform separation method:
and inverting the elastic water release coefficient and the inelastic water release coefficient of the aquifer by utilizing a pre-constructed ground subsidence soil layer compression model.
3. The aquifer deformation signal separation method according to claim 2, wherein the ground subsidence soil layer compression model is as follows:
calculating the settlement of the aquifer using the following formula:
ΔbInSAR＝-miΔtA；
wherein A is the thickness of the aquifer,. DELTA.bInSARIs the amount of compression of the aquifer over a time Δ t, SskAnd h is the water release coefficient of the aquifer, wherein the water release coefficient comprises an elastic water release coefficient and an inelastic water release coefficient, and h is the water level value of the aquifer water head.
4. The aquifer deformation signal separation method according to claim 1, wherein the acquiring of the surface time sequence deformation sequence of the region to be measured based on the SBAS-InSAR technique comprises:
acquiring all interferograms of a region to be detected acquired by the SBAS-InSAR technology;
unwrapping the wrapping phase of each interference pattern by using a three-dimensional unwrapping method to obtain phase information of each interference pattern;
and calculating corresponding average deformation speed according to the phase information of each interference pattern, and performing time domain high-pass filtering and space domain low-pass filtering to generate a surface time sequence deformation sequence of the region to be detected.
5. The aquifer deformation signal separation method according to claim 4, wherein after the acquiring of all interferograms of the area to be measured by the SBAS-InSAR technique, the method further comprises:
carrying out self-adaptive filtering on each interference pattern in space, and selecting pixels meeting conditions as initial high coherence points;
calculating the standard deviation of the noise phase of the high coherence point of each interferogram, and discarding the high coherence points which do not meet the deviation condition;
the high coherence points in each interferogram are error corrected using the SRTM3 DEM.
6. The aquifer deformation signal separation method according to any one of claims 1 to 5, wherein the separating the elastic deformation time sequence information and the inelastic deformation time sequence information of the surface time sequence deformation sequence by using a continuous wavelet transformation separation method comprises:
separating the earth surface time sequence deformation sequence into a long period and a short period by using a continuous wavelet change separation method;
and calculating the intensity range of each period by utilizing the Fourier change, and calculating the intensity value of each period of the corresponding basic wave according to each intensity range and the basic wave fitting metadata.
7. An aquifer deformation signal separation device, comprising:
the system comprises a surface time sequence deformation sequence acquisition module, a time sequence deformation sequence acquisition module and a time sequence deformation sequence acquisition module, wherein the surface time sequence deformation sequence acquisition module is used for acquiring a surface time sequence deformation sequence of a region to be detected based on an SBAS-InSAR technology;
and the signal separation module is used for separating the elastic deformation time sequence information and the inelastic deformation time sequence information of the surface time sequence deformation sequence by adopting a continuous wavelet change separation method so as to realize the deformation signal separation of the aquifer of the region to be detected.
8. The aquifer deformation signal separation device of claim 7, further comprising:
and the coefficient inversion module is used for inverting the elastic water release coefficient and the inelastic water release coefficient of the aquifer by utilizing a pre-constructed ground subsidence soil layer compression model.
9. An aquifer deformation signal separation apparatus comprising a processor for implementing the steps of the aquifer deformation signal separation method according to any one of claims 1 to 6 when executing a computer program stored in a memory.
10. A computer readable storage medium, having stored thereon an aquifer deformation signal separation program, which when executed by a processor, performs the steps of the aquifer deformation signal separation method according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910942684.XA CN110657742A (en) | 2019-09-30 | 2019-09-30 | Aquifer deformation signal separation method, aquifer deformation signal separation device, aquifer deformation signal separation equipment and readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910942684.XA CN110657742A (en) | 2019-09-30 | 2019-09-30 | Aquifer deformation signal separation method, aquifer deformation signal separation device, aquifer deformation signal separation equipment and readable storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110657742A true CN110657742A (en) | 2020-01-07 |
Family
ID=69040272
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910942684.XA Pending CN110657742A (en) | 2019-09-30 | 2019-09-30 | Aquifer deformation signal separation method, aquifer deformation signal separation device, aquifer deformation signal separation equipment and readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110657742A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113175949A (en) * | 2021-04-23 | 2021-07-27 | 首都师范大学 | Method and system for inverting water release coefficient by combining surface deformation and water level information |
| CN114066749A (en) * | 2021-10-21 | 2022-02-18 | 高速铁路建造技术国家工程实验室 | Phase-dependent anti-noise displacement estimation method, equipment and storage medium |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104459692A (en) * | 2014-11-27 | 2015-03-25 | 西安空间无线电技术研究所 | Quick data processing method for improving GEOSAR difference interference deformation measuring accuracy |
| CN108627832A (en) * | 2018-05-11 | 2018-10-09 | 电子科技大学 | A method of passway for transmitting electricity Ground Deformation is extracted based on multiple timings SAR image |
| CN108801211A (en) * | 2018-04-08 | 2018-11-13 | 山东农业大学 | A kind of sinking land in coalmining areas INSAR image method for quickly identifying |
| CN108957456A (en) * | 2018-08-13 | 2018-12-07 | 伟志股份公司 | Landslide monitoring and EARLY RECOGNITION method based on multi-data source SBAS technology |
| CN109100720A (en) * | 2018-09-14 | 2018-12-28 | 长安大学 | A kind of InSAR Ground Deformation monitoring method |
| CN109993336A (en) * | 2017-12-29 | 2019-07-09 | 广州潽蓝信息科技有限公司 | Financial investment data mutation analysis method and its system based on wavelet analysis |
| JP2019132707A (en) * | 2018-01-31 | 2019-08-08 | 株式会社パスコ | Ground surface displacement observation apparatus of tunnel path and ground surface displacement observation program of tunnel path |
-
2019
- 2019-09-30 CN CN201910942684.XA patent/CN110657742A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104459692A (en) * | 2014-11-27 | 2015-03-25 | 西安空间无线电技术研究所 | Quick data processing method for improving GEOSAR difference interference deformation measuring accuracy |
| CN109993336A (en) * | 2017-12-29 | 2019-07-09 | 广州潽蓝信息科技有限公司 | Financial investment data mutation analysis method and its system based on wavelet analysis |
| JP2019132707A (en) * | 2018-01-31 | 2019-08-08 | 株式会社パスコ | Ground surface displacement observation apparatus of tunnel path and ground surface displacement observation program of tunnel path |
| CN108801211A (en) * | 2018-04-08 | 2018-11-13 | 山东农业大学 | A kind of sinking land in coalmining areas INSAR image method for quickly identifying |
| CN108627832A (en) * | 2018-05-11 | 2018-10-09 | 电子科技大学 | A method of passway for transmitting electricity Ground Deformation is extracted based on multiple timings SAR image |
| CN108957456A (en) * | 2018-08-13 | 2018-12-07 | 伟志股份公司 | Landslide monitoring and EARLY RECOGNITION method based on multi-data source SBAS technology |
| CN109100720A (en) * | 2018-09-14 | 2018-12-28 | 长安大学 | A kind of InSAR Ground Deformation monitoring method |
Non-Patent Citations (1)
| Title |
|---|
| 张子文: "时序InSAR地面沉降监测与地下水-地面沉降预测模型参数反演", 《中国博士学位论文全文数据库（电子期刊）》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113175949A (en) * | 2021-04-23 | 2021-07-27 | 首都师范大学 | Method and system for inverting water release coefficient by combining surface deformation and water level information |
| CN113175949B (en) * | 2021-04-23 | 2022-08-02 | 首都师范大学 | Method and system for inverting water release coefficient by combining surface deformation and water level information |
| CN114066749A (en) * | 2021-10-21 | 2022-02-18 | 高速铁路建造技术国家工程实验室 | Phase-dependent anti-noise displacement estimation method, equipment and storage medium |
| CN114066749B (en) * | 2021-10-21 | 2024-08-02 | 高速铁路建造技术国家工程实验室 | Phase correlation anti-noise displacement estimation method, device and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106772342B (en) | Time sequence differential radar interference method suitable for large-gradient ground surface settlement monitoring | |
| Pascazio et al. | Multifrequency InSAR height reconstruction through maximum likelihood estimation of local planes parameters | |
| Fornaro et al. | Imaging of single and double scatterers in urban areas via SAR tomography | |
| Jiang et al. | Hybrid approach for unbiased coherence estimation for multitemporal InSAR | |
| Liu et al. | Estimating Spatiotemporal Ground Deformation With Improved Persistent-Scatterer Radar Interferometry $^\ast$ | |
| Liang et al. | Interferometry with ALOS-2 full-aperture ScanSAR data | |
| US20180011187A1 (en) | Synthetic-aperture radar signal processing apparatus | |
| CN104698457B (en) | InSAR (interferometric synthetic aperture radar) imaging and altitude estimation method by iterative curved surface prediction | |
| Pepe et al. | New advances of the extended minimum cost flow phase unwrapping algorithm for SBAS-DInSAR analysis at full spatial resolution | |
| CN103901442B (en) | A kind of single-frequency multipath effect is surveyed and remission method | |
| CN111239736B (en) | Single-baseline-based surface elevation correction method, device, equipment and storage medium | |
| CN111856459B (en) | Improved DEM maximum likelihood constraint multi-baseline InSAR phase unwrapping method | |
| CN110018476B (en) | Time difference baseline set time sequence interference SAR processing method | |
| WO2020240720A1 (en) | Synthetic aperture radar signal processing device and signal processing method | |
| CN113589287B (en) | Synthetic aperture radar sparse imaging method and device, electronic equipment and storage medium | |
| Mestre-Quereda et al. | An improved phase filter for differential SAR interferometry based on an iterative method | |
| CN110174673A (en) | A method of it is superimposed using timing relay interference pattern and efficiently weakens atmospheric phase influence | |
| CN110657742A (en) | Aquifer deformation signal separation method, aquifer deformation signal separation device, aquifer deformation signal separation equipment and readable storage medium | |
| CN111950438A (en) | Depth learning-based effective wave height inversion method for Tiangong No. two imaging altimeter | |
| CN114812491A (en) | Power transmission line earth surface deformation early warning method and device based on long-time sequence analysis | |
| CN105678716A (en) | Correction method and apparatus for atmospheric interference phase in ground-based SAR | |
| Virtanen et al. | Lag profile inversion method for EISCAT data analysis | |
| Karimi et al. | Mean sea surface and mean dynamic topography determination from CryoSat-2 data around Australia | |
| CN116148855B (en) | Method and system for removing atmospheric phase and resolving deformation of time sequence InSAR | |
| Nies et al. | Phase unwrapping using 2D-Kalman filter-potential and limitations |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200107 |
|
| WD01 | Invention patent application deemed withdrawn after publication |