CN101699236B - 深海热液口温度场原位在线声学检测方法 - Google Patents
深海热液口温度场原位在线声学检测方法 Download PDFInfo
- Publication number
- CN101699236B CN101699236B CN 200910154610 CN200910154610A CN101699236B CN 101699236 B CN101699236 B CN 101699236B CN 200910154610 CN200910154610 CN 200910154610 CN 200910154610 A CN200910154610 A CN 200910154610A CN 101699236 B CN101699236 B CN 101699236B
- Authority
- CN
- China
- Prior art keywords
- underwater acoustic
- delta
- sound wave
- acoustic transducer
- temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Measuring Temperature Or Quantity Of Heat (AREA)
Abstract
本发明涉及一种深海热液口温度场原位在线声学检测方法。现有的深海热液口测温的主要手段无法对热液流体的温度场准确测量。本发明方法首先将使同层的其中一个水声换能器作为声源发声,其余十五个水声换能器作为接受端接受声波信号,由此计算出每条路径上的声波传播时间;然后计算声波传播路径的平均温度值t;最后重建海底热液口温度场,具体是将所在平面划分为多个网格,利用数据层析还原算法计算每个网格内的平均温度,从而得到温度场。本发明采用一种非接触式的声学热液温度场原位测量,可以克服传统接触式测量的缺点,在高温、腐蚀、多悬浮颗粒的恶劣环境下连续实时测量,提高热液口热通量测量的准确度。
Description
技术领域
本发明属于深海装备技术领域,涉及一种深海热液口温度场原位在线声学检测方法。
背景技术
对现代热液活动的温度场分布原位监测是热液活动研究中的重要内容之一,也是研究深海热液喷口及其周围深部生物圈、探究海底热液成矿作用机理的基础研究工作。温度场分布又是热液系统输入海洋热通量计算的重要参数。但由于热液活动区的高温高压极端环境限制,现有的测定热液口温度的有效手段和技术非常有限。迄今为止已经开发了可以测量热液口点温度的温度计。现有的手段主要用温度传感器阵列进行测量,但传统的接触式温度测量装置具有很大的局限性:首先,接触式测量会干扰热液口的温度;其次,接触式温度测量只能测量热液口某些点的温度,无法得到热液区域二维和三维的温度场分布数据。综上所述,目前深海热液口测温的主要手段无法对热液流体的温度场准确测量。因此海洋科学家无法对脊-洋反应区内整个温洋脊剖面的热能通量进行准确估计,只能在现有资料基础上进行一个恰当的估计。目前,国内尚未有一种能够实现热液喷口温度场原位监测,给出定量温度信息的方法面市。
发明内容
本发明的目的在于克服传统的海底热液喷口温度测量的缺点,提供了一种运行稳定可靠、测量准确快速的深海热液口温度场原位在线声学检测方法。
本发明方法的步骤是:
步骤(1)、将三十二个水声换能器分为两层设置在支架上,层高为0.5m,每一层有十六个水声换能器,十六个水声换能器组成围合成正方形,正方形的每边上设置四个水声换能器,每条边上水声换能器的间距为0.3m;
步骤(2)、用扫频范围为18kHz~23kHz,声压级为69dB的扫频信号源发射卡驱动第一层任意一个水声换能器发射声波信号,同层其余15个水声换能器都作为接收水声换能器接收声信号,分别将这15个接收水声换能器接收到的声波信号与发射水声换能器发射的声波信号进行互相关函数分析
其中R(n)为互相关函数,N为采样点数,可根据所需的测量精度设定,f1(k)是发送水声换能器发送的声源信号,f2(k)是接收水声换能器采集到的离散声波信号,互相关函数取得最大值的延迟时间即为声波的传播时间。
再通过切换电路把按照逆时针方向确定的紧临上一个水声换能器的下一个水声换能器作为发射声波信号的水声换能器,同层其余十五个水声换能器作为接收声波信号的水声换能器,再次利用互相关函数计算声波的传播时间,以此类推,直到完成第二层中的最后一个水声换能器作为发射声波信号的水声换能器时的互相关函数分析。这样切换电路共切换三十一次形成一个完整的周期;
步骤(3)、计算声波传播路径的平均温度值t。
c(D,S,t)=c(0,S,t)+(16.23+0.253t)D
+(0.213-0.1t)D2+[0.016
+0.0002(S-35)](S-35)tD
式中:
c(0,S,t)=1449.05+45.7t-5.21t2
+0.23t3-(1.333-0.126t
+0.009t2)(S-35)
其中c(D,S,t)为声波在海水介质中的速度,D为被测平面的深度,单位为Km;t=T/10,T为海水温度,单位为℃,S为含盐量的千分数。在固定的热液口盐度和被测海平面的深度均可视为常数,因此声速变化主要与温度有关;
步骤(4)、利用数据层析还原算法重建海底热液口温度场。
建立同层平面中十六个水声换能器之间的九十六条声波传播路径的声学测量模型。将被测正方形平面划分为8×8共六十四个网格区域。假定每块区域上的温度分布相同,网格区域标记为m,Δs为声波通过每一块区域的长度,则可建立如下关系。
其中,m=64,n=96,Tn是第n条声波传播路径上的声波飞渡时间。ai是第i个网格区域上声速的倒数,Δsni是声波在第n条传播路径通过第i个网格区域的长度,ai与Δsni都是常数。ai是各路径上温度的特征函数,是要计算的值,写成矩阵形式,得到如下的方程组:
AX=t
其中
矩阵A是已知的,矩阵t中的元素为步骤(2)中得到的声波传播时间,对于此线性方程组,由总体最小二乘方法可得到变量X,进而求出声波在每一块网格区域上传播速度的倒数,利用步骤(3)中速度与温度的关系式,可重建出海底热液口截面温度场。
所述的总体最小二乘方法为:
(A+ΔA)X=t+Δt
其中,ΔA表示声波在网格中传播路径长度的噪声矩阵,Δt表示声波在网格中传播路径上飞渡时间测量数据的噪声向量。则(A+ΔA)X=t+Δt可以表述为:
或者
(W+ΔW)M=0
其中,W=(At),ΔW=(ΔA Δt),M=(X-1)T,W为增广输入数据矩阵,ΔW是叠加于增广数据矩阵的噪声,M是增广权向量。这样,(A+ΔA)X=t+Δt总体最小二乘解可以简单表示为求解一个解向量MTLS,使得min||ΔA Δt||F或min||ΔW||F为最小。对增广输入数据矩阵W进行奇异值分解,得到:
其中,σj是第j个奇异值,∑=diag(σ1,…,σn+1),则总体最小二乘解为:
得出解向量MTLS后再计算变量X,求出声波在每一块网格区域上传播速度的倒数,利用步骤(3)中速度与温度的关系式,可精确地重建出海底热液口截面温度场。
本发明采用一种非接触式的声学热液温度场原位测量,可以克服传统接触式测量的缺点,在高温、腐蚀、多悬浮颗粒的恶劣环境下连续实时测量,提高热液口热通量测量的准确度。为科学家研究洋壳散热量、建立大洋环流模型、研究热液对临近海域海洋物理特性影响提供帮助。
附图说明
图1为同层平面的声学路径和拓扑网络图。
具体实施方式
以下结合附图对本发明的原理、结构做进一步的说明。
本发明方法所采用的测温装置由水声换能器、机架、电缆和数据采集系统组成。安装在机架上的水声换能器通过电缆连接到驱动与前置放大腔、驱动与前置放大腔、信号采集腔和电池腔组成的数据采集系统。
如图1所示,是同层平面的声学路径和拓扑网络图。在平面的每一边放置四个水声换能器,四个位置处于每条边五等分处。这样一个平面中共有十六个个水声换能器。
本发明方法的具体步骤如下:
步骤(1)、在支架上共安装两层RHS-30水声换能器,每层RHS-30水声换能器都处于同一水平面,每层装有十六个水声换能器,层高0.5m,每层的支架围合成正方形,正方形的每边上设置四个标准水声换能器,每条边上水声换能器的间距为0.3m,按照逆时针的顺序从T1换能器到T16换能器均匀布置。
步骤(2).用扫频范围为18kHz~23kHz,声压级为69dB的扫频信号源发射卡驱动第T1个水声换能器发射声波信号,同层其余15个水声换能器(T2至T16)都作为接收水声换能器接收声信号,分别将T2至T16接收水声换能器接收到的声波信号与发射水声换能器T1发射的声波信号利用公式(1)进行互相关函数分析
其中R(n)为互相关函数,N为采样点数,可根据所需的测量精度设定,f1(k)是发送水声换能器发送的声源信号,f2(k)是接收水声换能器采集到的离散声波信号,互相关函数取得最大值的延迟时间即为声波的传播时间(TOF)。
再通过切换电路把第2个水声换能器T2作为发射声波信号的水声换能器,其余同层15个作为接收声波信号的水声换能器,再次利用公式(1)计算声波的传播时间,以此类推,直到完成第二层中的第16个水声换能器作为发射声波信号的水声换能器时的互相关函数分析。这样切换电路共切换31次才是一个完整的周期。每个测量周期的时间约为1分钟。
步骤(3)、计算声波传播路径的平均温度值t。已知声波在海水中的传播速度与温度、盐度和深度有关,关系式如下:
c(D,S,t)=c(0,S,t)+(16.23+0.253t)D
+(0.213-0.1t)D2+[0.016 (2)
+0.0002(S-35)](S-35)tD
其中:
c(0,S,t)=1449.05+45.7t-5.21t2
+0.23t3-(1.333-0.126t
+0.009t2)(S-35)
c(D,S,t)是声波在某种海水介质中的速度,单位为m/s;D为被测平面的深度,单位为Km;t=T/10,液口盐度和被测海平面的深度均可视为常数,声速变化主要与温度有关。利用式(2)中温度与声速的关系得出这些有效传播路径的平均温度值。
步骤(4)、利用数据层析还原算法重建海底热液口温度场。建立如图
(2)所示的同层16个水声换能器,96条声波传播路径的声学测量模型。水声换能器对称分布在正方形控制面上,被测空间划分为8×8共64个网格区域。假定每块区域上的温度分布相同,网格区域标记为m,Δs为声波通过每一块区域的长度,则可建立如下关系。
其中,m=64,n=96,Tn是第n条声波传播路径上的声波飞渡时间,ai是第i个网格区域上声速的倒数,Δsni是声波在第n条传播路径通过第i个网格区域的长度,它们都是常数,取决于测温系统装置和各水声换能器的安装位置。ai是各路径上温度的特征函数,是要计算的值,写成矩阵形式,得到如下的方程组:
AX=t (4)
其中,
方程(4)中,矩阵A是已知的,矩阵t中的元素为步骤(2)中利用公式(1)得到的声波传播时间,对于此线性方程组,由方法一即最小二乘法求解可得到变量X,求出声波在每一块网格区域上传播速度的倒数,利用式(2)中速度与温度的关系式,可重建出海底热液口截面温度场。
最小二乘作为一种数据拟合方法,在求解矩阵方程AX=t时,仅假设方程右端n条路径飞渡时间矩阵t中的测量数据存在误差。而在实际测量系统中,声波在网格中的传播路径长度矩阵A和飞渡时间矩阵t中同时存在误差,因此最小二乘方法具有局限性。更合理的方法是同时考虑A、t二者的扰动。本实施例中优先采用Golub G H,Loan C F V等人提出的总体最小二乘方法求解矩阵反问题,使得来自A、t的噪声扰动影响最小。总体最小二乘的求解方法如下,考虑矩阵方程:
(A+ΔA)X=t+Δt (5)
其中,ΔA表示声波在网格中传播路径长度的噪声矩阵,Δt表示声波在网格中传播路径上飞渡时间测量数据的噪声向量。式(5)可以表述为:
或者
(W+ΔW)M=0 (7)
其中,W=(A t),ΔW=(ΔA Δt),M=(X-1)T,W为增广输入数据矩阵,ΔW是叠加于增广数据矩阵的噪声,M是增广权向量。这样,矩阵方程(5)的总体最小二乘解可以简单表示为求解一个解向量MTLS,使得min||ΔA Δt||F或min||ΔW||F为最小。对增广输入数据矩阵W进行奇异值分解,得到:
得出解向量MTLS后再计算变量X,求出声波在每一块网格区域上传播速度的倒数,利用式(2)中速度与温度的关系式,可精确地重建出海底热液口截面温度场。
Claims (1)
1.深海热液口温度场原位在线声学检测方法,其特征在于该方法包括如下步骤:
步骤(1)、将三十二个水声换能器分为两层设置在支架上,层高为0.5m,每一层有十六个水声换能器,十六个水声换能器组成围合成正方形,正方形的每边上设置四个水声换能器,每条边上水声换能器的间距为0.3m;
步骤(2)、用扫频范围为18kHz~23kHz,声压级为69dB的扫频信号源发射卡驱动第一层任意一个水声换能器发射声波信号,同层其余15个水声换能器都作为接收水声换能器接收声信号,分别将这15个接收水声换能器接收到的声波信号与发射水声换能器发射的声波信号进行互相关函数分析
其中R(n)为互相关函数,N为采样点数,可根据所需的测量精度设定,f1(k)是发射水声换能器发送的声源信号,f2(k)是接收水声换能器采集到的离散声波信号,互相关函数取得最大值的延迟时间即为声波的传播时间;
再通过切换电路把按照逆时针方向确定的紧临上一个水声换能器的下一个水声换能器作为发射声波信号的水声换能器,同层其余十五个水声换能器作为接收声波信号的水声换能器,再次利用互相关函数计算声波的传播时间,以此类推,直到完成第二层中的最后一个水声换能器作为发射声波信号的水声换能器时的互相关函数分析;这样切换电路共切换三十一次形成一个完整的周期;
步骤(3)、计算声波传播路径的平均温度值t,
c(D,S,t)=c(0,S,t)+(16.23+0.253t)D
+(0.213-0.1t)D2+[0.016
+0.0002(S-35)](S-35)tD
式中:
c(0,S,t)=1449.05+45.7t-5.21t2
+0.23t3-(1.333-0.126t
+0.009t2)(S-35)
其中c(D,S,t)为声波在海水介质中的速度,D为被测平面的深度,单位为Km;t=T/10,T为海水温度,单位为℃,S为含盐量的千分数;在固定的热液口盐度和被测海平面的深度均可视为常数,因此声速变化主要与温度有关;
步骤(4)、利用数据层析还原算法重建海底热液口温度场;
建立同层平面中十六个水声换能器之间的九十六条声波传播路径的声学测量模型,将被测正方形平面划分为8×8共六十四个网格区域,假定每块区域上的温度分布相同,网格区域标记为m,Δs为声波通过每一块区域的长度,则可建立如下关系:
其中,m=64,n取1至96,Tn是第n条声波传播路径上的声波飞渡时间,ai是第i个网格区域上声速的倒数,Δsni是声波在第n条传播路径通过第i个网格区域的长度,ai与Δsni都是常数,ai是各路径上温度的特征函数,是要计算的值,写成矩阵形式,得到如下的方程组:
AX=t
其中
矩阵A是已知的,矩阵t中的元素为步骤(2)中得到的声波传播时间,对于此线性方程组,由总体最小二乘方法可得到变量X,进而求出声波在每一块网格区域上传播速度的倒数,利用步骤(3)中速度与温度的关系式,可重建出海底热液口截面温度场;
所述的总体最小二乘方法为:
(A+ΔA)X=t+Δt
其中,ΔA表示声波在网格中传播路径长度的噪声矩阵,Δt表示声波在网格中传播路径上飞渡时间测量数据的噪声向量,则(A+ΔA)X=t+Δt可以表述为:
或者
(W+ΔW)M=0
其中,W=(At),ΔW=(ΔAΔt),M=(X-1)T,W为增广输入数据矩阵,ΔW是叠加于增广数据矩阵的噪声,M是增广权向量;这样,(A+ΔA)X=t+Δt总体最小二乘解可以简单表示为求解一个解向量MTLS,使得min||ΔAΔt||F或min||ΔW||F为最小,对增广输入数据矩阵W进行奇异值分解,得到:
得出解向量MTLS后再计算变量X,求出声波在每一块网格区域上传播速度的倒数,利用步骤(3)中速度与温度的关系式,可精确地重建出海底热液口截面温度场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910154610 CN101699236B (zh) | 2009-11-19 | 2009-11-19 | 深海热液口温度场原位在线声学检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910154610 CN101699236B (zh) | 2009-11-19 | 2009-11-19 | 深海热液口温度场原位在线声学检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101699236A CN101699236A (zh) | 2010-04-28 |
CN101699236B true CN101699236B (zh) | 2013-04-24 |
Family
ID=42147699
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200910154610 Expired - Fee Related CN101699236B (zh) | 2009-11-19 | 2009-11-19 | 深海热液口温度场原位在线声学检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101699236B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2853818C (en) * | 2011-10-27 | 2020-02-25 | Paul LARKIN | Drive signal distribution for direct field acoustic testing |
CN104155030B (zh) * | 2014-08-13 | 2016-10-05 | 沈阳工业大学 | 一种考虑声线弯曲的声学ct温度场重建方法 |
CN104464855A (zh) * | 2014-12-25 | 2015-03-25 | 中国科学院合肥物质科学研究院 | 一种基于高频电磁力的液态重金属超声波测温装置 |
CN104502001B (zh) * | 2014-12-31 | 2017-10-13 | 中国科学院南海海洋研究所 | 一种海洋湍流热通量仪及其使用方法 |
CN107389220B (zh) * | 2017-07-06 | 2019-07-16 | 中国船舶重工集团公司第七一九研究所 | 一种锅炉内温度场分布的检测方法 |
CN110006550B (zh) * | 2019-02-02 | 2024-02-23 | 四川大学 | 一种声波传播路径自适应组网的温度场重建系统及其方法 |
CN112197884B (zh) * | 2020-10-29 | 2024-04-16 | 西安科技大学 | 基于声学法测量松散介质温度的实验装置及方法 |
CN113375828B (zh) * | 2021-05-10 | 2023-03-28 | 电子科技大学 | 一种基于超声波的温度场重建方法 |
CN113465770B (zh) * | 2021-05-17 | 2022-05-10 | 浙江大学 | 一种小尺度垂直剖面二维温度场观测方法 |
CN114216582B (zh) * | 2021-12-08 | 2022-09-06 | 电子科技大学长三角研究院(湖州) | 一种高精度快速温度场重建方法、系统、设备及终端 |
CN115248030B (zh) * | 2021-12-28 | 2023-09-22 | 浙江海洋大学 | 一种基于声层析的海洋涡旋涡量测量方法 |
CN115389046B (zh) * | 2022-10-26 | 2023-08-11 | 浙江大学 | 测温方法及装置 |
-
2009
- 2009-11-19 CN CN 200910154610 patent/CN101699236B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101699236A (zh) | 2010-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101699236B (zh) | 深海热液口温度场原位在线声学检测方法 | |
US11733381B2 (en) | Sound velocity profile inversion method based on inverted multi-beam echo sounder | |
CN105911551B (zh) | 一种基于加权集合卡尔曼滤波算法的声速剖面反演方法 | |
Lermusiaux et al. | Coupled ocean–acoustic prediction of transmission loss in a continental shelfbreak region: Predictive skill, uncertainty quantification, and dynamical sensitivities | |
Fan et al. | Experimental study on underwater acoustic imaging of 2-D temperature distribution around hot springs on floor of Lake Qiezishan, China | |
CN104678384B (zh) | 一种波束域的声压差互相关谱分析水下目标速度估计方法 | |
CN105403878A (zh) | 一种基于时延敏感核的海洋声层析方法 | |
CN103076594A (zh) | 一种基于互相关的水声脉冲信号双阵元定位的方法 | |
CN103217706B (zh) | 管理沿拖曳的声学线天线布置的声学节点网络的声学性能的方法和装置 | |
Zhu et al. | Method for rapid impact localization for subsea structures | |
Lin et al. | Merging multiple-partial-depth data time series using objective empirical orthogonal function fitting | |
Pritchard | Arctic Ocean background noise caused by ridging of sea ice | |
CN116660996B (zh) | 一种基于深度学习的漂流式浅海局部地声参数预测方法 | |
Chen et al. | Study on acoustic reflection characteristics of layered sea ice based on boundary condition method | |
CN101726626B (zh) | 深海热液口速度场原位在线声学检测方法 | |
CN105785350A (zh) | 一种基于最大熵原理的浅海流场声层析方法 | |
Miron-Morin et al. | The oceanographic sensitivity of the acoustic channel in shallow water | |
Kexin et al. | Physics-aided data-driven modal ocean acoustic propagation modeling | |
Skarsoulis et al. | Tomographic inversions from the Fram Strait 2008–9 experiment | |
Goodney et al. | Water temperature sensing with microtomography | |
Lunkov | Interference structure of low-frequency reverberation signals in shallow water | |
McCarthy et al. | Reduced-Order Machine-Learning Model for Transmission Loss Prediction in Underwater Acoustics | |
Ohta et al. | Inversion for seabed geoacoustic properties in shallow water experiments | |
Zhou et al. | Impact location estimation in anisotropic structures | |
Attarsharghi | Design criteria and practical insights into an underwater current measurement system along with simulation results of a real-case scenario in the Northwest Atlantic Ocean |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130424 Termination date: 20141119 |
|
EXPY | Termination of patent right or utility model |