CN109557333B - Method for inverting surface ocean current based on navigation buoy data - Google Patents
Method for inverting surface ocean current based on navigation buoy data Download PDFInfo
- Publication number
- CN109557333B CN109557333B CN201811209727.5A CN201811209727A CN109557333B CN 109557333 B CN109557333 B CN 109557333B CN 201811209727 A CN201811209727 A CN 201811209727A CN 109557333 B CN109557333 B CN 109557333B
- Authority
- CN
- China
- Prior art keywords
- buoy
- ocean current
- parameters
- ocean
- navigation
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 230000033001 locomotion Effects 0.000 claims abstract description 25
- 238000005259 measurement Methods 0.000 claims abstract description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 230000007547 defect Effects 0.000 abstract description 2
- 238000000691 measurement method Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000004873 anchoring Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007227 biological adhesion Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 229930002875 chlorophyll Natural products 0.000 description 1
- 235000019804 chlorophyll Nutrition 0.000 description 1
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P13/00—Indicating or recording presence, absence, or direction, of movement
- G01P13/02—Indicating direction only, e.g. by weather vane
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/14—Receivers specially adapted for specific applications
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Aviation & Aerospace Engineering (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
Abstract
The invention discloses a method for inverting surface ocean current based on navigation buoy data, which comprises the following steps: step one, arranging a navigation buoy in an observation sea area, and carrying a sea current data measuring system; collecting buoy position data according to an ocean current data measuring system, and correcting to obtain motion parameters of the buoy; thirdly, collecting buoy position data according to the ocean current data measurement system to establish a buoy stress model; analyzing the parameters of the ocean current resistance of the buoy in the buoy stress model through an ocean current inversion algorithm; calculating the movement parameters of the buoy in the step two and the ocean current resistance parameters of the buoy in the step three to obtain initial values of the flow velocity and the direction of the ocean current in the observation sea area; and step six, correcting the initial values of the flow velocity and the direction of the ocean current obtained in the step five to obtain a final ocean current flow velocity value. The invention realizes accurate measurement of the ocean current in a specific sea area through the navigation buoy, and overcomes the defects of the traditional measurement method.
Description
Technical Field
The invention relates to the technical field of ocean current observation, in particular to a method for inverting surface ocean current based on navigation buoy data.
Background
Ocean currents have important influence and restriction on various physical, chemical, biological and geological effects and processes of salinity, chlorophyll and ocean heat transport in the ocean, as well as the formation and changes of climate and weather above the ocean, and have great social and economic significance on ocean development, navigation, military, port construction, offshore fishing and protection against oceanographic disasters for human activities. The traditional method for measuring the current flow mainly measures by fixing a current meter in the sea or performs inversion by using height satellite data, compared with the traditional measuring method, the method for inverting the current flow by using the navigation buoy can reduce the cost of purchasing instruments, avoid errors caused by attachment of marine organisms to the surfaces of the instruments, and can realize fixed-point current measurement.
Disclosure of Invention
The invention aims to overcome the defects in the prior art, and provides a method for inverting surface ocean currents based on navigation buoy data, which has the following specific technical scheme:
a method for inverting surface ocean currents based on navigation buoy data comprises the following steps:
step one, arranging a navigation buoy in an observation sea area, and carrying a sea current data measuring system;
step two, collecting a buoy position value and a wind speed value according to the ocean current data measurement system, and correcting to obtain a motion parameter of the buoy; meanwhile, collecting actual parameters of the quality, the underwater depth and the anchor chain length of the buoy according to an ocean current data measuring system;
step three, establishing a buoy stress model according to the buoy motion parameters and the actual parameters;
analyzing the parameters of the ocean current resistance of the buoy in the buoy stress model through an ocean current inversion algorithm;
calculating the movement parameters of the buoy in the step two and the ocean current resistance parameters of the buoy in the step three to obtain initial values of the flow velocity and the direction of the ocean current in the observation sea area;
and step six, correcting the initial values of the flow velocity and the direction of the ocean current obtained in the step five to obtain a final ocean current flow velocity value.
The ocean current data measuring system comprises a navigation buoy, a GPS locator, an anemoscope, a data transmission unit and an anchor chain; the upper portion of the buoy is connected with an anemoscope, the gravity center position inside the buoy is provided with a GPS locator and a data transmission unit respectively, and the lower portion of the buoy is connected with an anchor chain.
The ocean current inversion algorithm in the fourth step comprises the following steps:
3.1, calculating to obtain the movement acceleration of the buoy according to the buoy position value in the buoy stress model;
3.2, calculating according to the buoy wind speed value in the buoy stress model and the formula (1) to obtain the buoy wind stress;
R2＝kqAd (1)
3.3, calculating the anchor chain tension borne by the buoy according to the length of the anchor chain of the buoy in the buoy stress model and a formula (2);
3.4, calculating the gravity borne by the buoy according to the mass of the buoy in the buoy stress model;
3.5, calculating buoyancy borne by the buoy according to the water entry depth value of the buoy in the buoy stress model;
3.6, obtaining an ocean current resistance parameter R1 borne by the buoy according to the known parameters of the buoy generated in the step 3.1 to 3.5 through a buoy motion N-S control equation established by a Newton second motion law;
and 3.7, obtaining the ocean current flow velocity V according to the formula (3) according to the known ocean current resistance parameter R1.
R1＝C(ρV^2)Ω/2 (3)
Advantageous effects
Compared with the prior art, the technical scheme of the invention has the following beneficial effects:
1. compared with the current traditional method for acquiring the ocean current flow velocity through an ocean current meter and a satellite data inversion method, the method for inverting the ocean current by installing the GPS locator in the navigation buoy and combining the buoy stress model is obviously capable of greatly reducing the cost for acquiring the ocean current flow velocity, and compared with the satellite inversion method, the method is simple in algorithm and small in data size.
2. The invention can effectively solve the measurement error caused by the problems of biological adhesion and the like. In addition, the anchor point of the buoy can be changed at any time according to the measurement requirement, so that the size and the direction of the ocean current in a specific sea area can be measured at any time and any place.
Drawings
FIG. 1 is a schematic view of instrument installation of a navigation buoy measurement system;
FIG. 2 is a flow chart of the ocean current inversion algorithm.
Reference numerals:
1-buoy 2-GPS locator 3-data transmission equipment 4-anemoscope 5-anchor chain
Detailed Description
The technical solution of the present patent will be described in further detail with reference to the following embodiments.
The invention introduces a method for inverting surface ocean current based on navigation buoy data, which is a method for acquiring data by carrying a high-precision GPS locator on an ocean buoy and analyzing the flow velocity and direction of ocean current by data, and the specific process is as follows:
step one, a navigation buoy 1 is arranged in an observation sea area and fixed at an anchor point by an anchor chain 5. The buoy is loaded with a high-precision GPS locator 2, an anemoscope 4 and data transmission equipment 3. All instruments are rigidly connected to the buoy 1, mounted and positioned as shown in fig. 1. The GPS locator 2 is used for collecting position parameters of the buoy, including information such as longitude, latitude and elevation of the buoy; the anemoscope 4 is used for measuring the wind power of the buoy on the sea; the data transfer device 3 is used to transfer data to the shore base.
And step two, collecting buoy position data. The GPS locator 2 is used for acquiring longitude and latitude information of the buoy and correcting the posture of the buoy to obtain a buoy movement track, and accordingly, movement parameters such as acceleration and the like of the buoy are acquired; and further establishing a buoy stress model, wherein the establishing method comprises the following steps:
buoys in water are mainly subject to the action of factors such as ocean current resistance, wind load, and the like, and in addition, are subject to the tension of the mooring cables.
At present, a formula for calculating the frictional resistance of ocean currents is more, and the resistance calculation formula is commonly used in China and is obtained based on flat plate tests and theoretical analysis, and the expression is as follows:
R1＝C(ρV^2)Ω/2 (3)
wherein V is the flow velocity of the seawater.
The wind pressure acting on the buoy mark body can be determined by a calculation formula given in Port engineering calculation Specifications:
R2＝kqAd (1)
according to the balance analysis of the smart force of the buoy of the H.0. Betto and the rock well, the cable is in a catenary curve shape, and the length and the force of the cable can be determined by the following calculation formula:
the force analysis of the buoy shows that the water flow force borne by the buoy has a certain correlation with the ocean current flow rate, the water flow force can be obtained through a buoy motion equation according to other known forces, then the ocean current flow rate is obtained by referring to an empirical formula and is corrected by using experimental data. The direction of the flow velocity may also be determined by the direction of the float's trajectory. The core idea of the method is to invert the flow velocity of the ocean current through the motion track of the anchoring buoy.
Step four, combining the winning motion parameters in step two and the algorithm of calculating the ocean current by the ocean current resistance parameters of the buoy in step three to obtain the initial values of the flow speed and the direction of the ocean current in the sea area; as shown in the flowchart of fig. 2, the motion trajectory of the buoy is obtained by correcting the data acquired by the buoy and the posture of the buoy, and the motion trajectory is embodied by the force received by the buoy; then combining the known forces obtained by stress analysis, namely buoyancy, wind power, anchor chain restraining force, gravity and buoyancy, the ocean current resistance borne by the buoy can be obtained by resolving according to Newton' S law of motion and motion N-S control equation, and finally the ocean current flow velocity is obtained according to the empirical formula of the ocean current resistance and the ocean current magnitude, namely R1 ═ C (rho V ^2)/2 omega. The specific method comprises the following steps:
a. and a GPS positioning instrument is used for obtaining the position value of the buoy, and an anemoscope is used for obtaining the wind speed value borne by the buoy. And obtaining the numerical values of the buoy, such as the mass, the water penetration depth, the anchor chain length and the like according to actual measurement.
b. The position value is used to obtain the moving acceleration of buoy, the wind speed value is used to obtain the wind stress applied to buoy according to the wind pressure formula R2-kqAd, and the anchor chain length is used to obtain the anchor chain tension formula according to H.0 Betto and rock well intelligence
Calculating the tension of an anchor chain borne by the buoy, calculating the gravity borne by the buoy by using the mass of the buoy, and calculating the buoyancy borne by the buoy by using the depth of water;
c. and establishing an N-S control equation of the buoy motion according to the Newton second motion law, and solving the N-S control equation according to the known acceleration, wind stress, anchor chain tension, gravity and buoyancy to obtain the water flow tension R1 borne by the buoy.
d. And solving the ocean current flow velocity V according to the obtained current drag force R1 by using an ocean current resistance calculation formula R1 ═ C (rho V ^2) omega/2 obtained based on flat plate test and theoretical analysis. Thus, the initial values of the magnitude and direction of the current are obtained.
And step five, correcting the initial values of the flow velocity and the direction of the ocean current obtained in the step four to obtain the final result. Since the result obtained in step four is obtained by an empirical formula, it is necessary to repeatedly perform experiments to correct the relevant parameters to obtain more accurate results.
The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for a person skilled in the art, several modifications and variations can be made without departing from the principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.
Claims (2)
1. A method for inverting surface ocean currents based on navigation buoy data comprises the following steps:
step one, arranging a navigation buoy in an observation sea area, and carrying a sea current data measuring system;
step two, collecting a buoy position value and a wind speed value according to the ocean current data measurement system, and correcting to obtain a motion parameter of the buoy; meanwhile, collecting actual parameters of the quality, the underwater depth and the anchor chain length of the buoy according to an ocean current data measuring system;
step three, establishing a buoy stress model according to the buoy motion parameters and the actual parameters;
analyzing the parameters of the ocean current resistance borne by the buoy in the buoy stress model through an ocean current inversion algorithm;
the ocean current inversion algorithm in the fourth step comprises the following steps:
3.1, calculating to obtain the movement acceleration of the buoy according to the buoy position value in the buoy stress model;
3.2, calculating according to the buoy wind speed value in the buoy stress model and the formula (1) to obtain the buoy wind stress;
R2＝kqAd (1)
3.3, calculating the anchor chain tension borne by the buoy according to the length of the anchor chain of the buoy in the buoy stress model and a formula (2);
3.4, calculating the gravity borne by the buoy according to the mass of the buoy in the buoy stress model;
3.5, calculating buoyancy borne by the buoy according to the water entry depth of the buoy in the buoy stress model;
3.6, obtaining an ocean current resistance parameter R1 borne by the buoy according to the known parameters of the buoy generated in the step 3.1 to 3.5 through a buoy motion N-S control equation established by a Newton second motion law;
3.7 obtaining the current flow velocity V according to the formula (3) according to the known current resistance parameter R1:
R1＝C(ρV^2)Ω/2 (3)
step five, calculating the movement parameters of the buoy in the step two and the ocean current resistance parameters borne by the buoy in the step four to obtain initial values of the flow velocity and the direction of the ocean current in the observation sea area;
and step six, correcting the initial values of the flow velocity and the direction of the ocean current obtained in the step five to obtain a final ocean current flow velocity value.
2. The method for inverting surface ocean currents based on navigation buoy data according to claim 1, wherein the ocean current data measuring system comprises a navigation buoy, a GPS locator, an anemometer, a data transmission unit and an anchor chain; the upper portion of the buoy is connected with an anemoscope, the gravity center position inside the buoy is provided with a GPS locator and a data transmission unit respectively, and the lower portion of the buoy is connected with an anchor chain.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811209727.5A CN109557333B (en) | 2018-10-17 | 2018-10-17 | Method for inverting surface ocean current based on navigation buoy data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811209727.5A CN109557333B (en) | 2018-10-17 | 2018-10-17 | Method for inverting surface ocean current based on navigation buoy data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109557333A CN109557333A (en) | 2019-04-02 |
| CN109557333B true CN109557333B (en) | 2022-04-15 |
Family
ID=65865070
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811209727.5A Expired - Fee Related CN109557333B (en) | 2018-10-17 | 2018-10-17 | Method for inverting surface ocean current based on navigation buoy data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109557333B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112433241A (en) * | 2020-11-09 | 2021-03-02 | 天津大学 | Submersible vehicle positioning method based on buoy |
| CN113405537A (en) * | 2021-07-20 | 2021-09-17 | 中国海洋大学 | Wave direction inversion method based on satellite navigation positioning |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101813476A (en) * | 2010-03-19 | 2010-08-25 | 天津大学 | Three-dimensional real-time monitoring system for offshore wave parameters |
| KR101360296B1 (en) * | 2013-11-15 | 2014-02-11 | 주식회사 범아엔지니어링 | Survey appratus of measurement system for visualizing 3d shape with measured data in the ocean |
| CN106940382A (en) * | 2017-03-22 | 2017-07-11 | 国家海洋局第二海洋研究所 | A kind of ocean current tester |
| CN107748360A (en) * | 2017-09-05 | 2018-03-02 | 浙江海洋大学 | Extra large table Wind-field Retrieval method and device |
-
2018
- 2018-10-17 CN CN201811209727.5A patent/CN109557333B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101813476A (en) * | 2010-03-19 | 2010-08-25 | 天津大学 | Three-dimensional real-time monitoring system for offshore wave parameters |
| KR101360296B1 (en) * | 2013-11-15 | 2014-02-11 | 주식회사 범아엔지니어링 | Survey appratus of measurement system for visualizing 3d shape with measured data in the ocean |
| CN106940382A (en) * | 2017-03-22 | 2017-07-11 | 国家海洋局第二海洋研究所 | A kind of ocean current tester |
| CN107748360A (en) * | 2017-09-05 | 2018-03-02 | 浙江海洋大学 | Extra large table Wind-field Retrieval method and device |
Non-Patent Citations (1)
| Title |
|---|
| 海流测量技术发展及应用;周庆伟 等;《海洋测绘》;20180531;第38卷(第3期);73-77 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109557333A (en) | 2019-04-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103017755B (en) | A kind of underwater navigation attitude measurement method | |
| McMillan et al. | Rates of dissipation of turbulent kinetic energy in a high Reynolds number tidal channel | |
| CN107585263B (en) | One kind is orientable to be tethered at test platform from sink-float deep-sea | |
| St. Laurent et al. | Measurements of near-surface turbulence and mixing from autonomous ocean gliders | |
| CN109557333B (en) | Method for inverting surface ocean current based on navigation buoy data | |
| CN110488334A (en) | A kind of underwater object locator, positioning system and its localization method | |
| CN104613906A (en) | Sound ray tracking based reservoir area deep water depth measurement method | |
| CN107462891B (en) | Three-point type deep sea submerged buoy positioning method | |
| CN109059746A (en) | A kind of bathymetric surveying method based on accurate POS | |
| Long et al. | Surface gas exchange determined from an aquatic eddy covariance floating platform | |
| RU2561229C1 (en) | Buoy for determination of characteristics of sea wind waves | |
| Feddersen et al. | The wavedrifter: a low-cost IMU-based Lagrangian drifter to observe steepening and overturning of surface gravity waves and the transition to turbulence | |
| Toma et al. | Performance of WAVY Ocean Lagrangian Drifters for Surface Characterization of Ocean Dynamic Structures | |
| CN109398615B (en) | Deep sea laying method of sinking block for buoy | |
| KR101301743B1 (en) | Oceanographic Observation Station for Height and Direction Measuring of a Wave | |
| CN110222304B (en) | Method for calculating amount of motion of floating structure based on mooring force | |
| Richardson | Drifters and floats | |
| CN102540257B (en) | Positioning method of earthquake signal receiving device | |
| McMillan | Turbulence Measurements in a High Reynolds Number Tidal Channel | |
| Prater et al. | Observations of the Faroe Bank Channel overflow using bottom-following RAFOS floats | |
| Sheng et al. | Wave-induced sediment resuspension and mixing in shallow waters | |
| Johnson | Measuring the orientation and movement of marine animals using inertial and magnetic sensors-a tutorial | |
| Chenbing et al. | A small buoy for flux measurement in air-sea boundary layer | |
| Lueck et al. | Hub-height time series measurements of velocity and dissipation of turbulence kinetic energy in a tidal channel | |
| McCullough | Problems in measuring currents near the ocean surface |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220415 |