A method for measuring river waves

Abstract

The present invention relates to the field of water level measurement technology, and specifically discloses a river wave measurement method, including 9 detection points preset according to an expert knowledge base on a water level monitoring rod and a pressure sensor corresponding to the 9 detection points, the pressure sensor is used to feed back the real-time monitored electrical measurement value Y; the water level monitoring rod is also provided with at least two electrode wires arranged along the axial direction, the two electrode wires respectively collect the first electrical measurement value G1 and the second electrical measurement value G2 corresponding to the 9 detection points; the measurement method includes the following steps: S01, real-time data of the 9 detection points are collected in real time, the real-time data includes the time series h1(t) of the water level values ​​measured by the first electrode wire corresponding to the 9 detection points, and the time series h2(t) of the water level values ​​measured by the second electrode wire. The invention realizes the measurement of waves with shallow water depth, small wave height and short period in the river.

Description

Technical Field

[0001] The invention relates to the technical field of water level measurement, and in particular to a river wave measurement method. Background Art

[0002] At present, the standard calculation method for river waves adopts the research results based on sea areas, reservoirs or lakes. Compared with sea areas, reservoirs or lakes, the wind zone length of rivers is limited, the coastline is tortuous and changeable, and the water flow and wave effects are complex. The calculation accuracy and reliability of the Putian formula suitable for open water wind and wave calculation in river environments have not been fully tested and need to be corrected and improved based on field observation data.

[0003] At present, wave buoys, bottom-mounted acoustic wave meters, and pressure wave meters (RBR) can be used for wave measurement. Due to the characteristics of river wave measurement such as shallow water depth, small wave height, and short period, wave buoys have poor wave-following properties for waves with short period and small wave height, resulting in inaccurate measurements; bottom-mounted acoustic wave meters use acoustic principles, with a measurement resolution of 5 cm and are unable to derive wave sequences, and are generally used for offshore wave measurements; pressure wave meters use the pressure conversion principle to measure waves, and can derive wave sequences with a sampling rate of up to 16 Hz, but for waves with small wave height and short period, the pressure conduction loss is serious, and there is a serious peak elimination phenomenon, resulting in distorted measurement results. Summary of the invention

[0004] The purpose of the present invention is to provide a river wave measurement method for solving the above-mentioned problem.

[0005] In order to achieve the above object, the present invention provides the following technical solutions: a river wave measurement method, comprising 9 detection points preset according to an expert knowledge base and pressure sensors arranged one-to-one corresponding to the 9 detection points on a water level monitoring rod, wherein the pressure sensor is used to feed back the electrical measurement value Y monitored in real time;

[0006] The water level monitoring rod is also provided with at least two electrode wires arranged along the axial direction, and the two electrode wires respectively collect the first electrical measurement value G1 and the second electrical measurement value G2 corresponding to the nine detection points;

[0007] The measuring method comprises the following steps:

[0008] S01, real-time data of 9 detection points are collected in real time, wherein the real-time data includes a time series h1(t) of water level values ​​measured by the first electrode wire corresponding to the 9 detection points, a time series h2(t) of water level values ​​measured by the second electrode wire, and a time series h3(t) of water level data measured by the pressure sensor;

[0009] S02, performing function calculation on the acquired real-time data to eliminate abnormal values ​​and obtain an optimal data set;

[0010] S03. Based on the obtained optimal data set, the normal river data of the measured values ​​are traversed according to time.

[0011] |h1(t)-h2(t)|≤5mm (1);

[0012] (h1(t)+h2(t)) / 2-h3(t)≤10mm (2);

[0013] Based on the above formula, we can make the following judgment:

[0014] If formula (1) cannot be satisfied, but formula (2) is satisfied, then there is only one abnormal electrode wire data, then h2(t) or h1(t) = h3(t), which is used as the water level output value;

[0015] If neither formula (1) nor formula (2) is satisfied, (h1(t)+h2(t)) / 2 is used as the water level output value;

[0016] If formula (1) is satisfied, (h1(t)+h2(t)) / 2 is used as the water level output value;

[0017] S04, arranging the continuous time to obtain the water level height output value to obtain a water level sequence;

[0018] S05. Calculate the wave height and wave period for the acquired water level sequence, wherein the wave height and wave period include the maximum wave height hmax, the average wave height ha, the effective wave height hs, the average period ta and the effective period ts.

[0019] Preferably, the expert knowledge base is that the monitoring length of the water level monitoring rod is the highest water level and the lowest water level in the historical water level data of the monitored water area;

[0020] The historical water level data include water level data of the currently monitored river channel for nearly 50 years.

[0021] Preferably, the deployment of the nine detection points includes:

[0022] S61, determining the required monitoring height of the water level monitoring rod for the monitored water area;

[0023] S62, taking the middle value based on the determined monitoring height as the starting monitoring origin f5;

[0024] S63, taking the monitoring origin f5 as the center, taking four detection points upwards along the axis of the water level monitoring rod, namely f1, f2, f3 and f4, and taking four detection points downwards, namely f6, f7, f8 and f9;

[0025] The real-time data collected at each of the above-mentioned detection points include G1, G2, Y and h, where h is the calibrated height of the detection point.

[0026] Preferably, performing function calculation on the acquired real-time data to remove abnormal values ​​in step S02 includes:

[0027] S21, when there are at least two detection points at a certain moment as a group to obtain measurement values, the remaining detection points are used as n groups to obtain measurement values ​​in the future time period as a data set of 9 detection points, and the measurement values ​​include G1, G2 and Y;

[0028] S22, performing calculations on the measurement values ​​of the two detection points of each group in the acquired data set based on the discrete function in Excel to obtain

[0029] G1=fn(h) (3);

[0030] G2 = fn(h) (4);

[0031] Y = fn(h) (5);

[0032] Where f is the detection point, n is the position of the detection point corresponding to 1, 2, ..., 9, and h is the calibration height of the corresponding detection point;

[0033] S23, performing least squares polynomial curve fitting method on the three functions of the 9 detection points obtained to obtain continuous functions of the 9 detection points of the data set;

[0034] S24. Then perform the inverse function operation in Excel on the continuous function to obtain the water level wave.

[0035] Preferably, the measurement value obtained at each detection point fn in the optimal data set obtained in step S02 is used to perform the judgment.

[0036] If G1=G2=Y, it is determined as the water level output value, and step S03 is skipped to execute step S04;

[0037] If G1=G2=Y, then execute step S03.

[0038] Preferably, the step S05 of calculating the wave height and wave period of the acquired water level sequence comprises:

[0039] S51. Take h1(t) as the default water level value, use frequency H, and the continuously collected water level sequence value is [h1(0), h1(1 / H), h1(2 / H), h1(3 / H), ... h1(n / H)]. Traverse the entire array and compare the value with the previous and next values. The condition for determining the peak value is:

[0040] h1(m / H)>h1(m-1 / H), and h1(m / H)>h1(m+1 / H) (6);

[0041] h1(n / H) is less than h1(n-1 / H), and h1(n / H)>h1(n+1 / H) (7);

[0042] Among them, h(m / H) is a peak value, and h(n / H) traverses the array in the pre-entered database to obtain the peak value array [h1(m1 / H), h1(m2 / H), h1(m3 / H)...] and the valley value array [h1(n1 / H), h1(n2 / H), h1(n3 / H)...] arranged by time;

[0043] S52. Calculate the wave height values ​​and cycle values ​​between the peak values ​​through the peak value array. The calculation method is as follows:

[0044] Then the three sets of wave heights when all data are obtained at the 9 detection points are calculated as follows:

[0045] b1=h1(m1 / f)-h1(n1 / f) (8);

[0046] b2=h1(m2 / f)-h1(n2 / f) (9);

[0047] b3=h1(m3 / f)-h1(n3 / f) (10);

[0048] Among them, b is the wave height, 1, 2, and 3 are the combinations of any two detection points among the 9 detection points;

[0049] The three sets of cycle calculations for all acquired data at 9 detection points are as follows:

[0050] t1=(n1-m1) / H (11);

[0051] t2=(n2-m2) / H (12);

[0052] t3=(n3-m3) / H (13);

[0053] Among them, t is the period, 1, 2, and 3 are the combinations of any two detection points among the 9 detection points;

[0054] Get the wave height array [b1, b2, b3, b4.....] and the period array [t1, t2, t3, t4,...];

[0055] S53, the wave height array [b1, b2, b3, b4.....] is arranged from large to small, the first maximum wave height is hmax, the average value is the average wave height ha, and the 13% wave height value in the wave height array is the effective wave height hH;

[0056] S54. Arrange the periodic array [t1, t2, t3, t4, ...] from large to small, the average value is the average wave height ta, and the 13% wave height is the effective wave height ts.

[0057] In the above technical scheme, a river wave measurement method provided by the present invention has the following beneficial effects: wave height data about river waves are obtained through 9 preset detection points, and the data obtained from every two groups of detection points constitute a complete wave height, and the wave height is obtained by performing calculations on any three groups of data obtained, thereby obtaining the wave height and wave period obtained at the 9 detection points. BRIEF DESCRIPTION OF THE DRAWINGS

[0058] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in the present invention. For ordinary technicians in this field, other drawings can also be obtained based on these drawings.

[0059] Figure 1 This is a logic flow chart provided for the first embodiment of the present invention. DETAILED DESCRIPTION

[0060] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0061] Embodiment 1:

[0062] like Figure 1 As shown, a river wave measurement method includes: 9 detection points preset according to an expert knowledge base and pressure sensors arranged one-to-one corresponding to the 9 detection points on a water level monitoring rod, the pressure sensor is used to feed back the electrical measurement value Y monitored in real time;

[0063] At least two electrode wires arranged along the axial direction are also provided on the water level monitoring rod, and the two electrode wires respectively collect the first electrical measurement value G1 and the second electrical measurement value G2 corresponding to the nine detection points.

[0064] In the implementation process of this embodiment, G1, G2 and Y of each detection point are actually collected, and the height values ​​h of the 9 detection points in the embodiment are fixed. The height values ​​h of the 9 detection points are determined by an expert knowledge base, and the expert knowledge base in the above embodiment can be data obtained by water level monitoring before the water level monitoring pole is installed, or can be divided equally according to the annual average value of the river water level height, etc., and any method known to those skilled in the art can be used, and no detailed description will be given here.

[0065] Therefore, the above embodiment is completed by obtaining the changes in the electrical measurement values ​​of the two electrode wires and the pressure sensor under a complete wave collected at 9 detection points.

[0066] The specific measurement method provided in the above embodiment includes the following steps:

[0067] S01, real-time data of 9 detection points are collected in real time, the real-time data includes the time series h1(t) of the water level values ​​measured by the first electrode wire corresponding to the 9 detection points, the time series h2(t) of the water level values ​​measured by the second electrode wire, and the time series h3(t) of the water level data measured by the pressure sensor;

[0068] S02, performing function calculation on the acquired real-time data to eliminate abnormal values ​​and obtain an optimal data set;

[0069] S03. Based on the obtained optimal data set, the normal river data of the measured values ​​are traversed according to time.

[0070] |h1(t)-h2(t)|≤5mm (1);

[0071] (h1(t)+h2(t)) / 2-h3(t)≤10mm (2);

[0072] Based on the above formula, we can make the following judgment:

[0073] If formula (1) cannot be satisfied, but formula (2) is satisfied, then there is only one abnormal electrode wire data, then h2(t) or h1(t) = h3(t), which is used as the water level output value;

[0074] If neither formula (1) nor formula (2) is satisfied, (h1(t)+h2(t)) / 2 is used as the water level output value;

[0075] If formula (1) is satisfied, (h1(t)+h2(t)) / 2 is used as the water level output value;

[0076] S04, sorting the continuous time to obtain the water level height output value to obtain the water level sequence;

[0077] S05. Calculate the wave height and wave period for the acquired water level sequence, where the wave height and wave period include the maximum wave height hmax, the average wave height ha, the effective wave height hs, the average period ta and the effective period ts.

[0078] The detailed steps for calculating the wave height and wave period in this step are as follows:

[0079] S51. Take h1(t) as the default water level value and the frequency as H. If the continuously collected water level sequence values ​​are [h1(0), h1(1 / H), h1(2 / H), h1(3 / H), ..... h1(n / H)], traverse the entire array and compare the value with the previous and next values. The condition for determining the peak value is:

[0080] h1(m / H)>h1(m-1 / H), and h1(m / H)>h1(m+1 / H) (6);

[0081] h1(n / H) is less than h1(n-1 / H), and h1(n / H)>h1(n+1 / H) (7);

[0082] Among them, h(m / H) is a peak value, and h(n / H) traverses the array in the pre-entered database to obtain the peak value array [h1(m1 / H), h1(m2 / H), h1(m3 / H)...] and the valley value array [h1(n1 / H), h1(n2 / H), h1(n3 / H)...] arranged by time;

[0083] S52. Calculate the wave height values ​​and cycle values ​​between the peak values ​​through the peak value array. The calculation method is as follows:

[0084] Then the three sets of wave heights when all data are obtained at the 9 detection points are calculated as follows:

[0085] b1=h1(m1 / f)-h1(n1 / f) (8);

[0086] b2=h1(m2 / f)-h1(n2 / f) (9);

[0087] b3=h1(m3 / f)-h1(n3 / f) (10);

[0088] Among them, b is the wave height, 1, 2, and 3 are the combinations of any two detection points among the 9 detection points;

[0089] The three sets of cycle calculations for all acquired data at 9 detection points are as follows:

[0090] t1=(n1-m1) / H (11);

[0091] t2=(n2-m2) / H (12);

[0092] t3=(n3-m3) / H (13);

[0093] Among them, t is the period, 1, 2, and 3 are the combinations of any two detection points among the 9 detection points;

[0094] Get the wave height array [b1, b2, b3, b4.....] and the period array [t1, t2, t3, t4,...];

[0095] S53, the wave height array [b1, b2, b3, b4.....] is arranged from large to small, the first maximum wave height is hmax, the average value is the average wave height ha, and the 13% wave height value in the wave height array is the effective wave height hH;

[0096] S54. Arrange the periodic array [t1, t2, t3, t4, ...] from large to small, the average value is the average wave height ta, and the 13% wave height is the effective wave height ts.

[0097] Specifically, the monitoring length of the water level monitoring rod under the expert knowledge base in the above embodiment is the highest water level and the lowest water level in the historical water level data of the monitored water area; and the historical water level data includes the water level data of the current monitored river in the past 50 years. In this embodiment, the highest water level and the lowest water level can be the highest water level and the lowest water level recorded in history, or the highest water level and the lowest water level of each year, and the average value of the highest water level and the average value of the lowest water level in the past 50 years are calculated.

[0098] Furthermore, in the actual operation process of the above embodiment, when the wave triggers any of the 9 detection points, the first electrode wire, the second electrode wire and the pressure sensor respectively acquire electrical signals, and then based on each detection point for acquiring data, the electrical signals are used as a set of data. If there are multiple detection points, multiple sets of data are acquired. It should be noted that in the above embodiment, a complete wave at least triggers two detection points, that is, a set of data. After the data is acquired, data denoising is performed in step S02, that is, outliers are removed. The implementation method includes but is not limited to the K nearest neighbor (KNN) algorithm, density-based clustering method and isolation forest and other signal denoising processing methods known in the art.

[0099] Normal river data refers to the river wave height data predicted by the water level monitoring poles. This data collection requires a collection of multiple sample data collected over at least one month or more.

[0100] In the above technology, wave height data about river waves are obtained through 9 preset detection points, and the data obtained from every two groups of detection points constitute a complete wave height. The wave height is obtained by performing calculations on any three groups of data obtained, thereby obtaining the wave height and wave period obtained at the 9 detection points.

[0101] Embodiment 2:

[0102] Based on the first embodiment, this embodiment aims to provide a method for controlling and processing 9 detection points, including:

[0103] S61, determining the required monitoring height of the water level monitoring pole for the monitored water area;

[0104] S62, taking the middle value based on the determined monitoring height as the starting monitoring origin f5;

[0105] S63, taking the monitoring origin f5 as the center, taking four detection points upwards along the axis of the water level monitoring rod, namely f1, f2, f3 and f4, and taking four detection points downwards, namely f6, f7, f8 and f9;

[0106] The real-time data collected at the above detection points include G1, G2, Y and h, where h is the calibrated height of the detection point.

[0107] Specifically, the monitoring height is determined based on the average water level of the river in each season and month of the year according to historical data, so as to achieve dynamic monitoring. In short, the 9 detection points will adapt to the changes in the monitoring height according to the season or month of the monitored river, so as to cooperate in monitoring to ensure the accuracy of the data collected throughout the year.

[0108] Furthermore, in the process of acquiring detection point data based on the above implementation, the lowest point of each wave corresponds to an electrical measurement value of any detection point f6, f7, f8 and f9, and the highest point corresponds to an electrical measurement value of any detection point f1, f2, f3 and f4, forming a data group, and there are at most four groups of data (excluding f5). In the detection method provided in the first embodiment of this implementation, it is necessary to obtain data at all 9 detection points before a complete data collection is completed. After a complete data collection is completed, it is necessary to determine three data groups according to the time sequence mechanism to perform the processing of steps S02 to S05.

[0109] Embodiment three:

[0110] Based on the first embodiment, this embodiment aims to perform function calculation on the acquired real-time data in step S02 to remove abnormal values, including:

[0111] S21, when there are at least two detection points at a certain moment as a group to obtain measurement values, the remaining detection points are used as n groups to obtain measurement values ​​in the future time period, as a data set of 9 detection points, the measurement values ​​include G1, G2 and Y;

[0112] S22, the measured values ​​of the two detection points of each group in the acquired data set are calculated based on the discrete function in Excel, and the obtained values ​​are obtained.

[0113] G1=fn(h) (3);

[0114] G2 = fn(h) (4);

[0115] Y = fn(h) (5);

[0116] Where f is the detection point, n is 1, 2, ..., 9 as the position of the corresponding detection point, and h is the calibration height of the corresponding detection point;

[0117] S23, performing least squares polynomial curve fitting method on the three functions of the 9 detection points obtained to obtain continuous functions of the 9 detection points of the data set;

[0118] S24. Then perform the inverse function operation in Excel on the continuous function to obtain the water level wave.

[0119] Preferably, the measurement value obtained at each detection point fn in the optimal data set obtained in step S02 is used to perform the judgment.

[0120] If G1=G2=Y, it is determined as the water level output value, and step S03 is skipped to execute step S04;

[0121] If G1=G2=Y, then execute step S03.

[0122] Specifically, when there are at least two detection points at a certain moment as a group to obtain measurement values, that is, there are two corresponding detection points f1-f4 and f6-f9 to obtain data, and the data is recorded, then a first data group is formed, and the measurement values ​​obtained from subsequent detection points are recorded. After the measurement values ​​are obtained at all 9 detection points, the multiple data obtained are recorded in chronological order, such as the difference between (f1, f9), the difference between (f2, f8), etc. After all 9 detection points have obtained data, the data is entered into Excel, and then each obtained data is processed through the built-in discrete function in Excel, so as to generate multiple scattered points under the preset coordinates.

[0123] Furthermore, based on the obtained multiple scattered points, least squares polynomial curve fitting is performed (executed using Excel preset algorithm) to obtain a continuous function of the 9 detection points in the data set, that is, multiple scattered points are connected to obtain a curve.

[0124] Furthermore, in the above embodiment, during the process of collecting the electrical values ​​of the 9 detection points, i.e., during the process of executing steps S21 to S24, a factor is given priority, i.e., whether the data obtained from each detection point is G1=G2=Y. If not, step S03 needs to be executed. This avoids the situation in which the pressure sensor and the two electrode wires are damaged during the data collection process, resulting in inconsistent collected data.

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

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

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

[0128] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0129] The present invention uses specific embodiments to illustrate the principles and implementation methods of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation on the present invention.

[0130] The embodiments of the present application also provide a specific implementation of an electronic device capable of implementing all steps of the method in the above embodiments, and the electronic device specifically includes the following contents:

[0131] Processor, memory, communications interface and bus;

[0132] Wherein, the processor, memory, and communication interface communicate with each other via the bus;

[0133] The processor is used to call the computer program in the memory, and when the processor executes the computer program, all the steps in the method in the above embodiment are implemented.

[0134] The embodiments of the present application also provide a computer-readable storage medium capable of implementing all the steps of the method in the above embodiments, wherein a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, all the steps of the method in the above embodiments are implemented.

[0135] Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the hardware + program embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment. Although the embodiment of this specification provides the method operation steps described in the embodiment or flow chart, more or less operation steps can be included based on conventional or non-creative means. The order of steps listed in the embodiment is only one way of executing the order of many steps, and does not represent the only execution order. When the device or terminal product in practice is executed, it can be executed in sequence or in parallel according to the method shown in the embodiment or the accompanying drawings (for example, a parallel processor or a multi-threaded processing environment, or even a distributed data processing environment). The term "include", "include" or any other variant thereof is intended to cover non-exclusive inclusion, so that the process, method, product or device including a series of elements includes not only those elements, but also includes other elements that are not explicitly listed, or also includes elements inherent to such a process, method, product or device. In the absence of more restrictions, it is not excluded that there are other identical or equivalent elements in the process, method, product or device including the elements. For the convenience of description, the above device is described by dividing it into various modules according to its functions. Of course, when implementing the embodiments of this specification, the functions of each module can be implemented in the same one or more software and / or hardware, or the module implementing the same function can be implemented by a combination of multiple sub-modules or sub-units. The device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the coupling or direct coupling or communication connection between each other shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms. The present invention is described with reference to the flowchart and / or block diagram of the method, device (system) and computer program product according to the embodiments of the present invention. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the process and / or box in the flowchart and / or block diagram can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the process Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0136] Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment. In the description of this specification, the description of the reference terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the embodiments of this specification.

[0137] In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradicting each other. The above is only an embodiment of the embodiment of this specification and is not intended to limit the embodiment of this specification. For those skilled in the art, the embodiment of this specification may have various changes and variations. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the embodiment of this specification shall be included in the scope of the claims of the embodiment of this specification.

Claims

1. A method for measuring river waves, characterized in that: It includes 9 detection points preset according to the expert knowledge base on the water level monitoring rod and pressure sensors corresponding to the 9 detection points, and the pressure sensors are used to feed back the electrical measurement value Y monitored in real time; The water level monitoring rod is also provided with at least two electrode wires arranged along the axial direction, and the two electrode wires respectively collect the first electrical measurement value G1 and the second electrical measurement value G2 corresponding to the nine detection points; The measuring method comprises the following steps: S01, real-time data of 9 detection points are collected in real time, wherein the real-time data includes a time series h1(t) of water level values ​​measured by the first electrode wire corresponding to the 9 detection points, a time series h2(t) of water level values ​​measured by the second electrode wire, and a time series h3(t) of water level data measured by the pressure sensor; S02, performing function calculation on the acquired real-time data to eliminate abnormal values ​​and obtain an optimal data set; S03. Based on the obtained optimal data set, the normal river data of the measured values ​​are traversed according to time. |h1(t)-h2(t)|≤5mm (1); (h1(t)+h2(t)) / 2-h3(t)≤10mm (2); Based on the above formula, we can make the following judgment: If formula (1) cannot be satisfied, but formula (2) is satisfied, then there is only one abnormal electrode wire data, then h2(t) or h1(t) = h3(t), which is used as the water level output value; If neither formula (1) nor formula (2) is satisfied, (h1(t)+h2(t)) / 2 is used as the water level output value; If formula (1) is satisfied, (h1(t)+h2(t)) / 2 is used as the water level output value; S04, arranging the continuous time to obtain the water level height output value to obtain a water level sequence; S05, calculating the wave height and wave period of the acquired water level sequence, wherein the wave height and wave period include the maximum wave height hmax, the average wave height ha, the effective wave height hs, the average period ta and the effective period ts; The step S02 of performing function calculation on the acquired real-time data to remove abnormal values ​​includes: S21, when there are at least two detection points at a certain moment as a group to obtain measurement values, the remaining detection points are used as n groups to obtain measurement values ​​in the future time period as a data set of 9 detection points, and the measurement values ​​include G1, G2 and Y; S22, performing calculations on the measurement values ​​of the two detection points of each group in the acquired data set based on the discrete function in Excel to obtain G1=fn(h) (3); G2 = fn(h) (4); Y = fn(h) (5); Where f is the detection point, n is the position of the detection point corresponding to 1, 2, ..., 9, and h is the calibration height of the corresponding detection point; S23, performing least squares polynomial curve fitting method on the three functions of the 9 detection points obtained to obtain continuous functions of the 9 detection points of the data set; S24, then perform an inverse function operation in Excel on the continuous function to obtain the water level wave; The step S05 of calculating the wave height and wave period of the acquired water level sequence comprises: S51. Take h1(t) as the default water level value, use frequency H, and the continuously collected water level sequence value is [h1(0), h1(1 / H), h1(2 / H), h1(3 / H), ... h1(n / H)]. Traverse the entire array and compare the value with the previous and next values. The condition for determining the peak value is: h1(m / H)>h1(m-1 / H), and h1(m / H)>h1(m+1 / H) (6); h1(n / H) is less than h1(n-1 / H), and h1(n / H)>h1(n+1 / H) (7); Among them, h(m / H) is a peak value, and h(n / H) traverses the array in the pre-entered database to obtain the peak value array [h1(m1 / H), h1(m2 / H), h1(m3 / H)...] and the valley value array [h1(n1 / H), h1(n2 / H), h1(n3 / H)...] arranged by time; S52. Calculate the wave height values ​​and cycle values ​​between the peak values ​​through the peak value array. The calculation method is as follows: The three sets of wave heights when all data are obtained at the 9 detection points are calculated as follows: b1=h1(m1 / f)-h1(n1 / f) (8); b2=h1(m2 / f)-h1(n2 / f) (9); b3=h1(m3 / f)-h1(n3 / f) (10); Among them, b is the wave height, 1, 2, and 3 are the combinations of any two detection points among the 9 detection points; The three sets of cycle calculations for all acquired data at 9 detection points are as follows: t1=(n1-m1) / H (11); t2=(n2-m2) / H (12); t3=(n3-m3) / H (13); Among them, t is the period, 1, 2, and 3 are the combinations of any two detection points among the 9 detection points; Get the wave height array [b1, b2, b3, b4.....] and the period array [t1, t2, t3, t4,...]; S53, the wave height array [b1, b2, b3, b4.....] is arranged from large to small, the first maximum wave height is hmax, the average value is the average wave height ha, and the 13% wave height value in the wave height array is the effective wave height hH; S54. Arrange the periodic array [t1, t2, t3, t4, ...] from large to small, the average value is the average wave height ta, and the 13% wave height is the effective wave height ts.

2. A river wave measurement method according to claim 1, characterized in that: The expert knowledge base is that the monitoring length of the water level monitoring rod is the highest water level and the lowest water level in the historical water level data of the monitored water area; The historical water level data include water level data of the currently monitored river channel for nearly 50 years.

3. A river wave measurement method according to claim 1, characterized in that: The deployment of the 9 detection points includes: S61, determining the required monitoring height of the water level monitoring rod for the monitored water area; S62, taking the middle value based on the determined monitoring height as the starting monitoring origin f5; S63, taking the monitoring origin f5 as the center, taking four detection points upwards along the axis of the water level monitoring rod, namely f1, f2, f3 and f4, and taking four detection points downwards, namely f6, f7, f8 and f9; The real-time data collected at each of the above-mentioned detection points include G1, G2, Y and h, where h is the calibrated height of the detection point.

4. A river wave measurement method according to claim 1, characterized in that: Performing judgment on the measurement value obtained at each detection point fn in the optimal data set obtained in step S02, If G1=G2=Y, it is determined as the water level output value, and step S03 is skipped to execute step S04; If G1=G2=Y, then execute step S03.

5. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the steps of the river wave measurement method according to any one of claims 1 to 4 are implemented.

6. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the river wave measurement method according to any one of claims 1 to 4 are implemented.

Images