摘要
海底地形的测量是海洋环境调查、资源开发及航行保障等工作的基础,主要任务是测量海底深度,并进一步绘制海底地形图。多波束测深系统利用广角度定向发射和多通道接收技术,形成条带状的高密度水深数据,是一项高技术海底地形测量手段。近五十年来,多波束测深技术不断完善,并持续朝着宽覆盖、高分辨以及高精度的方向发展。
多波束测深系统作为一种水声探测设备,其性能最终取决于测量过程背后的物理机制。系统接收的海底后向散射回波信号在空间和时间上都是连续分布的,相对于传统声纳研究的点目标而言,回波信号严重展宽,底部对应的回波信号特征不明显,信号的分辨和特征参数的估计都更为复杂,对底部回波信号检测/估计性能的分析更具困难。
多波束测深中的底检测处理属于含有未知参量的检测问题,应用时回波参量估计与水底检测是联合进行的,一方面估计中需用到检验为真的信号样本,另一方面检测判决需要有效参量估计的支撑。底检测处理中的目标为水底,面临的问题主要包括到达角度/到达时间估计、后向散射信号检测和条带测量模型变化检测。
论文以提高边缘波束测量精度和水底测深数据可靠性为目的,以国家863计划课题中的多波束浅海地形测量系统为研究平台,通过估计-检测-模型联动的方式,围绕着影响测线的三个主要因素(测点位置、测点深度以及测线模型),对多波束测深中的关键技术——波束形成处理、底检测处理和测深数据的模型检测展开研究工作。
首先针对多波束测深系统的应用环境,构建了回波信号的传播模型。对于浅海环境,信号在近距离情况时以球面波的形式传播,此时按照平面波模型处理,将会导致测点位置出现偏差。论文发展了近场聚焦处理与远场常规处理相结合的波束形成方法,以提高波束定向的精确度以及测点的定位精度;针对近场聚焦算法计算量大,很难实时处理的问题,提出了近似的聚焦波束形成算法,并仿真研究了近似聚焦的波束模式图与理想情况的差别,验证了近似算法的有效性。
测点的深度由回波信号的到达时间和到达角度两个待估参量决定,其中回波到达角度通常根据预成的波束角度获得,对于边缘的倾斜波束,脚印明显展宽,导致测深精度严重下降。针对边缘波束脚印内地形起伏导致回波信号表现出的点状特征,论文采用高分辨力方位估计算法分辨波束宽度内的多向回波信号,获得亚波束宽度的分辨性能。具体提出了波束域多指向的基于旋转不变技术的信号参数估计方法(Estimation of Signal Parameters via Rotational Invariance Techniques,简写为ESPRIT),降低了此类算法对信噪比以及样本点数的要求,且计算复杂度降低。试验数据处理结果表明,与传统的加权平均时间(Weighted Mean Time,简写为WMT)或者方位偏离指示(Bearing Deviation Indicator,简写为BDI)底检测算法相比,波束域多指向ESPRIT算法能够提供更好的细节分辨效果。
另一方面,系统覆盖范围/分辨力要求、安装空间限制等因素使得多波束系统多采用特殊形状的接收阵列,导致基于均匀线阵假设的高分辨力方位估计算法无法直接应用,同时回波信号的短时平稳性使得算法所需的协方差矩阵的估计较为困难,增加了高分辨力方位估计算法的应用复杂度。论文通过将虚拟阵列变换、多子阵空间前后向平滑方法与波束域多指向ESPRIT算法相结合,得到适用于任意形状平面阵的高分辨力底检测算法,计算机仿真以及试验数据分析结果均验证了算法的有效性。
回波信号的到达时间和到达角度可转化为水底的测点数据。由于各种干扰的影响,测点数据中存在着异常值,同时测点分布的不均匀性也不利于获得连续的水底测线。论文提出利用最小二乘支持向量机对条带的测点数据进行建模,根据模型对测深数据进行筛选,剔除异常值,进一步提高模型的拟合精度。试验数据分析表明,对水底地形建模既可以填补测量间隙,又能够实现多种算法测量结果的融合处理,提高测点数据的可靠性。
论文在理论研究的基础上完成了WMT和BDI两种底检测算法的信号处理机的实现。通过多次水池测试、湖上以及海上试验的实际使用,验证了信号处理机及其算法实现的有效性和可靠性。
Information on bottom topography is the foundation of numerous human ocean activities including marine environmental investigation, resource exploitation and navigation security. The goal of a bathymetry measurement system is to measure the water depth and further map the seabed topography. As such a high-tech instrument, multi-beam echo sounders (MBES) can generate high-density strip depth measurements by exploiting wide-swath directional transmitting and multi-channel receiving. During the past 50 years, the MBES technology has been improved significantly; it continues to improve toward even wider coverage, higher resolution and better precision.
As an underwater acoustic system, performance of the MBES is ultimately determined by the physical mechanism of the measurement process. The backscattering signal of the seabed is continuously distributed in both space and time. Compared to the case of point objects studied in conventional sonar applications, detection and estimation of signal parameters in MBES are more complex, and the associated performance analysis is even more challenging.
The bottom detection process in MBES is a hypothesis testing problem with unknown parameters, i.e., signal detection is jointly done with parameter estimation. On one hand, data samples, for which the hypothesis testing is true, are needed in parameter estimation; on the other hand, detection is supported by effective estimation results. Here, the target is the seabed, and the problem facing includes target direction-of-arrival (DOA) and time-of-arrival (TOA) estimation, backscattering signal detection, and model-based change detection.
This thesis is to improve the depth accuracy of the outer beams and the reliability of the measurement results, using the shallow-water MBES system funded by the national 863 program as the research platform. Through an interactive linkage among detection, estimation and model, considering main factors impacting the measurement lines (the position and depth of measuring points, and the seabed model), the thesis research focuses on three key aspects of MBES signal processing:beamforming. bottom detection, and model-based bottom change detection.
Firstly, a seabed backscattering signal model is developed. For a shallow-water environment, the signal in the near distance propagates in the form of a spherical wave; in this case, use of the plane wave model will result in some deviation of the measured position. As such, a combination of the near-field focused beamforming and plane wave beamforming is presented to make the beams pointed more precisely and improve the localization accuracy. Considering the large amount of calculation required for focused beamforming, the thesis proposes an approximate algorithm. The difference between the approximate and the ideal beam patterns is discussed; and simulations validate the effectiveness of the approximate algorithm.
The depth of measuring point is determined by DOA and TOA estimation. The estimation of DOA can be obtained from the beam angles; however, broadening of the footprints of the outer inclined beams reduces the resolution significantly and affects the accuracy of depth estimation. In fact, the topographic relief within the footprints of the inclined beams makes the echo present point-like features. Hence the high resolution algorithms can be used to resolve simultaneously echo signals and obtain a sub-beamwidth resolution. In this thesis, estimation of signal parameters via rotational invariance techniques (ESPRIT) with multiple-angle subarray beamforming is developed to reduce the requirement for the signal-to-noise ratio (SNR), the number of snapshots, and computational efforts. Experimental data processing results show that, compared to the traditional bottom detection algorithms such as weighted mean time (WMT) and bearing deviation indicator (BDI), ESPRIT with multiple-angle subarray beamforming can provide better fine-structure resolution.
On the other hand, because of the system coverage/resolution requirement and the space constraint of installation, an MBES usually employs a receiving array of special shape, which leads to that the high resolution algorithms based on a uniform linear array (ULA) can not be applied directly. Meanwhile the short-term stationary of the echo signals makes it difficult to estimate the covariance matrix required, further increasing the complexity in applying a high-resolution algorithm. In this thesis, a high resolution bottom detection algorithm is developed, which can be used on a planar array of any shape, by combining the ESPRIT with multiple-angle subarray beamforming with virtual array transformation and multiple subarray forward-backward spatial smoothing. Computer simulations and experimental data processing results verify the effectiveness of the developed algorithm.
The DOA and TOA estimation can be converted to the horizontal distance and depth of a measuring point. Due to the impact of operational errors and external interferences, there are outliers in the results of bottom detection, while the uneven distribution of measuring points is not conducive to obtaining a continuous bottom profile. In the thesis, the LS-SVM (least squares support vector machine) approach is studied to model the bottom profile. Based on this model, the outliers are removed to further improve the fitting precision. Experimental data processing results show that the model can not only fill the measurement gap, but also achieve the fusion of the results with a variety of algorithms, and improve the reliability of the results.
During the thesis research, two kinds of bottom detection algorithms—WMT and BDI have been designed on a digital signal processor for the MBES. Numerous tank tests, lake and sea trials have demonstrated that the signal processor and its algorithms can work effectively and reliably.
引文
[1]Hammerstad E, Asheim S, Nilsen K, et al. Advance in multibeam echosounder technology [C]//Proceedings of MTS/IEEE OCEANS Conference.1993: 482-487.
[2]黄谟涛.多波束测深技术研究进展与展望[J].海洋测绘,2000,78(2):2-7.
[3]赵会滨,徐新盛,吴英姿.多波束条带测深技术发展动态展望[J].哈尔滨工程大学学报,001,22(2):42-45.
[4]International Hydrographic Bureau. IHO standards for Hydrographic surveys [S]. 5th ed.2008. http://www.iho.int.
[5]中船重工第七一五研究所.国家高技术研究发展(863)计划海洋技术领域海洋环境立体监测技术专题/项目——“多波束海底地形测量系统”课题实施方案[R].中国船舶重工集团公司,2007年8月.
[6]SeaBeam Instrument Inc. SeaBeam 2000 series 2100 multibeam survey system [J/OL].1993. http://www.seabeamcanada.ca.
[7]ELAC Nautik, BottomChart compact dual frequency MK Ⅱ50/180kHz [J/OL]. 1998. http://www.elac-nautik.de.
[8]Simrad Inc. Simrad EM3000 multibeam echosounder:operator manual [J/OL]. 1996. http://www.simrad.com.
[9]ATLAS. Product information of ATLAS fansweep 20 shallow water multibeam sweeping echosounder systems [J/OL].1997. http://www.atlashydro.atlas-elektronik.com.
[10]RESON Inc. Seabat 8100 series new-generation multibeam echo sounders [J/OL].1997. http://www.reson.com.
[11]哈尔滨工程大学水声研究所.条带测深仪技术方案报告[R].哈尔滨工程大学,1998年.
[12]de Moustier C. State of the art in swath bathymetry survey system [D]. Int. Hyd. Rev., Monaco, LⅩⅤ (2), July 1988.
[13]赵建虎,刘经南.多波束测深及图像数据处理[M].武汉大学出版社,2008.
[14]Urick R J. Principles of Underwater Sound [M].3rd ed. McGraw Hill,1983.
[15]Time-Pekka Jantti. Trials and experiment results of the ECHOS XD multibeam echosounder [J]. IEEE Journal of Oceanic Engineering,1989,14(1):306-313.
[16]Ol'chevskii V V. Characteristics of Sea Reverberation [M]. Consultant Bureau, 1967.
[17]Boehme H, Chotiros N P. Acoustic backscattering at low grazing angles from the ocean bottom [J]. Journal of the Acoustical Society of America,1988,84(9): 1018-1029.
[18]Ogilvy J A. Theory of Wave Scattering from Random Rough Surfaces [M]. Adam Hilger,1991.
[19]Van Trees H L. Detection, Estimation and Modulation Theory, Part III [M]. Wiley,1971.
[20]杜选民,李海森,姚蓝.多波束条带测深系统中正交信号的获取技术[J].声学技术,1998,17(1).
[21]Maranda B. Efficient digital beamforming in the frequency domain [J]. Journal of the Acoustical Society of America,1989,85(5):1818-1819.
[22]Hacker P S, Schrank H E. Range distance requirements for measuring low and ultra low sidelobe antenna patterns [J]. IEEE Transactions on Antennas and Propagation,1982, AP-30(9):956-965.
[23]Hansen R C. Measurement distance effects on low sidelobe patterns [J]. IEEE Transactions on Antennas and Propagation,1984, AP-32(6):591-594.
[24]Silver S. Microwave Antenna Theory and Design [M]. McGraw Hill,1949.
[25]Rhodes D R. On minimum range for radiation patterns [J]. Proceedings of the Institute of Radio Engineers (IRE),1954,9:1409-1410.
[26]Jiang Y, Xu W, Chen L, et al. Near-field beamforming for a multibeam echo sounder:Approximation and error analysis [C]//Proceedings of MTS/IEEE OCEANS Conference.2010.
[27]陈龙永,梁兴东,丁赤飚.基于查表法的快速求浮点数平方根方法[J].微计算机信息,2009,25(2-3):205-206.
[28]Yan S F, Ma Y L, Sun C. Optimal beamforming for arbitrary arrays using second-order cone programming [J]. Chinese Journal of Acoustics.2005,1(24): 1-9.
[29]Yang L, Taxt T, Albregtsen F. Bottom detection for multibeam sonars with active contours [C]//Proceedings of MTS/IEEE OCEANS Conference.1997,2(10): 943-950.
[30]Lurton X. Swath bathymetry using phase difference:Theoretical analysis of acoustical measurement precision [J]. IEEE Journal of Oceanic Engineering, 2000,25(3):351-363.
[31]Manolakis D Q Ingle V K, Kogon S M. Statistical and adaptive signal processing:Spectral estimation, signal modeling, adaptive filtering and array processing [M]. McGraw Hill,2000.
[32]de Moustier C. Beyond bathymetry:Mapping acoustic backscattering from the deep seafloor with the SeaBeam [J]. Journal of the Acoustical Society of America, 1986,79(2):316-331.
[33]陈非凡,吴英姿.一种快速能量中心收敛方法[C]//水声学学术会议.1997:121-123.
[34]Alexandrou D, de Moustier C. Adaptive noise canceling applied to SeaBeam sidelobe interference rejection [J]. IEEE Journal of Oceanic Engineering,1988, 13(2):70-76.
[35]de Moustier C. Angular dependence of 12kHz seafloor acoustic backscatter [J]. Journal of the Acoustical Society of America,1991,90(1).
[36]陈非凡,吴英姿,卢逢春.一种多波束声呐回波信号时延的实时特征窗分析方法[J].海洋技术,1998,17(4).
[37]Kraeutner P H, Bird J S. Beyond interferometry:Resolving multiple angles-of-arrival in swath bathymetric imaging [C]//Proceedings of MTS/IEEE OCEANS Conference.1999,1:37-45.
[38]Stergiopoulos S, Ashley A T. An experimental of split-beam processing as a broadband bearing estimator for line array sonar systems [J]. Journal of the Acoustical Society of America,1997,102(6):3556-3563.
[39]Masnadi-Shirazi M A, de Moustier C, Cervenka P, et al. Differential phase estimation with the SeaMARCⅡ Bathymetric Sidescan sonar system [J]. IEEE Journal of Oceanic Engineering,1992,17(3):241.
[40]Yang L, Taxt T. Multibeam sonar bottom detection using multiple subarrays [C]// Proceedings of MTS/IEEE OCEANS Conference.1997,2(10):932-938.
[41]周天,朱志德,李海森等.多子阵幅度-相位联合检测法在多波束测深系统中的应用[J].海洋测绘,2004,24(4).
[42]王永良,陈辉,彭应宁等.空间谱估计理论与算法[M].清华大学出版社,2004.
[43]Van Trees H L. Optimum Array Processing [M]. Wiley,2002.
[44]张贤达.现代信号处理[M].清华大学出版社,2003.
[45]Scharf L L. Statistical signal processing:Detection, estimation and time series analysis [M]. Addison-Wesley Publishers,1991.
[46]Capon J. High resolution frequency-wavenumber spectrum analysis [J]. Proceedings of IEEE,1969,57(8):1408-1418.
[47]Hellequin L. Statistical characterization of multibeam echosounder data [C]// Proceedings of MTS/IEEE OCEANS Conference.1998,1(10):228-233.
[48]Roy R, Kailath T. ESPRIT-Estimation of signal parameters via rotational invariance techniques [J]. IEEE Transactions on Acoustic, Speech, and Signal Processing,1989,37(7):984-995.
[49]Richard R H. ESPRIT:Estimation of signal parameters via rotational invariance techniques [D]. PhD thesis. Stanford University,1987.
[50]Bienvenu G, Kopp L. Decreasing high resolution technique sensitivity by conventional beamformer preprocessing [C]//Proceedings of International Conference on Acoustics, Speech, and Signal Processing.1984,33(2):1-4.
[51]Xu X L, Buckley K. Reduced-dimension beamspace broadband source localization:Preprocessor design and evaluation [C]//Proceedings of IEEE 4th Annual ASSP Workshop on Spectrum Estimation and Modeling,2008: 1148-1152.
[52]Li J. Improving ESPRIT via beamforming [J]. IEEE Transactions on Aerospace and Electronic systems,1992,28(2):520-528.
[53]Xu G, Silverstern S D, Roy R, et al. Beamspace ESPRIT [J]. IEEE Transactions on Signal Processing,1994,42(2):349-356.
[54]Natterer F. The Mathematics of Computerized Tomography [M]. Wiley,1986.
[55]Golub G H, Loan C F V. Matrix Computations [M], Johns Hopkins University Press,1989.
[56]Moura J M F, Lourtie I M G,Eds. Signal processing for swath bathymetry and concurrent seafloor acoustic imaging [J]. Acoustic Signal Processing for Ocean Exploration,1993:329-354.
[57]Friedlander B. Direction finding using spatial smoothing with interpolated arrays [J]. IEEE Transactions on Aerospace and Electronic Systems,1992,28(2): 574-586.
[58]Weiss A J, Gavish M. Direction finding using ESPRIT with interpolated arrays [J]. IEEE Transactions on Signal Processing,1991,39(6):1473-1478.
[59]杨鹏,杨峰,聂在平.基于子阵分割和虚拟内插ESPRIT算法的圆柱共形阵DOA估计试[A].全国天线年会论文集(下),2009.
[60]Ruilin L, Dejun Z. Estimation via exploitation of circular array with the interpolation and spatial smoothing [C]//3rd International Conference on Signal Processing.1996,10:509-512.
[61]Ronhovde A. High resolution beamforming of SIMRAD EM3000 bathymetric multibeam sonar data [D]. Cand Scient thesis. Universitas Osloensis,1999.
[62]Rao B D, Hari K V S. Effect of spatial smoothing on state space methods/ESPRIT [J]. IEEE ASSP 5th Workshop on Spectrum Estimation and Modeling,1990:377-381.
[63]Akaike H. A new look at the statistical model identification [J]. IEEE Transactions on Automatic Control,1974,19(6):716-723.
[64]Achwartz G. Estimation the dimension of a model [J]. The Annals of Statistics, 1978,6:461-464.
[65]Zhao L C, Krishnaiah P R, Bai Z D. Remark on certain criteria for detection of numbers of signals [J]. IEEE Transactions on Acoustic, Speech, and Signal Processing,1987,35(1):129-132.
[66]李家彪等.多波束勘测原理技术方法[M].海洋出版社,1999.
[67]赵建虎.多波束深度及图像数据处理方法研究[D].博士论文.武汉大学,2002.
[68]宋玲玲.多波束测量数据处理研究[D].硕士论文.南京航空航天大学,2007.
[69]吴自银等.多波束测深边缘波束误差的综合校正[J].海洋学报,2005,27(4):88-94.
[70]吴英姿.多波束测深系统地形跟踪与数据处理技术研究[D].博士论文.哈尔滨工程大学,2001.
[71]Vapnik V N. The Nature of Statistical Learning Theory[M]. Springer-Verlag, 1995.
[72]Sundance Multiprocessor Technology Limited. User manual for SMT148FX [J/OL].2006. http://www.sundance.com.
[73]Rao C R. Information and accuracy attainable in the estimation of statistical parameters [J]. Bull Calcutta Math. Soc,1945,37(1).
[74]Kraeutner P H, Bird J S. Principle components array processing for swath bathymetric mapping [C]//Proceedings of MTS/IEEE OCEANS Conference. 1997:1264-1254.
[75]Bird J S, Mullins J K. Analysis of swath bathymetric sonar accuracy [J]. IEEE Journal of Oceanic Engineering,2005,30(2):372-390.
[76]Bird J S, Kraeutner P H. Cramer-Rao bound investigation of swath bathymetric accuracy [C]//Proceedings of MTS/IEEE OCEANS Conference,2001,3: 1640-1647.
[77]Abel J S. A bound on mean-square-estimate error [J]. IEEE Transactions on Information Theory,1993,39(5):1675-1680.
[78]Xu W. Cramer-Rao bound for bearing estimation with bias correction [C]// Proceedings of MTS/IEEE OCEANS Conference.2007.
[79]Naftali E, Makris N C. Necessary conditions for a maximum likelihood estimate to become asymptotically unbiased and attain the Cramer-Rao lower bound. Part I. General approach with an application to time delay and Doppler shift estimation [J]. Journal of the Acoustical Society of America,2001,110(4): 1917-1930.
[80]Thode A, Naftali E, Ingram I, et al. Necessary conditions for a maximum likelihood estimate to become asymptotically unbiased and attain the Cramer-Rao lower bound. Part Ⅱ.Range and depth localization of a sound source in an ocean waveguide [J]. Journal of the Acoustical Society of America, 2002,112(5):1890-1908.
[81]Jiang Y, Xu W, Pan X. Barankin bound for bearing estimation with bias correction [C]//Proceedings of MTS/IEEE OCEANS Conference.2008.
[82]Xu W, Chen Q, Jiang Y. Performance bound approximation for bearing estimation with bias correction [J]. IEEE Signal Processing Letters,2009,16(10): 833-836.