Hydrological flow reorganization and alignment method based on least square method

Abstract

The invention relates to the technical field of hydrological data reorganization, and discloses a hydrological flow reorganization and alignment method based on a least square method, comprising the following steps: S1, acquiring site water level flow data, and preprocessing the site water level flow data to generate standard site water level flow data; s2, analyzing an optimal solution of the curve based on a least square method, smoothing the curve, and controlling an error value to be minimum; s3, according to the drawn curve, carrying out three kinds of inspection according to hydrological specifications; and S4, calculating a water level and flow relation result table. According to the method, the water level flow is reorganized through the least square method, one alignment work can be completed in only one minute, the working efficiency is greatly improved, data can be processed more quickly, the alignment task can be completed more quickly, an efficient tool is provided for hydrological monitoring and analysis, and water resource management and hydraulic engineering decision making are assisted to be more timely and accurate.

Description

Technical Field

[0001] The invention relates to the technical field of hydrological data compilation, in particular to a hydrological flow compilation and routing method based on the least squares method. Background Art

[0002] Compilation of hydrological data is one of the important functions of hydrological work. The water level-flow relationship curve is an important achievement in the compilation of hydrological data. It is an important basis for hydrological forecasting, water resources management and water conservancy project design, and is of great significance for ensuring the safety and sustainable use of water resources. Hydrological alignment refers to drawing the water level-flow relationship curve and completing related calculation results. From the perspective of the scale and growth of the water conservancy industry, with the continuous growth of the demand for water resources utilization and management, the demand for hydrological alignment will continue to expand. At the same time, with the continuous advancement and development of science and technology, the technical means and level of hydrological alignment will continue to improve.

[0003] Due to the importance of the water level-flow relationship curve, the accuracy of the water level-flow relationship curve in hydrological data is very high. In the compilation of hydrological data, the traditional drawing of the water level-flow relationship curve is all done manually, which results in: first, it takes a mature worker about 4 to 6 hours to complete the work of collecting, sorting, analyzing and completing the line setting, which is inefficient; second, due to the subjectivity and experience of manual line setting, there is often a large degree of freedom, resulting in large differences in the line setting results and low curve accuracy; third, due to the lack of unified standards and specifications in the calculation and drawing process of the curve, different hydrological staff often have different drawing methods and drawing specifications, resulting in the lack of standardization of the line setting process, affecting the accuracy and availability of the drawn curve. Summary of the invention

[0004] In view of the shortcomings of the prior art, the present invention provides a hydrological flow compilation and routing method based on the least squares method, which solves the problems of low efficiency and accuracy in completing the routing work.

[0005] To achieve the above objectives, the present invention is implemented through the following technical solutions: A method for hydrologic flow alignment based on the least squares method, comprising the following steps:

[0006] S1: Obtain site water level and flow data, pre-process the site water level and flow data, and generate standard site water level and flow data;

[0007] S2: Analyze the optimal solution of the curve based on the least squares method, smooth the curve, and minimize the error value;

[0008] S3: Based on the drawn curve, three tests are carried out according to the hydrological specifications;

[0009] S4: Calculation results table of water level-flow relationship.

[0010] Preferably, the S1 comprises the following steps:

[0011] S11: sort the original data sequence S1 of the water level and flow series from small to large based on the water level to obtain a new water level and flow sequence S2;

[0012] S12: Format the water level sequence in the S2 sequence, that is, keep two decimal places for the water level; format the flow sequence in the S2 sequence, that is, according to the provisions of the "Specifications for the Compilation of Hydrological Data" (SL / T247-2020), the flow data is rounded up to the nearest even number, and three significant digits are retained. The formatted S2 sequence is taken as the new sequence S3;

[0013] S3=[(x1,y1),(x2,y2),(x3,y3)……(x n ,y n )]

[0014] Where X=(x1,x2,x3……x n ) is the flow series, Y=(y1,y2,y3……y n ) is the water level series;

[0015] The "rounding up to the nearest even number" method means:

[0016] (1) When retaining n significant digits, if the n+1th digit is less than or equal to 4, it is discarded;

[0017] (2) When n significant digits are retained, if the n+1th digit is greater than or equal to 6, the nth digit is incremented by 1;

[0018] (3) When n significant digits are retained, if the n+1th digit is equal to 5 and the following digits are 0, if the nth digit is even, discard the following digits; if the nth digit is odd, add 1; if the n+1th digit is equal to 5 and there is any non-zero digit following it, add 1 regardless of whether the nth digit is odd or even.

[0019] Preferably, S2 comprises the following steps:

[0020] S21, according to the flow series in S3 sequence (x1, x2, x3...x n ) Create a polynomial P m (x) = a0 + a1x + a2x 2 +……+a n x n

[0021] Among them, a0, a1, a2...an are polynomial coefficients, a0, a1, a2…a n The value is a random value;

[0022] S22, according to the S3 sequence, use the polynomial to fit n pairs of data (x i ,y i )(i=1,2,3……n) corresponds to it, then the S3 sequence corresponds to node x i The residual value at

[0023] R i =P m (x i )-y i (i=1,2,3……n)

[0024] All residual values ​​R i With S3 sequence node x i Correspondingly, a residual sequence R (R1, R2, R3...R n ), the sum of squares of the residual sequence is:

[0025]

[0026] S23, according to the residual sequence R calculated in step S22, adjust the polynomial coefficients a0, a1, a2...a n is b0, b1, b2...b n , using b0, b1, b2...b n Create a new polynomial and repeat step S22. Let a0, a1, a2...a n The sum of squares of the residual sequence of the polynomial formed by the coefficients is φ a , b0, b1, b2...b n The sum of squares of the residual sequence of the polynomial formed by the coefficients is φ b , let the sum of squares of the two residual sequences be Δd = φ b -φ a , when Δd≤0.001, it means that the residual sequence sum of squares has converged, and the final adjusted polynomial is the required polynomial:

[0027] P m (x) = c0 + c1x + c2x 2 +……+c n x n

[0028] Otherwise, it indicates that the sum of squares of the residual sequence has not converged, and step S23 is repeated until a polynomial that satisfies the condition of Δd≤0.001 appears;

[0029] S24, according to the polynomial obtained in step S23 and the lowest water level and the highest water level in the S3 sequence, calculate the 0.01 meter water level difference to obtain the fitting sequence:

[0030] L1=[(q1,z1),(q2,z2),(q3,z3)……(q n ,z n )]

[0031] Where, Q=(q1,q2,q3……q n ) is the flow series, Z=(z1,z2,z3……z n ) is the water level series;

[0032] S25, suppose the difference between the flow rates is: Δt = q i+1 -q i (i=1,2,3……n), the level difference between q1 and q2 is: Δt 1-2 =q2-q1, the level difference between q2 and q3 is: Δt 2-3 =q3-q2, and so on, we get the flow difference sequence: t=(Δt 1-2 ,Δt 2-3 ,Δt 3-4 ...Δt n-(n-1) );

[0033] S26, according to the flow difference sequence t, adjust the flow series in L1, smooth the curve and perform anti-curve control, specifically including the following steps:

[0034] S261, create a differential accumulation, k=0;

[0035] S262, create a new flow sequence W, and use q1 and q2 in the L1 flow series as the first two items of the W series, that is, W = (q1, q2);

[0036] S263. Starting from the third item, determine the size of the difference between two adjacent levels in the t sequence, which can be divided into the following situations:

[0037] ①When Δt n-(n-1) <Δt (n-1)-(n-2) When the nth term in the W series is W n =W n-1 +Δt (n-1)-(n-2) , and at the same time, the Δt in the t sequence n-(n-1) Replaced by Δt (n-1)-(n-2) , and Δ=Δt 9n-1)-(n-2) -Δt n-(n-1) Add to k, that is, k = k + Δ;

[0038] ②When Δt n-(n-1) >Δt (n-1)-(n-2)When k = 0, the nth term in the W series is W n =W n-1 +Δt n-(n-1) ;

[0039] ③When Δt n-(n-1) >Δt (n-1)-(n-2) When k>0, the nth term in the W series is W n =W n-1 +Δt (n-1)-(n-2) , and at the same time, the Δt in the t sequence n-(n-1) Replaced by Δt (n-1)-(n-2) ; and Δ=Δt n-(n-1) -Δt (n-1)-(n-2) Subtract it from k, that is, k = k-Δ;

[0040] ④When Δt n-(n-1) =Δt (n-1)-(n-2) When W n =W n-1 +Δt n-(n-1) ;

[0041] When calculating W3, when Δt 1-2 >Δt 2-3 , then W3=q2+Δt 1-2 , Δ=Δt 1-2 -Δt 2-3 , k=k+Δ, update t sequence to (Δt 1-2 ,Δt 1-2 ,Δt 3-4 ...Δt n-(n-1) );

[0042] When Δt 1-2 <Δt 2-3 , and k = 0, then W3 = q2 + Δt 2-3 , the t sequence remains unchanged as (Δt 1-2 ,Δt 2-3 ,Δt 3-4 ...Δt n-(n-1) );

[0043] When Δt 1-2 <Δt 2-3 , and k>0, then W3=q2+Δt 1-2 , Δ=Δt 2-3 -Δt 1-2 , k = k-Δ, update the t sequence to (Δt 1-2 ,Δt 1-2 ,Δt 3-4 ...Δt n-(n-1) );

[0044] When Δt 1-2 =Δt 2-3, W3=q2+Δt 2-3 , the t sequence does not change to:

[0045] (Δt 1-2 ,Δt 2-3 ,Δt 3-4 ...Δt n-(n-1) );

[0046] Through the above control and adjustment method of k value, the curve is realized as a one-way increasing curve, while the sum of squares of residuals is minimized. According to Q = (q1, q2, q3 ... q n ) The number of n items in the traffic sequence is calculated as W = (w1, w2, w3…w n );

[0047] S264, match the W series with the Z series to obtain the L2 sequence:

[0048] L2=[(w1,z1),(w2,z2),(w3,z3)……(w n ,z n )].

[0049] Preferably, S3 comprises the following steps:

[0050] S31. According to the water level sequence in the S3 sequence, find the corresponding water level value in L2 and form the H sequence: H = [(e1, r1), (e2, r2), (e3, r3) ... (e n ,r n )], where R = (r1, r2, r3 ... r n ) is the water level sequence, which is the same as Y=(y1,y2,y3……y n ) The water level sequence corresponds one to one, that is, R = Y; E = (e1, e2, e3...e n ) is the flow sequence;

[0051] S32 and H sequences are curve values, and S3 sequence is measured value. According to the provisions of the "Specifications for the Compilation of Hydrological Data" (SL / T247-2020), the sign test calculation of H sequence is performed:

[0052]

[0053] Among them, μ is the statistic, n is the total number of measurement points, k is the number of positive or negative signs, p and q are the probabilities of positive and negative signs, each is 0.5, and * is the number of continuous corrections;

[0054] S33. According to the provisions of the "Specifications for the Compilation of Hydrological Data" (SL / T 247-2020), the H series is subjected to a line fit test calculation:

[0055]

[0056] Wherein, μ is the statistic, n is the total number of measurement points, k is the number of times the sign is changed, a test is performed when k<0.5(n-1), otherwise no such test is performed, p and q are the probabilities of changing and not changing the sign, respectively, and are 0.5 respectively;

[0057] S34. According to the provisions of the "Specifications for the Compilation of Hydrological Data" (SL / T 247-2020), the deviation value test calculation of the H series is performed:

[0058]

[0059] Among them, t is the statistic, is the average relative deviation, for The standard deviation of n is the total number of measurement points, p i It is the relative deviation between the measuring point and the relationship curve.

[0060] Preferably, in said S4, a standard water level flow result table is prepared according to the L2 sequence.

[0061] The present invention provides a method for hydrologic flow alignment based on the least squares method, which has the following beneficial effects:

[0062] 1. The present invention uses the least squares method to compile the water level flow of the route and complete a route determination task in just one minute, which greatly improves work efficiency and can process data and complete route determination tasks more quickly. It provides an efficient tool for hydrological monitoring and analysis, and helps water resources management and water conservancy project decision-making to achieve more timely and accurate results.

[0063] 2. The solution obtained by the least squares method used in the present invention is not the optimal solution of the polynomial. The optimal solution of the polynomial is infinitely approximated by the trial calculation method. When the residual square sum converges, the optimal solution of the curve can be obtained. This method can ensure that the curve is not excessively curved, and has a good fitting effect on the measured water level flow data points, greatly improving the accuracy of the water level flow relationship line.

[0064] 3. The present invention adopts the "Specifications for the Compilation of Hydrological Data" (SL / T 247-2020) standards to carry out technical operations. The route determination results meet the requirements of the "Specifications for the Compilation of Hydrological Data" (SL / T 247-2020), and are commonly used in hydrological-related industries, greatly improving the availability of hydrological data. BRIEF DESCRIPTION OF THE DRAWINGS

[0065] Figure 1 A flow chart of a method for hydrologic flow alignment based on the least squares method of the present invention;

[0066] Figure 2This is the water level-flow relationship line of the Chenming Station of the Tangwang River in 2021, which is a hydrological flow compilation and routing method based on the least squares method of the present invention. The red line is the least squares method routing, and the purple line is a schematic diagram of manual routing. DETAILED DESCRIPTION

[0067] The technical solution of the present invention will be described clearly and completely below in conjunction with the accompanying drawings of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. 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.

[0068] Please refer to the attached Figure 1 -Attached Figure 2 The embodiment of the present invention provides a hydrological flow compilation and routing method based on the least squares method, which takes the measured flow of the Chenming Station of the Tangwang River in 2021 as a case for specific implementation, uses Python3.8 and PyQt5 as the platform basis of the technical solution, uses NumPy and Pandas as the data processing tools of the technical solution, and uses SciPy as the function construction tool of the technical solution, including the following steps:

[0069] S101: The original data sequence S1 of the water level and flow series of Chenming Station in 2021 is sorted from small to large based on water level to obtain a new water level and flow sequence S2

[0070] S102: Format the water level sequence in the S2 sequence, i.e., keep two decimal places for the water level; format the flow sequence in the S2 sequence, i.e., round the flow data to the nearest even number according to the "Specifications for the Compilation of Hydrological Data" (SL / T247-2020), and keep three significant figures. The formatted S2 sequence is taken as the new sequence S3.

[0071] S201: According to the flow series (x1, x2, x3...x n ) creates a polynomial,

[0072] P m (x) = a0 + a1x + a2x 2 +……+a n x n

[0073] S202: According to the S3 sequence, use the polynomial to fit the data and calculate the residual value at the node and the sum of squares of the residual sequence φ a

[0074] S203: Adjust the polynomial coefficients according to the residual values, refit the data, calculate the residual values ​​at the nodes, and the sum of squares of the residual sequence φ b , and then determine Δd=φ b -φ a Whether it has converged. If it has converged, the last adjusted polynomial is the required polynomial. If it has not converged, repeat the steps until a polynomial that satisfies the condition of Δd≤0.001 appears.

[0075] S204: According to the polynomial obtained in step S203 and the lowest and highest water levels in the S3 sequence, the water level difference of 0.01 meters is calculated to obtain the fitting sequence:

[0076] L1=[(q1,z1),(q2,z2),(q3,z3)……(q n ,z n )].

[0077] S205: Calculate the flow difference sequence t

[0078] S206: According to the flow extreme difference sequence t, adjust the flow series in L1, smooth the curve and perform anti-curve control. Get the new flow sequence W

[0079] S207: Match the W series with the Z series in L1 to obtain the L2 sequence:

[0080] L2=[(w1,z1),(w2,z2),(w3,z3)……(w n ,z n )].

[0081] S301: According to the water level sequence in the S3 sequence, find the corresponding water level value in L2 and form the H sequence H = [(e1, r1), (e2, r2), (e3, r3) ... (e n ,r n )]

[0082] S302: According to the provisions of the "Specifications for the Compilation of Hydrological Data" (SL / T 247-2020), the sign check calculation of the H sequence is performed

[0083]

[0084] S303: According to the provisions of the "Specifications for the Compilation of Hydrological Data" (SL / T 247-2020), the H series is tested for line fit calculation

[0085]

[0086] S304: According to the provisions of the "Specifications for the Compilation of Hydrological Data" (SL / T 247-2020), the deviation value test calculation of the H series is performed

[0087]

[0088] S04: Prepare a standard water level and flow results table based on the L2 sequence.

[0089] Table 1 shows the measured water level and flow at the Chenming Station of Tangwang River:

[0090]

[0091]

[0092] Table 2 shows the three test tables calculated at the Tangwang River Chenming Station in 2021:

[0093]

[0094] Table 3 is the water level-flow relationship table calculated at the Chenming Station of the Tangwang River in 2021:

[0095]

[0096] Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the appended claims and their equivalents.

Claims

1. A method for hydrologic flow alignment based on the least squares method, characterized in that: The following steps are involved: S1: Obtain site water level and flow data, pre-process the site water level and flow data, and generate standard site water level and flow data; S2: Analyze the optimal solution of the curve based on the least squares method, smooth the curve, and minimize the error value; S3: Based on the drawn curve, three tests are carried out according to the hydrological specifications; S4: Calculation results table of water level-flow relationship.

2. The method for hydrologic flow alignment based on least squares method according to claim 1 is characterized in that: The S1 comprises the following steps: S11: sort the original data sequence S1 of the water level and flow series from small to large based on the water level to obtain a new water level and flow sequence S2; S12: Format the water level sequence in the S2 sequence, i.e., keep two decimal places for the water level; format the flow sequence in the S2 sequence, i.e., according to the provisions of the "Specifications for the Compilation of Hydrological Data" (SL / T247-2020), the flow data is rounded up to the nearest even number, and three significant digits are retained. The formatted S2 sequence is taken as the new sequence S3; <h2 style=";text-align:left;direction:ltr">S3=[(x1,y1),(x2,y2),(x3,y3)……(x<h2 style=";text-align:left;direction:ltr"> n <h2 style=";text-align:left;direction:ltr"> ,y<h2 style=";text-align:left;direction:ltr"> n <h2 style=";text-align:left;direction:ltr"> )] Where X=(x1,x2,x3……x n ) is the flow series, Y=(y1,y2,y3……y n ) is the water level series; The "rounding up to even numbers" method means: (1) When retaining n significant digits, if the n+1th digit is less than or equal to 4, it is discarded; (2) When n significant digits are retained, if the n+1th digit is greater than or equal to 6, the nth digit is incremented by 1; (3) When n significant digits are retained, if the n+1th digit is equal to 5 and the following digits are 0, if the nth digit is even, discard the following digits; if the nth digit is odd, add 1; if the n+1th digit is equal to 5 and there is any non-zero digit following it, add 1 regardless of whether the nth digit is odd or even.

3. The method for hydrologic flow alignment based on least square method according to claim 1 is characterized in that: The S2 comprises the following steps: S21, according to the flow series in S3 sequence (x1, x2, x3...x n ) Create a polynomial P m (x) = a0 + a1x + a2x 2 +……+a n x n Among them, a0, a1, a2…a n are polynomial coefficients, a0, a1, a2…a n The value is a random value; S22, according to the S3 sequence, use the polynomial to fit n pairs of data (x i ,y i )(i=1,2,3……n) corresponds to it, then the S3 sequence corresponds to node x i The residual value at is R i =P m (x i )-y i (i=1,2,3……n) All residual values ​​R i With S3 sequence node x i Correspondingly, a residual sequence R (R1, R2, R3...R n ), the sum of squares of the residual sequence is: S23, according to the residual sequence R calculated in step S22, adjust the polynomial coefficients a0, a1, a2...a n is b0, b1, b2...b n , using b0, b1, b2...b n Create a new polynomial and repeat step S22. Let a0, a1, a2...a n The sum of squares of the residual sequence of the polynomial formed by the coefficients is φ a , b0, b1, b2...b n The sum of squares of the residual sequence of the polynomial formed by the coefficients is φ b , let the sum of squares of the two residual sequences be Δd = φ b -φ a , when Δd≤0.001, it means that the residual sequence sum of squares has converged, and the final adjusted polynomial is the required polynomial: P m (x)=c0+c1x+c2x 2 +……+c n x n Otherwise, it indicates that the sum of squares of the residual sequence has not converged, and step S23 is repeated until a polynomial that satisfies the condition of Δd≤0.001 appears; S24, according to the polynomial obtained in step S23 and the lowest water level and the highest water level in the S3 sequence, calculate the 0.01 meter water level difference to obtain the fitting sequence: <h2 style=";text-align:left;direction:ltr">L1 = [(q1,z1),(q2,z2),(q3,z3)……(q<h2 style=";text-align:left;direction:ltr"> n <h2 style=";text-align:left;direction:ltr"> ,z<h2 style=";text-align:left;direction:ltr"> n <h2 style=";text-align:left;direction:ltr"> )] Where, Q=(q1,q2,q3……q n ) is the flow series, Z=(z1,z2,z3……z n ) is the water level series; S25, suppose the difference between the flow rates is: Δt = q i+1 -q i (i=1,2,3……n), the level difference between q1 and q2 is: Δt 1-2 =q2-q1, the level difference between q2 and q3 is: Δt 2-3 =q3-q2, and so on, we get the flow difference sequence: t=(Δt 1-2 ,Δt 2-3 ,Δt 3-4 ...Δt n-(n-1) ); S26, according to the flow difference sequence t, adjust the flow series in L1, smooth the curve and perform anti-curve control, specifically including the following steps: S261, create a differential accumulation, k=0; S262, create a new flow sequence W, and use q1 and q2 in the L1 flow series as the first two items of the W series, that is, W = (q1, q2); S263. Starting from the third item, determine the size of the difference between two adjacent levels in the t sequence, which can be divided into the following situations: ①When Δt n-(n-1) <Δt (n-1)-(n-2) When the nth term in the W series is W n =W n-1 +Δt (n-1)-(n-2) , and at the same time, the Δt in the t sequence n-(n-1) Replaced by Δt (n-1)-(n-2) , and Δ=Δt (n-1)-(n-2) -Δt n-(n-1) Add to k, that is, k = k + Δ; ②When Δt n-(n-1) >Δt (n-1)-(n-2) When k = 0, the nth term in the W series is W n =W n-1 +Δt n-(n-1) ; ③When Δt n-(n-1) >Δt (n-1)-(n-2) When k>0, the nth term in the W series is W n =W n-1 +Δt (n-1)-(n-2) , and at the same time, the Δt in the t sequence n-(n-1) Replaced by Δt (n-1)-(n-2) ; and Δ=Δt n-(n-1) -Δt (n-1)-(n-2) Subtract it from k, that is, k = k-Δ; ④When Δt n-(n-1) =Δt (n-1)-(n-2) When W n =W n-1 +Δt n-(n-1) ; When calculating W3, when Δt 1-2 >Δt 2-3 , then W3=q2+Δt 1-2 , Δ=Δt 1-2 -Δt 2-3 , k=k+Δ, update t sequence to (Δt 1-2 ,Δt 1-2 ,Δt 3-4 ...Δt n-(n-1) ); When Δt 1-2 <Δt 2-3 , and k = 0, then W3 = q2 + Δt 2-3 , the t sequence remains unchanged as (Δt 1-2 ,Δt 2-3 ,Δt 3-4 ...Δt n-(n-1) ); When Δt 1-2 <Δt 2-3 , and k>0, then W3=q2+Δt 1-2 , Δ=Δt 2-3 -Δt 1-2 , k = k-Δ, update the t sequence to (Δt 1-2 ,Δt 1-2 ,Δt 3-4 ...Δt n-(n-1) ); When Δt 1-2 =Δt 2-3 , W3=q2+Δt 2-3 , the t sequence does not change to: (Δt 1-2 ,Δt 2-3 ,Δt 3-4 ……Δt n-(n-1) ); Through the above control and adjustment method of k value, the curve is realized as a one-way increasing curve, while the sum of squares of residuals is minimized. According to Q = (q1, q2, q3 ... q n ) The number of n items in the traffic sequence is calculated as W = (w1, w2, w3…w n ); S264, match the W series with the Z series to obtain the L2 sequence: <h2 style=";text-align:left;direction:ltr">L2 = [(w1,z1),(w2,z2),(w3,z3)……(w<h2 style=";text-align:left;direction:ltr"> n <h2 style=";text-align:left;direction:ltr"> ,z<h2 style=";text-align:left;direction:ltr"> n <h2 style=";text-align:left;direction:ltr"> )]。 4. The method for hydrologic flow alignment based on least squares method according to claim 1 is characterized in that: The S3 comprises the following steps: S31. According to the water level sequence in the S3 sequence, find the corresponding water level value in L2 and form the H sequence: H = [(e1, r1), (e2, r2), (e3, r3) ... (e n ,r n )], where R = (r1, r2, r3 ... r n ) is the water level sequence, which is the same as Y=(y1,y2,y3……y n ) The water level sequence corresponds one to one, that is, R = Y; E = (e1, e2, e3...e n ) is the flow sequence; S32 and H sequences are curve values, and S3 sequence is measured value. According to the provisions of the "Specifications for the Compilation of Hydrological Data" (SL / T247-2020), the sign test calculation of H sequence is performed: Among them, μ is the statistic, n is the total number of measurement points, k is the number of positive or negative signs, p and q are the probabilities of positive and negative signs, each is 0.5, and * is the number of continuous corrections; S33. According to the provisions of the "Specifications for the Compilation of Hydrological Data" (SL / T 247-2020), the H series is subjected to a line fit test calculation: Wherein, μ is the statistic, n is the total number of measurement points, k is the number of times the sign is changed, a test is performed when k<0.5(n-1), otherwise no such test is performed, p and q are the probabilities of changing and not changing the sign, respectively, and are 0.5 respectively; S34. According to the provisions of the "Specifications for the Compilation of Hydrological Data" (SL / T 247-2020), the deviation value test calculation of the H series is performed: Among them, t is the statistic, is the average relative deviation, for The standard deviation of n is the total number of measurement points, p i It is the relative deviation between the measuring point and the relationship curve.

5. The method for hydrologic flow alignment based on least square method according to claim 1 is characterized in that: In the above S4, a standard water level flow result table is prepared according to the L2 sequence.

Images