CN113296137A - Interferometric deformation monitoring method and device and receiver - Google Patents

Abstract

The application provides an interference type deformation monitoring method, an interference type deformation monitoring device and a receiver, and relates to the technical field of remote sensing measurement and control, wherein the method comprises the following steps: the method comprises the steps of firstly obtaining a first interference signal and then obtaining a second interference signal, wherein the first interference signal and the second interference signal are both formed by a direct signal and a reflected signal corresponding to the direct signal, the reflected signal corresponding to the direct signal is a signal obtained by reflecting the direct signal by a measured surface, the direct signal is a Global Navigation Satellite System (GNSS) signal, the emission angle of the direct signal forming the first interference signal is the same as that of the direct signal forming the second interference signal, the emission time is different, and finally, the deformation quantity of the measured surface is determined according to a first carrier-to-noise ratio sequence of the first interference signal and a second carrier-to-noise ratio sequence of the second interference signal. The method and the device do not need to continuously detect interference signals, so that the reliability of the receiver can be effectively improved.

Description

Interferometric deformation monitoring method, device and receiver

technical field

The application belongs to the technical field of remote sensing measurement and control, and in particular relates to an interferometric deformation monitoring method, device and receiver.

Background technique

The deformation monitoring technology based on Global Navigation Satellite System-Reflectometry (GNSS-R) has the advantages of high measurement accuracy, low cost and high operational safety, and can be applied to geological surveys, landslides and buildings detection and other fields.

GNSS-R deformation monitoring technology belongs to a remote sensing technology. The receiver receives the direct signal emitted by the satellite and the reflected signal generated by the reflection of the direct signal from the measured surface, and integrates the direct signal and the reflected signal to obtain the measured signal. shape of the surface. However, the signal strength of the reflected signal generated by the reflection of the measured surface is usually relatively low, so that the receiver needs to have a strong signal tracking ability to continuously receive the reflected signal of good quality. If the signal tracking ability of the receiver is insufficient, it is easy to cause the problem that the integral calculation cannot be processed due to the loss of the reflected signal, and the deformation of the measured surface cannot be monitored, so the reliability is not high enough.

SUMMARY OF THE INVENTION

In view of this, embodiments of the present application provide an interferometric deformation monitoring method, device, and receiver, which are used to improve the reliability of the receiver.

In a first aspect, an embodiment of the present application provides an interference-type deformation monitoring method, including:

obtain the first interference signal;

Obtain the second interference signal, wherein the first interference signal and the second interference signal are both formed by the direct signal and the reflected signal corresponding to the direct signal, and the reflected signal corresponding to the direct signal is the signal after the direct signal is reflected by the measured surface, The direct signal is a GNSS signal of the Global Navigation Satellite System, and the transmission angle of the direct signal forming the first interference signal is the same as that of the direct signal forming the second interference signal, but the transmission time is different;

The deformation amount of the measured surface is determined according to the first carrier-to-noise ratio sequence of the first interference signal and the second carrier-to-noise ratio sequence of the second interference signal.

Optionally, determining the deformation amount of the measured surface according to the first carrier-to-noise ratio sequence of the first interference signal and the second carrier-to-noise ratio sequence of the second interference signal, including:

determining a first phase difference according to the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence, where the first phase difference is the phase difference between the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence;

Determine the equivalent height angle of the measured surface according to the preset navigation message and the preset measured surface data, and the equivalent height angle is the complementary angle of the incident angle of the direct signal on the measured surface;

The deformation amount of the measured surface is determined according to the first phase difference and the equivalent height angle of the measured surface.

Optionally, determine the equivalent altitude angle of the measured surface according to the preset navigation message and the preset measured surface data, including:

Determine the satellite altitude and satellite azimuth of the measured surface according to the navigation message and the receiving position of the direct signal;

According to the satellite altitude angle, satellite azimuth angle and the measured surface data of the measured surface, the equivalent altitude angle of the measured surface is determined.

Optionally, the measured surface data includes the inclination angle of the measured surface, the azimuth angle of the measured surface, and the vertical distance from the position where the reflected signal is received to the measured surface.

Optionally, according to the satellite altitude angle, satellite azimuth angle and the measured surface data of the measured surface, determine the equivalent altitude angle of the measured surface, including:

Use the following formula to determine the equivalent height angle of the measured surface:

Among them, β represents the equivalent elevation angle of the measured surface, β″ represents the projection angle of the equivalent elevation angle of the measured surface on the first reference plane, θ′ represents the projection angle of the satellite elevation angle on the first reference plane, γ represents the inclination angle of the measured surface, α represents the difference between the satellite azimuth angle and the measured surface azimuth angle, α′ represents the projection angle of the difference between the satellite azimuth angle and the measured surface azimuth angle on the second reference plane, θ represents the satellite elevation angle, α _s represents the satellite azimuth angle, α _r represents the azimuth angle of the measured surface, the first reference plane is perpendicular to the measured surface and the ground, and the reflected signal receiving position is located in the first reference plane, the second reference plane The planes are respectively perpendicular to the measured surface and the first reference plane, and the reflection position of the measured surface is located in the second reference plane.

Optionally, determine the deformation amount of the measured surface according to the first phase difference and the equivalent height angle of the measured surface, including:

Use the following formula to determine the deformation of the measured surface:

Among them, d _def represents the deformation of the measured surface, β represents the equivalent height angle of the measured surface, ΔΦ represents the first phase difference, and λ represents the wavelength in the direct signal.

In a second aspect, an embodiment of the present application provides a communication device for a receiver, including:

an obtaining module, configured to obtain a first interference signal and obtain a second interference signal, wherein the first interference signal and the second interference signal are both formed by the direct signal and the reflected signal corresponding to the direct signal, and the reflected signal corresponding to the direct signal It is the signal after the direct signal is reflected by the measured surface. The direct signal is the GNSS signal of the Global Navigation Satellite System. The emission angle of the direct signal forming the first interference signal is the same as the emission angle of the direct signal forming the second interference signal, but the emission time is different. ;

The determining module is configured to determine the deformation amount of the measured surface according to the first carrier-to-noise ratio sequence of the first interference signal and the second carrier-to-noise ratio sequence of the second interference signal.

Optionally, determine that the module is specifically used for:

determining a first phase difference according to the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence, where the first phase difference is the phase difference between the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence;

Determine the equivalent height angle of the measured surface according to the preset navigation message and the preset measured surface data, and the equivalent height angle is the complementary angle of the incident angle of the direct signal on the measured surface;

The deformation amount of the measured surface is determined according to the first phase difference and the equivalent height angle of the measured surface.

Optionally, determine that the module is specifically used for:

Determine the satellite altitude and satellite azimuth of the measured surface according to the navigation message and the receiving position of the direct signal;

According to the satellite altitude angle, satellite azimuth angle and the measured surface data of the measured surface, the equivalent altitude angle of the measured surface is determined.

Optionally, the measured surface data includes the inclination angle of the measured surface, the azimuth angle of the measured surface, and the vertical distance from the position where the reflected signal is received to the measured surface.

Optionally, determine that the module is specifically used for:

Use the following formula to determine the equivalent height angle of the measured surface:

Among them, β represents the equivalent elevation angle of the measured surface, β″ represents the projection angle of the equivalent elevation angle of the measured surface on the first reference plane, θ′ represents the projection angle of the satellite elevation angle on the first reference plane, γ represents the inclination angle of the measured surface, α represents the difference between the satellite azimuth angle and the measured surface azimuth angle, α′ represents the projection angle of the difference between the satellite azimuth angle and the measured surface azimuth angle on the second reference plane, θ represents the satellite elevation angle, α _s represents the satellite azimuth angle, α _r represents the azimuth angle of the measured surface, the first reference plane is perpendicular to the measured surface and the ground, and the reflected signal receiving position is located in the first reference plane, the second reference plane The planes are respectively perpendicular to the measured surface and the first reference plane, and the reflection position of the measured surface is located in the second reference plane.

Optionally, determine that the module is specifically used for:

Use the following formula to determine the deformation of the measured surface:

Among them, d _def represents the deformation of the measured surface, β represents the equivalent height angle of the measured surface, ΔΦ represents the first phase difference, and λ represents the wavelength in the direct signal.

In a third aspect, an embodiment of the present application provides a receiver, including a memory, a processor, and a computer program stored in the memory and running on the processor, and the processor implements the first aspect or the first aspect when the processor executes the computer program. The method of any embodiment of the aspect.

In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium, including a computer-readable storage medium storing a computer program, and when the computer program is executed by a processor, the first aspect or any implementation manner of the first aspect is implemented Methods.

The present application provides an interferometric deformation monitoring method, device and receiver. First, a first interference signal can be obtained, and then a second interference signal can be obtained. The reflected signal corresponding to the signal is formed by the reflected signal corresponding to the direct signal. The reflected signal corresponding to the direct signal is the signal after the direct signal is reflected by the measured surface, and the direct signal is the GNSS signal of the Global Navigation Satellite System. The direct signals of the two interference signals have the same emission angle and different emission times. Finally, the deformation amount of the measured surface is determined according to the first carrier-to-noise ratio sequence of the first interference signal and the second carrier-to-noise ratio sequence of the second interference signal. The present application can determine the deformation of the measured surface through the relationship between the carrier-to-noise ratio sequence of the interference signal and the carrier phase difference between the direct signal and the reflected signal. Since it is not necessary to continuously detect the interference signal, the reliability of the receiver can be effectively improved. sex.

Description of drawings

In order to illustrate the technical solutions in the embodiments of the present application more clearly, the following briefly introduces the accompanying drawings that need to be used in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only for the present application. In some embodiments, for those of ordinary skill in the art, other drawings can also be obtained according to these drawings without any creative effort.

FIG. 1 is a schematic flow chart of an interference-type deformation monitoring method provided by an embodiment of the present application;

2 is a schematic diagram of a receiver application scenario provided by an embodiment of the present application;

3 is a two-dimensional schematic diagram of a signal reflection model provided by an embodiment of the present application;

4 is a three-dimensional schematic diagram of a GNSS reflection geometric model provided by an embodiment of the present application;

5 is a partial enlarged view of a GNSS reflection geometric model provided by an embodiment of the present application;

Fig. 6 is the flow chart of determining the shape variable of the measured surface provided by the embodiment of the present application;

7 is a structural block diagram of an interferometric deformation monitoring device provided by an embodiment of the present application;

FIG. 8 is a schematic structural diagram of a receiver provided by an embodiment of the present application.

Detailed ways

In the following description, for the purpose of illustration rather than limitation, specific details such as a specific system structure and technology are set forth in order to provide a thorough understanding of the embodiments of the present application. However, it will be apparent to those skilled in the art that the present application may be practiced in other embodiments without these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present application with unnecessary detail.

It is to be understood that, when used in this specification and the appended claims, the term "comprising" indicates the presence of the described feature, integer, step, operation, element and/or component, but does not exclude one or more other The presence or addition of features, integers, steps, operations, elements, components and/or sets thereof.

It will also be understood that, as used in this specification and the appended claims, the term "and/or" refers to and including any and all possible combinations of one or more of the associated listed items.

In addition, in the description of the specification of the present application and the appended claims, the terms "first", "second", "third", etc. are only used to distinguish the description, and should not be construed as indicating or implying relative importance.

The Global Navigation Satellite System (GNSS) is a space-based radio navigation and positioning system that can provide users with all-weather three-dimensional coordinates, speed and time information at any location on the earth's surface or near-Earth space. Real-time and high precision and many other advantages. The GNSS system may specifically be a Global Positioning System (Global Positioning System, GPS), a GLONASS navigation system (GLONASS), a Galileo (Galileo) system, a Beidou satellite navigation system, or the like. In this embodiment of the present application, the receiver may receive a GNSS signal of any GNSS system, which is not particularly limited in this embodiment.

GNSS-R deformation monitoring technology belongs to a remote sensing technology. The main working method is that the receiver receives the direct signal transmitted by the satellite. Since the direct signal transmitted by the satellite covers a wide area, some of the direct information will be reflected by the measured surface to form a reflection. The signal and reflected signal will be affected by the measured surface to produce some changes, and then the receiver receives the reflected signal. Finally, by comparing the direct signal and the reflected signal, the information of the different parts of the reflected signal and the direct signal is obtained, and the information is analyzed to obtain the measured surface. certain information. However, the signal strength of the reflected signal generated by the reflection of the measured surface is usually relatively low, so that the receiver needs to have a strong signal tracking ability to continuously receive the reflected signal of good quality. If the signal tracking ability of the receiver is insufficient, it is easy to cause the problem that the integral calculation cannot be processed due to the loss of the reflected signal, and the deformation of the measured surface cannot be monitored, so the reliability is not high enough.

Therefore, the embodiments of the present application provide a technical solution that does not require continuous detection of direct signals and reflected signals, and can effectively improve the reliability of the receiver.

The receiver provided in this embodiment of the present application may be a software receiver, a hardware receiver, a single-frequency receiver, or a dual-frequency receiver, and other equipment having a function of receiving GNSS signals.

The technical solutions of the present application will be described in detail below with specific examples. The following specific embodiments may be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments.

FIG. 1 is a schematic flowchart of an interference-type deformation monitoring method provided by an embodiment of the present application. As shown in FIG. 1 , the method may include the following steps:

S110. Obtain a first interference signal.

In this embodiment of the present application, the receiver may receive the direct signal of the GNSS through the right-handed antenna, receive the reflected signal of the GNSS through the left-handed antenna, and then pass the direct signal and the reflected signal through a combiner to form an interference signal, wherein the reflected signal is the direct signal The signal reflected by the measured surface, the direct signal is the GNSS signal of the global navigation satellite system. The receiver may determine the first-obtained interference signal as the first interference signal.

In one embodiment, the receiver may receive the direct signal of the GNSS through the right-handed antenna, and receive the reflected signal of the GNSS through the left-handed antenna. FIG. 2 is a schematic diagram of an application scenario of a receiver provided by an embodiment of the present application. As shown in FIG. 2 , two direct signals are sent by the same satellite, and one of the direct signals is emitted to the measured surface of the building 1, and then passes through the measured surface 1. reflection, producing a reflected signal. In the case of obtaining the direct signal and the reflected signal, the receiver can combine the direct signal and the reflected signal into an interference signal through a combiner.

It should be noted that the receiver can select the frequency components it needs from the numerous electromagnetic waves existing in the air, suppress or filter out unwanted signals, noise or interference signals, and then obtain the original useful information through amplification and demodulation. Therefore, in this embodiment of the present application, the receiver can parse the acquired GNSS signal to obtain information such as carrier information, pseudo-range code information, and navigation message, and the specific parsing steps will not be repeated in this application.

S120. Obtain a second interference signal.

After obtaining the first interference signal, the receiver may obtain the second interference signal when the emission angle of the direct signal forming the first interference signal and the emission angle of the direct signal forming the second interference signal are the same.

In the embodiment of the present application, since the reflected signal is generated by the reflection of the measured surface, the deformation of the measured surface will cause corresponding changes in the reflected signal, so the receiver can obtain the deformation change of the measured surface through the interference signal.

Specifically, in the monitoring model composed of the satellite, the measured surface and the receiver, since the satellite is too far away from the measured surface and the receiver, the position change of the satellite relative to the measured surface can be approximately considered as the direct signal emission angle relative to the Changes in the measured surface. Further, the satellite moves around the earth in a fixed period, so the change of the direct signal emission angle relative to the measured surface also changes periodically, so at the same time in different cycles, the direct signal emission angle relative to the measured surface changes. Also the same. Among them, in different satellite systems, the movement period of the satellites is also different: the revisit period of the GLONASS satellite is 8 days, the revisit period of the Galileo satellite is 10 days, and the revisit period of the inclined geosynchronous orbit (IGSO) satellite of the Beidou system. It is 1 day, and the medium circular orbit (MEO) satellite revisit period is 7 days.

For example, if the revisit period of the satellite is 24 hours, and the receiver obtains the first interference signal between 12:00 and 12:30 on the first day, the receiver can obtain the first interference signal between 12:00 and 12:30 on the third day. The second interference signal is obtained between 12:00 and 12:30 on the first day, because the emission angle of the direct signal relative to the measured surface between 12:00 and 12:30 on the first day is different from that between 12:00 and 12:30 on the third day. The emission angle to the measured surface is the same.

It should be noted that, the first interference signal and the second interference signal described in the embodiments of the present application are defined with respect to the timing relationship of signal acquisition in a deformation monitoring process.

In the embodiment of the present application, the carrier phase of the interference signal is not used, so it is not necessary to continuously obtain the interference signal, so that the receiver can stop receiving the signal after obtaining enough interference signal, thereby reducing the signal tracking of the receiver. capacity requirements and monitoring time, effectively improving the reliability of the receiver.

S130. Determine the deformation amount of the measured surface according to the first carrier-to-noise ratio sequence of the first interference signal and the second carrier-to-noise ratio sequence of the second interference signal.

In order to facilitate understanding, the principle of determining the measured surface deformation in this step is introduced first.

First, the GNSS signal T(t) emitted by each satellite can be expressed by the following formula:

表示GNSS信号的载波，y(t)表示伪距码，d(t)表示导航电文。Among them, A _T represents the amplitude of the transmitted signal, f represents the carrier frequency of the GNSS signal,

Represents the carrier of the GNSS signal, y(t) represents the pseudo-range code, and d(t) represents the navigation message.

After the GNSS signal transmitted by the satellite passes through the atmosphere and is finally acquired by the receiver (here refers to the direct signal R _d (t)), it can be expressed by the following formula:

表示接收到的直射信号的载波，A_d表示直射信号的幅度，n(t)表示观测噪声。Among them, τ _d is the transmission time of the direct signal from the satellite to the receiver, f(t-τ _d ) is the carrier frequency of the direct signal, f _Dd (t-τ _d ) is the Doppler frequency shift of the direct signal,

represents the carrier of the received direct signal, Ad represents the amplitude of the direct signal, and n( _t ) represents the observation noise.

The direct signal is reflected by the measured surface to generate a reflected signal, and the reflected signal R _r (t) can be expressed by the following formula:

表示反射信号的载波相位，A_r表示反射信号的幅度。where τ _r is the transit time of the reflected signal from the satellite to the receiver, f(t-τ _r ) is the carrier frequency of the reflected signal, f _Dr (t-τ _r ) is the Doppler frequency shift of the reflected signal,

represents the carrier phase of the reflected signal, and _Ar represents the amplitude of the reflected signal.

The direct signal and the reflected signal can generate an interference signal through the combiner. Both the direct signal and the reflected signal come from the same satellite, so the frequency between the direct signal and the reflected signal can be regarded as the same, so according to the principle of sine wave superposition, the interference The signal can still be represented as a single sine wave. The interference signal R _i (t) can be expressed by the following formula:

表示干涉信号的载波相位，且d(t)＝d(t-τ_d)d(t-τ_r)，y(t)＝y(t-τ_d)y(t-τ_r)。Among them, A _i represents the amplitude of the interference signal, f(t-τ _i ) represents the carrier frequency of the interference signal,

represents the carrier phase of the interference signal, and d(t)=d(t-τ _d )d(t-τ _r ), y(t)=y(t-τ _d )y(t-τ _r ).

Further, the amplitude A _i of the interference signal can be expressed by the following formula:

表示直射信号和反射信号的载波相位差。where C _Δ (τ _d -τ _r ) represents the autocorrelation function,

Indicates the carrier phase difference between the direct signal and the reflected signal.

的大小与

(A_rC_Δ(τ_d-τ_r))²以及2A_dA_rC_Δ(τ_d-τ_r)有关系，而幅值

的震荡波形则与载波相位差

有关系，由于在短时间内

(A_rC_Δ(τ_d-τ_r))²以及2A_dA_rC_Δ(τ_d-τ_r)可看作是常量，因此接收机可以根据幅值

的震荡波形确定载波相位差

Due to the movement of the satellite, the carrier phase difference between the direct signal and the reflected signal will change all the time. According to formula (5), it can be known that the amplitude of the interference signal

size with

(A _r C _Δ (τ _d -τ _r )) ² and 2A _d A _r C _Δ (τ _d -τ _r ) are related, and the amplitude

The oscillation waveform is out of phase with the carrier

related, due to the short period of time

(A _r C _Δ (τ _d -τ _r )) ² and 2A _d A _r C _Δ (τ _d -τ _r ) can be regarded as constants, so the receiver can

The oscillator waveform determines the carrier phase difference

Further, the amplitude of the signal in the receiver is also the strength of the signal, which is usually represented by the carrier-to-noise ratio C/N ₀ , and the unit is dB-Hz. ratio. Specifically, the receiver has a certain distance from the measured surface, so the pseudo-range codes of the direct signal and the reflected signal are different. Therefore, the receiver can distinguish the pseudo-range code phase of the direct signal and the reflected signal according to the autocorrelation function. The computer can determine the peak value of the interference signal through the phase-locked loop, and calculate the carrier-to-noise ratio sequence, that is, the oscillation waveform of the amplitude.

Further, according to the relationship between phase, distance and wavelength, the receiver can subtract the distance traveled by the reflected signal from the distance traveled by the direct signal to determine the distance traveled by the reflected signal more than the direct signal. FIG. 3 is a two-dimensional schematic diagram of a signal reflection model provided by an embodiment of the present application. As shown in Figure 3, point A is the position of the receiver, the direct signal is parallel to the reflected signal, DO represents the measured surface, the direct signal reaches point A after being reflected by the measured surface at point O, and point B is the mirror point of point A , so the COA is the distance that the reflected signal travels more than the direct signal (it can be called the propagation path difference).

The propagation path difference Δd can be expressed by the following formula:

where λ represents the wavelength of the GNSS signal.

Further, the receiver can model the GNSS reflection geometry. FIG. 4 is a three-dimensional schematic diagram of a GNSS reflection geometric model provided by an embodiment of the present application, and FIG. 5 is a partial enlarged view of the GNSS reflection geometric model provided by an embodiment of the present application. As shown in FIG. 4 and FIG. There is a certain geometric relationship between the measuring surfaces, and the change of the propagation path difference Δd can also be calculated by analyzing the geometric relationship between the satellite, the receiver and the measured surface. Among them, d represents the vertical distance from the receiver to the measured surface, β represents the equivalent altitude angle of the measured surface, β″ represents the projection angle of the equivalent altitude angle of the measured surface on the first reference plane, and θ′ represents the satellite The projection angle of the elevation angle on the first reference plane, γ represents the inclination angle of the measured surface, α represents the difference between the azimuth angle of the satellite and the azimuth angle of the measured surface, and α′ represents the difference between the azimuth angle of the satellite and the azimuth angle of the measured surface. The projection angle of the difference on the second reference plane, θ represents the satellite elevation angle, α _s represents the satellite azimuth angle, α _r represents the azimuth angle of the measured surface, the first reference plane is perpendicular to the measured surface and the ground, and the reflected signal is The receiving position is located in the first reference plane, the second reference plane is perpendicular to the measured surface and the first reference plane, respectively, and the measured surface reflection position is located in the second reference plane.

Specifically, in conjunction with Figure 3, Figure 4 and Figure 5, the following formula can be obtained through geometric analysis:

Among them, the satellite altitude angle θ and satellite azimuth angle α _s can be calculated from the navigation text in the GNSS signal and the position of the receiver, the inclination angle γ of the measured surface, the azimuth angle α _r of the measured surface and the vertical distance d can be calculated. measured in advance.

Therefore, the propagation path difference Δd can be determined by formula (7), and the following formula can be obtained by combining formulas (6) and (7):

并且对于同一个被测面，当卫星等效高度角相同时，直射信号和反射信号的载波相位差

也相同，所以干涉信号的载噪比序列也相同。The equivalent altitude angle also affects the carrier phase difference between the direct signal and the reflected signal

And for the same measured surface, when the equivalent altitude angle of the satellite is the same, the carrier phase difference between the direct signal and the reflected signal is

is also the same, so the carrier-to-noise ratio sequence of the interference signal is also the same.

When the measured surface is deformed, the emission angle of the direct signal is the same as when the measured surface is not deformed, then the formula (8) can be deformed into the following formula:

表示在被测面发生形变后的直射信号和反射信号的载波相位差，d_def表示被测面的形变量。in,

Represents the carrier phase difference between the direct signal and the reflected signal after the measured surface is deformed, and d _def represents the deformation of the measured surface.

Further, by subtracting formula (8) and formula (9), and combining formula (5), the following formula can be obtained:

Among them, ΔΦ is the phase difference between the carrier-to-noise ratio sequence before the measured surface is deformed and the carrier-to-noise ratio sequence after the measured surface is deformed. It should be noted that when the receiver obtains two carrier-to-noise ratio sequences, the phase difference of the two carrier-to-noise ratio sequences can be directly obtained, which further simplifies the calculation steps.

Further, by transforming formula (10), the following formula can be obtained:

After the above theoretical analysis, the receiver can obtain the carrier-to-noise ratio sequence of the interference signal under the condition of obtaining the interference signal with the same emission angle of the two direct signals, and obtain the phase difference by comparing the two carrier-to-noise ratio sequences. The phase difference, the equivalent height angle and the wavelength of the direct signal are input into formula (11) to obtain the deformation of the measured surface.

Specifically, FIG. 6 is a flowchart of determining the deformation variable of the measured surface provided by the embodiment of the present application. As shown in FIG. 6 , the receiver can determine the deformation variable of the measured surface through the following steps:

S131. Determine the first phase difference according to the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence.

The receiver may determine the peak value of the first interference signal through a phase-locked loop under the condition of acquiring the first interference signal, and obtain the first carrier-to-noise ratio sequence by calculation. Similarly, the receiver can also obtain the second carrier-to-noise ratio sequence. Then the receiver can determine the first phase difference according to the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence, where the first phase difference is the phase difference between the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence.

Specifically, because the carrier-to-noise ratio sequence has relatively large noise, the receiver can first use the orthogonal fitting method to fit a sine wave to the two carrier-to-noise ratio sequences, and then obtain the phase spectrum of the sine wave through Fourier transform, and finally do The difference is the phase difference, that is, the first phase difference.

S132: Determine the equivalent height angle of the measured surface according to the preset navigation message and the preset measured surface data.

Usually, the satellite will carry the latest navigation message in the direct signal when transmitting the direct signal, but the latest navigation message is often not corrected, so there is a certain error with the real situation. After receiving the navigation message, the satellite-related agency will correct it to improve the accuracy of the navigation message. Therefore, the use of the navigation message issued by the satellite-related agency can improve the accuracy of the monitoring results.

In the embodiment of the present application, since there is no requirement for the continuity of the received signal, the receiver can obtain a more accurate navigation message through the network or other means after obtaining the interference signal for several days.

Specifically, the receiver can determine the satellite altitude angle and satellite azimuth angle of the measured surface according to the receiving position of the navigation message and the direct signal. Among them, in practical applications, the direct signal receiving position is the actual position of the receiver.

Further, the receiver can determine the equivalent altitude angle of the measured surface according to the satellite altitude angle, satellite azimuth angle and measured surface data of the measured surface. The measured surface data includes the inclination angle of the measured surface, the azimuth angle of the measured surface, and the vertical distance from the direct signal receiving position to the measured surface.

Specifically, the receiver can input the satellite altitude angle, satellite azimuth angle and the measured surface data of the measured surface into formula (7) to obtain the equivalent altitude angle of the measured surface. Complementary to the angle of incidence of the face.

It should be noted that, when the equivalent height angle is specifically determined, the specific angle value of the equivalent height angle can be determined according to formula (7), or the sin value of the equivalent height angle can be directly determined without calculating the angle value.

S133: Determine the deformation amount of the measured surface according to the first phase difference and the equivalent height angle of the measured surface.

After obtaining the first phase difference and the equivalent height angle of the measured surface, the receiver can input the above data into formula (11) to obtain the deformation value of the measured surface.

In the interferometric deformation monitoring method provided in the embodiment of the present application, a first interference signal can be obtained, and then a second interference signal can be obtained, wherein the first interference signal and the second interference signal are both formed by the direct signal and the reflected signal corresponding to the direct signal. The reflected signal corresponding to the direct signal is the signal after the direct signal is reflected by the measured surface, the direct signal is the GNSS signal of the global navigation satellite system, the emission angle of the direct signal forming the first interference signal and the direct signal forming the second interference signal The emission angles of the 1 and 2 are the same, but the emission times are different. Finally, the deformation amount of the measured surface is determined according to the first carrier-to-noise ratio sequence of the first interference signal and the second carrier-to-noise ratio sequence of the second interference signal. The present application can determine the deformation of the measured surface through the relationship between the carrier-to-noise ratio sequence of the interference signal and the carrier phase difference between the direct signal and the reflected signal. Since it is not necessary to continuously detect the interference signal, the reliability of the receiver can be effectively improved. sex.

It should be understood that the size of the sequence numbers of the steps in the above embodiments does not mean the sequence of execution, and the execution sequence of each process should be determined by its function and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present application.

FIG. 7 is a structural block diagram of an interferometric deformation monitoring device provided by an embodiment of the present application. As shown in FIG. 7 , the device may include:

The obtaining module 110 is configured to obtain a first interference signal and obtain a second interference signal, wherein the first interference signal and the second interference signal are both formed by the direct signal and the reflected signal corresponding to the direct signal, and the reflection corresponding to the direct signal The signal is the signal of the direct signal reflected by the measured surface, and the direct signal is the GNSS signal of the global navigation satellite system. The emission angle of the direct signal forming the first interference signal is the same as that of the direct signal forming the second interference signal, and the emission time is the same. different;

The determining module 110 is configured to determine the deformation amount of the measured surface according to the first carrier-to-noise ratio sequence of the first interference signal and the second carrier-to-noise ratio sequence of the second interference signal.

Optionally, the determining module 110 is specifically used for:

determining a first phase difference according to the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence, where the first phase difference is the phase difference between the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence;

Determine the equivalent height angle of the measured surface according to the preset navigation message and the preset measured surface data, and the equivalent height angle is the complementary angle of the incident angle of the direct signal on the measured surface;

The deformation amount of the measured surface is determined according to the first phase difference and the equivalent height angle of the measured surface.

Optionally, the determining module 110 is specifically used for:

Determine the satellite altitude and satellite azimuth of the measured surface according to the navigation message and the receiving position of the direct signal;

According to the satellite altitude angle, satellite azimuth angle and the measured surface data of the measured surface, the equivalent altitude angle of the measured surface is determined.

Optionally, the measured surface data includes the inclination angle of the measured surface, the azimuth angle of the measured surface, and the vertical distance from the position where the reflected signal is received to the measured surface.

Optionally, the determining module 110 is specifically used for:

Use the following formula to determine the equivalent height angle of the measured surface:

Among them, β represents the equivalent elevation angle of the measured surface, β″ represents the projection angle of the equivalent elevation angle of the measured surface on the first reference plane, θ′ represents the projection angle of the satellite elevation angle on the first reference plane, γ represents the inclination angle of the measured surface, α represents the difference between the satellite azimuth angle and the measured surface azimuth angle, α′ represents the projection angle of the difference between the satellite azimuth angle and the measured surface azimuth angle on the second reference plane, θ represents the satellite elevation angle, α _s represents the satellite azimuth angle, α _r represents the azimuth angle of the measured surface, the first reference plane is perpendicular to the measured surface and the ground, and the reflected signal receiving position is located in the first reference plane, the second reference plane The planes are respectively perpendicular to the measured surface and the first reference plane, and the reflection position of the measured surface is located in the second reference plane.

Optionally, the determining module 110 is specifically used for;

Use the following formula to determine the deformation of the measured surface:

Among them, d _def represents the deformation of the measured surface, β represents the equivalent height angle of the measured surface, ΔΦ represents the first phase difference, and λ represents the wavelength in the direct signal.

FIG. 8 is a schematic structural diagram of a receiver provided by an embodiment of the present application. As shown in FIG. 8 , the receiver of this embodiment includes: at least one processor 20 (only one is shown in FIG. 8 ), a memory 21 , and storage in the memory A computer program 22 in 21 that can be run on at least one processor 20. When the processor 20 executes the computer program 22, the steps in any of the foregoing embodiments of the receiver control method are implemented.

The receiver can be a software receiver, a hardware receiver, a single-frequency receiver or a dual-frequency receiver. Those skilled in the art can understand that FIG. 8 is only an example of a receiver, and does not constitute a limitation on the receiver, and may include more or less components than those shown in the figure, or combine some components, or different components, such as It can also include input and output devices, network access devices, and the like.

The so-called processor 20 may be a central processing unit (Central Processing Unit, CPU), and the processor 20 may also be other general-purpose processors, digital signal processors (Digital Signal Processors, DSP), application specific integrated circuits (Application Specific Integrated Circuits) , ASIC), off-the-shelf programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.

The memory 21 may in some embodiments be an internal storage unit of the receiver, such as a hard disk or memory of the receiver. In other embodiments, the memory 21 may also be an external storage device of the receiver, such as a plug-in hard disk equipped on the receiver, a Smart Media Card (SMC), a Secure Digital (SD) card, a flash memory Card (Flash Card) and so on. Further, the memory 21 may also include both an internal storage unit of the receiver and an external storage device. The memory 21 is used to store an operating system, an application program, a boot loader (Boot Loader), data, and other programs, such as program codes of computer programs, and the like. The memory 21 can also be used to temporarily store data that has been output or is to be output.

Embodiments of the present application further provide a computer-readable storage medium, where a computer program is stored in the computer-readable storage medium, and when the computer program is executed by a processor, the steps in the foregoing method embodiments can be implemented.

It should be noted that the information exchange, execution process and other contents between the above-mentioned devices/units are based on the same concept as the method embodiments of the present application. For specific functions and technical effects, please refer to the method embodiments section. It is not repeated here.

Those skilled in the art can clearly understand that, for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example. Module completion, that is, dividing the internal structure of the device into different functional units or modules to complete all or part of the functions described above. Each functional unit and module in the embodiment may be integrated in one processing unit, or each unit may exist physically alone, or two or more units may be integrated in one unit, and the above-mentioned integrated units may adopt hardware. It can also be realized in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing from each other, and are not used to limit the protection scope of the present application. For the specific working processes of the units and modules in the above-mentioned system, reference may be made to the corresponding processes in the foregoing method embodiments, which will not be repeated here.

Those of ordinary skill in the art can realize that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Skilled artisans may implement the described functionality using different methods for each particular application, but such implementations should not be considered beyond the scope of this application.

The above-mentioned embodiments are only used to illustrate the technical solutions of the present application, but not to limit them; although the present application has been described in detail with reference to the above-mentioned embodiments, those of ordinary skill in the art should understand that: it can still be used for the above-mentioned implementations. The technical solutions described in the examples are modified, or some technical features thereof are equivalently replaced; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions in the embodiments of the application, and should be included in the within the scope of protection of this application.

Claims (10)

1. an interference-type deformation monitoring method, is characterized in that, comprises: obtain the first interference signal; A second interference signal is obtained, wherein the first interference signal and the second interference signal are both formed by a direct signal and a reflected signal corresponding to the direct signal, and the reflected signal corresponding to the direct signal is obtained by the direct signal. The signal reflected by the surface to be measured, the direct signal is a GNSS signal of the Global Navigation Satellite System, and the emission angle of the direct signal forming the first interference signal is the same as the emission angle of the direct signal forming the second interference signal. different time; The deformation amount of the measured surface is determined according to the first carrier-to-noise ratio sequence of the first interference signal and the second carrier-to-noise ratio sequence of the second interference signal. 2 . The interferometric deformation monitoring method according to claim 1 , wherein, according to the first carrier-to-noise ratio sequence of the first interference signal and the second carrier-to-noise ratio sequence of the second interference signal, 2 . Determine the deformation of the measured surface, including: Determine a first phase difference according to the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence, where the first phase difference is the sum of the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence phase difference between Determine the equivalent height angle of the measured surface according to the preset navigation message and the preset measured surface data, and the equivalent height angle is the complementary angle of the incident angle of the direct signal on the measured surface ; The deformation amount of the measured surface is determined according to the first phase difference and the equivalent height angle of the measured surface. 3. The interferometric deformation monitoring method according to claim 2, wherein, determining the equivalent height angle of the measured surface according to a preset navigation message and preset measured surface data, comprising: Determine the satellite altitude and satellite azimuth of the measured surface according to the navigation message and the receiving position of the direct signal; The equivalent altitude angle of the measured surface is determined according to the satellite altitude angle, satellite azimuth angle and the measured surface data of the measured surface. 4 . The interferometric deformation monitoring method according to claim 3 , wherein the measured surface data includes the inclination angle of the measured surface, the azimuth angle of the measured surface, and the received reflected signal. 5 . The vertical distance from the position to the measured surface. 5 . The method for monitoring interferometric deformation according to claim 4 , wherein the measurement of the measured surface is determined according to the satellite altitude angle, satellite azimuth angle of the measured surface and the measured surface data. 6 . Equivalent elevation angles, including: Use the following formula to determine the equivalent height angle of the measured surface:

Wherein, β represents the equivalent elevation angle of the measured surface, β″ represents the projection angle of the equivalent elevation angle of the measured surface on the first reference plane, and θ′ represents the satellite elevation angle on the first reference plane The projection angle on the plane, γ represents the inclination angle of the measured surface, α represents the difference between the azimuth angle of the satellite and the azimuth angle of the measured surface, and α′ represents the azimuth angle of the satellite and the measured surface. The projection angle of the difference between the azimuth angles on the second reference plane, θ represents the satellite elevation angle, α _s represents the satellite azimuth angle, α _r represents the azimuth angle of the measured surface, and the first reference plane are respectively perpendicular to the measured surface and the ground, and the reflected signal receiving position is located in the first reference plane, the second reference plane is respectively perpendicular to the measured surface and the first reference plane, and The measured surface reflection position is located in the second reference plane. 6 . The interferometric deformation monitoring method according to any one of claims 2 to 5 , wherein the measured measured surface is determined according to the first phase difference and the equivalent height angle of the measured surface. 7 . Surface deformation variables, including: The deformation of the measured surface is determined by the following formula:

Wherein, d _def represents the deformation amount of the measured surface, β represents the equivalent height angle of the measured surface, ΔΦ represents the first phase difference, and λ represents the wavelength in the direct signal. 7. An interferometric deformation monitoring device, wherein the device comprises: an obtaining module, configured to obtain a first interference signal and obtain a second interference signal, wherein the first interference signal and the second interference signal are both formed by the direct signal and the reflected signal corresponding to the direct signal, the The reflected signal corresponding to the direct signal is the signal after the direct signal is reflected by the measured surface, the direct signal is the GNSS signal of the global navigation satellite system, and the emission angle of the direct signal forming the first interference signal and forming the first interference signal The emission angle of the direct signal of the second interference signal is the same, and the emission time is different; A determination module, configured to determine the deformation amount of the measured surface according to the first carrier-to-noise ratio sequence of the first interference signal and the second carrier-to-noise ratio sequence of the second interference signal. 8. The deformation monitoring device according to claim 7, wherein the determining module is specifically used for: Determine a first phase difference according to the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence, where the first phase difference is the sum of the first carrier-to-noise ratio sequence and the second carrier-to-noise ratio sequence phase difference between Determine the equivalent height angle of the measured surface according to the preset navigation message and the preset measured surface data, and the equivalent height angle is the complementary angle of the incident angle of the direct signal on the measured surface ; The deformation amount of the measured surface is determined according to the first phase difference and the equivalent height angle of the measured surface. 9. A receiver comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the computer program as claimed in the claims when executing the computer program The method of any one of 1 to 6. 10 . A computer-readable storage medium storing a computer program, wherein when the computer program is executed by a processor, the method according to any one of claims 1 to 6 is implemented. 11 .

Images