The filtering and compressing of outer beams to multibeam bathymetric data

Some errors and noises are often present in multibeam swath bathymetric data. Echo detection error (EDE) is one of the main errors. It causes the depth error to become bigger in outer beams and looks like sound refraction. But depth errors due to EDEs have a trumpet-shaped appearance, instead of a curved appearance that is caused by the sound refraction errors. EDEs, including systematic acoustic signal detection errors and internal noises, cannot be removed during the correction of sound refraction. It causes depth inconsistencies between adjacent swaths and degrades precision of outer beams. Sometimes, the bathymetric errors caused by EDEs do not even meet the requirements of IHO (International Hydrographic Organization). Therefore, a post-processing method is presented to minimize the EDEs by filtering outliers and compressing outer beams of multibeam bathymetric data. The outliers caused by internal noises are removed by an automatic filter algorithm first. Then the outer beams are compressed to reduce systematic acoustic signal detection errors according to their depths, the calculated depth line and standard deviations (SDs). The automatic filter process is important for calculating the depth line. The selection of inner beams to calculate the average SD of beam depths is crucial to achieving compressing goals. The quality of final bathymetric data in outer beams can be improved by these steps. The method is verified by a field test.     

深海主跃层环境对会聚区影响的数值分析

应用BELLHOP模式,对声速剖面的声跃层结构变化引起会聚区偏移特征进行了分析。结果表明,声速垂直结构的变化可导致会聚区位置出现不同程度的偏移:跃层强度增加0.01 s-1将使会聚区向远离声源方向偏移1.5~2.0km;跃层厚度增大50m将使会聚区向靠近声源方向偏移0.3~0.5km;跃层位置加深100m将使会聚区向远离声源方向偏移0.5~1.0km。在跃层的三个特征量中,跃层强度起主导作用。跃层强度变化引起的声线在海洋次表层的偏折差异,进而导致进入深海等温层的入射角差异,是使会聚区发生偏移的决定性因素。     

A model for variations in the range and depth dependence of thesound speed and attenuation induced by bubble clouds under wind-drivensea surfaces

When modeling sound propagation through the uppermost layers of the ocean, the presence of bubble clouds cannot be ignored. Their existence can convert a range-independent sound propagation problem into a range-dependent one. Measurements show that strong changes in sound speed and attenuation are produced by the presence of swarms of microbubbles which can be depicted as patchy clouds superimposed on a very weak background layer. While models suitable for use in acoustic calculations are available for the homogeneous bubble layer (which results from long time averages of the total bubble population), no similar parameterizations are available for the more realistic inhomogeneous bubble layer. Based on available information and within the framework of a classification scheme for bubble plumes proposed by Monahan, a model for the range and depth dependence of the bubbly environment is developed to fill this void. This model, which generates a possible realization of the bubbly environment, is then used to calculate the frequency-dependent change in the sound speed and attenuation induced by the presence of the bubble plumes. Time evolution is not addressed in this work     

一种提高深远海全深度声速剖面重构精度的方法

开展多波束水深测量应同步进行声速剖面探测。因海上作业条件恶劣、作业时间受限及设备性能局限等影响,在深远海海域常获取不到全深度的实测声速剖面。尽管利用温盐场模型可将声速剖面直接延拓至实地水深的最大深度,但这种气候态平均声速剖面与实际的声速剖面间存在不可控的系统性偏差,会给声速改正及水深测量成果带来质量隐患。给出了一种提高深远海全深度声速剖面重构精度的方法,即利用有效探测深度附近的实测温度盐度值,对大于有效探测深度的各水层的模型温度盐度值施加程度不一的约束控制。结果表明,经优化后全深度声速剖面的重构精度得到明显提高,其中2个XCTD站点声速剖面的互差SSPD分别由-2.5～1.0 m/s优化为0.0～1.0 m/s、0.0～2.6 m/s优化为-1.5～0.0 m/s; 2个CTD站点声速剖面的互差SSPD分别由-0.5～1.7 m/s优化为-0.4～0.3 m/s、-2.15～0.8 m/s优化为-1.4～0.8 m/s。     

潮间带低速层中高速埋伏岩体计算模式的建立

在某些沿海潮间带区域的低速层中常存在一定数量的巨大高速滚石埋伏岩体。这些埋伏岩体会给工程地质勘察带来很大的困难 ,甚至无法进行工程地质勘察。根据浅层折射波方法原理利用浅层折射波相遇时距曲线与折射波射线路径关系建立了求取浅层高速埋伏岩体的埋深及其厚度的计算模式 ,利用这些公式可以较准确地确定高速埋伏岩体的位置、埋深及其厚度 ,对工程地质勘察是十分有利的。     

浅海双跃层脉冲声传播

利用在东海测量的双跃层声速剖面和修改的单跃层声速剖面,数值模拟了2种跃层条件下不同收发深度声脉冲传播的波形。模拟结果表明,当声源或接收器位于上混合层时,信号波形在2种条件下都出现梳状多途结构。当声源和接收器都位于下混合层时,信号波形在2种条件下均相似。当声源位于中间均匀层时,信号波形在除上混合层以外的4层都有显著差异。用简正波的深度-简正波号域的幅度和相应的群速度解释了双跃层和单跃层声速剖面条件下信号波形特点以及异同的原因。     

Seasonal variations of sound speed in the Arabian Gulf

Sound-speed computations from CTD casts in the Arabian Gulf during 1992, reveal spatial and temporal variations in acoustic properties. Hydrographic conditions affecting sound speed propagation were seasonally investigated. A monotonic decrease in sound speed profiles with depth was commonly observed at almost all the stations in the Gulf. However, an exception occurred at Hormuz strait during winter. The water exchange pattern between the Gulf of Oman and the Arabian Gulf seems to influence the sound-speed structure, especially in the southern part of the latter. Winter profiles along the Gulf axis showed almost vertically homogenous sound speed. Maximum speeds are observed in summer, with a strong gradient associated with the development of the summer thermocline layer. Horizontal distributions in both winter and summer show a decreasing trend in sound speed from the Strait of Hormuz to the head of the Gulf. The resultant profiles provide a more comprehensive and reliable data set than any that have been reported in the literature. Shallowness and multiple refraction and reflection in the Arabian Gulf may cause sound speed energy to be trapped. No sound channel was detected inside the Gulf. A correlation analysis shows that sound speed is closely correlated with temperature throughout the Gulf, except in winter in the southern half where salinity effects, as a result of inversion and water exchange at the entrance, are found to be dominant.     

基于二次多项式拟合的超短基线声线修正方法

超短基线定位解算中的距离观测值是指换能器与水下应答器之间的直线距离,而海水声速的不均匀分布导致声波在海水中的实际传播路径为连续弯曲的曲线,需要结合实测声速剖面进行声线修正。根据声速在分层介质中的传播特性,本文提出了一种基于二次多项式拟合的声线跟踪算法,采用线性插值方法对声速剖面数据进行合理加密并按等深度进行分层,设定每层声速梯度是不断变化的,用二次多项式拟合声速,基于运动学原理建立了完整的数学解算模型。仿真结果表明,该方法修正后的水下目标分布具有明显的收敛性,且优于等梯度声线跟踪算法和等效声速剖面法,显著提高了超短基线水声定位系统的定位精度。     

多波束测深残余折射处理技术的对比研究

针对多波束条幅存在的残余声线折射影响,加拿大Caris HIPS多波束资料处理软件除了支持声速剖面改正外,还提供了一个折射改正工具。该折射改正工具简单易用,效果所见即所得,能够快速消除声线折射产生的假地形。通过实测数据对这两种改正方法的结果进行了比较。     

Basement imaging with Side-Looking Seismics

The multistreamer Side-Looking Seismic system presented in this paper makes a sonograph of uncovered or buried crustal topography, thus revealing the structural fabric of the oceanic basement, even when this is covered with a sedimentary layer. Major elements of the system are an airgun as a sound source, five single-channel parallel streamers and two minicomputers for signal capture and processing.The system is used simultaneously for enhanced single-channel seismic profiling and for side-looking seismics. A vertical section with an improved signal-to-noise ratio and a suppression of side-echoes is produced on a digital seismic recorder. Primary side-looking seismic output in the form of 5 profiles with different angles of incidence is obtained within 10 seconds. This part of the processing can be done in real time.In sediment-covered areas the low frequencies used cause the slanted profiles (the side beams in the primary output) to be side-looking sonar images of buried topography. The projection process yielding final side-looking output corrects for slant range deformation caused by the water column and, if necessary, for deformation caused by refraction within the sedimentary column. The result approaches a conformal map of the structure of the traversed basement. Swath width is mainly determined by water depth and refraction effects in the sediment. In Madeira abyssal plain a swath width of 8000 m was attained in a water depth of 5000 m.Within the swath, oceanic basement structures are recognized in the form of elongate more or less parallel reflectors. They are interpreted as buried spreading topography. The lack of side-echoes within fracture zones combined with typical wall signatures can be used to trace fracture zones. These features are demonstrated for an area in Madeira abyssal plain.     

菲律宾海的声速剖面结构特征及季节性变化

应用Argo资料研究了菲律宾海的声速剖面结构特征。通过统计分析选取了合理的跃层标准，分析了主跃层、季节性跃层和表面正梯度层的区域性分布及季节性变化。结果表明，菲律宾海主要受赤道流系和北太平洋西边界流系的支配，其环流结构和水团配置对声场结构影响很大；主跃层的经向差异显著，但季节性变化较小，其平均位置由南向北逐渐加深，强度逐渐减弱；季节性跃层的分布及变化主要受混合层的季节性变化以及北部海区冬季温跃层通风过程的影响，夏季较强较厚，冬季较弱较薄；深海声道轴季节性变化较小，南极中层水和北太平洋中层水的温盐差异是其经向分布差异的主要原因。综合考虑海区声速结构区域性和季节性特征，将其归纳为6种典型结构，得出了各类声速剖面的模态特征及垂直结构参数的统计特征值。     

Correction for depth biases to shallow water multibeam bathymetric data

Vertical errors often present in multibeam swath bathymetric data. They are mainly sourced by sound refraction, internal wave disturbance, imperfect tide correction, transducer mounting, long period heave, static draft change, dynamic squat and dynamic motion residuals, etc. Although they can be partly removed or reduced by specific algorithms, the synthesized depth biases are unavoidable and sometimes have an important influence on high precise utilization of the final bathymetric data. In order to confidently identify the decimeter-level changes in seabed morphology by MBES, we must remove or weaken depth biases and improve the precision of multibeam bathymetry further. The fixed-interval profiles that are perpendicular to the vessel track are generated to adjust depth biases between swaths. We present a kind of postprocessing method to minimize the depth biases by the histogram of cumulative depth biases. The datum line in each profile can be obtained by the maximum value of histogram. The corrections of depth biases can be calculated according to the datum line. And then the quality of final bathymetry can be improved by the corrections. The method is verified by a field test.     

黄海、渤海温度垂直结构三维数值预报

文中将海洋分为具有密度阶跃的两层，设各层内各物理量与深度无关，用推广的ＡＤＩ方法进行包括潮流和风海流的二层流场计算，待流场计算稳定后，与温度扩散方程和上均匀层深度预报方程相耦合，对黄渤海区表、底层温度和上均匀层深度作了为期４ｄ的试报。然后，通过温度垂直剖面自模函数预报出三维的水温场。试报结果表明，整个模式的预报性能及试报结果与实测的吻合程度是令人满意的。     

海底沉积物的基本地声结构与地声模型

为准确建立海底地声模型,本文探讨地声模型的基本组成和基本结构。通过样品实验室测量,分析南海海底表层沉积物的密度、孔隙度与声速随着埋深变化的关系,得出海底实际存在的低声速表面–声速缓慢变化类型、低声速表面–声速增大类型、高声速表面–声速缓慢变化类型和高声速表面–声速增大类型4种典型地声结构;对比钻探测量,分析黄海海底沉积物的密度、孔隙度与声速随埋深变化关系,得出海底地声模型分层特征与地声结构组合特征。研究表明,地声模型可以归结为4种基本地声结构的组合,通过与底层海水声速、同层内声速剖面以及与上层海底沉积物下表面声速的比较,可以建立各种海底地声模型;基于实验室测量法建立的地声模型可以作为参考地声模型,但需要考虑实际海底温度和压力梯度以及海底沉积物的频散特性等,借助于声速比校正法和频散性理论模型进行计算及修正。     

Acoustic environment and physico‐mechanical properties of the shelf seabed off the pearl river mouth

Abstract

The vast shallow sea off the Pearl River mouth in the northern South China Sea is an important prospecting area for offshore oil development. In recent years, the authors have investigated acoustic and geotechnical characteristics of marine sediments in this area. An intercalated layer of low sound velocity and low compressive strength has been found within the seabed, in which the median diameter of sediment grains is fine and the sound velocity is 100–200 m/s lower than that of the overlying and underlying layers. The minimum unconfined compressive strength of this layer is 0.075 kg/cm², which is lower than that of the over‐ and underlying layers by an order of magnitude. Such an intercalation often constitutes a threat to the stability of shallow foundation soil. In case of overloading, the layer may be weakened, and seafloor sliding between different sediment layers may occur. The regional distribution of these kinds of weak intercalations of low sound velocity may be traced by a subbottom profiler and by means of sediment acoustical investigations.

Correlation between the gray level of a layer on subbottom profile records and physical properties of the layer (including sound velocity and reflectivity) suggests that the layer of sufficient bearing capacity must be searched by means of sediment acoustics at least to the depth of a high‐velocity substratum of stronger reflection.     

由正入射遥测海底介质参数方法初探

通过数值模拟计算发现，对于分层有吸收海底，在层厚已知的条件下，除最下面的介质以外．有可能利用垂直入射的脉冲声波遥测各层的密度、声速与吸收系数。此法的要点是，根据入射波和层参数的假设值算得反射信号，与实际反射信号比较；然后改变参数再作计算，直到误差最小为止．此时对应的层参数即为介质参数值实际的估计值。指出根据这种方法有可能实现以自适应方法遥测分层海底的介质参数。     

Using matched-field processing to estimate shallow-water bottomproperties from shot data taken in the Mediterranean Sea

It is extremely difficult to determine shallow ocean bottom properties (such as sediment layer thicknesses, densities, and sound speeds). However, when acoustic propagation is affected by such environmental parameters, it becomes possible to use acoustic energy as a probe to estimate them. Matched-field processing (MFP) which relies on both field amplitude and phase can be used as a basis for the inversion of experimental data to estimate bottom properties. Recent inversion efforts applied to a data set collected in October 1993 in the Mediterranean Sea north of Elba produce major improvements in MFP power, i.e., in matching the measured field by means of a model using environmental parameters as inputs, even using the high-resolution minimum variance (MV) processor that is notoriously sensitive and usually results in very low values. The inversion method applied to this data set estimates water depth, sediment thickness, density, and a linear sound-speed profile for the first layer, density and a linear sound-speed profile for a second layer, constant sound speed for the underlying half space, array depth, and source range and depth. When the inversion technique allows for the array deformations in range as additional parameters (to be estimated within fractions of a wavelength, e.g., 0.1 m), the MFP MV peak value for the Med data at 100 Hz can increase from 0.48 (using improved estimates of environmental parameters and assuming a vertical line array) to 0.68 (using improved estimates of environmental parameters PLUS improved phone coordinates). The ideal maximum value would be 1.00 (which is achieved for the less sensitive Linear processor). However, many questions remain concerning the reliability of these inversion results and of inversion methods in general     

Heterogeneous distribution of organochlorine contaminants in the blubber of baleen whales: implications for sampling procedures

The concentration of DDTs and PCBs was determined in two strata (inner and outer) of the blubber of fin and sei whales from the eastern North Atlantic. In the two species residue levels of all organochlorines in the outer layer were significantly higher than in the inner one. The DDE/tDDT ratio did not show significant variation between strata, but the tDDT/PCB ratio was higher in the outer layer than in the inner one. Such differences between strata cannot be explained by variation in the quantity of lipids present in the blubber and probably indicate a constant organochlorine distribution pattern in the blubber of large whales. The differential role of the blubber strata in the fattening cycle, the heterogeneous lipid composition of blubber throughout its depth and the different turnover rates of pollutants in the layers are the probable factors responsible for these observed differences in pollutant distribution between blubber strata. Because of this heterogeneous distribution, blubber samples collected from cetaceans for pollutant analyses should include equal representation of all layers in order to be representative of the individual's pollutant load. This is especially important when sampling stranded specimens or taking biopsy samples from large free-ranging whales.     

台风“天秤”对两种深海声道下声传播的影响

将Argo浮标资料与卫星遥感再分析数据相结合,调用基于抛物方程算法的RAM(Range-dependent Acoustic Model)模型,研究了2012年第14号台风“天秤”对不完整深海声道(3 000 m)和完整深海声道(5 500 m)两种水深条件下声传播特性的影响。结果显示:台风对海水的影响局限于表层水体,具体为混合层加厚,混合层内温度梯度接近于零,声速在混合层内正梯度分布;混合层下方一定深度的水体增温,相应的声速也增大。声源在混合层内时,主要对海表层的声传播产生影响,两种水深条件下均出现表面波导声传播模式以及泄漏模式。声源在混合层以下时,不完整深海声道条件下台风使得会聚区向着声源方向靠近;完整深海声道条件下台风对会聚区的位置影响不明显,但声波的翻转深度增加近500 m。     

两种基于贝叶斯点估计理论的多声源定位方法研究

海洋环境参数失配是制约匹配场定位性能的主要因素之一。为了克服环境失配，本文基于贝叶斯理论，将环境参数与声源的距离和深度一起作为未知量进行反演。然而在进行多声源定位时，反演参数的维数几何增长，极大地增加了反演问题的复杂性和计算量。为此本文将声源强度和噪声方差表示成其极大似然估计值，从而将这些参数进行隐式采样，大大降低了反演的维数和难度。文章比较了两种贝叶斯点估计方法，最大后验概率密度方法和最大边缘后验概率密度方法。最大后验概率密度方法的解是令后验概率密度取得最大值的参数组合，可以利用优化算法快速获得。最大边缘后验概率密度法将其他参数积分，得到目标参数的一维边缘概率分布，分布的最大值为反演结果。该方法得到最优估计值的同时可以获取参数估计的不确定信息。在环境参数和声源参数都未知的情况下，利用蒙特卡洛法在不同信噪比情况下对两种声源定位方法进行分析，实验结果表明：(1)对于敏感参数，如声源距离、水深和海水声速，最大边缘后验概率密度法比最大边缘后验概率密度方法的性能好。(2)对于较不敏感的参数，如海底声速、海底密度和海底声衰减，当信噪比较低时，最大边缘后验概率密度方法能较好地平滑噪声，从而比最大边缘后验概率密度法具有更好的性能。由于声源距离和深度是敏感参数，研究表明最大边缘后验概率密度法提供了一种在不确知环境下更可靠的多声源定位方法。