CN110146895B - 基于倒置式多波束回声仪的声速剖面反演方法 - Google Patents
基于倒置式多波束回声仪的声速剖面反演方法 Download PDFInfo
- Publication number
- CN110146895B CN110146895B CN201910410164.4A CN201910410164A CN110146895B CN 110146895 B CN110146895 B CN 110146895B CN 201910410164 A CN201910410164 A CN 201910410164A CN 110146895 B CN110146895 B CN 110146895B
- Authority
- CN
- China
- Prior art keywords
- sound velocity
- sound
- formula
- echo
- depth
- 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
Links
Images
Classifications
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H5/00—Measuring propagation velocity of ultrasonic, sonic or infrasonic waves, e.g. of pressure waves
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于倒置式多波束回声仪的声速剖面反演方法,该方法包括以下步骤:将多波束回声仪倒置式安装于水下潜标或者固定于水底,通过发射换能器阵向水面发射波束,通过接收换能器阵接收回波信号;多波束回声仪根据接收到的回波信号获得回波的到达角和到达时间;根据声速剖面先验信息求解EOF,得到声速剖面的降维基函数描述方法;结合EOF、声线跟踪算法、表面声速以及多波束数据,建立优化模型;根据建立的优化模型、接收到的回波的到达时间和到达角以及表面声速,利用优化算法获得测量区域声速剖面的估计结果。进一步,利用声速剖面的估计值计算测量区域的水温剖面。本发明能够快速、准确地跟踪声速剖面及温度剖面的起伏变化。
Description
技术领域
本发明涉及的是水声测量和信号信息处理技术领域,具体涉及一种基于倒置式多波束回声仪的声速剖面反演方法。
背景技术
典型的海洋***是一个由众多机理且相互作用的物理、声学、光学、生物、化学及地质过程组成的复杂***,其内部存在着形式多样的动力过程。这些动力过程以及它们之间的耦合,加上与海床的相互作用,导致海洋内部的不均匀性以及时间与空间上的显著变化。在诸多环境参数中,声速剖面具有明显的时变特性,对声信号传播有重要影响,进而影响水声通信、定位、测绘等声纳***的工作性能。海洋环境的有效测量技术一直受到关注,是当今水声研究的一个主要方向。
目前获得声速剖面的方式主要有两种,即测量方法与反演方法。直接测量法利用温盐深仪(CTD)通过逐点测量的方式获取声速剖面。这种方法虽然具有非常高的测量精度,但直接测量不仅费时费力,而且无法获取大范围海域动态变化的观测数据,难以满足实际应用的需求。近年来,声学反演技术通过将声传播模型与观测数据相结合,为获取感兴趣的环境参数提供了一种长期、快速、高效的方法,与传统方法相比具有很大的优势。
多波束回声仪在测量过程中,利用特殊的发射和接收波束结构,在一个较宽的范围内形成许多较窄的测量波束,从而一次性获得一个较宽范围内水体回声信息。多波速回声仪采用正交的“T”型的发射接收阵结构,发射阵沿船航行方向布置,形成纵向(沿航行方向)较窄横向(垂直航行方向)较宽的发射波束,形成一个测量条带;接收阵则横向布置,通过接收波束形成处理,形成许多横向较窄而纵向较宽的波束,接收波束与发射波束相交形成许多测量点(脚印),最后利用各个脚印回波的到达时间与角度计算出整个测量值。近几十年来,利用多波束回声仪的测深技术得到了快速发展,已成为水下勘测的主要手段之一,在海洋工程建设、资源开发、海图绘制、海洋研究与权益维护等方面起到了极为重要的作用。应用于测深时为保证测量精度,数据需经声速剖面修正,同时数据中也包含了声速剖面信息。利用多波束回声仪数据反演声速及温度等水文参数,在改进地形测量精度的同时获取水体环境信息,体现了垂直方向声场-动力的深度耦合关系。
倒置式回声仪长期固定布防于海底,通过测量声波在垂直方向传播时间的变化,反演温度、盐度等物理海洋学参量,可应用于大洋环流、中尺度涡旋、内波等深海观测。目前,倒置式回声仪采用单波束方式,通过最严经验模态方法(GEM-Gravest EmpiricalMode)得到海洋动力参数的剖面分布。在深海情况下,单波束方式的倒置式回声仪测量误差较大。同时,考虑到浅海复杂的温跃层结构与时空变化环境,单波束方式的倒置式回声仪应用受限。
寻找能够准确、快速地获得声速剖面时变特性的反演方法成为目前研究的重点。
发明内容
本发明的目的在于针对现有技术的不足,提供了一种基于倒置式多波束回声仪的声速剖面反演方法和基于倒置式多波束回声仪的水温剖面反演方法。
为实现上述目的,本发明所采取的其中一种技术方案是:基于倒置式多波束回声仪的声速剖面反演方法包括以下步骤:
(1)将多波束回声仪倒置式安装于水下潜标或者固定于水底,通过发射换能器阵向水面发射波束,通过接收换能器阵接收回波信号;
(2)多波束回声仪根据接收到的回波信号,获得回波的到达角和到达时间;
(3)根据声速剖面先验信息,求解经验正交函数(Empirical OrthogonalFunction-EOF),得到声速剖面的降维基函数描述方法:
某一时刻t,声速剖面可以表示为:
式(1)中,z表示深度;c0(z)表示深度z处的平均声速,Δc(z,t)表示扰动量;i表示EOF的阶数,1≤i≤I,I表示EOF的总阶数,αi(t)表示第i阶EOF系数,ψi(z)表示第i阶EOF;
(4)结合EOF、声线跟踪算法、表面声速以及多波束数据,建立优化模型:
式(2)中,K为权重系数,取[0,1];Toameas(t)为t时刻所有波束到达时间的实际测量数据,Toaest(t)为t时刻所有波束到达时间的估计值,c0meas(t)为t时刻多波束回声仪接收换能器阵所在位置声速的测量数据,c0est(t)为t时刻多波束回声仪接收换能器阵所在位置声速的估计值;为测量算子,可根据声线跟踪算法进行计算,Doameas(t)为t时刻所有波束的到达方向测量值;
(5)根据步骤(4)建立的优化模型、步骤(2)中接收到的回波的到达角和到达时间以及接收换能器阵所在位置声速的测量值,利用优化算法获得EOF系数的反演值,进而利用式(1)获得声速剖面的估计结果。
进一步地,所述步骤(5)中,所述优化算法为模拟退火算法、爬山算法、粒子群算法、遗传算法、蚁群算法。
本发明所采取的另一种技术方案是:基于倒置式多波束回声仪的水温剖面反演方法包括如下步骤:
(1)将多波束回声仪倒置式安装于水下潜标或者固定于水底,通过发射换能器阵向水面发射波束,通过接收换能器阵接收回波信号;
(2)多波束回声仪根据接收到的回波信号,获得回波的到达角和到达时间;
(3)根据声速剖面先验信息,求解经验正交函数(Empirical OrthogonalFunction-EOF),得到声速剖面的降维基函数描述方法:
某一时刻t,声速剖面可以表示为:
式(3)中,z表示深度;c0(z)表示深度z处的平均声速,Δc(z,t)表示扰动量;i表示EOF的阶数,1≤i≤I,I表示EOF的总阶数,αi(t)表示第i阶EOF系数,ψi(z)表示第i阶EOF;
(4)结合EOF、声线跟踪算法、表面声速以及多波束数据,建立优化模型:
式(4)中,K为权重系数,取[0,1];Toameas(t)为t时刻所有波束到达时间的实际测量数据,Toaest(t)为t时刻所有波束到达时间的估计值,c0meas(t)为t时刻多波束回声仪接收换能器阵所在位置声速的测量数据,c0est(t)为t时刻多波束回声仪接收换能器阵所在位置声速的估计值;为测量算子,可根据声线跟踪算法进行计算,Doameas(t)为t时刻所有波束的到达方向测量值;
(5)根据步骤(4)建立的优化模型、步骤(2)中接收到的回波的到达角和到达时间以及接收换能器阵所在位置声速的测量值,利用优化算法获得EOF系数的反演值,进而利用式(3)获得声速剖面的估计结果。
(6)根据步骤(5)的声速剖面的估计值计算测量区域的水温剖面。
进一步地,所述步骤(5)中,所述优化算法为模拟退火算法、爬山算法、粒子群算法、遗传算法、蚁群算法。
进一步地,所述步骤(6)中,采用Mackenzie声速公式、Chen-Millero-Li声速公式、Dell Grosso声速公式、Leroy声速公式、Wilson精确声速公式、Wilson简化声速公式或EM分层简化声速公式计算测量区域的水温剖面。
与现有技术相比,本发明的有益效果是:
(1)本发明方法直接利用倒置式多波束回声仪的测量数据(到达时间与到达角度)实现测量区域声速剖面估计,进一步可以获得测量区域的水温剖面,相比于直接投放CTD的方法更为快速;
(2)本发明采用倒置式多波束回声仪能够同时获得数十甚至上百个窄波束,相比于单波束回声仪,具有更高的海洋环境估计精度;
(3)本发明采用倒置式多波束回声仪实现环境参数(声速剖面及水温剖面)估计,可以坐底安装或安装于潜标、浮标上,在满足安全性的同时,可以实现对海洋环境的大范围、长时序、动态监测;
(4)本发明采用倒置式多波束回声仪安装于潜标上,可以实现潜标安装位置至水面环境参数估计,解决传统潜标无法测量海水表层水文参数的问题,对于潜标观测应用有极大助益。
附图说明
图1是多波束回声仪的结构框图;
图2是本发明反演方法的流程示意图;
图3是常声速声线跟踪法示意图;
图4是利用本发明反演方法所获得的声速剖面估计结果与真实值的比较图;
图5是利用本发明反演方法所获得的水温剖面估计结果与真实值的比较图。
具体实施方式
下面结合附图和实施例对本发明作进一步具体描述:
如图1所示,多波束回声仪包括:
发射电路:用于产生发射电信号。
发射换能器阵:用于将产生的发射电信号转化为声信号发射。
接收换能器阵:用于接收经过水体及水面反射、散射的回波,并将回波声信号转化为电信号。
接收电路:用于将接收换能器阵转换得到的电信号进行放大、滤波等预处理。
姿态传感器:用于实现多波束回声仪的姿态测定。
压力传感器:用于实现多波束回声仪的压力测定。
表面声速头:用于实现多波束回声仪安装位置的声速测定。
***控制与信号处理***:用于实现信号的高速采集和实时处理,并将采集的数据和处理结果通过网络传给显控与后处理平台,其中采集数据包括姿态传感器数据、压力传感器数据及表面声速数据,处理结果包括回波的到达角和到达时间。
显控与后处理平台:用于下达指令给***控制与信号处理***,并接收其处理结果和采集数据;在此基础上利用所建立的优化模型,实现声速剖面的反演,进一步,利用声速剖面的估计值计算测量区域的水温剖面。
本发明基于倒置式多波束回声仪的声速剖面反演方法的流程如图2所示,包括如下步骤:
(一)将多波束回声仪倒置式安装于水下潜标或者固定于水底,通过发射换能器阵向水面发射波束,通过接收换能器阵接收回波信号;
(二)多波束回声仪根据接收到的回波信号,获得回波的到达角和到达时间;
(三)根据声速剖面先验信息,求解EOF,得到声速剖面的降维基函数描述方法:
EOF可以根据测量海域的已知声速剖面样本集合进行计算。假定有M次实测声速剖面数据,如
式(5)中zi表示深度值。则可得到声速的确定量为
协方差矩阵为
(四)结合EOF、声线跟踪算法、表面声速以及多波束数据,建立优化模型:
由于海水的不均匀性以及各向异性,声波在穿越不同声速层时会发生折射,从而导致声线弯曲,如图3所示。在非均匀介质条件下,假设声速只在垂直方向上变化,而在水平方向上是均匀的,声线传播遵循Snell定理:
式(8)中,c0表示换能器基阵位置处的声速,称之为表面声速;θ表示波束回波到达接收阵时与垂直方向的夹角;c(z)表示深度z处的声速,θ(z)表示深度z处声线切线与垂直方向的夹角。假设在很小的一段时间dt内,声线传播的水平距离dr、垂直距离dz以及路径长度ds满足
则自换能器基阵发出的声线到达深度z处的传播时间为:
且声线传播的水平距离为:
式(11)中,z0为换能器基阵深度。
由于海洋环境复杂多变,难以采用简单的连续函数表示声速剖面c(z)。根据这些实测声速剖面测量点的深度将整个水体分成若干层,如图3所示,并采用层内等梯度假设进行声线跟踪计算。
首先,将水体沿深度方向分为N层,假设分层处的深度与声速值分别为z=[z0,z1,…,zN]T和c=[c0,c1,…,cN]T,并假设层内声速常梯度,即:
c(z)=ci+gi(z-zi). (12)
式(12)中,i表示层数(i=0,2,…,N-1),ci表示第i层分界处的声速值,zi表示第i层的深度,gi表示第i层内的声速梯度,根据式(10)可知,在层内常梯度假设下,第i层内声线的双程传播时间可表示为:
根据式(11)可知,第i层声线传播的水平距离可表示为:
则整个深度范围内,声线总的传播时间、水平距离与深度的关系可表示为:
(五)根据上述建立的优化模型,通过网络将回波的到达角与到达时间数据以及采集的姿态、压力、表面声速数据传给显控与后处理平台,利用优化算法获得声速剖面基函数EOF的系数反演。
通常情况下,可采用模拟退火算法、爬山算法、粒子群算法、遗传算法、蚁群算法等优化算法实现待估参数的估计。图4给出了Matlab软件自带遗传算法的仿真实例。EOFs阶数为3阶。遗传算法的参数估计下限为[-100 -100 -100],估计上限为[100 100 100],种群数为200,迭代次数为100。由图4可以看出,采用本发明的方法估计出的声速剖面与真实剖面十分接近,整个剖面的均方误差为0.9726m/s,可见本发明方法估计效果良好。
(六)根据步骤(五)的声速剖面的估计值计算测量区域的水温剖面。
采用适用于测量区域的声速经验公式,根据步骤(五)中的声速剖面的估计结果,可进一步计算测量区域的水温剖面。目前被普遍认可的声速经验公式主要有Mackenzie声速公式、Chen-Millero-Li声速公式、Dell Grosso声速公式、Leroy声速公式、Wilson精确声速公式、Wilson简化声速公式或EM分层简化声速公式。其中,EM分层简化声速公式结构简单,计算精度高,且能够适用于海水盐度、深度或温度变化范围大的测量区域,是本发明的优选实施方式。
以下以EM分层简化声速公式为例进行水温剖面计算。
EM分层简化声速公式可以表示为:
表层声速计算公式为
淡水中深度达到200米、海水中深度达到1000米时声速计算公式为
c(T,z,S)=c(T,0,S)+16.5z. (18)
淡水中深度达到2000米、海水中深度达到11000米时声速计算公式为
深度大于5000米时,声速计算公式应考虑纬度改正:
根据压力传感器数据即可知测量***安放位置的深度,根据式(17)-(20)的海水深度适用范围,选择其中合适的一个进行水温剖面估计。由于仿真实例中***安放在200米深度左右,所以选取式(17)进行测量区域水温剖面的估计,如图5所示。由图4可以看出,采用本发明的方法估计出的水温剖面与真实剖面十分接近,可见本发明方法估计效果良好。
上述实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
Claims (7)
1.一种基于倒置式多波束回声仪的声速剖面反演方法,其特征是,包括如下步骤:
(1)将多波束回声仪倒置式安装于水下潜标或者固定于水底,通过发射换能器阵向水面发射波束,通过接收换能器阵接收回波信号;
(2)多波束回声仪根据接收到的回波信号,获得回波的到达角和到达时间;
(3)根据声速剖面先验信息,求解经验正交函数EOF,得到声速剖面的降维基函数描述方法:
某一时刻t,声速剖面可以表示为:
式(1)中,z表示深度;c0(z)表示深度z处的平均声速,Δc(z,t)表示扰动量;i表示EOF的阶数,1≤i≤I,I表示EOF的总阶数,αi(t)表示第i阶EOF系数,ψi(z)表示第i阶EOF;
(4)声线跟踪算法:
由于海水的不均一性以及各向异性,声波在穿越不同声速层时会发生折射,从而导致声线弯曲;在非均匀介质条件下,假设声速只在垂直方向上变化,而在水平方向上是均匀的,声线传播遵循Snell定理:
式(2)中,c0表示换能器基阵位置处的声速,称之为表面声速;θ表示波束回波到达接收阵时与垂直方向的夹角;c(z)表示深度z处的声速,θ(z)表示深度z处声线切线与垂直方向的夹角;假设在很小的一段时间dt内,声线传播的水平距离dr、垂直距离dz以及路径长度ds满足
则自换能器基阵发出的声线到达深度z处的传播时间t为:
且声线传播的水平距离r为:
式(5)中,z0为换能器基阵深度;
由于海洋环境复杂多变,难以采用简单的连续函数表示声速剖面c(z);实测的声速剖面一般是由一些离散的测量点所组成,可以根据这些测量点的深度将整个水体分成若干层;具体为假设层内声速为常梯度变化进行声线跟踪计算;
首先,将水体沿深度方向分为N层,假设分层处的深度与声速值分别为z=[z0,z1,…,zN]T和c=[c0,c1,…,cN]T,并假设层内声速常梯度,即:
c(z)=ci+gi(z-zi), (6)
其中,i表示层数(i=0,2,…,N-1),ci表示第i层分界处的声速值,zi表示第i层的深度,gi表示第i层内的声速梯度;根据式(4)可知,在层内常梯度假设下,第i层内声线的双程传播时间可表示为:
根据式(5)可知,第i层声线传播的水平距离可表示为:
则整个深度范围内,声线总的传播时间、水平距离与深度的关系可表示为:
(5)结合EOF、声线跟踪算法、表面声速以及多波束数据,建立优化模型:
式(11)中,K为权重系数,Toameas(t)为t时刻所有波束到达时间的实际测量数据,Toaest(t)为t时刻所有波束到达时间的估计值,c0meas(t)为t时刻多波束回声仪接收换能器阵所在位置声速的测量数据,c0est(t)为t时刻多波束回声仪接收换能器阵所在位置声速的估计值;为测量算子,可根据声线跟踪算法进行计算,Doameas(t)为t时刻所有波束的到达方向测量值;
(6)根据步骤(5)建立的优化模型、步骤(2)中接收到的回波的到达角和到达时间以及接收换能器阵所在位置声速的测量值,利用优化算法获得EOF系数的反演值,进而利用式(1)获得声速剖面的估计结果。
2.根据权利要求1所述的基于倒置式多波束回声仪的声速剖面反演方法,其特征是,所述步骤(5)中,所述优化算法为模拟退火算法、爬山算法、粒子群算法、遗传算法、蚁群算法。
3.一种基于倒置式多波束回声仪的水温剖面反演方法,其特征是,包括如下步骤:
(1)将多波束回声仪倒置式安装于水下潜标或者固定于水底,通过发射换能器阵向水面发射波束,通过接收换能器阵接收回波信号;
(2)多波束回声仪根据接收到的回波信号,获得回波的到达角和到达时间;
(3)根据声速剖面先验信息,求解经验正交函数EOF,得到声速剖面的降维基函数描述方法:
某一时刻t,声速剖面可以表示为:
式(12)中,z表示深度;c0(z)表示深度z处的平均声速,Δc(z,t)表示扰动量;i表示EOF的阶数,1≤i≤I,I表示EOF的总阶数,αi(t)表示第i阶EOF系数,ψi(z)表示第i阶EOF;
(4)声线跟踪算法:
由于海水的不均一性以及各向异性,声波在穿越不同声速层时会发生折射,从而导致声线弯曲;在非均匀介质条件下,假设声速只在垂直方向上变化,而在水平方向上是均匀的,声线传播遵循Snell定理:
式(13)中,c0表示换能器基阵位置处的声速,称之为表面声速;θ表示波束回波到达接收阵时与垂直方向的夹角;c(z)表示深度z处的声速,θ(z)表示深度z处声线切线与垂直方向的夹角;假设在很小的一段时间dt内,声线传播的水平距离dr、垂直距离dz以及路径长度ds满足
则自换能器基阵发出的声线到达深度z处的传播时间为t:
且声线传播的水平距离r为:
式(16)中,z0为换能器基阵深度;
由于海洋环境复杂多变,难以采用简单的连续函数表示声速剖面c(z);实测的声速剖面一般是由一些离散的测量点所组成,可以根据这些测量点的深度将整个水体分成若干层;具体为假设层内声速为常梯度变化进行声线跟踪计算;
首先,将水体沿深度方向分为N层,假设分层处的深度与声速值分别为z=[z0,z1,…,zN]T和c=[c0,c1,…,cN]T,并假设层内声速常梯度,即:
c(z)=ci+gi(z-zi), (17)
式(17)中,i表示层数(i=0,2,…,N-1),ci表示第i层分界处的声速值,zi表示第i层的深度,gi表示第i层内的声速梯度;根据式(15)可知,在层内常梯度假设下,第i层内声线的双程传播时间可表示为:
根据式(16)可知,第i层声线传播的水平距离可表示为:
则整个深度范围内,声线总的传播时间、水平距离与深度的关系可表示为:
(5)结合EOF、声线跟踪算法、表面声速以及多波束数据,建立优化模型:
式(22)中,为EOF系数向量的估计结果,α为待估计的EOF系数向量,K为权重系数,取[0,1];Toameas(t)为t时刻所有波束到达时间的实际测量数据,Toaest(t)为t时刻所有波束到达时间的估计值,c0meas(t)为t时刻多波束回声仪接收换能器阵所在位置声速的测量数据,c0est(t)为t时刻多波束回声仪接收换能器阵所在位置声速的估计值;为测量算子,可根据声线跟踪算法进行计算,Doameas(t)为t时刻所有波束的到达方向测量值;
(6)根据步骤(5)建立的优化模型、步骤(2)中接收到的回波的到达角和到达时间以及接收换能器阵所在位置声速的测量值,利用优化算法获得EOF系数的反演值,进而利用式(3)获得声速剖面的估计结果;
(7)根据步骤(6)的声速剖面的估计值计算测量区域的水温剖面。
4.根据权利要求3所述的基于倒置式多波束回声仪的水温剖面反演方法,其特征是,所述步骤(5)中,所述优化算法为模拟退火算法、爬山算法、粒子群算法、遗传算法、蚁群算法。
5.根据权利要求3或4所述的基于倒置式多波束回声仪的水温剖面反演方法,其特征是,所述步骤(6)中,采用Mackenzie声速公式、Chen-Millero-Li声速公式、Dell Grosso声速公式、Leroy声速公式、Wilson精确声速公式、Wilson简化声速公式或EM分层简化声速公式计算测量区域的水温剖面。
6.一种倒置式多波束回声仪测量***,其特征包括:采用权利要求1或权利要求2所述的声速剖面反演方法。
7.一种倒置式多波束回声仪测量***,其特征包括:采用权利要求3、权利要求4或权利要求5所述的水温剖面反演方法。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910410164.4A CN110146895B (zh) | 2019-05-16 | 2019-05-16 | 基于倒置式多波束回声仪的声速剖面反演方法 |
PCT/CN2020/088233 WO2020228547A1 (zh) | 2019-05-16 | 2020-04-30 | 基于倒置式多波束回声仪的声速剖面反演方法 |
US17/232,163 US11733381B2 (en) | 2019-05-16 | 2021-04-16 | Sound velocity profile inversion method based on inverted multi-beam echo sounder |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910410164.4A CN110146895B (zh) | 2019-05-16 | 2019-05-16 | 基于倒置式多波束回声仪的声速剖面反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110146895A CN110146895A (zh) | 2019-08-20 |
CN110146895B true CN110146895B (zh) | 2021-04-20 |
Family
ID=67595603
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910410164.4A Active CN110146895B (zh) | 2019-05-16 | 2019-05-16 | 基于倒置式多波束回声仪的声速剖面反演方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US11733381B2 (zh) |
CN (1) | CN110146895B (zh) |
WO (1) | WO2020228547A1 (zh) |
Families Citing this family (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019010022A1 (en) * | 2017-07-03 | 2019-01-10 | R2Sonic, Llc | SYSTEM AND METHOD FOR MULTI-PERSPECTIVE INSONIFICATION |
CN109696235B (zh) * | 2018-12-27 | 2020-11-06 | 国家深海基地管理中心 | 一种深海有效声速确定方法 |
CN110146895B (zh) | 2019-05-16 | 2021-04-20 | 浙江大学 | 基于倒置式多波束回声仪的声速剖面反演方法 |
CN110837791B (zh) * | 2019-11-02 | 2023-04-07 | 山东科技大学 | 一种基于过完备字典的声速剖面反演方法 |
CN110865333B (zh) * | 2019-11-19 | 2021-08-24 | 浙江大学 | 洋流影响下水下滑翔机单信标无源声学定位方法 |
CN111142113A (zh) * | 2020-01-08 | 2020-05-12 | 青岛深海云海洋科技有限公司 | 一种实时化深海逆式回声观测*** |
CN112348280B (zh) * | 2020-11-23 | 2022-08-12 | 山东科技大学 | 一种利用改进型Argo浮标数据进行EOF分析的方法 |
CN112598113A (zh) * | 2020-12-15 | 2021-04-02 | 广东海洋大学 | 一种基于自组织竞争型神经网络的海洋声速剖面获取方法 |
CN113218493B (zh) * | 2021-04-08 | 2023-12-08 | 中国人民解放军国防科技大学 | 一种基于经验正交函数法的声速剖面反演方法 |
CN113281757B (zh) * | 2021-04-22 | 2023-07-11 | 浙江大学 | 一种倒置式多波束回声仪及海面波浪测量方法 |
CN113406645B (zh) * | 2021-05-10 | 2022-08-05 | 山东科技大学 | 一种平均声速水下声呐定位方法 |
CN113984178B (zh) * | 2021-09-18 | 2024-03-15 | 中国人民解放军国防科技大学 | 一种构建中尺度涡水下声速扰动普遍结构的方法 |
CN114061730B (zh) * | 2022-01-19 | 2023-09-19 | 中国船舶工业***工程研究院 | 目标散射回波变步长快速自适应估计方法 |
CN114646304B (zh) * | 2022-03-11 | 2022-11-08 | 广州海洋地质调查局 | 一种基于多波束数据的海洋内波识别方法 |
CN114782745B (zh) * | 2022-04-19 | 2024-04-23 | 广东海洋大学 | 一种基于机器学习的海洋声速剖面分类方法及装置 |
CN114978287B (zh) * | 2022-05-23 | 2023-03-28 | 北京理工大学 | 基于多波束协作的天基安全调制方法 |
CN115164944B (zh) * | 2022-09-09 | 2022-12-02 | 长江水利委员会水文局 | 一种多波束测深异常的角度和传播时间相关改正方法 |
CN115963445B (zh) * | 2023-01-04 | 2023-06-06 | 中国科学院声学研究所 | 水下传感节点的声学定位方法、装置、设备和介质 |
CN116044385B (zh) * | 2023-01-09 | 2024-04-26 | 西南石油大学 | 一种同位素示踪流量测井吸水剖面解释方法 |
CN116643061A (zh) * | 2023-04-13 | 2023-08-25 | 广东海洋大学 | 一种基于声学多普勒流速剖面仪的近岸海水层次划分方法 |
CN116202611B (zh) * | 2023-05-06 | 2023-07-14 | 中国海洋大学 | 一种基于元学习的小样本声速剖面反演方法 |
CN116774197B (zh) * | 2023-08-24 | 2023-11-03 | 交通运输部天津水运工程科学研究所 | 基于声呐扫测深度变化的固定目标物的声速剖面求取方法 |
CN117151198B (zh) * | 2023-09-06 | 2024-04-09 | 广东海洋大学 | 基于自组织竞争型神经网络的水声被动定位方法及装置 |
CN117665782B (zh) * | 2024-02-02 | 2024-04-26 | 山东科技大学 | 一种多探头多波束的实时声线弯曲改正方法 |
CN118226427A (zh) * | 2024-05-23 | 2024-06-21 | 山东科技大学 | 面向全海深环境的海量声呐波束高效声线跟踪方法及*** |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5608689A (en) * | 1995-06-02 | 1997-03-04 | Seabeam Instruments Inc. | Sound velocity profile signal processing system and method for use in sonar systems |
DK2435806T3 (en) * | 2009-05-27 | 2016-08-29 | Teledyne Instruments Inc | A method and system for remotely measuring the speed of sound |
JP5447323B2 (ja) * | 2010-10-06 | 2014-03-19 | 新日鐵住金株式会社 | 音速分布測定装置、音速分布測定方法及びプログラム |
CN102749622B (zh) * | 2012-07-03 | 2013-10-09 | 杭州边界电子技术有限公司 | 基于多波束测深的声速剖面及海底地形的联合反演方法 |
US9702973B2 (en) * | 2012-11-27 | 2017-07-11 | Nec Corporation | Environment measurement system and environment measurement method |
CN104181523B (zh) * | 2013-05-21 | 2017-12-29 | 中国科学院声学研究所 | 一种基于横摇稳定策略的多波束测深方法及*** |
CN105403860B (zh) * | 2014-08-19 | 2017-10-31 | 中国科学院声学研究所 | 一种基于支配相关的多稀疏声源定位方法 |
CN105487046A (zh) * | 2015-11-19 | 2016-04-13 | 中国石油大学(华东) | 大入射角声线跟踪定位方法 |
CN105911551B (zh) * | 2016-05-09 | 2018-05-08 | 浙江大学 | 一种基于加权集合卡尔曼滤波算法的声速剖面反演方法 |
CN106950568B (zh) * | 2017-02-23 | 2018-08-03 | 山东科技大学 | 一种自适应多节点等效声速剖面的构建方法 |
CN110146895B (zh) * | 2019-05-16 | 2021-04-20 | 浙江大学 | 基于倒置式多波束回声仪的声速剖面反演方法 |
-
2019
- 2019-05-16 CN CN201910410164.4A patent/CN110146895B/zh active Active
-
2020
- 2020-04-30 WO PCT/CN2020/088233 patent/WO2020228547A1/zh active Application Filing
-
2021
- 2021-04-16 US US17/232,163 patent/US11733381B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US11733381B2 (en) | 2023-08-22 |
CN110146895A (zh) | 2019-08-20 |
US20210231800A1 (en) | 2021-07-29 |
WO2020228547A1 (zh) | 2020-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110146895B (zh) | 基于倒置式多波束回声仪的声速剖面反演方法 | |
CN102749622B (zh) | 基于多波束测深的声速剖面及海底地形的联合反演方法 | |
CN106154276B (zh) | 基于海底混响和传播损失的深海海底参数反演方法 | |
CN110208812A (zh) | 半潜无人航行器海底三维地形探测装置及方法 | |
US7768874B2 (en) | System and method of horizontal wave measurement | |
CN109900256B (zh) | 一种自适应海洋移动声层析***和方法 | |
CN107132520A (zh) | 一种基于水声超短基线定位***的声线修正方法及*** | |
CN105258684A (zh) | 基于激光点云为约束的多波束低掠射角波束归位方法 | |
CN113156413B (zh) | 一种基于双程声路径的海底基准校准方法 | |
CN103217706B (zh) | 管理沿拖曳的声学线天线布置的声学节点网络的声学性能的方法和装置 | |
CN113281757B (zh) | 一种倒置式多波束回声仪及海面波浪测量方法 | |
CN108931777A (zh) | 海底冷泉输出气体运移通量的测量方法 | |
CN112540348A (zh) | 一种基于空间尺度上的声线修正算法在长基线水声定位***上的应用 | |
CN111220146B (zh) | 一种基于高斯过程回归学习的水下地形匹配定位方法 | |
CN112147578B (zh) | 一种高精度深水发射阵及多元垂直接收阵阵元定位***与方法 | |
CN114838808B (zh) | 一种基于浅海声信道响应的打桩冲击噪声源级测量***及方法 | |
US7161871B2 (en) | Method for determining ocean current and associated device | |
Didier et al. | Real-time correction of sound refraction errors in bathymetric measurements using multiswath multibeam echosounder | |
CN113009417B (zh) | 利用声场干涉特性的海底声学阵列阵形估计方法 | |
RU2769550C1 (ru) | Способ достроения измеренной от поверхности моря части профиля вертикального распределения скорости звука до дна | |
CN215340310U (zh) | 一种水下垂直阵阵元位置测量*** | |
CN111982343B (zh) | 一种适用于河流入海区域的水温声层析方法及其装置 | |
CN114675331A (zh) | 水面走航式探测海底气泡型浅层气的装置及方法 | |
RU2477498C1 (ru) | Метод мониторинга вертикального распределения скорости звука в условиях мелководных акваторий | |
Morozov et al. | Methodical aspects of the application of acoustic doppler current profilers in the black sea |
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 |