专利名称：声多普勒海流剖面仪测量海流剖面流速的方法和系统的制作方法
技术领域：
本发明涉及海流剖面流速的测量，更具体地说，涉及一种声多普勒海流剖面仪测量流速的方法和宽带声多普勒海流剖面仪。
背景技术：
宽带声多普勒海流剖面仪是测量海流剖面流速的声学设备。目前，宽带声多普勒海流剖面仪的信号处理方法是相关脉冲对协方差法，它发射的脉冲对由两个子脉冲组成，每个子脉冲包括多个脉冲元，求两个子脉冲中对应脉冲元之间的相关函数的相位，由此求出流速。(参见文献[1]Zhu Weiqing，etc，″Spectral moment estimatesof broadband backscattering acoustic wave in moving medium″，Chinese Journal ofAcoustics，Vol.15，No.2，1996，106-114)它的速度估计V的标准偏差由下式表示σ(V)=Va[116π5/2Lρ2(1-ρ2+2NS+N2S2)]1/2---(1)]]>其中，λ为剖面仪工作的中心频率对应的波长，T和ρ分别为相关脉冲对中两个子脉冲的时间间距和相关系数，L为子脉冲的脉冲元的个数，S/N为信噪比。特别是Va＝λ/2T，称为模糊速度，当载体与海流间的相对速度超过Va后，就会发生相位�：�。这里所说的相位�：侵赶辔豢梢允铅铡�2mπ中的任一个值，其中m＝0，±1……±N。因此在发生相位�：保辔徊獠蛔迹俣纫膊獠蛔�。
在B.H.Brumley等人的美国专利5208785、5483499和5615173“Broadbandacoustic Doppler current profiler”中介绍了锚系安装的宽带声多普勒海流剖面仪，海底有一沉块，用系缆与宽带声多普勒海流剖面仪相连，剖面仪的四个波束向上，测量剖面内的流速。剖面仪测量的是载体与海流的相对速度，是载体速度与流速的矢量和。通过锚系安装方式将载体速度限制得相当低，测量的主要是流速，而流速值比较�。梢员苊夥⑸辔荒：菀谆竦米既返牧魉僦�。
但是，锚系安装的宽带声多普勒海流剖面仪存在不足。在实际工作中，常常希望将宽带声多普勒海流剖面仪安装在船上，在走航的过程中测量大面积海域的流�。耸逼拭嬉遣饬康氖谴俣扔肓魉俚氖噶亢停捎诖俦攘魉俅螅３Ｊ勾怂俣仁噶亢痛笥谀：俣萔a，产生相位模糊，速度测不准。
因此就需要有一种宽带声多普勒海流剖面仪测量流速的方法，在该方法中，声多普勒海流剖面仪可以安装在航行中的载体上。还需要有一种宽带声多普勒海流剖面仪测量流速的方法，该方法能克服相位�：佣既返夭獬隽魉�。
发明内容
本发明的目的是为安装在航行中载体上的声多普勒海流剖面仪提供一种测量海流剖面流体速度的方法，并对。
为了实现上述目的，本发明提供一种声多普勒海流剖面仪测量海流剖面流速的方法，声多普勒海流剖面仪装载于载体上，包括如下步骤(1)预估一船速和流速的矢量和V的估计值VE，该估计值VE大于真实的船速和流速的矢量和V。
(2)根据估计值VE设置相关脉冲对中两个子脉冲之间的第一时间间隔T1，其中，T1≤λ/(2VE)，λ为剖面仪的工作中心频率对应的波长，优选0.7λ/(2VE)≤T1≤λ/(2VE)。
(3)根据步骤(2)中的第一时间间隔T1选择相关脉冲对发射编码，声多普勒海流剖面仪向待测的流体介质发射相关脉冲对信号，并根据回波信号采用相关脉冲对协方差法计算得到船速和流速的矢量和V的粗测值VS。
(4)根据粗测值VS设置相关脉冲对中两个子脉冲之间的第二时间间隔T2，其中，T2≤λ/(2VS)，优选0.9λ/(2VS)≤T2≤λ/(2VS)。
(5)根据步骤(4)中的第二时间间隔T2选择相关脉冲对发射编码，声多普勒海流剖面仪向待测的流体介质发射相关脉冲对信号，并根据回波信号采用相关脉冲对协方差法计算得到船速和流速的矢量和V的精测值VP。
(6)在精测值VP中扣除船速，获得流速。
本发明还提供一种声多普勒海流剖面仪，包括计算机406，该计算机406内包括速度估计�？椋糜诖娲⒃ど璧幕蛲獠渴淙氲拇俸土魉俚氖噶亢偷腣估计值VE；第一设置�？椋糜诟菟俣裙兰颇？橹械墓兰浦礦E设置相关脉冲对中两个子脉冲之间的第一时间间隔T1，其中，T1≤λ/(2VE)，λ为剖面仪的工作中心频率对应的波长；第二设置模块，用于根据速度计算�？榛竦玫拇植庵礦S设置相关脉冲对中两个子脉冲之间的第二时间间隔T2，其中，T2≤λ/(2VS)，λ为剖面仪的工作中心频率对应的波长；速度计算模块，用于根据第一设置�？橹械牡谝皇奔浼涓鬞1或第二设置模块中的第二时间间隔T2选择相关脉冲对发射编码，并根据剖面仪采集的回波信号采用相关脉冲对协方差法计算船速和流速的矢量和V的粗测值VS或精测值VP。
本发明优点在于1.本发明的声多普勒海流剖面仪测量流速的方法适用于航行中的载体。
2.本发明的自动选择子脉冲时间间距的方法，在宽广的范围内适应载体的船速和流速的变化，同时达到标定的作用距离和测深精度。
3.本发明的声多普勒海流剖面仪测量流速的方法，同样也适用锚系的载体，尤其是流速大的水域。
图1是声多普勒海流剖面仪工作示意图；图2是声多普勒海流剖面仪装置示意图；图3是本发明的声多普勒海流剖面仪软件流程图。
图面说明载体100 湿端200 换能器阵300 电子分机400 传感器500 干端600波束101 海底102 流体10具体实施方式
下面结合附图与具体方式对本发明作进一步详细描述。
如图1所示，声多普勒海流剖面仪测量流体速度时，它装于载体(船)100上，声多普勒海流剖面仪总体上包括湿端200和干端600。所述的湿端包括换能器阵300、电子分机400和传感器500；所述的干端600包括终端计算机601和若干外设，它们装在水面上。声多普勒海流剖面仪通过换能器阵300向水下发射波束101，并接受流体103和/或海底102的反射信号，经电子分机400处理后得到流体103的速度。
声多普勒海流剖面仪的装置示意图如图2所示，它包括由换能器阵列组成的换能器阵300、电子分机400、传感器500和干端600。电子分机包括发射机401、接收机402、数据采集器403、高速数字信号处理器404、声纳控制器405、PC104计算机406以及PC104总线407。干端600包括终端计算机601、电罗经602、打印机603和显示器604，其中电罗经602为选配件。
电子分机400中发射机401、接收机402与换能器阵300连接，发射机401分别与声纳控制器405和高速数字信号处理器404连接，数据采集器403与接收机402连接，高速数字信号处理器404与数据采集器403连接，声纳控制器405与高速数字信号处理器404连接，PC104计算机406经由PC104总线407与数据采集器403、高速数字信号处理器404和声纳控制器405连接。传感器500包括温度传感器501和姿态传感器502，分别与声纳控制器405连接。
干端600中终端计算机601与电子分机400中的PC104计算机406连接，终端计算机601与电罗经602、打印机603和显示器604连接。
该系统的专用测速程序装载在计算机406的存储器中，在本发明中，该程序包括速度估计�？�、第一设置�？�、第二设置�？楹退俣燃扑隳？椋赐�3所示流程图中的步骤执行。
在图3中，701是开始步骤，由终端计算机601发出指令，传输到计算机406，启动计算机406存储器中的程序，使声纳处于工作状态。步骤702和703中，对软件和系统硬件进行初始化。
在步骤704中，根据海底的深度选择发射编码，这里的发射编码为相关脉冲对发射编码。
在根据步骤704的发射编码发射声信号之前，需要设置该相关脉冲对中两个子脉冲的时间间距T。已经知道，�：俣萔a＝λ/2T，其中，λ为剖面仪工作的中心频率对应的波长。由Va的表式可知，当λ固定，T要短些，以使Va变大，则剖面仪适用的速度范围大些。但是T短了，子脉冲的脉冲元的个数L和信噪比S/N就�。饔镁嗬氡涠蹋馑倬缺洳�。特别是水深时，T要长，使L和S/N变大，保持足够的信噪比，以达到标定的作用距离和测速精度。所以在选择发射编码时，相关脉冲对中两个子脉冲的时间间距T的合理选取是十分重要的。这也是本发明的特征所在。
在步骤705中，速度估计模块预设一个船速和流速的矢量和V的估计值VE，或者由该速度估计�？榻邮沼刹僮魅嗽笔淙氲囊桓龉兰浦礦E。该估计值VE可以根据测量船和待测量流体区域的现场情况来确定，为了保证不发生相位�：霉兰浦礦E的取值要比实际船速和流速的矢量和V要大，这一点是本领域的一般技术人员很容易判断并作出估计的。例如，当V大约为2m/s的时候，可以将估计值VE设为2.6m/s。
在步骤706中，第一设置�？楦莨兰浦礦E设置相关脉冲对中两个子脉冲之间的第一时间间隔T1，其中，T1≤λ/(2VE)。当T1≤λ/(2VE)时，可得到VE≤λ/(2T1)。在选定T1后，根据�：俣萔a1＝λ/2T1可知VE≤Va1，又因为估计值VE的取值要比实际船速和流速的矢量和V要大，这样就可以保证在速度测量时不会发生相位的�：５牵捎赥1的取值越�。饬烤仍讲�。在实际操作中优选0.7λ/(2VE)≤T1≤λ/(2VE)，(2)这样既保证不出现相位�：挚梢员Ｖひ欢ǖ牟饬烤取＠纾诓街�(1)的例子中，将VE设为2.6m/s后，在λ为0.5cm时，T1应该满足T1≤0.96ms，优选0.67ms≤T1≤ 0.96ms。
在步骤707中，根据步骤704选择的发射编码和步骤706得到的第一时间间隔T1，计算机406经PC104总线407使高速数字信号处理器404和声纳控制器405产生发射信号，然后驱动发射机401，再驱动换能器阵300，发出声脉冲对到流体介质中。在步骤708中，换能器阵300接收从流体介质和海底介质反向散射的回波信号，经接收机402馈送给数据采集器403采样，再馈送给高速数字信号处理器405。在步骤709中，对回波信号进行调解滤波，得到回波的基波信号。
在步骤710中，速度计算�？楦莶街�709中获得的回波信号，采用常规的相关脉冲对协方差法计算速度，即根据回波信号求延时为T1的相关函数R(T1)，并由R(T1)求得速度船速和流速的矢量和VS。在这里，将计算得到的VS称为粗测值VS。在步骤705和706中，T1是按照船速和流速的矢量和V的估计值VS来选取的，由于VE的设定值要比实际速度V大，使得T1的取值要比其最佳的取值要�。佣沟貌饬烤炔患�。因此，将该步骤中测得的VS称为粗测值，但是该粗测值VS要比估计值VE更接近于实际的船速和流速矢量和V。
在步骤711中，第二设置�？楦荽植庵礦S设置相关脉冲对中两个子脉冲之间的第二时间间隔T2，其中，T2≤λ/(2VS)。在选定T2后，根据模糊速度Va2＝λ/2T2可知VS≤Va2，这样，就可以保证在速度测量时不会发生相位的�：�。特别是，由于粗测值VS已经基本上能够反映船速和该流域流速的情况，因此，在选择T2时，则可以在不出现相位�：那榭鱿陆渚×垦≡竦谋冉洗螅佣岣卟饬烤�。在实际操作中优选0.9λ/(2VS)≤T2≤λ/(2VS)，(3)例如，在步骤710测得的粗测值VS为2.1m/s，λ为5cm时，T2应该满足T2≤1.19ms，优选1.07ms≤T2≤1.19ms。
最后，返回步骤707，根据步骤711确定的两个子脉冲之间的第二时间间隔T2重复进行步骤707～710，获得最后的船速和流速矢量和精测值VP。再根据外输入的准确船速，求出准确的流速Vc。
最后还可以将温度传感器501和姿态传感器502的数据经声纳主控器405馈送给PC104计算机，再馈送给终端计算机601，融合电罗经602的数据，给出校正过的流速数据。
声多普勒海流剖面仪装在航行的载体上测量时，在标定的作用距离和测深精度下，采用本发明的方法可避免发生相位�：�。
此外，在一般情况下，载体的速度是大量而流速是小量，所以在本发明的方法中求出的第二时间间隔T2主要决定于船速。另一方面，由于载体的航速(即船速)是一般是恒定的，因此T2在相当长的时间内是稳定的。这样，在相当长的时间内可按照第二时间间隔T2发射同一个信号，这为数据后处理提供了方便。
本发明的方法中采取的避免相位�：拇胧┮彩视糜诿档脑靥澹绕涫窃诹魉俳洗蟮乃�。
权利要求
1.一种声多普勒海流剖面仪测量海流剖面流速的方法，声多普勒海流剖面仪装载于载体上，包括如下步骤(1)预估一船速和流速的矢量和V的估计值VE，该估计值VE大于真实的船速和流速的矢量和V；(2)根据估计值VE设置相关脉冲对中两个子脉冲之间的第一时间间隔T1，其中，T1≤λ/(2VE)，λ为剖面仪的工作中心频率对应的波长；(3)根据步骤(2)中的第一时间间隔T1选择相关脉冲对发射编码，声多普勒海流剖面仪向待测的流体介质发射相关脉冲对信号，并根据回波信号采用相关脉冲对协方差法计算得到船速和流速的矢量和V的粗测值VS；(4)根据粗测值VS设置相关脉冲对中两个子脉冲之间的第二时间间隔T2，其中，T2≤λ/(2VS)；(5)根据步骤(4)中的第二时间间隔T2选择相关脉冲对发射编码，声多普勒海流剖面仪向待测的流体介质发射相关脉冲对信号，并根据回波信号采用相关脉冲对协方差法计算得到船速和流速的矢量和V的精测值VP；(6)在精测值VP中扣除船速，获得流速。
2.根据权利要求1所述的声多普勒海流剖面仪测量流速的方法，其特征在于，在步骤(2)中，0.7λ/(2VE)≤T1≤λ/(2VE)。
3.根据权利要求1或2所述的声多普勒海流剖面仪测量流速的方法，其特征在于，在步骤(4)中，0.9λ/(2VS)≤T2≤λ/(2VS)。
4.一种实施权利要求1所述方法的声多普勒海流剖面仪，包括计算机(406)；其特征在于，该计算机(406)内包括速度估计�？椋糜诖娲⒃ど璧幕蛲獠渴淙氲拇俸土魉俚氖噶亢偷腣估计值VE；第一设置模块，用于根据速度估计模块中的估计值VE设置相关脉冲对中两个子脉冲之间的第一时间间隔T1，其中，T1≤λ/(2VE)，λ为剖面仪的工作中心频率对应的波长；第二设置�？椋糜诟菟俣燃扑隳？榛竦玫拇植庵礦S设置相关脉冲对中两个子脉冲之间的第二时间间隔T2，其中，T2≤λ/(2VS)，λ为剖面仪的工作中心频率对应的波长；速度计算模块，用于根据第一设置�？橹械牡谝皇奔浼涓鬞1或第二设置模块中的第二时间间隔T2选择相关脉冲对发射编码，并根据剖面仪采集的回波信号采用相关脉冲对协方差法计算船速和流速的矢量和V的粗测值VS或精测值VP。
全文摘要
本发明公开了一种声多普勒海流剖面仪测量流速的方法，首先预估船速和流速的矢量和V的估计值V
文档编号G01S7/52GK1632592SQ200310122440
公开日2005年6月29日 申请日期2003年12月23日 优先权日2003年12月23日
发明者朱维庆, 潘锋, 朱敏, 汪玉玲 申请人:中国科学院声学研究所
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
