CN106840362B - Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoplasticity method - Google Patents
Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoplasticity method Download PDFInfo
- Publication number
- CN106840362B CN106840362B CN201710166704.XA CN201710166704A CN106840362B CN 106840362 B CN106840362 B CN 106840362B CN 201710166704 A CN201710166704 A CN 201710166704A CN 106840362 B CN106840362 B CN 106840362B
- Authority
- CN
- China
- Prior art keywords
- vibration
- hifu
- vibration measurement
- point
- viscoplasticity
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring 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 damage viscoplasticity method based on acoustic radiation force impulse response, comprising the following steps: 1) controls HIFU energy converter generation driving pulse and act on the point of the excitation inside tested tissue；Vibration measurement with laser equipment tracks the vibration of tested tissue inside test point；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 motivate hyaline tissue inside to generate 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 motion information, that is, shear wave velocity and the frequency of each mark point, and then realizes the real time monitoring and evaluation damaged to hyaline tissue inside HIFU.
Description
Technical field
The invention belongs to biomedical measurement technical field, in particular to a kind of monitoring high intensity focused ultrasound (HIFU) damage
Hurt viscoplasticity method.
Background technique
High intensity focused ultrasound (High Intensity Focused Ultrasound, HIFU) is a kind of noninvasive, peace
Complete and efficient treatment method, it focuses high energy ultrasound beam by the HIFU energy converter of placement in vitro, and in focal zone
Biological modification is generated in domain to realize the treatment of disease.In order to effectively guide treatment focal regions, the control ultrasonic therapy of HIFU
Dosage guarantees that HIFU can be applied to clinic to effective and safe, needs that the therapeutic process of HIFU is monitored in real time and evaluated.?
During HIFU Treatment, due to the dissipation effect of tissue, so that supersonic beam kinetic energy in communication process is converted into thermal energy, it is located at coke
The histocyte in area is destroyed, protein denaturation to forming irreversible coagulation necrosis, cause tissue hardness variation from
And change its mechanical characteristic.Therefore, it can be supervised in real time by the especially viscoelastic property of the mechanical characteristic to damage field
Control is 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 not only separate viscosity and elasticity, but also available they are each
From quantized value, be a kind of very effective detection method, have been applied in in vitro and living tissue.
Local organization is activated to steady-state response using the continuous wave of amplitude modulation by existing SDUV method, to generate harmonic wave
Wave motion is sheared, then measures the shear wave velocity under different frequency respectively, then the shear viscosity and bullet of tissue is obtained by calculation
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, this method is adopted
Excitation and motion detection are carried out to tissue respectively with a ultrasonic transducer.By the continuous wave excitation tissue for applying certain frequency
It generates harmonic wave and shears wave motion, then obtain the harmonic wave shearing wave of two different locations on transonic direction using coherent plane wave
Motor message, the harmonic wave shearing wave motion information such as amplitude and phase etc. by analyzing different frequency can carry out viscoplasticity and estimate
Meter.In the method by modification Kalman filtering parameter come and meanwhile determine multiple harmonic motion frequencies, and during filtering will
Background motion is modeled as random process to compensate the movement 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 be motivated using continuous wave to tissue in interleaved push and detection pulses ".Unlike
By pumping signal transmitting alternate with detection signal in this method, so as to obtain in biggish vibration displacement and detection process
Secure threshold of the ultrasonic radiation dosage without departing from FDA.However, existing SDUV method needs 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.
It generallys use ultrasonic imaging in SDUV to be monitored shearing wave, but there are certain limitations for the method for ultrasonic imaging
Property, as image resolution ratio is low, echo vulnerable to noise and organism displacement interference and can be generated with HIFU sound field mutual
Effect etc., so that this method is unfavorable for monitoring and detecting HIFU Treatment process in real time.Vibration measurement with laser method can achieve sub-micron
Precision, and can be applied to the measurement in tissue surface and hyaline tissue, since laser signal does not influence ultrasonic sound field, because
This can be used to the various vibrations during real time monitoring HIFU Treatment.The patent of invention of James F.Greenleaf et al.
US9044192B2, patent name " System and method for non-invasively measuring tissue
The system using surface wave measurement tissue viscoelasticity and side are provided in viscoelasticity using surfacewaves "
Method generates vibration in tissue surface using energy converter or mechanical vibrator in the method, then using laser vibration measurer to group
The surface wave for knitting surface carries out tracking to estimate its viscoplasticity information.However, existing vibration measurement with laser method is confined to tissue
The detection of surface vibration effectively cannot be tracked and be detected to organization internal vibration.
Summary of the invention
The purpose of the present invention is to provide a kind of, and the vibration measurement with laser monitoring HIFU damage based on acoustic radiation force impulse response is viscous
Elastic method acts on tissue using pulse, and vibration measurement with laser technology is applied to the monitoring of hyaline tissue's internal shear wave
In, it is transparent to overcome the shortcomings of to act on tissue using continuous wave in existing SDUV, and realize for vibration measurement with laser technology to be applied to
The breakthrough of organization internal.
To achieve the above object, the present invention adopts the following technical scheme that:
Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoplasticity method, comprising the following steps:
1) HIFU energy converter generates driving pulse and acts on the point of the excitation 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, it in step 1), adjusts vibration measurement with laser equipment and HIFU energy converter is confocal, make vibration measurement with laser device signal
The direction of the launch and the HIFU transducer ultrasonic direction of the launch are in the same direction；Synchronization signal control system generates synchronization signal and respectively drives arbitrarily
Waveform generator and data collecting card, the pulse train that arbitrary waveform generator generates pass through power amplifier and impedance matching net
Network acts on HIFU energy converter, and control HIFU energy converter generates driving pulse and acts on the point of the excitation inside tested tissue；Data
The signal that capture card exports vibration measurement with laser equipment is acquired.
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 includes the following steps:
Step S1, by HIFU transducers focus at a certain position in tested tissue inside, which is referred to as excitation point, and HIFU is changed
Energy device generates pulse and acts on point induction generation shearing wave signal；
Step S2, direction of wave travel to excitation point and is sheared up to the vibration of few two test points using vibration measurement with laser equipment
It is monitored and records its oscillating curve.
Further, the method for shear wave velocity and frequency is calculated specifically:
Mark point position is fixed and spacing is it is known that indicate the shearing wave between two o'clock with the delay of same-phase between two o'clock
Excitation point is denoted as 0 in the time of a certain phase by speed, 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 that two o'clock reaches same-phase is Δ T1, then shearing wave is detecting
The speed of point are as follows:
Time difference in shearing wave oscillating curve between the minimum and maximum two o'clock of relaxation amplitude is considered π/ω, therefore
Its frequency can be calculated；Remember that this two o'clock time difference is Δ T2, then shearing wave vibration frequency be
In the case where obtaining local organization shear wave velocity and frequency distribution, according to the Voigt model of biological tissue and its
Derivation formula carries out viscoplasticity fitting；
Method for carrying out viscoplasticity fitting are as follows:
Relationship between plane shear velocity of wave and tissue viscoelasticity are as follows:
In formula: ρ --- density/kgm of medium-3；Angular frequency/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 arbitrary waveform generator generates is long pulse signal.
Further, the HIFU energy converter for being 1MHz for centre frequency, pulse width is between 50 μ s-500 μ s.
Compared with the existing technology, the present invention has following significant advantage:
The present invention provides a kind of, and the vibration measurement with laser based on acoustic radiation force impulse response evaluates high intensity focused ultrasound
(HIFU) system and method for damaging viscoelastic property, overcomes the deficiency of existing SDUV method.The present invention is made using impulse wave
For organizing, not only process is simple, the short, strong interference immunity of time-consuming, 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, can be effectively used for the evaluation of hyaline tissue's internal injury；It will swash
Flash ranging vibration technology is combined with the method for acoustic radiation force impulse response, makes measurement response faster, precision is higher, overcomes existing side
Method is difficult to the deficiency being monitored in speed and viscoplasticity separation to damage simultaneously；Vibration measurement with laser apparatus structure is simple,
Without additional fixation, and the vibration information of measurement point can be directly obtained, later data processing is simple；Laser signal and HIFU without
Interaction can be used to monitor in real time and detect the vibration of local organization in therapeutic process completely.
Detailed description of the invention
The present invention is described in further details with reference to the accompanying drawings and detailed description.
Fig. 1 is a kind of vibration measurement with laser monitoring HIFU damage viscoplasticity method institute based on acoustic radiation force impulse response of the present invention
Using the block diagram of monitoring system.Wherein: 1 is synchronization signal control system, and 2 be arbitrary waveform generator, and 3 be power amplifier, 4
It is high-intensity focusing ultrasonic transducer (HIFU) for impedance matching network, 5,6 be water tank, and 7 be tested tissue, and 8 acquire for data
Card, 9 be vibration measurement with laser equipment, and 10 be PC machine, and 11 be three-dimensional moving device.
Fig. 2 is driving pulse schematic diagram in the present invention.
Fig. 3 is the schematic diagram that shear wave velocity method is calculated in the present invention.
Fig. 4 is the schematic diagram that shearing wave vibration frequency method is calculated in the present invention.
Fig. 5 is the target damage and mark point distribution map after treatment.Wherein (a) is that LDV monitors lower target damage and label
The form that point has just been completed is (b) that LDV monitors lower target damage and the form after mark point placement one hour, is (c) lateral damage
Wound and label dot shape are (d) longitudinal damage and label dot shape.
Fig. 6 is HIFU energy converter acoustic pressure distribution map.
Fig. 7 is the vibrational waveform time domain and frequency domain distribution comparison diagram of one group of excitation point and test point.It (a) is one group of measurement point
Characteristic shear wave amplitude and time graph comparison；(b) (a) two pulse frequency domain comparisons are corresponded to；It (c) is filtered through band logical
The time domain waveform of monochromatic after wave；(d) two pulse frequency domains comparison of (c) is corresponded to.
Fig. 8 is point of the shear wave velocity that the pulse excitation of distinct pulse widths generates in imitative body and relaxation curve vibration angular frequency
Butut.Wherein (a) is the shear wave velocity distribution that distinct pulse widths generate；It (b) is that the relaxation curve that distinct pulse widths generate vibrates angular frequency
Rate distribution；It (c) is the average clearance method positioned at 1~4 mark point with pulse width variations curve；It (d) is the pine for being located at 1-4 mark point
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 distribution.
Specific embodiment
The known polyacrylamide gel containing bovine serum albumin (BSA) imitates the shear viscosity coefficient μ of body1, shearing elasticity
Modulus μ2, shear-wave velocity c in the imitative bodysIt is as follows with the relationship of the angular frequency of shearing wave vibration:
In formula: ρ --- density/kgm of medium-3；Angular frequency/rads of ω --- vibration-1；μ1--- shear viscosity
Coefficient/Pas；μ2--- the coefficient of rigidity/kPa.
From the above equation, we can see that shearing elasticity and viscous modulus can pass through multiple marks under the premise of known to BSA imitates volume density
The shear wave velocity variation of note point and corresponding vibration angular frequency are fitted to obtain.
For obtain medium in shear wave velocity, present invention utilizes same-phase point (usually selecting Oscillation Amplitude peak point) when
Between it is poor: excitation point is reached into time of amplitude peak and is denoted as 0, then it is exactly shearing wave that test point reaches the time of amplitude peak immediately
The time of propagation.Since excitation point and the location interval of test point are fixed in the present invention, by two o'clock distance interval divided by
Average propagation rate of the shearing wave between this two o'clock can be obtained in time interval.
For the frequency information for obtaining shearing wave vibration, present invention utilizes on the vibration of test point and time curve two
The time difference of extreme point: when shearing wave amplitude arrival peak value, the mechanical characteristic of local organization especially viscoplasticity is determined
The process of relaxation, due to the influence of soft tissue viscosity, detected wave-shape amplitude vibrates decaying near equilbrium position, decline
The time difference of vibration amplitude maximum point and smallest point may be considered π/ω during subtracting, therefore can calculate vibration frequency accordingly
The value of rate ω.
Following steps are taken to obtain vibration and time curve, the present invention of the excitation point with test point: adjusting laser
Vibration measuring equipment and HIFU energy converter are confocal, keep optical signal launch direction and HIFU ultrasound emission direction in the same direction.Two three-dimensionals are set
Mobile device is respectively used to mobile vibration measurement with laser equipment 9 and tested tissue 7.Before experiment, using HIFU energy converter 5 in target group
It knits inside and equally spacedly does several mark points, be located on same straight line, while doing target damage on this line, mark point
Size and degree of injury be much smaller than target damage.HIFU energy transducer focus is adjusted to a mark point, remembers the point for excitation point.Together
Whistle control system 1 is walked to generate synchronization signal feeding arbitrary waveform generator 2 and then drive HIFU energy converter 5 to complete shearing wave
Excitation.Trigger data acquisition card 8 while signal occurs, the signal exported to vibration measurement with laser equipment 9 are acquired.Remember mark point
Between spacing be Δ d, make excitation point and test point be all located in mark point.Vibration at tracking excitation point obtains its waveform song
Line is then put adjacent point as test point using the excitation, is recorded the oscillating curve at the point using the reference as subsequent detection.
Again using Δ d as the mobile three-dimensional moving device for carrying imitative body of step-length, repeats above step and obtained until testing all mark points
Vibrational waveform to multiple groups excitation point with test point.
The present invention will be described in detail with reference to the accompanying drawing.It should be appreciated that specific embodiment described herein only to
It explains the present invention, is not intended to limit the present invention.
Refering to Figure 1, a kind of vibration measurement with laser monitoring HIFU damage based on acoustic radiation force impulse response of the present invention is viscous
Monitoring system used in elastic method, including synchronization signal control system 1, arbitrary waveform generator 2, power amplifier 3, resistance
Anti- matching network 4, high-intensity focusing ultrasonic transducer (HIFU) 5, water tank 6, tested tissue 7, data collecting card 8, vibration measurement with laser
Equipment 9, PC machine 10 and two three-dimensional moving devices 11.Here tested tissue is that BSA imitates body, and vibration measurement with laser equipment is more using laser
General Le vialog (LDV).
Synchronization signal control system 1 connects arbitrary waveform generator 2 and data collecting card 8, arbitrary waveform generator 2 pass through
Power amplifier 3, impedance matching network 4 connect HIFU5；HIFU5 is set to outside water tank 6, and tested tissue 7 is set to water tank 6
Internal and one three-dimensional moving device of connection；Vibration measurement with laser equipment 9 is set to 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 machine 10 through data collecting card 8.
Synchronization signal control system 1 export synchronization signal drive all the way arbitrary waveform generator 2 (AWG420,
Tektronix Inc., US) pulse train is generated, by (AG1017, the T&C Power of radio-frequency power amplifier 3
Conversion Inc., US) with impedance matching network 4 amplify after motivate HIFU energy converter 5, to the tested tissue 7 in water tank 6
(BSA imitates body) application effect, HIFU energy converter used in experiment are spherical surface list element transducer, centre frequency 1.06MHz,
Focal length is 100mm, and burnt section length is 18mm (Chongqing Haifu Medical Technology Co., Ltd, CN)；
The vibration that LDV9 generates shearing wave tracks, and another way trigger data acquisition card 8 adopts LDV signal in synchronization signal
Collection.The image collected is saved and is handled by PC machine 10.Realized to LDV in experiment using three-dimensional moving device 11 and
BSA imitates the movement of body.
It is that bovine serum albumin (BSA) polyacrylamide gel that mass fraction is 7% imitates body that BSA, which imitates body 6, and table 1 is imitative for BSA
The formula of body, wherein TRIS is used to adjust the PH to 8 of imitative body, and for acrylamide for increasing imitative body hardness, TEMED is imitative for being catalyzed
Body condensation.
1 100mlBSA of table imitates liquid solution formula
A kind of vibration measurement with laser monitoring HIFU damage viscoplasticity method based on acoustic radiation force impulse response of the present invention, including with
Lower step: synchronization signal control system drives arbitrary waveform generator to generate pulse train for generating shearing wave excitation.Driving
Signal is long pulse signal, and all pulses once emit, due to the HIFU energy converter centre frequency that is used in this experiment for
1.06MHz, therefore selecting pulse width (PD) is that 50 μ s-500 μ s then may be used according to the HIFU energy converter of other centre frequencies
It is appropriate to change pulse width value, and using 50 μ s as step-length, the pulse of each width is made of 10 group pulses, and pulse repeats frequency
Rate (PRF) is 100Hz.It is spaced 30ms to prevent temperature excessively high, between every group pulse string to guarantee to radiate, the sound function of HIFU excitation
Rate is fixed on 12W.The vibration that shearing wave generates is tracked by LDV, while synchronization signal control system generates pulse train
Trigger data acquisition card is acquired LDV signal.The output power of power amplifier is 50W.It is adjusted first before experiment
LDV and HIFU energy converter is confocal, keeps the laser signal direction of the launch of LDV and the HIFU transducer ultrasonic direction of the launch in the same direction.It is three-dimensional
Mobile device is used to move LDV and imitative body, completes the vibration acquisition of multiple positions.In the present embodiment, used mark point
It is all thermal damage with target damage, HIFU energy converter is driven by continuous sine wave, and acoustical power is fixed on 60W.Damage passes through in LDV
It is camera supervised, lesion size controls by adjusting treatment time.When treatment end, mark point is equally spacedly distributed in
Around target damage, and it is located at same depth, mark point is enough reflection laser beam, and will not be to the tissue characteristics of background
Generation significantly affects.
In experiment, 13 mark points are shared as shown in Figure 5, and wherein mark point 7 is located at target damage center.The ruler of mark point
Very little and degree of injury is much smaller than target damage.It is first vibration at excitation point detection excitation point as subsequent inspection using mark point 1
The reference of survey, test point is located at mark point 2 at this time, and the three-dimensional moving device of mobile carrying LDV carries out the vibration of test point 2
Detection.Carrying BSA is moved again imitates the three-dimensional moving device of body in the center for making HIFU energy transducer focus be located at test point 2, with
Test point 2 is excitation point, repeats the above detection process until testing all mark points and obtains 12 groups of oscillating curves.It is every kind wide
The pulse detection step of degree is identical.
The magnitude peak for taking same group of excitation point and test point respectively remembers that its time interval is Δ T1, due to the present embodiment
Middle excitation point and test point location interval are fixed as 1mm, therefore can be obtained and cut divided by time interval by the distance interval of two o'clock
Cut average propagation rate of the wave between this two o'clock are as follows:
The vibration of analysis detection point and the relation curve of time, the time difference of two extreme points of note are Δ T2, then vibration frequency
The value of ω is
Angular frequencies are vibrated with corresponding difference according to the shear wave velocity variation between multiple mark points obtained as above, in conjunction with
Shear viscosity and elasticity modulus are fitted, the viscoplasticity distribution of local organization can be obtained.
Interpretation of result:
As shown in figure 5, mark point is equally spacedly distributed in target damage two sides, diameter is about 0.8mm, the ruler of target damage
Very little and degree of injury is much larger than mark point, and it is laterally about 2mm that axial, which is about 7mm,.By the sound intensity point of HIFU energy converter in Fig. 6
Butut is it is found that burnt area axial dimension is about 20mm, and lateral dimension is 3mm or so.Therefore, if transducers focus is damaging
Center, then vibration caused by pulse excitation can cover entire damage field.It is being propagated through by the plane shear wave that excitation generates
The vibration for causing mark point and target damage to generate in journey can by laser vibration measured by mark point and target damage Lai
It indicates.
The vibrational waveform time domain and frequency domain distribution pair of one group of excitation point and test point when being illustrated in figure 7 PD=100 μ s
Than figure.It can be seen that the time domain waveform curve at this two o'clock has similitude, and there is apparent phase delay.Each wavy curve is by two
Part forms: first part is that curve is vibrated to amplitude peak；Second part be displacement by amplitude peak vibration decay to it is flat
Weighing apparatus position.The Nomenclature Composition and Structure of Complexes of this and typical curve in Fig. 4 is consistent.On frequency domain, point is motivated it can be seen from frequency domain figure
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 frequency variations, and bandpass filtering has effectively filtered original waveform into monochromatic.
The relaxation curve vibration angular frequency and shear wave velocity of distinct pulse widths generation are calculated separately according to the wavy curve, 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 be it is similar, in the target damage area coefficient of rigidity and shearing
Viscous modulus has obvious rising compared with mark zone, therefore can significantly tell damage field and mark point region.
Therefore above examples describe a specific embodiment of the invention, and the description thereof is more specific and detailed, but can not be
And it is interpreted as limitations on the scope of the patent of the present invention.It should be pointed out that make without departing from the inventive concept of the premise
Several modifications and improvements belong to protection scope of the present invention.Therefore, the scope of protection of the patent of the present invention should be wanted with appended right
Subject to asking.
Claims (5)
1. based on acoustic radiation force impulse response vibration measurement with laser monitoring HIFU damage viscoplasticity method, which is characterized in that including with
Lower step:
1) control HIFU energy converter generates driving pulse and acts on the point of the excitation 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；
In step 1), adjust vibration measurement with laser equipment and HIFU energy converter it is confocal, make the vibration measurement with laser device signal direction of the launch and
The HIFU transducer ultrasonic direction of the launch is in the same direction；Synchronization signal control system generates synchronization signal and respectively drives arbitrary waveform generator
And data collecting card, the pulse train that arbitrary waveform generator generates are acted on by power amplifier and impedance matching network
HIFU energy converter, control HIFU energy converter generate driving pulse and act on the excitation point in tested tissue inside；Data collecting card
The signal of vibration measurement with laser equipment output is acquired；
The pulse train that arbitrary waveform generator generates is long pulse signal；
Step 1) specifically includes the following steps:
Step S1, by HIFU transducers focus at a certain position in tested tissue inside, which is referred to as excitation point, HIFU energy converter
It generates pulse and acts on point induction generation shearing wave signal；
Step S2, excitation point and shearing direction of wave travel are carried out up to the vibration of few two test points using vibration measurement with laser equipment
It monitors and records its oscillating curve；
The method for calculating shear wave velocity and frequency specifically:
Mark point position is fixed and spacing is it is known that indicate the shear wave velocity between two o'clock with the delay of same-phase between two o'clock,
Excitation point is denoted 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
Encouraging a little is Δ d with the distance between test point, and the time interval that two o'clock reaches same-phase is Δ T1, then shearing wave is in test point
Speed are as follows:
Time difference in shearing wave oscillating curve between the minimum and maximum two o'clock of relaxation amplitude is considered π/ω, therefore can
Calculate its frequency；Remember that this two o'clock time difference is Δ T2, then shearing wave vibration frequency be
In the case where obtaining local organization shear wave velocity and frequency distribution, according to the Voigt model of biological tissue and its derivation
Formula carries out viscoplasticity fitting；
Method for carrying out viscoplasticity fitting are as follows:
Relationship between plane shear velocity of wave and tissue viscoelasticity are as follows:
In formula: ρ --- density/kgm of medium-3；Angular frequency/rads of ω --- vibration-1；μ1--- shear viscosity system
Number/Pas；μ2--- the coefficient of rigidity/kPa.
2. the vibration measurement with laser monitoring HIFU according to claim 1 based on acoustic radiation force impulse response damages viscoplasticity side
Method, which is characterized in that the tested tissue is transparent prosthese or transparent in vitro tissue.
3. the vibration measurement with laser monitoring HIFU according to claim 1 based on acoustic radiation force impulse response damages viscoplasticity side
Method, which is characterized in that vibration measurement with laser equipment is laser doppler vibrometer.
4. the vibration measurement with laser monitoring HIFU according to claim 1 based on acoustic radiation force impulse response damages viscoplasticity side
Method, which is characterized in that tested tissue and vibration measurement with laser equipment are respectively connected with a three-dimensional moving device.
5. the vibration measurement with laser monitoring HIFU according to claim 1 based on acoustic radiation force impulse response damages viscoplasticity side
Method, which is characterized in that the HIFU energy converter for being 1MHz for centre frequency, pulse width is between 50 μ s-500 μ s.
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 CN106840362A (en) | 2017-06-13 |
| CN106840362B true 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) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| CN110954207B (en) * | 2019-11-15 | 2021-12-10 | 重庆医科大学 | Detection device and detection method for focused ultrasound focus sound wave structure |
| CN111307269B (en) * | 2020-03-11 | 2021-01-01 | 北京理工大学 | Laser confocal/differential confocal Raman spectrum vibration parameter measuring method |
| CN114173670A (en) * | 2020-04-24 | 2022-03-11 | 深圳迈瑞生物医疗电子股份有限公司 | Viscoelasticity measuring method and ultrasonic measuring system |
| 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 |
| CN115753502B (en) * | 2022-11-14 | 2023-08-18 | 西安交通大学 | Device and method for testing micro-nano rheological characteristics of biological tissue |
Citations (5)
| 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 |
-
2017
- 2017-03-20 CN CN201710166704.XA patent/CN106840362B/en not_active Expired - Fee Related
Patent Citations (5)
| 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)
| Title |
|---|
| 基于正交频率脉冲编码激励的软组织粘弹性定量测量;林浩铭 等;《中国生物医学工程学报》;20140831;第33卷(第4期);第392-401页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106840362A (en) | 2017-06-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106840362B (en) | Vibration measurement with laser monitoring HIFU based on acoustic radiation force impulse response damages viscoplasticity method | |
| CN110248606B (en) | Air pocket positioning | |
| CN101784234B (en) | Method and device for measuring a mean value of visco-elasticity of a region of interest | |
| CN107550458B (en) | Biological tissue multi-characteristic imaging method based on acoustoelectric effect and acoustic radiation force | |
| US10332250B2 (en) | Three-dimensional cavitation quantitative imaging method for microsecond-resolution cavitation spatial-temporal distribution | |
| BRPI0313214B1 (en) | DEVICE AND PROCEDURE FOR MEASURING THE ELASTICITY OF A HUMAN OR ANIMAL ORGAN | |
| Zhao et al. | Shear wave speed measurement using an unfocused ultrasound beam | |
| Maleke et al. | Single-element focused ultrasound transducer method for harmonic motion imaging | |
| CN109745077A (en) | Based on the elastic characteristic detection method for focusing ultrasonic sound and vibration signal | |
| CN109077754A (en) | A kind of method and apparatus measuring organization mechanics characterisitic parameter | |
| JP7295129B2 (en) | Hybrid elastography method, probes and devices for hybrid elastography | |
| CN106232013B (en) | Ultrasonic method and apparatus for characterizing a weakly anisotropic soft medium, and ultrasonic probe assembly for use in the characterizing apparatus | |
| 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 | |
| CN109717904A (en) | Elastogram system | |
| JP2017164559A (en) | Ultrasonic device | |
| CN111970972A (en) | Acoustic field mapped with ultrasonic particle velocity estimator | |
| CN110547825B (en) | Ultrasonic elastography technology based on multi-frequency shear wave amplitude analysis | |
| KR101411099B1 (en) | Method for generating mechanical waves by creating an interfacial acoustic radiation force | |
| CN105078412B (en) | Tissue elasticity analysis method and device based on optoacoustic spectrum analysis | |
| JP5435455B2 (en) | Focused vibration exciter | |
| KR102608323B1 (en) | Methods, systems and computer program products for augmentative shear wave ultrasound imaging | |
| CN106510766A (en) | Cartilage tissue elasticity measuring device and method thereof based on shear wave propagation | |
| Kim | Principles of ultrasound and applied ultrasound physics relevant for advanced sonographers | |
| 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 |