超声治癌与无损测温研究

至今已有很多关于超声治癌的研究。热疗已上升为治疗癌肿的主要方法之一。文中简述了热疗的历史发展，生物物理学基础和超声热疗的技术面貌，报道了对肿瘤或前列腺烧疗的高强度相控聚焦超声技术，它可在１０ｓ以内使肿瘤内升温超过７０℃。     

激光超声无损检测技术

 传统的无损检测方法具有一定的局限性,例如:在高温高压、高湿、有毒等某些恶劣检测环境中,或被测工件具有放射性或腐蚀性,以及被测工件具有较快的运动速度时,已不能完全满足要求。因此,人们一直都在努力寻找适合于上述情况的无损检测方法。这里介绍一种激光超声检测技术,它是利用激光脉冲激发超声从而实现无损检测的一种方法,已逐渐成为材料无损检测的一种重要手段。本文详细介绍了激光超声的检测原理及应用,并指出目前尚存在的问题及解决措施。激光超声的激发原理激光超声是指用脉冲激光在介质中所产生的超声波或利用激光来产生超声这一物理过程。激光可以在固体中产生超声,也可以在气体和液体中产生超声。     

超声红外热像技术及其在无损评价中的应用

本介绍在光热红外检测的基础上发展的超声红外热像技术。该技术利用超声脉冲作为激发源，当超声脉冲在试样上传播的过程中遇到裂纹等缺陷时，缺陷引起超声附加衰减而局部升温。利用红外照相机获取试样表面的温度分布，可显示裂纹等缺陷。超声红外热像技术发挥了超声和红外技术的优点，可实时地检测裂纹等缺陷，在无损评价和检测中有广泛应用。     

三次谐波超声脉冲回波法对厚板的无损检测

在同一厚板表面的两个不同位置处分别进行超声回波脉冲测量，一处有细微裂纹而另一处无细微裂纹。实验结果表明，即使两个位置回波的基波频谱非常接近，但相应的三次谐波回波频谱仍存在较显著的区别，这一结果可望用于厚板具有微缺陷时的超声无损检测。     

激光超声无损检测技术及其在岩体物性测试中的应用

激光超声无损检测技术是材料无损检测的一种重要手段,是目前人们所关注的一种重要的检测方法.文章阐述了激光超声的激发原理和激光超声信号检测方法,详细介绍其在对岩体物性测试中的应用,并且针对目前存在的一些技术难题指出了今后的解决方向.     

用于铝板缺陷无损检测的激光超声有限元模拟研究

激光超声技术能够在材料表面形成超声波，是实现材料缺陷无损检测的重要环节。借助Abaqus有限元分析工具，基于激光超声热弹机制建立了轴对称铝板的表面缺陷模型，模拟了激光激发产生的表面波在材料中的传播特性及其与铝板缺陷的相互作用过程。数值模拟实验表明，铝板表面缺陷的分布深度值越大，反射波越强，并且当缺陷深度达到一定程度时，反射波的幅值趋于稳定；但缺陷的分布宽度对于反射波的影响则十分有限。所得结论为基于激光超声的材料缺陷的定量检测及识别提供有效的理论依据。     

基于遗传算法的超声信号LMS自适应时延估计

为了克服ＬＭＳ自适应时延估计（ＬＭＳＴＤＥ）算法计算量大的问题,引进遗传算法进行ＬＭＳＴＤＥ的寻优规划,并采取了克服过早收敛的措施。对超声信号进行时延估计的实验表明,该方法大大减少了计算量,并有较高的时延估计精度。     

Estimation of speed of sound in dual-layered media using medical ultrasound image deconvolution

The speed of sound in soft tissues is usually assumed to be 1540 m/s in medical pulse-echo ultrasound imaging systems. When the true speed is different, the mismatch can lead to distortions in the acquired images, and so reduce their clinical value. Previously we reported a new method of sound-speed estimation in the context of image deconvolution. Unlike most other sound-speed estimation methods, this enables the use of unmodified ultrasound machines and a normal scanning pattern. Our approach was validated for largely homogeneous media with single sound speeds. In this article, we demonstrate that sound speeds of dual-layered media can also be estimated through image deconvolution. An ultrasound simulator has been developed for layered media assuming that, for moderate speed differences, the reflection at the interface may be neglected. We have applied our dual-layer algorithm to simulations and in vitro phantoms. The speed of the top layer is estimated by our aforesaid method for homogeneous media. Then, when the layer boundary position is known, a series of deconvolutions are carried out with dual-layered point-spread functions having different lower-layer speeds. The best restoration is selected using a correlation metric. The error level (e.g., a mean error of −9 m/s with a standard deviation of 16 m/s) for in vitro phantoms is found to be not as good as that of our single-speed algorithm, but is comparable to other local speed estimation methods where the data acquisition may not be as simple as in our proposed method.     

A new approach for the estimation of the axial velocity using ultrasound

The most used estimation method for calculating the blood velocity in commercial scanners is the autocorrelation approach. The calculation of the mean velocity used in this method depends on the center frequency of the interacting ultrasound pulse which downshifts as a function of depth, introducing a bias. A new velocity estimator for the mean axial velocity is presented. The estimation principle is based on the 2D Fourier transform and the Radon transform. The input data are a sequence of RF data forming a 2D data input, one column for each pulse emission. A 2D segment is selected for a specific depth. This data segment is first transformed by a 2D Fourier transform, and the result is then transformed by a Radon transform. The center of gravity for the angles of the lines intersecting the origin of the R-theta coordinate system in the Radon domain gives the mean axial velocity for the data segment. The benefit of this method is an estimate of the mean axial velocity which is independent of the center frequency of the propagating ultrasound pulse. The estimate will only depend on f(s) and f(pr f). Results of the estimation method is presented based on both simple generated RF harmonic data for different signal/noise ratios and simulated acoustic RF responses from a 3D measurement situation with an array transducer and a tube with plug flow. The new method shows improvement with a factor of 1.5-4 on the standard deviation on the estimated mean velocity for the simulated case.     

基于最优估计理论、联合星载主被动传感器资料的液态云微物理特性反演研究

针对液态云微物理特性精确反演的迫切需求, 综合主被动传感器的探测优势, 联合CloudSat雷达反射率和Aqua光学厚度资料, 提出基于最优估计理论的液态云微物理参数反演算法.通过假设粒子谱服从对数正态分布, 基于前向物理模式建立测量变量与反演变量的函数关系, 借助谱分布参数的先验信息、通过算法迭代得到谱参数的最优解, 进而利用前向物理模式反演液态云微物理参数, 并根据误差传递理论计算反演不确定度.通过设计反演方案, 基于实测个例数据并与CloudSat官方发布产品和经验算法反演结果对比验证.结果表明: 基于最优估计理论、联合主被动传感器资料的液态云微物理参数反演结果与官方发布产品一致性较好, 弥补了经验算法误差大、扩展性差的不足, 对于开展国内星载和机载W波段毫米波雷达液态云微物理参数反演研究具有重要的借鉴意义. 关键词： CloudSat Aqua 液态云 最优估计理论     

Simultaneous determination of the ultrasound velocity and the thickness of solid plates from the analysis of thickness resonances using air-coupled ultrasound

A method that combines transmission of air-coupled ultrasound pulses through solid plates and amplitude and phase spectral analysis is presented. In particular, the method analyzes the first thickness resonance of the plates. The purpose is to determine, simultaneously, velocity and attenuation coefficient of the ultrasounds in the material and the thickness of the plate. This is especially useful when thickness can not be measured independently. The method is successfully applied to soft membranes, biological samples and FRP composites.     

Optimal estimation of the diffusion coefficient from non-averaged and averaged noisy magnitude data

The magnitude operation changes the signal distribution in MRI images from Gaussian to Rician. This introduces a bias that must be taken into account when estimating the apparent diffusion coefficient. Several estimators are known in the literature. In the present paper, two novel schemes are proposed. Both are based on simple least squares fitting of the measured signal, either to the median (MD) or to the maximum probability (MP) value of the Probability Density Function (PDF). Fitting to the mean (MN) or a high signal-to-noise ratio approximation to the mean (HS) is also possible. Special attention is paid to the case of averaged magnitude images. The PDF, which cannot be expressed in closed form, is analyzed numerically. A scheme for performing maximum likelihood (ML) estimation from averaged magnitude images is proposed. The performance of several estimators is evaluated by Monte Carlo (MC) simulations. We focus on typical clinical situations, where the number of acquisitions is limited. For non-averaged data the optimal choice is found to be MP or HS, whereas uncorrected schemes and the power image (PI) method should be avoided. For averaged data MD and ML perform equally well, whereas uncorrected schemes and HS are inadequate. MD provides easier implementation and higher computational efficiency than ML. Unbiased estimation of the diffusion coefficient allows high resolution diffusion tensor imaging (DTI) and may therefore help solving the problem of crossing fibers encountered in white matter tractography.     

Investigation of transmit and receive performance at the fundamental and third harmonic resonance frequency of a medical ultrasound transducer

In this study, the phenomenon of higher harmonic thickness resonance of a piezoelectric transducer was used to investigate potentially additional sensitivity at the third harmonic frequency for conventional medical transducers. The motivation for this research is that some applications in medical ultrasound (e.g. third harmonic transmit phasing and contrast imaging) need probes which are sensitive around both the fundamental and third harmonic frequencies, and that these higher harmonic thickness modes, although often considered as undesired, might be used beneficially. The novelty aspect in this study is the presented transmit and receive potential at both the fundamental and third harmonic of a conventional cardiac probe with modified electrical tuning. Elements of an experimental PZT-based phased-array probe (f_c = 3 MHz, 64 elements, element width = 0.3 mm, elevation aperture = 13 mm) were electrically retuned with series inductors around the third harmonic resonance frequency at 10 MHz. Hydrophone measurements with 10-MHz-tuned elements showed that, as compared to a conventionally tuned element, the transmit transfer function at the third harmonic increased more than 23 dB, while the sensitivity at the fundamental frequency was only 6 dB lower. Pulse-echo measurements showed that the two-way transfer function of a 10-MHz-tuned element resulted in 20 dB increased sensitivity around the third harmonic as compared to an untuned element. Simulated transfer functions, from both a 1D KLM and 2D finite element model of an element of the experimental array transducer, confirmed the measured sensitivity peaks at the fundamental and third harmonic. In conclusion, this study demonstrated the effect of changing the electrical tuning on a conventional array transducer which increased the sensitivity around the third harmonic resonance frequency, while maintaining good sensitivity at the fundamental frequency.     

The estimation of lower refractivity uncertainty from radar sea clutter using the Bayesian-MCMC method

The estimation of lower atmospheric refractivity from radar sea clutter(RFC) is a complicated nonlinear optimization problem.This paper deals with the RFC problem in a Bayesian framework.It uses the unbiased Markov Chain Monte Carlo(MCMC) sampling technique,which can provide accurate posterior probability distributions of the estimated refractivity parameters by using an electromagnetic split-step fast Fourier transform terrain parabolic equation propagation model within a Bayesian inversion framework.In contrast to the global optimization algorithm,the Bayesian-MCMC can obtain not only the approximate solutions,but also the probability distributions of the solutions,that is,uncertainty analyses of solutions.The Bayesian-MCMC algorithm is implemented on the simulation radar sea-clutter data and the real radar seaclutter data.Reference data are assumed to be simulation data and refractivity profiles are obtained using a helicopter.The inversion algorithm is assessed(i) by comparing the estimated refractivity profiles from the assumed simulation and the helicopter sounding data;(ii) the one-dimensional(1D) and two-dimensional(2D) posterior probability distribution of solutions.     

从相关测量的数据恢复超短激光脉冲的波形

本文证明了,只要测量了脉冲光强的二阶相关函数G~(2)(τ)和单延迟的三阶相关函数G~(3)(τ),便可以恢复超短激光脉冲的强度波形。并提出一种恢复脉冲波形的计算方法,编制了计算机程序,对实验数据作了处理,证实这种计算方法是可行的。在此基础上建议一种通过相关测量确定超短激光脉冲波形的新型仪器。     

Investigation of low energy pion-nucleus interaction using electrofission data for heavy nuclei at the pion threshold

A theoretical approach to describe electro- and photofission processes at intermediate energies was worked out. Photopion reabsorption mechanisms by two and three nucleons (2NA and 3NA) were incorporated in the calculations. The comparison with electrofission data for preactinides showed that a substantial 3NA component should be added in the pion absorption to fit the data near photopion threshold, in accordance with a recent theoretical estimate. It was also shown that the shape of the fission response curve is sensitive to detailed aspects of the pion mean free path. Received: 10 March 1998 / Revised version: 12 March 1999     

An exact inversion method for the determination of spin–orbit potentials from scattering data

A generalized Newton–Sabatier inversion method which permits extraction from scattering data of central and spin-orbit potentials is presented. The inversion method originally developed by Sabatier and further elaborated by Hooshyar and Richardson, has been reformulated to lead to physically reasonable solutions and to allow for its numerical implementation. Numerical problems due to the occurrence of singularities in the transformation kernel are discussed and a successful application using schematic scattering data is reported. Received: 30 April 1997 / Revised version: 20 October 1997