CN106618496A - All-optical photoacoustic Doppler transverse flow speed measuring method and device - Google Patents
All-optical photoacoustic Doppler transverse flow speed measuring method and device Download PDFInfo
- Publication number
- CN106618496A CN106618496A CN201611127381.5A CN201611127381A CN106618496A CN 106618496 A CN106618496 A CN 106618496A CN 201611127381 A CN201611127381 A CN 201611127381A CN 106618496 A CN106618496 A CN 106618496A
- Authority
- CN
- China
- Prior art keywords
- doppler
- photoacoustic signal
- laser
- optoacoustic
- signal
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000001228 spectrum Methods 0.000 claims abstract description 35
- 238000005259 measurement Methods 0.000 claims abstract description 25
- 239000000835 fiber Substances 0.000 claims description 25
- 230000003287 optical effect Effects 0.000 claims description 24
- 230000005622 photoelectricity Effects 0.000 claims description 20
- 238000001514 detection method Methods 0.000 claims description 15
- 230000003111 delayed effect Effects 0.000 claims description 13
- 239000002245 particle Substances 0.000 claims description 12
- 241000931526 Acer campestre Species 0.000 claims description 6
- 239000013307 optical fiber Substances 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000009738 saturating Methods 0.000 claims 1
- 238000010408 sweeping Methods 0.000 claims 1
- 230000017531 blood circulation Effects 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 6
- 230000033001 locomotion Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 210000001519 tissue Anatomy 0.000 description 5
- 239000007788 liquid Substances 0.000 description 3
- 239000006096 absorbing agent Substances 0.000 description 2
- 210000001367 artery Anatomy 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 210000003743 erythrocyte Anatomy 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000004087 circulation Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 210000002615 epidermis Anatomy 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 230000007102 metabolic function Effects 0.000 description 1
- 230000004089 microcirculation Effects 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000010186 staining Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 210000005239 tubule Anatomy 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0093—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
- A61B5/0095—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy by applying light and detecting acoustic waves, i.e. photoacoustic measurements
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Acoustics & Sound (AREA)
- Hematology (AREA)
- Cardiology (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The invention discloses an all-optical photoacoustic Doppler transverse flow speed measuring method and device. The method comprises the following steps: radiating laser to a measured object to generate a photoacoustic signal; detecting the photoacoustic signal of the measured object by a laser beam through a michelson interferometer; converting interference signals detected by laser into electric signals by using a photoelectric balance detector, acquiring the converted electric signals by using a digital acquiring card, and storing the acquired photoacoustic signal; and carrying out hilbert transform on the photoacoustic signal acquired by an scanner A, defining a complex function, calculating the photoacoustic signal as Doppler frequency spectrum standard difference, establishing a measurement model of blood flow velocity to obtain the relation between the Doppler frequency spectrum standard difference and a Doppler frequency spectrum band spread, and calculating the speed of a flowing object in a measurement region. The problems of low reflection, deep layer and low flow velocity are solved, and furthermore, the problem that the Doppler frequency shift cannot be measured when a Doppler angle is close to 90 degrees is further solved.
Description
Technical field
The present invention relates to a kind of flow-speed measurement method and device, especially a kind of full optics optoacoustic Doppler transverse flow speed are surveyed
Amount method and device, belongs to optoacoustic flow-speed measurement field.
Background technology
When the relative position of wave source and observer changes, observer receives the frequency of ripple and also can change,
This phenomenon is called Doppler effect.From after 19th century Christian Doppler discovery Doppler effect, Doppler technology
Paid close attention to by the lasting of people all the time, be successively found that ultrasonic Doppler technique, laser doppler technique and optoacoustic are more
General Le technology, is widely used in the field of testing the speed.And optoacoustic Doppler technology is used as a kind of new means that test the speed, light
Sound Doppler measurement flow velocity makes up the defect of the ultrasonic Doppler deficiency that fathoms insensitive to low flow velocity and laser-Doppler.
Optoacoustic Doppler technology is developed on the basis of optoacoustic conversion, the luminous energy of modulation when material absorbing, because what is absorbed
Luminous energy causes absorber local moment to produce temperature rise, so as to cause medium because of the generation ultrasonic signal that expands with heat and contract with cold, if absorber
Be motion, due to Doppler effect, the ultrasonic signal for detecting also can occurrence frequency movement, by the signal to collecting
The computing of Doppler frequency shift, and then the movement velocity of the body that is absorbed, the technology has had the measurement depth that ultrasonic Doppler tests the speed concurrently
The advantage of degree and LDV accuracy.The measurement of optoacoustic Doppler blood flow velocity provides abundant physiology, pathology
The information such as information and metabolic function, has broad application prospects in biomedical sector.
Ultrasonic Doppler velocity measuring technique is larger compared with deep and speed because fathoming, and the imaging of pipe arteries and veins is mainly used at present, special
It is not the measurement of this deep layer high speed of arteries, and to the low speed flow velocity ultrasonic Doppler velocity measuring technique in capillary just not
It is sensitive.And laser Doppler vibration relies primarily on the scattering of light realizing what is tested the speed, because during light transmission biological tissue
Absorbed by biological tissue, luminous energy is rapidly reduced, so the height that laser Doppler vibration is mainly used in epidermis is anti-
The measurement of blood flow in the biological tissue penetrated.When measuring low reflection, microcirculation blood flow velocity at deep tissues, ultrasonic Doppler is surveyed
Fast technology and laser Doppler vibration all can not be realized well, and should include abundant erythrocyte in blood, can
Light acoustical Doppler effect is produced with the absorption by the erythrocyte to moving to light to realize low reflection, it is micro- at deep tissues
The measurement of circulation VPV.
The content of the invention
The invention aims to solve the defect of above-mentioned prior art, there is provided a kind of full optics optoacoustic Doppler is horizontal
To flow-speed measurement method.
Another object of the present invention is to provide a kind of full optics optoacoustic Doppler transverse flow speed for realizing said method survey
Amount device.
The purpose of the present invention can be reached by adopting the following technical scheme that：
Full optics optoacoustic Doppler's transverse flow speed measuring method, the method comprising the steps of：
S1, laser is irradiated on testee produces photoacoustic signal；
S2, the photoacoustic signal that laser beam is detected testee by Michelson's interferometer；
S3, using photoelectricity balanced detector by laser acquisition to interference signal change into electric signal, then using numeral
Capture card is acquired to the electric signal for changing into, and the photoacoustic signal for collecting is stored；
S4, the photoacoustic signal to A sweep collection carry out Hilbert transform, then define complex function；
S5, the photoacoustic signal to A sweep collection do the calculating of Doppler frequency spectrum standard deviation；
S6, the measurement model for setting up VPV, obtain Doppler frequency spectrum standard deviation and Doppler frequency spectrum band spread it
Between relation, calculate the size of scattering particles flowing velocity；
S7, the photoacoustic signal repeat step S4～S6 to each A sweep collection, calculate flowing object speed in measured zone
Degree.
Further, step S4 includes：
S41, the photoacoustic signal to A sweep collection carry out Hilbert transform, such as following formula：
Wherein, H represents Hilbert transform, and * represents convolution, and p (t) represents the photoacoustic signal of A sweep collection,Represent
Signal after Hilbert transform, p (t) withIt is orthogonal；
S42, such as definition complex function, following formula：
Further, in step S5, the calculating of the Doppler frequency spectrum standard deviation, such as following formula：
Wherein, S (ω) represents Doppler power spectra,Average Doppler angular frequency is represented,Represent that jth time A sweep is adopted
The complex function of the photoacoustic signal of collection, T represents the time interval of adjacent A sweep.
Further, step S6 includes：
S61, Doppler frequency spectrum band spread are the limit Doppler frequency shifts caused by two rim rays of detection light beam
Difference determine that the corresponding average Doppler frequency shift of two edge optical fibers is：
S62, simultaneous formula (4) and formula (5), extrapolate Doppler frequency spectrum band spread, such as following formula：
Wherein, BdFor Doppler frequency spectrum band spread, f0For the centre frequency of the photoacoustic signal of A sweep collection, c is ultrasonic
Ripple spread speed in media as well, V is scattering particles flowing velocity, and θ is doppler angle,For effective numerical aperture angle, NAeff
For the effective numerical aperture of lens；
S63, simultaneous formula (3) and formula (6), according toThe size of scattering particles flowing velocity V is calculated, it is as follows
Formula：
Further, the optical maser wavelength in step S1 is 532nm.
Further, optical maser wavelength 1310nm in step S2.
Another object of the present invention can be reached by adopting the following technical scheme that：
The full optics optoacoustic Doppler transverse flow speed measurement apparatus of said method are realized, described device is visited including photoacoustic signal
Survey source, optical circulator, fiber coupler, delayed sweep line, collimating mirror, photoacoustic signal excitaton source, dichroscope, the pole of photoelectricity two
Pipe, photoelectricity balanced detector, Data Acquisition Card and computer；
The photoacoustic signal excitaton source adopts first laser device, the photoacoustic signal excitaton source, dichroscope and collimating mirror
It is sequentially connected, the laser that first laser device sends is irradiated on testee by dichroscope and collimating mirror and produces optoacoustic letter
Number；
The optical circulator, fiber coupler, delayed sweep line and collimating mirror constitute Michelson's interferometer, the optoacoustic
Signal detection source adopts second laser, and the photoacoustic signal detection source, optical circulator and fiber coupler are sequentially connected, described
Fiber coupler is connected respectively with delayed sweep line and collimating mirror, and the laser beam that second laser sends passes sequentially through light and goes in ring
Device, fiber coupler be divided into two-beam：A branch of being radiated on testee by collimating mirror is used as detection light, and another beam passes through
Delayed sweep line produces interference as reference light, the light that two beams are reflected in fiber coupler, detects the optoacoustic of testee
Signal；
The optical circulator and fiber coupler are connected respectively with photoelectricity balanced detector, the dichroscope, photoelectricity two
Pole pipe and Data Acquisition Card are sequentially connected；
The photoelectricity balanced detector, Data Acquisition Card and computer are sequentially connected, and photoelectricity balanced detector visits laser
The interference signal for measuring changes into electric signal, and then Data Acquisition Card is acquired to the electric signal for changing into, by what is collected
Photoacoustic signal storage is to computer.
Further, there are the MATLAB programs of LABVIEW data acquisition control platforms and image reconstruction in the computer.
The present invention has following beneficial effect relative to prior art：
1st, the photoacoustic signal that the present invention first collects Data Acquisition Card carries out Hilbert transform, and real number signal is opened up
Open up as complex analytic signal (defining complex function), then carry out the calculating of Doppler frequency spectrum standard deviation；Then simultaneous Doppler frequency
Spectrum standard deviation formula and Doppler frequency spectrum band spread formula, by the speed connection of the intensity of the photoacoustic signal of collection and material flows
System gets up, and the size of speed is obtained by the relevance between the photoacoustic signal that adjacent A sweep is gathered, can be with by the method
Solve the problems, such as that doppler angle cannot detect speed at 90 °, by relevance between the photoacoustic signal that adjacent A sweep is gathered
Carry out calculating speed, rather than photoacoustic signal frequency displacement is asked come calculating speed, the method for simplifying measurement with locking mutually.
2nd, the present invention does not rely on ultrasonic probe, but using the complete optical method detection ultrasound letter of michelson interferometer
Number, realize the measurement of flow velocity.
Description of the drawings
Fig. 1 is the structured flowchart of the full optics optoacoustic Doppler transverse flow speed measurement apparatus of the embodiment of the present invention 1.
Fig. 2 is the flow chart of the full optics optoacoustic Doppler transverse flow speed measuring method of the embodiment of the present invention 1.
Fig. 3 a～Fig. 3 c are the full optics optoacoustic Doppler transverse flow speed measuring method of the embodiment of the present invention 2 in fluid simulation
Measurement result figure in test.
Wherein, 1- photoacoustic signals detection source, 2- optical circulators, 3- fiber couplers, 4- delayed sweep lines, 5- collimating mirrors,
6- photoacoustic signal excitaton sources, 7- dichroscopes, 8- testees, 9- photodiodes, 10- photoelectricity balanced detectors, 11- is digital
Capture card, 12- computers.
Specific embodiment
With reference to embodiment and accompanying drawing, the present invention is described in further detail, but embodiments of the present invention are not limited
In this.
Embodiment 1：
As shown in figure 1, full optics optoacoustic Doppler's transverse flow speed measurement apparatus of the present embodiment are detected including photoacoustic signal
Source 1, optical circulator 2, fiber coupler 3, delayed sweep line 4, collimating mirror 5, photoacoustic signal excitaton source 6, dichroscope 7, photoelectricity
Diode 9, photoelectricity balanced detector 10, Data Acquisition Card 11 and computer 12.
The photoacoustic signal excitaton source 6 adopts first laser device, the photoacoustic signal excitaton source 6, dichroscope 7 and collimation
Mirror 5 is sequentially connected, and the laser that first laser device sends is irradiated on testee 8 and produces light by dichroscope 7 and collimating mirror 5
Acoustical signal；Wherein, the optical maser wavelength that first laser device sends is 532nm；
The optical circulator 2, fiber coupler 3, delayed sweep line 4 and collimating mirror 5 constitute Michelson's interferometer, described
Photoacoustic signal detection source 1 adopts second laser, and the photoacoustic signal detection source 1, optical circulator 2 and fiber coupler 3 are successively
Connection, the fiber coupler 3 is connected respectively with delayed sweep line 4 and collimating mirror 5, the laser beam that second laser sends according to
It is secondary that two-beam is divided into by optical circulator 2, fiber coupler 3：A branch of being radiated on testee 8 by collimating mirror 5 is used as spy
Light-metering, another beam is used as reference light by delayed sweep line 4, and the light that two beams are reflected produces interference in fiber coupler 3, visits
Survey the photoacoustic signal of testee 8；Wherein, the optical maser wavelength that second laser sends is 1310nm, and the fiber coupler 3 is adopted
Use 2*2 fiber couplers；
The optical circulator 2 and fiber coupler 3 are connected respectively with photoelectricity balanced detector 10；The dichroscope 7, light
Electric diode 9 and Data Acquisition Card 11 are sequentially connected, and a normalized (energy of laser output is unstable) are done, as energy
Amount correction；
The photoelectricity balanced detector 10, Data Acquisition Card 11 and computer 12 are sequentially connected, photoelectricity balanced detector 10
By laser acquisition to interference signal change into electric signal, then 11 pairs of electric signals for changing into of Data Acquisition Card are acquired,
By the photoacoustic signal for collecting storage to computer 12；Wherein, the Data Acquisition Card 11 adopts NI (National
Instruments, Nat Instr Corp.) Data Acquisition Card, have LABVIEW data acquisitions in the computer 12
The MATLAB programs of control platform and image reconstruction.
As shown in Fig. 2 the present embodiment additionally provides a kind of full optics optoacoustic Doppler transverse flow speed measuring method, the method
Realize on the basis of said apparatus, comprise the following steps：
S1, the laser by wavelength for 532nm are irradiated on testee and produce photoacoustic signal；
S2, the photoacoustic signal that wavelength is detected testee for 1310nm laser beams by Michelson's interferometer；
S3, using photoelectricity balanced detector by laser acquisition to interference signal change into electric signal, then using numeral
Capture card is acquired to the electric signal for changing into, and the photoacoustic signal for collecting is stored；
S4, the photoacoustic signal to A sweep collection carry out Hilbert transform, then define complex function, including：
S41, the photoacoustic signal to A sweep collection carry out Hilbert transform, such as following formula：
Wherein, H represents Hilbert transform, and * represents convolution, and p (t) represents the photoacoustic signal of A sweep collection,Represent
Signal after Hilbert transform, p (t) withIt is orthogonal；
S42, such as definition complex function, following formula：
S5, the photoacoustic signal to A sweep collection do the calculating of Doppler frequency spectrum standard deviation, such as following formula：
Wherein, S (ω) represents Doppler power spectra,Average Doppler angular frequency is represented,Represent that jth time A sweep is adopted
The complex function of the photoacoustic signal of collection, T represents the time interval of adjacent A sweep.
S6, the measurement model for setting up VPV, obtain Doppler frequency spectrum standard deviation and Doppler frequency spectrum band spread it
Between relation, calculate the size of scattering particles flowing velocity；
When the particle of motion passes through supersonic beam at an angle, in addition to producing Doppler frequency shift f, while also depositing
Effect is widened in the frequency spectrum of transition time：The back scattered signal of single scattering object is not single frequency, its frequency spectrum widened for
Bd；
Affected by video stretching effect, when the angle between wave beam and flow velocity is 90 °, the Doppler frequency shift of measurement
Be not as being generally considered to is zero, but the symmetrical horizontal Doppler frequency spectrum with certain bandwidth centered on zero frequency；
When angle of the conventional Doppler technology between wave beam and flow velocity is more than 75 °, the movement of frequency cannot be just detected, cause survey
The blind area of amount, and the angle that laterally Doppler technology can be between wave beam and flow velocity measures Doppler frequency spectrum broadening at 90 °
To realize testing the speed.
S61, Doppler frequency spectrum band spread are the limit Doppler frequency shifts caused by two rim rays of detection light beam
Difference determine that the corresponding average Doppler frequency shift of two edge optical fibers is：
S62, simultaneous formula (4) and formula (5), extrapolate Doppler frequency spectrum band spread, such as following formula：
Wherein, BdFor Doppler frequency spectrum band spread, f0For the centre frequency of the photoacoustic signal of A sweep collection, c is ultrasonic
Ripple spread speed in media as well, V is scattering particles flowing velocity, and θ is doppler angle,For effective numerical aperture angle, NAeff
For the effective numerical aperture of lens；
S63, simultaneous formula (3) and formula (6), according toThe size of scattering particles flowing velocity V is calculated, it is as follows
Formula：
S7, the photoacoustic signal repeat step S4～S6 to each A sweep collection, calculate flowing object speed in measured zone
Degree.
S8, calculate in measured zone and terminate flow process after flowing object speed.
The photoacoustic imaging system of the present embodiment does not rely on ultrasonic probe, and utilizes the complete optical side of michelson interferometer
Method detects ultrasonic signal, realizes the measurement of flow velocity.
Embodiment 2：
The present embodiment is specific application example, and the laser instrument of 532nm is radiated at tested by dichroscope and collimating mirror
Photoacoustic signal is produced on object, the laser beam that the laser instrument of 1310nm sends passes sequentially through optical circulator 2, fiber coupler point
Into two-beam：A branch of being radiated on testee by collimating mirror is used as detection light, and another beam is used as ginseng by delayed sweep line
Light is examined, the light that two beams are reflected produces interference in fiber coupler, the photoacoustic signal of testee is detected, finally by calculating
Machine obtains the size of flowing object speed in measured zone to photoacoustic signal analytical calculation.
It is below fluid simulation experiment：
Tested tubule (bore：300um, external diameter 500um) in be filled with the red staining particle (0.4um of diameter 5；Absorb
Peak wavelength：512nm), particle is controlled by syringe pump in the range of 0-3mm/s (speed at intervals of 0.5mm/s) at the uniform velocity
Flowing, with the laser illumination testee of 532nm photoacoustic signal is produced, and then detects optoacoustic with the laser instrument of 1310nm
Signal, computer is sent to by photoelectricity balanced detector by the photoacoustic signal for detecting, and carries out data processing.
As shown in Figure 3 a, it is imaged cross-section figure in pipe in pickup area, is calculated by the photoacoustic signal to gathering
Process obtains flow velocity image in region as shown in Figure 3 b, and Fig. 3 c are in the range of 0-3mm/s (speed at intervals of 0.5mm/s)
Graph of a relation between the liquid velocity that obtains of measurement and actual liquid speed, wherein scatterplot value are laboratory, straight line representation theory
Value, as can be seen that for liquid different in flow rate, the result measured using the method for the present embodiment is non-with theoretical value from image
Very close to it can be seen that the method for the present embodiment has the higher degree of accuracy.
In sum, the present invention is provided full optics optoacoustic Doppler's transverse flow speed measuring method and device, will survey first
Photoacoustic signal does Hilbert transform and obtains complex function, calculate Doppler frequency spectrum standard deviation；Then set up the measurement of flow velocity
Model, and Doppler frequency spectrum standard deviation and dopplerbroadening formula are connected, with the object speed that flows in computation and measurement region
Degree.
The above, patent preferred embodiment only of the present invention, but the protection domain of patent of the present invention is not limited to
This, any those familiar with the art in the scope disclosed in patent of the present invention, according to the skill of patent of the present invention
Art scheme and its inventive concept equivalent or change in addition, belong to the protection domain of patent of the present invention.
Claims (8)
1. full optics optoacoustic Doppler transverse flow speed measuring method, it is characterised in that：The method comprising the steps of：
S1, laser is irradiated on testee produces photoacoustic signal；
S2, the photoacoustic signal that laser beam is detected testee by Michelson's interferometer；
S3, using photoelectricity balanced detector by laser acquisition to interference signal change into electric signal, then using digital collection
Block the electric signal to changing into be acquired, the photoacoustic signal for collecting is stored；
S4, the photoacoustic signal to A sweep collection carry out Hilbert transform, then define complex function；
S5, the photoacoustic signal to A sweep collection do the calculating of Doppler frequency spectrum standard deviation；
S6, the measurement model for setting up VPV, obtain between Doppler frequency spectrum standard deviation and Doppler frequency spectrum band spread
Relation, calculates the size of scattering particles flowing velocity；
S7, the photoacoustic signal repeat step S4～S6 to each A sweep collection, calculate flowing object speed in measured zone.
2. full optics optoacoustic Doppler transverse flow speed measuring method according to claim 1, it is characterised in that：Step S4 bag
Include：
S41, the photoacoustic signal to A sweep collection carry out Hilbert transform, such as following formula：
Wherein, H represents Hilbert transform, and * represents convolution, and p (t) represents the photoacoustic signal of A sweep collection,Represent Xi Er
Bert conversion after signal, p (t) withIt is orthogonal；
S42, such as definition complex function, following formula：
3. full optics optoacoustic Doppler transverse flow speed measuring method according to claim 1, it is characterised in that：Step S5
In, the calculating of the Doppler frequency spectrum standard deviation, such as following formula：
Wherein, S (ω) represents Doppler power spectra,Average Doppler angular frequency is represented,Represent jth time A sweep collection
The complex function of photoacoustic signal, T represents the time interval of adjacent A sweep.
4. full optics optoacoustic Doppler transverse flow speed measuring method according to claim 3, it is characterised in that：Step S6 bag
Include：
S61, Doppler frequency spectrum band spread are the differences of the limit Doppler frequency shift caused by two rim rays of detection light beam
Determine, the corresponding average Doppler frequency shift of two edge optical fibers is：
S62, simultaneous formula (4) and formula (5), extrapolate Doppler frequency spectrum band spread, such as following formula：
Wherein, BdFor Doppler frequency spectrum band spread, f0The centre frequency of the photoacoustic signal gathered for A sweep, c exists for ultrasonic wave
Spread speed in medium, V is scattering particles flowing velocity, and θ is doppler angle,For effective numerical aperture angle, NAeffFor saturating
The effective numerical aperture of mirror；
S63, simultaneous formula (3) and formula (6), according toCalculate the size of scattering particles flowing velocity V, such as following formula：
5. the full optics optoacoustic Doppler transverse flow speed measuring method according to any one of claim 1-4, it is characterised in that：
Optical maser wavelength in step S1 is 532nm.
6. the full optics optoacoustic Doppler transverse flow speed measuring method according to any one of claim 1-4, it is characterised in that：
Optical maser wavelength 1310nm in step S2.
7. the full optics optoacoustic Doppler transverse flow speed measurement apparatus of claim 1-6 any one methods described, its feature are realized
It is：Described device includes photoacoustic signal detection source, optical circulator, fiber coupler, delayed sweep line, collimating mirror, optoacoustic letter
Number excitaton source, dichroscope, photodiode, photoelectricity balanced detector, Data Acquisition Card and computer；
The photoacoustic signal excitaton source adopts first laser device, and the photoacoustic signal excitaton source, dichroscope and collimating mirror are successively
Connection, the laser that first laser device sends is irradiated on testee by dichroscope and collimating mirror and produces photoacoustic signal；
The optical circulator, fiber coupler, delayed sweep line and collimating mirror constitute Michelson's interferometer, the photoacoustic signal
Detection source adopts second laser, the photoacoustic signal detection source, optical circulator and fiber coupler to be sequentially connected, the optical fiber
Coupler is connected respectively with delayed sweep line and collimating mirror, the laser beam that second laser sends pass sequentially through optical circulator,
Fiber coupler be divided into two-beam：A branch of being radiated on testee by collimating mirror is used as detection light, and another beam is by sweeping
Delay line is retouched as reference light, the light that two beams are reflected produces interference in fiber coupler, detect the optoacoustic letter of testee
Number；
The optical circulator and fiber coupler are connected respectively with photoelectricity balanced detector, the dichroscope, photodiode
It is sequentially connected with Data Acquisition Card；
The photoelectricity balanced detector, Data Acquisition Card and computer are sequentially connected, and photoelectricity balanced detector arrives laser acquisition
Interference signal change into electric signal, then Data Acquisition Card is acquired to the electric signal for changing into, by the optoacoustic for collecting
Signal storage is to computer.
8. described full optics optoacoustic Doppler transverse flow speed measurement apparatus according to claim 7, it is characterised in that：Institute
Stating has the MATLAB programs of LABVIEW data acquisition control platforms and image reconstruction in computer.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611127381.5A CN106618496A (en) | 2016-12-09 | 2016-12-09 | All-optical photoacoustic Doppler transverse flow speed measuring method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611127381.5A CN106618496A (en) | 2016-12-09 | 2016-12-09 | All-optical photoacoustic Doppler transverse flow speed measuring method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106618496A true CN106618496A (en) | 2017-05-10 |
Family
ID=58825581
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611127381.5A Pending CN106618496A (en) | 2016-12-09 | 2016-12-09 | All-optical photoacoustic Doppler transverse flow speed measuring method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106618496A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108444896A (en) * | 2018-01-10 | 2018-08-24 | 深圳先进技术研究院 | A kind of optoacoustic detection device |
| CN109975578A (en) * | 2019-04-09 | 2019-07-05 | 武汉新烽光电股份有限公司 | Runoff velocity method for real-time measurement, system and control device |
| CN111110193A (en) * | 2019-12-27 | 2020-05-08 | 华南师范大学 | Device and method for evaluating pipe diameter narrow fraction based on photoacoustic flow velocity measurement |
| CN111207820A (en) * | 2020-01-09 | 2020-05-29 | 哈尔滨工程大学 | Method for calibrating array elements of buoy hydrophone array in reverberation pool |
| CN111351959A (en) * | 2018-12-20 | 2020-06-30 | 罗伯特·博世有限公司 | Method and device for determining the absolute value of the flow velocity of a particulate conveying medium |
| CN114152775A (en) * | 2021-11-05 | 2022-03-08 | 安徽大学 | Deep sea hydrothermal and cold spring nozzle flow velocity measuring instrument based on optical Doppler and use method of instrument |
| CN115128299A (en) * | 2022-08-31 | 2022-09-30 | 之江实验室 | Photoacoustic particle image speed measurement system and method for measuring non-transparent flow field |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103018171A (en) * | 2012-11-29 | 2013-04-03 | 华南师范大学 | Wide-frequency-band optical-acoustic and fluorescent double-imaging device without energy converter and detection method thereof |
| CN103054610A (en) * | 2012-11-29 | 2013-04-24 | 华南师范大学 | Photoacoustic imaging device free of limitation of ultrasonic transducer frequency bands and detection method of photoacoustic imaging device |
| US20130151188A1 (en) * | 2008-07-06 | 2013-06-13 | Or-Nim Medical Ltd. | Method and system for non-invasively monitoring fluid flow in a subject |
| CN105411627A (en) * | 2015-12-28 | 2016-03-23 | 武汉大学 | Device and method for measuring transverse velocity of blood streams |
| WO2016125157A1 (en) * | 2015-02-05 | 2016-08-11 | Or-Nim Medical Ltd. | Method and system for noninvasively monitoring conditions of a subject |
-
2016
- 2016-12-09 CN CN201611127381.5A patent/CN106618496A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130151188A1 (en) * | 2008-07-06 | 2013-06-13 | Or-Nim Medical Ltd. | Method and system for non-invasively monitoring fluid flow in a subject |
| CN103018171A (en) * | 2012-11-29 | 2013-04-03 | 华南师范大学 | Wide-frequency-band optical-acoustic and fluorescent double-imaging device without energy converter and detection method thereof |
| CN103054610A (en) * | 2012-11-29 | 2013-04-24 | 华南师范大学 | Photoacoustic imaging device free of limitation of ultrasonic transducer frequency bands and detection method of photoacoustic imaging device |
| WO2016125157A1 (en) * | 2015-02-05 | 2016-08-11 | Or-Nim Medical Ltd. | Method and system for noninvasively monitoring conditions of a subject |
| CN105411627A (en) * | 2015-12-28 | 2016-03-23 | 武汉大学 | Device and method for measuring transverse velocity of blood streams |
Non-Patent Citations (3)
| Title |
|---|
| HONGWU REN ET AL: "Imaging and quantifying transverse flow velocity with the Doppler bandwidth in a phase-resolved functional optical coherence tomography", 《OPTICS LETTERS》 * |
| ZHONGJIANG CHEN ET AL: "Optically integrated trimodality imaging system: combined all-optical photoacoustic microscopy, optical coherence tomography, and fluorescence imaging", 《OPTICS LETTERS》 * |
| 卢涛等: "基于多普勒频移和频带展宽的自相关光声流速矢量测量", 《激光生物学报》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108444896A (en) * | 2018-01-10 | 2018-08-24 | 深圳先进技术研究院 | A kind of optoacoustic detection device |
| CN111351959A (en) * | 2018-12-20 | 2020-06-30 | 罗伯特·博世有限公司 | Method and device for determining the absolute value of the flow velocity of a particulate conveying medium |
| CN109975578A (en) * | 2019-04-09 | 2019-07-05 | 武汉新烽光电股份有限公司 | Runoff velocity method for real-time measurement, system and control device |
| CN109975578B (en) * | 2019-04-09 | 2021-12-31 | 武汉新烽光电股份有限公司 | Runoff flow velocity real-time measuring method, system and control device |
| CN111110193A (en) * | 2019-12-27 | 2020-05-08 | 华南师范大学 | Device and method for evaluating pipe diameter narrow fraction based on photoacoustic flow velocity measurement |
| CN111207820A (en) * | 2020-01-09 | 2020-05-29 | 哈尔滨工程大学 | Method for calibrating array elements of buoy hydrophone array in reverberation pool |
| CN111207820B (en) * | 2020-01-09 | 2021-10-26 | 哈尔滨工程大学 | Method for calibrating array elements of buoy hydrophone array in reverberation pool |
| CN114152775A (en) * | 2021-11-05 | 2022-03-08 | 安徽大学 | Deep sea hydrothermal and cold spring nozzle flow velocity measuring instrument based on optical Doppler and use method of instrument |
| CN114152775B (en) * | 2021-11-05 | 2024-01-12 | 安徽大学 | Deep sea hot liquid cold spring nozzle flow velocity measuring instrument based on optical Doppler and using method of instrument |
| CN115128299A (en) * | 2022-08-31 | 2022-09-30 | 之江实验室 | Photoacoustic particle image speed measurement system and method for measuring non-transparent flow field |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106618496A (en) | All-optical photoacoustic Doppler transverse flow speed measuring method and device | |
| CN105559756B (en) | Based on the compound microangiography method and system of total space modulation spectrum segmentation angle | |
| CN1279345C (en) | Method and systems using field-based light scattering spectroscopy | |
| JP5234470B2 (en) | Laser Doppler blood flow measurement method and apparatus | |
| JP3433508B2 (en) | Scattering absorber measurement method and scattering absorber measuring device | |
| US8263947B2 (en) | Method and device for localising fluorophores or absorbers in a surrounding medium | |
| CN103308142B (en) | A kind of speed of ultrasonic travelling wave in liquid and method and device of frequency measured | |
| Norgia et al. | Low-cost optical flowmeter with analog front-end electronics for blood extracorporeal circulators | |
| CN108801439A (en) | A kind of sound field measuring device and measurement method | |
| US5954658A (en) | Method and apparatus for measuring blood flow at precise depths in tissue and skin | |
| CN110179446A (en) | A kind of multi-modality imaging equipment for combining optoacoustic and laser speckle | |
| CN102258386A (en) | Display data obtaining apparatus and display data obtaining method | |
| CN106680159A (en) | Particle fast detection method and system based on dynamic light scattering sample ensemble analysis | |
| CN113367660B (en) | Photoacoustic Doppler flow velocity measuring device and method | |
| CN104777077A (en) | Liquid viscous coefficient measuring device and measuring method based on optical trap effect | |
| CN108535194A (en) | A kind of opto-acoustic microscopic imaging system and method based on surface plasma resonance | |
| Bauer-Marschallinger et al. | Fiber-optic annular detector array for large depth of field photoacoustic macroscopy | |
| JP5183406B2 (en) | Biological information processing apparatus and biological information processing method | |
| CN102788678B (en) | Oceanic turbulence parameter non-contact measuring method and system | |
| US4397550A (en) | Laser doppler velocimeter | |
| CN104977275B (en) | A kind of OCT image device and method under the ultrasonic field excitation of SBS | |
| CN209264563U (en) | A kind of refractive index micrometering system | |
| CN103908239A (en) | Non-contact imaging system and imaging method thereof | |
| CN206740252U (en) | Optical fiber image transmission beam velocity interferometer | |
| Savo et al. | Mean path length invariance in multiple light scattering |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170510 |
|
| RJ01 | Rejection of invention patent application after publication |