骨小梁中超声背散射信号的频率分析及其微结构的估计

骨质疏松症是一种骨强度下降的全身性骨骼疾病,骨强度的下降是骨量减少和骨微结构退化的共同结果。相比于传统的超声透射方法,超声背散射法可提供更多的骨微结构信息,而对于松质骨结构的建模能有助于结构信息的获取。本文将骨小梁简化为单圆柱模型(圆柱状的单根骨小梁浸于骨髓中),并基于此模型对超声背散射与频率的关系进行分析。用铝线代替骨小梁做仿体实验,通过实验与理论结果的比较来验证单圆柱模型的可行性。     

Computer simulation model on ultrasonic myocardial backscatter

Myocardial ultrasonic tissue characterization with integrated backscatter (IB) has been studied in recent years. To evaluate the dynamic characters of the myocardium, dynamic tracing of each unit volume of myocardium is a key point. Thus, a 2D CVIB imaging algorithm based on an auto-tracing method was developed in our previous work. Where, an auto-tracking method was used in CVIB-weighted images so that the POI in each frame throughout the cardiac cycle is automatically traced to reflect the real myocardial tissue movement. It is necessary to validate the auto-tracing method. However, at present there are no appropriate models to support the process. In order to validate the auto-tracing method, two myocardium models have been established to simulate the short and long axis view of echocardiography. The motions of partially ischemic myocardium have been simulated. The models have been used to validate the 2D CVIB imaging methods in the detection of myocardial ischemia. The simulation results show that the auto-tracing algorithm is effective. The model developed in this work provides a tool for studying and developing technologies in myocardial behavior.     

Anisotropy of the ultrasonic backscatter of myocardial tissue: I. Theory and measurements in vitro

This research addresses the variations in the ultrasonic backscatter from specimens consisting of a suspension of approximately aligned cylindrical scatterers in a fluid medium as a function of the angle of propagation in the sample. Predictions of the angular dependence of backscatter based on the time-domain Born approximation described by Rose and Richardson [J. H. Rose and J. M. Richardson, J. Nondestr. Eval. 3, 45-53 (1982)] were compared with experimental measurements of the backscatter from both tissue-mimicking phantoms consisting of graphite fibers suspended in gelatin and from canine myocardial tissue. The angular dependence of the backscatter was predicted and measured to be maximum for propagation perpendicular to the cylinder axes and minimum for propagation parallel to the axes. Maximum to minimum (i.e., perpendicular to parallel) changes in the integrated backscatter were predicted to be between 5 and 10 dB in the phantom. The corresponding quantity measured in both the phantom and in canine myocardial tissue was approximately 6 dB.     

Frequency dependence of ultrasonic backscatter from human trabecular bone: theory and experiment

A model describing the frequency dependence of backscatter coefficient from trabecular bone is presented. Scattering is assumed to originate from the surfaces of trabeculae, which are modeled as long thin cylinders with radii small compared with the ultrasonic wavelength. Experimental ultrasonic measurements at 500 kHz, 1 MHz, and 2.25 MHz from a wire target and from trabecular bone samples from human calcaneus in vitro are reported. In both cases, measurements are in good agreement with theory. For mediolateral insonification of calcaneus at low frequencies, including the typical diagnostic range (near 500 kHz), backscatter coefficient is proportional to frequency cubed. At higher frequencies, the frequency response flattens out. The data also suggest that at diagnostic frequencies, multiple scattering effects on the average are relatively small for the samples investigated. Finally, at diagnostic frequencies, the data suggest that absorption is likely to be a larger component of attenuation than scattering.     

Anisotropy of the ultrasonic backscatter of myocardial tissue: II. Measurements in vivo

The purpose of this investigation was to determine the angular dependence of the backscatter from canine myocardial tissue in vivo and to compare it with the variation of backscatter over the cardiac cycle that has been recognized and reported previously. The backscatter was measured from regions of left ventricular wall in canine hearts in which the fibers of the muscle lay parallel to the surface of the heart and were oriented predominantly in a circumferential fashion. Because of technical considerations, the angle of insonification was varied systematically through two cycles in which the angle relative to the muscle fiber axes ranged from 60 degrees-120 degrees. Backscatter was maximum at angles of interrogation perpendicular to the myocardial fibers and minimum at those most acute (60 degrees) relative to the orientation of the fibers. The previously observed variation of integrated backscatter over the heart cycle was evident at each angle of interrogation. At end systole, the average maximum-to-minimum angular variation of integrated backscatter as 5.0 +/- 0.4 dB. At end diastole, the average maximum-to-minimum angular variation was 3.2 +/- 0.4 dB. Thus, even though angular dependence of the backscatter from tissues with directionally oriented structures is substantial, the anisotropy does not account for cardiac-cycle-dependent variation of backscatter. Accordingly, the angular dependence should be incorporated in approaches to quantitative tissue characterization with ultrasound.     

Statistics of the envelope of ultrasonic backscatter from human trabecular bone

The paper describes the investigations intended to compare the results of experimental measurements of backscattering properties of the trabecular bone with the results of computer simulations. Ultrasonic RF echoes were collected using two bone scanners operating at 0.58 and 1.3 MHz. The simulations of the backscattered RF echoes were performed using the scattering model of the trabecular bone that consisted of cylindrical and spherical elements uniformly distributed in water-like medium. For each measured or simulated RF backscatter the statistical properties of the signal envelope were determined. Experimental results suggest deviations of the backscattering properties from the Rayleigh distribution. The results of simulation suggest that deviation from Rayleigh distribution depends on the variation of trabeculae diameters and the number of thin trabeculae. Experimentally determined deviations corresponded well to the deviations calculated from simulated echoes assuming trabeculae thickness variation equaled to the earlier published histomorphometric study results.     

Fundamental precision limitations for measurements of frequency dependence of backscatter: applications in tissue-mimicking phantoms and trabecular bone.

Various models for ultrasonic scattering from trabecular bone have been proposed. They may be evaluated to a certain extent by comparison with experimental measurements. In order to appreciate limitations of these comparisons, it is important to understand measurement precision. In this article, an approach proposed by Lizzi and co-workers is adapted to model precision of estimates of frequency-dependent backscatter for scattering targets (such as trabecular bone) that contain many scatterers per resolution cell. This approach predicts uncertainties in backscatter due to the random nature of the interference of echoes from individual scatterers as they are summed at the receiver. The model is validated in experiments on a soft-tissue-mimicking phantom and on 24 human calcaneus samples interrogated in vitro. It is found that while random interference effects only partially explain measured variations in the magnitude of backscatter, they are virtually entirely responsible for observed variations in the frequency dependence (exponent of a power law fit) of backscatter.     

The effect of grain size distribution on the frequency dependence of the ultrasonic attenuation in polycrystalline materials

An expression has been derived for the effect of the grain size distribution in polycrystalline materials on the frequency dependence of the ultrasonic attenuation. It has been shown that two specimens with the same mean grain size can have significantly different ultrasonic attenuations if their grain size distributions are different and that no unique solution in terms of the grain size exists for a particular ultrasonic attenuation against frequency curve. Qualitative agreement has been found between the theory and some of the experimental data available.     

Measurement of dependence of backscatter coefficient from cylinders on frequency and diameter using focused transducers--with applications in trabecular bone

A theory for the elastic scattering response from a cylinder insonified by a plane wave was previously derived by Faran. In the present paper, the empirical relationship between Faran's theory and measurements of backscatter coefficient from cylindrical targets using focused transducers is investigated. Experimental measurements of dependence of backscatter coefficient on frequency and diameter for nylon wires are reported. It is found that, under certain conditions (including weak, incoherent scattering), backscatter coefficient measurements from collections of cylindrical scatterers may be meaningfully compared with Faran's model predictions. At low frequencies, the theory and experimental measurements exhibit similar dependences on frequency and diameter, provided that the scatterers are not too densely packed. At higher frequencies, the fine structure of Faran's predictions becomes difficult to reproduce experimentally with a focused transducer. Implications regarding applications to characterization of trabecular bone are discussed.     

In vitro ultrasonic characterization of human cancellous femoral bone using transmission and backscatter measurements: relationships to bone mineral density

Thirty-eight slices of pure trabecular bone 1-cm thickness were extracted from human proximal femurs. A pair of 1-MHz central frequency transducers was used to measure quantitative ultrasound (QUS) parameters in transmission [normalized broadband ultrasound attenuation (nBUA), speed of sound (SOS)] and in backscatter [broadband ultrasound backscatter (BUB)]. Bone mineral density (BMD) was measured using clinical x-ray quantitative computed tomography. Site-matched identical region of interest (ROIs) of 7 x 7 mm2 were positioned on QUS and QCT images. This procedure resulted in 605 ROIs for all the specimens data pooled together. The short-term precision of the technique expressed in terms of CV was found to be 2.3% for nBUA, 0.3% for SOS and 4.5% for BUB. Significant linear correlation between QUS and BMD were found for all the 605 ROIs pooled, with r2 values of 0.73, 0.77, and 0.58 for nBUA, SOS, and BUB, respectively (all p < 0.05). For the BUB, the best regression was obtained with a polynomial fit of second order (r2 = 0.63). An analysis of measurements errors was developed. It showed that the residual variability of SOS is almost completely predicted by measurements errors, which is not the case for BUA and BUB, suggesting a role for micro-architecture in the determination of BUA and BUB.     

Laser ultrasonic propagation imaging method in the frequency domain based on wavelet transformation

A wavelet-transformed ultrasonic propagation imaging method capable of ultrasonic propagation imaging in the frequency domain was developed and applied as a new structural damage or flaw visualization algorithm. Since the wavelet-transformed ultrasonic propagation imaging method has strong frequency selectivity, it can visualize the propagation of ultrasonic waves of a specific frequency (for example, to isolate ultrasonic mode of interest and a damage-related ultrasonic wave). The strong frequency selectivity of the wavelet-transformed ultrasonic propagation imaging method was demonstrated, isolating only the zeroth-order asymmetrical mode of the fundamental Lamb wave modes in an anisotropic carbon fiber-reinforced plastic plate with a thickness of 5 mm. The wavelet-transformed ultrasonic propagation imaging method can also convert a complex time domain multiple wavefield into a simple frequency domain single wavefield. This feature enables easy interpretation of the results, and facilitates the precise evaluation of the location and size of structural damage or flaws. We demonstrated this capability by detecting a disbond in a sandwich structure made of Al-alloy skins and a foam core. A disbond with a diameter of 20 mm, which is representative of a common manufacturing flaw, was successfully detected, localized, and evaluated. Since a method to determine the allowable maximum pulse repetition frequency depending on target materials and structures was found by investigating the residual wave caused from the previous laser impinging, our laser ultrasonic system can scan rapidly the target with an optimal pulse repetition rate. In addition, the proposed wavelet-transformed ultrasonic propagation imaging method can visualize damage or flaw without the need for reference data from the intact state of the structure. Hence, we propose the wavelet-transformed ultrasonic propagation imaging approach for automatic inspection of in-service engineering structures, or in-process quality inspection in manufacturing.     

The dependence of time-domain speed-of-sound measurements on center frequency, bandwidth, and transit-time marker in human calcaneus in vitro

Time-domain speed-of-sound (SOS) measurements in calcaneus are effective predictors of osteoporotic fracture risk. High attenuation and dispersion in bone, however, produce severe distortion of transmitted pulses that leads to ambiguity of time-domain SOS measurements. An equation to predict the effects of system parameters (center frequency and bandwidth), algorithm parameters (pulse arrival-time marker), and bone properties (attenuation coefficient and thickness) on time-domain SOS estimates is derived for media with attenuation that varies linearly with frequency. The equation is validated using data from a bone-mimicking phantom and from 30 human calcaneus samples in vitro. The data suggest that the effects of dispersion are small compared with the effects of frequency-dependent attenuation. The equation can be used to retroactively compensate data. System-related variations in SOS are shown to decrease as the pulse-arrival-time marker is moved toward the pulse center. Therefore, compared with other time-domain measures of SOS, group velocity exhibits the minimum system dependence.     

High frequency ultrasonic characterization of human vocal fold tissue

Recently, endolaryngeal sonography at frequencies ranging from 10 to 30 MHz has been found to be useful in diagnosing diseases of the vocal folds (VFs). However, image resolution can be further improved by ultrasound at higher frequencies, necessitating the measurement of high-frequency acoustic properties of VF tissue. The ultrasonic parameters of integrated backscatter, sound velocity, and frequency-dependent attenuation coefficient were measured in both the lamina propria (LP) and vocalis muscle (VM) of human VFs using a 47 MHz high-frequency ultrasonic transducer. The integrated backscatter was -173.44+/-6.14 (mean+/-s.d.) and -195.13+/-3.58 dB in the LP and VM, respectively, the sound velocity was 1667.68+/-44.9 and 1595.07+/-39.33 ms, and the attenuation coefficient at 47 MHz was 8.28+/-1.72 and 7.17+/-1.30 dBmm. The difference between these ultrasonic parameters may be attributed to variations in the structure and fiber concentrations in VF tissue. These results could serve as a useful clinical reference for the further development of high-frequency ultrasound devices for endolarynx sonography applications.     

Dependence of the cavitation threshold on the ultrasonic frequency

A theoretical prediction of the dependence of the ultrasonic cavitation threshold for sonoluminescence on the acoustic frequency is presented. Data were obtained from the numerical solution of nonlinear oscillations of a single isolated gas-filled bubble in a viscous compressible liquid. Principal reasons for the increase of cavitation threshold with acoustic frequency and liquid viscosity are also briefly discussed.The author wishes to thank professor Ivo Hrazdira for permanent support of this work. All calculations were run on ICL 2950/10 computer of the Regional Computing Centre of the Czechoslovak Academy of Sciences, Brno.     

Stress-dependent changes in the diffuse ultrasonic backscatter coefficient in steel: experimental results

In this article, the effects of uniaxial compressive loading on the ultrasonic scattering from polycrystalline grains are shown for 10 MHz ultrasound in annealed, 1018 steel. The results show a decreasing value of the stress-dependent backscatter coefficient for normal incident ultrasound when the compression loading is perpendicular to the scattering direction. The change due to scattering is about 2 orders of magnitude greater than changes observed by others using ultrasonic wavespeed measurements. It is anticipated that this research can serve as the basis for many methods associated with nondestructive determination of stress in structural materials.     

简易超声波换能器频率测试仪的研究

对超声波换能器频率测试系统进行了探讨，并在此基础上提出了改装设想，分析了所研究的简易超声波换能器频率测试仪的电路原理及实验结果。     

Influence of ultrasonic frequency on multibubble sonoluminescence

Computer simulations of bubble oscillations are performed under conditions of multibubble sonoluminescence (MBSL) in water for various ultrasonic frequencies. The range of the ambient bubble radius for sonoluminescing bubbles narrows as the ultrasonic frequency increases; at 20 kHz it is 0.1-100 microm while at 1 MHz it is 0.1-3 microm. At 1 MHz, any sonoluminescing bubble disintegrates into a mass of smaller bubbles in a few or a few tens of acoustic cycles, while at 20 kHz and 140 kHz some sonoluminescing bubbles are shape stable. The mechanism of the light emission also depends on the ultrasonic frequency. As the ultrasonic frequency increases, the amount of water vapor trapped inside bubbles at the collapse decreases. As a result, MBSL originates mainly in plasma emissions at 1 MHz while it originates in chemiluminescence of OH radicals and plasma emissions at 20 kHz.     

Ultrasonic backscatter from flowing whole blood. II: Dependence on frequency and fibrinogen concentration

Earlier studies showed that ultrasonic backscatter from erythrocytes suspended in saline is a function of hematocrit and frequency and that it can be affected by flow disturbance. The experimental data agree well with the theories. Recently, studies have been extended to flowing whole blood. The results indicated that ultrasonic backscatter from flowing whole blood differs from that from saline suspensions of erythrocytes in that it is shear-rate dependent and species dependent. In the present article, data on the dependence of ultrasonic backscatter from flowing whole blood on frequency and on fibrinogen concentration are reported. It was found that ultrasonic backscatter from flowing whole blood also depends on fibrinogen concentration when red blood cell (RBC) aggregation exists. Moreover, when the blood is under conditions that favor RBC aggregation such as low shear rates, high fibrinogen concentration, or high hematocrits, Rayleigh scattering apparently is no longer sufficient to describe its scattering behavior.     

Origin of the bell-like dependence of the DPOAE amplitude on primary frequency ratio.

For low and medium sound pressure levels (SPLs), the amplitude of the distortion product otoacoustic emission (DPOAE) recorded from guinea pigs at the 2f1-f2 frequency is maximal when f2/f1 approximately 1.23 and decreases for lower and higher f2/f1 ratios. The high-ratio slope of the DPOAE dependence on the ratio of the primary frequencies might be anticipated since the f1 amplitude at the f2 place is expected to decrease for higher f2/f1 ratios. The low-ratio slope of the dependence at low and medium SPLs of the primaries is actually one slope of a notch. The DPOAE amplitude recovers from the notch when the f2/f1 ratio is further reduced. In two-dimensional space formed by the f2/f1 ratio, and the levels of the primaries, the notch is continuous and has a level-dependent phase transition. The notch is identical to that seen in DPOAE growth functions. Similar notches and phase transitions were observed for high-order and high-frequency DPOAEs. Theoretical analysis reveals that a single saturating nonlinearity is capable of generating similar amplitude notch and phase transition when the f2/f1 ratio is decreased because of the increase in f1 amplitude at the DPOAE generation place (f2 place). The difference between the DPOAE recorded from guinea pigs and humans is discussed in terms of different position of the operating point of the DPOAE generating nonlinearity.