CN106840362A - Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoelastic method - Google Patents

Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoelastic method Download PDF

Info

Publication number
CN106840362A
CN106840362A CN201710166704.XA CN201710166704A CN106840362A CN 106840362 A CN106840362 A CN 106840362A CN 201710166704 A CN201710166704 A CN 201710166704A CN 106840362 A CN106840362 A CN 106840362A
Authority
CN
China
Prior art keywords
vibration measurement
hifu
point
vibration
laser
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710166704.XA
Other languages
Chinese (zh)
Other versions
CN106840362B (en
Inventor
陆明珠
韩丹
张灵璐
马风超
李玉娇
刘妍杉
关宇波
万明习
王晓东
骆志坚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201710166704.XA priority Critical patent/CN106840362B/en
Publication of CN106840362A publication Critical patent/CN106840362A/en
Application granted granted Critical
Publication of CN106840362B publication Critical patent/CN106840362B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H9/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means

Abstract

The present invention discloses a kind of vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response and damages viscoelastic method, comprises the following steps:1) control HIFU transducers produce driving pulse to act on the excitation point inside tested tissue;Vibration of the vibration measurement with laser equipment to tested tissue inside test point is tracked;Record the oscillating curve of each test point;2) according to oscillating curve, shear wave velocity and frequency are calculated, then carries out viscoplasticity fitting.The present invention uses pulse to enter row energization to hyaline tissue inside to produce shearing wave motion, the vibration put in shearing direction of wave travel is tracked using vibration measurement with laser equipment, viscoplasticity estimation is carried out by analyzing the movable information i.e. shear wave velocity and frequency of each mark point, and then is realized the monitor in real time damaged to hyaline tissue inside HIFU and is evaluated.

Description

Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoelastic Method
Technical field
The invention belongs to biomedical measurement technical field, more particularly to a kind of monitoring high intensity focused ultrasound (HIFU) is damaged Hinder viscoelastic method.
Background technology
High intensity focused ultrasound (High Intensity Focused Ultrasound, HIFU), is a kind of noninvasive, peace Complete and efficient treatment method, it focuses on high energy ultrasound beam by the HIFU transducers placed in vitro, and in focal zone Biological modification is produced to realize the treatment of disease in domain.In order to effectively guide treatment focal regions, the control ultrasonic therapy of HIFU Dosage, guarantee HIFU can be applied to clinic, it is necessary to carry out monitor in real time and evaluation to the therapeutic process of HIFU effective and safe. During HIFU Treatment, due to the dissipation effect organized so that supersonic beam kinetic energy in communication process is converted into heat energy, positioned at Jiao The histocyte in area is destroyed, protein denaturation so as to form irreversible coagulation necrosis, cause tissue hardness change from And change its mechanical characteristic.Therefore, it can by the mechanical characteristic to damage field that especially viscoelastic property is supervised in real time Control to evaluate therapeutic effect.
Shearing scattering of wave ultrasonic vibration measurement (Shear wave Dispersion Ultrasound Vibrometry, SDUV) method is a kind of method based on acoustic radiation force shearing wave, can both separate viscosity and elasticity, them can be obtained again each From quantized value, be a kind of very effective detection method, in having been applied in vitro and biological tissue.
Local organization is activated to steady-state response by existing SDUV methods using the continuous wave of amplitude modulation(PAM), to produce harmonic wave Shearing wave motion, then measures the shear wave velocity under different frequency respectively, then by being calculated the shear viscosity and bullet of tissue Property modulus.The patent of invention US8659975B2 of James F.Greenleaf et al., patent name " Vibration generation and detection in shear wave dispersion ultrasound vibrometry with In large background motions ", the system and method that a kind of coherent imaging is applied to SDUV is given, the method is adopted With a ultrasonic transducer respectively to organizing into row energization and motion detection.By the continuous wave excitation tissue for applying certain frequency Harmonic wave shearing wave motion is produced, then two harmonic wave shearing waves of diverse location on transonic direction are obtained using coherent plane wave Motor message, just can carry out viscoplasticity and estimate by harmonic wave shearing wave the movable information such as amplitude and phase etc. of analyzing different frequency Meter.In the method by changing Kalman filtering parameter come while determine multiple harmonic motion frequencies, and will during filtering Background motion is modeled as random process to compensate the motion of the overall background in SDUV.In the patent of Shigao Chen et al. US20110263978A1, patent name " Method for shear wave ultrasound vibrometry with It is also to organizing into row energization using continuous wave in interleaved push and detection pulses ".Unlike By pumping signal transmitting alternate with detection signal in the method, such that it is able in obtaining larger vibration displacement and detection process Secure threshold of the ultrasonic radiation dosage without departing from FDA.However, existing SDUV methods need to use continuous wave by tissue excitation To steady-state response, this process needs to consume the plenty of time, is unfavorable for using in real time.
Shearing wave is monitored using ultrasonic imaging generally in SDUV, but there is certain limitation in the method for ultrasonic imaging Property, such as image resolution ratio is low, echo easily by noise and organism displacement interference and can be produced with HIFU sound fields mutual Effect etc. so that the method is unfavorable for monitor in real time and detection HIFU Treatment process.Vibration measurement with laser method can reach sub-micron Precision, and can apply to the measurement in tissue surface and hyaline tissue, because laser signal does not influence ultrasonic sound field, because This can be used to the various vibrations during monitor in real time HIFU Treatment.The patent of invention of James F.Greenleaf et al. US9044192B2, patent name " System and method for non-invasively measuring tissue System and the side using surface wave measurement tissue viscoelasticity are provided in viscoelasticity using surfacewaves " Method, produces vibration, then using laser vibration measurer to group using transducer or mechanical vibrator in tissue surface in the method The surface wave for knitting surface is tracked so as to estimate its viscoplasticity information.However, existing vibration measurement with laser method is confined to tissue The detection of surface vibration, it is impossible to effectively organization internal vibration is tracked and is detected.
The content of the invention
Damage viscous it is an object of the invention to provide a kind of vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response The method of elasticity, acts on tissue, and vibration measurement with laser technology is applied to the prison of hyaline tissue's internal shear ripple using pulse In survey, to overcome the shortcomings of to act on tissue using continuous wave in existing SDUV, and be applied to for vibration measurement with laser technology by realization The breakthrough of bright organization internal.
To achieve the above object, the present invention is adopted the following technical scheme that:
Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoelastic method, comprises the following steps:
1) HIFU transducers produce driving pulse to act on the excitation point inside tested tissue;Vibration measurement with laser equipment is to tested The vibration of the test point of organization internal is tracked;Record the oscillating curve of each test point;
2) according to oscillating curve, shear wave velocity and frequency are calculated, then carries out viscoplasticity fitting.
Further, step 1) in, regulation vibration measurement with laser equipment is confocal with HIFU transducers, makes vibration measurement with laser device signal The direction of the launch is in the same direction with the HIFU transducer ultrasonic direction of the launch;Synchronizing signal control system produces synchronizing signal to drive respectively arbitrarily Waveform generator and data collecting card, the pulse train that AWG is produced is by power amplifier and impedance matching net Network acts on HIFU transducers, and control HIFU transducers produce driving pulse to act on the excitation point inside tested tissue;Data Capture card is acquired to the signal that vibration measurement with laser equipment is exported.
Further, the tested tissue is transparent prosthese or transparent in vitro tissue.
Further, vibration measurement with laser equipment is laser doppler vibrometer.
Further, step 1) specifically include following steps:
Step S1, by HIFU transducers focus at a certain position in tested tissue inside, the point be called to encourage point, HIFU to change Energy device produces pulse to act on point induction and produces shearing wave signal;
Step S2, using vibration measurement with laser equipment to excitation point and shearing direction of wave travel up to lack two vibrations of test point It is monitored and records its oscillating curve.
Further, calculate shear wave velocity and the method for frequency is specially:
Mark point position is fixed and spacing with the time delay of same-phase between 2 points, it is known that represent the shearing wave between 2 points Speed, 0 is designated as by excitation point in the time of a certain phase, and test point is in the time that the time of same distance is exactly that shearing wave is propagated; The distance between note excitation point and test point are Δ d, and the time interval of 2 points of arrival same-phases is Δ T1, then shearing wave is in detection Point speed be:
Time difference in shearing wave oscillating curve between minimum and maximum 2 points of relaxation amplitude is considered π/ω, therefore Its frequency can be calculated;Remember that this 2 time differences are Δ T2, then shearing wave vibration frequency be
In the case where local organization shear wave velocity and frequency distribution is drawn, Voigt models according to biological tissue and its Derivation formula carries out viscoplasticity fitting;
Method for carrying out viscoplasticity fitting is:
Relation between plane shear velocity of wave and tissue viscoelasticity is:
In formula:Density/the kgm of ρ --- medium-3;Angular frequency/the rads of ω --- vibration-1;μ1--- shear viscosity Coefficient/Pas;μ2--- the coefficient of rigidity/kPa.
Further, tested tissue and vibration measurement with laser equipment are respectively connected with a three-dimensional moving device.
Further, the pulse train that AWG is produced is long pulse signal.
Further, it is the HIFU transducers of 1MHz for centre frequency, pulse width is between 50 μ s-500 μ s.
Relative to prior art, the present invention has following significant advantage:
High intensity focused ultrasound is evaluated the invention provides a kind of vibration measurement with laser based on acoustic radiation force impulse response (HIFU) system and method for viscoelastic property is damaged, the deficiency of existing SDUV methods is overcome.The present invention is made using impulse wave For organizing, not only process is simple, take short, strong interference immunity, moreover it is possible to effectively reduce ultrasonic radiation dosage;Vibration measurement with laser technology Shear wave velocity is calculated by tracking the vibration of organization internal point, the evaluation of hyaline tissue's internal injury is can be effectively used for;To swash Light vibration measuring technology is combined with the method for acoustic radiation force impulse response, measurement is responded faster, and precision is higher, overcomes existing side Method is difficult to be separated to damaging the deficiency being monitored in speed and viscoplasticity simultaneously;Vibration measurement with laser apparatus simple structure, Without extra fixation, and the vibration information of measurement point can be directly obtained, later data treatment is simple;Laser signal and HIFU without Interact, the vibration of local organization in monitor in real time and detection therapeutic process is can be used to completely.
Brief description of the drawings
The present invention is described in further details with reference to the accompanying drawings and detailed description.
Fig. 1 is that a kind of vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response of the present invention damages viscoelastic method The block diagram of used monitoring system.Wherein:1 is synchronizing signal control system, and 2 is AWG, and 3 is power amplification Device, 4 is impedance matching network, and 5 is high-intensity focusing ultrasonic transducer (HIFU), and 6 is water tank, and 7 is tested tissue, and 8 is data Capture card, 9 is vibration measurement with laser equipment, and 10 is PC, and 11 is three-dimensional moving device.
Fig. 2 is driving pulse schematic diagram in the present invention.
Fig. 3 is the schematic diagram of calculating shear wave velocity method in the present invention.
Fig. 4 is the schematic diagram of calculating shearing wave vibration frequency method in the present invention.
Fig. 5 is the target damage and mark point distribution map after treatment.Wherein (a) is LDV monitoring lower target damage and mark The form that point has just been completed, (b) is the lower target damage of LDV monitoring and the form after mark point is placed a hour, and (c) is horizontal damage Wound and mark point shape, (d) is longitudinal damage and mark point shape.
Fig. 6 is HIFU transducer acoustic pressure distribution maps.
Fig. 7 is the vibrational waveform time domain and frequency domain distribution comparison diagram of one group of excitation point and test point.A () is one group of measurement point Characteristic shear ripple amplitude and time graph contrast;B () corresponds to (a) two pulse frequency domain contrasts;C () is filtered through band logical The time domain waveform of the monochromatic after ripple;D () corresponds to two pulse frequency domains contrast of (c).
Fig. 8 be the shear wave velocity that is produced in imitative body of pulse excitation and the relaxation curve vibration angular frequency of distinct pulse widths point Butut.Wherein (a) is the shear wave velocity distribution that distinct pulse widths are produced;B () is the relaxation curve vibration angular frequency that distinct pulse widths are produced Rate is distributed;C () is the average clearance method for being located at 1~4 mark point with pulse width variations curve;D () is the pine for being located at 1-4 mark points The Mean Oscillation angular frequency of relaxation curve is with pulse width variations curve.
Fig. 9 (a) is the coefficient of rigidity distribution of target damage and mark point, and Fig. 9 (b) is target damage and mark point Shear viscosity is distributed.
Specific embodiment
The shear viscosity coefficient μ of the imitative body of the known polyacrylamide gel containing bovine serum albumin (BSA)1, shearing elasticity Modulus μ2, shear-wave velocity c in the imitative bodysThe relation of the angular frequency vibrated with shearing wave is as follows:
In formula:Density/the kgm of ρ --- medium-3;Angular frequency/the rads of ω --- vibration-1;μ1--- shear viscosity Coefficient/Pas;μ2--- the coefficient of rigidity/kPa.
From above formula, under the premise of known to the imitative volume densitys of BSA, shearing elasticity and viscous modulus can be marked by multiple The shear wave velocity change and corresponding vibration angular frequency for remembering point are fitted and obtain.
Be obtain medium in shear wave velocity, present invention utilizes same-phase point (generally selecting Oscillation Amplitude peak point) when Between it is poor:The time that point will be encouraged to reach amplitude peak is designated as 0, then the time that test point reaches amplitude peak immediately is exactly shearing wave The time of propagation.Due in the present invention encourage point fixed with the location interval of test point, therefore by 2 points distance be spaced divided by Time interval is the average propagation speed that can obtain shearing wave between this 2 points.
To obtain the frequency information of shearing wave vibration, present invention utilizes on the vibration of test point and time curve two The time difference of extreme point:When wave amplitude arrival peak value is sheared, especially viscoplasticity is determined the mechanical characteristic of local organization Lax process, due to the influence of soft tissue viscosity, detected wave-shape amplitude is vibrated decay near equilbrium position, is declined The time difference of vibration amplitude maximum point and smallest point may be considered π/ω during subtracting, therefore can accordingly calculate vibration frequency The value of rate ω.
To obtain vibration of the excitation point with test point and time curve, the present invention takes following steps:Regulation laser Vibration measuring equipment is confocal with HIFU transducers, makes optical signal launch direction in the same direction with HIFU ultrasound emissions direction.Two three-dimensionals are set Mobile device, is respectively used to move vibration measurement with laser equipment 9 and imitative body 7.Before experiment, using HIFU transducers 5 in destination organization Several mark points are equally spacedly done in portion, are located on same straight line, while target damage is done on this line, the chi of mark point Very little and degree of injury is much smaller than target damage.Regulation HIFU energy transducer focus remember that the point is excitation point to a mark point.Synchronous letter Number control system 1 produces synchronizing signal feeding AWG 2 and then drives HIFU transducers 5 to swash to complete shearing wave Encourage.Trigger data acquisition card 8 while signal occurs, the signal to the output of vibration measurement with laser equipment 9 is acquired.Note mark point it Between spacing be Δ d, excitation point and test point is all located in mark point.Vibration at tracking excitation point obtains its wavy curve As the reference of subsequent detection, then using the adjacent point of the excitation point as test point, the oscillating curve at the point is recorded.Again The three-dimensional moving device of imitative body is carried as step-length movement with Δ d, above step is repeated, until having tested all mark points, obtained Multigroup vibrational waveform of the excitation point with test point.
The present invention is described in detail below in conjunction with the accompanying drawings.It should be appreciated that specific embodiment described herein is only used to The present invention is explained, is not intended to limit the present invention.
Refer to shown in Fig. 1, a kind of vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response of the present invention damages viscous The monitoring system that is used of method of elasticity, including synchronizing signal control system 1, AWG 2, power amplifier 3, Impedance matching network 4, high-intensity focusing ultrasonic transducer (HIFU) 5, water tank 6, tested tissue 7, data collecting card 8, Laser Measuring Equipment of shaking 9, PC 10 and two three-dimensional moving devices 11.Here tested tissue is that BSA imitates body, and vibration measurement with laser equipment uses laser Doppler vibrometer (LDV).
Synchronizing signal control system 1 connects AWG 2 and data collecting card 8, and AWG 2 passes through Power amplifier 3, the connection of impedance matching network 4 HIFU5;HIFU5 is arranged at outside water tank 6, and tested tissue 7 is arranged at water tank 6 Internal and one three-dimensional moving device of connection;Vibration measurement with laser equipment 9 is arranged at outside water tank 6, and it is mobile to connect another three-dimensional Device;The output end of vibration measurement with laser equipment 9 connects PC 10 through data collecting card 8.
Synchronizing signal control system 1 output synchronizing signal drive all the way AWG 2 (AWG420, Tektronix Inc., US) pulse train is produced, by (AG1017, the T&C Power of radio-frequency power amplifier 3 Conversion Inc., US) amplify with impedance matching network 4 after encourage HIFU transducers 5, the imitative bodies 7 of BSA in water tank 6 are applied Plus effect, the HIFU transducers used in experiment are sphere list element transducer, and centre frequency is 1.06MHz, and focal length is 100mm, burnt section length is 18mm (Chongqing Haifu Medical Technology Co., Ltd, CN);LDV9 is to cutting The vibration for cutting ripple generation is tracked, and another road trigger data acquisition card 8 is acquired to LDV signals in synchronizing signal.Collection The image for obtaining is preserved and processed by PC 10.Realized to the imitative bodies of LDV and BSA using three-dimensional moving device 11 in experiment It is mobile.
The imitative bodies 6 of BSA are the imitative body of bovine serum albumin (BSA) polyacrylamide gel that mass fraction is 7%, and table 1 is imitative for BSA The formula of body, wherein TRIS are used to adjust the PH to 8 of imitative body, and acrylamide is used to increase imitative body hardness, and TEMED is imitative for being catalyzed Body condenses.
The imitative liquid solution formulas of the 100mlBSA of table 1
A kind of vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response of the present invention damages viscoelastic method, including Following steps:Synchronizing signal control system drives AWG to produce pulse train to be used for generation and shears wave excitation.Drive Dynamic signal is long pulse signal, and all pulses are once launched, due to the HIFU transducers centre frequency that is used in this experiment for 1.06MHz, therefore be 50 μ s-500 μ s from pulse width (PD), according to the HIFU transducers of other centre frequencies, then may be used It is appropriate to change pulse width value, and be step-length with 50 μ s, the pulse of each width is made up of 10 group pulses, and pulse repeats frequency Rate (PRF) is 100Hz.To prevent temperature too high, radiating, the sound work(of HIFU excitations are ensured per 30ms is spaced between group pulse string Rate is fixed on 12W.The vibration that shearing wave is produced is tracked by LDV, while synchronizing signal control system produces pulse train Trigger data acquisition card, is acquired to LDV signals.The power output of power amplifier is 50W.Adjusted first before experiment LDV and HIFU transducers are confocal, make the laser signal direction of the launch of LDV in the same direction with the HIFU transducer ultrasonic direction of the launch.It is three-dimensional Mobile device is used for moving LDV and imitative body, completes the vibration acquisition of multiple positions.In the present embodiment, the mark point for being used All it is fire damage with target damage, HIFU transducers are driven by continuous sine wave, and acoustical power is fixed on 60W.Damage by LDV It is camera supervised, lesion size is controlled by adjusting treatment time.When treatment end, mark point is equally spacedly distributed in Around target damage, and positioned at same depth, mark point is enough to reflection laser beam, and will not be to the tissue characteristics of background Generation is significantly affected.
In experiment, 13 mark points are had as shown in Figure 5, wherein mark point 7 is located at target damage center.The chi of mark point Very little and degree of injury is much smaller than target damage.It is the vibration at excitation point detection excitation point as follow-up inspection first using mark point 1 The reference of survey, now positioned at mark point 2, the mobile three-dimensional moving device for carrying LDV, the vibration to test point 2 is carried out test point Detection.Move again and carry the three-dimensional moving device of the imitative bodies of BSA in causing HIFU energy transducer focus to be located at the center of test point 2, with Test point 2 is excitation point, repeats above detection process, until having tested all mark points, obtains 12 groups of oscillating curves.Every kind of width The pulse detection step of degree is identical.
The magnitude peak of same group of excitation point and test point is taken respectively, remembers its time at intervals of Δ T1, due to the present embodiment Middle excitation point is fixed as 1mm with test point location interval, therefore it is to be cut to be spaced divided by time interval by 2 points of distance Cutting average propagation speed of the ripple between this 2 points is:
The vibration of analysis test point and the relation curve of time, the time difference for remembering two extreme points is Δ T2, then vibration frequency The value of ω is
Angular frequency is vibrated according to the shear wave velocity change between multiple mark points obtained as above and corresponding difference, with reference to
Shear viscosity and elastic modelling quantity are fitted, the viscoplasticity distribution of local organization is just can obtain.
Interpretation of result:
As shown in figure 5, mark point is equally spacedly distributed in target damage both sides, diameter is about 0.8mm, the chi of target damage Very little and degree of injury is much larger than mark point, and it is axially about 7mm, is laterally about 2mm.By the sound intensity point of HIFU transducers in Fig. 6 Butut understands that burnt area's axial dimension is about 20mm, and lateral dimension is 3mm or so.Therefore, if transducers focus are being damaged Center, the then vibration produced by pulse excitation can cover whole damage field.The plane shear ripple produced by excitation is being propagated through The vibration for causing mark point and target damage to produce in journey can be as measured by laser in mark point and target damage vibration Represent.
The vibrational waveform time domain and frequency domain distribution pair of one group of excitation point when being illustrated in figure 7 PD=100 μ s and test point Than figure.It can be seen that the time domain waveform curve at this 2 point has similitude, and there is obvious phase delay.Each wavy curve is by two Part constitutes:Part I is that curve is vibrated to amplitude peak;Part II is displacement decayed to by amplitude peak vibration it is flat Weighing apparatus position.This is consistent with the Nomenclature Composition and Structure of Complexes of typical curve in Fig. 4.On frequency domain, by frequency domain figure as can be seen that excitation point As being with the shearing wave frequency spectrum of test point, it can be considered that communication process of the shearing wave between excitation point and test point In there is no occurrence frequency to change, and bandpass filtering has effectively filtered into monochromatic original waveform.
Angular frequency and shear wave velocity are vibrated according to the relaxation curve that the wavy curve calculates distinct pulse widths generation respectively, as a result As shown in Figure 7.
It is illustrated in figure 9 the fitting result figure of the coefficient of rigidity and shear viscosity at target damage and mark point.The figure Show, the distribution curve of the coefficient of rigidity and shear viscosity is similar, in the target damage area coefficient of rigidity and shearing Viscous modulus have substantially rising compared with mark zone, therefore can significantly tell damage field and mark point region.
Specific embodiment of the invention is above examples described, therefore its description is more specific and detailed, but can not be And it is interpreted as the limitation to the scope of the claims of the present invention.It should be pointed out that make without departing from the inventive concept of the premise Some deformations and improvement, belong to protection scope of the present invention.Therefore, the protection domain of patent of the present invention should be wanted with appended right Ask and be defined.

Claims (9)

1. the vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoelastic method, it is characterised in that including Following steps:
1) control HIFU transducers produce driving pulse to act on the excitation point inside tested tissue;Vibration measurement with laser equipment is to tested The vibration of organization internal test point is tracked;Record the oscillating curve of each test point;
2) according to oscillating curve, shear wave velocity and frequency are calculated, then carries out viscoplasticity fitting.
2. the vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response according to claim 1 damages viscoelastic side Method, it is characterised in that step 1) in, regulation vibration measurement with laser equipment is confocal with HIFU transducers, sends out vibration measurement with laser device signal Penetrate direction in the same direction with the HIFU transducer ultrasonic direction of the launch;Synchronizing signal control system produces synchronizing signal to drive any ripple respectively Shape generator and data collecting card, the pulse train that AWG is produced is by power amplifier and impedance matching network HIFU transducers are acted on, control HIFU transducers produce driving pulse to act on the excitation point in tested tissue inside;Data Capture card is acquired to the signal that vibration measurement with laser equipment is exported.
3. the vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response according to claim 1 damages viscoelastic side Method, it is characterised in that the tested tissue is transparent prosthese or transparent in vitro tissue.
4. the vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response according to claim 1 damages viscoelastic side Method, it is characterised in that vibration measurement with laser equipment is laser doppler vibrometer.
5. the vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response according to claim 1 damages viscoelastic side Method, it is characterised in that step 1) specifically include following steps:
Step S1, by HIFU transducers focus at a certain position in tested tissue inside, the point be called to encourage point, HIFU transducers Produce pulse to act on point induction and produce shearing wave signal;
Step S2, using vibration measurement with laser equipment to excitation point and shearing direction of wave travel up to lack two vibrations of test point carry out Monitor and record its oscillating curve.
6. the vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response according to claim 5 damages viscoelastic side Method, it is characterised in that the method for calculating shear wave velocity and frequency is specially:
Mark point position is fixed and spacing with the time delay of same-phase between 2 points, it is known that represent the shear wave velocity between 2 points, Excitation point is designated as 0 in the time of a certain phase, test point is in the time that the time of same distance is exactly that shearing wave is propagated;Note swashs It is a little Δ d with the distance between test point to encourage, and the time intervals of 2 points of arrival same-phases are Δ T1, then shearing wave is in test point Speed is:
c s = Δ d ΔT 1 - - - ( 1 )
Time difference in shearing wave oscillating curve between minimum and maximum 2 points of relaxation amplitude is considered π/ω, therefore, it is possible to Calculate its frequency;Remember that this 2 time differences are Δ T2, then shearing wave vibration frequency be
ω = π ΔT 2 - - - ( 2 )
In the case where local organization shear wave velocity and frequency distribution is drawn, Voigt models and its derivation according to biological tissue Formula carries out viscoplasticity fitting;
Method for carrying out viscoplasticity fitting is:
Relation between plane shear velocity of wave and tissue viscoelasticity is:
c s = 2 ( μ 1 2 + ω 2 μ 2 2 ) ρ ( μ 1 + ( μ 1 2 + ω 2 μ 2 2 ) ) - - - ( 3 )
In formula:Density/the kgm of ρ --- medium-3;Angular frequency/the rads of ω --- vibration-1;μ1--- shear viscosity system Number/Pas;μ2--- the coefficient of rigidity/kPa.
7. the vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response according to claim 1 damages viscoelastic side Method, it is characterised in that tested tissue and vibration measurement with laser equipment are respectively connected with a three-dimensional moving device.
8. the vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response according to claim 2 damages viscoelastic side Method, it is characterised in that the pulse train that AWG is produced is long pulse signal.
9. the vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response according to claim 8 damages viscoelastic side Method, it is characterised in that for the HIFU transducers that centre frequency is 1MHz, pulse width is between 50 μ s-500 μ s.
CN201710166704.XA 2017-03-20 2017-03-20 Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoplasticity method Expired - Fee Related CN106840362B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710166704.XA CN106840362B (en) 2017-03-20 2017-03-20 Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoplasticity method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710166704.XA CN106840362B (en) 2017-03-20 2017-03-20 Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoplasticity method

Publications (2)

Publication Number Publication Date
CN106840362A true CN106840362A (en) 2017-06-13
CN106840362B CN106840362B (en) 2019-08-23

Family

ID=59130053

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710166704.XA Expired - Fee Related CN106840362B (en) 2017-03-20 2017-03-20 Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoplasticity method

Country Status (1)

Country Link
CN (1) CN106840362B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108732240A (en) * 2018-04-26 2018-11-02 西安交通大学 Laser Doppler monitoring pulse acoustic radiation force shearing wave quantitative predications HIFU damages viscoelastic system and method
CN110251159A (en) * 2019-05-31 2019-09-20 深圳先进技术研究院 A kind of tissue testing method, device and storage medium
CN110954207A (en) * 2019-11-15 2020-04-03 重庆医科大学 Detection device and detection method for focused ultrasound focus sound wave structure
CN111307269A (en) * 2020-03-11 2020-06-19 北京理工大学 Laser confocal/differential confocal Raman spectrum vibration parameter measuring method
CN111751200A (en) * 2020-06-29 2020-10-09 中国科学院声学研究所 System and method for testing dynamic modulus of viscoelastic material
WO2021212578A1 (en) * 2020-04-24 2021-10-28 深圳迈瑞生物医疗电子股份有限公司 Viscoelasticity measuring method, elasticity measuring method and ultrasonic measuring system
CN114838802A (en) * 2022-04-26 2022-08-02 南京工业职业技术大学 Touch slider sub-nanoscale harmonic detection system based on Hilbert-Huang transform
CN115753502A (en) * 2022-11-14 2023-03-07 西安交通大学 Device and method for testing micro-nano rheological characteristics of biological tissues

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009140607A1 (en) * 2008-05-15 2009-11-19 Mayo Foundation For Medical Education And Research Vibration generation and detection in shear wave dispersion ultrasound vibrometry with large background motions
CN102151152A (en) * 2011-03-01 2011-08-17 深圳市一体医疗科技股份有限公司 Measurement probe, system and method for measuring elasticity of viscoelastic medium
CN103006274A (en) * 2012-12-21 2013-04-03 深圳大学 Method and system for ultrasonic detection of cornea viscoelasticity
CN106037816A (en) * 2016-05-19 2016-10-26 深圳大学 Ultrasonic viscoelasticity measurement method and system
CN106175831A (en) * 2016-09-28 2016-12-07 重庆理工大学 A kind of tissue viscoelasticity based on shearing wave amplitude and phase-detection measuring method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009140607A1 (en) * 2008-05-15 2009-11-19 Mayo Foundation For Medical Education And Research Vibration generation and detection in shear wave dispersion ultrasound vibrometry with large background motions
CN102151152A (en) * 2011-03-01 2011-08-17 深圳市一体医疗科技股份有限公司 Measurement probe, system and method for measuring elasticity of viscoelastic medium
CN103006274A (en) * 2012-12-21 2013-04-03 深圳大学 Method and system for ultrasonic detection of cornea viscoelasticity
CN106037816A (en) * 2016-05-19 2016-10-26 深圳大学 Ultrasonic viscoelasticity measurement method and system
CN106175831A (en) * 2016-09-28 2016-12-07 重庆理工大学 A kind of tissue viscoelasticity based on shearing wave amplitude and phase-detection measuring method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
林浩铭 等: "基于正交频率脉冲编码激励的软组织粘弹性定量测量", 《中国生物医学工程学报》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108732240A (en) * 2018-04-26 2018-11-02 西安交通大学 Laser Doppler monitoring pulse acoustic radiation force shearing wave quantitative predications HIFU damages viscoelastic system and method
CN108732240B (en) * 2018-04-26 2020-06-19 西安交通大学 System and method for quantitatively estimating HIFU damage viscoelasticity by laser Doppler monitoring pulsed acoustic radiation force shear wave
CN110251159A (en) * 2019-05-31 2019-09-20 深圳先进技术研究院 A kind of tissue testing method, device and storage medium
CN110954207A (en) * 2019-11-15 2020-04-03 重庆医科大学 Detection device and detection method for focused ultrasound focus sound wave structure
CN111307269A (en) * 2020-03-11 2020-06-19 北京理工大学 Laser confocal/differential confocal Raman spectrum vibration parameter measuring method
WO2021212578A1 (en) * 2020-04-24 2021-10-28 深圳迈瑞生物医疗电子股份有限公司 Viscoelasticity measuring method, elasticity measuring method and ultrasonic measuring system
CN111751200A (en) * 2020-06-29 2020-10-09 中国科学院声学研究所 System and method for testing dynamic modulus of viscoelastic material
CN111751200B (en) * 2020-06-29 2021-06-15 中国科学院声学研究所 System and method for testing dynamic modulus of viscoelastic material
CN114838802A (en) * 2022-04-26 2022-08-02 南京工业职业技术大学 Touch slider sub-nanoscale harmonic detection system based on Hilbert-Huang transform
CN115753502A (en) * 2022-11-14 2023-03-07 西安交通大学 Device and method for testing micro-nano rheological characteristics of biological tissues
CN115753502B (en) * 2022-11-14 2023-08-18 西安交通大学 Device and method for testing micro-nano rheological characteristics of biological tissue

Also Published As

Publication number Publication date
CN106840362B (en) 2019-08-23

Similar Documents

Publication Publication Date Title
CN106840362B (en) Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoplasticity method
CN106730424B (en) Hundred microsecond pulse ultrasonic tissue of confocal harmonic superposition damages mode control method
CN110248606B (en) Air pocket positioning
CN107550458B (en) Biological tissue multi-characteristic imaging method based on acoustoelectric effect and acoustic radiation force
CN101784234A (en) Method and device for measuring a mean value of visco-elasticity of a region of interest
Maleke et al. Single-element focused ultrasound transducer method for harmonic motion imaging
WO2005002447A1 (en) Harmonic motion imaging
Zhao et al. Shear wave speed measurement using an unfocused ultrasound beam
CN110243454B (en) Micro-bubble synergy dynamics experiment system and method based on dual-frequency superposed ultrasonic pulses in viscoelastic medium
CN109077754A (en) A kind of method and apparatus measuring organization mechanics characterisitic parameter
Hazard et al. Integration of crawling waves in an ultrasound imaging system. Part 1: System and design considerations
CN109745077A (en) Based on the elastic characteristic detection method for focusing ultrasonic sound and vibration signal
JP7295129B2 (en) Hybrid elastography method, probes and devices for hybrid elastography
CN109730722A (en) Based on the elastograph imaging method for focusing ultrasonic sound and vibration signal
CN108732240B (en) System and method for quantitatively estimating HIFU damage viscoelasticity by laser Doppler monitoring pulsed acoustic radiation force shear wave
CA2667527C (en) Method for generating mechanical waves by creating an interfacial acoustic radiation force
Lu et al. Enhanced-cavitation heating protocols in focused ultrasound surgery with broadband split-focus approach
CN110243950B (en) Acoustic radiation force pulse excitation method based on amplitude modulation
Zheng et al. Shear Wave Propagation in Soft Tissue with Ultrasound Vibrometry
Miller et al. Bubble-induced color Doppler feedback for histotripsy tissue fractionation
JP2017164559A (en) Ultrasonic device
KR101411099B1 (en) Method for generating mechanical waves by creating an interfacial acoustic radiation force
CN105246416A (en) Ultrasonic medical device
CN112023284A (en) Focus position real-time monitoring method for focused acoustic dynamic therapy
Amador et al. Effects of phase aberration on acoustic radiation force-based shear wave generation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190823