High-precision decontamination judgment method and system

Abstract

The invention relates to a high-precision sewage disposal judgment method and system, and the method comprises the steps: obtaining an average thickness value in a plurality of cross sections of a whole river channel, and carrying out the processing of the obtained average thickness value; performing curve fitting by taking the processed average thickness value as a vertical coordinate and the spacing distance of each cross section as a horizontal coordinate to obtain a sludge average thickness fitting curve; selecting a length value x according to the size of a required sludge judgment region, calculating an area value S under a fitting curve at a certain distance x by utilizing calculus, further obtaining the sludge volume of the length section x, and correcting V according to the position of the region needing to be judged; comparing the V value with preset alarm threshold values V low, V high and V risk; and selecting a desilting device and a desilting method corresponding to the V value according to the alarm grade. According to the method, the sludge volume is adopted to represent the sludge amount, judgment is more accurate, and the corrected volume value can better reflect the actual sludge amount; the labor cost is saved, and the judgment precision of the sludge amount is greatly improved.

Description

Technical Field

[0001] The present invention belongs to the technical field of river silt removal, and particularly relates to a high-precision sewage cleaning determination method and system. Background Technique

[0002] In practical engineering, to prevent silt from blocking the river channel and affecting the water flow velocity in the river channel, it is necessary to regularly clean the silt in the river channel. The existing silt cleaning technology first requires manual regular detection of the silt volume in the river channel. If cleaning is required, the cleaning program is started. It mostly relies on manual determination, with high labor costs, and the determination accuracy is subject to the experience of the operator, resulting in low accuracy.

[0003] The existing method for detecting the underwater silt thickness value usually collects the thickness value at a single point and determines whether cleaning is required based on the one-dimensional silt thickness value. This determination method has a lot of noise. For example, if the silt accumulates too high at a certain point due to a stone or a branch, but there is no silt accumulation before and after that point. If the detection point is set here, the detected data will inevitably affect the judgment and cause misjudgment of sewage cleaning. Summary of the Invention

[0004] The purpose of the present invention is to propose a high-precision sewage cleaning determination method and system in view of the defects and deficiencies of the existing technology, automatically detect and determine the sediment deposition situation in the river channel, achieve high-precision determination of the sediment deposition volume, and can perform high-precision determination of the silt volume at key areas regularly, set the alarm level, set the sewage cleaning measures based on the alarm level, set a timing module to regularly determine the silt volume in the relevant area, and utilize the automatic collection and automatic determination technology to save labor costs and greatly improve the determination accuracy of the silt volume.

[0005] To achieve the purpose of this invention, the following technical solutions are adopted in the present invention: A high-precision sewage cleaning determination method specifically includes the following steps: (1) Set a number of cross-sections si at intervals of n meters along the length direction of the river channel to be measured, where i is a natural number greater than 1. The cross-section si+1 is spaced n meters from the cross-sections si and si+2 respectively in the length direction of the river channel, and n < D / 5, where D is the total length value of the river channel to be measured; (2) Uniformly set a number of silt thickness measurement units in each cross-section si of the river channel, and obtain a number of silt thickness values a1, a2... ai according to the number of silt thickness measurement units, where i is a natural number greater than 1; (3) Calculate the average thickness avgi in the cross-section si according to the thickness values a1, a2... ai, avgi = (a1 + a2 +... + ai) / i; (4) Repeat steps (2) and (3) to obtain the average thickness values avg1, avg2... avgi in a number of cross-sections on the entire river channel; (5) Process the above average thickness values obtained in step (4). The specific processing method is as follows: If a certain average thickness value is 5 times or more of the average of its previous and next average thickness values, then delete this value; if the average thickness value of the cross-section is missing or for the above deleted average thickness values, the Kriging interpolation algorithm is used for interpolation and the obtained interpolation is corrected; (6) Use the processed average thickness values in step (5) as the ordinate and the interval distance of each cross-section as the abscissa for curve fitting to obtain the fitting curve of the average silt thickness; (7) Based on the above fitting curve of the average silt thickness, according to the required size of the silt determination area, select the length value x, and use calculus to calculate the area value S under the fitting curve within a certain distance x. , where x is a length value containing several n, x ≤ D, and avg is the average thickness value contained in the river channel of length x; (8) Determine the silt volume of this length segment x, V = S × d, where d is the average width of the river channel; (9) Correct V according to the position of the area to be determined; (10) Determine the alarm level, and compare the V value with the preset alarm thresholds V_low, V_high, and V_danger: If V ≤ V_low, it is the safe level and no cleaning is required; if V_low < V ≤ V_high, it is the general level and cleaning is required but not urgently; if V_high < V, it is the emergency level and immediate cleaning is required; (11) According to the alarm level, select the appropriate silt cleaning device and cleaning method corresponding to the V value; (12)定时循环执行步骤(3)至(10)。 (12) Timingly and repeatedly execute steps (3) to (10).

[0006] As an optimal solution, in step (2), the silt thickness measurement unit uses ultrasonic waves to measure the silt thickness, which includes a switching power supply, a main processor, a communication circuit, a boost circuit, an underwater transducer, a signal clamping circuit, a band-pass filter circuit, a logarithmic detection circuit, an amplification circuit, and a hardware watchdog. The switching power supply, the communication circuit, the boost circuit, the amplification circuit, and the hardware watchdog are all connected to the main processor. The boost circuit is connected to the signal clamping circuit through the underwater transducer, the signal clamping circuit is connected to the logarithmic detection circuit through the band-pass filter circuit, and the logarithmic detection circuit is connected to the amplification circuit; the switching power supply is used to provide DC power, the underwater transducer is a sensor for ultrasonic wave emission and echo, the signal clamping circuit is used to clamp the signal below 0.9v, the band-pass filter circuit is used to perform band-pass filtering on the echo signal, and the logarithmic detection circuit is used to perform envelope detection on the echo signal.

[0007] As an optimal solution, in step (5), the interpolation using the Kriging interpolation algorithm and the correction of the obtained interpolation are as follows: 5.1 Set the inland river section to be tested as A, and set the measured value of the river layer to be constructed as Z(avg), where {Z(avg)∈A}, Z(avg) is a second-order stable function, and the values ​​of Z(avg) in space are set as Z(avg1), Z(avg2), ..., Z(avg n ), where avg represents the spatial position of the measurement point; according to the principle of the ordinary kriging algorithm, the Z(avg0) estimate of the missing or deleted point avg0 is the weighted sum of multiple known measurement points: Z(avg0) , where Z(avg i )(i=0,1,2,...,n) is the average thickness of the silt in the cross section si avg i The function value of avg0 is the preliminary interpolation value obtained, and the other avg i is a known point; i (i=1,2,...,n) is the weight, which is determined by the value calculated by the variation function; The ordinary Kriging equations are:

[0008] According to the minimum variance estimated by ordinary Kriging, it is transformed into: =C - +u; in, represents covariance, C represents covariance function, E is the mathematical expectation value, u is the Lagrange multiplier; 5.2 Obtain the historical average thickness value vector corresponding to each moment of the cross section s0 , , , …, , t is an integer greater than 1; Using preliminary interpolation The correlation between the above historical average thickness value vectors is used to correct the above preliminary interpolation data, specifically: , is the feature vector after processing by the convolution model, (i∈[1,t]), , ,Pick ,and Close in size , (k∈[1,t]), The corresponding historical average thickness value vector is , calculate the preliminary interpolation To vector The Euclidean distance of each element in, vector The element corresponding to the jth row and the pth column in The Euclidean distance between is: , j and p are 0, 1, 2, ..., n, and the minimum Euclidean distance is obtained , the minimum Euclidean distance The corresponding element value The final interpolation data is ; As a preferred solution, the curve fitting method in step (6) adopts the spline method to perform curve fitting; As a preferred solution, step (7) targets a key area in the river channel, selects a length value xi corresponding to the length of the area at the horizontal coordinate position corresponding to the area, and determines the S value in the key area at a regular interval.

[0009] As a preferred solution, in step (9), the specific correction method is: if the measurement area is upstream of the gate, the corrected volume is V, the measurement area is downstream of the gate, then the corrected volume is V, the measurement area setting includes the gate, then the corrected volume is V1+ V2, V1 is the sludge volume corresponding to the part located upstream of the gate, and V2 is the sludge volume corresponding to the part located downstream of the gate.

[0010] As a more preferred solution, the correction parameters in step (9) are and The calculation formula is as follows: (1) The part located upstream of the gate The calculation formula is: When the cross-sections s1, s2, ..., si correspond to an upward trend in the bottom water: ; When the cross-section s1, s2, ..., si corresponds to a downward trend in the bottom water: ; in, , is a correction factor related to the trend of the water bottom in the cross section, and < , h is the average water level height corresponding to the cross sections s1, s2, ..., si, λ i is the weight corresponding to each cross section, i is a natural number greater than 1, T is the temperature data, is the water flow velocity corresponding to the cross section si; (2) Portion located downstream of the gate The calculation formula is: When the cross-sections s1, s2, ..., si correspond to an upward trend in the bottom water: ; The bottom trend corresponding to the cross sections s1, s2, ..., si is a downward trend ; in, , is a correction factor related to the trend of the water bottom in the cross section, and < , h is the average water level height corresponding to the cross sections s1, s2, ..., si, λ i is the weight corresponding to each cross section, i is a natural number greater than 1, T is the temperature data, is the water flow velocity corresponding to the cross section si.

[0011] A determination system for a high-precision pollution removal determination method, the determination system comprising: A parameter monitoring module is used to measure the silt thickness value, wherein a plurality of silt thickness measuring units are arranged at the same intervals on the cross section in the length direction of the river channel, and the silt thickness values ​​in the plurality of cross sections are measured; The first calculation module is used to obtain the average silt thickness value, calculate the average value of the silt thickness values ​​in the same cross section, and repeat the calculation to obtain the average silt thickness values ​​in multiple cross sections; The data processing module is used to process the obtained average silt thickness values. If a certain average thickness value is 5 times or more than the average of the average thickness values ​​before and after it, the value is deleted; if the average thickness value avg0 of a certain cross section is missing or for the above-mentioned deleted average thickness value, the Kriging interpolation algorithm is used to interpolate and correct the obtained interpolation; A curve fitting module is used to obtain a fitting curve of average sludge thickness, and to perform curve fitting with the average thickness value after data processing as the ordinate and the interval distance of each cross section as the abscissa, so as to obtain a fitting curve of average sludge thickness; The second calculation module is used to obtain the silt volume. Based on the above curve and the required silt determination area, the appropriate length value x is selected and the area value S under the fitting curve within a certain distance x is calculated by calculus. , x is a length value containing several n, x≤D, avg is the average thickness value contained in the river channel of length x, and the silt volume of the length segment x is obtained: V=S×d, where d is the average width of the river channel; The correction module is used to correct V. If the measurement area is upstream of the gate, the corrected volume is V1. If the measurement area is downstream of the gate, the corrected volume is V2. If the measurement area setting includes the gate, the corrected volume is V1+ V2, (V1 is the sludge volume corresponding to the part located upstream of the gate, and V2 is the sludge volume corresponding to the part located downstream of the gate); Correction parameters and The calculation formulas are as follows: (1) The part located upstream of the gate The calculation formula is: When the bottom trends corresponding to the cross-sections s1, s2,..., si are upward trends: ; When the bottom trends corresponding to the cross-sections s1, s2,..., si are downward trends: ; Among them, , are correction factors related to the walking trends of the bottoms corresponding to the cross-sections, and < , h is the average value of the water level heights corresponding to the cross-sections s1, s2,..., si, λ i is the weight corresponding to each cross-section, and i is a natural number greater than 1, T is the temperature data, is the water flow velocity corresponding to the cross-section si; (2) The part located downstream of the gate The calculation formula is: When the bottom trends corresponding to the cross-sections s1, s2,..., si are upward trends: ; When the bottom trends corresponding to the cross-sections s1, s2,..., si are downward trends ; Among them, , are correction factors related to the walking trends of the bottoms corresponding to the cross-sections, and < , h is the average value of the water level heights corresponding to the cross-sections s1, s2,..., si, λ i is the weight corresponding to each cross-section, and i is a natural number greater than 1, T is the temperature data, is the water flow velocity corresponding to the cross-section si.

[0012] The alarm level determination module compares the S value with the preset alarm thresholds S low, S high, and S danger: If S ≤ S low, it is the safe level and no cleaning is required; If S low < S ≤ S high, it is the general level and cleaning is required, but not urgently; If S high < S, it is the emergency level and immediate cleaning is required; The main control module selects the appropriate dredging device and dredging method corresponding to the S value according to the alarm level; The timing module sets the time interval and repeats the dredging determination.

[0013] The present invention has the following beneficial effects: The present invention uses the volume of silt to characterize the amount of silt, which can more accurately reflect the underwater silt amount compared with the existing single-point silt thickness value, and the determination is more precise. The determination standard is upgraded from a one-dimensional measurement method to a three-dimensional silt volume to characterize, which can achieve high-precision determination of the silt amount and reduce the probability of misjudgment. Specifically, the average value of the silt thickness values on the same cross-section is taken, and the average values of the silt thicknesses of multiple cross-sections are processed, including interpolation and correction of the interpolation, so that the silt thickness value is closer to the actual thickness value. The values after data processing are curve-fitted, and the area corresponding to the fitting curve in the relevant section is calculated based on the fitted curve, and the characteristic value representing the silt volume is further obtained by using the area. Among them, the interpolation is first based on the average thickness values corresponding to the front and rear cross-sections at the same moment for preliminary interpolation, and then the correlation between the preliminary interpolation and the historical data of the cross-section is determined. Further, based on the correlation and the Euclidean distance, the preliminary interpolation is further corrected. The interpolation data obtained by the above means is considered from both horizontal and vertical aspects, which can better reflect the actual silt thickness value and has high data accuracy. After calculating the volume value representing the silt amount, correction parameters are set based on the different positions of the section to be determined. Specifically, factors such as whether the section is upstream or downstream of the gate are considered, and the flow velocity, water level, bottom trend, temperature, and the weights of the corresponding cross-sections in the section are combined, which can filter out the influence of factors such as the environment, flow velocity, water level, and water potential. The corrected volume value can better reflect the actual silt amount. In addition, an alarm level is set, and a cleaning measure is set based on the alarm level, and a timing module is set to regularly determine the silt amount in the relevant area. By using automatic collection and automatic determination, the labor cost is saved, and the determination accuracy of the silt amount is greatly improved. Detailed implementation mode

[0014] To better understand the present invention, the content of the present invention will be further clarified below in conjunction with embodiments. However, in order to enable those of ordinary skill in the art to fully understand the technical solution and beneficial effects of the present invention, the following is further described in conjunction with specific embodiments. The embodiments are only simple examples of the present invention and do not represent or limit the scope of the patent protection of the present invention. The scope of protection of the present invention is subject to the claims.

[0015] A high-precision cleaning determination method specifically includes the following steps: (1) A plurality of cross-sections si are set at intervals of n meters along the length direction of the river to be measured, where i is a natural number greater than 1. The cross-section si + 1 is spaced n meters from the cross-sections si and si + 2 in the river length direction, and n < D / 5, where D is the total length value of the river to be measured; (2) A plurality of silt thickness measurement units are evenly arranged in each cross-section si of the river, and a plurality of silt thickness values a1, a2... ai are obtained according to the plurality of silt thickness measurement units, where i is a natural number greater than 1; (3) Calculate the average thickness avg_i within the cross-section s_i based on the thickness values a1, a2,..., ai: avg_i = (a1 + a2 +... + ai) / i; (4) Repeat steps (2) and (3) to obtain the average thickness values avg1, avg2,..., avg_i for multiple cross-sections across the entire river channel; (5) Process the above average thickness values obtained in step (4). Specific processing method: When a certain average thickness value is 5 times or more the average of its adjacent average thickness values, delete this value; if the average thickness value of a cross-section is missing or for the above-deleted average thickness values, use the Kriging interpolation algorithm for interpolation and correct the obtained interpolation; (6) Perform curve fitting with the processed average thickness values in step (5) as the ordinate and the interval distance of each cross-section as the abscissa to obtain the fitted curve of the average silt thickness; (7) Based on the above-mentioned fitted curve of the average silt thickness, select the length value x according to the required size of the silt determination area, and use calculus to calculate the area value S under the fitted curve within a certain distance x interval, , where x is a length value containing several n, x ≤ D, and avg is the average thickness value avg contained in the river channel of length x; (8) Determine the silt volume of this length segment x, V = S × d, where d is the average width of the river channel; (9) Correct V according to the location of the area to be determined; (10) Determine the alarm level by comparing the V value with the preset alarm thresholds V_low, V_high, and V_critical: If V ≤ V_low, it is the safe level and no cleaning is required; if V_low < V ≤ V_high, it is the general level and cleaning is required but not urgently; if V_high < V, it is the emergency level and immediate cleaning is required; (11) Select the appropriate silt removal device and silt removal method corresponding to the V value according to the alarm level; (12)定时循环执行步骤(3)至(10)。 (12) Timely and repeatedly execute steps (3) to (10).

[0016] As a preferred solution, the sludge thickness measuring unit in step (2) uses ultrasound to measure the sludge thickness, which includes a switching power supply, a main processor, a communication circuit, a boost circuit, an underwater transducer, a signal clamping circuit, a bandpass filter circuit, a logarithmic detection circuit, an amplifier circuit, and a hardware watchdog. The switching power supply, the communication circuit, the boost circuit, the amplifier circuit, and the hardware watchdog are all connected to the main processor. The boost circuit is connected to the signal clamping circuit through the underwater transducer, the signal clamping circuit is connected to the logarithmic detection circuit through the bandpass filter circuit, and the logarithmic detection circuit is connected to the amplifier circuit; the switching power supply is used to provide a DC power supply, the underwater transducer is a sensor that transmits and echoes ultrasonic waves, the signal clamping circuit is used to clamp the signal below 0.9V, the bandpass filter circuit is used to perform bandpass filtering on the echo signal, and the logarithmic detection circuit is used to perform envelope detection on the echo signal.

[0017] As a preferred solution, in step (5), the Kriging interpolation algorithm is used for interpolation as follows: 5.1 Set the inland river section to be tested as A, and set the measured value of the river layer to be constructed as Z(avg), where {Z(avg)∈A}, Z(avg) is a second-order stable function, and the values ​​of Z(avg) in space are set as Z(avg1), Z(avg2), ..., Z(avg n ), where avg represents the spatial position of the measurement point; according to the principle of the ordinary kriging algorithm, the Z(avg0) estimate of the missing or deleted point avg0 is the weighted sum of multiple known measurement points: Z(avg0) , where Z(avg i )(i=0,1,2,...,n) is the average thickness of the silt in the cross section si avg i The function value of avg0 is the preliminary interpolation value obtained, and the other avg i is a known point; i (i=1,2,...,n) is the weight, which is determined by the value calculated by the variation function; The ordinary Kriging equations are:

[0018] According to the minimum variance estimated by ordinary Kriging, it is transformed into: =C - +u; Where C represents the covariance function, E is the mathematical expectation value, u is the Lagrange multiplier; 5.2 Obtain the historical average thickness value vector corresponding to each moment of the cross section s0 , , , …, , t is an integer greater than 1, using preliminary interpolation The correlation between the above historical average thickness value vectors is used to correct the above preliminary interpolation data, specifically: , is the feature vector after processing by the convolution model, (i∈[1,t]), , ,Pick ,and The size closest to , (k∈[1,t]), The corresponding historical average thickness value vector is , calculate the preliminary interpolation To vector The Euclidean distance of each element in, vector The element corresponding to the jth row and the pth column in The Euclidean distance between is: , j and p are 0, 1, 2, ..., n, and the minimum Euclidean distance is obtained , the minimum Euclidean distance The corresponding element value The final interpolation data is ; As a preferred solution, the curve fitting method in step (6) adopts the spline method to perform curve fitting; As a preferred solution, step (7) targets a key area in the river channel, selects a length value xi corresponding to the length of the area at the horizontal coordinate position corresponding to the area, and determines the S value in the key area at a regular interval.

[0019] As a preferred solution, in step (9), the specific correction method is: if the measurement area is upstream of the gate, the corrected volume is V1; if the measurement area is downstream of the gate, the corrected volume is V2; if the measurement area setting includes the gate, the corrected volume is V1+ V2, V1 is the sludge volume corresponding to the part located upstream of the gate, and V2 is the sludge volume corresponding to the part located downstream of the gate.

[0020] As a more preferred solution, the correction parameters in step (9) are and The calculation formula is as follows: (1) The part located upstream of the gate The calculation formula is: When the cross-sections s1, s2, ..., si correspond to an upward trend in the bottom water: ; When the cross-section s1, s2, ..., si corresponds to a downward trend in the bottom water: ; in, , is a correction factor related to the trend of the water bottom in the cross section, and < , h is the average water level height corresponding to the cross sections s1, s2, ..., si, λ i is the weight corresponding to each cross section, i is a natural number greater than 1, T is the temperature data, is the water flow velocity corresponding to the cross section si; (2) Portion located downstream of the gate The calculation formula is: When the cross-sections s1, s2, ..., si correspond to an upward trend in the bottom water: ; The bottom trend corresponding to the cross sections s1, s2, ..., si is a downward trend ; in, , is a correction factor related to the trend of the water bottom in the cross section, and < , h is the average water level height corresponding to the cross sections s1, s2, ..., si, λ i is the weight corresponding to each cross section, i is a natural number greater than 1, T is the temperature data, is the water flow velocity corresponding to the cross section si.

[0021] A determination system for a high-precision pollution removal determination method, the determination system comprising: A parameter monitoring module is used to measure the silt thickness value, wherein a plurality of silt thickness measuring units are arranged at the same intervals on the cross section in the length direction of the river channel, and the silt thickness values ​​in the plurality of cross sections are measured; The first calculation module is used to obtain the average silt thickness value, calculate the average value of the silt thickness values ​​in the same cross section, and repeat the calculation to obtain the average silt thickness values ​​in multiple cross sections; The data processing module is used to process the obtained average silt thickness values. If a certain average thickness value is 5 times or more than the average of the average thickness values ​​before and after it, the value is deleted; if the average thickness value avg0 of a certain cross section is missing or for the deleted average thickness value, the Kriging interpolation algorithm is used to interpolate and correct the obtained interpolation; A curve fitting module is used to obtain a fitting curve of average sludge thickness, and to perform curve fitting with the average thickness value after data processing as the ordinate and the interval distance of each cross section as the abscissa, so as to obtain a fitting curve of average sludge thickness; The second calculation module is used to obtain the silt volume. Based on the above curve and the required silt determination area, the appropriate length value x is selected and the area value S under the fitting curve within a certain distance x is calculated by calculus. , x is a length value containing several n, x≤D, avg is the average thickness value contained in the river channel of length x, and the silt volume of the length segment x is obtained: V=S×d, where d is the average width of the river channel; The correction module is used to correct V. If the measurement area is upstream of the gate, the corrected volume is V1. If the measurement area is downstream of the gate, the corrected volume is V2. If the measurement area setting includes the gate, the corrected volume is V1+ V2, (V1 is the sludge volume corresponding to the part located upstream of the gate, and V2 is the sludge volume corresponding to the part located downstream of the gate); Correction parameters and The calculation formula is as follows: (1) The part located upstream of the gate The calculation formula is: When the cross-sections s1, s2, ..., si correspond to an upward trend in the bottom water: ; When the cross-section s1, s2, ..., si corresponds to a downward trend in the bottom water: ; in, , is a correction factor related to the trend of the water bottom in the cross section, and < , h is the average water level height corresponding to the cross sections s1, s2, ..., si, λ i is the weight corresponding to each cross section, i is a natural number greater than 1, T is the temperature data, is the water flow velocity corresponding to the cross section si; (2) Portion located downstream of the gate The calculation formula is: When the cross-sections s1, s2, ..., si correspond to an upward trend in the bottom water: ; The bottom trend corresponding to the cross sections s1, s2, ..., si is a downward trend ; in, , is a correction factor related to the trend of the water bottom in the cross section, and < , h is the average water level height corresponding to the cross sections s1, s2, ..., si, λ iis the weight corresponding to each cross-section, where i is a natural number greater than 1, and T is temperature data. is the water flow velocity corresponding to cross-section si.

[0022] An alarm level determination module compares the S value with preset alarm thresholds S low, S high, and S critical: If S ≤ S low, it is a safe level and no cleaning is required; if S low < S ≤ S high, it is a general level and cleaning is required but not urgently; if S high < S, it is an emergency level and immediate cleaning is required. A main control module selects a suitable silt cleaning device and method corresponding to the S value according to the alarm level. A timing module sets a time interval and repeats the silt cleaning determination.

[0023] The specific implementation is as follows: Table 1: The average thickness value avgi (unit: meter) corresponding to each cross-section si and the weight corresponding to each cross-section (unit: 1), and the distance between adjacent cross-sections is 1 meter.

[0024] It can be seen from Table 1 that the thickness value of 3m corresponding to cross-section s3 is significantly greater than 5 times the average value ((0.1 + 0.7) / 2) of the thickness value of 0.1 of the previous cross-section s2 and the thickness value of 0.7 of the subsequent cross-section s4, so it is excluded. Assume Z(avg) = , then from Z(avg) = Z(0.6)×0.1 + Z(0.1)×0.3 + Z(0.7)×0.3 + Z(1.2)×0.1 + Z(2.2)×0.2 = 1.298, the preliminary interpolation of avg3 is , and further correction is made using the correlation with the historical average data of this cross-section s3, and finally the interpolation of avg3 is = 1.036; Measure the volume of silt in the area where s1 to s6 are located, where s1 - s3 is upstream of the gate and s4 - s6 is downstream of the gate. Table 2: The volume of silt V1 (unit: cubic meter) corresponding to the upstream and the volume of silt V2 (unit: cubic meter) corresponding to the downstream

[0025] Thus, the corrected volume of silt is 1.637×1.2 + 0.32×0.8 = 3.0389 , It can be seen that the volume of silt in the area where cross-sections s1 to s6 are located is greater than the emergency level threshold (V high = 3 ), so immediate silt cleaning is required. The technical solution of the present invention is explained through the above embodiments, but the present invention is not limited to the above embodiments, that is, it does not mean that the present invention must rely on the above specific embodiments to be implemented. Any improvements made by technicians in the relevant field based on the present invention, or equivalent replacement of the materials selected by the present invention, etc., fall within the scope of protection of the patent.

Claims

1. A high-precision pollution determination method, characterized in that: Specifically, it includes the following steps: (1) Set a number of cross-sections si at intervals of n meters along the length direction of the river to be measured, where i is a natural number greater than 1. The cross-section si+1 is spaced n meters from the cross-sections si and si+2 in the length direction of the river, and n < D / 5, where D is the total length value of the river to be measured; (2) Uniformly set multiple silt thickness measurement units in each cross-section si of the river, and obtain multiple silt thickness values a1, a2... ai according to the multiple silt thickness measurement units, where i is a natural number greater than 1; (3) Calculate the average thickness avg i in the cross-section si according to the thickness values a1, a2... ai, where avg i = (a1 + a2 +... + ai) / i; (4) Repeat steps (2) and (3) to obtain the average thickness values avg1, avg2... avg i in multiple cross-sections of the entire river; (5) Process the above average thickness values obtained in step (4). The specific processing method is as follows: when a certain average thickness value is 5 times or more of the average value of its previous and next average thickness values, then delete this value; if the average thickness value of the cross-section is missing or for the above deleted average thickness values, the Kriging interpolation algorithm is used for interpolation and the obtained interpolation is corrected; (6) Perform curve fitting with the average thickness values processed in step (5) as the ordinate and the interval distance of each cross-section as the abscissa to obtain the silt average thickness fitting curve; (7) Based on the above average silt thickness fitting curve, according to the required silt determination area size, select the length value x, and use calculus to calculate the area value S under the fitting curve within a certain distance x. , x is a length value containing several n, x≤D, avg is the average thickness value avg contained in the river channel of length x; (8) Determine the silt volume V of this length segment x, where V = S × d, and d is the average width of the river; (9) Correct V according to the position of the area to be determined; (10) Determine the alarm level, and compare the V value with the preset alarm thresholds V low, V high, and V danger: if V ≤ V low, it is the safe level and no cleaning is required; if V low < V ≤ V high, it is the general level and cleaning is required, but not in a hurry; if V high < V, it is the emergency level and immediate cleaning is required; (11) According to the alarm level, select a suitable silt removal device and silt removal method corresponding to the V value; (12)定时循环执行步骤(3)至(10)。定时循环执行步骤(3)至(10)。 2. The high-precision cleaning determination method according to claim 1, characterized in that: In step (2), the silt thickness measurement unit uses ultrasonic waves to measure the silt thickness.

3. The high-precision cleaning determination method according to claim 2, characterized in that: The silt thickness measurement unit in step (2) includes a switching power supply, a main processor, a communication circuit, a boost circuit, an underwater transducer, a signal clamping circuit, a band-pass filter circuit, a logarithmic detection circuit, an amplifier circuit, and a hardware watchdog. The switching power supply, the communication circuit, the boost circuit, the amplifier circuit, and the hardware watchdog are all connected to the main processor. The boost circuit is connected to the signal clamping circuit through the underwater transducer. The signal clamping circuit is connected to the logarithmic detection circuit through the band-pass filter circuit. The logarithmic detection circuit is connected to the amplifier circuit. The switching power supply is used to provide DC power. The underwater transducer is a sensor for ultrasonic wave emission and echo. The signal clamping circuit is used to clamp the signal below 0.9v. The band-pass filter circuit is used to perform band-pass filtering on the echo signal. The logarithmic detection circuit is used to perform envelope detection on the echo signal.

4. The high-precision cleaning determination method according to claim 1, characterized in that: In step (5), the specific method of using the Kriging interpolation algorithm for interpolation is as follows: Step 5.1 Set the inland river section to be tested as A, and set the measured value of the river layer to be constructed as Z(avg), where {Z(avg)∈A}, Z(avg) is a second-order stable function, and the values ​​of Z(avg) in space are set as Z(avg1), Z(avg2), ..., Z(avg n ), where avg represents the spatial position of the measurement point; according to the principle of the ordinary kriging algorithm, the Z(avg0) estimate of the missing or deleted point avg0 is the weighted sum of multiple known measurement points: Z(avg0) , where Z(avg i )(i=0,1,2,...,n) is the average thickness of the silt in the cross section si avg i The function value of avg0 is the preliminary interpolation value obtained, and the other avg i is a known point; i (i=1,2,...,n) is the weight, which is determined by the value calculated by the variation function; The ordinary Kriging equations are: ; According to the minimum variance estimated by ordinary Kriging, it is transformed into: =C - +u; in, represents covariance, C represents covariance function, E is the mathematical expectation value, u is the Lagrange multiplier; Step 5.2 Obtain the historical average thickness value vector corresponding to each moment of the cross section s0 , , , …, , t is an integer greater than 1; Using preliminary interpolation The correlation between the above historical average thickness value vectors is used to correct the above preliminary interpolation data, specifically: , is the feature vector after processing by the convolution model, (i∈[1,t]), , ,Pick ,and The size closest to , (k∈[1,t]), The corresponding historical average thickness value vector is , calculate the preliminary interpolation To vector The Euclidean distance of each element in, vector The element corresponding to the jth row and the pth column in The Euclidean distance between is: , j and p are 0, 1, 2, ..., n, and the minimum Euclidean distance is obtained , the minimum Euclidean distance The corresponding element value The final interpolation data is .

5. The high-precision cleaning determination method according to claim 1, characterized in that: In step (6), the curve fitting method uses the spline method for curve fitting.

6. The high-precision cleaning determination method according to claim 1, characterized in that: In step (7), for the key areas in the river channel, at the abscissa position corresponding to this area, a length value xi corresponding to the length of this area is selected, and the S value in this key area is determined at regular intervals.

7. The high-precision cleaning determination method according to claim 1, characterized in that: In step (9), the specific correction method is: if the measurement area is upstream of the gate, the corrected volume is V1; if the measurement area is downstream of the gate, the corrected volume is V2; if the measurement area setting includes the gate, the corrected volume is V1+ V2, V1 is the sludge volume corresponding to the part located upstream of the gate, and V2 is the sludge volume corresponding to the part located downstream of the gate. and To correct the parameters.

8. The high-precision cleaning determination method according to claim 7, characterized in that: The part located upstream of the gate The calculation formula is: When the bottom trends corresponding to the cross-sections s1, s2, …, si are upward trends: ； When the bottom trends corresponding to the cross-sections s1, s2, …, si are downward trends: ； in, , is a correction factor related to the trend of the water bottom in the cross section, and < , h is the average water level height corresponding to the cross sections s1, s2, ..., si, λ i is the weight corresponding to each cross section, i is a natural number greater than 1, T is the temperature data, is the water flow velocity corresponding to the cross section si.

9. The high-precision cleaning determination method according to claim 7, characterized in that: The part downstream of the gate The calculation formula is: When the bottom trends corresponding to the cross-sections s1, s2, …, si are upward trends: ； When the bottom trends corresponding to the cross-sections s1, s2, …, si are downward trends: ； in, , is a correction factor related to the trend of the water bottom in the cross section, and < , h is the average water level height corresponding to the cross sections s1, s2, ..., si, λ i is the weight corresponding to each cross section, i is a natural number greater than 1, T is the temperature data, is the water flow velocity corresponding to the cross section si.

10. A high-precision pollution determination system, characterized in that: The determination system includes: A parameter monitoring module, which is used to measure the silt thickness value. A plurality of silt thickness measurement units are arranged at the same interval on the cross-sections in the length direction of the river channel, and the silt thickness values in a plurality of cross-sections are measured; A first calculation module, which is used to obtain the average silt thickness value, calculate the average value of the silt thickness values in the same cross-section, and repeat the calculation to obtain the average silt thickness values in a plurality of cross-sections; A data processing module, which is used to process the obtained average silt thickness value. If a certain average thickness value is 5 times or more of the average value of its previous and subsequent average thickness values, then this value is deleted; if the average thickness value avg0 of a certain cross-section is missing or for the above-deleted average thickness values, the Kriging interpolation algorithm is used for interpolation and the obtained interpolation is corrected; A curve fitting module, which is used to obtain the silt average thickness fitting curve. Using the processed average thickness value as the ordinate and the interval distance of each cross-section as the abscissa for curve fitting, the silt average thickness fitting curve is obtained; The second calculation module is used to obtain the silt volume. Based on the above curve and the required silt determination area, the appropriate length value x is selected and the area value S under the fitting curve within a certain distance x is calculated by calculus. , x is a length value containing several n, x≤D, avg is the average thickness value contained in the river channel of length x, and the silt volume of the length segment x is obtained: V=S×d, where d is the average width of the river channel; The correction module is used to correct V. If the measurement area is upstream of the gate, the corrected volume is V1. If the measurement area is downstream of the gate, the corrected volume is V2. If the measurement area setting includes the gate, the corrected volume is V1+ V2, V1 is the sludge volume corresponding to the part located upstream of the gate, and V2 is the sludge volume corresponding to the part located downstream of the gate; and is the correction parameter; An alarm level determination module, which compares the S value with the preset alarm thresholds S low, S high, and S danger: If S ≤ S low, it is a safe level and no cleaning is required; if S low < S ≤ S high, it is a general level and cleaning is required, but not urgently; if S high < S, it is an emergency level and immediate cleaning is required; A main control module, which selects a suitable silt removal device and silt removal method corresponding to the S value according to the alarm level; A timing module, which sets a time interval and repeats the silt removal determination.

Images