CN103090855A - Method for determining arrival of tidal bore based on water velocity - Google Patents
Method for determining arrival of tidal bore based on water velocity Download PDFInfo
- Publication number
- CN103090855A CN103090855A CN201310019205XA CN201310019205A CN103090855A CN 103090855 A CN103090855 A CN 103090855A CN 201310019205X A CN201310019205X A CN 201310019205XA CN 201310019205 A CN201310019205 A CN 201310019205A CN 103090855 A CN103090855 A CN 103090855A
- Authority
- CN
- China
- Prior art keywords
- speed data
- flow velocity
- flow speed
- flow
- state
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method for determining the arrival of tidal bore based on water velocity. The method specifically comprises the following steps of: correctly mounting a radar wave flow velocity meter on the surface of a river; collecting the flow velocity data; determining whether or not to switch the flow velocity data from state 1 to state 2 according to whether the flow velocity increases; determining whether or not to switch the flow velocity data from state 2 to state 3 according to whether the flow velocity decreases; determining whether or not to switch the flow velocity data from state 3 to state 4 according to whether the flow velocity is stable; determining that whether the reason causing the flow velocity change is interference or arrival of the tidal bore according to the difference between the stabilized flow velocity after changing and the flow velocity before changing; and determining the level of the tidal bore according to the magnitude of the difference between the maximal flow velocity in the change process and the flow velocity before changing. The method provided by the invention not only can determine the arrival of a tidal bore, but also can determine the flow velocity of the tidal bore, thereby providing necessary data for forecasting the arrival time of the tidal bore to the next station, and comprehensively realizing real-time automatic tidal bore detection and warning of rivers.
Description
Technical field
The invention belongs to technical field of automation, specifically a kind ofly judge based on flow rate of water flow the method whether tidal bore arrives, can be used for river tidal bore monitoring and real-time prediction early warning.
Background technology
The Qiantang River flows into the East Sea at Hangzhou Wan, and wide interior narrow outside the river mouth, wide place reaches 100 kms, and narrow place only has several kms.Seawater has just entered the Haikou, and the water surface is wide, and more more have to gush by the constraint of two sides, river landform, seawater and amass, and more long-pending higher, the river of Qiantang River outflow is stopped in addition, promotes again water level rise, as one upright waterwall.When the seawater of back constantly pours in, just form back rise and caught up with front wave, the grand scene of a folded wave of wave, surging, spectacular, become the whole world the most grand maximum morning and evening tides.Meanwhile, Qianjiang tidewater volume people event was almost never stopped over the past thousands of years, and just have more than ten people because of the wrecked number of tidewater every year.For preventing Qiantang River tidewater volume people, local government and personage among the people have organized " crying out damp team " specially, and establish the warning poster at dangerous area, warn the visitor to take care.Although certain effect is arranged, Qiantang River tidewater " is eaten people's event " and is never stopped.Therefore, the early warning of river tidal bore monitoring and prediction is for just become very meaningful of the safety of protection visitor and people across the Straits' life and property.
General, in river tidal bore monitoring and prediction early warning technology, whether needs contact tidewater according to checkout equipment, are divided into contact measurement and non-contact detection.For contact measurement, because marine corrosion can cause the equipment loss, on the other hand, for the large tidal bore rivers as the Qiantang River, the impulsive force of its morning and evening tides makes for the installation requirement of equipment very high on the one hand.Therefore present general employing non-contact detection, the existence limitation separately but existing technology all remains unchanged on using.For example, for the river tidal surge detection method based on audio frequency characteristics identification, in the situation that tidewater is less, audio frequency characteristics is obvious not, can't reach the purpose of identification tidal bore.Be unwell to night or river surface has the situations such as mist based on the river tidal surge detection method of video identification.In addition, existing technology only is confined to judge whether tidal bore arrives, and can't judge that the rank of tidal bore is big or small, is unfavorable for forecasting that tidal bore arrives the time of next website, makes real time automatic detection and the early warning of river tidal bore comprehensive not.
Summary of the invention
For the limitation of prior art, the invention provides and a kind ofly judge based on flow rate of water flow the method whether tidal bore arrives.The flow rate of water flow data that the present invention utilizes the radar wave current meter to record judge whether tidal bore arrives, and judge simultaneously the rank size of tidal bore, realize real time automatic detection and the alarm of river tidal bore.
The present invention includes following steps:
Step 1: select suitable place that the radar wave current meter is installed and measure flow rate of water flow, install on the bank iron shelf with the radar frame on river surface, make radar over against water (flow) direction.
Step 2: gather the instantaneous velocity data, two variablees are set are used for the flow speed data that record collects, be respectively flow speed data old and flow speed data new, to flow speed data old, the 2nd data assignment is to flow speed data new with the 1st data assignment collecting.
Step 3: calculate the difference of flow velocity data n ew and flow speed data old, relatively the size of the threshold value threshold1 of this difference and setting is to judge whether current state changes state two over to by state one; If this difference judges that less than or equal to the threshold value threshold1 that sets current state as state one, enters step 4; If this difference begins to judge whether to get the hang of two greater than the threshold value threshold1 that sets, give original flow velocity original with flow speed data old assignment, wherein the original flow velocity original of variable is used for the flow velocity of recording status one, enters step 5.
Step 4: flow speed data new assignment to flow speed data old, is then read the new data of radar wave current meter and assignment to flow speed data new, go back to step 3.
Step 5: for Peak Flow Rate max initialize is flow speed data new, wherein Peak Flow Rate max will be used for recording the maximum data that occurs in the tidal bore arrival process; Be 1 for being used for recording the integer variable increase assignment that has had several data to increase, expression has had 1 data to increase, wherein the size of integer variable increase can be used for characterizing two the confidence level of getting the hang of, the larger confidence level of numerical value is higher, when increasing to designated value n1, be judged to be and got the hang of two.
Step 6: read the new data of radar wave current meter and assignment to flow speed data new, calculate the difference of flow velocity data n ew and original flow velocity original, the size of the threshold value threshold1 of this difference and setting relatively, be divided into situation 1. with situation 2.:
Situation is 1.: this difference is greater than the threshold value threshold1 that sets, and judges that flow speed data new is whether greater than Peak Flow Rate max: if greater than Peak Flow Rate max, with its assignment to Peak Flow Rate max, otherwise Peak Flow Rate max remains unchanged; Integer variable increase is added 1, then judge whether integer variable increase has arrived designated value n1: if not yet arrive designated value n1, again change step 6 over to, if arrived designated value n1, judge and got the hang of two, therefore next begin to judge whether to get the hang of three, the integer variable decrease assignment that several data reduce has been arranged is 0 in order to be used for recording, and enters step 7; Wherein the size of integer variable decrease is used for characterizing three the confidence level of having got the hang of, and the larger confidence level of numerical value is higher, is judged to be when increasing to designated value n2 and has got the hang of three.
Situation is 2.: this difference judges that whether integer variable increase is greater than 1: if greater than 1, integer variable increase is subtracted 1, reenter step 6 less than or equal to the threshold value threshold1 that sets; If be not more than 1, judge not get the hang of two, namely be still state one, so integer variable increase assignment is 0, return to step 4.
Step 7: read the new data of radar wave current meter and assignment to flow speed data new, judge that flow speed data new is whether greater than Peak Flow Rate max: if greater than Peak Flow Rate max, with its assignment to Peak Flow Rate max, otherwise Peak Flow Rate max remains unchanged; Calculate the difference of Peak Flow Rate max and flow speed data new, relatively whether the size of the threshold value threshold2 of this difference and setting changes state three over to by state two with the residing state of decision data, is divided into situation I and situation II; Wherein threshold value threshold2 can select according to Peak Flow Rate max.
The situation I: this difference adds 1 greater than threshold value threshold2 with integer variable decrease, judges then whether integer variable decrease has arrived designated value n2: if not yet arrive designated value n2, again change step 7 over to; If arrived designated value n2, judge to have got the hang of three, therefore next begin to judge whether to get the hang of four, the integer variable steady of several data stabilizations assignment has been arranged is 0 in order to be used for recording, flow speed data new assignment to flow speed data old, is entered step 8; Wherein the size of integer variable steady is used for characterizing four the confidence level of having got the hang of, and the larger confidence level of numerical value is higher, is judged to be when increasing to designated value n3 and has got the hang of four.
The situation II: this difference judges that whether integer variable decrease is greater than 1: if greater than 1, integer variable decrease is subtracted 1, reenter step 7 less than or equal to threshold value threshold2; If be not more than 1, judge not get the hang of three, be still namely state two, be therefore 0 with integer variable increase assignment, then return to step 6.
Step 8: read the new data of radar wave current meter and assignment to flow speed data new, and calculate the absolute value of the difference of flow speed data old and flow speed data new; Relatively whether the size of the threshold value threshold3 of this absolute value and setting changes state four over to by state three with the residing state of decision data, is divided into situation a and situation b.
Situation a: this absolute value judges that whether integer variable steady is greater than 1: if greater than 1, integer variable steady is subtracted 1, reenter step 8 more than or equal to threshold value threshold3; If be not more than 1, judge not get the hang of four, be still namely state three, be therefore 0 with integer variable decrease assignment, then return to step 7.
Situation b: this absolute value is less than threshold value threshold3, integer variable steady is added 1, then judge whether integer variable steady has arrived designated value n3: if not yet arrive designated value n3, with flow speed data new assignment to flow speed data old, again change step 8 over to; If arrived designated value n3, judge to have got the hang of four, calculate the difference of flow velocity data n ew and original flow velocity original, relatively the size of this difference and threshold value threshold4, be divided into situation b-1 and situation b-2; The flow velocity of flow speed data new recording status four wherein, the flow velocity of original flow velocity original recording status one.
Situation b-1: if this difference is greater than threshold value threshold4, show that causing the whole variation of this time flow speed data is not to be caused by interference, therefore judge that tidal bore arrives, and according to the rank size of the size judgement tidal bore of the difference of Peak Flow Rate max and original flow velocity original, wherein max has recorded the Peak Flow Rate that reaches in whole change in flow process, and original flow velocity original has recorded the flow velocity of state one.
Situation b-2: if this difference is less than or equal to threshold value threshold4, show that causing the whole variation of this time flow speed data is to be caused by interference, therefore judge the state one that is in, therefore with flow speed data new assignment to flow speed data old, read the new data of radar wave current meter and assignment to flow speed data new, return to step 3.
Described state one not yet arrives for tidal bore, and state two increases for flow velocity, and state three is for flow velocity reduces, and state four is flow speed stability.
The invention has the beneficial effects as follows:
1. not only can judge whether tidal bore arrives, and obtain simultaneously the flow velocity of tidal bore, the time that arrives next website for the forecast tidal bore provides necessary data, has realized more all sidedly real time automatic detection and the early warning of river tidal bore.
2. due to the existence of disturbing, may cause erroneous judgement if directly continue to judge greater than certain threshold value whether tidal bore arrives according to flow velocity, so whether change to eliminate erroneous judgement according to the front and back data at last.
3. introduce and be used for recording the integer variable increase that had several data to increase, be used for recording and the integer variable decrease that several data reduce arranged and be used for recording the integer variable steady that several data stabilizations have been arranged, characterize respectively with its size and got the hang of two, the confidence level of state three and state four.So can effectively reject the interference of accidentalia, make for whether got the hang of two, the judgement of state three and state four is more rigorous, increased the accuracy of judgement.
Description of drawings
The change procedure of flow velocity when Fig. 1 is tidal bore arrival.
Fig. 2 is the inventive method process flow diagram.
Embodiment
Below in conjunction with accompanying drawing, method of the present invention is further described.
As shown in Figure 1, analyze the situation of flow velocity when in many experiments, tidal bore arrives, when the most obvious feature is tidal bore arrival as can be known, at first the unexpected moment change of flow velocity is large, begin after a period of time to diminish, keep for a long time a less stable flow velocity (flow speed value of this value when not yet arriving after being varied down to a certain degree, can eliminate with this and disturb the erroneous judgement that causes, if because interference causes change in flow, disturbing the flow restoration after disappearing so is original flow velocity).The present invention judges by this feature whether tidal bore arrives.Be that in the tidal bore arrival process, flow speed data is divided into one of four states: 1. tidal bore not yet arrives, and 2. flow velocity increases, and 3. flow velocity reduces, 4. flow speed stability.The present invention is in when judgement, if find that flow speed data has intactly experienced this one of four states, and the flow velocity after vary stable is judged as tidal bore and arrives greater than the flow velocity before changing, and judges simultaneously the rank size of this time tidal bore.
As shown in Figure 2, a kind ofly judge based on flow rate of water flow the method whether tidal bore arrives, comprise the following steps:
Step 1: select suitable place that the radar wave current meter is installed and measure flow rate of water flow, install on the bank iron shelf with the radar frame on river surface (T-shaped dam and beach are avoided in attention), make radar over against water (flow) direction.
Step 2: gather instantaneous velocity data (these instantaneous velocity data are to issue microprocessor by serial ports by the RS232 interface by the radar wave current meter), two variablees are set are used for the flow speed data that record collects, be respectively flow speed data old and flow speed data new, to flow speed data old, the 2nd data assignment is to flow speed data new with the 1st data assignment collecting.
Step 3: the difference (new-old) of calculating flow velocity data n ew and flow speed data old.Relatively the size of the threshold value threshold1 of this difference and setting is to judge whether current state changes state two (flow velocity increase) over to by state one (tidal bore not yet arrives).If this difference judges that less than or equal to the threshold value threshold1 that sets current state as state one (tidal bore not yet arrives), enters step 4; If this difference is greater than the threshold value threshold1 that sets, begin to judge whether to get the hang of two (flow velocity increases), give original flow velocity original with flow speed data old assignment, flow velocity when wherein the original flow velocity original of variable is used for recording status one (tidal bore not yet arrives) enters step 5.
Step 4: flow speed data new assignment to flow speed data old, is then read the new data of radar wave current meter and assignment to flow speed data new, go back to step 3.
Step 5: for Peak Flow Rate max initialize is flow speed data new, wherein Peak Flow Rate max will be used for recording the maximum data that occurs in the tidal bore arrival process.Be 1 for being used for recording the integer variable increase assignment that has had several data to increase, expression has had 1 data to increase, wherein the size of integer variable increase can be used for characterizing the confidence level of got the hang of two (flow velocity increases), the larger confidence level of numerical value is higher, when increasing to designated value n1, be judged to be two (the flow velocity increases) that got the hang of.
Step 6: read the new data of radar wave current meter and assignment to flow speed data new, calculate the difference (new-original) of flow velocity data n ew and original flow velocity original, the size of the threshold value threshold1 of this difference and setting relatively, be divided into situation 1. with situation 2..
Situation is 1.: this difference is greater than the threshold value threshold1 that sets, and judges that flow speed data new is whether greater than Peak Flow Rate max: if greater than Peak Flow Rate max, with its assignment to Peak Flow Rate max, otherwise Peak Flow Rate max remains unchanged; Integer variable increase is added 1, then judge whether integer variable increase has arrived designated value n1: if not yet arrive designated value n1, again change step 6 over to, if arrived designated value n1, judge two (the flow velocity increases) that got the hang of, therefore next begin to judge whether to get the hang of three (flow velocity reduces), the integer variable decrease assignment that several data reduce has been arranged is 0 in order to be used for recording, and enters step 7.Wherein the size of integer variable decrease is used for characterizing the confidence level of got the hang of three (flow velocity reduces), and the larger confidence level of numerical value is higher, is judged to be three (flow velocity reduces) that got the hang of when increasing to designated value n2.
Situation is 2.: this difference judges that whether integer variable increase is greater than 1: if greater than 1, integer variable increase is subtracted 1, reenter step 6 less than or equal to the threshold value threshold1 that sets; If be not more than 1, judge two (the flow velocity increases) that do not get the hang of, namely be still state one (tidal bore not yet arrives), so integer variable increase assignment is 0, return to step 4.
Step 7: read the new data of radar wave current meter and assignment to flow speed data new, judge that flow speed data new is whether greater than Peak Flow Rate max: if greater than Peak Flow Rate max, with its assignment to Peak Flow Rate max, otherwise Peak Flow Rate max remains unchanged.Calculate the difference (max-new) of Peak Flow Rate max and flow speed data new, relatively whether the size of the threshold value threshold2 of this difference and setting changes state three (flow velocity reduces) over to by state two (flow velocity increase) with the residing state of decision data, is divided into situation I and situation II.Wherein threshold value threshold2 can select according to Peak Flow Rate max, as gets max/3, max/2 etc.
The situation I: this difference adds 1 greater than threshold value threshold2 with integer variable decrease, judges then whether integer variable decrease has arrived designated value n2: if not yet arrive designated value n2, again change step 7 over to; If arrived designated value n2, judge three (flow velocity reduces) that got the hang of, therefore next begin to judge whether to get the hang of four (flow speed stabilities), the integer variable steady of several data stabilizations assignment has been arranged is 0 in order to be used for recording, flow speed data new assignment to flow speed data old, is entered step 8.Wherein the size of integer variable steady is used for characterizing the confidence level of got the hang of four (flow speed stabilities), and the larger confidence level of numerical value is higher, is judged to be four (flow speed stabilities) that got the hang of when increasing to designated value n3.
The situation II: this difference judges that whether integer variable decrease is greater than 1: if greater than 1, integer variable decrease is subtracted 1, reenter step 7 less than or equal to threshold value threshold2; If be not more than 1, judge three (flow velocity reduces) that do not get the hang of, be still namely state two (flow velocity increase), be therefore 0 with integer variable increase assignment, then return to step 6.
Step 8: read the new data of radar wave current meter and assignment to flow speed data new, and calculate the absolute value of the difference of flow speed data old and flow speed data new | old-new|.Relatively whether the size of the threshold value threshold3 of this absolute value and setting changes state four (flow speed stability) over to by state three (flow velocity reduces) with the residing state of decision data, is divided into situation a and situation b.
Situation a: this absolute value judges that whether integer variable steady is greater than 1: if greater than 1, integer variable steady is subtracted 1, reenter step 8 more than or equal to threshold value threshold3; If be not more than 1, judge four (flow speed stabilities) that do not get the hang of, namely still be state three (flow velocity reduces), be therefore 0 with integer variable decrease assignment, then return to step 7.
Situation b: this absolute value is less than threshold value threshold3, integer variable steady is added 1, then judge whether integer variable steady has arrived designated value n3: if not yet arrive designated value n3, with flow speed data new assignment to flow speed data old, again change step 8 over to; If arrived designated value n3, judge four (flow speed stabilities) that got the hang of, calculate the difference (new-original) of flow velocity data n ew and original flow velocity original, relatively the size of this difference and threshold value threshold4, be divided into situation b-1 and situation b-2.The flow velocity of flow speed data new recording status four (flow speed stability) wherein, the flow velocity of original flow velocity original recording status one (tidal bore not yet arrives).
Situation b-1: if this difference is greater than threshold value threshold4, show that causing the whole variation of this time flow speed data is not to be caused by interference, therefore judge that tidal bore arrives, and according to the rank size of the size judgement tidal bore of the difference (max-original) of Peak Flow Rate max and original flow velocity original, wherein max has recorded the Peak Flow Rate that reaches in whole change in flow process, and original flow velocity original has recorded the flow velocity of state one (tidal bore not yet arrives).
Situation b-2: if this difference is less than or equal to threshold value threshold4, show that causing the whole variation of this time flow speed data is to be caused by interference, therefore judge the state one (tidal bore not yet arrives) that is in, therefore with flow speed data new assignment to flow speed data old, read the new data of radar wave current meter and assignment to flow speed data new, return to step 3.
Be below the embodiment that takes according to the inventive method.
On August 21st, 2012 (lunar calendar the fifth day of a lunar month in July), near bank, Qiantang River crossing, No. 27, No. ten street of husky Higher Education Zone under the Jianggan District of Hangzhou (being positioned at approximately 3000 meters, bridge upstream, the east of a river), the radar wave current meter is installed, and is made radar over against the current inflow direction.
Threshold value Threshold1 gets 2, and threshold value threshold2 gets 2, and threshold value threshold3 gets 0.5, and threshold value threshold4 gets 0.5; Judgement enters the required designated value n1 that reaches of NextState, n2, n3 all get 15, and the resulting judged result of this method is as shown in the table when adopting respectively.
Claims (1)
1. one kind judges based on flow rate of water flow the method whether tidal bore arrives, and it is characterized in that the method comprises the following steps:
Step 1: select suitable place that the radar wave current meter is installed and measure flow rate of water flow, install on the bank iron shelf with the radar frame on river surface, make radar over against water (flow) direction;
Step 2: gather the instantaneous velocity data, two variablees are set are used for the flow speed data that record collects, be respectively flow speed data old and flow speed data new, to flow speed data old, the 2nd data assignment is to flow speed data new with the 1st data assignment collecting;
Step 3: calculate the difference of flow velocity data n ew and flow speed data old, relatively the size of the threshold value threshold1 of this difference and setting is to judge whether current state changes state two over to by state one; If this difference judges that less than or equal to the threshold value threshold1 that sets current state as state one, enters step 4; If this difference begins to judge whether to get the hang of two greater than the threshold value threshold1 that sets, give original flow velocity original with flow speed data old assignment, wherein the original flow velocity original of variable is used for the flow velocity of recording status one, enters step 5;
Step 4: flow speed data new assignment to flow speed data old, is then read the new data of radar wave current meter and assignment to flow speed data new, go back to step 3;
Step 5: for Peak Flow Rate max initialize is flow speed data new, wherein Peak Flow Rate max will be used for recording the maximum data that occurs in the tidal bore arrival process; Be 1 for being used for recording the integer variable increase assignment that has had several data to increase, expression has had 1 data to increase, wherein the size of integer variable increase can be used for characterizing two the confidence level of getting the hang of, the larger confidence level of numerical value is higher, when increasing to designated value n1, be judged to be and got the hang of two;
Step 6: read the new data of radar wave current meter and assignment to flow speed data new, calculate the difference of flow velocity data n ew and original flow velocity original, the size of the threshold value threshold1 of this difference and setting relatively, be divided into situation 1. with situation 2.;
Situation is 1.: this difference is greater than the threshold value threshold1 that sets, and judges that flow speed data new is whether greater than Peak Flow Rate max: if greater than Peak Flow Rate max, with its assignment to Peak Flow Rate max, otherwise Peak Flow Rate max remains unchanged; Integer variable increase is added 1, then judge whether integer variable increase has arrived designated value n1: if not yet arrive designated value n1, again change step 6 over to, if arrived designated value n1, judge and got the hang of two, therefore next begin to judge whether to get the hang of three, the integer variable decrease assignment that several data reduce has been arranged is 0 in order to be used for recording, and enters step 7; Wherein the size of integer variable decrease is used for characterizing three the confidence level of having got the hang of, and the larger confidence level of numerical value is higher, is judged to be when increasing to designated value n2 and has got the hang of three;
Situation is 2.: this difference judges that whether integer variable increase is greater than 1: if greater than 1, integer variable increase is subtracted 1, reenter step 6 less than or equal to the threshold value threshold1 that sets; If be not more than 1, judge not get the hang of two, namely be still state one, so integer variable increase assignment is 0, return to step 4;
Step 7: read the new data of radar wave current meter and assignment to flow speed data new, judge that flow speed data new is whether greater than Peak Flow Rate max: if greater than Peak Flow Rate max, with its assignment to Peak Flow Rate max, otherwise Peak Flow Rate max remains unchanged; Calculate the difference of Peak Flow Rate max and flow speed data new, relatively whether the size of the threshold value threshold2 of this difference and setting changes state three over to by state two with the residing state of decision data, is divided into situation I and situation II; Wherein threshold value threshold2 can select according to Peak Flow Rate max;
The situation I: this difference adds 1 greater than threshold value threshold2 with integer variable decrease, judges then whether integer variable decrease has arrived designated value n2: if not yet arrive designated value n2, again change step 7 over to; If arrived designated value n2, judge to have got the hang of three, therefore next begin to judge whether to get the hang of four, the integer variable steady of several data stabilizations assignment has been arranged is 0 in order to be used for recording, flow speed data new assignment to flow speed data old, is entered step 8; Wherein the size of integer variable steady is used for characterizing four the confidence level of having got the hang of, and the larger confidence level of numerical value is higher, is judged to be when increasing to designated value n3 and has got the hang of four;
The situation II: this difference judges that whether integer variable decrease is greater than 1: if greater than 1, integer variable decrease is subtracted 1, reenter step 7 less than or equal to threshold value threshold2; If be not more than 1, judge not get the hang of three, be still namely state two, be therefore 0 with integer variable increase assignment, then return to step 6;
Step 8: read the new data of radar wave current meter and assignment to flow speed data new, and calculate the absolute value of the difference of flow speed data old and flow speed data new; Relatively whether the size of the threshold value threshold3 of this absolute value and setting changes state four over to by state three with the residing state of decision data, is divided into situation a and situation b;
Situation a: this absolute value judges that whether integer variable steady is greater than 1: if greater than 1, integer variable steady is subtracted 1, reenter step 8 more than or equal to threshold value threshold3; If be not more than 1, judge not get the hang of four, be still namely state three, be therefore 0 with integer variable decrease assignment, then return to step 7;
Situation b: this absolute value is less than threshold value threshold3, integer variable steady is added 1, then judge whether integer variable steady has arrived designated value n3: if not yet arrive designated value n3, with flow speed data new assignment to flow speed data old, again change step 8 over to; If arrived designated value n3, judge to have got the hang of four, calculate the difference of flow velocity data n ew and original flow velocity original, relatively the size of this difference and threshold value threshold4, be divided into situation b-1 and situation b-2; The flow velocity of flow speed data new recording status four wherein, the flow velocity of original flow velocity original recording status one;
Situation b-1: if this difference is greater than threshold value threshold4, show that causing the whole variation of this time flow speed data is not to be caused by interference, therefore judge that tidal bore arrives, and according to the rank size of the size judgement tidal bore of the difference of Peak Flow Rate max and original flow velocity original, wherein max has recorded the Peak Flow Rate that reaches in whole change in flow process, and original flow velocity original has recorded the flow velocity of state one;
Situation b-2: if this difference is less than or equal to threshold value threshold4, show that causing the whole variation of this time flow speed data is to be caused by interference, therefore judge the state one that is in, therefore with flow speed data new assignment to flow speed data old, read the new data of radar wave current meter and assignment to flow speed data new, return to step 3;
Described state one not yet arrives for tidal bore, and state two increases for flow velocity, and state three is for flow velocity reduces, and state four is flow speed stability.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310019205.XA CN103090855B (en) | 2013-01-17 | 2013-01-17 | Method for determining arrival of tidal bore based on water velocity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310019205.XA CN103090855B (en) | 2013-01-17 | 2013-01-17 | Method for determining arrival of tidal bore based on water velocity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103090855A true CN103090855A (en) | 2013-05-08 |
| CN103090855B CN103090855B (en) | 2014-12-17 |
Family
ID=48203722
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310019205.XA Expired - Fee Related CN103090855B (en) | 2013-01-17 | 2013-01-17 | Method for determining arrival of tidal bore based on water velocity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103090855B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103424106A (en) * | 2013-08-27 | 2013-12-04 | 杭州市水文水资源监测总站 | River tidal bore automatic detection and real-time forecasting method based on water level and flow velocity |
| CN110388899A (en) * | 2019-06-04 | 2019-10-29 | 浙江省水利河口研究院 | The personal module calculated for tidal bore flow velocity vertical characteristics |
| CN110412312A (en) * | 2019-06-04 | 2019-11-05 | 浙江省水利河口研究院 | A kind of calculation method of tidal bore flow velocity vertical characteristics |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101403636A (en) * | 2008-10-21 | 2009-04-08 | 杭州电子科技大学 | River tidal surge detection method based on audio characteristic recognition |
| CN101625732A (en) * | 2009-08-03 | 2010-01-13 | 杭州电子科技大学 | Forecasting method of water level of potamic tidewater |
| CN201600384U (en) * | 2010-02-02 | 2010-10-06 | 杭州电子科技大学 | Non-contact tidal current measuring device |
| CN102156803A (en) * | 2011-03-18 | 2011-08-17 | 杭州电子科技大学 | Video recognition based river tidal bore detection method |
| CN202648661U (en) * | 2012-06-21 | 2013-01-02 | 杭州市水文水资源监测总站 | Flow velocity-based river eagre detection device |
| CN202676162U (en) * | 2012-06-21 | 2013-01-16 | 杭州电子科技大学 | Automatic river tidal bore real-time detection device with two radar wave flow meters |
-
2013
- 2013-01-17 CN CN201310019205.XA patent/CN103090855B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101403636A (en) * | 2008-10-21 | 2009-04-08 | 杭州电子科技大学 | River tidal surge detection method based on audio characteristic recognition |
| CN101625732A (en) * | 2009-08-03 | 2010-01-13 | 杭州电子科技大学 | Forecasting method of water level of potamic tidewater |
| CN201600384U (en) * | 2010-02-02 | 2010-10-06 | 杭州电子科技大学 | Non-contact tidal current measuring device |
| CN102156803A (en) * | 2011-03-18 | 2011-08-17 | 杭州电子科技大学 | Video recognition based river tidal bore detection method |
| CN202648661U (en) * | 2012-06-21 | 2013-01-02 | 杭州市水文水资源监测总站 | Flow velocity-based river eagre detection device |
| CN202676162U (en) * | 2012-06-21 | 2013-01-16 | 杭州电子科技大学 | Automatic river tidal bore real-time detection device with two radar wave flow meters |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103424106A (en) * | 2013-08-27 | 2013-12-04 | 杭州市水文水资源监测总站 | River tidal bore automatic detection and real-time forecasting method based on water level and flow velocity |
| CN103424106B (en) * | 2013-08-27 | 2015-06-03 | 杭州市水文水资源监测总站 | River tidal bore automatic detection and real-time forecasting method based on water level and flow velocity |
| CN110388899A (en) * | 2019-06-04 | 2019-10-29 | 浙江省水利河口研究院 | The personal module calculated for tidal bore flow velocity vertical characteristics |
| CN110412312A (en) * | 2019-06-04 | 2019-11-05 | 浙江省水利河口研究院 | A kind of calculation method of tidal bore flow velocity vertical characteristics |
| CN110412312B (en) * | 2019-06-04 | 2021-05-11 | 浙江省水利河口研究院 | Calculation method for tidal bore flow velocity vertical distribution |
| CN110388899B (en) * | 2019-06-04 | 2022-08-09 | 浙江省水利河口研究院 | Special assembly for calculating tidal bore flow velocity vertical distribution |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103090855B (en) | 2014-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Stockdon et al. | A simple model for the spatially-variable coastal response to hurricanes | |
| Douglas et al. | Long-term shoreline position prediction and error propagation | |
| Trung et al. | Vulnerability to flood in the Vietnamese Mekong Delta: mapping and uncertainty assessment | |
| Mel et al. | Addressing the effect of the Mo. SE barriers closure on wind setup within the Venice lagoon | |
| Miller et al. | A hydrological assessment of the November 2009 floods in Cumbria, UK | |
| Liao et al. | The cause of the 2008 cold disaster in the Taiwan Strait | |
| CN103090855B (en) | Method for determining arrival of tidal bore based on water velocity | |
| CN103546319B (en) | The alarming flow method and system of the network equipment | |
| Tseng et al. | Validation of the Kuroshio current system in the dual-domain Pacific Ocean model framework | |
| CN114510765B (en) | Dynamic dam-breaking submergence forecasting method for sea dam under typhoon storm surge action | |
| Kirkpatrick et al. | Modelling the effects of climate change on urban coastal-fluvial flooding | |
| CN202648661U (en) | Flow velocity-based river eagre detection device | |
| Zheng et al. | Numerical simulation of typhoon-induced storm surge along Jiangsu coast, Part I: Analysis of tropical cyclone | |
| Ko et al. | Impacts of tides and Typhoon Fanapi (2010) on seas around Taiwan | |
| CN101441078B (en) | River tidal bore subsection real time early warning method | |
| CN103424106A (en) | River tidal bore automatic detection and real-time forecasting method based on water level and flow velocity | |
| Hasegawa et al. | JMA’s storm surge prediction for the WMO Storm Surge Watch Scheme (SSWS) | |
| Yuk et al. | Modelling of storm-induced seawater flooding in the Suyeong River area, South Korea: A case study due to the storm surge and waves during Typhoon Sanba | |
| CN203204010U (en) | Eagre propagation speed measurement device based on step jump identification | |
| Waterhouse et al. | Influence of two tropical storms on the residual flow in a subtropical tidal inlet | |
| Tanaka et al. | Numerical model for predicting the seasonal migration of a river mouth | |
| Chaiongkarn et al. | Characteristics of Seasonal Wind and Wind-driven current in the Gulf of Thailand | |
| Lawless et al. | Development of the European tidal database and its potential application to marine renewables | |
| Kawai et al. | Uncertainty of extreme storm surge estimation by high wind sea surface drag coefficient and future typhoon change | |
| Chen et al. | Fast and minimally intrusive method for measuring tidal-stream discharge |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141217 Termination date: 20170117 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |