海洋声层析研究的两个热点问题

综述了目前海洋声层析研究的两个热点问题，归纳了用声学方法监测气候变化和反演内波两个领域所取得的研究进展，总结了研究过程中存在的问题。     

射线稳定性参数与波导不变量

基于简正波理论的波导不变量在声传播特性描述中已得到广泛应用,而当环境与距离无关或存在轻微扰动时,基于射线理论的声场特性可以用一个新的物理量—射线稳定性参数来描述。本文通过简正波干涉声场数值模拟,通过分析简正波相长干涉形成的能量包传播路径与射线传播路径相互关系,利用简正波水平干涉结构距离上的周期性与射线水平跨距相等的特点,采用简洁方式证明了射线稳定性参数与简正波波导不变量间的等价性,从而在射线和简正波理论中建立了一个新的联系。对远距离射线模型声传播仿真结果的分析表明,与射线传播时间相比,射线稳定性参数作为射线理论的一个新的参数,可有效用于描述声传播特性,并能够反映海洋环境变化,应用于海洋环境参数声学反演。     

海洋声层析研究的两个热点问题

综述了目前海洋声层析研究的两个热点问题 ,归纳了用声学方法监测气候变化和反演内波两个领域所取得的研究进展 ,总结了研究过程中存在的问题     

浅海内波及声场起伏数值研究

为探讨声学方法监测(反演)海洋内波的可行性,分析ASIAEX2001亚洲海海洋声学实验定点温度链实验数据,结果表明：实验期间存在M2内潮波、周期10～30min非线性内波现象,并观测到双温跃层同相传播“双非线性内波”。分别对垂直和平行于内波传播方向传播声场进行了数值计算分析。垂直于内波传播方向的声传播计算结果表明：简正波波数时间演化函数的频谱与内波引起的声速剖面变化频谱基本一致;简正波系数(模值)时间演化函数的频谱与内波引起的声源位置处的声速演化函数的频谱基本一致;高阶简正波波数差随时间的变化要比低阶简正波波数差小约1个量级。平行于内波传播方向的声传播计算结果表明：内潮波主要导致相邻简正波间耦合;非线性内波包能够导致跨号甚至垮多号简正波间耦合。讨论了2种海洋内波声学监测方法。     

胶州湾流速场的声层析反演研究

采用结合匹配法和经验正交函数法的射线声层析反演方法,针对流场水平分布不均匀特点,使用距离等效分段方法对反演算法进行改进,利用三个断面声层析数据,对胶州湾口潮流场的垂直和水平分布进行反演。与传统方法比较,大大降低与实测流速垂直分布间的偏差,平均偏差小于0.02 m/s,流场垂直分布间的相关系数提高到0.85以上。分析了声传播断面上不同的流场结构,以及涨潮期间在胶州湾团岛附近出现的涡流现象,并计算得到湾口海水流量变化。结果表明,改进的射线声层析反演方法可有效地用于水平分布不均匀流场监测,仅采用少量声学观测站位,即可获得大范围的复杂流场三维信息,有利于近海海洋资源开发、海洋环境保护和船舶航行安全。     

海洋声学测温及声层析法

海洋声层析法作为检测海洋内部变化一有力手段，正面临着海洋声学大尺度观测新挑战。本文就海洋声传播的绝热近似；声传播时间的非线性偏移；声速剖面的线性及非线性反演等与声学大尺度观测相关的一些水声学问题的发展现状及存在问题等进行了综述性论述。文中还介绍了赫德岛预演实验等利用声层析法观测海洋内部变化的一些具体应用例     

BDRM理论在双轴海洋声道中的应用

由于表面声道与深海声道间的耦合效应,声波在双轴海洋声道中的传播比较复杂,因此求解双轴海洋声道中的声场就比较困难。在 WKBZ 本征函数的基础上,推导出了参考界面相位修正的一致表达式,并将浅海声传播的波束位移射线简正波(BDRM)理论应用于计算双轴海洋声道中的声场,进行了数值模拟并与传统简正波方法进行比较,结果表明应用 BDRM 理论计算的传播损失具有很高的精度和速度,可以对双轴海洋声道内声传播问题进行分析和预报。     

海洋声层析观测技术和方法

叙述海洋声层析观测系统,以声线传播时间层析为重点概括了海洋声层析的基本原理和其他主要方法,共6个方法。对运用海洋声层析观测来反演海洋状态问题的建立、求解及其误差来源作了分析和讨论。以测量声线传播时间为例介绍了海洋声层析观测系统主要设计技术。     

射线声学法在会聚区调查中的应用研究

会聚区调查是海洋声传播损失调查的一项重要内容,在调查前大致了解海区中的会聚区分布情况,将给会聚区调查带来极大的便利.针对深海会聚区调查,提出了一种利用射线声学法对会聚区进行预估,从而大大提高了工作效率,降低了调查成本,保证了调查数据的质量.     

声学与海洋沉积学交叉领域的研究

声学与海洋沉积学交叉领域研究可分为“沉积层声学特性的研究”、“海底高频声散射或低频声反射与底质类型之间关系的研究”、“回声参数反演海底类型技术”和“海底回声图象识别海底沉积类型技术”4个方面的研究,较详细论述了4个方面的国内外研究现状,最后提出了该领域进一步研究的战略思路。     

The coupled oceanographic–tomographic analysis and prediction system: A practical means for utilizing high frequency tomography in the world's oceans

Details are presented of a methodology that utilizes acoustic travel time information in an ocean circulation model. Recent developments of this model-oriented tomography are discussed, representing some significant improvements over earlier formulations. More accurate means of determining the arrival times of specific ray paths are detailed, along with a means of estimating possible errors in the calculated travel times. The assimilation of the observed arrival time information into an ocean model is achieved using a Kalman gain, and more advanced expressions for calculating the Kalman gain are presented. A formulation to account for errors in the stated positions of a source and receiver is also presented. It is shown that the methodology performs fairly well in reproducing observed travel time anomalies. However, the model-predicted anomalies along a specific ray path may not always track the observed anomalies for that path when assimilating multiple ray path data. Results indicate that additional work is required to determine a means of handling observed arrival time data without having prior knowledge of the magnitude of errors in the observations. Results from simulation experiments provide estimates of: (1) potential errors when the travel times for ray paths are only sampled at discreet intervals as opposed to continuously and (2) to what degree acoustic data can be expected to “correct” model-predicted fields.     

海洋声学层析研究现状与展望

海洋声层析是1979年美国科学家W.Munk等人提出的,通过测量传播时间等声传播信号有关参量反演声波所穿越的海洋特征,得到大面积海域中的海洋动力学状态及其变化的估计。声层析方法自提出以来受到各国的重视。在最初的20 a间,以全球测温计划(ATOC)为代表的一系列验证实验推动了声层析的理论研究、设备研制和应用。步入21世纪后,声层析与海洋动力学的数据同化、利用滑翔机等设备的移动声层析以及沿海声层析得到重视,并进行了相关的理论与实验研究。纳入成像旗下的声层析希尔伯特方法、有效低功耗的分布式传感网络声层析,以及被动声层析,或将成为声层析未来研究的重要方向。     

基于声传播时间的二维流场反演数值仿真研究

基于南海1998年夏季调查航次诊断计算的流函数场,选取越南以东偶极子发生海域,进行不同的声层析观测站位设置实验。模拟计算声线传播时间信息,然后应用基函数重建方法进行了流函数场的模拟反演研究,讨论了不同随机观测误差对反演结果的影响。研究结果表明该方法是可行的,在所选取的约833 km×833 km海域内,在观测海域外围配置19个声层析观测站位就能够很好地重构原流函数场,空间分辨率约为63.7 km,可以分辨模拟海域中尺度涡场结构,这是传统的基于海流计(CM)或声学多普勒剖面仪(ADCP)等观测方法所不能企及的。该观测技术和方法可以实际应用于近海大范围流场结构的遥感实时监测,为近海污染物等的扩散研究,海洋环境变化等提供实时观测流场资料。     

Acoustic tomography at basin scales and clock errors

A basin-scale acoustic tomography simulation is carried out for the northeast Pacific Ocean to determine the accuracy with which time must be kept at the sources when clocks at the receivers are accurate. A sequential Kalman filter is used to estimate sound-speed fluctuations and clock errors. Sound-speed fluctuations in the simulated ocean are estimated from an eddy-resolving hydrodynamic model of the Pacific forced by realistic wind fields at daily resolution from 1981-1993. The model output resembles features associated with El Nino and the Southern Oscillation, as well as many other features of the ocean's circulation. Using a Rossby-wave resolving acoustic array of four fixed sources and twenty drifting receivers, the authors find that the percentage of the modeled ocean's sound-speed variance accounted for with tomography is 92% at 400-km resolution, regardless of the accuracy of the clocks. Clocks which drift up to hundreds of seconds of error or more for a year do not degrade tomographic images of the model ocean. Tomographic reconstructions of the sound-speed field are insensitive to clock error primarily because of the wide variety of distances between the receivers from each source. Every receiver “sees” the same clock error from each source, regardless of section length, but the sound-speed fluctuations in the modeled ocean cannot yield travel times which lead to systematic changes in travel time that are independent of section length. The Kalman filter is thus able to map the sound-speed field accurately in the presence of large errors at the source's clocks     

Antipodal acoustic thermometry: 1960, 2004

On 21 March 1960, sounds from three 300-lb depth charges deployed at 5.5-min intervals off Perth, Australia were recorded by the SOFAR station at Bermuda. The recorded travel time of these signals, about 13,375 s, is a historical measure of the ocean temperature averaged across several ocean basins. The 1960 travel time measurement has about 3-s precision. High-resolution global ocean state estimates for 2004 from the “Estimating the Circulation and Climate of the Ocean, Phase II” (ECCO2) project were combined with ray tracing to determine the paths followed by the acoustic signals. The acoustic paths are refracted geodesics that are slightly deflected by either small-scale topographic features in the Southern Ocean or the coast of Brazil. The refractive influences of intense, small-scale oceanographic features, such as Agulhas Rings or eddies in the Antarctic Circumpolar Current, greatly reduce the necessary topographic deflection and cause the acoustic paths to meander in time. The ECCO2 ocean state estimates, which are constrained by model dynamics and available data, were used to compute present-day travel times. Measured and computed arrival coda were in good agreement. Based on recent estimates of warming of the upper ocean, the travel-time change over the past half-century was nominally expected to be about −9 s, but little difference between measured (1960) and computed (2004) travel times was found. Taking into account uncertainties in the 1960 measurements, the 2004 ocean state estimates, and other approximations, the ocean temperature averaged along the sound channel axis over the antipodal paths has warmed at a rate less than about 4.6 m °C yr⁻¹ (95% confidence). At this time, the estimated uncertainties are comparable in size to the expected warming signal, however.     

Oceanic sound speed profile inversion

A modal (full-wave) method has been developed to predict ocean sound speed profiles from propagated acoustic field data. The method assumes a point source of sound in the ocean and uses as data the values of the transmitted acoustic field at an array. The formalism for depth-dependent sound speeds consists of the standard Hankel integral transform of the depth solution. In the travel length coordinate, the latter is written exactly, using the Green's function, in terms of an integral equation whose kernel includes the sound speed profile correction. A Born approximation to this equation is used. This is just the WKB solution, and permits the use of a nontrivial input (or guess) profile, here chosen as bilinear. The use of asymptotic methods enables us to write the data as an integral transform over the profile correction. The transform can be inverted. An example is presented for full-bandwidth inversion.     

基于正交匹配追踪的声层析方法

声速剖面的变化会对声传播产生较大的影响，经验正交函数模型经常用来实现对声速剖面数据的简化描述。然而在内波、湍流等海水不均匀性存在时，这种正则化操作会造成声速重构精度的大幅降低。本文利用字典学习生成声速剖面的非正交原子，在稀疏编码时采用正交匹配追踪（OMP，Orthogonal Matching Pursuit）算法，更新字典则使用KSVD （Kernel Singular Value Decomposition）的字典更新算法。由于字典学习不需要强制使用正交条件，对于训练数据更加灵活，从而可以使用少数的原子组合达到更高的重构精度。利用一次浅海声学实验多次测量的声速剖面研究了海水声速剖面的经验正交函数表示和字典学习，研究表明：相比于正交函数表示，学习字典可以利用少数原子（甚至一个原子）更好的解释声速剖面扰动。字典学习可以提高声速剖面的稀疏性，从而提高声速剖面的反演精度。     

A fast inversion method for ocean parameters based on dispersion curves with a single hydrophone

The dispersion characteristics of shallow water can be described by the dispersion curves, which contain substantial ocean parameter information. A fast ocean parameter inversion method based on dispersion curves with a single hydrophone is presented in this paper. The method is achieved through Bayesian theory. Several sets of dispersion curves extracted from measured data are used as the input function. The inversion is performed by matching a replica calculated with a dispersion formula. The bottom characteristics can be described by the bottom reflection phase shift parameter P. The propagation range and the depth can be inverted quickly when the seabed parameters are represented by on parameter P. The inversion results improve the inversion efficiency of the seabed parameters. Consequently, the inversion efficiency and accuracy are improved while the number of inversion parameters is decreased and the computational speed of replica is increased. The inversion results have lower error than the reference values, and the dispersion curves calculated with inversion parameters are also in good agreement with extracted curves from measured data; thus, the effectiveness of the inversion method is demonstrated.