CN117825745A - Typhoon wind speed real-time observation method and equipment based on single vector hydrophone - Google Patents
Typhoon wind speed real-time observation method and equipment based on single vector hydrophone Download PDFInfo
- Publication number
- CN117825745A CN117825745A CN202410158331.1A CN202410158331A CN117825745A CN 117825745 A CN117825745 A CN 117825745A CN 202410158331 A CN202410158331 A CN 202410158331A CN 117825745 A CN117825745 A CN 117825745A
- Authority
- CN
- China
- Prior art keywords
- typhoon
- wind speed
- vector hydrophone
- real
- sea surface
- 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000001228 spectrum Methods 0.000 claims abstract description 55
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims description 44
- 230000006870 function Effects 0.000 claims description 36
- 238000004891 communication Methods 0.000 claims description 25
- 230000003595 spectral effect Effects 0.000 claims description 15
- 238000005457 optimization Methods 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 230000002068 genetic effect Effects 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 6
- 230000005284 excitation Effects 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 16
- 238000013500 data storage Methods 0.000 description 9
- 238000004088 simulation Methods 0.000 description 9
- 229910052741 iridium Inorganic materials 0.000 description 8
- GKOZUEZYRPOHIO-UHFFFAOYSA-N iridium atom Chemical compound [Ir] GKOZUEZYRPOHIO-UHFFFAOYSA-N 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 6
- 239000002131 composite material Substances 0.000 description 5
- 230000005693 optoelectronics Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 241001386813 Kraken Species 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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
- G01P5/24—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/10—Devices for predicting weather conditions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Aviation & Aerospace Engineering (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本申请公开了一种基于单矢量水听器的台风风速实时观测方法及设备,包括:将观测设备布放在台风风场外围海上,使用单矢量水听器采集台风环境下矢量水听器布放深度处的噪声声场数据,根据噪声声场数据实时计算噪声复声强垂向分量实部,根据经验公式,建立风成噪声源谱强度模型,根据风成噪声源谱强度模型,建立简正波水声传播模型,根据简正波水声传播模型和噪声复声强垂向分量实部,建立矢量水听器上方海表面风速的代价函数,估计矢量水听器上方海表面的风速,实时接收台风眼气压、台风眼壁半径信息,结合矢量水听器上方海表面估计的实时风速,动态计算海表面风速,降低了台风观测和预报的实施难度和成本,提高了时效性和精度。
The present application discloses a typhoon wind speed real-time observation method and equipment based on a single vector hydrophone, comprising: deploying the observation equipment on the sea outside the typhoon wind field, using the single vector hydrophone to collect noise sound field data at the deployment depth of the vector hydrophone in a typhoon environment, calculating the real part of the vertical component of the complex sound intensity of the noise in real time according to the noise sound field data, establishing a wind-induced noise source spectrum intensity model according to an empirical formula, establishing a simple normal wave water acoustic propagation model according to the wind-induced noise source spectrum intensity model, establishing a cost function of the sea surface wind speed above the vector hydrophone according to the simple normal wave water acoustic propagation model and the real part of the vertical component of the complex sound intensity of the noise, estimating the wind speed of the sea surface above the vector hydrophone, receiving the typhoon eye pressure and typhoon eye wall radius information in real time, combining the real-time wind speed estimated on the sea surface above the vector hydrophone, and dynamically calculating the sea surface wind speed, thereby reducing the implementation difficulty and cost of typhoon observation and forecasting, and improving timeliness and accuracy.
Description
技术领域Technical Field
本申请涉及海洋观测和探测技术领域,特别是涉及一种基于单矢量水听器的台风风速实时观测方法、系统、台风观测设备及计算机可读存储介质。The present application relates to the field of ocean observation and detection technology, and in particular to a method and system for real-time observation of typhoon wind speed based on a single vector hydrophone, a typhoon observation device, and a computer-readable storage medium.
背景技术Background technique
台风是严重威胁人类生活活动的灾害性天气之一,观测和预报台风风速等参数可以为人类提供及时的信息,进而有效制定应急措施,对降低台风的灾害影响具有重要的意义。Typhoon is one of the disastrous weather types that seriously threatens human life and activities. Observing and forecasting typhoon wind speed and other parameters can provide humans with timely information, thereby effectively formulating emergency measures, which is of great significance to reducing the impact of typhoon disasters.
在现有技术中,台风观测和预报主要是通过卫星遥感技术和气象侦察飞机,利用卫星遥感技术虽然可以大范围有效监测台风形成及移动路径,但对于台风风速等参数的预报精度仍待提高;而利用侦察飞机实地观测台风实施难度大,且成本高昂,飞机和人员安全都存在很高的风险,因此,侦察飞机也很少使用。In the existing technology, typhoon observation and forecasting are mainly carried out through satellite remote sensing technology and meteorological reconnaissance aircraft. Although satellite remote sensing technology can effectively monitor the formation and movement path of typhoons over a large area, the forecast accuracy of parameters such as typhoon wind speed still needs to be improved. On-site observation of typhoons using reconnaissance aircraft is difficult and costly, and there are high risks to the safety of aircraft and personnel. Therefore, reconnaissance aircraft are rarely used.
发明内容Summary of the invention
有鉴于此,本申请提供了一种单基于矢量水听器的台风风速实时观测方法、系统、台风观测设备及计算机可读存储介质,以解决现有技术中通过卫星遥感技术进行台风观测和预报,带来的观测和预报精度差,通过气象侦察飞机进行台风观测和预报带来的实施难度大,成本高,安全性低的问题。In view of this, the present application provides a typhoon wind speed real-time observation method, system, typhoon observation equipment and computer-readable storage medium based solely on vector hydrophones to solve the problems in the prior art of poor observation and forecasting accuracy caused by typhoon observation and forecasting through satellite remote sensing technology, and the difficulty, high cost and low safety of typhoon observation and forecasting through meteorological reconnaissance aircraft.
本申请提出一种基于单矢量水听器的台风风速实时观测方法,所述基于单矢量水听器的台风风速实时观测方法包括:The present application proposes a typhoon wind speed real-time observation method based on a single vector hydrophone, and the typhoon wind speed real-time observation method based on a single vector hydrophone includes:
将台风观测设备布放在台风风场外围海上,使用台风观测设备的单个矢量水听器采集台风环境下所述矢量水听器布放深度处的噪声声场数据;Deploy typhoon observation equipment on the sea outside the typhoon wind field, and use a single vector hydrophone of the typhoon observation equipment to collect noise sound field data at the deployment depth of the vector hydrophone in a typhoon environment;
根据所述噪声声场数据实时计算噪声复声强垂向分量实部;Calculate the real part of the vertical component of the complex sound intensity of the noise in real time according to the noise sound field data;
根据经验公式,建立风成噪声源谱强度模型;According to the empirical formula, the spectrum intensity model of wind-generated noise source is established;
根据所述风成噪声源谱强度模型,建立简正波水声传播模型;According to the wind-generated noise source spectrum intensity model, a simple normal wave underwater sound propagation model is established;
根据所述简正波水声传播模型和所述噪声复声强垂向分量实部,建立矢量水听器上方海表面风速的代价函数,估计所述矢量水听器上方海表面的风速;According to the simple normal wave underwater acoustic propagation model and the real part of the vertical component of the complex sound intensity of the noise, a cost function of the wind speed on the sea surface above the vector hydrophone is established to estimate the wind speed on the sea surface above the vector hydrophone;
通过卫星实时接收岸基中心发送的台风眼气压和台风眼壁半径信息,结合所述矢量水听器上方海表面估计的实时风速,动态计算距离台风中心不同距离上的海表面风速。The typhoon eye pressure and typhoon eye wall radius information sent by the shore-based center are received in real time via satellite, and the sea surface wind speed at different distances from the typhoon center is dynamically calculated in combination with the real-time wind speed estimated on the sea surface above the vector hydrophone.
可选地,所述噪声声场数据包括声压分量和垂直振速分量;所述根据所述噪声声场数据实时计算噪声复声强垂向分量实部,具体包括:Optionally, the noise sound field data includes a sound pressure component and a vertical vibration velocity component; and the real part of the vertical component of the complex sound intensity of the noise is calculated in real time according to the noise sound field data, specifically comprising:
获取第一时间长度的噪声声场数据,对所述声压分量和所述垂直振速分量进行傅里叶变换,得到第一时间长度的声压复频谱和垂直振速复频谱;Acquire noise sound field data of a first time length, perform Fourier transform on the sound pressure component and the vertical vibration velocity component, and obtain a sound pressure complex spectrum and a vertical vibration velocity complex spectrum of the first time length;
根据所述声压复频谱和所述垂直振速复频谱,计算所述矢量水听器的实测噪声复声强垂向分量实部;Calculating the real part of the vertical component of the complex sound intensity of the measured noise of the vector hydrophone according to the complex sound pressure spectrum and the vertical vibration velocity complex spectrum;
其中,所述实测噪声复声强垂向分量实部的计算公式为:The calculation formula of the real part of the vertical component of the measured noise complex sound intensity is:
其中,Iz实测(f,zr)为实测噪声复声强垂向分量实部;f为信号的中心频率;Re为复数量的实部;*为复共轭符号;<·>表示滑动窗平均周期图;zr为矢量水听器布放水深,r为矢量水听器距离台风中心水平距离,Sp(f,zr)为声压复频谱,Svz(f,zr)为垂直振速复频谱。Among them, Izmeasured (f, zr ) is the real part of the vertical component of the measured noise complex sound intensity; f is the center frequency of the signal; Re is the real part of the complex quantity; * is the complex conjugate symbol; <·> represents the sliding window average periodogram; zr is the water depth of the vector hydrophone, r is the horizontal distance between the vector hydrophone and the typhoon center, Sp (f, zr ) is the complex spectrum of sound pressure, and Svz (f, zr ) is the complex spectrum of vertical vibration velocity.
可选地,所述根据经验公式,建立风成噪声源谱强度模型,具体包括:Optionally, establishing a wind-induced noise source spectrum intensity model according to an empirical formula specifically includes:
根据海洋环境风成噪声源的频率和风速,以及风成噪声源经验公式,建立台风环境下海表面的风成噪声源谱强度模型;According to the frequency and wind speed of wind-induced noise sources in the marine environment and the empirical formula of wind-induced noise sources, a spectrum intensity model of wind-induced noise sources on the sea surface in a typhoon environment is established.
所述风成噪声源谱强度模型计算公式为:The calculation formula of the wind-generated noise source spectrum intensity model is:
SIw(U,f,zs)=C×fq×Un;SI w (U,f,z s ) = C × f q × U n ;
zs=cw/4f;z s = c w /4f;
其中,SIw(U,f,zs)为风成噪声源谱强度;U为海表面风速,U的参数寻优范围为4-16;zs为噪声源设置水深;cw为声源处声波在水体中的传播速度;C为第一参数;q为第二参数;n为第三参数,zs为台风激发下的各个风成噪声源均匀分布在大海平面的深度。Among them, SI w (U,f,z s ) is the spectral intensity of the wind-induced noise source; U is the sea surface wind speed, and the parameter optimization range of U is 4-16; z s is the water depth of the noise source; c w is the propagation speed of the sound wave in the water body at the sound source; C is the first parameter; q is the second parameter; n is the third parameter, and z s is the depth at which each wind-induced noise source excited by the typhoon is evenly distributed on the sea surface.
可选地,所述根据所述风成噪声源谱强度模型,建立简正波水声传播模型,具体包括:Optionally, establishing a normal wave underwater sound propagation model according to the wind-induced noise source spectrum intensity model specifically includes:
根据台风环境下海表面的所述风成噪声源谱强度模型和简正波理论,建立矢量水听器接收位置处台风激发下风成噪声源的水下简正波水声传播模型;According to the spectral intensity model of wind-induced noise source on the sea surface under typhoon environment and the simple normal wave theory, an underwater simple normal wave acoustic propagation model of wind-induced noise source under typhoon excitation at the receiving position of the vector hydrophone is established;
所述水下简正波水声传播模型的计算公式为:The calculation formula of the underwater normal wave acoustic propagation model is:
km=αm+iβm; km = αm + iβm ;
其中,Iz(f,zr)为水下简正波水声传播数据,ρ=1.0,ρ为水体密度,ψm为噪声源信号中心频率为f时第m个垂直模态特征函数,km为相应的模态水平传播特征值,αm为km的实部,βm为km的虚部,ψm′为特征函数导数,i为虚数单位,Re为复数量的实部,SIw(U,f,zs,zr)为水听器接收点上方海表面风速对应的风成噪声源谱强度,Ur为水听器上方海表面风速。Among them, Iz (f, zr ) is the underwater simple normal wave acoustic propagation data, ρ=1.0, ρ is the water density, ψm is the mth vertical mode characteristic function when the center frequency of the noise source signal is f, km is the corresponding modal horizontal propagation eigenvalue, αm is the real part of km , βm is the imaginary part of km , ψm ′ is the derivative of the characteristic function, i is the imaginary unit, Re is the real part of the complex number, SIw (U,f, zs , zr ) is the spectral intensity of the wind-induced noise source corresponding to the sea surface wind speed above the hydrophone receiving point, and Ur is the sea surface wind speed above the hydrophone.
可选地,所述根据所述简正波水声传播模型和所述噪声复声强垂向分量实部,建立矢量水听器上方海表面风速的代价函数,估计所述矢量水听器上方海表面的风速,具体包括:Optionally, establishing a cost function of the sea surface wind speed above the vector hydrophone according to the simple normal wave underwater acoustic propagation model and the real part of the vertical component of the complex sound intensity of the noise, and estimating the sea surface wind speed above the vector hydrophone specifically includes:
根据所述简正波水声传播模型和所述噪声复声强垂向分量实部,利用多频联合反演的方式,采用多参数并行粒子群遗传优化算法,建立所述矢量水听器上方海表面风速的关于所述第一参数、所述第二参数和所述第三参数的代价函数,得到最优的所述第一参数、所述第二参数和所述第三参数;According to the simple normal wave underwater acoustic propagation model and the real part of the vertical component of the complex sound intensity of the noise, a cost function of the sea surface wind speed above the vector hydrophone with respect to the first parameter, the second parameter and the third parameter is established by using a multi-frequency joint inversion method and a multi-parameter parallel particle swarm genetic optimization algorithm to obtain the optimal first parameter, the second parameter and the third parameter;
所述代价函数的计算公式为:The calculation formula of the cost function is:
其中,E(C,q,n)为代价函数,J表示用于反演计算的频率的个数,n的参数寻优范围为2.5-3.5。Among them, E(C,q,n) is the cost function, J represents the number of frequencies used for inversion calculation, and the parameter optimization range of n is 2.5-3.5.
可选地,所述通过卫星实时接收岸基中心发送的台风眼气压和台风眼壁半径信息,结合所述矢量水听器上方海表面估计的实时风速,动态计算距离台风中心不同距离上的海表面风速,具体包括:Optionally, the receiving of typhoon eye pressure and typhoon eye wall radius information sent by a shore-based center in real time via a satellite, combined with the estimated real-time wind speed of the sea surface above the vector hydrophone, and dynamically calculating the sea surface wind speed at different distances from the typhoon center specifically includes:
获取最优的所述第一参数、所述第二参数和所述第三参数,得到所述矢量水听器上方海表面的风速;Obtaining the optimal first parameter, second parameter and third parameter to obtain the wind speed of the sea surface above the vector hydrophone;
通过卫星实时接收所述岸基中心发送的台风眼气压和台风眼壁半径信息;receiving in real time via satellite the typhoon eye pressure and typhoon eye wall radius information sent by the shore-based center;
结合所述矢量水听器上方海表面估计的实时风速、所述台风眼气压和所述台风眼壁半径,使用台风风场计算模型,得到距离台风中心不同距离的海表面风速的动态计算表达式,得到所述海表面风速。The real-time wind speed estimated on the sea surface above the vector hydrophone, the typhoon eye pressure and the typhoon eyewall radius are combined, and a typhoon wind field calculation model is used to obtain a dynamic calculation expression for the sea surface wind speed at different distances from the typhoon center to obtain the sea surface wind speed.
可选地,所述海表面风速的动态计算表达式为:Optionally, the dynamic calculation expression of the sea surface wind speed is:
Rw=A1/B; Rw = A1 /B ;
其中,U(d)为矢量水听器处海表面风速,pn为台风眼气压,pc为台风风场外围气压,Rw为台风眼壁半径,e≈2.71828,为自然常数,ρa=1.15,为空气密度,A为第一中间参数,B为第二中间参数,d为被观测风速处与台风中心的水平距离。Among them, U(d) is the sea surface wind speed at the vector hydrophone, pn is the typhoon eye pressure, pc is the periphery of the typhoon wind field, Rw is the radius of the typhoon eyewall, e≈2.71828 is a natural constant, ρa =1.15 is the air density, A is the first intermediate parameter, B is the second intermediate parameter, and d is the horizontal distance between the observed wind speed and the typhoon center.
本申请还提出一种基于单矢量水听器的台风风速实时观测系统,所述基于单矢量水听器的台风风速实时观测系统包括:The present application also proposes a typhoon wind speed real-time observation system based on a single vector hydrophone, the typhoon wind speed real-time observation system based on a single vector hydrophone comprises:
噪声声场数据采集模块,用于将台风观测设备布放在台风风场外围海上,使用台风观测设备的单个矢量水听器采集台风环境下所述矢量水听器布放深度处的噪声声场数据;A noise sound field data acquisition module is used to deploy typhoon observation equipment on the sea outside the typhoon wind field, and use a single vector hydrophone of the typhoon observation equipment to collect noise sound field data at the deployment depth of the vector hydrophone in a typhoon environment;
噪声复声强垂向分量实部计算模块,用于根据所述噪声声场数据实时计算噪声复声强垂向分量实部;A noise complex sound intensity vertical component real part calculation module, used for calculating the noise complex sound intensity vertical component real part in real time according to the noise sound field data;
风成噪声源谱强度模型建立模块,用于根据经验公式,建立风成噪声源谱强度模型;A module for establishing a spectrum intensity model of wind-generated noise sources is used to establish a spectrum intensity model of wind-generated noise sources based on empirical formulas;
简正波水声传播模型建立模块,用于根据所述风成噪声源谱强度模型,建立简正波水声传播模型;A simple normal wave underwater acoustic propagation model establishment module is used to establish a simple normal wave underwater acoustic propagation model according to the wind-induced noise source spectrum intensity model;
矢量水听器上方海表面风速观测模块,用于根据所述简正波水声传播模型和所述噪声复声强垂向分量实部,建立矢量水听器上方海表面风速的代价函数,估计所述矢量水听器上方海表面的风速;The sea surface wind speed observation module above the vector hydrophone is used to establish a cost function of the sea surface wind speed above the vector hydrophone according to the simple normal wave underwater acoustic propagation model and the real part of the vertical component of the complex sound intensity of the noise, and estimate the sea surface wind speed above the vector hydrophone;
海表面风速实时观测模块,用于通过卫星实时接收岸基中心发送的台风眼气压和台风眼壁半径信息,结合所述矢量水听器上方海表面估计的实时风速,动态计算距离台风中心不同距离上的海表面风速。The sea surface wind speed real-time observation module is used to receive the typhoon eye pressure and typhoon eye wall radius information sent by the shore-based center in real time via satellite, and dynamically calculate the sea surface wind speed at different distances from the typhoon center in combination with the real-time wind speed estimated on the sea surface above the vector hydrophone.
本申请还提出一种台风观测设备,所述台风观测设备包括:单个矢量水听器、卫星通信模块、存储器、处理器及存储在所述存储器上并可在所述处理器上运行的基于单矢量水听器的台风风速实时观测程序,所述基于单矢量水听器的台风风速实时观测程序被所述处理器执行时实现如上所述的基于单矢量水听器的台风风速实时观测方法的步骤。The present application also proposes a typhoon observation device, which includes: a single vector hydrophone, a satellite communication module, a memory, a processor, and a typhoon wind speed real-time observation program based on a single vector hydrophone stored in the memory and executable on the processor. When the typhoon wind speed real-time observation program based on a single vector hydrophone is executed by the processor, the steps of the typhoon wind speed real-time observation method based on a single vector hydrophone as described above are implemented.
本申请还提出一种计算机可读存储介质,所述计算机可读存储介质存储有基于单矢量水听器的台风风速实时观测程序,所述基于单矢量水听器的台风风速实时观测程序被处理器执行时实现如上所述的基于单矢量水听器的台风风速实时观测方法的步骤。The present application also proposes a computer-readable storage medium, which stores a typhoon wind speed real-time observation program based on a single vector hydrophone. When the typhoon wind speed real-time observation program based on a single vector hydrophone is executed by a processor, the steps of the typhoon wind speed real-time observation method based on a single vector hydrophone as described above are implemented.
本申请的有益效果是:区别于现有技术,本申请通过将所述台风观测设备布放在台风风场外围海上,使用台风观测设备的单矢量水听器采集台风环境下矢量水听器布放深度处的噪声声场数据,并利用岸基遥感卫星数据实时接收噪声数据,动态优化台风风速模型参数,实现不同距离台风风场风速的动态计算和实时回传,有效提高台风风速观测和预报的精度及效率;其次,本申请根据噪声声场数据实时计算噪声复声强垂向分量实部,根据经验公式,建立风成噪声源谱强度模型,根据风成噪声源谱强度模型,建立简正波水声传播模型,根据简正波水声传播模型和噪声复声强垂向分量实部,建立矢量水听器上方海表面风速的代价函数,估计矢量水听器上方海表面的风速,集中获取水听器海表面噪声特性,有效滤除了其他噪声干扰,降低了实施难度和成本,提高了观测和预报精度;最后,本申请通过通过卫星实时接收岸基中心发送的台风眼气压、台风眼壁半径信息,结合矢量水听器上方海表面估计的实时风速,动态计算距离台风中心不同距离上的海表面风速,并将计算结果通过铱星通信实时发送至岸基,实现台风风场风速信息的实时观测和预报,精准高效且易于工程应用。The beneficial effects of the present application are as follows: Different from the prior art, the present application deploys the typhoon observation equipment at sea outside the typhoon wind field, uses the single vector hydrophone of the typhoon observation equipment to collect the noise sound field data at the depth of the vector hydrophone deployment in the typhoon environment, and uses shore-based remote sensing satellite data to receive noise data in real time, dynamically optimizes the typhoon wind speed model parameters, realizes the dynamic calculation and real-time feedback of wind speeds in typhoon wind fields at different distances, and effectively improves the accuracy and efficiency of typhoon wind speed observation and forecasting; secondly, the present application calculates the real part of the vertical component of the complex sound intensity of the noise in real time according to the noise sound field data, establishes a spectral intensity model of the wind-induced noise source according to the empirical formula, and establishes a simple normal wave water acoustic propagation model according to the spectral intensity model of the wind-induced noise source. The simple normal wave underwater acoustic propagation model and the real part of the vertical component of the complex sound intensity of the noise are used to establish the cost function of the sea surface wind speed above the vector hydrophone, estimate the sea surface wind speed above the vector hydrophone, and centrally obtain the sea surface noise characteristics of the hydrophone, effectively filter out other noise interference, reduce the implementation difficulty and cost, and improve the observation and forecasting accuracy; finally, this application receives the typhoon eye pressure and typhoon eyewall radius information sent by the shore-based center via satellite in real time, combines the real-time wind speed estimated on the sea surface above the vector hydrophone, dynamically calculates the sea surface wind speed at different distances from the typhoon center, and sends the calculation results to the shore in real time via Iridium communication, so as to realize real-time observation and forecasting of typhoon wind field wind speed information, which is accurate, efficient and easy for engineering application.
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,而非限制本申请。It is to be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the present application.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.
图1是本申请基于单矢量水听器的台风风速实时观测方法的较佳实施例的流程图;FIG1 is a flow chart of a preferred embodiment of a method for real-time observation of typhoon wind speed based on a single vector hydrophone of the present application;
图2是本申请台风观测设备的结构示意图;FIG2 is a schematic diagram of the structure of the typhoon observation equipment of the present application;
图3是本申请基于单矢量水听器的台风风速实时观测方法的海洋空间二维坐标系示意图;FIG3 is a schematic diagram of a two-dimensional coordinate system in ocean space of the real-time observation method of typhoon wind speed based on a single vector hydrophone of the present application;
图4是本申请基于单矢量水听器的台风风速实时观测方法的台风风场海表面风速信息实时观测和预报计算步骤图;4 is a diagram of the real-time observation and forecast calculation steps of the typhoon wind speed real-time observation method of the present application based on a single vector hydrophone for typhoon wind speed information on the sea surface in a typhoon wind field;
图5是本申请基于单矢量水听器的台风风速实时观测方法的台风风场海表面气压仿真结果图;FIG5 is a diagram showing the simulation results of the sea surface pressure in a typhoon wind field according to the method for real-time observation of typhoon wind speed based on a single vector hydrophone of the present application;
图6是本申请基于单矢量水听器的台风风速实时观测方法的台风风场海表面风速仿真结果图;FIG6 is a diagram showing the simulation results of the sea surface wind speed in a typhoon wind field according to the real-time observation method of typhoon wind speed based on a single vector hydrophone of the present application;
图7是本申请基于单矢量水听器的台风风速实时观测系统的较佳实施例的原理示意图;FIG7 is a schematic diagram of the principle of a preferred embodiment of the typhoon wind speed real-time observation system based on a single vector hydrophone of the present application;
图8为本申请台风观测设备的较佳实施例的运行环境示意图。FIG8 is a schematic diagram of the operating environment of a preferred embodiment of the typhoon observation equipment of the present application.
其中,图中各附图标记:10、浮标组件;11、浮标本体;12、数据存储与处理模块;13、卫星通信模块;20、矢量水听器;30、光电复合缆;40、锚块;50、声学释放器;60、岸基中心;70、存储器;80、处理器;90、显示器;100、基于单矢量水听器的台风风速实时观测程序。Among them, the figure reference numerals are: 10, buoy assembly; 11, buoy body; 12, data storage and processing module; 13, satellite communication module; 20, vector hydrophone; 30, optoelectronic composite cable; 40, anchor block; 50, acoustic releaser; 60, shore-based center; 70, memory; 80, processor; 90, display; 100, typhoon wind speed real-time observation program based on single vector hydrophone.
具体实施方式Detailed ways
为使本领域的技术人员更好地理解本申请的技术方案,下面结合附图和具体实施方式对本申请所提供的基于单矢量水听器的台风风速实时观测方法、系统、台风观测设备及计算机可读存储介质做进一步详细描述。可以理解的是,所描述的实施例仅仅是本申请一部分实施例,而不是全部实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性的劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。In order to enable those skilled in the art to better understand the technical solution of the present application, the typhoon wind speed real-time observation method, system, typhoon observation equipment and computer-readable storage medium based on a single vector hydrophone provided by the present application are further described in detail below in conjunction with the accompanying drawings and specific implementation methods. It can be understood that the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in the field without making creative work are within the scope of protection of this application.
本申请中的术语“第一”、“第二”等是用于区别不同对象,而不是用于描述特定顺序。此外,术语“包括”和“具有”以及它们任何变形,意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元,而是可选地还包括没有列出的步骤或单元,或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。The terms "first", "second", etc. in this application are used to distinguish different objects, rather than to describe a specific order. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not limited to the listed steps or units, but optionally includes steps or units that are not listed, or optionally includes other steps or units inherent to these processes, methods, products or devices.
本申请提供一种基于单矢量水听器的台风风速实时观测方法、系统、台风观测设备及计算机可读存储介质,以解决现有技术中通过卫星遥感技术进行台风观测和预报,带来的观测和预报精度差,通过气象侦察飞机进行台风观测和预报带来的实施难度大,成本高,安全性低的问题。The present application provides a typhoon wind speed real-time observation method, system, typhoon observation equipment and computer-readable storage medium based on a single vector hydrophone to solve the problems in the prior art of poor observation and forecasting accuracy caused by typhoon observation and forecasting through satellite remote sensing technology, and the difficulty, high cost and low safety caused by typhoon observation and forecasting through meteorological reconnaissance aircraft.
请参阅图1至图6,图1是本申请基于单矢量水听器的台风风速实时观测方法的较佳实施例的流程图;图2是本申请台风观测设备的结构示意图;图3是本申请基于单矢量水听器的台风风速实时观测方法的海洋空间二维坐标系示意图;图4是本申请基于单矢量水听器的台风风速实时观测方法的台风风场海表面风速信息实时观测和预报计算步骤图;Please refer to Figures 1 to 6, Figure 1 is a flow chart of a preferred embodiment of the typhoon wind speed real-time observation method based on a single vector hydrophone of the present application; Figure 2 is a schematic diagram of the structure of the typhoon observation equipment of the present application; Figure 3 is a schematic diagram of the ocean space two-dimensional coordinate system of the typhoon wind speed real-time observation method based on a single vector hydrophone of the present application; Figure 4 is a diagram of the real-time observation and forecast calculation steps of the typhoon wind field sea surface wind speed information of the typhoon wind speed real-time observation method based on a single vector hydrophone of the present application;
图5是本申请基于单矢量水听器的台风风速实时观测方法的台风风场海表面气压仿真结果图;图6是本申请基于单矢量水听器的台风风速实时观测方法的台风风场海表面风速仿真结果图。Figure 5 is a diagram showing the simulation results of the sea surface air pressure in a typhoon wind field according to the real-time observation method of typhoon wind speed based on a single vector hydrophone of the present application; Figure 6 is a diagram showing the simulation results of the sea surface wind speed in a typhoon wind field according to the real-time observation method of typhoon wind speed based on a single vector hydrophone of the present application.
本申请提出一种基于单矢量水听器的台风风速实时观测方法,所述基于单矢量水听器的台风风速实时观测方法使用台风观测设备进行观测,如图2所示,台风观测设备包括浮标组件10、矢量水听器20、光电复合缆30和锚块40。其中,浮标组件10可设置在海面上,浮标组件10可包括数据存储与处理模块12和卫星通信模块13,数据存储与处理模块12通信连接于卫星通信模块13,卫星通信模块13通过卫星通信连接于岸基中心60;矢量水听器20可设置在浮标组件10的下方,矢量水听器20可用于监测和采集台风的噪声声场数据,噪声声场数包括但不限于噪声声压分量和垂直振速分量;光电复合缆30可连接在浮标组件10和矢量水听器20之间,光电复合缆30可用于将噪声声场数据传输至数据存储与处理模块12;锚块40可设置在海底,且锚块40可连接在矢量水听器20远离浮标组件10的一端,锚块40可用于系泊矢量水听器20;卫星通信模块13可用于与岸基中心60进行双向通信,数据存储与处理模块12可用于对所述噪声声场数据进行存储和计算,得到所述矢量水听器上方海表面的风速。The present application proposes a real-time observation method for typhoon wind speed based on a single vector hydrophone. The real-time observation method for typhoon wind speed based on a single vector hydrophone uses typhoon observation equipment for observation. As shown in FIG2 , the typhoon observation equipment includes a buoy assembly 10, a vector hydrophone 20, an optoelectronic composite cable 30 and an anchor block 40. Wherein, the buoy assembly 10 can be set on the sea surface, and the buoy assembly 10 may include a data storage and processing module 12 and a satellite communication module 13, the data storage and processing module 12 is communicatively connected to the satellite communication module 13, and the satellite communication module 13 is connected to the shore-based center 60 through satellite communication; the vector hydrophone 20 can be set below the buoy assembly 10, and the vector hydrophone 20 can be used to monitor and collect the noise field data of the typhoon, and the noise field data includes but is not limited to the noise sound pressure component and the vertical vibration velocity component; the optoelectronic composite cable 30 can be connected between the buoy assembly 10 and the vector hydrophone 20, and the optoelectronic composite cable 30 can be used to transmit the noise field data to the data storage and processing module 12; the anchor block 40 can be set on the seabed, and the anchor block 40 can be connected to the end of the vector hydrophone 20 away from the buoy assembly 10, and the anchor block 40 can be used to moor the vector hydrophone 20; the satellite communication module 13 can be used to communicate with the shore-based center 60 in a two-way manner, and the data storage and processing module 12 can be used to store and calculate the noise field data to obtain the wind speed of the sea surface above the vector hydrophone.
在一些实施例中,卫星通信模块13可以选择铱星通信,也可以选择北斗通信等其他通信方式,能够满足需要即可。In some embodiments, the satellite communication module 13 may select Iridium communication, or other communication methods such as Beidou communication, as long as the needs are met.
在一些实施例中,台风观测设备还可包括声学释放器50,声学释放器50可设置在锚块40和矢量水听器20之间,声学释放器50可用于矢量水听器20的二次回收。In some embodiments, the typhoon observation equipment may further include an acoustic releaser 50 . The acoustic releaser 50 may be disposed between the anchor block 40 and the vector hydrophone 20 . The acoustic releaser 50 may be used for secondary recovery of the vector hydrophone 20 .
在一些实施例中,浮标组件10包括浮标本体11,数据存储与处理模块12和卫星通信模块13均设置在浮标本体11内,浮标本体11用于支撑数据存储与处理模块12和卫星通信模块13。In some embodiments, the buoy assembly 10 includes a buoy body 11 , in which the data storage and processing module 12 and the satellite communication module 13 are both disposed. The buoy body 11 is used to support the data storage and processing module 12 and the satellite communication module 13 .
在一些实施例中,矢量水听器20可以设置为单个,也可根据需要设置有多个,多个矢量水听器20阵列设置,可以使得检测结果更准确。In some embodiments, the vector hydrophone 20 may be provided as a single one, or may be provided in plurality as required. An array arrangement of a plurality of vector hydrophones 20 may make the detection result more accurate.
在一些实施例中,矢量水听器20可以设计为自容式,在台风发生后回收观测设备,在岸基中心60利用采集数据按照同样的反演计算方法获得台风风速信息。In some embodiments, the vector hydrophone 20 can be designed to be self-contained, and the observation equipment can be recovered after a typhoon occurs. The shore-based center 60 uses the collected data to obtain typhoon wind speed information according to the same inversion calculation method.
在一些实施例中,如图1所示,所述基于单矢量水听器的台风风速实时观测方法包括步骤:In some embodiments, as shown in FIG1 , the method for real-time observation of typhoon wind speed based on a single vector hydrophone comprises the steps of:
步骤S100:将台风观测设备布放在台风风场外围海上,使用台风观测设备的单个矢量水听器采集台风环境下所述矢量水听器布放深度处的噪声声场数据。Step S100: deploying typhoon observation equipment on the sea outside the typhoon wind field, and using a single vector hydrophone of the typhoon observation equipment to collect noise sound field data at the deployment depth of the vector hydrophone in a typhoon environment.
具体的,将台风观测设备布放在台风风场外围远离台风中心的海上,也可以是距离台风中心第一距离的海上,其中,第一距离可以选择500到1000公里的范围内,使用单矢量水听器采集台风环境下矢量水听器布放深度处的噪声声场数据,通过光电复合缆将噪声声场数据传输至数据存储与处理模块,数据存储与处理模块将噪声声场数据传递给卫星通信模块,卫星通信模块将噪声声场数据传输给岸基中心。Specifically, the typhoon observation equipment is deployed on the sea outside the typhoon wind field, away from the typhoon center, or on the sea at a first distance from the typhoon center, wherein the first distance can be selected within the range of 500 to 1000 kilometers. A single vector hydrophone is used to collect noise sound field data at the depth of the vector hydrophone deployment in a typhoon environment, and the noise sound field data is transmitted to the data storage and processing module through an optoelectronic composite cable. The data storage and processing module transmits the noise sound field data to the satellite communication module, and the satellite communication module transmits the noise sound field data to the shore-based center.
在海洋空间建立二维坐标系,如图3所示,台风激发下的各个风成噪声源均匀分布在无限大海平面深度zs处,zr为水听器布放水深,r为水听器距离台风中心水平距离。A two-dimensional coordinate system is established in the ocean space, as shown in Figure 3. The various wind-generated noise sources excited by the typhoon are evenly distributed at the infinite sea surface depth zs , where zr is the water depth of the hydrophone and r is the horizontal distance between the hydrophone and the typhoon center.
步骤S200:根据所述噪声声场数据实时计算噪声复声强垂向分量实部。Step S200: Calculate the real part of the vertical component of the noise complex sound intensity in real time according to the noise sound field data.
具体的,根据噪声声场数据实时计算噪声复声强垂向分量实部,为建立矢量水听器上方海表面风速的代价函数做准备。Specifically, the real part of the vertical component of the complex sound intensity of the noise is calculated in real time according to the noise sound field data, in preparation for establishing a cost function of the sea surface wind speed above the vector hydrophone.
其中,所述噪声声场数据包括声压分量和垂直振速分量;所述步骤S200:根据所述噪声声场数据实时计算噪声复声强垂向分量实部,具体包括:The noise sound field data includes a sound pressure component and a vertical vibration velocity component; the step S200: calculating the real part of the vertical component of the noise complex sound intensity in real time according to the noise sound field data specifically includes:
获取第一时间长度的噪声声场数据,对所述噪声声场数据的所述声压分量和所述垂直振速分量进行傅里叶变换,得到第一时间长度的声压复频谱和垂直振速复频谱;Acquire noise sound field data of a first time length, perform Fourier transform on the sound pressure component and the vertical velocity component of the noise sound field data, and obtain a sound pressure complex spectrum and a vertical velocity complex spectrum of the first time length;
根据所述声压复频谱和所述垂直振速复频谱,计算所述矢量水听器的实测噪声复声强垂向分量实部;Calculating the real part of the vertical component of the complex sound intensity of the measured noise of the vector hydrophone according to the complex sound pressure spectrum and the vertical vibration velocity complex spectrum;
其中,所述实测噪声复声强垂向分量实部的计算公式为:The calculation formula of the real part of the vertical component of the measured noise complex sound intensity is:
其中,Iz实测(f,zr)为实测噪声复声强垂向分量实部;f为信号的中心频率;Re为复数量的实部;*为复共轭符号;<·>表示滑动窗平均周期图;zr为矢量水听器布放水深,r为矢量水听器距离台风中心水平距离,Sp(f,zr)为声压复频谱,Svz(f,zr)为垂直振速复频谱。 Among them, Izmeasured (f, zr ) is the real part of the vertical component of the measured noise complex sound intensity; f is the center frequency of the signal; Re is the real part of the complex quantity; * is the complex conjugate symbol; <·> represents the sliding window average periodogram; zr is the water depth of the vector hydrophone, r is the horizontal distance between the vector hydrophone and the typhoon center, Sp (f, zr ) is the complex spectrum of sound pressure, and Svz (f, zr ) is the complex spectrum of vertical vibration velocity.
具体的,实时获取矢量水听器位置处的噪声声场数据,噪声声场数据包括但不限于声压分量p(zr)和垂直振速分量vz(zr),获取第一时间长度tn的相同时间段不同分量噪声声场数据,对噪声声场数据的声压分量和垂直振速分量进行傅里叶变换,得到第一时间长度的声压复频谱和垂直振速复频谱,根据声压复频谱均值和垂直振速复频谱,计算矢量水听器的实测噪声复声强垂向分量实部为:Specifically, noise sound field data at the position of the vector hydrophone is acquired in real time, the noise sound field data including but not limited to a sound pressure component p(z r ) and a vertical vibration velocity component v z (z r ), noise sound field data of different components in the same time period of a first time length t n are acquired, the sound pressure component and the vertical vibration velocity component of the noise sound field data are Fourier transformed to obtain a sound pressure complex spectrum and a vertical vibration velocity complex spectrum of the first time length, and the real part of the vertical component of the complex sound intensity of the measured noise of the vector hydrophone is calculated according to the mean value of the complex sound pressure spectrum and the complex vertical vibration velocity spectrum:
其中,tn为傅里叶变化时间长度;f为信号的中心频率;Re为复数量的实部;*为复共轭符号;<·>表示滑动窗平均周期图。Where tn is the length of Fourier transform time; f is the center frequency of the signal; Re is the real part of the complex quantity; * is the complex conjugate symbol; <·> represents the sliding window average periodogram.
步骤S300:根据经验公式,建立风成噪声源谱强度模型。Step S300: Establishing a wind-induced noise source spectrum intensity model according to an empirical formula.
具体的,根据根据经验公式、海洋环境噪声与海表风速,建立风成噪声源谱强度模型,为建立简正波水声传播模型做准备。Specifically, based on empirical formulas, ocean environmental noise and sea surface wind speed, a spectral intensity model of wind-induced noise sources is established to prepare for the establishment of a simple normal wave underwater acoustic propagation model.
其中,所述步骤S300:根据经验公式,建立风成噪声源谱强度模型,具体包括:The step S300: establishing a wind-induced noise source spectrum intensity model according to an empirical formula, specifically includes:
根据海洋环境风成噪声源的频率和风速,以及风成噪声源经验公式,建立台风环境下海表面风成噪声源谱强度模型;According to the frequency and wind speed of wind noise sources in the marine environment and the empirical formula of wind noise sources, a spectrum intensity model of wind noise sources on the sea surface in a typhoon environment is established.
所述风成噪声源谱强度模型计算公式为:The calculation formula of the wind-generated noise source spectrum intensity model is:
SIw(U,f,zs)=C×fq×Un;SI w (U,f,z s ) = C × f q × U n ;
zs=cw/4f;z s = c w /4f;
其中,SIw(U,f,zs)为风成噪声源谱强度;U为海表面风速,U的参数寻优范围为4-16;zs为噪声源设置水深;cw为声源处声波在水体中的传播速度;C为第一参数;q为第二参数;n为第三参数,zs为台风激发下的各个风成噪声源均匀分布在大海平面的深度。Among them, SI w (U,f,z s ) is the spectral intensity of the wind-induced noise source; U is the sea surface wind speed, and the parameter optimization range of U is 4-16; z s is the water depth of the noise source; c w is the propagation speed of the sound wave in the water body at the sound source; C is the first parameter; q is the second parameter; n is the third parameter, and z s is the depth at which each wind-induced noise source excited by the typhoon is evenly distributed on the sea surface.
具体的,风成噪声源由频率和风速共同决定,根据台风环境下海洋环境风成噪声源的频率和风速,以及Wilson和Piggott风成噪声源经验公式,建立台风环境下海表面风成噪声源谱强度模型SIw(U,f,zs),其表达式为:Specifically, the wind-induced noise source is determined by the frequency and wind speed. According to the frequency and wind speed of the wind-induced noise source in the ocean environment under the typhoon environment, and the empirical formula of Wilson and Piggott wind-induced noise source, the spectrum intensity model SI w (U, f, z s ) of the wind-induced noise source on the sea surface under the typhoon environment is established, and its expression is:
SIw(U,f,zs)=C×fq×Un;SI w (U,f,z s ) = C × f q × U n ;
zs=cw/4f;z s = c w /4f;
其中,U为海表面风速,U的参数寻优范围为4-16,为待反演参数,单位为节(kn);zs为噪声源设置水深;cw为声源处声波在水体中的传播速度;C为第一参数;q为第二参数;n为第三参数,均为待求参数。Among them, U is the sea surface wind speed, the parameter optimization range of U is 4-16, it is the parameter to be inverted, and the unit is knot (kn); z s is the water depth of the noise source; c w is the propagation speed of the sound wave in the water body at the sound source; C is the first parameter; q is the second parameter; n is the third parameter, all of which are parameters to be determined.
步骤S400:根据所述风成噪声源谱强度模型,建立简正波水声传播模型。Step S400: establishing a simple normal wave underwater sound propagation model according to the wind-induced noise source spectrum intensity model.
具体的,根据风成噪声源谱强度模型,建立简正波水声传播模型,为建立矢量水听器上方海表面风速的代价函数做准备。Specifically, a simple normal wave underwater acoustic propagation model is established based on the spectral intensity model of wind-generated noise sources, in preparation for establishing a cost function of the sea surface wind speed above the vector hydrophone.
其中,所述步骤S400:根据所述风成噪声源谱强度模型,建立简正波水声传播模型,具体包括:The step S400: establishing a normal wave underwater acoustic propagation model according to the wind-induced noise source spectrum intensity model, specifically includes:
根据台风环境下海表面的所述风成噪声源谱强度模型和简正波理论,建立矢量水听器接收位置处台风激发下风成噪声源的水下简正波水声传播模型;According to the spectral intensity model of wind-induced noise source on the sea surface under typhoon environment and the simple normal wave theory, an underwater simple normal wave acoustic propagation model of wind-induced noise source under typhoon excitation at the receiving position of the vector hydrophone is established;
所述水下简正波水声传播模型的计算公式为:The calculation formula of the underwater normal wave acoustic propagation model is:
km=αm+iβm; km = αm + iβm ;
其中,Iz(f,zr)为水下简正波水声传播数据,ρ=1.0,ρ为水体密度,ψm为噪声源信号中心频率为f时第m个垂直模态特征函数,km为相应的模态水平传播特征值,αm为km的实部,βm为km的虚部,ψm和km为通过简正波声场模型仿真获取的参数,ψm′为特征函数导数,i为虚数单位,SIw(U,f,zs,zr)为水听器接收点上方海表面风速对应的风成噪声源谱强度,Ur为水听器上方海表面风速。Among them, Iz (f, zr ) is the underwater simple normal wave acoustic propagation data, ρ=1.0, ρ is the water density, ψm is the mth vertical mode characteristic function when the center frequency of the noise source signal is f, km is the corresponding modal horizontal propagation eigenvalue, αm is the real part of km , βm is the imaginary part of km , ψm and km are parameters obtained by simulation of the simple normal wave acoustic field model, ψm ′ is the derivative of the characteristic function, i is the imaginary unit, SIw (U,f, zs , zr ) is the spectral intensity of the wind-induced noise source corresponding to the sea surface wind speed above the hydrophone receiving point, and Ur is the sea surface wind speed above the hydrophone.
具体的,根据台风环境下海表面的建立的风成噪声源谱强度模型和简正波理论,假设各个风成噪声源为均匀分布的点源,分布在深度为zs处与海面平行的无限平面上,海洋环境垂直分层且不随距离变化,当接收水听器距离台风中心较远时,假设台风激发的各点噪声源互不相干,则在台风激发的风成噪声源作用下,根据简正波理论,建立矢量水听器接收位置处台风激发下风成噪声源的水下简正波水声传播模型:Specifically, according to the established spectral intensity model of wind-induced noise sources on the sea surface under typhoon environment and the simple normal wave theory, it is assumed that each wind-induced noise source is a uniformly distributed point source, distributed on an infinite plane parallel to the sea surface at a depth of zs , and the marine environment is vertically stratified and does not change with distance. When the receiving hydrophone is far away from the typhoon center, it is assumed that the point noise sources excited by the typhoon are irrelevant to each other. Then, under the action of the wind-induced noise source excited by the typhoon, according to the simple normal wave theory, an underwater simple normal wave acoustic propagation model of the wind-induced noise source under typhoon excitation at the receiving position of the vector hydrophone is established:
其中,ρ=1.0,ρ为水体密度;ψm为噪声源信号中心频率为f时第m个垂直模态特征函数,km=αm+iβm为相应的模态水平传播特征值,ψm和km通过简正波声场模型Kraken仿真计算获取;ψm′为特征函数导数;i为虚数单位;SIw(U,f,zs,zr)为水听器接收点上方海表面风速对应的风成噪声源谱强度;Ur为水听器上方海表面风速。Among them, ρ = 1.0, ρ is the water density; ψ m is the mth vertical mode characteristic function when the center frequency of the noise source signal is f, km = α m + iβ m is the corresponding modal horizontal propagation eigenvalue, ψ m and km are obtained by Kraken simulation calculation of the simple normal wave acoustic field model; ψ m ′ is the derivative of the characteristic function; i is the imaginary unit; SI w (U,f,z s ,z r ) is the spectral intensity of the wind-induced noise source corresponding to the sea surface wind speed above the hydrophone receiving point; U r is the sea surface wind speed above the hydrophone.
步骤S500:根据所述简正波水声传播模型和所述噪声复声强垂向分量实部,建立矢量水听器上方海表面风速的代价函数,估计所述矢量水听器上方海表面的风速。Step S500: establishing a cost function of the sea surface wind speed above the vector hydrophone according to the simple normal wave underwater acoustic propagation model and the real part of the vertical component of the noise complex sound intensity, and estimating the sea surface wind speed above the vector hydrophone.
具体的,根据简正波水声传播模型和噪声复声强垂向分量实部,建立矢量水听器上方海表面风速的代价函数,通过代价函数模型仿真结果与实测计算数据,利用多参数并行粒子群遗传优化算法反演计算关系模未知参量,估计矢量水听器上方海表面的风速。Specifically, according to the simple normal wave underwater acoustic propagation model and the real part of the vertical component of the complex sound intensity of the noise, a cost function of the sea surface wind speed above the vector hydrophone is established. Through the simulation results of the cost function model and the measured calculation data, a multi-parameter parallel particle swarm genetic optimization algorithm is used to inversely calculate the unknown parameters of the relational norm, and the wind speed on the sea surface above the vector hydrophone is estimated.
其中,所述步骤S500:根据所述简正波水声传播模型和所述噪声复声强垂向分量实部,建立矢量水听器上方海表面风速的代价函数,估计所述矢量水听器上方海表面的风速,具体包括:The step S500: establishing a cost function of the wind speed on the sea surface above the vector hydrophone according to the simple normal wave underwater acoustic propagation model and the real part of the vertical component of the complex sound intensity of the noise, and estimating the wind speed on the sea surface above the vector hydrophone, specifically includes:
根据所述简正波水声传播模型和所述噪声复声强垂向分量实部,利用多频联合反演的方式,采用多参数并行粒子群遗传优化算法,建立所述矢量水听器上方海表面风速的关于所述第一参数、所述第二参数和所述第三参数的代价函数,得到最优的所述第一参数、所述第二参数和所述第三参数;According to the simple normal wave underwater acoustic propagation model and the real part of the vertical component of the complex sound intensity of the noise, a cost function of the sea surface wind speed above the vector hydrophone with respect to the first parameter, the second parameter and the third parameter is established by using a multi-frequency joint inversion method and a multi-parameter parallel particle swarm genetic optimization algorithm to obtain the optimal first parameter, the second parameter and the third parameter;
所述代价函数的计算公式为:The calculation formula of the cost function is:
其中,E(C,q,n)为代价函数,J表示用于反演计算的频率的个数,n的参数寻优范围为2.5-3.5。Among them, E(C,q,n) is the cost function, J represents the number of frequencies used for inversion calculation, and the parameter optimization range of n is 2.5-3.5.
具体的,根据简正波水声传播模型和噪声复声强垂向分量实部,利用多频联合反演的方式,计算过程采用多参数并行粒子群遗传优化算法,计算J个频率下的累积最小平方和误差,得到目标函数E(C,q,n)的最小值,建立矢量水听器上方海表面风速的关于第一参数、第二参数和第三参数的代价函数,得到最优的第一参数、第二参数和第三参数;Specifically, according to the simple normal wave underwater acoustic propagation model and the real part of the vertical component of the complex sound intensity of the noise, the multi-frequency joint inversion method is used, and the calculation process adopts a multi-parameter parallel particle swarm genetic optimization algorithm to calculate the cumulative minimum square error under J frequencies, and obtain the minimum value of the objective function E(C, q, n). The cost function of the sea surface wind speed above the vector hydrophone with respect to the first parameter, the second parameter and the third parameter is established to obtain the optimal first parameter, the second parameter and the third parameter;
其中,代价函数的计算公式为:The calculation formula of the cost function is:
其中,J表示用于反演计算的频率的个数,n的参数寻优范围为2.5-3.5。Wherein, J represents the number of frequencies used for inversion calculation, and the parameter optimization range of n is 2.5-3.5.
步骤S600:通过卫星实时接收岸基中心发送的台风眼气压和台风眼壁半径信息,结合所述矢量水听器上方海表面估计的实时风速,动态计算距离台风中心不同距离上的海表面风速。Step S600: receiving the typhoon eye pressure and typhoon eyewall radius information sent by the shore-based center in real time via satellite, combining the real-time wind speed estimated on the sea surface above the vector hydrophone, and dynamically calculating the sea surface wind speed at different distances from the typhoon center.
具体的,通过卫星实时接收岸基中心发送的台风眼气压、台风眼壁半径信息,结合所述矢量水听器上方海表面估计的实时风速,动态计算距离台风中心不同距离上的海表面风速,将反演计算的台风风速通过铱星通信传输至岸基中心,实现台风风场风速信息的实时观测和预报,精准高效且易于工程应用;Specifically, the typhoon eye pressure and typhoon eye wall radius information sent by the shore-based center are received in real time by satellite, and the sea surface wind speed at different distances from the typhoon center is dynamically calculated in combination with the real-time wind speed estimated on the sea surface above the vector hydrophone. The inverse calculated typhoon wind speed is transmitted to the shore-based center via Iridium communication, so as to realize the real-time observation and forecast of the typhoon wind field wind speed information, which is accurate, efficient and easy for engineering application.
其中,所述步骤S600:通过卫星实时接收岸基中心发送的台风眼气压和台风眼壁半径信息,结合所述矢量水听器上方海表面估计的实时风速,动态计算距离台风中心不同距离上的海表面风速,具体包括:The step S600 includes receiving the typhoon eye pressure and typhoon eye wall radius information sent by the shore-based center in real time via a satellite, combining the real-time wind speed estimated on the sea surface above the vector hydrophone, and dynamically calculating the sea surface wind speed at different distances from the typhoon center, specifically including:
获取最优的所述第一参数、所述第二参数和所述第三参数,得到所述矢量水听器上方海表面的风速;Obtaining the optimal first parameter, second parameter and third parameter to obtain the wind speed of the sea surface above the vector hydrophone;
通过卫星实时接收所述岸基中心发送的台风眼气压和台风眼壁半径信息;receiving in real time via satellite the typhoon eye pressure and typhoon eye wall radius information sent by the shore-based center;
结合所述矢量水听器上方海表面估计的实时风速、所述台风眼气压和所述台风眼壁半径,使用台风风场计算模型,得到距离台风中心不同距离的海表面风速的动态计算表达式,得到所述海表面风速。The real-time wind speed estimated on the sea surface above the vector hydrophone, the typhoon eye pressure and the typhoon eyewall radius are combined, and a typhoon wind field calculation model is used to obtain a dynamic calculation expression for the sea surface wind speed at different distances from the typhoon center to obtain the sea surface wind speed.
其中,所述海表面风速的动态计算表达式为:The dynamic calculation expression of the sea surface wind speed is:
Rw=A1/B; Rw = A1 /B ;
其中,U(d)为矢量水听器处海表面风速,pn为台风眼气压,pc为台风风场外围气压,Rw为台风眼壁半径,e≈2.71828,为自然常数,ρa=1.15,为空气密度,A为第一中间参数,B为第二中间参数,d为被观测风速处与台风中心的水平距离。Among them, U(d) is the sea surface wind speed at the vector hydrophone, pn is the typhoon eye pressure, pc is the periphery of the typhoon wind field, Rw is the radius of the typhoon eyewall, e≈2.71828 is a natural constant, ρa =1.15 is the air density, A is the first intermediate parameter, B is the second intermediate parameter, and d is the horizontal distance between the observed wind speed and the typhoon center.
具体的,获取最优的第一参数、第二参数和第三参数,得到所述矢量水听器上方海表面的风速,通过卫星实时接收所述岸基中心发送的台风眼气压和台风眼区半径信息,根据矢量水听器上方海表面的风速、台风眼气压和台风眼区半径,使用Holland台风风场计算模型,给出距离台风中心r处矢量水听器处海表面风速实时计算表达式,得到海表面风速,如图5所示,台风风场海表面气压仿真结果图,如图6所示,台风风场海表面风速仿真结果图。Specifically, the optimal first parameter, second parameter and third parameter are obtained to obtain the wind speed on the sea surface above the vector hydrophone, and the typhoon eye pressure and typhoon eye radius information sent by the shore-based center are received in real time via satellite. According to the wind speed on the sea surface above the vector hydrophone, the typhoon eye pressure and the typhoon eye radius, the Holland typhoon wind field calculation model is used to give a real-time calculation expression for the sea surface wind speed at the vector hydrophone at a distance r from the typhoon center to obtain the sea surface wind speed, as shown in FIG5 , a simulation result diagram of the sea surface pressure in a typhoon wind field, and as shown in FIG6 , a simulation result diagram of the sea surface wind speed in a typhoon wind field.
矢量水听器处海表面风速实时计算表达式为:The real-time calculation expression of sea surface wind speed at the vector hydrophone is:
Rw=A1/B; Rw = A1 /B ;
其中,pn为台风眼气压,单位为帕(Pa),pc为台风风场外围气压,单位为帕(Pa),Rw为台风眼壁半径,e≈2.71828为自然常数,ρa=1.15为空气密度,A为第一中间参数,B为第二中间参数,均为待求参数,根据反演计算获取的Ur实时计算获取。Among them, pn is the typhoon eye pressure, in Pascal (Pa), pc is the typhoon wind field periphery pressure, in Pascal (Pa), Rw is the typhoon eyewall radius, e≈2.71828 is a natural constant, ρa =1.15 is the air density, A is the first intermediate parameter, and B is the second intermediate parameter, all of which are parameters to be determined and are obtained in real time based on Ur obtained by inversion calculation.
Rw、pn、pc和浮标距离台风中心的距离通过浮标安装的铱星通信从岸基中心远程获取,进而对距离台风中心d在0-1000km范围取值,获得不同距离台风激发下的海表面风速U(d),将反演计算的台风风速通过铱星通信传输至岸基中心,实现台风的实时观测和预报。R w , p n , pc and the distance from the buoy to the typhoon center are remotely acquired from the shore-based center via the Iridium satellite communication installed on the buoy. Then, the distance d from the typhoon center is taken in the range of 0-1000km to obtain the sea surface wind speed U(d) under typhoon excitation at different distances. The inverted typhoon wind speed is transmitted to the shore-based center via the Iridium satellite communication to realize real-time observation and forecast of the typhoon.
基于单矢量水听器的台风风速实时观测方法的原理,如图4所示,首先矢量水听器获取噪声声压分量和垂直振速分量,计算噪声复声强垂向分量实部;然后利用KRAKEN声场计算程序仿真获取简正波特征函数和特征值;再利用构建的关系模型计算海表面风成噪声源谱强度;再利用多参数并行粒子群遗传优化算法求取关系模型未知参量;再计算获取水听器上方海表面风速;再通过铱星通信接收岸基中心发送的台风气压、眼壁半径等参数;再动态计算距离台风中心不同距离海表风速;最后将计算结果通过铱星通信实时发送至岸基。The principle of the real-time observation method of typhoon wind speed based on a single vector hydrophone is shown in Figure 4. First, the vector hydrophone obtains the noise sound pressure component and the vertical vibration velocity component, and calculates the real part of the vertical component of the noise complex sound intensity; then the KRAKEN sound field calculation program is used to simulate and obtain the simple normal wave characteristic function and eigenvalue; then the constructed relationship model is used to calculate the spectrum intensity of the wind-induced noise source on the sea surface; then the multi-parameter parallel particle swarm genetic optimization algorithm is used to obtain the unknown parameters of the relationship model; then the sea surface wind speed above the hydrophone is calculated; then the typhoon pressure, eyewall radius and other parameters sent by the shore-based center are received through Iridium communication; then the sea surface wind speed at different distances from the typhoon center is dynamically calculated; finally, the calculation results are sent to the shore in real time through Iridium communication.
请参阅图7至图8,图7是本申请基于单矢量水听器的台风风速实时观测系统的较佳实施例的原理示意图;图8为本申请台风观测设备的较佳实施例的运行环境示意图。Please refer to Figures 7 to 8. Figure 7 is a schematic diagram of the principle of a preferred embodiment of the typhoon wind speed real-time observation system based on a single vector hydrophone of the present application; Figure 8 is a schematic diagram of the operating environment of a preferred embodiment of the typhoon observation equipment of the present application.
在一些实施例中,如图7所示,基于上述基于单矢量水听器的台风风速实时观测方法,本申请还提出一种基于单矢量水听器的台风风速实时观测系统,所述基于单矢量水听器的台风风速实时观测系统包括:In some embodiments, as shown in FIG. 7 , based on the above-mentioned typhoon wind speed real-time observation method based on a single vector hydrophone, the present application further proposes a typhoon wind speed real-time observation system based on a single vector hydrophone, and the typhoon wind speed real-time observation system based on a single vector hydrophone includes:
噪声声场数据采集模块51,用于将台风观测设备布放在台风风场外围海上,使用台风观测设备的单个矢量水听器采集台风环境下所述矢量水听器布放深度处的噪声声场数据;The noise sound field data collection module 51 is used to deploy the typhoon observation equipment on the sea outside the typhoon wind field, and use a single vector hydrophone of the typhoon observation equipment to collect the noise sound field data at the deployment depth of the vector hydrophone in the typhoon environment;
噪声复声强垂向分量实部计算模块52,用于根据所述噪声声场数据实时计算噪声复声强垂向分量实部;A noise complex sound intensity vertical component real part calculation module 52, used for calculating the noise complex sound intensity vertical component real part in real time according to the noise sound field data;
风成噪声源谱强度模型建立模块53,用于根据经验公式,建立风成噪声源谱强度模型;A wind-generated noise source spectrum intensity model establishment module 53 is used to establish a wind-generated noise source spectrum intensity model according to an empirical formula;
简正波水声传播模型建立模块54,用于根据所述风成噪声源谱强度模型,建立简正波水声传播模型;A simple normal wave underwater acoustic propagation model establishing module 54 is used to establish a simple normal wave underwater acoustic propagation model according to the wind-induced noise source spectrum intensity model;
矢量水听器上方海表面风速观测模块55,用于根据所述简正波水声传播模型和所述噪声复声强垂向分量实部,建立矢量水听器上方海表面风速的代价函数,估计所述矢量水听器上方海表面的风速;The sea surface wind speed observation module 55 above the vector hydrophone is used to establish a cost function of the sea surface wind speed above the vector hydrophone according to the simple normal wave underwater acoustic propagation model and the real part of the vertical component of the complex sound intensity of the noise, and estimate the sea surface wind speed above the vector hydrophone;
海表面风速实时观测模块56,用于通过卫星实时接收岸基中心发送的台风眼气压和台风眼壁半径信息,结合所述矢量水听器上方海表面估计的实时风速,动态计算距离台风中心不同距离上的海表面风速。The sea surface wind speed real-time observation module 56 is used to receive the typhoon eye pressure and typhoon eye wall radius information sent by the shore-based center in real time via satellite, and dynamically calculate the sea surface wind speed at different distances from the typhoon center in combination with the real-time wind speed estimated on the sea surface above the vector hydrophone.
在一些实施例中,如图8所示,基于上述基于单矢量水听器的台风风速实时观测方法和系统,本申请还相应提出一种台风观测设备,所述台风观测设备包括:单个矢量水听器20、卫星通信模块13、存储器70、处理器80、显示器90,图8仅示出了台风观测设备的部分组件,但是应理解的是,并不要求实施所有示出的组件,可以替代的实施更多或者更少的组件。In some embodiments, as shown in FIG8 , based on the above-mentioned typhoon wind speed real-time observation method and system based on a single vector hydrophone, the present application also proposes a typhoon observation device accordingly, and the typhoon observation device includes: a single vector hydrophone 20, a satellite communication module 13, a memory 70, a processor 80, and a display 90. FIG8 only shows some components of the typhoon observation device, but it should be understood that it is not required to implement all the components shown, and more or fewer components may be implemented instead.
所述存储器70在一些实施例中可以是所述台风观测设备的内部存储单元,例如台风观测设备的硬盘或内存。所述存储器70在另一些实施例中也可以是所述台风观测设备的外部存储设备,例如所述台风观测设备上配备的插接式硬盘,智能存储卡(Smart MediaCard,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,所述存储器70还可以既包括所述台风观测设备的内部存储单元也包括外部存储设备。所述存储器70用于存储安装于所述台风观测设备的应用软件及各类数据,例如所述安装台风观测设备的程序代码等。所述存储器70还可以用于暂时地存储已经输出或者将要输出的数据。In some embodiments, the memory 70 may be an internal storage unit of the typhoon observation device, such as a hard disk or memory of the typhoon observation device. In other embodiments, the memory 70 may also be an external storage device of the typhoon observation device, such as a plug-in hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), etc. equipped on the typhoon observation device. Further, the memory 70 may also include both an internal storage unit and an external storage device of the typhoon observation device. The memory 70 is used to store application software and various types of data installed in the typhoon observation device, such as the program code of the typhoon observation device. The memory 70 may also be used to temporarily store data that has been output or is to be output.
在一实施例中,存储器70上存储有基于单矢量水听器的台风风速实时观测程序100,该基于单矢量水听器的台风风速实时观测程序100可被处理器80所执行,从而实现本申请中基于单矢量水听器的台风风速实时观测方法。In one embodiment, a typhoon wind speed real-time observation program 100 based on a single vector hydrophone is stored in the memory 70, and the typhoon wind speed real-time observation program 100 based on a single vector hydrophone can be executed by the processor 80, thereby realizing the typhoon wind speed real-time observation method based on a single vector hydrophone in the present application.
所述处理器80在一些实施例中可以是一中央处理器(Central Processing Unit,CPU),微处理器或其他数据处理芯片,用于运行所述存储器70中存储的程序代码或处理数据,例如执行所述基于单矢量水听器的台风风速实时观测方法等。In some embodiments, the processor 80 may be a central processing unit (CPU), a microprocessor or other data processing chip, used to run the program code or process data stored in the memory 70, such as executing the real-time observation method of typhoon wind speed based on a single vector hydrophone.
所述显示器90在一些实施例中可以是LED显示器、液晶显示器、触控式液晶显示器以及OLED(Organic Light-Emitting Diode,有机发光二极管)触摸器等。所述显示器90用于显示在所述台风观测设备的信息以及用于显示可视化的用户界面。所述台风观测设备的部件10-30通过系统总线相互通信。In some embodiments, the display 90 may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, an OLED (Organic Light-Emitting Diode) touch device, etc. The display 90 is used to display information on the typhoon observation device and to display a visual user interface. The components 10-30 of the typhoon observation device communicate with each other via a system bus.
本申请还提出一种计算机可读存储介质,所述计算机可读存储介质存储有基于单矢量水听器的台风风速实时观测程序,所述基于单矢量水听器的台风风速实时观测程序被处理器执行时实现如上所述的基于单矢量水听器的台风风速实时观测方法的步骤。The present application also proposes a computer-readable storage medium, which stores a typhoon wind speed real-time observation program based on a single vector hydrophone. When the typhoon wind speed real-time observation program based on a single vector hydrophone is executed by a processor, the steps of the typhoon wind speed real-time observation method based on a single vector hydrophone as described above are implemented.
综上所述,本申请通过将所述台风观测设备布放在台风风场外围海上,使用台风观测设备的单矢量水听器采集台风环境下矢量水听器布放深度处的噪声声场数据,并利用岸基遥感卫星数据实时接收噪声数据,动态优化台风风速模型参数,实现不同距离台风风场风速的动态计算和实时回传,有效提高台风风速观测和预报的精度及效率;其次,本申请根据噪声声场数据实时计算噪声复声强垂向分量实部,根据经验公式,建立风成噪声源谱强度模型,根据风成噪声源谱强度模型,建立简正波水声传播模型,根据简正波水声传播模型和噪声复声强垂向分量实部,建立矢量水听器上方海表面风速的代价函数,估计矢量水听器上方海表面的风速,集中获取水听器海表面噪声特性,有效滤除了其他噪声干扰,降低了实施难度和成本,提高了观测和预报精度;最后,本申请通过卫星通过实时接收岸基中心发送的台风眼气压、台风眼壁半径信息,结合矢量水听器上方海表面估计的实时风速,动态计算距离台风中心不同距离上的海表面风速,并将计算结果通过铱星通信实时发送至岸基,实现台风风场风速信息的实时观测和预报,精准高效且易于工程应用。In summary, the present application deploys the typhoon observation equipment on the sea outside the typhoon wind field, uses the single vector hydrophone of the typhoon observation equipment to collect the noise sound field data at the deployment depth of the vector hydrophone in the typhoon environment, and uses the shore-based remote sensing satellite data to receive the noise data in real time, dynamically optimizes the typhoon wind speed model parameters, realizes the dynamic calculation and real-time feedback of the wind speed in the typhoon wind field at different distances, and effectively improves the accuracy and efficiency of typhoon wind speed observation and forecasting; secondly, the present application calculates the real part of the vertical component of the complex sound intensity of the noise in real time according to the noise sound field data, establishes a wind-induced noise source spectrum intensity model according to the empirical formula, and establishes a simple normal wave water acoustic propagation model according to the wind-induced noise source spectrum intensity model. The model and the real part of the vertical component of the complex sound intensity of the noise are used to establish a cost function of the sea surface wind speed above the vector hydrophone, estimate the wind speed on the sea surface above the vector hydrophone, centrally obtain the sea surface noise characteristics of the hydrophone, effectively filter out other noise interference, reduce the implementation difficulty and cost, and improve the observation and forecast accuracy; finally, the application uses a satellite to receive the typhoon eye pressure and typhoon eyewall radius information sent by the shore-based center in real time, combined with the real-time wind speed estimated on the sea surface above the vector hydrophone, dynamically calculates the sea surface wind speed at different distances from the typhoon center, and sends the calculation results to the shore in real time through Iridium communication, so as to realize real-time observation and forecasting of typhoon wind field wind speed information, which is accurate, efficient and easy for engineering application.
需要说明的是,本申请实施例中介绍的多种可选的实施方式,彼此可以相互结合实现,也可以单独实现,对此本申请实施例不作限定。It should be noted that the various optional implementations introduced in the embodiments of the present application can be implemented in combination with each other or can be implemented separately, and the embodiments of the present application are not limited to this.
在本申请的描述中,需要理解的是,术语“上”、“下”、“左”、“右”等指示方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以及特定的方位构造和操作。因此,不能理解为对本申请的限制。此外,“第一”、“第二”仅由于描述目的,且不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。因此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者多个该特征。本申请的描述中,除非另有说明,“多个”的含义是两个或两个以上。In the description of the present application, it should be understood that the terms "upper", "lower", "left", "right", etc. indicate orientations or positional relationships based on the orientations or positional relationships shown in the accompanying drawings, and are only for the convenience of describing the present application and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, and a specific orientation structure and operation. Therefore, it cannot be understood as a limitation on the present application. In addition, "first" and "second" are only for descriptive purposes, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Therefore, the features defined as "first" and "second" may explicitly or implicitly include one or more of the features. In the description of the present application, unless otherwise specified, "multiple" means two or more.
在本申请的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”等应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接连接,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本申请中的具体含义。In the description of this application, it should be noted that, unless otherwise clearly specified and limited, the terms "installed", "connected", "connected", etc. should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or an indirect connection through an intermediate medium, or it can be the internal communication of two components. For ordinary technicians in this field, the specific meanings of the above terms in this application can be understood according to specific circumstances.
上述实施例是参考附图来描述的,其他不同的形式和实施例也是可行而不偏离本申请的原理,因此本申请不应被建构成为在此所提出实施例的限制。更确切地说,这些实施例被提供以使得本申请会是完善又完整,且会将本申请范围传达给本领域技术人员。在附图中,组件尺寸及相对尺寸也许基于清晰起见而被夸大。在此所使用的术语只是基于描述特定实施例目的,并无意成为限制。术语“包含”及/或“包括”在使用于本说明书时,表示所述特征、整数、构件及/或组件的存在,但不排除一或更多其他特征整数、构件、组件及/或其族群的存在或增加。除非另有所示,陈述时,数值范围包含该范围的上下限及其间的任何子范围。The above embodiments are described with reference to the accompanying drawings, and other different forms and embodiments are also feasible without departing from the principles of the present application, so the present application should not be constructed as a limitation of the embodiments proposed herein. More specifically, these embodiments are provided so that the present application will be perfect and complete, and the scope of the present application will be conveyed to those skilled in the art. In the accompanying drawings, the component sizes and relative sizes may be exaggerated for clarity. The terms used herein are only based on the purpose of describing specific embodiments and are not intended to be limiting. The terms "comprising" and/or "including" when used in this specification indicate the presence of the features, integers, components and/or components, but do not exclude the presence or increase of one or more other feature integers, components, components and/or their groups. Unless otherwise indicated, when stated, the numerical range includes the upper and lower limits of the range and any sub-ranges therebetween.
以上所述仅为本申请的部分实施例,并非因此限制本申请的保护范围,凡是利用本申请说明书及附图内容所作的等效装置或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本申请的专利保护范围内。The above descriptions are only some embodiments of the present application, and do not limit the protection scope of the present application. Any equivalent device or equivalent process transformation made using the contents of the present application specification and drawings, or directly or indirectly used in other related technical fields, are also included in the patent protection scope of the present application.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410158331.1A CN117825745B (en) | 2024-02-04 | 2024-02-04 | A real-time observation method and device for typhoon wind speed based on single vector hydrophone |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410158331.1A CN117825745B (en) | 2024-02-04 | 2024-02-04 | A real-time observation method and device for typhoon wind speed based on single vector hydrophone |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117825745A true CN117825745A (en) | 2024-04-05 |
CN117825745B CN117825745B (en) | 2024-07-26 |
Family
ID=90505932
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410158331.1A Active CN117825745B (en) | 2024-02-04 | 2024-02-04 | A real-time observation method and device for typhoon wind speed based on single vector hydrophone |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117825745B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117991254A (en) * | 2024-04-07 | 2024-05-07 | 南京信息工程大学 | A method and system for estimating typhoon moving velocity vector based on synthetic aperture radar in quasi-real-time monitoring |
CN118393607A (en) * | 2024-06-25 | 2024-07-26 | 南方海洋科学与工程广东省实验室(广州) | Typhoon intensity monitoring method, system and terminal based on optical fiber hydrophone buoy |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2376653C1 (en) * | 2008-05-20 | 2009-12-20 | Виктор Сергеевич Аносов | Device of hydrometeorological surveys of water area of sea polygon |
CN109443516A (en) * | 2018-12-25 | 2019-03-08 | 西北工业大学 | A kind of passive acquisition methods of Bottom sound speed based on the vertical vibration velocity signal of noise field |
CN109815942A (en) * | 2019-03-18 | 2019-05-28 | 西北工业大学 | Normal wave feature extraction method based on marine environmental noise signal |
CN115951328A (en) * | 2023-03-10 | 2023-04-11 | 中国人民解放军国防科技大学 | Wind speed estimation method and device for wind lidar based on probability density constraints |
CN116203651A (en) * | 2023-01-18 | 2023-06-02 | 浙江大学 | Typhoon wind speed inversion method based on extremely low frequency noise |
WO2023202008A1 (en) * | 2022-04-19 | 2023-10-26 | 中国科学院声学研究所 | Marine environment noise forecasting method, computer device, and storage medium |
-
2024
- 2024-02-04 CN CN202410158331.1A patent/CN117825745B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2376653C1 (en) * | 2008-05-20 | 2009-12-20 | Виктор Сергеевич Аносов | Device of hydrometeorological surveys of water area of sea polygon |
CN109443516A (en) * | 2018-12-25 | 2019-03-08 | 西北工业大学 | A kind of passive acquisition methods of Bottom sound speed based on the vertical vibration velocity signal of noise field |
CN109815942A (en) * | 2019-03-18 | 2019-05-28 | 西北工业大学 | Normal wave feature extraction method based on marine environmental noise signal |
WO2023202008A1 (en) * | 2022-04-19 | 2023-10-26 | 中国科学院声学研究所 | Marine environment noise forecasting method, computer device, and storage medium |
CN116203651A (en) * | 2023-01-18 | 2023-06-02 | 浙江大学 | Typhoon wind speed inversion method based on extremely low frequency noise |
CN115951328A (en) * | 2023-03-10 | 2023-04-11 | 中国人民解放军国防科技大学 | Wind speed estimation method and device for wind lidar based on probability density constraints |
Non-Patent Citations (6)
Title |
---|
徐朋豪 等: "利用噪声谱进行风浪声场特征反演", 测控技术, no. 08, 31 August 2016 (2016-08-31), pages 41 - 44 * |
李跃金 等: "海底原位分层声学测量模拟装置的设计与实验", 广东海洋大学学报, no. 1, 31 January 2024 (2024-01-31), pages 111 - 117 * |
李风华 等: "台风激发水下噪声场的建模及其在台风风速反演中的应用", 声学学报, no. 05, 30 September 2016 (2016-09-30), pages 210 - 217 * |
李风华;刘姗琪;王琰;: "台风激发水下噪声场的建模及其在台风风速反演中的应用", 声学学报, no. 05, 15 September 2016 (2016-09-15) * |
王琰 等: "深海海洋环境噪声垂直相关性研究", 声学技术, no. 02, 30 April 2016 (2016-04-30), pages 25 - 29 * |
邹吉武 等: "海洋矢量噪声场试验研究", 哈尔滨工程大学学报, no. 01, 31 January 2011 (2011-01-31), pages 16 - 20 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117991254A (en) * | 2024-04-07 | 2024-05-07 | 南京信息工程大学 | A method and system for estimating typhoon moving velocity vector based on synthetic aperture radar in quasi-real-time monitoring |
CN117991254B (en) * | 2024-04-07 | 2024-06-21 | 南京信息工程大学 | Estimation method and system for typhoon moving speed vector based on synthetic aperture radar quasi-real-time monitoring |
CN118393607A (en) * | 2024-06-25 | 2024-07-26 | 南方海洋科学与工程广东省实验室(广州) | Typhoon intensity monitoring method, system and terminal based on optical fiber hydrophone buoy |
Also Published As
Publication number | Publication date |
---|---|
CN117825745B (en) | 2024-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN117825745B (en) | A real-time observation method and device for typhoon wind speed based on single vector hydrophone | |
US20240093989A1 (en) | System and method for wave prediction | |
CN102446367B (en) | Method for constructing three-dimensional terrain vector model based on multi-beam sonar submarine measurement data | |
CN114067530B (en) | Marine information perception and early warning method, system and storage medium based on optical fiber sensing | |
Rabault et al. | Measurements of waves in landfast ice using inertial motion units | |
CN109556829B (en) | A real-time detection method of ocean internal wave acoustics based on line spectrum recognition | |
CN108414002A (en) | A kind of environmental parameter test device of offshore wind farm unit | |
CN116628559A (en) | Underwater big data calculation comprehensive experiment classification system and training method of classification model | |
CN115128617A (en) | High-precision seafloor imaging method suitable for deep-sea mineral resource exploration areas | |
CN110068309A (en) | One kind being based on New High Order disperse relationship deep water offshore measurement method and device | |
CN209638791U (en) | A kind of submarine pipeline leakage monitoring system | |
CN201672919U (en) | A mid-range and long-range automatic self-adaptive wave measuring buoy | |
CN113156393A (en) | Airborne laser sounding broken wind wave sea surface model construction method | |
CN109426253A (en) | A kind of control method and system for determining region sail unmanned boat at a distance | |
KR102558955B1 (en) | System and method for measuring water depth using radar image processing | |
KR102111637B1 (en) | System and method for predicting wave using observation data | |
CN116958435A (en) | Wave information inversion system based on X-band radar image | |
CN115220026A (en) | A method for cross-ice localization in matching field based on horizontal array of ice surface | |
CN113126029B (en) | A multi-sensor pulsed sound source localization method suitable for deep-sea reliable sound path environment | |
US20230167796A1 (en) | Method for predicting a characteristic resulting from a swell on the basis of a spectral model of the swell | |
CN1687705A (en) | System for observing wave in waterway | |
Ikpekha | Underwater acoustics propagation analysis and modelling of sound emitting devices | |
CN211714048U (en) | Transmission tower foundation dimension measurement transmission system | |
Watanabe et al. | Data assimilation of the stereo reconstructed wave fields to a nonlinear phase resolved wave model | |
CN217501868U (en) | Floating wind power generation system |
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 |