A method for correcting vertical displacement results of annual regulation reservoir

Abstract

The present invention discloses a method for correcting the vertical displacement results of an annual regulating reservoir, which is applicable to the annual regulating reservoir. The existing vertical line observation does not consider the influence of reservoir water level change on the vertical line deviation and the stability of the deep anchor point. Therefore, on the basis of the original vertical line system, the present invention arranges GNSS receivers at the same position to carry out continuous synchronous observation, and uses the GNSS monitoring results to correct the vertical line.

Description

Technical Field

[0001] The invention relates to the field of dam monitoring, and in particular to a method for correcting vertical displacement results of an annual regulating reservoir. Background Art

[0002] With the in-depth development of my country's hydropower construction, many super-high dams have been built, such as the Xiaowan Dam, which is a hyperbolic arch dam with a height of more than 290 meters. Dam safety monitoring usually uses the inverted vertical system as the monitoring benchmark. The inverted vertical connection is set under the dam foundation, and the other elevations are all based on the positive vertical line to observe the horizontal displacement of the dam.

[0003] The current method for calculating the horizontal displacement of a dam by vertical observation is basically to first obtain the radial and tangential displacement increments of each observation point by radial and tangential measurements at the location; then directly algebraically superimpose the radial and tangential displacement increments obtained from observation points at different elevations to obtain the radial and tangential displacements of the corresponding measurement point. For annual regulation reservoirs, the existing method for calculating the results of vertical observations does not take into account the impact of reservoir water level changes on the vertical deviation and the stability of deep anchor points. Summary of the invention

[0004] The present invention aims to provide a method for correcting the vertical displacement results of an annual regulating reservoir, by synchronously observing the vertical line observation of a GNSS and a measuring point, and correcting the measuring point using the GNSS monitoring results.

[0005] The present invention is achieved through the following technical solutions:

[0006] A method for correcting the vertical displacement results of an annual regulating reservoir comprises the following steps:

[0007] S1: After arranging the GNSS reference points, the measurement point data set is obtained through GNSS, and the displacement results of the measurement points are obtained by performing single-base station differential positioning with the known reference station as the starting point;

[0008] S2: Evaluate the accuracy of the displacement results of the measuring points, perform weighted fusion on the displacement results that meet the requirements, and output the GNSS displacement sequence;

[0009] S3: Based on the output result of the above step S2, the difference is calculated with the vertical observation data of the measuring point, and the vertical-GNSS difference measurement sequence is output. After filtering and decomposition, the sequence is corrected. The corrected vertical correction number sequence is combined and analyzed using the reservoir water level data to determine whether its correlation meets the preset value. If it meets the preset value, the displacement deviation of the measuring point is corrected after the parameters are corrected; if it does not meet the preset value, the filtering and decomposition are continued and the sequence correction is continued.

[0010] As an optional method of the present invention, the GNSS reference point arrangement includes the following steps:

[0011] At least two GNSS reference points are arranged outside the influence range of the water level of the reservoir behind the dam, and basic observation piers are constructed at the same location where each GNSS reference point is arranged.

[0012] As an optional method of the present invention, in the above step S1, obtaining the measurement point data set through GNSS includes the following steps:

[0013] The time dimension direction is determined, and based on the time dimension direction, the data of multiple measurement points are continuously detected by GNSS to generate a detection data set, and the detection data set is preprocessed to obtain a measurement point data set.

[0014] As an optional method of the present invention, preprocessing the detection data set includes the following steps:

[0015] Calibrate the reference frame data, and identify the difference between the reference frame data and its adjacent data;

[0016] Traverse the detection data of all frames, extract all difference information, and select valid detection data that meets the preset difference requirements;

[0017] The effective detection data are merged and processed to obtain the measurement point data set.

[0018] As an optional method of the present invention, in the above step S2, the accuracy assessment includes arranging the displacement results according to the time series, eliminating errors and outputting the observation accuracy, which includes the following steps:

[0019] Select any displacement results as annotation points, connect the annotation points, fit all displacement results to output the deformation process curve, move the annotation point connection line to the deformation process curve, remove the points of displacement results whose deviation values ​​exceed the preset value, and calculate the observation accuracy of the deformation process curve.

[0020] As an optional method of the present invention, weighted fusion includes merging the displacement result and the observation accuracy to generate the displacement accuracy, and then normalizing it with the number of measuring points to establish its accuracy weight.

[0021] As an optional method of the present invention, in the above step S3, the difference is calculated by statistically analyzing the difference between the displacement results of the GNSS observation data and the measurement points over the years; the vertical correction number sequence is obtained by analyzing the difference and extracting a curve with obvious annual periodic fluctuations.

[0022] As an optional method of the present invention, the correlation is obtained by statistically analyzing the adjustment of reservoir water level changes over the years and performing correlation analysis with the vertical line correction number sequence through grey correlation and cluster analysis.

[0023] In addition, to achieve the above-mentioned purpose, this embodiment also provides an electronic device, including a memory and a processor, the memory stores a computer program that can be run on the processor, and when the processor executes the computer program, the steps of the above-mentioned method for correcting the vertical displacement results of the annual regulation reservoir are implemented.

[0024] In addition, to achieve the above purpose, this embodiment also provides a storage medium on which a computer program is stored. When the computer program is executed by a processor, the steps in the above-mentioned method for correcting the vertical displacement results of the annual regulation reservoir are implemented.

[0025] Compared with the prior art, the present invention has the following advantages and beneficial effects:

[0026] The method for correcting the vertical displacement results of the annual regulating reservoir provided by the present invention is suitable for the annual regulating reservoir. The existing vertical line observation does not take into account the influence of the reservoir water level change on the vertical line deviation and the stability of the deep anchor point. Therefore, on the basis of the original vertical line system, GNSS receivers are arranged at the same position to carry out continuous synchronous observation, and the vertical line is corrected using the GNSS monitoring results. BRIEF DESCRIPTION OF THE DRAWINGS

[0027] The drawings described herein are used to provide a further understanding of the embodiments of the present invention, constitute a part of this application, and do not constitute a limitation of the embodiments of the present invention. In the drawings:

[0028] Figure 1 A schematic flow chart of a method for correcting vertical displacement results of an annual regulating reservoir provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0029] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with embodiments and drawings. The exemplary embodiments of the present invention and their description are only used to explain the present invention and are not intended to limit the present invention.

[0030] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0031] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0032] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0033] Therefore, the following detailed description of the embodiments of the present invention provided in the embodiments is not intended to limit the scope of the invention claimed for protection, but merely represents selected embodiments of the present invention. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0034] Therefore, the detailed description of the embodiments of the present invention provided below is not intended to limit the scope of the claimed invention, but merely represents selected embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0035] Example

[0036] The change of reservoir water level reflects the change of water volume in the reservoir. For large reservoirs, the change of volume will cause the change of local gravity field (including gravity value and gravity direction). The inverted vertical system based on gravity will affect the change of vertical line direction and generate vertical line deviation error due to the change of gravity direction, which will directly affect the measured value of the vertical line. This influence will vary depending on the distance from the dam to the reservoir, the reservoir type and the depth of the vertical anchor block.

[0037] Δ=ΔuL / p

[0038] Where: Δu is the vertical line deviation angle, and its size is related to factors such as reservoir type, vertical line position, and change in reservoir water level; L is the length of the inverted vertical line.

[0039] For large reservoirs with significant changes in water level, the vertical deviation has a certain value. Taking Longyangxia as an example, when the reservoir water level drops from 2600m to 2540m, the calculated vertical deviation angle Δu=0.51". If the total length of the vertical line is 90m, it can be calculated that

[0040]

[0041] It can be seen from this that the impact of changes in the local gravity field on the accuracy of vertical line observation will also reach a certain order of magnitude and increase with the increase in the length of the vertical line.

[0042] Different types of dam bodies have different effects on the stability of the inverted anchor points due to their own weight loads, different engineering geological conditions of the dam foundation, and especially the changes in the reservoir water level during operation. Due to the large fluctuation of the reservoir water level, the high dam has a certain impact on the downstream displacement of the dam foundation layer at different depths caused by the changes in the dam foundation and reservoir basin load. If the inverted anchor block is not buried deep enough, its stability cannot be guaranteed. According to the current research results, when the reservoir water level changes from 25m to 80m, in order to make the stability of the inverted anchor block not less than 0.1mm, the anchor block should be buried at a depth of more than 70m, and the corresponding applied buoyancy is 1143MPa. When the reservoir water level changes by 60m, the error of the inverted datum may reach 0.52mm, and its stability is not ideal. Further burying the inverted anchor will greatly increase the technical difficulty and cost, and will also affect its sensitivity.

[0043] GNSS baseline solution uses relative positioning technology, which is also called differential positioning. The most basic situation is that the receivers placed at both ends of the baseline perform synchronous observations to obtain observation data of the same epoch and satellite, and determine the relative position of the two end points of the baseline through data solution processing. Baseline solution is a key step in high-precision GNSS measurement and positioning data processing. It uses the carrier phase observation data obtained by synchronous observation of the receivers at both ends of the baseline, establishes the error equation through the differential positioning principle, and then uses the least squares principle to solve it, and finally obtains the three-dimensional coordinate vector of the baseline.

[0044] GNSS deformation monitoring requires the use of a reference station with known coordinates as the starting point. When solving the GNSS baseline, the reference point is used as a known quantity, and the coordinates of the monitoring point are obtained by calculating the baseline coordinate vector difference. Taking the GPS constellation system as an example, the coordinate frame used by the GPS system is the WGS-84 coordinate system. During the relative positioning process, the result is the displacement value under the reference ellipsoid, which is not affected by the change of the gravity vertical line deviation. At the same time, the GNSS reference point can be arranged outside the range affected by the reservoir water level change to avoid being affected by the rise and fall of the water level. Therefore, this embodiment is implemented as follows:

[0045] See also Figure 1 This embodiment provides a method for correcting the vertical displacement result of an annual regulating reservoir, comprising the steps of:

[0046] S1: After arranging the GNSS reference points, the measurement point data set is obtained through GNSS, and the displacement results of the measurement points are obtained by single base station differential positioning with the known reference station as the starting point. The reference station with known coordinates is used as the starting point. When solving the GNSS baseline, the reference point is used as a known quantity, and the coordinates of the monitoring point are obtained by calculating the baseline coordinate vector difference. The GNSS reference point arrangement process is to arrange two or more GNSS reference points outside the influence range of the reservoir water level behind the dam. Since the GNSS observation accuracy and reliability are lower than the vertical line, a multi-base station solution method is used to improve the accuracy and reliability of the GNSS results. At the same time, in order to ensure the consistency of the vertical line and the GNSS results, a basic observation pier is constructed at the same location where each GNSS reference point is arranged, and it needs to be arranged in the same position for point construction.

[0047] Among them, in order to ensure the validity of the measurement point data set, it includes preprocessing after the collection is completed. In this embodiment, the time dimension direction is first determined, that is, the passage of time and the interval of collection. And based on the time dimension direction, the data of multiple measurement points are continuously detected by GNSS to generate a detection data set, and the detection data set is preprocessed to obtain the measurement point data set. That is, following the passage of time, the vertical line information of the measurement points at different time points is collected, and multiple points with different vertical line angles following the time change are generated according to the order of the passage of time to achieve the monitoring effect.

[0048] In addition, in order to ensure the effectiveness of the monitoring results and to continue to participate in subsequent optimization, the preprocessing of the detection data set includes the following steps: calibrating the reference frame data, which can be customized according to the actual use environment, and this embodiment does not limit it. The reference frame data and its adjacent data are used to identify the difference; traverse the detection data of all frames, extract all the difference information, and select the valid detection data that meets the preset difference requirements; merge the valid detection data to obtain the measurement point data set.

[0049] S2: Evaluate the accuracy of the displacement results of the measuring points, perform weighted fusion on the displacement results that meet the conditions, and output the GNSS displacement sequence. However, since the observation error is an accidental error, it obeys the normal distribution. At the same time, the deformation of the building is also a dynamic and continuous process, which is relatively stable. Therefore, the accuracy evaluation includes arranging the displacement results according to the time series, eliminating errors and outputting the observation accuracy, which includes the following steps:

[0050] Select any displacement results as annotation points, connect the annotation points, fit all displacement results to output the deformation process curve, move the annotation point connection line to the deformation process curve, remove the points of displacement results whose deviation value exceeds the preset value, and calculate the observation accuracy σ of the deformation process curve ij .

[0051] After that, weighted fusion includes merging the displacement results and observation accuracy to generate displacement accuracy, and then normalizing it with the number of measurement points to establish its accuracy weight p ij

[0052]

[0053] Where n represents the number of base stations; p ij Represents weight, the value range is (0, 1), the higher the precision, the greater the weight.

[0054] S3: Based on the output result of the above step S2, the difference is calculated with the vertical observation data of the measuring point, and the vertical-GNSS difference measurement sequence is output. After filtering and decomposition, the sequence is corrected. The corrected vertical correction number sequence is combined and analyzed using the reservoir water level data to determine whether its correlation meets the preset value. If it meets the preset value, the displacement deviation of the measuring point is corrected after the parameters are corrected; if it does not meet the preset value, the filtering and decomposition are continued and the sequence correction is continued.

[0055] In this embodiment, synchronous observation of the vertical line and GNSS displacement results refers to the statistical deformation results of the dam vertical line and GNSS horizontal displacement during the dam operation period, water storage, water level drop, and high water level operation over the years. The difference calculation is to obtain the time-dependent vertical line-GNSS difference sequence by statistically analyzing the difference between the GNSS observation data and the displacement results of the measuring points over the years; the vertical line correction number sequence is obtained by analyzing the difference and extracting a curve with obvious annual cycle fluctuations. Among them, the correlation is obtained by statistically analyzing the adjustment of the reservoir water level changes over the years, and correlating it with the vertical line correction number sequence through grey correlation and cluster analysis. Please refer to it again. Figure 1 In this embodiment, the correlation should be greater than 0.8, then the correction is started; otherwise, the above filtering decomposition and sequence correction are repeated. In the subsequent correction, finite element statistical analysis is used to take the reservoir water level as the dependent variable, calculate its influence coefficient, C, and correct the vertical line observation result.

[0056]

[0057] Where α is the influence coefficient, and ΔH is the water level fluctuation (current reservoir water level - minimum reservoir water level).

[0058] In addition, to achieve the above-mentioned purpose, this embodiment also provides an electronic device, including a memory and a processor, the memory stores a computer program that can be run on the processor, and when the processor executes the computer program, the steps of the above-mentioned method for correcting the vertical displacement results of the annual regulation reservoir are implemented.

[0059] In addition, to achieve the above purpose, this embodiment also provides a storage medium on which a computer program is stored. When the computer program is executed by a processor, the steps in the above-mentioned method for correcting the vertical displacement results of the annual regulation reservoir are implemented.

[0060] This scheme is applicable to annual regulation reservoirs. The existing vertical line observation does not take into account the impact of reservoir water level changes on the vertical line deviation and the stability of deep anchor points. Therefore, on the basis of the original vertical line system, GNSS receivers are arranged at the same location for continuous synchronous observation, and the vertical line is corrected using the GNSS monitoring results.

[0061] The specific implementation methods described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation method of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for correcting the vertical displacement results of an annual regulating reservoir, characterized in that: Includes steps: S1: After arranging the GNSS reference points, the measurement point data set is obtained through GNSS, and the displacement results of the measurement points are obtained by performing single-base station differential positioning with the known reference station as the starting point; S2: Evaluate the accuracy of the displacement results of the measuring points, perform weighted fusion on the displacement results that meet the requirements, and output the GNSS displacement sequence; In step S2, the accuracy assessment includes arranging the displacement results in time series, eliminating errors and outputting the observation accuracy, which includes the following steps: Select any displacement results as marked points, connect the marked points, fit all displacement results to output a deformation process curve, move the marked point connection line to the deformation process curve, remove the points of the displacement results whose deviation value exceeds the preset value, and calculate the observation accuracy of the deformation process curve; The weighted fusion includes merging the displacement result and the observation accuracy to generate the displacement accuracy, and then normalizing it with the number of measuring points to establish its accuracy weight; S3: According to the output result of the above step S2, the difference is calculated with the vertical observation data of the measuring point, and the vertical-GNSS difference measurement sequence is output. After filtering and decomposing, the sequence is corrected. The corrected vertical correction number sequence is combined and analyzed using the reservoir water level data to determine whether its correlation meets the preset value. If it meets the preset value, the displacement deviation of the measuring point is corrected after the parameters are corrected; if it does not meet the preset value, the filtering and decomposition are continued and the sequence correction is continued; The difference calculation is to statistically calculate the difference between the displacement results of the GNSS observation data and the measurement points over the years; the vertical correction number sequence is obtained by analyzing the difference and extracting the curve with obvious annual cycle fluctuations.

2. The method for correcting the vertical displacement results of an annual regulating reservoir according to claim 1, characterized in that: The GNSS benchmark point layout includes the following steps: At least two GNSS reference points are arranged outside the influence range of the water level of the reservoir behind the dam, and basic observation piers are constructed at the same location where each GNSS reference point is arranged.

3. The method for correcting the vertical displacement results of an annual regulating reservoir according to claim 1, characterized in that: In the above step S1, obtaining the measurement point data set through GNSS includes the following steps: The time dimension direction is determined, and based on the time dimension direction, the data of multiple measurement points are continuously detected by GNSS to generate a detection data set, and the detection data set is preprocessed to obtain a measurement point data set.

4. A method for correcting the vertical displacement results of an annual regulating reservoir according to claim 3, characterized in that: The preprocessing of the detection data set includes the following steps: Calibrate the reference frame data, and identify the difference between the reference frame data and its adjacent data; Traverse the detection data of all frames, extract all difference information, and select valid detection data that meets the preset difference requirements; The effective detection data are merged to obtain the measurement point data set.

5. The method for correcting the vertical displacement result of the annual regulating reservoir according to claim 1 is characterized in that: The correlation is obtained by statistically analyzing the adjustment of reservoir water level changes over the years and correlating it with the vertical line correction number sequence through grey correlation and cluster analysis.

6. An electronic device comprising a memory and a processor, wherein the memory stores a computer program that can be run on the processor, characterized in that: When the processor executes the computer program, the steps of the method for correcting the vertical displacement results of the annual regulating reservoir as described in any one of claims 1 to 5 are implemented.

7. A storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps in the method for correcting the vertical displacement results of the annual regulating reservoir as described in any one of claims 1 to 5 are implemented.

Images