A high-precision pollution removal determination method and system

Abstract

The present invention relates to a high-precision pollution cleaning determination method and system, which obtains the average thickness value of multiple cross-sections on the entire river channel, processes the obtained average thickness value, and performs curve fitting with the processed average thickness value as the vertical coordinate and the interval distance of each cross-section as the horizontal coordinate to obtain the average silt thickness fitting curve; selects a length value x according to the required silt determination area size, uses calculus to calculate the area value S under the fitting curve within a certain interval x, further obtains the silt volume of the length segment x, and corrects V according to the location of the area to be determined; compares the V value with preset alarm thresholds V low, V high, and V dangerous; and selects a suitable silt cleaning device and silt cleaning method corresponding to the V value according to the alarm level. The present invention uses silt volume to characterize the silt amount, making the determination more accurate, and the corrected volume value can better reflect the actual silt amount; saves labor costs, and greatly improves the accuracy of silt amount determination.

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 Art

[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 detection of the silt volume in the river channel at regular intervals. If cleaning is required, the cleaning program is started. It mostly relies on manual judgment, resulting in high labor costs. The judgment accuracy is subject to the experience of the operator and is low.

[0003] The existing detection of 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, due to stones or branches at a certain place, the silt accumulation at that point is too high, but there is no silt accumulation before and after that point. If the detection point is set here, the detected data will surely 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 for 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 the timing module, regularly determine the silt volume in the relevant area, and utilize the automatic collection and automatic determination technology, saving labor costs and greatly improving 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:

[0006] A high-precision sewage cleaning determination method specifically includes the following steps:

[0007] (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 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;

[0008] (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;

[0009] (3) Calculate the average thickness avgi in the cross-section si according to the thickness values a1, a2... ai, where avgi = (a1 + a2 +... + ai) / i;

[0010] (4) Repeat steps (2) and (3) to obtain the average thickness values avg1, avg2... avgi in multiple cross-sections of the entire river channel;

[0011] (5) Process the above average thickness values obtained in step (4). Specific processing means: when 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, use the Kriging interpolation algorithm for interpolation and correct the obtained interpolation;

[0012] (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 fitting curve of the average silt thickness;

[0013] (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 interval, , where x is the length value including several n, x ≤ D, and avg is the average thickness value avg included in the river channel with length x;

[0014] (8) Calculate the silt volume of this length segment x, V = S × d, where d is the average width of the river channel;

[0015] (9) Correct V according to the position of the area to be determined;

[0016] (10) Determine the alarm level, and compare 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;

[0017] (11) According to the alarm level, select the appropriate silt cleaning device and silt cleaning method corresponding to the V value;

[0018] (12)定时循环执行步骤(3)至(10)。 (12) Timely and circularly execute steps (3) to (10).

[0019] As a preferred solution, the silt thickness measuring unit in step (2) uses ultrasound 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 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.

[0020] As a preferred solution, in step (5), the Kriging interpolation algorithm is used to interpolate and the obtained interpolation is corrected as follows:

[0021] 5.1 Set the inland river section to be tested as A, and the river layer measurement value to be constructed as Z(avg), with {Z(avg)∈A}, Z(avg) is a second-order stationary function, and the values ​​of Z(avg) in space are set as Z(avg1), Z(avg2), ..., Z(avg n ), where avg represents the spatial location of the measurement point; according to the principle of the ordinary kriging algorithm, the estimated value of Z(avg0) 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 silt in cross section si avg i Function value, avg0 is the preliminary interpolation value obtained, and 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;

[0022] The ordinary Kriging equations are:

[0023] According to the minimum variance estimated by ordinary Kriging, the transformation is: =C - +u;

[0024] in, represents covariance, C represents covariance function, E is the mathematical expectation value, u is the Lagrange multiplier;

[0025] 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;

[0026] Using preliminary interpolation The correlation between the above-mentioned historical average thickness value vectors is used to correct the above-mentioned 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 row j and column p in The Euclidean distance between them is: , j and p are 0, 1, 2, ..., n, to obtain the minimum Euclidean distance , the minimum Euclidean distance Corresponding element value The final interpolation data is ;

[0027] As a preferred solution, the curve fitting method in step (6) adopts the spline method to perform curve fitting;

[0028] 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 regularly determines the S value in the key area.

[0029] 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, and the corrected volume is V, the measurement area setting includes the gate, then the corrected volume is V1+ V2, V1 is the silt volume corresponding to the part upstream of the gate, and V2 is the silt volume corresponding to the part downstream of the gate.

[0030] As a more preferred solution, the correction parameters in step (9) and The calculation formula is as follows:

[0031] (1) The part located upstream of the gate The calculation formula is:

[0032] When the bottom trend corresponding to the cross sections s1, s2, ..., si is an upward trend: ;

[0033] When the bottom trend of cross sections s1, s2, ..., si is a downward trend: ;

[0034] in, 、 is a correction factor related to the trend of the bottom of the cross section, and < , h is the average water level height of 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;

[0035] (2) The part located downstream of the gate The calculation formula is:

[0036] When the bottom trend corresponding to the cross sections s1, s2, ..., si is an upward trend: ;

[0037] The bottom trend corresponding to cross sections s1, s2, ..., si is a downward trend ;

[0038] in, 、 is a correction factor related to the trend of the bottom of the cross section, and < , h is the average water level height of 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.

[0039] A high-precision pollution removal determination method determination system, the determination system comprising:

[0040] A parameter monitoring module is used to measure the silt thickness value, wherein a plurality of silt thickness measuring units are arranged at equal intervals on the cross section in the longitudinal direction of the river channel, and the silt thickness values ​​in the plurality of cross sections are measured;

[0041] The first calculation module is used to obtain an 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;

[0042] 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 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.

[0043] A curve fitting module is used to obtain a fitting curve of average silt thickness, and to perform curve fitting using the average thickness value after data processing as the ordinate and the interval distance of each cross section as the abscissa to obtain a fitting curve of average silt thickness;

[0044] 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 using 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 x is obtained: V=S×d, where d is the average width of the river channel;

[0045] Correction module, 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 silt volume corresponding to the part located upstream of the gate, and V2 is the silt volume corresponding to the part located downstream of the gate);

[0046] Correction parameters and The calculation formula is as follows:

[0047] (1) The part located upstream of the gate The calculation formula is:

[0048] When the bottom trend corresponding to the cross sections s1, s2, ..., si is an upward trend: ;

[0049] When the bottom trend of cross sections s1, s2, ..., si is a downward trend: ;

[0050] in, 、 is a correction factor related to the walking trend of the bottom corresponding to the cross-section, 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 temperature data, is the water flow velocity corresponding to the cross-section si;

[0051] (2) The part located downstream of the gate The calculation formula is:

[0052] When the bottom trends corresponding to the cross-sections s1, s2, …, si are upward trends: ;

[0053] When the bottom trends corresponding to the cross-sections s1, s2, …, si are downward trends ;

[0054] Among them, 、 are correction factors related to the walking trend of the bottom corresponding to the cross-section, 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 temperature data, is the water flow velocity corresponding to the cross-section si.

[0055] The alarm level determination module compares the S value with the 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;

[0056] The main control module selects a suitable dredging device and dredging method corresponding to the S value according to the alarm level;

[0057] The timing module sets a time interval and repeats the dredging determination.

[0058] The present invention has the following beneficial effects:

[0059] The present invention uses the volume of silt to characterize the amount of silt, which can more accurately reflect the underwater silt amount compared to 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 two aspects, horizontal and vertical, 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 upstream and downstream of the gate in this section and the flow velocity, water level, bottom trend, temperature, and the weights of each cross-section in this section are considered, 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 manners

[0060] In order 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 solutions and beneficial effects of the present invention, the following will be 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.

[0061] A high-precision cleaning determination method specifically includes the following steps:

[0062] (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;

[0063] (2) Uniformly arrange multiple silt thickness measurement units in each cross-section si of the river channel, 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;

[0064] (3) Calculate the average thickness avg i in the cross-section si based on the thickness values a1, a2 …… ai: avg i = (a1 + a2 + …… + ai) / i;

[0065] (4) Repeat steps (2) and (3) to obtain the average thickness values avg1, avg2 …… avg i in multiple cross-sections of the entire river channel;

[0066] (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, use the Kriging interpolation algorithm for interpolation and correct the obtained interpolation;

[0067] (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;

[0068] (7) Based on the above silt average thickness fitting curve, 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 fitting 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;

[0069] (8) Determine the silt volume of this length segment x, V = S × d, where d is the average width of the river channel;

[0070] (9) Correct V according to the position of the area to be determined;

[0071] (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;

[0072] (11) Select a suitable silt removal device and silt removal method corresponding to the V value according to the alarm level;

[0073] (12)定时循环执行步骤(3)至(10)。 (12)定时循环执行步骤(3)至(10)。 (12) Timingly and repeatedly execute steps (3) to (10).

[0074] As a preferred solution, the silt thickness measuring unit in step (2) uses ultrasound 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 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.

[0075] As a preferred solution, in step (5), the Kriging interpolation algorithm is used for interpolation as follows:

[0076] 5.1 Set the inland river section to be tested as A, and the river layer measurement value to be constructed as Z(avg), with {Z(avg)∈A}, Z(avg) is a second-order stationary function, and the values ​​of Z(avg) in space are set as Z(avg1), Z(avg2), ..., Z(avg n ), where avg represents the spatial location of the measurement point; according to the principle of the ordinary kriging algorithm, the estimated value of Z(avg0) 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 silt in cross section si avg i Function value, avg0 is the preliminary interpolation value obtained, and 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;

[0077] The ordinary Kriging equations are:

[0078] According to the minimum variance estimated by ordinary Kriging, the transformation becomes: =C - +u;

[0079] Among them, C represents the covariance function, E is the mathematical expectation value, u is the Lagrange multiplier;

[0080] 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-mentioned historical average thickness value vectors is used to correct the above-mentioned preliminary interpolation data, specifically: , is the feature vector after processing by the convolution model, (i∈[1,t]), , ,Pick ,and The closest 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 row j and column p in The Euclidean distance between them is: , j and p are 0, 1, 2, ..., n, to obtain the minimum Euclidean distance , the minimum Euclidean distance Corresponding element value The final interpolation data is ;

[0081] As a preferred solution, the curve fitting method in step (6) adopts the spline method to perform curve fitting;

[0082] 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 regularly determines the S value in the key area.

[0083] 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 silt volume corresponding to the part upstream of the gate, and V2 is the silt volume corresponding to the part downstream of the gate.

[0084] As a more preferred solution, the correction parameters in step (9) and The calculation formula is as follows:

[0085] (1) The part located upstream of the gate The calculation formula is:

[0086] When the bottom trend corresponding to the cross sections s1, s2, ..., si is an upward trend: ;

[0087] When the bottom trend of cross sections s1, s2, ..., si is a downward trend: ;

[0088] in, 、 is a correction factor related to the trend of the bottom of the cross section, and < , h is the average water level height of 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;

[0089] (2) The part located downstream of the gate The calculation formula is:

[0090] When the bottom trend corresponding to the cross sections s1, s2, ..., si is an upward trend: ;

[0091] The bottom trend corresponding to cross sections s1, s2, ..., si is a downward trend ;

[0092] in, 、 is a correction factor related to the trend of the bottom of the cross section, and < , h is the average water level height of 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.

[0093] A high-precision pollution removal determination method determination system, the determination system comprising:

[0094] A parameter monitoring module is used to measure the silt thickness value, wherein a plurality of silt thickness measuring units are arranged at equal intervals on the cross section in the longitudinal direction of the river channel, and the silt thickness values ​​in the plurality of cross sections are measured;

[0095] The first calculation module is used to obtain an 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;

[0096] 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 the average of the average thickness values ​​before and after it, the value will be 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.

[0097] A curve fitting module is used to obtain a fitting curve of average silt thickness, and to perform curve fitting using the average thickness value after data processing as the ordinate and the interval distance of each cross section as the abscissa to obtain a fitting curve of average silt thickness;

[0098] 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 using 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 x is obtained: V=S×d, where d is the average width of the river channel;

[0099] Correction module, 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 silt volume corresponding to the part located upstream of the gate, and V2 is the silt volume corresponding to the part located downstream of the gate);

[0100] Correction parameters and The calculation formula is as follows:

[0101] (1) The part located upstream of the gate The calculation formula is:

[0102] When the bottom trend corresponding to the cross sections s1, s2, ..., si is an upward trend: ;

[0103] When the bottom trend of cross sections s1, s2, ..., si is a downward trend: ;

[0104] in, 、 is a correction factor related to the walking trend of the bottom corresponding to the cross-section, 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;

[0105] (2) The part located downstream of the gate The calculation formula is:

[0106] When the bottom trend corresponding to the cross-sections s1, s2, …, si is an upward trend: ;

[0107] When the bottom trend corresponding to the cross-sections s1, s2, …, si is a downward trend ;

[0108] Among them, , are correction factors related to the walking trend of the bottom corresponding to the cross-section, 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.

[0109] 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 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 cleaning is required immediately;

[0110] The main control module selects a suitable dredging device and dredging method corresponding to the S value according to the alarm level;

[0111] The timing module sets a time interval and repeats the dredging determination.

[0112] The specific implementation is as follows:

[0113] 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), the distance between adjacent cross-sections is 1 meter;

[0114]

[0115] From Table 1, we can see that the thickness value 3m corresponding to section s3 is significantly greater than 5 times the average of the thickness value 0.1 of the previous section s2 and the thickness value 0.7 of the subsequent section s4 ((0.1+0.7) / 2), so it is eliminated;

[0116] Assume Z(avg)= , then 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 value of avg3 is , further corrected by using the correlation with the historical average data of the horizontal interface s3, and finally obtained the interpolation value of avg3 = 1.036;

[0117] Measure the silt volume in the area from s1 to s6, where s1-s3 are upstream of the gate and s4-s6 are downstream of the gate. Table 2: Silt volume V1 (unit: cubic meter) corresponding to the upstream and V2 (unit: cubic meter) corresponding to the downstream

[0118]

[0119] The corrected sludge volume is 1.637×1.2+0.32×0.8=3.0389 ,

[0120] It can be seen that the volume of silt in the area where sections s1 to s6 are located is greater than the emergency level threshold (V high = 3 ), then sludge cleaning is required immediately.

[0121] The above embodiments are used to explain the technical solutions of the present invention. However, 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 those skilled in the art based on the present invention, or equivalent substitutions of materials used in the present invention, etc., fall within the scope of protection of the patent.

Claims

1. A high-precision pollution determination method, characterized by: 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 channel to be measured. 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, where n < D / 5 and D is the total length value of the river channel to be measured; (2) Uniformly set multiple silt thickness measurement units in each cross-section si of the river channel, and obtain multiple silt thickness values a1, a2... ai according to the multiple silt thickness measurement units. i is a natural number greater than 1; (3) Calculate the average thickness avg i = (a1 + a2 +... + ai) / i within the cross-section si based on the thickness values a1, a2... ai; (4) Repeat steps (2) and (3) to obtain the average thickness values avg1, avg2... avg i within multiple cross-sections of the entire river channel; (5) Process the above average thickness values obtained in step (4). The specific processing method: when 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) 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 the 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 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 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 a suitable dredging device and dredging method corresponding to the V value according to the alarm level; (12)定时循环执行步骤(3)至(10)。 2. The high-precision pollution determination method according to claim 1, characterized in that: (12) Timely and circularly execute steps (3) to (10).

3. The high-precision pollution determination method according to claim 2, characterized in that: In step (2), the silt thickness measurement unit uses ultrasonic waves to measure the silt thickness.

4. The high-precision pollution determination method according to claim 3, characterized in that: In step (2), the silt thickness measurement unit uses ultrasonic waves to measure the silt thickness. It 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. In step (5), the specific method of using the Kriging interpolation algorithm for interpolation is as follows: 5.1 Set the inland river section to be tested as A, and the river layer measurement value to be constructed as Z(avg), with {Z(avg)∈A}, Z(avg) is a second-order stationary function, and the values ​​of Z(avg) in space are set as Z(avg1), Z(avg2), ..., Z(avg n ), where avg is the average thickness avg of the river channel of length x; according to the principle of the ordinary kriging algorithm, the estimated value of Z(avg0) 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 silt in cross section si avg i Function value, avg0 is the preliminary interpolation value obtained, and 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, the transformation becomes: =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-mentioned historical average thickness value vectors is used to correct the above-mentioned preliminary interpolation data, specifically: , is the feature vector after processing by the convolution model, (i∈[1,t]), , ,Pick ,and The closest 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 row j and column p in The Euclidean distance between them is: , j and p are 0, 1, 2, ..., n, to obtain the minimum Euclidean distance , the minimum Euclidean distance Corresponding element value The final interpolation data is .

5. The high-precision pollution 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 pollution 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, select the length value xi corresponding to the length of this area, and regularly determine the S value in this key area.

7. The high-precision pollution 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 silt volume corresponding to the part upstream of the gate, and V2 is the silt volume corresponding to the part downstream of the gate.

8. The high-precision pollution determination method according to claim 7, characterized in that: Correct the parameters in step (9) and The calculation formula is as follows: (1) The part located upstream of the gate The calculation formula is: When the corresponding bottom trends of the cross-sections s1, s2,..., si are upward trends: ； When the corresponding bottom trends of the cross-sections s1, s2,..., si are downward trends: ； in, 、 is a correction factor related to the trend of the bottom of the cross section, and < , h is the average water level height of 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) The part located downstream of the gate The calculation formula is: When the corresponding bottom trends of the cross-sections s1, s2,..., si are upward trends: ； When the corresponding bottom trends of the cross-sections s1, s2,..., si are downward trends: ； in, 、 is a correction factor related to the trend of the bottom of the cross section, and < , h is the average water level height of 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. A high-precision pollution cleaning determination system, characterized in that: The determination system includes: A parameter monitoring module, which is used to measure the silt thickness value, set multiple silt thickness measurement units at the same interval on the cross-section in the length direction of the river channel, and measure the silt thickness values in multiple cross-sections; 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 multiple 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 delete this value; if the average thickness value avg0 of a certain cross-section is missing or for the above-deleted average thickness values, then use the Kriging interpolation algorithm for interpolation and correct the obtained interpolation; A curve fitting module, which is used to obtain the silt average thickness fitting curve, perform curve fitting with the processed average thickness value as the ordinate and the interval distance of each cross-section as the abscissa, and obtain the silt average thickness fitting curve; 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 using 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 x is obtained: V=S×d, where d is the average width of the river channel; Correction module, 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 silt volume corresponding to the part located upstream of the gate, and V2 is the silt 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 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 in a hurry; if S high < S, it is the emergency level and immediate cleaning is required; A main control module, which selects the appropriate dredging device and dredging method corresponding to the S value according to the alarm level; A timing module, which sets the time interval and repeats the dredging determination.

Images