CN109557333A - A method of based on navigation buoy data inversion sea surface current - Google Patents

A method of based on navigation buoy data inversion sea surface current Download PDF

Info

Publication number
CN109557333A
CN109557333A CN201811209727.5A CN201811209727A CN109557333A CN 109557333 A CN109557333 A CN 109557333A CN 201811209727 A CN201811209727 A CN 201811209727A CN 109557333 A CN109557333 A CN 109557333A
Authority
CN
China
Prior art keywords
buoy
ocean current
parameter
current
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.)
Granted
Application number
CN201811209727.5A
Other languages
Chinese (zh)
Other versions
CN109557333B (en
Inventor
张安民
宁伟
宁一伟
张豪
欧阳知雨
刘禹
刘帅
刘荣霞
侯泽北
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201811209727.5A priority Critical patent/CN109557333B/en
Publication of CN109557333A publication Critical patent/CN109557333A/en
Application granted granted Critical
Publication of CN109557333B publication Critical patent/CN109557333B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P13/00Indicating or recording presence, absence, or direction, of movement
    • G01P13/02Indicating direction only, e.g. by weather vane
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/14Receivers 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 kind of methods based on navigation buoy data inversion sea surface current, comprising the following steps: Step 1: laying navigation buoy in observation sea area, carries ocean current data measurin system;Its kinematic parameter is obtained Step 2: collecting buoy position data according to ocean current data measurin system and being modified;Step 3: collecting buoy position data according to ocean current data measurin system establishes buoy stress model;Step 4: by buoy in ocean current inversion algorithm analysis buoy stress model by current drag parameter;Step 5: buoy in buoy movement parameter in step 2 and step 3 is carried out the initial value that operation obtains observation sea area ocean current flow velocity size and Orientation by current drag parameter;Final ocean current flow velocity numerical value is obtained Step 6: being modified to the initial value of the ocean current flow velocity size and Orientation obtained in step 5.The present invention accurately measures the ocean current of specified sea areas by navigation buoy realization, improves the deficiency of traditional measurement method.

Description

A method of based on navigation buoy data inversion sea surface current
Technical field
The present invention relates to Current Observations technical fields, are based on navigation buoy data inversion sea surface current more particularly to one kind Method.
Background technique
Ocean current is to physics a variety of in ocean, chemistry, biology and geologic process and conveying salinity, chlorophyll and ocean heat mistake The formation and variation of the climate and weather in journey and ocean overhead, play the role of important influence and restriction, simultaneously for people For class activity, ocean current all has ocean development, navigation, military affairs, port construction, marine fishing and defence maritime meteorology disaster There are great social and economic implications.The method of traditional sea flow velocity measuring in marine fixed current meter mainly by being surveyed Amount can both have been subtracted compared to traditional measurement method using navigation buoy inverting ocean current using the progress inverting of altitude satellite data The cost for buying instrument less also can avoid error caused by surface attachment of the marine organisms to instrument, and can be realized fixed point and survey Measure ocean current.
Summary of the invention
Purpose of the invention is to overcome the shortcomings in the prior art, proposes a kind of based on navigation buoy data inversion The method of sea surface current, specific technical solution are as follows:
A method of based on navigation buoy data inversion sea surface current, comprising the following steps:
Step 1: laying navigation buoy in observation sea area, ocean current data measurin system is carried;
It is obtained Step 2: being modified according to ocean current data measurin system collection buoy position value and air speed value
Kinematic parameter;Meanwhile quality, underwater penetration and the anchor of buoy are collected according to ocean current data measurin system
The actual parameter of chain length;
Step 3: establishing buoy stress model according to buoy movement parameter and actual parameter;
Step 4: by buoy in ocean current inversion algorithm analysis buoy stress model by current drag parameter;
Step 5: buoy in buoy movement parameter in step 2 and step 3 is carried out operation acquisition by current drag parameter Observe the initial value of sea area ocean current flow velocity size and Orientation;
Final ocean current is obtained Step 6: being modified to the initial value of the ocean current flow velocity size and Orientation obtained in step 5 Flow velocity numerical value.
The ocean current data measurin system includes navigation buoy, GPS positioning instrument, anemobiagraph, data transmission unit and anchor chain; The buoy top is connected with anemobiagraph, and position of centre of gravity is respectively arranged with GPS positioning instrument inside the buoy and data transmission is single Member, the buoy lower part are connected with anchor chain.
Ocean current inversion algorithm described in step 4 includes the following steps:
3.1, the acceleration of motion for obtaining buoy is calculated according to buoy position value in buoy stress model;
3.2, according to buoy air speed value in buoy stress model according to formula (1) calculate obtain buoy by wind-stress;
R2=kqAd (1)
3.3, anchor chain pulling force suffered by buoy is calculated according to formula (2) according to buoy mooring line length in buoy stress model;
3.4, the gravity according to suffered by the Mass Calculation buoy of buoy in buoy stress model;
3.5, the buoyancy according to suffered by the underwater penetration value calculating buoy of buoy in buoy stress model;
3.6, the parameter according to known to the buoy that step 3.1~3.5 generate is established floating by Newton interpolation algorithm Mark movement N-S governing equation obtains current drag parameter R1 suffered by buoy;
3.7, ocean current flow velocity V is obtained according to formula (3) according to known current drag parameter R1.
R1=C (ρ V^2)/2 Ω (3)
Beneficial effect
Compared with prior art, the beneficial effects brought by the technical solution of the present invention are as follows:
1, the inverting obtained by installing GPS positioning instrument combination buoy stress model in navigation buoy that the present invention uses The method of ocean current with it is traditional at present ocean current flow velocity is obtained by current meter and satellite data inversion method compared with, it is clear that the party Method can substantially reduce the cost for obtaining ocean current flow velocity, and the data for simply needing to use compared to satellite Retrieval this method algorithm Measure it is small,.
2, the present invention can effectively solve the measurement error occurred due to the problems such as biological attachment.In addition to this, the anchor of buoy Fixed point can need to change at any time according to measurement, to realize ocean current size and the side for measuring specified sea areas whenever and wherever possible To.
Detailed description of the invention
Fig. 1 is navigation buoy measuring system instrument scheme of installation;
Fig. 2 is the resolving flow chart of the ocean current inversion algorithm.
Appended drawing reference:
1- buoy 2-GPS position indicator 3- data transmission set 4- anemobiagraph 5- anchor chain
Specific embodiment
The technical solution of the patent is explained in further detail With reference to embodiment.
Invention describes a kind of methods based on navigation buoy data inversion sea surface current, refer to by oceanographic buoy The upper method for carrying high-precision GPS position indicator and obtaining data and analyzing ocean current flow velocity size and Orientation by data, detailed process It is as follows:
Step 1: laying navigation buoy 1 in observation sea area, and anchor point is fixed in anchor chain 5.Wherein carried in buoy High-precision GPS position indicator 2, anemobiagraph 4 and data transmission set 3.All appts are all rigidly connected with buoy 1, installation Mode is as shown in Figure 1 with position.GPS positioning instrument 2 is used to acquire the location parameter of buoy, longitude, latitude and height including buoy The information such as journey;Anemobiagraph 4 is for measuring the wind-force size that buoy is at sea subject to;Data transmission set 3 is for transmitting data To bank base.
Step 2: collecting buoy position data.The latitude and longitude information of buoy is obtained by GPS positioning instrument 2 and corrects buoy Buoy movement track is obtained after posture, obtains the kinematic parameter of buoy, such as acceleration accordingly;And then establish buoy stress mould Type, method for building up are as follows:
In addition to this buoy in water, is also tethered at cable mainly by the effect of the factors such as current drag, wind load The pulling force of rope.
Currently, the formula for calculating ocean current frictional resistance is more, the most commonly used is be based on treadmill test and theory analysis institute in China Obtained resistance calculation formulae, expression formula are as follows:
R1=C (ρ V^2)/2 Ω (3)
Wherein, V is seawater velocity.
The wind pressure acted in buoy standard type can be determined as the calculation formula given by " port engineering calculates specification ":
R2=kqAd (1)
According to the H.0. equilibrium analysis of shellfish support and rock well acute hearing to buoy stress, cable is in catenary curve shape, and its length It can be determined by following calculating formula with stress:
By the force analysis to buoy it is found that flow force and ocean current flow velocity that buoy is subject to have certain correlativity, Flow force can be found out according to the other known buoy movement equation that is forced through, obtain ocean current flow velocity referring next to empirical equation Size and be modified with experimental data.The direction of flow velocity can also be determined by the motion profile direction of buoy.The core of method Thought is exactly the flow velocity by the motion profile inverting ocean current of mooring buoys.
Step 4: in conjunction with buoy in step 2 acceptance of the bid kinematic parameter and step 3 by the calculation of current drag coaptation ocean current Method obtains the initial value of the sea area ocean current flow velocity size and Orientation;As shown in Fig. 2 flow chart, pass through the collected number of buoy first According to obtain the motion profile of buoy after buoy attitude rectification, which is embodiment strong suffered by buoy;So Buoyancy, wind-force and anchor chain restraining force, gravity, the buoyancy that the known force obtained afterwards in conjunction with force analysis has buoy to be subject to, according to Newton's laws of motion and movement N-S governing equation can resolve the current drag for showing that buoy is subject to, last according to current drag Ocean current flow velocity is obtained with empirical equation R1=C (ρ V^2)/2 Ω of ocean current size.The specific method is as follows:
A. the positional value that buoy is obtained with GPS positioning instrument obtains air speed value suffered by buoy with anemobiagraph.According to actual measurement Obtain the numerical value such as quality, underwater penetration, the length of chain cable of buoy.
B. the acceleration of motion that buoy is acquired with positional value acquires buoy according to wind pressure formula R2=kqAd with air speed value Suffered wind-stress, the anchor chain pulling force formula provided with the length of chain cable of buoy according to H.0. shellfish support and rock well acute hearingAnchor chain pulling force suffered by buoy is calculated, gravity suffered by the Mass Calculation buoy with buoy, with entering water Buoyancy suffered by depth calculation buoy;
C. buoy movement N-S governing equation is established according to Newton interpolation algorithm, is answered according to known acceleration, wind Power, anchor chain pulling force, gravity, buoyancy acquire water flow pulling force R1 suffered by buoy by solving N-S governing equation.
D. it is utilized according to the water flow pulling force R1 acquired and is calculated based on treadmill test and the obtained current drag of theory analysis Formula R1=C (ρ V^2)/2 Ω solve ocean current flow velocity V.Thus the initial value of the size and Orientation of ocean current flow velocity is obtained.
Step 5: being modified to the initial value of the ocean current flow velocity size and Orientation obtained in step 4, obtain to the end As a result.Due to obtained in step 4 the result is that being obtained by empirical equation, it is therefore necessary to repeat experiment and join to correlation Number is corrected to obtain more accurate result.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art For member, without departing from the principle of the present invention, several improvements and changes can also be made, such modifications and variations are also answered It is considered as the protection scope of invention.

Claims (3)

1. a kind of method based on navigation buoy data inversion sea surface current, comprising the following steps:
Step 1: laying navigation buoy in observation sea area, ocean current data measurin system is carried;
Its kinematic parameter is obtained Step 2: being modified according to ocean current data measurin system collection buoy position value and air speed value; Meanwhile the quality of buoy, the actual parameter of underwater penetration and length of chain cable are collected according to ocean current data measurin system;
Step 3: establishing buoy stress model according to buoy movement parameter and actual parameter;
Step 4: by buoy in ocean current inversion algorithm analysis buoy stress model by current drag parameter;
It is observed Step 5: buoy in buoy movement parameter in step 2 and step 3 is carried out operation by current drag parameter The initial value of sea area ocean current flow velocity size and Orientation;
Final ocean current flow velocity is obtained Step 6: being modified to the initial value of the ocean current flow velocity size and Orientation obtained in step 5 Numerical value.
2. according to claim 1 be based on navigation buoy data inversion sea surface current method, which is characterized in that the ocean current Data measurin system includes navigation buoy, GPS positioning instrument, anemobiagraph, data transmission unit and anchor chain;The buoy top connection There is anemobiagraph, position of centre of gravity is respectively arranged with GPS positioning instrument and data transmission unit inside the buoy, and the buoy lower part connects It is connected to anchor chain.
3. the method for navigation buoy data inversion sea surface current according to claim 1, which is characterized in that
Ocean current inversion algorithm described in step 4 includes the following steps:
3.1, the acceleration of motion for obtaining buoy is calculated according to buoy position value in buoy stress model;
3.2, according to buoy air speed value in buoy stress model according to formula (1) calculate obtain buoy by wind-stress;
R2=kqAd (1)
3.3, anchor chain pulling force suffered by buoy is calculated according to formula (2) according to buoy mooring line length in buoy stress model;
3.4, the gravity according to suffered by the Mass Calculation buoy of buoy in buoy stress model;
3.5, the buoyancy according to suffered by the underwater penetration value calculating buoy of buoy in buoy stress model;
3.6, the parameter according to known to the buoy that step 3.1~3.5 generate is transported by the buoy that Newton interpolation algorithm is established Dynamic N-S governing equation obtains current drag parameter R1 suffered by buoy;
3.7, ocean current flow velocity V is obtained according to formula (3) according to known current drag parameter R1;
R1=C (ρ V^2)/2 Ω (3).
CN201811209727.5A 2018-10-17 2018-10-17 Method for inverting surface ocean current based on navigation buoy data Expired - Fee Related CN109557333B (en)

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 true CN109557333A (en) 2019-04-02
CN109557333B 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)

Cited By (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
周庆伟 等: "海流测量技术发展及应用", 《海洋测绘》 *

Cited By (2)

* Cited by examiner, † Cited by third party
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

Also Published As

Publication number Publication date
CN109557333B (en) 2022-04-15

Similar Documents

Publication Publication Date Title
CN108362281B (en) Long-baseline underwater submarine matching navigation method and system
Stevens et al. Linking glacially modified waters to catchment-scale subglacial discharge using autonomous underwater vehicle observations
CN107585263B (en) One kind is orientable to be tethered at test platform from sink-float deep-sea
CN203785669U (en) High-precision deep sea depth detection device
CN106768187B (en) Fixing device for sea-way topography measuring water level self-recording instrument of multi-tide station
CN105738869B (en) A kind of deep water beacon search localization method suitable for single hydrophone
CN104613906A (en) Sound ray tracking based reservoir area deep water depth measurement method
CN109557333A (en) A method of based on navigation buoy data inversion sea surface current
CN104280024B (en) Device and method for integrated navigation of deepwater robot
CN107462891B (en) Three-point type deep sea submerged buoy positioning method
Caldwell Small‐scale physics of the ocean
Toma et al. Performance of WAVY Ocean Lagrangian Drifters for Surface Characterization of Ocean Dynamic Structures
Roiha et al. Estimating currents from argo trajectories in the bothnian sea, baltic sea
CN111427010A (en) ROV underwater positioning system and positioning method
KR101301743B1 (en) Oceanographic Observation Station for Height and Direction Measuring of a Wave
Prater et al. Observations of the Faroe Bank Channel overflow using bottom-following RAFOS floats
Shitashima et al. Strategies for detection and monitoring of CO2 leakage in sub-seabed CCS
Purba et al. Performance of lagriangan drifter oceanography coverage area (RHEA): second phase
Menard et al. Underwater mapping by diving geologists
Richardson Drifters and floats
Alexander et al. Multi-modal sensing drifters as a tool for repeatable glacial hydrology flow path measurements
CN111427011A (en) Submarine asset position calibration method and system
Shi Underwater Terrain Measurement Based on Unmanned Measurement Ship Multi-Beam Bathymeter
CN103791900A (en) Method for determining position of towed body under water
CN104677414B (en) CTD data processing methods based on AUV

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