CN1257695A - Dual-ultrasonic Doppler method for measuring blood flow speed - Google Patents

Dual-ultrasonic Doppler method for measuring blood flow speed Download PDF

Info

Publication number
CN1257695A
CN1257695A CN 99127206 CN99127206A CN1257695A CN 1257695 A CN1257695 A CN 1257695A CN 99127206 CN99127206 CN 99127206 CN 99127206 A CN99127206 A CN 99127206A CN 1257695 A CN1257695 A CN 1257695A
Authority
CN
China
Prior art keywords
ultrasonic
blood flow
doppler
flow rate
array element
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
CN 99127206
Other languages
Chinese (zh)
Other versions
CN1119974C (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.)
Fudan University
Original Assignee
Fudan 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 Fudan University filed Critical Fudan University
Priority to CN 99127206 priority Critical patent/CN1119974C/en
Publication of CN1257695A publication Critical patent/CN1257695A/en
Application granted granted Critical
Publication of CN1119974C publication Critical patent/CN1119974C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

A dual-ultrasonic Doppler method for measuring blood flow speed features that two groups of transducer in Doppler probe are used to emit ultrasonic beams and after the proper included angle between ultrasonic beam and blood flow speed is chosen, the calculation equation between blood flow speed V and the deviation of blood vessel to skin is found out. In the condition of pulse Doppler effect, said transducers work in time division mode. Its advantages are simple operation and high accuracy. The real blood flow direction can be detected.

Description

Dual-ultrasonic Doppler method for measuring blood flow speed
The invention belongs to the ultrasonic technique field, is the method for a kind of pair of ultrasonic beam Doppler measurement blood flow rate.
Ultrasonic technique is used medically existing widely because of its AT characteristics.Blood flow velocity is measured and the Type B ultrasonography is that ultrasonic technique is in two broad aspect of medically using.Type B ultrasonography equipment is to the utmost for popularizing clinically, all have Doppler range rate measurement part (being called duplex system) in its high-grade image-forming instrument concurrently, can diagnose the illness in conjunction with blood flow information with morphologic information, especially at cardiovascular, the diagnostic field of diseases such as tumor.Also there is simultaneously independently the instrument of ultrasonic doppler diagnostic widely use clinically, for example diagnostic instruments such as cerebral blood flow, Umbilical Blood Flow.At present, the instrument majority of simple blood flow measurement is based on continuous wave Doppler technique.
Doppler technique also can be applicable to industrial circle except that being applied to medical science.For example, can be used for flow rate of measuring general fluid etc.
The fundamental formular of ultrasonic Doppler technique is:
Figure A9912720600031
, v is for treating velocity measurement in the formula, c is blood flow (fluid) velocity of sound, f 0Be emission supersonic frequency, f dBe detected ultrasonic doppler frequency displacement, θ is the angle between acoustic beam and the blood flow rate.Under the θ known case, can be by f dObtain speed v.But in fact the θ angle is often unpredictable, can only suppose a value, or approximate measure, so the v of gained can only be qualitatively as a result.Wang Weiqi etc. once proposed the technology (Chinese patent application number: 85100528) obtained application in continuous wave Doppler ultrasonic measurement instrument of the accurate penalty method continuous wave Doppler of alliteration bundle quantitative measurement flow velocity.But this method has adopted mimic accurate compensation method, and some shortcomings are arranged: 1. this method can only be applied to the continuous wave mode, and can not be used for pulsed mode, and what adopt in the duplexing image-forming instrument must be the pulse Doppler working method.Therefore, the scope of application has very big limitation.Must seek a best measurement position when 2. this method is measured, promptly will rotate probe and reaching a certain angle and make θ angle as far as possible little (the simulation indication is arranged), to guarantee accurate compensation.But can cause probe this moment and survey the bad bad result such as skin of maybe need pressing of position contact skin, using has quite inconvenience.3. the actual flow of blood trend can't be known.
The objective of the invention is to propose a kind of deficiency of avoiding said method, easy to use, the degree of accuracy height can record the actual flow of blood flow velocity, especially applicable to the method for two ultrasonic beam Doppler measurement blood flow rate of pulse Doppler mode.
The basic ideas of the method for two ultrasonic beam Doppler measurement blood flow rate that the present invention proposes are as follows, and two groups of transducer A, B in the ultrasonic Doppler probe are launched ultrasonic beam respectively, and frequency is designated as f respectively 01And f 02, select the two suitable ultrasonic beams and the angle of blood flow velocity, be designated as θ respectively 1, θ 2, the Doppler frequency shift f in the ultrasound echo signal that then receives D1, f D2For:
Fig. 1 measures diagram for it, and the angle of two groups of transducers is designated as 2 θ 0, 2 θ are then arranged 01+ θ 2θ 0<90 °, C is the velocity of sound of ultrasonic beam in tissue, generally gets 1540 meter per seconds.
In addition, note
Figure A9912720600042
, the calculating formula of bias angle theta that then obtains blood flow rate v and blood vessel and skin is as follows:
Figure A9912720600044
The above results has provided the size and Orientation of blood flow rate, has realized the quantitative measurement to blood flow rate.In the situation of continuous wave, concrete measuring process is as follows:
The detected frequency displacement of two pass doppler system (Jie transfers acquisition through echo) f D1And f D2Signal becomes two paths of digital signals through analog digital conversion, sends into the real-time computing of microprocessor MC or digital signal processor DSP, draws flow velocity v and bias angle theta, and real-time continuous shows or output.Also digital signal directly can be sent into the post processing that PC carries out flow velocity analysis demonstration etc.The structured flowchart of system as shown in Figure 2.
What the fundamental difference of the CN85100528 patent of the present invention and Wang Weiqi etc. was that the latter's analog compensation obtains is the approximation of flow velocity, and there is not angle information, and the inconvenience on using, inaccurate, and the present invention utilizes the Direct Digitalization of signal, can obtain flow speed value and angle message accurately with digital technology, and it is easy to use, unrestricted.
In the pulse Doppler situation, measuring process of the present invention is as follows: utilize dipulse Doppler acoustic beam to realize and the same quantitative measurement of CW (continuous wave) mode, substantially a difficulty is arranged, promptly be that the phase mutual interference can take place two acoustic beams because pulse mode is to be operated in broadband state (for CW).And the CW mode is arrowband work, and the frequency in two acoustic beams differs (as difference 500KHz) slightly, can by wave filter two bundle echoes separately handle respectively, and pulse mode be not all right.
The present invention is according to the quasi-steady state characteristic of human bloodstream: promptly several milliseconds in ten several milliseconds, human bloodstream can be considered stationary process, promptly can think less than this time blood flow state (comprising direction velocity magnitude etc.) not change.The ultransonic emission simultaneously of two-beam is received, changes time-sharing work into, promptly the A bundle is ultrasonic finish reception after, carry out B again and restraint ultransonic emission and receive, both intervals should be less than several milliseconds, even difference hundred microseconds only, for example interval is: 100 microseconds~14 millisecond.In selected same area, use the echo data of this two acoustic beam, carry out computing, obtain the size and Orientation of blood flow rate, system block diagram is as shown in Figure 3.
Above-mentioned impulse type alliteration Shu Dingliang tests the speed, and its calculating process is basic identical with the continuous wave mode.Difference is that the mid frequency of wafer is identical in the pulse mode, and f is promptly arranged 01=f 02=f 0So, have:
Figure A9912720600051
The computing formula of closing v and θ at last is:
Figure A9912720600052
Further introducing the present invention below utilizes the array transducer to carry out the measurement of the blood flow rate of pulse Doppler.
A large amount of array transducers (for example linear array, protruding battle array, and ring battle array, face battle array etc.) that use carry out ultrasonography and blood flow measurement in medical ultrasonic.Utilize linear array and convex array transducer, but the present invention's quantitative measurement blood flow rate (comprising v and θ).Detailed process is as follows:
1. system is by aforesaid time-sharing format work, two bundle ultrasonic front and back emission/receptions, and with the data of same area foundation as computing, and then obtain the result by MC or DSP etc.
2. work process (is example with the linear array, protruding battle array all fours) as shown in Figure 4.(1) sets search coverage, comprise the degree of depth and the size of direction (move show with sample line position D) and sample area (demonstration of sampling volume).(2) press above-mentioned parameter, select two groups of array element A and B, and the time delay of the reception/emission of two each array element of array element group, make the sound beam direction of this array element group point to sampled targets respectively, and make two acoustic beams intersect at sample area.Allow two acoustic beams all focus on sample area simultaneously.(3), and carry out computing with the sample area echo data by calculating formula (3), (4) and obtain the result with array element group A, B timesharing emission/reception ultrasonic beam.
The selection principle that array element group A, B are two groups is:
(a) be positioned at sample line lateral symmetry zone, if depth d is less, along with the increase of d, the array element group selection can be to probe two side shiftings, move to a determining deviation after (example is to the end of popping one's head on one side) can strengthen θ gradually 0Numerical value, certain θ 0<90 °.
(b) θ in ultrasonic device 0Selection, be the transmitting-receiving time delay that changes different chips in the array element group choose, concern as shown in Figure 5 time delay.If array element width is b, one group of array element is N altogether, then should be the time delay of i array element
Figure A9912720600053
, c is the velocity of sound in tissue in the formula, generally gets 1540 meter per seconds.
(c) selection of array element group and θ 0Selection, can import according to the parameter of setting (sample line position D depth d), carry out fast to read the software register case form, need only and in advance related data put into form.Certainly also can obtain by real-time operation.
Description of drawings
Fig. 1 measures the diagram of blood flow rate for the two ultrasonic beams of the present invention.
Fig. 2 is for measuring the system block diagram of blood flow rate in the present invention of continuous wave situation.
Fig. 3 is for measuring the system block diagram of blood flow rate in the present invention of pulse Doppler situation.
Fig. 4 utilizes the array transducer to carry out the diagram of pulse Doppler blood flow velocity measurement.
Fig. 5 is different chips sending and receiving graph of a relation time delay in the array element group.
Fig. 6 is sample line and a sampling volume relation diagram in the duplex system.
Number in the figure: 1 is ultrasonic probe, and 2 is tested blood vessel, and 3 is human body skin, and 4 is sample line, and 5 is sampling volume, and A, B are two groups of transducers in the probe 1.
The present invention be directed to the quantitative measurment of VPV and propose, but also can be applicable to the pipeline inner fluid quantitatively Measure, can be applied industrial. The quantitative measurment of VPV and direction is the urgent need of clinical medicine Want. The present invention will make the measurement of supersonic blood be developed into quantitatively by qualitative, for just making a definite diagnosis of disease provides accurately Foundation. The present invention not only goes for the continuous wave situation, also applicable to pulsed mode, has overcome existing The limitation of ultrasonic Doppler technique, and, easy and simple to handle, the accuracy height.
Provide the application example of the present invention in duplexing ultrasonic system below.
In duplexing ultrasonic system, Type B imaging and Doppler range rate measurement often need carry out simultaneously, its principle is to utilize the time-sharing work of the two to realize, scanning of image for several times normally, after obtaining the data of several scanning of image lines, the several scanning of image is carried out in the array element group that jumps to sample line position correspondence in the Doppler range rate measurement mode scanning of once testing the speed again, test the speed again, go round and begin again.When constituting b mode image and Doppler signal, the data of gained show (Doppler signal is often worked with sound spectrogram mode and Doppler's voice output mode).When looking for required speed measuring point, fixedly behind the position and size of sample line position and sampling volume, can change pure Doppler's mode over to, promptly carry out pulse Doppler scanning merely and test the speed.
What duplex system was used all is arrayed ultrasonic transducer (for example linear array, protruding battle array etc.), hardware system has functions such as the grouping focusing, scan deflection of array element group, and the realization of these functions is that the emission by different array elements in the suitable arrangement array element group receives to have and realizes different time delay.System all has the choice of location indication (sample line indication) and the sampling volume indication of Doppler range rate measurement.Utilize the select progressively of array element group can realize that the sequential scanning of B ultrasonic image and doppler position select (the corresponding array element group that sample line is indicated).Sample line and sampling volume relation are as shown in Figure 6.
The present invention can utilize the existence conditions of duplexing B ultrasonic, and the steady principle of the standard of studying the human bloodstream signal that draws according to us is (in number ms to tens ms, the blood flow signal characteristic can be considered stationary process), change by systems soft ware can realize that alliteration bundle Doppler flow velocity detects, to realize the quantitative Doppler range rate measurement in the duplex system, its angle information is a new extremely useful information.
Details are as follows for the specific implementation process:
(1) in array (for example linear array, protruding battle array, phased array etc.) system, former should have deflection delay and focusing to postpone (but do not have deflection to postpone in the general linear array, hardware will have change slightly for this reason), can realize the diplex operation of b mode image and pulse Doppler flow-speed measurement.
(2) method is changed, change the working method of (back of b mode image scanning is for several times inserted a Doppler and worked) in original duplex demonstration, change the back of b mode image scanning for several times into and insert secondary Doppler work successively, and the gained signal is carried out by formula (3) and (4) carry out computing, obtain v and θ.
(3) Doppler's array element group of working was exactly place, sample line position during original duplex showed, but must change now, should select by the position of sample line and depth location and big or small this Several Parameters of sampling volume.Usually these parameters all are manually to be provided with by keying is manual, use when of the present invention, needn't change, and be by the parameter that gained is set is determined suitable θ in software 0With two groups of array element group combinations selecting Doppler to detect, these two groups is the both sides that are symmetrical in position shown in the sample line.Selection principle should be at near-field region, and these two groups of array elements can be approaching as far as possible, and θ 0Desirable less, strengthen with the degree of depth, the array element group selection can be to probe two side shiftings, move to a determining deviation after, can strengthen θ gradually 0Numerical value.Certain θ 0<90 °.θ in ultrasonic device 0Selection, essence is the time delay that changes the different chips transmitting-receiving in the array element group chosen.
(4) selection of array element group and θ 0Selection, can do input according to the parameter that manual keying is set and carry out fast, as long as in advance relevant data are put into form to read the software register case form.Can certainly obtain by real-time operation,, not influence executed in real time (for example being placed on interior calculating flyback time etc.) as long as speed is enough.
Flow rate result can show record with digital form on screen, the drift angle can be on b mode image label orientation.

Claims (3)

1. the method for two ultrasonic beam Doppler measurement blood flow rate, it is characterized in that launching two groups of transducer A, B in the ultrasonic Doppler probe ultrasonic respectively, select the two suitable ultrasonic beams and the angle of blood flow rate, then the bias angle theta of blood flow rate v and blood vessel and skin is provided by following formula:
Figure A9912720600021
Wherein:
Figure A9912720600022
Figure A9912720600023
f 01, f 02Be respectively transducer A, B ultrasonic waves transmitted frequency, f D1, f D2Be respectively the Doppler frequency shift in the ultrasound echo signal that receives, 2 θ 0Be the angle of transducer A and B, θ 0<90,2 θ 01+ θ 2, θ 1, θ 2Be respectively the angle of two bundle ultrasonic beams and blood flow rate, c is the velocity of sound of ultrasonic beam in tissue.
2. the method for according to claim 1 pair of ultrasonic beam Doppler measurement blood flow rate is characterized in that in the pulse Doppler situation, the ultrasonic time-sharing work of two-beam: A bundle is ultrasonic finish reception after, carry out B again and restraint ultransonic emission/reception, both intervals are 100 microseconds to 14 millisecond, at this moment, and f 01=f 02=f 0, and have:
Blood flow rate v and bias angle theta are provided by following formula:
Figure A9912720600025
3. the method for according to claim 2 pair of ultrasonic beam Doppler measurement blood flow rate is characterized in that described transducer utilizes the array transducer, and detailed process is as follows:
(1) sets search coverage, comprise the degree of depth and the size of direction and sample area;
(2) select two groups of array element A and B by above-mentioned parameter, and the time delay of the reception/emission of each array element of two array element groups, make the sound beam direction of this array element group point to sampled targets respectively, and make two acoustic beams intersect at sample area, allow two acoustic beams all focus on sample area simultaneously;
(3) array element A, B are launched/receive ultrasonic beam respectively, and carry out computing, obtain the result with the sample area echo data.
CN 99127206 1999-12-30 1999-12-30 Dual-ultrasonic doppler method for measuring blood flow speed Expired - Fee Related CN1119974C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 99127206 CN1119974C (en) 1999-12-30 1999-12-30 Dual-ultrasonic doppler method for measuring blood flow speed

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 99127206 CN1119974C (en) 1999-12-30 1999-12-30 Dual-ultrasonic doppler method for measuring blood flow speed

Publications (2)

Publication Number Publication Date
CN1257695A true CN1257695A (en) 2000-06-28
CN1119974C CN1119974C (en) 2003-09-03

Family

ID=5284794

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 99127206 Expired - Fee Related CN1119974C (en) 1999-12-30 1999-12-30 Dual-ultrasonic doppler method for measuring blood flow speed

Country Status (1)

Country Link
CN (1) CN1119974C (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102401896A (en) * 2011-08-26 2012-04-04 湖南湘依铁路机车电器股份有限公司 Method for calculating Doppler signal spectrum
CN101375803B (en) * 2007-08-31 2012-06-13 奥林巴斯医疗株式会社 Ultrasonic diagnosis device
CN102988092A (en) * 2012-11-02 2013-03-27 中国人民解放军总医院 Quantitative limiting device of blood flow
CN103932737A (en) * 2014-04-28 2014-07-23 刘树英 Cardiovascular blood flow velocity sensor
US9295446B2 (en) 2010-07-29 2016-03-29 Shenzhen Mindray Bio-Medical Electronics Co., Ltd Methods and systems for pulse scanning and simultaneously displaying a blood flow image and a B-mode image
CN107080558A (en) * 2017-03-27 2017-08-22 北京大学 A kind of local pulse wave velocity device and its measuring method
CN110958858A (en) * 2017-07-28 2020-04-03 皇家飞利浦有限公司 Intraluminal imaging device with multiple center frequencies
CN115517708A (en) * 2022-11-24 2022-12-27 苏州圣泽医疗科技有限公司 Method for determining blood flow velocity by using double-crystal group and Doppler blood flow detection device

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101375803B (en) * 2007-08-31 2012-06-13 奥林巴斯医疗株式会社 Ultrasonic diagnosis device
US9295446B2 (en) 2010-07-29 2016-03-29 Shenzhen Mindray Bio-Medical Electronics Co., Ltd Methods and systems for pulse scanning and simultaneously displaying a blood flow image and a B-mode image
CN102401896A (en) * 2011-08-26 2012-04-04 湖南湘依铁路机车电器股份有限公司 Method for calculating Doppler signal spectrum
CN102988092A (en) * 2012-11-02 2013-03-27 中国人民解放军总医院 Quantitative limiting device of blood flow
CN102988092B (en) * 2012-11-02 2014-10-08 中国人民解放军总医院 Quantitative limiting device of blood flow
CN103932737A (en) * 2014-04-28 2014-07-23 刘树英 Cardiovascular blood flow velocity sensor
CN107080558A (en) * 2017-03-27 2017-08-22 北京大学 A kind of local pulse wave velocity device and its measuring method
CN110958858A (en) * 2017-07-28 2020-04-03 皇家飞利浦有限公司 Intraluminal imaging device with multiple center frequencies
US11576652B2 (en) 2017-07-28 2023-02-14 Philips Image Guided Therapy Corporation Intraluminal imaging devices with multiple center frequencies
CN110958858B (en) * 2017-07-28 2023-05-05 皇家飞利浦有限公司 Intraluminal imaging device having multiple center frequencies
CN115517708A (en) * 2022-11-24 2022-12-27 苏州圣泽医疗科技有限公司 Method for determining blood flow velocity by using double-crystal group and Doppler blood flow detection device

Also Published As

Publication number Publication date
CN1119974C (en) 2003-09-03

Similar Documents

Publication Publication Date Title
US4217909A (en) Directional detection of blood velocities in an ultrasound system
US20190200961A1 (en) Motion detection using ping-based and multiple aperture doppler ultrasound
US4257278A (en) Quantitative volume blood flow measurement by an ultrasound imaging system featuring a Doppler modality
EP0008517B1 (en) Duplex ultrasonic imaging system with repetitive excitation of common transducer in doppler modality
EP0168565B1 (en) Ultrasonic measurement method and apparatus therefor
EP0952462A2 (en) Method and apparatus for improving visualization of biopsy needle in ultrasound imaging
JPH05506371A (en) Ultrasonic testing method and device for determining tubular body position and dimensions
JPS60249946A (en) Ultrasonic tissue diagnostic method and apparatus
CN107690313B (en) Ultrasound transducer array probe for shear wave imaging
KR20190088165A (en) Ultrasound probe and manufacturing method for the same
US6293914B1 (en) Ultrasonic system and method for measurement of fluid flow
CN110301938A (en) Probe and tissue elasticity detection system
JP2021522011A (en) Systems and methods for ultrasound screening
CN1119974C (en) Dual-ultrasonic doppler method for measuring blood flow speed
JP3182419B2 (en) Blood flow measurement and display device
JPH03155843A (en) Ultrasonic diagnostic device
JPS6253182B2 (en)
Teague et al. A combined ultrasonic linear array scanner and pulsed Doppler velocimeter for the estimation of blood flow in the foetus and adult abdomen-I: Technical aspects
JPH10314167A (en) Ultrasonic diagnostic device
JPH0499566A (en) Ultrasonic diagnostic device
JPH0368694B2 (en)
JPH10314171A (en) Ultrasonic diagnostic device
WO2023278445A9 (en) System and method for non-invasive determination of pressure in a biological compartment
JPS6066732A (en) Ultrasonic pulse doppler apparatus
He et al. Doppler beam steering for blood flow velocity vector imaging

Legal Events

Date Code Title Description
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C06 Publication
PB01 Publication
C14 Grant of patent or utility model
GR01 Patent grant
C19 Lapse of patent right due to non-payment of the annual fee
CF01 Termination of patent right due to non-payment of annual fee