A forward model and conjugate gradient inversion technique for low-frequency ultrasonic imaging

Emerging methods of hyperthermia cancer treatment require noninvasive temperature monitoring, and ultrasonic techniques show promise in this regard. Various tomographic algorithms are available that reconstruct sound speed or contrast profiles, which can be related to temperature distribution. The requirement of a high enough frequency for adequate spatial resolution and a low enough frequency for adequate tissue penetration is a difficult compromise. In this study, the feasibility of using low frequency ultrasound for imaging and temperature monitoring was investigated. The transient probing wave field had a bandwidth spanning the frequency range 2.5-320.5 kHz. The results from a forward model which computed the propagation and scattering of low-frequency acoustic pressure and velocity wave fields were used to compare three imaging methods formulated within the Born approximation, representing two main types of reconstruction. The first uses Fourier techniques to reconstruct sound-speed profiles from projection or Radon data based on optical ray theory, seen as an asymptotical limit for comparison. The second uses backpropagation and conjugate gradient inversion methods based on acoustical wave theory. The results show that the accuracy in localization was 2.5 mm or better when using low frequencies and the conjugate gradient inversion scheme, which could be used for temperature monitoring.     

适用于循环平稳声场的基于波叠加法的近场声全息技术

在原有的平面循环平稳近场声全息基础上，提出一种基于波叠加法的循环平稳近场声全息技术，可以对具有复杂表面的声源进行全息重建，重建的声源表面声压谱相关密度函数能反映出调制信号的信息.声源表面声压谱相关密度函数全息图形象地反映了调制信号在表面的强弱分布情况，可由此确定调制信号源的产生位置.仿真分析和实验验证表明，基于波叠加法的循环平稳近场声全息技术可以更准确地反映循环平稳声场的调制特性.该方法继承了波叠加法的优点，无需计算边界奇异积分，计算效率高、精度好.     

适用于循环平稳声场的基于波叠加法的近场声全息技术

在原有的平面循环平稳近场声全息基础上,提出一种基于波叠加法的循环平稳近场声全息技术,可以对具有复杂表面的声源进行全息重建,重建的声源表面声压谱相关密度函数能反映出调制信号的信息.声源表面声压谱相关密度函数全息图形象地反映了调制信号在表面的强弱分布情况,可由此确定调制信号源的产生位置.仿真分析和实验验证表明,基于波叠加法的循环平稳近场声全息技术可以更准确地反映循环平稳声场的调制特性.该方法继承了波叠加法的优点,无需计算边界奇异积分,计算效率高、精度好. 关键词： 近场声全息 循环平稳信号 波叠加     

An effective quiescent medium for sound propagating through an inhomogeneous,moving fluid

The idea of similarity between acoustic fields in a moving fluid and in a certain "effective" quiescent medium, first put forward by Lord Rayleigh, proved very helpful in understanding and modeling sound propagation in an atmosphere with winds and in an ocean with currents, as well as in other applications involving flows with small velocity compared to sound speed. Known as effective sound speed approximation, the idea is routinely utilized in the contexts of the ray theory, normal mode representation of the sound field, and the parabolic approximation. Despite the wide use of the concept of effective sound speed in acoustics of moving media, no theoretical justification of Rayleigh's idea was published that would be independent of the chosen representation of the sound field and uniformly apply to distinct propagation regimes. In this paper, we present such a justification by reducing boundary conditions and a wave equation governing sound fields in the inhomogeneous moving fluid with a slow flow to boundary conditions and a wave equation in a quiescent fluid with effective sound speed and density. The derivation provides insight into validity conditions of the concept of effective quiescent fluid. Introduction of effective density in conjunction with effective sound speed is essential to ensure accurate reproduction of acoustic pressure amplitude in the effective medium. Effective parameters depend on sound speed, flow velocity, and density of the moving fluid as well as on sound propagation direction. Conditions are discussed under which the dependence on the propagation direction can be avoided or relaxed.     

三维复杂温度场高精度声学测量方法

声学温度场检测技术通过多路径声波传播时间数据,反演被测区域的温度分布。提供了一种高精度的三维复杂温度场的声学测量方法。首先从射线声学角度给出了三维非均匀温度场中声波传播路径的数学模型。在此基础上,将三维温度场的重建问题转化为声波传播路径的求解和温度场的反演问题,建立了基于多项式修正径向基函数(RBF-PR)和改进的Tikhonov正则化三维温度场重建算法。采用两种典型的炉膛三维温度场模型,在信噪比SNR=35 dB下进行了数值模拟,分析了声波传播路径在非均匀温度场中的弯曲特性、算法的重建质量和抗噪性,同时进行了实际炉膛内二维温度场的重建。结果表明了提出的考虑声线弯曲的温度场重建算法具有精度高,抗噪性强、适用性好的特点,为实现高精度的复杂温度场的声学测量提供了有效方法。     

Multi-step geoacoustic inversion with a towed line array

Geoacoustic inversion is important for acoustic field predictions and matched field localizations in shallow water.Combing the Matched-field inversion(MFI) and the Reflectionloss inversion,a multi-step Bayesian inversion for geoacoustic parameters was presented.This method applied the posterior probability density(PPD) or inversion results from one inversion as prior information for subsequent inversion.First,the sensitive parameters were determined by the MFI.Second,the insensitive parameters were determined by the Reflection-loss inversion based on the PPD and inversion results from the MFI.The PPD results indicated that the multi-step inversion method was performed better than direct matched-field inversion,and the inversion results of some parameters were improved significantly.To demonstrate the advantages of the multi-step inversion method on the sound field prediction,the statistical properties of transmission loss based on the posterior probability were introduced.The transmission loss distribution showed that the predicted acoustic fields based on the multi-step inversion method had smaller errors.     

舰船辐射噪声贝叶斯方法反演浅海底层结构及地声参数EI北大核心CSCD

针对以舰船辐射噪声为参考声源的浅海海底分层结构及地声参数反演问题,研究了一种基于贝叶斯理论的浅海多层海底地声参数反演方法。反演中以舰船辐射噪声的线谱成分为研究对象,进而采用非线性贝叶斯反演方法反演浅海底层结构、层中声速、声速衰减和密度,并对反演结果的不确定性进行分析。反演结果的最大后验概率估计值和边缘概率分布分别通过拨正模拟退火算法和Metropolis-Hastings采样法在各参数先验区间内计算获得,并根据贝叶斯信息准则确定最佳海底分层结构。海上实验表明:根据该方法反演获得海底分层结构及地声参数,计算得到的声压场与实测舰船辐射噪声传播损失误差不超过10%,反演结果能够准确表征实验海区海底特征。反演结果不确定性分析表明:海底纵波声速、横波声速以及密度的不确定性更小,对声压场变化更加敏感,反演结果更有效、准确。     

Full wave-field reflection coefficient inversion

This paper develops a Bayesian inversion for recovering multilayer geoacoustic (velocity, density, attenuation) profiles from a full wave-field (spherical-wave) seabed reflection response. The reflection data originate from acoustic time series windowed for a single bottom interaction, which are processed to yield reflection coefficient data as a function of frequency and angle. Replica data for inversion are computed using a wave number-integration model to calculate the full complex acoustic pressure field, which is processed to produce a commensurate seabed response function. To address the high computational cost of calculating short range acoustic fields, the inversion algorithms are parallelized and frequency averaging is replaced by range averaging in the forward model. The posterior probability density is interpreted in terms of optimal parameter estimates, marginal distributions, and credibility intervals. Inversion results for the full wave-field seabed response are compared to those obtained using plane-wave reflection coefficients. A realistic synthetic study indicates that the plane-wave assumption can fail, producing erroneous results with misleading uncertainty bounds, whereas excellent results are obtained with the full-wave reflection inversion.     

气泡体积分数对沙质沉积物低频声学特性的影响

由于有机物质分解等原因,实际的海底沉积物中存在气泡,气泡的存在会显著影响沉积物低频段的声学特性,因此研究气泡对沉积物低频段声速的影响机理具有重要意义.考虑到外场环境的不可控性,在室内水池中搭建了大尺度含气非饱和沙质沉积物声学特性获取平台,在有界空间中应用多水听器反演方法首次获取了含气非饱和沙质沉积物300—3000 Hz频段内的声速数据(79—142 m/s),并同时利用双水听器法获取了同一频段的数据(112—121 m/s).在声波频率远低于沉积物中最大气泡的共振频率时,根据等效介质理论,将孔隙水和气泡等效为一种均匀流体,改进了水饱和等效密度流体近似模型.模型揭示了气泡对沉积物低频段声学特性的影响规律,理论上解释了沉积物中声速的降低.通过分析模型预报声速对模型参数的敏感性,根据测量得到的声速分布反演得到了沉积物不同区域的气泡体积分数,气泡体积分数从1.07%变化到2.81%.改进的模型为沉积物中气泡体积分数估计提供了一种新方法.     

谐振管谐振频率计算方法的研究

为了准确地计算声谐振管的谐振频率,提出了基于对管内声压分布进行理论模拟的声压模拟法,以及基于驻波最小点位置与管两端声阻抗关系的驻波最小点法。结果表明这两种方法因为考虑了声谐振管两端的声阻对谐振频率的影响,所得谐振频率更加精确、严格,所以特别适用于要求谐振管的谐振频率严格控制的声学领域。同时表明:声压模拟法与驻波最小点法虽然出发点不同,但谐振频率的计算结果完全相同。驻波最小点法计算比较方便,而声压模拟法在计算出谐振频率的同时,还可以给出管中的声压随位置的分布情况。     

A wave field synthesis approach to reproduction of spatially correlated sound fields

This article discusses an open-loop wave field synthesis (WFS) approach for the reproduction of spatially correlated sound fields. The main application concerns laboratory reproduction of turbulent boundary layer wall pressure on aircraft fuselages and measurement of their sound transmission loss. The problem configuration involves reconstruction of random sound pressure distributions on a planar reproduction surface using a planar array of reproduction monopoles parallel to the reproduction plane. In this paper, the WFS formulation is extended to sound fields with imposed time and spatial correlation properties (or equivalently imposed cross-spectral density in the frequency and wave number domains). Numerical examples are presented for the reproduction of a propagating plane wave, diffuse acoustic field and wall pressure in subsonic or supersonic turbulent boundary layers. The reproduction accuracy is examined in terms of the size of the source plane and reproduction plane, their separation, and the number of reproduction sources required per acoustic wavelength. While the reproduction approach cannot reconstruct sub-wavelength correlation scales of subsonic turbulent boundary layers, it effectively reconstructs correlation scales larger than the acoustic wavelength, making it appropriate for diffuse acoustic field and supersonic turbulent layers.     

Adjoint modeling for acoustic inversion

The use of adjoint modeling for acoustic inversion is investigated. An adjoint model is derived from a linearized forward propagation model to propagate data-model misfit at the observation points back through the medium to the medium perturbations not being accounted for in the model. This adjoint model can be used to aid in inverting for these unaccounted medium perturbations. Adjoint methods are being applied to a variety of inversion problems, but have not drawn much attention from the underwater acoustic community. This paper presents an application of adjoint methods to acoustic inversion. Inversions are demonstrated in simulation for both range-independent and range-dependent sound speed profiles using the adjoint of a parabolic equation model. Sensitivity and error analyses are discussed showing how the adjoint model enables calculations to be performed in the space of observations, rather than the often much larger space of model parameters. Using an adjoint model enables directions of steepest descent in the model parameters (what we invert for) to be calculated using far fewer modeling runs than if a forward model only were used.     

匹配垂直质点振速的浅海地声参数反演

海底声学参数对海洋波导中的声场特性研究和相关应用具有重要意义。针对一次夏季黄海声传播实验,分析了浅海负跃层环境下垂直质点振速的传播特性和简正波结构,说明当声源和接收器均位于负跃层下时,除海底附近外的大部分深度上垂直质点振速能量较高,且与声压相比,号数高的简正波对垂直质点振速的贡献更大,利用垂直质点振速进行匹配场反演能获得更高的海底参数敏感性。分析了海底吸收系数对垂直质点振速匹配场反演的影响,结果表明只有当进行匹配场反演时设置的海底吸收系数接近真实值时,才能获得准确的海底声速、密度和海深反演结果。利用实验中矢量水听器获取的垂直质点振速信号进行匹配场反演,将海底吸收系数在变化范围内取不同值对海底声速、密度和平均海深进行全局搜索,根据代价函数值最大确定了海底声速、密度及平均海深的反演结果,并利用不同距离上的声压传播损失反演出不同频率下的海底吸收系数。根据反演得到的海底声学参数计算声压传播损失,与实验中声压水听器测量结果符合较好。     

The effects of ignored seabed variability in geoacoustic inversion

Geoacoustic inversion using a matched-field inversion algorithm is a well-established technique for estimating the geoacoustic parameters of the seabed. This paper demonstrates how parameter estimation can be affected by unknown or wishfully ignored random range dependence of the true environment when the inversion model is--for practicality--assumed to be range independent. Simulations with controlled statistics were carried out using a simple shallow water model: an isospeed water column over a homogeneous elastic halfspace. The inversion parameters included water depth, compressional speed in the seabed, seabed density, and compressional wave attenuation. On average the environment is range independent: some parameters are constant while other parameters are random with range-independent means and variances. A Parabolic Equation underwater acoustic propagation model is used to calculate the simulated data fields for the range-dependent environment as well as to calculate the model fields for the range-independent inversion model. The Adaptive Simplex Simulated Annealing inversion algorithm is used to estimate the best-fit solution. It is found that ignoring the variability of even a single geoacoustic parameter leads to significant and correlated uncertainty (bias and variance) in the estimation of all inverted parameters. Results are presented for range variation of compressional sound speed and water depth.     

Acoustic waves in random fields

The effect of space- and time-dependent random mass density, velocity, and pressure fields on frequencies and amplitudes of acoustic waves is considered by means of the analytical perturbative method. The analytical results, which are valid for weak fluctuations and long wavelength sound waves, reveal frequency and amplitude alteration, the effect of which depends on the type of random field. In particular, the effect of a random mass density field is to increase wave frequencies. Space-dependent random velocity and pressure fields reduce wave frequencies. While space-dependent random fields attenuate wave amplitudes, their time-dependent counterparts lead to wave amplification. In another example, sound waves that are trapped in the vertical direction but are free to propagate horizontally are affected by a space-dependent random mass density field. This effect depends on the direction along which the field is varying. A random field, which varies along the horizontal direction, does not couple vertically standing modes but increases their frequencies and attenuates amplitudes. These modes are coupled by a random field which depends on the vertical coordinate, but the dispersion relation remains the same as in the case of the deterministic medium.     

Geoacoustic inversion for bottom parameters in a thermocline environment in the northern area of the South China Sea

Bottom acoustic parameters have important influence on the application of underwater acoustic propagation and source location.The acoustic parameters of the seabed in the northern of the South China Sea(SCS) were inversed using the experiment data from an acoustic experiment in 2015.Based on the comprehensive analysis of the influence of the sound speed fluctuation and the geoacoustic model on seabed inversion,the multi-parameter hybrid acoustic inversion scheme is improved by selecting the equi...     

偶极随钻声波测井声压和径向位移响应的差异*

声波测井仪接收到的电信号通常是多个压电片响应的叠加,它主要是由声压还是径向位移响应转化而来,或是两种响应兼有目前未有定论。该文通过实轴积分法和复变函数法计算并对比分析随钻声波测井的声压和径向位移场,发现这两种响应特性有着显著的差异。首先,软地层的偶极随钻测井时,声压信号包含钻铤波和舒尔特波两个波群,而径向位移信号仅有钻铤波波群;其次,单极声源情况下,声压和径向位移信号的钻铤波能量分别集中在钻铤内、外壁,而偶极情况恰好相反,可见,钻铤按照单极情况的分析结果进行刻槽后,高频时的拖尾现象会影响偶极信号中舒尔特波对横波速度的反演。因此,阐明两类信号的差异对横波速度的反演和钻铤波的压制都具有重要意义。     

空间声场全息重建的波叠加方法研究

提出了基于波叠加法的近场声场全息技术，并将其用于任意形状物体的声辐射分析.在声辐射计算问题中，边界元法是通过离散边界面上的声学和位置变量来实现，而波叠加方法则通过叠加辐射体内部若干个简单源产生的声场来完成.因而，基于波叠加法的声全息就不存在边界面上的参数插值和奇异积分等问题，而这些问题是基于边界元法的声全息所固有的.与基于边界元法的声全息相比较，基于波叠加法的声全息在原理上更易于理解，在计算机上更容易实现.实验结果表明：该种全息技术在重建声场时，具有令人满意的重建精度. 关键词： 声全息 逆问题 波叠加方法 正则化方法     

Sound-absorbing media with two types of acoustic losses

Sound-absorbing media with two types of acoustic losses, namely, inertial and deformation are considered. Such a medium is described by a complex density and a complex compressibility. It is shown that a certain relation between the latter quantities imparts some specific properties to the medium: a real wave impedance, complete absorption of sound waves at the boundary with the external medium, the properties of a viscous or “superviscous” medium, etc. The results are generalized to the case of a layered inhomogeneous medium with two types of losses.     

Direct analysis of dispersive wave fields from near-field pressure measurements

Flexural waves play a significant role for the radiation of sound from plates. The analysis of flexural wave fields enables the detection of sources and transmission paths in plate-like structures. The measurement of these wave fields can be carried out indirectly by means of near-field acoustic holography, which determines the vibrational wave field from pressure information measured in a plane close to the plate under investigation. The reconstruction of the plate vibration is usually obtained by inverting the forward radiation problem, i.e., by inversion of an integral operator. In this article, it is shown that a pressure measurement taken in the extreme near-field of a vibrating plate can directly be used for the approximate analysis of the dispersive flexural wave field. The inversion step of near-field acoustic holography is not necessarily required if such an approximate solution is sufficient. The proposed method enables fast and simple analysis of dispersion characteristics. Application of dispersion compensation to the measured field allows for visualizations of propagating wavefronts, such that sources and scatterers in the plate can be detected. The capabilities of the described approach are demonstrated on several measurements.