[go: up one dir, main page]

CN115327669B - 海面风场反演方法及装置 - Google Patents

海面风场反演方法及装置 Download PDF

Info

Publication number
CN115327669B
CN115327669B CN202211002563.5A CN202211002563A CN115327669B CN 115327669 B CN115327669 B CN 115327669B CN 202211002563 A CN202211002563 A CN 202211002563A CN 115327669 B CN115327669 B CN 115327669B
Authority
CN
China
Prior art keywords
sequence
spectrum
parameter
buoy
wave
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.)
Active
Application number
CN202211002563.5A
Other languages
English (en)
Other versions
CN115327669A (zh
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.)
National Ocean Technology Center
Original Assignee
National Ocean Technology Center
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 National Ocean Technology Center filed Critical National Ocean Technology Center
Priority to CN202211002563.5A priority Critical patent/CN115327669B/zh
Publication of CN115327669A publication Critical patent/CN115327669A/zh
Application granted granted Critical
Publication of CN115327669B publication Critical patent/CN115327669B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/02—Instruments for indicating weather conditions by measuring two or more variables, e.g. humidity, pressure, temperature, cloud cover or wind speed
    • G01W1/04—Instruments for indicating weather conditions by measuring two or more variables, e.g. humidity, pressure, temperature, cloud cover or wind speed giving only separate indications of the variables measured
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63B—SHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B22/00—Buoys
    • 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
    • G01P5/001—Full-field flow measurement, e.g. determining flow velocity and direction in a whole region at the same time, flow visualisation
    • 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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/52—Determining velocity
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Ocean & Marine Engineering (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Chemical & Material Sciences (AREA)
  • Mechanical Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Combustion & Propulsion (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)
  • Other Liquid Machine Or Engine Such As Wave Power Use (AREA)

Abstract

本申请提供一种海面风场反演方法及装置，涉及海洋波浪测量技术领域，该方法获取目标海面中浮标的运动参数序列，根据浮标的运动参数序列，确定浮标在笛卡尔坐标系的第二横向加速度序列、第二纵向加速度序列、第二竖直加速度序列，根据第二竖直加速度序列，确定目标海面的海浪能量谱，根据海浪能量谱确定目标海面的风速，根据第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列，确定目标海面的海浪方向谱，根据海浪方向谱确定目标海面的风向。本申请根据浮标的运动参数序列得到了目标海面的海浪方向谱和能量谱，再根据海浪方向谱和能量谱确定海面的风速和方向，实现了对海浪、风要素小范围、高时空分辨率的观测。

Description

海面风场反演方法及装置
技术领域
本申请涉及海洋波浪测量技术领域，具体涉及一种海面风场反演方法及装置。
背景技术
海面的风要素测量对于海洋工程研究、海洋气象预报、海洋灾害预警等有着重要意义。传统的海面风要素测量一般采用多点锚系浮标组网、卫星遥感、船基观测等方式。但是锚系浮标造价和维护成本高，多点锚系浮标组网测量方式需要耗费大量的成本。卫星遥感测量方式对海浪、风要素的获取是通过间接观测得到的，其数据观测精度和数据质量高度依赖于定标、现场真实性检验以及经验模型的选取。因此，传统的卫星遥感测量方式得到的风要素数据准确性低。传统的船基观测方式主要依靠船舶航行过程中获取海浪、风要素数据。然而，各类船舶航线相对比较固定，且大多数航线贴近大陆边缘，观测范围难以有效覆盖全球深海大洋，只能对区域性的海洋环境进行基本描述，因此，船基观测方式观测范围小、获取的海浪、风要素数据有限。
综上，传统采用多点锚系浮标组网、卫星遥感、船基观测等风速测量方式具有成本高昂、观测要素单一、观测准确性低、时空分辨率不足等问题，无法实现对海浪、风要素小范围、高时空分辨率的观测。
发明内容
本申请实施例提供一种海面风场反演方法及装置，实现了对海浪、风要素小范围、高时空分辨率的观测。
一方面，本申请实施例提供了一种海面风场反演方法，包括：
获取目标海面中浮标的运动参数序列，运动参数序列包括浮标角度参数序列和浮标在东北天坐标系的东北天坐标参数序列，浮标角度参数序列包括横摇角参数序列、纵摇角参数序列、艏摇角参数序列，东北天坐标参数序列包括第一横向加速度序列、第一纵向加速度序列、第一竖直加速度序列；
根据浮标的运动参数序列，确定浮标在笛卡尔坐标系的笛卡尔坐标参数序列，笛卡尔坐标参数序列包括第二横向加速度序列、第二纵向加速度序列、第二竖直加速度序列；
根据第二竖直加速度序列，确定目标海面的海浪能量谱；
根据海浪能量谱确定目标海面的风速；
根据第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列，确定目标海面的海浪方向谱；
根据海浪方向谱确定目标海面的风向。
可选地，在本申请一些可能的实现方式中，根据第二竖直加速度序列，确定目标海面的海浪能量谱，包括：
对第二竖直加速度序列进行积分，得到浮标在竖直方向上的位移序列；
将浮标在竖直方向上的位移序列分段，得到多个位移子序列；
确定多个位移子序列的能量谱密度；
根据多个位移子序列的能量谱密度，确定目标海面的海浪能量谱。
可选地，在本申请一些可能的实现方式中，根据第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列，确定目标海面的海浪方向谱，包括：
对第二横向加速度序列进行积分，得到浮标在横向上的位移序列；
对第二纵向加速度序列进行积分，得到浮标在纵向上的位移序列；
对第二竖直加速度序列进行积分，得到浮标在竖直方向的位移序列；
根据浮标在横向上的位移序列、浮标在纵向上的位移序列，以及浮标在竖直方向的位移序列，确定海浪方向谱的海浪方向参数及海浪方向谱。
可选地，在本申请一些可能的实现方式中，根据浮标在横向上的位移序列、浮标在纵向上的位移序列，以及浮标在竖直方向的位移序列，确定海浪方向谱的海浪方向参数及海浪方向谱，包括：
根据浮标在横向上的位移序列以及浮标在纵向上的位移序列，得到目标海面的海浪的第一交叉谱，根据浮标在横向上的位移序列以及浮标在竖直方向的位移序列，得到海浪的第二交叉谱，根据浮标在纵向上的位移序列以及浮标在竖直方向的位移序列，得到海浪的第三交叉谱，根据两个浮标在横向上的位移序列，得到海浪的第四交叉谱，根据两个浮标在纵向上的位移序列，得到海浪的第五交叉谱，根据两个浮标在竖直方向的位移序列，得到海浪的第六交叉谱；
根据第一交叉谱、第二交叉谱、第三交叉谱、第四交叉谱、第五交叉谱以及第六交叉谱，确定海浪方向谱的海浪方向参数及海浪方向谱。
可选地，在本申请一些可能的实现方式中，海浪方向交叉谱包括第一交叉谱、第二交叉谱、第三交叉谱、第四交叉谱、第五交叉谱以及第六交叉谱，海浪方向交叉谱包括实部和虚部；根据第一交叉谱、第二交叉谱、第三交叉谱、第四交叉谱、第五交叉谱以及第六交叉谱，确定海浪方向谱的海浪方向参数及海浪方向谱，包括：
根据第一交叉谱实部的实部参数、第四交叉谱实部的实部参数、第五交叉谱实部的实部参数、第六交叉谱实部的实部参数、第二交叉谱虚部的虚部参数，以及第三交叉谱虚部的虚部参数，确定海浪方向谱的海浪方向参数；
根据海浪方向参数，确定海浪方向谱。
可选地，在本申请一些可能的实现方式中，海浪方向谱的海浪方向参数包括第一海浪方向参数、第二海浪方向参数、第三海浪方向参数、第四海浪方向参数以及第五海浪方向参数；根据第一交叉谱实部的实部参数、第四交叉谱实部的实部参数、第五交叉谱实部的实部参数、第六交叉谱实部的实部参数、第二交叉谱虚部的虚部参数，以及第三交叉谱虚部的虚部参数，确定海浪方向谱的海浪方向参数，包括：
根据第六交叉谱实部的实部参数，确定第一海浪方向参数；
根据第四交叉谱实部的实部参数、第五交叉谱实部的实部参数、第六交叉谱实部的实部参数、第二交叉谱虚部的虚部参数，确定第二海浪方向参数；
根据第四交叉谱实部的实部参数、第五交叉谱实部的实部参数、第六交叉谱实部的实部参数、第三交叉谱虚部的虚部参数，确定第三海浪方向参数；
根据所述第四交叉谱实部的实部参数和所述第五交叉谱实部的实部参数，确定所述第四海浪方向参数；
根据所述第四交叉谱实部的实部参数、所述第五交叉谱实部的实部参数，以及所述第一交叉谱实部的实部参数，确定所述第五海浪方向参数。
可选地，在本申请一些可能的实现方式中，根据海浪方向谱确定目标海面的风向，包括：
根据第二海浪方向参数以及第三海浪方向参数，确定目标海面的风向。
可选地，在本申请一些可能的实现方式中，根据第二海浪方向参数以及第三海浪方向参数，确定目标海面的风向，包括：
根据第二海浪方向参数以及第三海浪方向参数，确定目标海面多个频率对应的海浪波向；
根据多个海浪波向，确定海浪波向的平均值；
根据海浪波向的平均值，确定目标海面的风向。
可选地，在本申请一些可能的实现方式中，根据海浪能量谱确定目标海面的风速，包括：
基于海浪作用方程、目标海面的波浪与风场耦合关系，以及海浪能量谱确定目标海面的风速。
一方面，本申请实施例提供了海面风场反演装置，包括：
获取模块，用于获取目标海面中浮标的运动参数序列，运动参数序列包括浮标角度参数序列和浮标在东北天坐标系的东北天坐标参数序列，浮标角度参数序列包括横摇角参数序列、纵摇角参数序列、艏摇角参数序列，东北天坐标参数序列包括第一横向加速度序列、第一纵向加速度序列、第一竖直加速度序列；
第一确定模块，用于根据浮标的运动参数序列，确定浮标在笛卡尔坐标系的笛卡尔坐标参数序列，笛卡尔坐标参数序列包括第二横向加速度序列、第二纵向加速度序列、第二竖直加速度序列；
第二确定模块，用于根据第二竖直加速度序列，确定目标海面的海浪能量谱；
第三确定模块，用于根据海浪能量谱确定目标海面的风速；
第四确定模块，用于根据第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列，确定目标海面的海浪方向谱；
第五确定模块，用于根据海浪方向谱确定目标海面的风向。
一方面，本申请实施例提供了一种电子设备，电子设备包括存储器，处理器及存储在存储器中并可在处理器上运行的计算机程序，处理器执行计算机程序时实现如上述海面风场反演方法中的步骤。
一方面，本申请实施例还提供了一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序被处理器执行时实现上述海面风场反演方法中的步骤。
本申请实施例提供一种海面风场反演方法及装置，该方法首先获取目标海面中浮标的运动参数序列，根据浮标的运动参数序列，确定浮标在笛卡尔坐标系的第二横向加速度序列、第二纵向加速度序列、第二竖直加速度序列，再根据第二竖直加速度序列，确定目标海面的海浪能量谱，根据海浪能量谱确定目标海面的风速，再根据第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列，确定目标海面的海浪方向谱，最后根据海浪方向谱确定目标海面的风向。由于本申请只需要通过传感器等设备获取浮标的运动参数序列就能反演出目标海面的风速和风向，节省了海洋观测的成本。且本申请是先根据浮标的运动参数序列得到了目标海面的海浪方向谱和能量谱，再根据海浪方向谱和能量谱计算得到了海面的风速和方向，因此，实现了对海浪、风要素小范围、高时空分辨率的观测。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的海面风场反演系统的场景示意图。
图2是本申请实施例提供的海面风场反演方法的流程示意图。
图3a是本申请实施例提供的浮标的主视图。
图3b是本申请实施例提供的浮标的剖视图。
图4是本申请实施例提供的海面风场反演装置的结构示意图。
图5是本申请实施例提供的电子设备的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
本申请实施例提供海面风场反演方法及装置。其中，该海面风场反演装置可以集成在电子设备中，该电子设备可以是服务器，也可以是终端等设备，其中，该终端可以包括是智能手机、平板电脑、笔记本电脑、个人计算机（PC，Personal Computer）、微型处理盒子、或者其他设备等。
本申请说明书、权利要求书和附图中出现的术语“包括”和“具有”以及它们任何变形，意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可选地还包括没有列出的步骤或单元，而是可选地还包括没有列出的步骤或单元，或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。此外，术语“第一”、“第二”和“第三”等是用于区别不同的对象，而非用于描述特定的顺序。
在本申请中，对海面风场进行反演主要是针对风场中的风速和风向进行反演。
在本申请中，海浪谱（ocean wave spectrum）包括海浪的能量谱和方向谱，海浪谱是描述海浪内部能量相对于频率和方向分布的图谱。又称海洋能量谱，是研究海浪的重要概念。海浪谱不仅表明海浪内部构成，还能给出海浪的外部特征。海浪可视作由无限多个振幅不同、频率不同、方向不同、位相杂乱的组成波组成，这些组成波便构成海浪谱。此谱描述海浪能量相对于各组成波的分布，故又名能量谱、功率谱、方向谱。它是随机海浪的一个重要统计性质，它不仅包含着海浪的二阶信息，而且还直接给出海浪组成波能量相对于频率和方向的分布。它用于描述海浪内部能量相对于频率和方向的分布。
在本申请中，浮标角度参数包括横摇角参数、纵摇角参数，以及艏摇角参数，其中，横摇角参数指的是指浮标横向与水平面的夹角大小，纵摇角参数指的是浮标纵轴与水平面的夹角大小，艏摇角参数指浮标垂直方向与水平面的夹角大小。
请参阅图1，图1是本申请实施例所提供的海面风场反演系统的场景示意图，请参阅图1，仅以浮标10和电子设备11为例进行说明，电子设备11获取目标海面中浮标10的运动参数序列，运动参数序列包括浮标角度参数序列和浮标在东北天坐标系的东北天坐标参数序列，浮标角度参数序列包括横摇角参数序列、纵摇角参数序列、艏摇角参数序列，东北天坐标参数序列包括第一横向加速度序列、第一纵向加速度序列、第一竖直加速度序列，再根据浮标的运动参数序列，确定浮标在笛卡尔坐标系的笛卡尔坐标参数序列，笛卡尔坐标参数序列包括第二横向加速度序列、第二纵向加速度序列、第二竖直加速度序列，再根据浮标在笛卡尔坐标系中的第二竖直加速度序列，确定目标海面的海浪能量谱，根据海浪能量谱确定目标海面的风速，再根据浮标在笛卡尔坐标系中的第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列，确定目标海面的海浪方向谱，最后，根据海浪方向谱确定目标海面的风向。
需要说明的是，图1所示的海面风场反演的场景示意图仅仅是一个示例，本申请实施例描述的电子设备、浮标是为了更加清楚的说明本申请实施例的技术方案，并不构成对于本申请实施例提供的技术方案的限定，本领域普通技术人员可知，随着系统的演变和新业务局景的出现，本申请实施例提供的技术方案对于类似的技术问题，同样适用。以下分别进行详细说明。需说明的是，以下实施例的描述顺序不作为对实施例优选顺序的限定。
图2是本申请实施例提供的海面风场反演方法的流程示意图。如图2所示，海面风场反演方法应用于电子设备，该海面风场反演方法包括以下几个步骤：
步骤201、获取目标海面中浮标的运动参数序列。
在本申请实施例中，可以通过传感器获取目标海面中浮标的运动参数，其中传感器包括但不局限于MEMS运动传感器，MEMS传感器即微机电系统（Microelectro MechanicalSystems），MEMS传感器是采用微电子和微机械加工技术制造出来的新型传感器。与传统的传感器相比，它具有体积小、重量轻、成本低、功耗低、可靠性高、适于批量化生产、易于集成和实现智能化的特点。同时，在微米量级的特征尺寸使得它可以完成某些传统机械传感器所不能实现的功能。
其中，运动参数序列包括浮标角度参数序列和浮标在东北天坐标系的东北天坐标参数序列，浮标角度参数序列包括横摇角参数序列、纵摇角参数序列、艏摇角参数序列，东北天坐标参数序列包括第一横向加速度序列、第一纵向加速度序列、第一竖直加速度序列。其中，横摇角参数序列由浮标在不同的时间点的横摇角参数组成，纵摇角参数序列由浮标在不同的时间点的纵摇角参数组成，艏摇角参数序列由浮标在不同时间点的艏摇角参数组成。第一横向加速度序列由浮标在东北天坐标系的X轴方向上不同时间点的加速度组成，第一纵向加速度序列由浮标在东北天坐标系的Y轴方向上不同时间点的加速度组成，第一竖直加速度序列由浮标在东北天坐标系的Z轴方向上不同时间点的加速度组成。
步骤202、根据浮标的运动参数序列，确定浮标在笛卡尔坐标系的笛卡尔坐标参数序列，笛卡尔坐标参数序列包括第二横向加速度序列、第二纵向加速度序列、第二竖直加速度序列。
电子设备根据坐标转换公式和浮标角度参数将浮标在东北天坐标中的东北天坐标参数序列转换成东北天坐标系的东北天坐标参数序列。
步骤203、根据第二竖直加速度序列，确定目标海面的海浪能量谱。
可以通过神经网络模型自动得到海浪能量谱，将第二竖直加速度序列输入到训练好的能量谱确定模型中，得到目标海面的海浪能量谱。由于海浪能量谱是波高谱，是通过竖直方向的位移计算的能量谱密度，因此，也可以通过浮标在Z轴方向的第二竖直加速度序列计算求得目标海面的海量能量谱。
在一种实施例中，根据第二竖直加速度序列，确定目标海面的海浪能量谱，包括：对第二竖直加速度序列进行积分，得到浮标在竖直方向上的位移序列；将浮标在竖直方向上的位移序列分段，得到多个位移子序列；确定多个位移子序列的能量谱密度；根据多个位移子序列的能量谱密度，确定目标海面的海浪能量谱。
首先对浮标在z轴方向的第二竖直加速度序列az进行二次积分得到竖直方向的位移序列dz[n]，运用Wench方法将dz[n]分成L段，得到多个位移子序列，每个位移子序列有M个点，每段可以有部分重叠，分别求取每段的能量谱密度，然后将平均能量谱密度求和得到海浪能量谱，海浪能量谱计算公式如下：
（公式1）
其中，第二竖直加速度序列az有多个序列，n为第二竖直加速度序列的个数，f为海浪对应的海波频率，S(f)为海浪能量谱，j=为虚数单位，W[n]为窗函数（如海明窗函数、汉宁窗函数或平顶窗函数等），为归一化因子，l为第n个第二竖直加速度序列的段数，L为将第n个第二竖直加速度序列分成L段，M为位移子序列的点数。
步骤204、根据海浪能量谱确定目标海面的风速。
可以通过自适应的方式根据海浪的能量谱自动确定目标海面的风速，也可以根据公式计算得到目标海面的风速。
在一种实施例中，将海浪能量谱输入到训练好的风速确定模型中进行处理，得到目标海面的风速。
其中，风速确定模型是通过大量的样本数据进行训练得到的，可以将大量海浪能量谱作为深度神经网络的训练样本，同时获取标注的海浪能量谱对应的风速，将训练样本输入到风速确定模型中，得到预测的风速，根据预测风速和标注的风速，采用深度学习方法进行迭代训练，直至风速确定模型收敛时，得到训练好的风速确定模型。
在一种实施例中，根据海浪能量谱确定目标海面的风速，包括：基于海浪作用方程、目标海面的波浪与风场耦合关系，以及海浪能量谱确定目标海面的风速。
其中，海浪作用方程为：
（公式2）。
其中Sin为能量输入项，即风通过多种海洋运动形式向海洋内部输入的能量。SnL是非线性波-波相互作用，其物理意义是在频谱内波分量的共振机制交换和重新分配能量，浅水时以三波非线性相互作用为主，能量从低频传递到高频。深水时，以四波非线性相互作用为主，能量由峰值波频率传递至低频再到高频。Sdiss为能量耗散项，是波浪能量的整体耗散，主要包括白冠破碎、底摩擦、浅化引起的波浪破碎这三个物理过程。在频谱的高频端，三个过程之间会迅速形成平衡，即：。
利用海浪的波浪与驱动波场的风之间的固有耦合，假设风是平稳的，在平衡范围内存在：
（公式3）。
其中，且，其中，E 0为常数，为经验常数，I为方向扩展函数，g为重力加速度，U为风摩擦速度，S(f)为公式1中的海浪能量谱。
进一步地，通常风速计测量海平面以上4-5米的风速和风向，为获取不同高度的风信息，使用海表面风速的竖直剖面方程：
（公式4）。
其中，Uz为风速，K为冯·卡曼常数，Z 0为海表面粗糙度长度，，为经验常数，z为距离海表面高度，U为风摩擦速度。
联立公式1、公式3、公式4，求得目标海面的风速Uz。
步骤205：根据第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列，确定目标海面的海浪方向谱。
可以通过神经网络模型自动得到海浪能量谱，将第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列同时输入到训练好的方向谱确定模型中，得到目标海面的海浪方向谱。由于海浪是在三维空间上传播的，因此，可以通过第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列计算求得目标海面的海量方向谱。
在一种实施例中，根据第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列，确定目标海面的海浪方向谱，包括：对第二横向加速度序列进行积分，得到浮标在横向上的位移序列；对第二纵向加速度序列进行积分，得到浮标在纵向上的位移序列；对第二竖直加速度序列进行积分，得到浮标在竖直方向的位移序列；根据浮标在横向上的位移序列、浮标在纵向上的位移序列，以及浮标在竖直方向的位移序列，确定海浪方向谱的海浪方向参数及海浪方向谱。
对浮标在x、y轴方向的加速度序列ax、ay分别进行二次积分得到水平方向的位移序列d x[n]、纵向的位移序列d y[n]，再结合竖直方向的位移序列d z[n]。再通过计算三维位移交叉谱得到海浪方向谱，其中，交叉谱定义为交叉相关函数的傅里叶变换。
其中，根据浮标在横向上的位移序列、浮标在纵向上的位移序列，以及浮标在竖直方向的位移序列，确定海浪方向谱的海浪方向参数及海浪方向谱，包括：根据浮标在横向上的位移序列以及浮标在纵向上的位移序列，得到目标海面的海浪的第一交叉谱，根据浮标在横向上的位移序列以及浮标在竖直方向的位移序列，得到海浪的第二交叉谱，根据浮标在纵向上的位移序列以及浮标在竖直方向的位移序列，得到海浪的第三交叉谱，根据两个浮标在横向上的位移序列，得到海浪的第四交叉谱，根据两个浮标在纵向上的位移序列，得到海浪的第五交叉谱，根据两个浮标在竖直方向的位移序列，得到海浪的第六交叉谱；根据第一交叉谱、第二交叉谱、第三交叉谱、第四交叉谱、第五交叉谱以及第六交叉谱，确定海浪方向谱的海浪方向参数及海浪方向谱。
交叉函数的计算公式为：
（公式5）
其中，i,j=1,2,3…式中A i(t), A j(t)表示浮标基于MEMS运动传感器输出参数计算的位移序列d[n]，t表示时间，表示时间延时，T表示周期，表示浮标的运动函数，交叉谱计算公式为：（公式6）
首先，将浮标在横向上的位移序列dx[n]以及浮标在纵向上的位移序列dy[n]代入上述公式5中，得到关于τ的第一运动函数，再将第一运动函数代入公式6中求解积分得到海浪的第一交叉谱。将浮标在横向上的位移序列dx[n]以及浮标在竖直方向的位移序列dZ[n]代入上述公式5中，得到关于τ的第二运动函数，再将第二运动函数代入公式6中求解积分得到海浪的第二交叉谱。将浮标在纵向上的位移序列dy[n]以及浮标在竖直方向的位移序列dZ[n]代入上述公式5中，得到关于τ的第三运动函数，再将第三运动函数代入公式6中求解积分得到海浪的第三交叉谱。将两个浮标在横向上的位移序列dx[n]代入上述公式5中，得到关于τ的第四运动函数，再将第四运动函数代入公式6中求解积分得到海浪的第四交叉谱。将两个浮标在纵向上的位移序列dy[n]代入上述公式5中，得到关于τ的第五运动函数，再将第五运动函数代入公式6中求解积分得到海浪的第五交叉谱。将两个浮标在竖直方向的位移序列dZ[n]代入上述公式5中，得到关于τ的第六运动函数，再将第六运动函数代入公式6中求解积分得到海浪的第六交叉谱。
最后再根据第一交叉谱、第二交叉谱、第三交叉谱、第四交叉谱、所第五交叉谱以及第六交叉谱，确定海浪方向谱的海浪方向参数及海浪方向谱。
海浪方向谱的计算公式为：
（公式7）
其中，海浪方向谱的海浪方向参数是根据第一交叉谱、第二交叉谱、第三交叉谱、第四交叉谱、第五交叉谱以及第六交叉谱中交叉谱的实部和虚部的参数计算得到的，因此，需要先分别计算第一交叉谱、第二交叉谱、第三交叉谱、第四交叉谱、第五交叉谱以及第六交叉谱对应的实部和虚部参数。
在一种实施例中，海浪方向交叉谱包括第一交叉谱、第二交叉谱、第三交叉谱、第四交叉谱、第五交叉谱以及第六交叉谱，海浪方向交叉谱包括实部和虚部；根据第一交叉谱、第二交叉谱、第三交叉谱、第四交叉谱、第五交叉谱以及第六交叉谱，确定海浪方向谱的海浪方向参数及海浪方向谱，包括：根据第一交叉谱实部的实部参数、第四交叉谱实部的实部参数、第五交叉谱实部的实部参数、第六交叉谱实部的实部参数、第二交叉谱虚部的虚部参数，以及第三交叉谱虚部的虚部参数，确定海浪方向谱的海浪方向参数；根据海浪方向参数，确定海浪方向谱。
通常交叉谱可由复数表示，其中实部为同相谱，虚部为正交谱，o，p为三个方向(x，y，z)中的任意位移序列。假设第一交叉谱表示为,第二交叉谱表示为,第三交叉谱表示为,第四交叉谱表示为,第五交叉谱表示为,第六交叉谱表示为。
海浪方向参数为：
（公式8）。
将上述海浪方向交叉谱中的实部和虚部参数代入公式8中，得到海浪方向谱的海浪方向参数。
在一种实施例中，海浪方向谱的海浪方向参数包括第一海浪方向参数、第二海浪方向参数、第三海浪方向参数、第四海浪方向参数以及第五海浪方向参数；根据第一交叉谱实部的实部参数、第四交叉谱实部的实部参数、第五交叉谱实部的实部参数、第六交叉谱实部的实部参数、第二交叉谱虚部的虚部参数，以及第三交叉谱虚部的虚部参数，确定海浪方向谱的海浪方向参数，包括：根据第六交叉谱实部的实部参数，确定第一海浪方向参数；根据第四交叉谱实部的实部参数、第五交叉谱实部的实部参数、第六交叉谱实部的实部参数、第二交叉谱虚部的虚部参数，确定第二海浪方向参数；根据第四交叉谱实部的实部参数、第五交叉谱实部的实部参数、第六交叉谱实部的实部参数、第三交叉谱虚部的虚部参数，确定第三海浪方向参数；根据第四交叉谱实部的实部参数和第五交叉谱实部的实部参数，确定第四海浪方向参数；根据第四交叉谱实部的实部参数、第五交叉谱实部的实部参数，以及第一交叉谱实部的实部参数，确定第五海浪方向参数。
206：根据海浪方向谱确定目标海面的风向。
在实际确定方向的过程中，可以通过自适应的方式根据海浪的方向谱自动确定目标海面的方向，也可以根据公式计算得到目标海面的方向。
在一种实施例中，将海浪方向谱输入到训练好的方向确定模型中进行处理，得到目标海面的方向。
其中，方向确定模型是通过大量的样本数据进行训练得到的，可以将大量海浪方向谱作为深度神经网络的训练样本，同时获取标注的海浪方向谱对应的方向，将训练样本输入到方向确定模型中，得到预测的方向，根据预测方向和标注的方向，采用深度学习方法进行迭代训练，直至方向确定模型收敛时，得到训练好的方向确定模型。
由于从风到波的大部分能量转换发生在海浪谱的高频段，故可以合理将波传播方向近似为风向。使用谱分析方法、傅里叶级数法，来估算方向的分布函数：
（公式9），
其中，等于海浪的波浪角度，f为波浪频率，波浪在频率f下对应的角度，b1(f)在数值上等于公式8中的b1，即b1(f)等于第三海浪方向参数，a1(f)在数值上等于公式8中的a1，即a1(f)等于第二海浪方向参数。由于海浪的方向计算公式为：
（公式10）。
其中，为海浪的方向，为海浪的平均波向，N为高频段的个数。因此，首先需要先计算得到波浪的平均角度，可以取高频段中多个频率f1、f2、f3……fN分别对应的海浪波向，则。最后将代入公式10中，得到目标海面的风向。
在一种实施例中，根据海浪方向谱确定目标海面的风向，包括：根据第二海浪方向参数以及第三海浪方向参数，确定目标海面的风向。
其中，根据第二海浪方向参数以及第三海浪方向参数，确定目标海面的风向，包括：根据第二海浪方向参数以及第三海浪方向参数，确定目标海面多个频率对应的海浪波向；根据多个海浪波向，确定海浪波向的平均值；根据海浪波向的平均值，确定目标海面的风向。
本申请实施例提供一种海面风场反演方法，该方法首先获取目标海面中浮标的运动参数序列，根据浮标的运动参数序列，确定浮标在笛卡尔坐标系的第二横向加速度序列、第二纵向加速度序列、第二竖直加速度序列，再根据第二竖直加速度序列，确定目标海面的海浪能量谱，根据海浪能量谱确定目标海面的风速，再根据第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列，确定目标海面的海浪方向谱，最后根据海浪方向谱确定目标海面的风向。由于本申请只需要通过传感器等设备获取浮标的运动参数序列就能反演出目标海面的风速和方向，节省了海洋观测的成本。且本申请是先根据浮标的运动参数序列得到了目标海面的海浪方向谱和能量谱，再根据海浪方向谱和能量谱计算得到了海面的风速和方向，因此，实现了对海浪、风要素小范围、高时空分辨率的观测。
如图3a所示为本申请实施例提供的浮标的主视图，如图3b所示为本申请实施例提供的浮标的剖视图。
浮标主要包括天线罩2、柔性太阳能板3、仪器舱9和配重底托12。卫星通信模组1置于天线罩2内，仪器舱9与天线罩2通过О型圈4密封连接。系统采集和处理电路6、MEMS运动传感器采集与调理电路8、电池组11置于仪器舱9内，通过支撑板7与仪器舱9固定，泡沫浮体5通过仪器舱9底端配置双头螺杆10与配重底托12配合夹紧固定。
在上述实施例方法的基础上，本实施例将从海面风场反演装置的角度进一步进行描述，请参阅图4，图4具体描述了本申请实施例提供的海面风场反演装置的结构示意图，其可以包括：
获取模块401，用于获取目标海面中浮标的运动参数序列，运动参数序列包括浮标角度参数序列和浮标在东北天坐标系的东北天坐标参数序列，浮标角度参数序列包括横摇角参数序列、纵摇角参数序列、艏摇角参数序列，东北天坐标参数序列包括第一横向加速度序列、第一纵向加速度序列、第一竖直加速度序列；
第一确定模块402，用于根据浮标的运动参数序列，确定浮标在笛卡尔坐标系的笛卡尔坐标参数序列，笛卡尔坐标参数序列包括第二横向加速度序列、第二纵向加速度序列、第二竖直加速度序列；
第二确定模块403，用于根据第二竖直加速度序列，确定目标海面的海浪能量谱；
在一种实施例中，第二确定模块403，包括：
第一积分子模块，用于对第二竖直加速度序列进行积分，得到浮标在竖直方向上的位移序列；
分段子模块，用于将浮标在竖直方向上的位移序列分段，得到多个位移子序列；
第一确定子模块，用于确定多个位移子序列的能量谱密度；
第二确定子模块，用于根据多个位移子序列的能量谱密度，确定目标海面的海浪能量谱。
第三确定模块404，用于根据海浪能量谱确定目标海面的风速；
在一种实施例中，第三确定模块404，包括：
反演子模块，用于基于海浪作用方程、目标海面的波浪与风场耦合关系，以及海浪能量谱确定目标海面的风速。
第四确定模块405，用于根据第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列，确定目标海面的海浪方向谱；
在一种实施例中，第四确定模块405，包括：
第二积分子模块，用于对第二横向加速度序列进行积分，得到浮标在横向上的位移序列；
第三积分子模块，用于对第二纵向加速度序列进行积分，得到浮标在纵向上的位移序列；
第四积分子模块，用于对第二竖直加速度序列进行积分，得到浮标在竖直方向的位移序列；
第三确定子模块，用于根据浮标在横向上的位移序列、浮标在纵向上的位移序列，以及浮标在竖直方向的位移序列，确定海浪方向谱的海浪方向参数及海浪方向谱。
第五确定模块406，用于根据海浪方向谱确定目标海面的风向。
相应的，本申请实施例还提供一种电子设备，如图5所示，该电子设备可以包括射频电路501、包括有一个或一个以上计算机可读存储介质的存储器502、输入单元503、显示单元504、传感器505、音频电路506、WiFi模块507、包括有一个或者一个以上处理核心的处理器508、以及电源509等部件。本领域技术人员可以理解，图5中示出的电子设备结构并不构成对电子设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。其中：
射频电路501可用于收发信息或通话过程中，信号的接收和发送，特别地，将基站的下行信息接收后，交由一个或者一个以上处理器508处理；另外，将涉及上行的数据发送给基站。存储器502可用于存储软件程序以及模块，处理器508通过运行存储在存储器502的软件程序以及模块，从而执行各种功能应用以及数据处理。输入单元503可用于接收输入的数字或字符信息，以及产生与用户设置以及功能控制有关的键盘、鼠标、操作杆、光学或者轨迹球信号输入。
显示单元504可用于显示由用户输入的信息或提供给用户的信息以及电子设备的各种图形用户接口，这些图形用户接口可以由图形、文本、图标、视频和其任意组合来构成。
电子设备还可包括至少一种传感器505，比如光传感器、运动传感器以及其他传感器。音频电路506包括扬声器，扬声器可提供用户与电子设备之间的音频接口。
WiFi属于短距离无线传输技术，电子设备通过WiFi模块507可以帮助用户收发电子邮件、浏览网页和访问流式媒体等，它为用户提供了无线的宽带互联网访问，虽然图5示出了WiFi模块507，但是可以理解的是，其并不属于电子设备的必须构成，完全可以根据需要在不改变申请的本质的范围内而省略。
处理器508是电子设备的控制中心，利用各种接口和线路连接整个手机的各个部分，通过运行或执行存储在存储器502内的软件程序和/或模块，以及调用存储在存储器502内的数据，执行电子设备的各种功能和处理数据，从而对手机进行整体监控。
电子设备还包括给各个部件供电的电源509（比如电池），优选的，电源可以通过电源管理系统与处理器508逻辑相连，从而通过电源管理系统实现管理充电、放电、以及功耗管理等功能。
尽管未示出，电子设备还可以包括摄像头、蓝牙模块等，在此不再赘述。具体在本实施例中，电子设备中的处理器508会按照如下的指令，将一个或一个以上的应用程序的进程对应的可执行文件加载到存储器502中，并由处理器508来运行存储在存储器502中的应用程序，从而实现以下功能：
获取目标海面中浮标的运动参数序列，运动参数序列包括浮标角度参数序列和浮标在东北天坐标系的东北天坐标参数序列，浮标角度参数序列包括横摇角参数序列、纵摇角参数序列、艏摇角参数序列，东北天坐标参数序列包括第一横向加速度序列、第一纵向加速度序列、第一竖直加速度序列；
根据浮标的运动参数序列，确定浮标在笛卡尔坐标系的笛卡尔坐标参数序列，笛卡尔坐标参数序列包括第二横向加速度序列、第二纵向加速度序列、第二竖直加速度序列；
根据第二竖直加速度序列，确定目标海面的海浪能量谱；
根据海浪能量谱确定目标海面的风速；
根据第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列，确定目标海面的海浪方向谱；
根据海浪方向谱确定目标海面的风向。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见上文的详细描述，此处不再赘述。
本领域普通技术人员可以理解，上述实施例的各种方法中的全部或部分步骤可以通过指令来完成，或通过指令控制相关的硬件来完成，该指令可以存储于一计算机可读存储介质中，并由处理器进行加载和执行。
为此，本申请实施例提供一种存储介质，其中存储有多条指令，该指令能够被处理器进行加载，以实现以下功能：
获取目标海面中浮标的运动参数序列，运动参数序列包括浮标角度参数序列和浮标在东北天坐标系的东北天坐标参数序列，浮标角度参数序列包括横摇角参数序列、纵摇角参数序列、艏摇角参数序列，东北天坐标参数序列包括第一横向加速度序列、第一纵向加速度序列、第一竖直加速度序列；
根据浮标的运动参数序列，确定浮标在笛卡尔坐标系的笛卡尔坐标参数序列，笛卡尔坐标参数序列包括第二横向加速度序列、第二纵向加速度序列、第二竖直加速度序列；
根据第二竖直加速度序列，确定目标海面的海浪能量谱；
根据海浪能量谱确定目标海面的风速；
根据第二横向加速度序列、第二纵向加速度序列以及第二竖直加速度序列，确定目标海面的海浪方向谱；
根据海浪方向谱确定目标海面的风向。
以上各个操作的具体实施可参见前面的实施例，在此不再赘述。
其中，该存储介质可以包括：只读存储器（ROM，Read Only Memory）、随机存取记忆体（RAM，Random Access Memory）、磁盘或光盘等。
由于该存储介质中所存储的指令，可以执行本申请实施例所提供的任一种海面风场反演方法中的步骤，因此，可以实现本申请实施例所提供的任一种海面风场反演方法所能实现的有益效果，详见前面的实施例，在此不再赘述。
以上对本申请实施例所提供的一种海面风场反演方法及装置进行了详细介绍，本文中应用了具体个例对本申请的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本申请的技术方案及其核心思想；本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例的技术方案的范围。

Claims (8)

1.一种海面风场反演方法，其特征在于，海面风场反演系统包括浮标和电子设备，所述浮标与电子设备通信连接，所述风场反演方法应用于所述电子设备；
所述浮标包括MEMS运动传感器；
所述海面风场反演方法包括：
获取目标海面中浮标的运动参数序列，所述运动参数序列是由所述MEMS运动传感器采集的；
所述运动参数序列包括浮标角度参数序列和浮标在东北天坐标系的东北天坐标参数序列；
所述浮标角度参数序列包括横摇角参数序列、纵摇角参数序列、艏摇角参数序列；
所述东北天坐标参数序列包括第一横向加速度序列、第一纵向加速度序列、第一竖直加速度序列；
根据所述浮标的运动参数序列，确定所述浮标在笛卡尔坐标系的笛卡尔坐标参数序列，所述笛卡尔坐标参数序列包括第二横向加速度序列、第二纵向加速度序列、第二竖直加速度序列；
根据所述第二竖直加速度序列，确定所述目标海面的海浪能量谱，其中，所述根据所述第二竖直加速度序列，确定所述目标海面的海浪能量谱，包括：对所述第二竖直加速度序列进行积分，得到所述浮标在竖直方向上的位移序列；将所述浮标在竖直方向上的位移序列分段，得到多个位移子序列；确定所述多个位移子序列的能量谱密度；根据所述多个位移子序列的能量谱密度，确定所述目标海面的海浪能量谱；
根据所述海浪能量谱确定所述目标海面的风速；
根据所述第二横向加速度序列、所述第二纵向加速度序列以及所述第二竖直加速度序列，确定所述目标海面的海浪方向谱；
根据所述海浪方向谱确定所述目标海面的风向，其中，所述根据所述海浪能量谱确定所述目标海面的风速，包括：基于海浪作用方程、所述目标海面的波浪与风场耦合关系，以及所述海浪能量谱确定所述目标海面的风速。
2.如权利要求1所述的海面风场反演方法，其特征在于，所述根据所述第二横向加速度序列、所述第二纵向加速度序列以及所述第二竖直加速度序列，确定所述目标海面的海浪方向谱，包括：
对所述第二横向加速度序列进行积分，得到所述浮标在横向上的位移序列；
对所述第二纵向加速度序列进行积分，得到所述浮标在纵向上的位移序列；
对所述第二竖直加速度序列进行积分，得到所述浮标在竖直方向的位移序列；
根据所述浮标在横向上的位移序列、所述浮标在纵向上的位移序列，以及所述浮标在竖直方向的位移序列，确定所述海浪方向谱的海浪方向参数及所述海浪方向谱。
3.如权利要求2所述的海面风场反演方法，其特征在于，所述根据所述浮标在横向上的位移序列、所述浮标在纵向上的位移序列，以及所述浮标在竖直方向的位移序列，确定所述海浪方向谱的海浪方向参数及所述海浪方向谱，包括：
根据所述浮标在横向上的位移序列以及所述浮标在纵向上的位移序列，得到所述目标海面的海浪的第一交叉谱，根据所述浮标在横向上的位移序列以及所述浮标在竖直方向的位移序列，得到海浪的第二交叉谱，根据所述浮标在纵向上的位移序列以及所述浮标在竖直方向的位移序列，得到海浪的第三交叉谱，根据两个所述浮标在横向上的位移序列，得到海浪的第四交叉谱，根据两个所述浮标在纵向上的位移序列，得到所述海浪的第五交叉谱，根据两个所述浮标在竖直方向的位移序列，得到所述海浪的第六交叉谱；
根据所述第一交叉谱、所述第二交叉谱、所述第三交叉谱、所述第四交叉谱、所述第五交叉谱以及所述第六交叉谱，确定所述海浪方向谱的海浪方向参数及所述海浪方向谱。
4.如权利要求3所述的海面风场反演方法，其特征在于，海浪方向交叉谱包括所述第一交叉谱、所述第二交叉谱、所述第三交叉谱、所述第四交叉谱、所述第五交叉谱以及所述第六交叉谱，
所述海浪方向交叉谱包括实部和虚部；所述根据所述第一交叉谱、所述第二交叉谱、所述第三交叉谱、所述第四交叉谱、所述第五交叉谱以及所述第六交叉谱，确定所述海浪方向谱的海浪方向参数及所述海浪方向谱，包括：
根据所述第一交叉谱实部的实部参数、所述第四交叉谱实部的实部参数、所述第五交叉谱实部的实部参数、所述第六交叉谱实部的实部参数、所述第二交叉谱虚部的虚部参数，以及所述第三交叉谱虚部的虚部参数，确定所述海浪方向谱的海浪方向参数；
根据所述海浪方向参数，确定所述海浪方向谱。
5.如权利要求4所述的海面风场反演方法，其特征在于，海浪方向谱的海浪方向参数包括第一海浪方向参数、第二海浪方向参数、第三海浪方向参数、第四海浪方向参数以及第五海浪方向参数；所述根据所述第一交叉谱实部的实部参数、所述第四交叉谱实部的实部参数、所述第五交叉谱实部的实部参数、所述第六交叉谱实部的实部参数、所述第二交叉谱虚部的虚部参数，以及所述第三交叉谱虚部的虚部参数，确定所述海浪方向谱的海浪方向参数，包括：
根据所述第六交叉谱实部的实部参数，确定所述第一海浪方向参数；
根据所述第四交叉谱实部的实部参数、所述第五交叉谱实部的实部参数、所述第六交叉谱实部的实部参数、所述第二交叉谱虚部的虚部参数，确定所述第二海浪方向参数；
根据所述第四交叉谱实部的实部参数、所述第五交叉谱实部的实部参数、所述第六交叉谱实部的实部参数、所述第三交叉谱虚部的虚部参数，确定所述第三海浪方向参数；
根据所述第四交叉谱实部的实部参数和所述第五交叉谱实部的实部参数，确定所述第四海浪方向参数；
根据所述第四交叉谱实部的实部参数、所述第五交叉谱实部的实部参数，以及所述第一交叉谱实部的实部参数，确定所述第五海浪方向参数。
6.如权利要求5所述的海面风场反演方法，其特征在于，所述根据所述海浪方向谱确定所述目标海面的风向，包括：
根据所述第二海浪方向参数以及所述第三海浪方向参数，确定所述目标海面的风向。
7.如权利要求6所述的海面风场反演方法，其特征在于，所述根据所述第二海浪方向参数以及所述第三海浪方向参数，确定所述目标海面的风向，包括：
根据所述第二海浪方向参数以及所述第三海浪方向参数，确定所述目标海面多个频率对应的海浪波向；
根据多个所述海浪波向，确定所述海浪波向的平均值；
根据所述海浪波向的平均值，确定所述目标海面的风向。
8.一种海面风场反演装置，其特征在于，海面风场反演系统包括浮标和电子设备，所述浮标与电子设备通信连接，所述风场反演装置应用于所述电子设备；
所述浮标包括MEMS运动传感器；
所述装置包括：
获取模块，用于获取目标海面中浮标的运动参数序列，所述运动参数序列是由所述MEMS运动传感器采集的；
所述运动参数序列包括浮标角度参数序列和浮标在东北天坐标系的东北天坐标参数序列，所述浮标角度参数序列包括横摇角参数序列、纵摇角参数序列、艏摇角参数序列，所述东北天坐标参数序列包括第一横向加速度序列、第一纵向加速度序列、第一竖直加速度序列；
第一确定模块，用于根据所述浮标的运动参数序列，确定所述浮标在笛卡尔坐标系的笛卡尔坐标参数序列，所述笛卡尔坐标参数序列包括第二横向加速度序列、第二纵向加速度序列、第二竖直加速度序列；
第二确定模块，用于根据所述第二竖直加速度序列，确定所述目标海面的海浪能量谱，其中，所述根据所述第二竖直加速度序列，确定所述目标海面的海浪能量谱，包括：对所述第二竖直加速度序列进行积分，得到所述浮标在竖直方向上的位移序列；将所述浮标在竖直方向上的位移序列分段，得到多个位移子序列；确定所述多个位移子序列的能量谱密度；根据所述多个位移子序列的能量谱密度，确定所述目标海面的海浪能量谱；
第三确定模块，用于根据所述海浪能量谱确定所述目标海面的风速；
第四确定模块，用于根据所述第二横向加速度序列、所述第二纵向加速度序列以及所述第二竖直加速度序列，确定所述目标海面的海浪方向谱；
第五确定模块，用于根据所述海浪方向谱确定所述目标海面的风向，其中，所述根据所述海浪能量谱确定所述目标海面的风速，包括：基于海浪作用方程、所述目标海面的波浪与风场耦合关系，以及所述海浪能量谱确定所述目标海面的风速。
CN202211002563.5A 2022-08-22 2022-08-22 海面风场反演方法及装置 Active CN115327669B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211002563.5A CN115327669B (zh) 2022-08-22 2022-08-22 海面风场反演方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211002563.5A CN115327669B (zh) 2022-08-22 2022-08-22 海面风场反演方法及装置

Publications (2)

Publication Number Publication Date
CN115327669A CN115327669A (zh) 2022-11-11
CN115327669B true CN115327669B (zh) 2024-10-15

Family

ID=83926755

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211002563.5A Active CN115327669B (zh) 2022-08-22 2022-08-22 海面风场反演方法及装置

Country Status (1)

Country Link
CN (1) CN115327669B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117664187A (zh) * 2024-02-02 2024-03-08 国家海洋技术中心 一种基于水声通讯的潜标实时健康状态监测系统及方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101813476A (zh) * 2010-03-19 2010-08-25 天津大学 近海波浪参数立体实时监测系统
CN114620186A (zh) * 2022-03-15 2022-06-14 国家海洋技术中心 一种小型化海洋测波浮标装置与测波及风要素反演方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006284242A (ja) * 2005-03-31 2006-10-19 Furuno Electric Co Ltd 潮流測定装置及び潮流測定方法
US10137972B2 (en) * 2014-07-28 2018-11-27 Furuno Electric Co., Ltd. Vessel characteristic estimation device and automatic steering device
CN209938881U (zh) * 2019-04-11 2020-01-14 南京信息工程大学 基于六维加速度传感器的海洋波浪浮标

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101813476A (zh) * 2010-03-19 2010-08-25 天津大学 近海波浪参数立体实时监测系统
CN114620186A (zh) * 2022-03-15 2022-06-14 国家海洋技术中心 一种小型化海洋测波浮标装置与测波及风要素反演方法

Also Published As

Publication number Publication date
CN115327669A (zh) 2022-11-11

Similar Documents

Publication Publication Date Title
CN112394349B (zh) 海洋浮标智能安防监控方法、装置、系统及存储介质
CN114580223B (zh) 风致漂移预测方法、装置、设备及存储介质
CN115327669B (zh) 海面风场反演方法及装置
CN114580309B (zh) 漂移速度预测方法、装置、计算设备及存储介质
CN117113798A (zh) 漂浮物的漂移速度预测方法、装置、设备及存储介质
CN112762935A (zh) 一种基于船体姿态监测的定位方法及系统
CN114620186A (zh) 一种小型化海洋测波浮标装置与测波及风要素反演方法
CN103323817A (zh) 一种机载合成孔径雷达海面风矢量反演方法
CN112415467B (zh) 一种基于神经网络的单矢量潜标目标定位实现方法
CN106628006A (zh) 一种智能手机波浪浮标
US20230167796A1 (en) Method for predicting a characteristic resulting from a swell on the basis of a spectral model of the swell
CN119047177B (zh) 一种用于水声探测仿真的数据处理方法及装置
CN118050541B (zh) 一种基于悬空动力机构的加速度计校准方法、系统及装置
Ding A ship-motion prediction algorithm based on modified covariance method and neural networks
Wang et al. Research on a perceptual system of small unmanned surface vehicle
Zheng et al. Design and test verification of integrated navigation system based on SINS/DVL for underwater search and rescue robot
Pimenov et al. Investigation of the Features of Radio Wave Propagation over the Sea Surface
CN119164425B (zh) 场景适应性多导航源动态优选与评估系统及方法
CN117592381B (zh) 一种大气波导参数反演模型训练方法、装置、设备及介质
Sun et al. A Sensor-Based Shape Estimation Method for Towed Array
CN117805877B (zh) 一种船用定位方法、装置、电子设备及存储介质
CN117554920B (zh) 一种水面探测方法、装置、存储介质及电子设备
CN116539918A (zh) 一种基于悬空动力机构的加速度计校准方法、系统及装置
CN109141521A (zh) 一种适用于船舶姿态和升沉及振动测量的一体化检测装置和方法
Chen et al. Kalman Filter Based Integrated Underwater Acoustic Localization Using SINS and Long Baseline Technique

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