[go: up one dir, main page]

WO2013140583A1 - Fluid evaluation device and method - Google Patents

Fluid evaluation device and method Download PDF

Info

Publication number
WO2013140583A1
WO2013140583A1 PCT/JP2012/057379 JP2012057379W WO2013140583A1 WO 2013140583 A1 WO2013140583 A1 WO 2013140583A1 JP 2012057379 W JP2012057379 W JP 2012057379W WO 2013140583 A1 WO2013140583 A1 WO 2013140583A1
Authority
WO
WIPO (PCT)
Prior art keywords
blood flow
flow rate
fluid
period
average value
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.)
Ceased
Application number
PCT/JP2012/057379
Other languages
French (fr)
Japanese (ja)
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.)
Pioneer Corp
Original Assignee
Pioneer Corp
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 Pioneer Corp filed Critical Pioneer Corp
Priority to PCT/JP2012/057379 priority Critical patent/WO2013140583A1/en
Publication of WO2013140583A1 publication Critical patent/WO2013140583A1/en
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024Measuring pulse rate or heart rate
    • A61B5/02416Measuring pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/026Measuring blood flow
    • A61B5/0261Measuring blood flow using optical means, e.g. infrared light

Definitions

  • the present invention relates to a technical field of a fluid evaluation apparatus and method for evaluating a fluid by detecting a flow rate of the fluid such as a blood flow.
  • Patent Document 1 there is a device that detects a blood flow using a laser blood flow meter (see, for example, Patent Document 1).
  • the laser blood flow meter disclosed in Patent Document 1 irradiates a living body with laser light and detects the blood flow rate of the living body based on a change in wavelength caused by Doppler shift that occurs during reflection or scattering. More specifically, when the living body is irradiated with laser light, scattered light is generated due to the flow of blood in the blood vessels of the living body (that is, movement of red blood cells that are scatterers). The frequency of the scattered light is changed by the laser Doppler action corresponding to the moving speed of the red blood cells that are the scatterers, as compared with the frequency of the original laser light. A signal component corresponding to a so-called frequency difference signal is obtained by mutual interference between the scattered light and the reflected light of the laser light from the living body. By analyzing this signal component, the blood flow of the living body is detected.
  • the blood flow generally varies greatly according to the heart beat (for example, changes in a sine curve). For this reason, when the detected blood flow is output as it is (for example, output to a display connected to the outside of the laser Doppler blood flow meter), the blood flow that changes every moment along the time axis is output. It will be. In some cases, the blood flow volume (that is, the detection result) output in this way may be difficult to utilize. Therefore, it is conceivable to try to make it easier to utilize the output blood flow volume (that is, the detection result) by calculating the time average value of the detected blood flow volume and outputting the calculated time average value.
  • the calculated time average value of blood flow may deviate from the time average value of blood flow that should be calculated (that is, ideal). This is clarified by experiments by the inventors of the present application. For example, suppose that the time average value of the blood flow volume for every predetermined fixed period is calculated. On the other hand, the pulsation cycle of blood flow can always vary according to the pulsation of the heart. For this reason, when the time average value of the blood flow rate for each predetermined fixed period is calculated, the calculated blood is caused by the difference between the pulsation cycle of the blood flow and the fixed cycle for calculating the time average value of the blood flow rate.
  • time average value of the flow rate may become larger or smaller than the time average value of blood flow that should be originally calculated (that is, ideal). It has become. Such a change (that is, error) in the time average value of the blood flow rate is not preferable from the viewpoint of highly accurate detection of the blood flow rate.
  • the above-described technical problem is not limited to the case where the time average value of the blood flow rate (blood flow rate) is calculated, but may similarly occur when the time average value of the flow rate of any fluid is calculated. .
  • the present invention has been made in view of the above-mentioned problems, for example, and an object of the present invention is to provide a fluid evaluation apparatus and method capable of suitably calculating an average value of fluid flow rates.
  • a fluid evaluation apparatus for solving the above-described problem is to detect an irradiation unit that irradiates a laser beam to a measurement target in which a fluid flows, and the laser beam irradiated to the measurement target.
  • First detection means for detecting the flow rate of the fluid based on a change in frequency resulting from a Doppler shift of the laser light, and a second detection for detecting a pulsation cycle of the flow rate based on the change in the flow rate.
  • a fluid evaluation method for solving the above-described problem includes an irradiation step of irradiating a measurement target with a fluid flowing therein to a laser beam, and detecting the laser beam irradiated to the measurement target.
  • the fluid evaluation apparatus is obtained by detecting an irradiation unit that irradiates laser light to a measurement target in which a fluid flows, and the laser light irradiated to the measurement target.
  • First detection means for detecting the flow rate of the fluid based on a frequency change resulting from the Doppler shift of the laser light; and second detection means for detecting a pulsation cycle of the flow rate based on the change in the flow rate; Calculating means for calculating an average value of the flow rate for each period synchronized with the pulsation cycle.
  • the irradiation means irradiates the measurement target with laser light.
  • a fluid flows inside the object to be measured.
  • the laser light is preferably applied to the fluid.
  • the measurement target and the fluid include living bodies such as humans and animals and blood. Or a water pipe and water are mention
  • the first detection means detects the flow rate of the fluid based on the frequency change caused by the Doppler shift of the laser light. Specifically, for example, when laser light is irradiated on the measurement target, scattered light is generated due to the flow of fluid inside the measurement target. The frequency of the scattered light is changed by the laser Doppler action corresponding to the moving speed of the fluid as compared with the frequency of the original laser light. A signal component corresponding to a so-called frequency difference signal is obtained by the mutual interference between the scattered light and the reflected light of the laser light from the object to be measured. By analyzing this signal component, the flow rate of the fluid is detected.
  • the second detection means detects the pulsation cycle of the flow rate based on the change in the flow rate detected by the first detection means.
  • the “pulsation cycle” means, for example, a flow rate change cycle that changes so as to pulsate (for example, a period from the peak value of the flow rate to the next peak value).
  • the pulsation cycle substantially coincides with the pulse cycle (typically, the cycle of one pulse).
  • the calculating means calculates an average value of the flow rate detected by the first detecting means (for example, a time average value that is an average value of the sum of the sample values of the flow rate every predetermined time). Particularly in the present embodiment, the calculation means calculates an average value of the flow rate for each period synchronized with the pulsation cycle (for example, a time average value that is an average value of the sum of the sample values of the flow rate for each period synchronized with the pulsation cycle). To do.
  • the calculation means calculates the average value of the flow rate for each period synchronized with the pulsation cycle (for example, for each period synchronized with the pulsation cycle, the sum of the sample values of the flow rate detected during the period) Calculate the average time value).
  • the calculation means calculates an average value of the flow rate in synchronization with the pulsation cycle (for example, calculates a time average value that is an average value of the sum of the sample values of the flow rate in synchronization with the pulsation cycle).
  • the pulsation cycle detected by the second detection means is 1 second.
  • the first detection means detects the flow rate once every 0.1 seconds.
  • the calculation means calculates an average value of the flow rate per second, which is a pulsation cycle.
  • the calculation means can obtain 10 flow rate samples detected during the 1 second, which is the pulsation cycle. The average value of the sum of values is calculated.
  • the pulsation cycle detected by the second detection unit is changed to 1.2 seconds, the calculation unit calculates the sum of the sample values of 12 flow rates detected during the pulsation cycle of 1.2 seconds. The average value of is calculated.
  • the average value of the fluid flow rate is calculated in accordance with the pulsation cycle of the fluid flow rate that may change due to some factor. That is, the calculation cycle for calculating the average value of the fluid flow rate is appropriately adjusted according to the pulsation cycle of the fluid flow rate. Therefore, according to the fluid evaluation device of the present embodiment, the pulsation cycle of the fluid flow rate and the fluid are compared with the fluid evaluation device of the comparative example in which the calculation cycle when calculating the average value of the fluid flow rate is fixed. The possibility of a deviation from the calculation cycle for calculating the average value of the flow rates of is reduced.
  • the fluid evaluation apparatus of the present embodiment there is little or no deviation between the pulsation cycle of the fluid flow rate and the calculation cycle for calculating the average value of the fluid flow rate. Therefore, according to the fluid evaluation apparatus of the present embodiment, the flow rate of the fluid calculated by the calculation unit as compared with the fluid evaluation apparatus of the comparative example in which the calculation cycle when calculating the average value of the flow rate of the fluid is fixed. Is less likely to deviate from the average value of the flow rate of the fluid to be originally calculated (that is, ideal). Alternatively, according to the fluid evaluation apparatus of the present embodiment, the average value of the flow rate of the fluid calculated by the calculation unit may deviate from the average value of the flow rate of the fluid that should be originally calculated (that is, ideal). Little or no. For this reason, the fluid evaluation apparatus of this embodiment can calculate the average value of the flow volume of a fluid suitably.
  • the calculation unit is configured to perform the calculation every the same period as the pulsation period or every N times (where N is an integer of 2 or more) the pulsation period. Calculate the average flow rate.
  • the calculation means preferably uses an average value of the flow rate for each period synchronized with the pulsation cycle (for example, a time average value that is an average value of the sum of the sample values of the flow rate for each period synchronized with the pulsation cycle). Can be calculated.
  • the fluid evaluation apparatus further includes an output unit that outputs the average value calculated by the calculation unit at a desired timing.
  • the average value of the flow rate of the fluid calculated by the calculation unit is output to another device (for example, a display or an analysis circuit) provided in the fluid evaluation device or attached to the outside of the fluid evaluation device. Is done.
  • another device for example, a display or an analysis circuit
  • the desired timing is a fixed timing that appears periodically.
  • the average value of the flow rate of the fluid calculated by the calculation means is output at a constant cycle. For this reason, even if the period during which the calculation unit calculates the average value (that is, the period synchronized with the pulsation cycle) fluctuates, the cycle in which the average value of the fluid flow rate calculated by the calculation unit is output from the output unit is constant. Kept. Therefore, compared with the case where the average value of the flow rate of the fluid calculated by the calculation unit is output in a random cycle, the average value can be easily used for other apparatuses that acquire the average value from the output unit.
  • the desired timing is a timing designated by the user.
  • the average value of the flow rate of the fluid calculated by the calculation means is output at a timing designated by the user. For this reason, even if the period during which the calculation unit calculates the average value (that is, the period synchronized with the pulsation cycle) fluctuates, the cycle in which the average value of the fluid flow rate calculated by the calculation unit is output from the output unit is It is adjusted appropriately according to the instructions. Therefore, compared with the case where the average value of the flow rate of the fluid calculated by the calculation means is output in a cycle unrelated to the user's instruction, for other devices that acquire the average value from the output means, the average value is It becomes easy to use.
  • the fluid evaluation device further includes third detection means for detecting the amplitude of the flow rate, and the calculation means (i) when the amplitude is equal to or greater than a predetermined amount, An average value of the flow rate for each period synchronized with the pulsation cycle detected by the detection unit is calculated, and (ii) the pulsation detected by the first detection unit in the past when the amplitude does not exceed the predetermined amount. The average value of the flow rate for each period synchronized with the cycle is calculated.
  • the calculating means replaces the pulsation cycle with relatively low reliability that the second detecting means detects in real time.
  • the second detection means having relatively high reliability calculates the average value of the flow rate for each period synchronized with the pulsation cycle detected in the past. For this reason, the reliability of the average value calculated by the calculation means is relatively compared with the case where the average value of the flow rate for each period synchronized with the pulsation cycle detected in real time by the second detection means is always calculated. Can be increased.
  • the calculation means is a period synchronized with the pulsation period itself of the flow rate detected by the second detection means in the past (that is, the same period as the pulsation period or N times the pulsation period (where N is an integer of 2 or more).
  • the average value of the flow rate for each 1 / N times may be calculated.
  • the calculation unit may calculate an average value of the flow rate for each period synchronized with the time average value of the pulsation cycle of the flow rate detected in the past by the second detection unit.
  • the fluid is a blood flow flowing inside the living body as the measurement target, and the pulsation cycle is a pulse wave cycle of the blood flow.
  • the calculation means has a blood flow rate for each period synchronized with the pulse wave cycle of blood (that is, the pulse cycle, typically one pulse cycle).
  • An average value of that is, blood flow volume
  • the fluid evaluation method of the present embodiment is obtained by irradiating a laser beam to a measurement target in which a fluid flows inside, and detecting the laser beam irradiated to the measurement target.
  • the fluid evaluation method of the present embodiment may adopt various aspects.
  • the fluid evaluation apparatus includes an irradiation unit, a first detection unit, a second detection unit, and a calculation unit.
  • the fluid evaluation method of this embodiment includes an irradiation process, a first detection process, a second detection process, and a calculation process. Therefore, the average value of the flow rate of the fluid can be suitably calculated.
  • the fluid evaluation device is applied to a blood flow rate detection device that detects a flow rate (that is, blood flow rate) of blood flowing in a blood vessel of a living body
  • a flow rate that is, blood flow rate
  • the fluid evaluation device can detect any flow rate of any fluid other than blood (for example, water flowing in a water pipe or blood flowing in a tube constituting a blood flow circuit of an artificial dialysis device).
  • the present invention may be applied to a flow rate detection device.
  • FIG. 1 is a block diagram illustrating a configuration of a blood flow rate detection device 1 according to the first embodiment.
  • the blood flow rate detection device 1 includes a laser element 11, a light receiving element 12, an amplifier 13, an analog / digital (A / D) converter 14, an arithmetic circuit 15, A pulse wave period detection circuit 16, an averaging circuit 17, a timing generator 18, and a sampling circuit 19 are provided.
  • the laser element 11 constitutes a specific example of “irradiation means” and irradiates the living body 100 such as a human being or an animal with the laser beam LB. At this time, it is preferable that the laser element 11 irradiates the blood vessel in the living body 100 with the laser beam LB.
  • the light receiving element 12 constitutes a specific example of “first detection means”, and beat signals generated by mutual interference between the reflected light of the laser light LB from the living body 100 and the scattered light of the laser light LB from the living body 100. Receives light. The light receiving element 12 generates a detection current obtained by converting the received beat signal light into an electrical signal.
  • the amplifier 13 constitutes a specific example of “first detection means”, which amplifies the detection current output from the light receiving element 12 after converting it into a voltage signal.
  • the A / D converter 14 constitutes a specific example of “first detection means”, and an A / D converter 14 outputs A / D with respect to the output of the amplifier 13 (that is, the voltage signal corresponding to the beat signal light received by the light receiving element 12). D conversion processing (that is, quantization processing) is performed. As a result, the A / D converter 14 outputs the sample value of the voltage signal (that is, the quantized voltage signal) corresponding to the beat signal light received by the light receiving element 12 to the arithmetic circuit 15.
  • the arithmetic circuit 15 constitutes a specific example of “first detection means”, and outputs the output of the A / D converter 14 (that is, the sample value of the voltage signal corresponding to the beat signal light received by the light receiving element 12).
  • frequency analysis using FFT is performed.
  • the arithmetic circuit 15 calculates a blood flow rate Q (k (where k is a variable assigned in time series in accordance with the blood flow calculation timing)).
  • the pulse wave cycle detection circuit 16 constitutes a specific example of “second detection means”, and by analyzing the temporal change of the blood flow rate Q (k) calculated by the calculation circuit 15, the pulse wave cycle (that is, , The period of blood pulsation, and substantially the period of pulse) Tr. Thereafter, the pulse wave cycle detection circuit 16 counts the number of samples n of the blood flow rate Q (k) corresponding to the pulse wave cycle Tr (that is, the number of samples of the blood flow rate Q (k) detected in the same period as the pulse wave cycle. n) is calculated.
  • the averaging circuit 17 constitutes a specific example of “calculation means”, and calculates an average value of n blood flow volumes Q (k) (hereinafter referred to as “average blood flow volume Qh” as appropriate).
  • the timing generator 18 generates a timing at which the sampling circuit 19 outputs the average blood flow rate Qh.
  • the sampling circuit 19 constitutes a specific example of “output means”, and the blood flow detection device 1 synchronizes the average blood flow Qh calculated by the averaging circuit 17 with the timing generated by the timing generator 18. (Or to a processing block (not shown) included in the blood flow rate detection device 1).
  • the blood flow detection device 1 may output the average blood flow Qh calculated by the averaging circuit 17 as it is. In this case, the blood flow rate detection device 1 may not include the timing generator 18 and the sampling circuit 19.
  • FIG. 2 is a flowchart showing the operation flow of the blood flow rate detection device 1 of the first embodiment.
  • FIG. 3 is a timing chart showing waveforms of various signals observed during the operation of the blood flow rate detection device 1 of the first embodiment.
  • the laser element 11 irradiates a living body 100 such as a human being or an animal with the laser beam LB (step S11). At this time, it is preferable that the laser element 11 irradiates the blood vessel in the living body 100 with the laser beam LB.
  • the light receiving element 12 receives the beat signal light generated by the mutual interference between the reflected light of the laser light LB from the living body 100 and the scattered light of the laser light LB from the living body 100 (step S12). Specifically, when the living body 100 is irradiated with the laser light LB, scattered light is generated due to blood flow inside the blood vessel in the living body 100 (that is, movement of red blood cells that are scatterers). The frequency of the scattered light is changed by the laser Doppler action corresponding to the moving speed of the blood as compared with the frequency of the original laser light. The light receiving element 12 receives beat signal light (so-called frequency difference signal) generated by mutual interference between the scattered light and the reflected light of the laser light LB from the living body 100.
  • beat signal light syn-called frequency difference signal
  • the scattered light that generates the beat signal light forward scattered light corresponding to the transmitted light of the laser light LB irradiated on the living body 100 may be used, or the reflected light of the laser light LB irradiated on the living body 100. Backscattered light corresponding to may be used.
  • the light receiving element 12 generates a detection current obtained by converting the received beat signal light into an electric signal.
  • the light receiving element 12 outputs the generated detection current to the amplifier 13.
  • the amplifier 13 converts the detection current output from the light receiving element 12 (that is, the detection current corresponding to the beat signal light received by the light receiving element 12) into a voltage signal and then amplifies it.
  • the amplifier outputs a voltage signal to the A / D converter 14.
  • the A / D converter 14 performs an A / D conversion process (that is, a quantization process) on the output of the amplifier 13 (that is, a voltage signal corresponding to the beat signal light received by the light receiving element 12) (steps). S13).
  • the A / D converter 14 outputs the sample value of the voltage signal (that is, the quantized voltage signal) corresponding to the beat signal light received by the light receiving element 12 to the arithmetic circuit 15.
  • the sampling period of the A / D converter 14 is Ta
  • the A / D converter 14 samples the voltage signal sample value (corresponding to the beat signal light received by the light receiving element 12 for each period Ta ( That is, a quantized voltage signal) is output.
  • the arithmetic circuit 15 performs frequency analysis using FFT (Fast Fourier Transform) on the output of the A / D converter 14 (that is, the sample value of the voltage signal corresponding to the beat signal light received by the light receiving element 12). .
  • the arithmetic circuit 15 calculates a blood flow rate Q (k) (step S14: refer to the first waveform in FIG. 3).
  • the arithmetic circuit 15 performs FFT on the sample value of the voltage signal corresponding to the beat signal light and the Nf sample values.
  • the arithmetic circuit 15 calculates the blood flow rate Q (k) using a first moment that is a result of multiplication of the power spectrum and the frequency vector obtained by performing the FFT.
  • the arithmetic circuit 15 outputs the calculated blood flow volume Q (k) to each of the pulse wave period detection circuit 16 and the averaging circuit 17.
  • the period in which the blood flow volume Q (k) is calculated (that is, the period in which one (in other words, one sample) blood flow volume Q (k) is calculated) Tq is the first and second stages in FIG. As shown in each waveform of the stage, Nf ⁇ Ta (that is, the number of FFT points Nf ⁇ the sampling period Ta of the A / D converter 14).
  • the pulse wave cycle detection circuit 16 analyzes the time change of the blood flow rate Q (k) calculated by the arithmetic circuit 15 (that is, the time change of the blood flow rate Q (k) of a plurality of samples) to thereby calculate the pulse wave cycle. Tr is calculated (step S15). For example, the pulse wave period detection circuit 16 binarizes the blood flow rate Q (k) with reference to the amplitude center of the blood flow rate Q (k) as shown in the first and third waveforms in FIG. And the pulse wave period Tr is calculated for a period from the rise to the fall of the binarized blood flow rate Q (k) (that is, a combination of a pulse with an amplitude greater than or equal to a pulse less than the amplitude center). May be.
  • the pulse wave period detection circuit 16 determines that the blood flow rate Q (k) is the amplitude center of the blood flow rate Q (k), as shown in the first and third waveforms in FIG.
  • the pulse wave period Tr may be calculated by detecting a crossing point (so-called zero cross point).
  • the pulse wave period detection circuit 16 detects the peak value of the blood flow rate Q (k) and calculates the interval between two peak values that are continuous on the time axis as the pulse wave period Tr. Also good.
  • the pulse wave cycle detection circuit 16 detects the bottom value of the blood flow rate Q (k) and calculates the interval between two bottom values that are continuous on the time axis as the pulse wave cycle Tr. Also good.
  • the pulse wave period detection circuit 16 calculates the sample number n of the blood flow rate Q (k) corresponding to the pulse wave period Tr as shown in the third waveform in FIG. 3 (step S16). Specifically, the pulse wave cycle detection circuit 16 calculates the number of samples n of the blood flow rate Q (k) detected in the same period as the pulse wave cycle Tr. In other words, the pulse wave period detection circuit 16 calculates the number of samples n of the blood flow rate Q (k) included in the same period as the pulse wave period Tr.
  • the pulse wave period detection circuit 16 may carry up the digits after the decimal point when the calculated sample number n has digits after the decimal point. Alternatively, digits after the decimal point may be carried down. Alternatively, digits after the decimal point may be rounded off.
  • the number of samples n calculated by the pulse wave period detection circuit 16 changes as appropriate according to the time variation of the blood flow rate Q (n). That is, the number of samples n calculated by the pulse wave cycle detection circuit 16 changes appropriately according to the pulse wave cycle Tr.
  • the third waveform in FIG. 3 shows an example in which the number of samples n1 is calculated in a certain period while the number of samples n2 is calculated in the subsequent period.
  • the pulse wave period detection circuit 16 outputs the calculated sample number n to the averaging circuit 17.
  • the averaging circuit 17 determines (Q (0) + Q (1) +... In a certain period (a period in which the number of samples n1 is calculated). Then, + (Q (n1-1)) / n1 is calculated as the average blood flow rate Qh, and then the averaging circuit 17 calculates (Q (n11) in another period (period in which the number of samples n2 is calculated). ) + Q (n1 + 1) +... + (Q (n1 + n2-1)) / n2 is calculated as the average blood flow Qh, and the averaging circuit 17 outputs the calculated average blood flow Qh to the sampling circuit 19. .
  • the sampling circuit 19 synchronizes the average blood flow Qh calculated by the averaging circuit 17 with the timing generated by the timing generator 18 (corresponding to the resampling timing shown in the fifth waveform in FIG. 3),
  • the data is output to the outside of the blood flow detection device 1 (or to a processing block (not shown) included in the blood flow detection device 1) (step S18).
  • the timing generated by the timing generator 18 is preferably a fixed timing (that is, a timing that appears in a fixed cycle).
  • the output rate (in other words, transfer rate) of the average blood flow Qh output from the blood flow detection device 1 is constant.
  • the sampling circuit 19 may output a time average value Qh ′ calculated by further averaging the average blood flow rate Qh every predetermined time. In this case, fluctuations in the value output from the blood flow rate detection device 1 (in this case, the time average value Qh ′) can be reduced.
  • FIG. 4 shows a blood flow detection of a comparative example in which an average value (average blood flow Qh) of blood flow Q (k) of m fixed samples is calculated without calculating the number n of samples corresponding to the pulse wave period Tr. It is a timing chart which shows the waveform of various signals observed at the time of operation of an apparatus.
  • the blood flow rate detection device of the comparative example calculates (Q (0) + Q (1) +... + (Q (m ⁇ 1)) / m as the average blood flow rate Qh.
  • the blood flow rate detection device of the comparative example adds 21 sample values of the blood flow rate Q (k) sequentially and sequentially calculated along the time axis.
  • the pulse wave cycle Tr varies. This is because the blood flow pulsates with the heartbeat, and the pulse wave period Tr naturally varies with the change of the heartbeat. Therefore, a deviation occurs between the pulse wave period Tr and the calculation period for calculating the average blood flow rate Qh (that is, the period corresponding to the fixed sample number m). Therefore, depending on the timing at which the average blood flow Qh is calculated, m blood flows Q (k) for calculating the average blood flow Qh, as shown on the relatively left side of the waveform in the first stage of FIG.
  • the number of blood flows Q (k) distributed on the side larger than the amplitude center may be larger than the number of blood flows Q (k) distributed on the side smaller than the amplitude center.
  • the average blood flow Qh calculated from the m blood flows Q (k) is ideally calculated as shown in the vicinity of the center of the fourth waveform in FIG. Is likely to become bigger than.
  • the timing of calculating the average blood flow Qh as shown on the relatively right side of the waveform in the first stage of FIG.
  • m blood flows Q (k) for calculating the average blood flow Qh the number of blood flows Q (k) distributed on the side larger than the amplitude center may be smaller than the number of blood flows Q (k) distributed on the side smaller than the amplitude center.
  • the average blood flow Qh calculated from such m blood flow volumes Q (k) is an ideal average blood to be originally calculated, as shown on the right side of the waveform in the fourth row in FIG. The possibility of becoming smaller than the flow rate Qh is high. Therefore, if the number of samples m when calculating the average blood flow Qh is fixed, the average blood flow Qh may include an error.
  • the fluctuation of the average blood flow Qh having a frequency lower than the pulse wave period Tr is a fluctuation involving the activity of the autonomic nervous system of the living body 100.
  • This fluctuation is an original fluctuation as a reaction of the living body 100 and is preferably detected.
  • fluctuations that is, fluctuations in the average blood flow Qh
  • the average blood flow Qh are generated in the average blood flow Qh output from the blood flow detection device of the comparative example, as shown in the waveform at the fourth stage in FIG. .
  • this fluctuation cycle is an erroneous fluctuation due to the above-described error in the average blood flow Qh (that is, the deviation of the blood flow Q (k) used for calculating the average blood flow Qh from the amplitude center).
  • the reaction of the living body 100 is different from the original fluctuation.
  • such an erroneous fluctuation also has a frequency (that is, a long period) lower than the pulse wave period Tr of the blood flow rate Q (k). Therefore, in the blood flow rate detection device of the comparative example, the original fluctuation as the reaction of the living body 100 and the error of the above average blood flow rate Qh (that is, the blood flow rate Q (k) used when calculating the average blood flow rate Qh). Incorrect fluctuations due to (bias with respect to the amplitude center) are mixed. For this reason, there is a possibility that the original fluctuation is not suitably detected as a reaction of the living body 100.
  • the average value (average blood flow volume Qh) of the blood flow volume Q (k) of n samples calculated in accordance with the pulse wave cycle Tr is calculated. That is, the calculation cycle (that is, the cycle corresponding to the number of samples n) when calculating the average blood flow Qh is appropriately adjusted according to the pulsation cycle Tr of the blood flow Q (k). Therefore, according to the blood flow rate detection device 1 of the present embodiment, there is a possibility that a deviation occurs between the pulse wave cycle Tr and the calculation cycle for calculating the average blood flow rate Qh, as compared with the blood flow rate detection device of the comparative example. Is reduced.
  • the blood flow detection device 1 of the present embodiment there is little or no deviation between the pulse wave cycle Tr and the calculation cycle for calculating the average blood flow Qh.
  • the n blood flows Q (k) for calculating the average blood flow Qh is larger than the amplitude center. The possibility that the number of blood flows Q (k) distributed on the side is larger than the number of blood flows Q (k) distributed on the side smaller than the amplitude center is reduced.
  • the blood flow volume Q (k) distributed on the side larger than the amplitude center among the n blood flow volumes Q (k) for calculating the average blood flow volume Qh is hardly or not at all.
  • the amplitude center of the n blood flow volumes Q (k) for calculating the average blood flow volume Qh is compared with the blood flow volume detection device of the comparative example. The possibility that the number of blood flow volumes Q (k) distributed on the larger side is smaller than the number of blood flow volumes Q (k) distributed on the smaller side than the amplitude center is reduced.
  • the average blood flow Qh that should be originally calculated (that is, the ideal) average blood flow Qh is compared with the blood flow detection device of the comparative example.
  • the possibility of shifting from is reduced.
  • the average blood flow rate Qh is hardly or completely deviated from the originally calculated (that is, ideal) average blood flow rate Qh. For this reason, the blood flow rate detection device 1 of the present embodiment can suitably calculate the average blood flow rate Qh.
  • the blood flow detection device 1 of the present embodiment compared to the blood flow detection device of the comparative example, it is used when calculating the error of the average blood flow Qh (that is, the average blood flow Qh).
  • the error of the average blood flow Qh that is, the average blood flow Qh.
  • the above-described error in the average blood flow Qh that is, the deviation of the blood flow Q (k) used for calculating the average blood flow Qh from the amplitude center.
  • the resulting false fluctuations will occur little or no in the mean blood flow Qh.
  • the original fluctuation is suitably detected as the reaction of the living body 100.
  • the blood flow rate detection device 1 of the present embodiment that calculates the number of samples n according to the pulse wave period Tr is very advantageous in practice.
  • the average blood flow Qh is output in synchronization with the timing generated by the timing generator 18. Therefore, the output rate (in other words, transfer rate) of the average blood flow Qh output from the blood flow detection device 1 can be made constant. Therefore, it is relatively easy to construct an interface between the blood flow rate detection device 1 and the external device.
  • the blood flow rate detection device 1 can output the averaged blood flow rate Qh in a cycle longer than the cycle Tq in which the blood flow rate Q (k) is calculated. For this reason, compared with the case where the blood flow rate Q (k) is output as it is, the output data amount per unit time can be reduced. Furthermore, the blood flow rate detection device 1 of the present embodiment can output an average blood flow rate Qh that is the result of averaging after the blood flow rate Q (k) is averaged. For this reason, compared with the case where the blood flow rate Q (k) is output as it is, the output data amount per unit time can be reduced.
  • FIG. 5 is a block diagram showing the configuration of the blood flow rate detection device 2 of the second embodiment.
  • the blood flow rate detection device 2 of the second embodiment does not include the timing generator 18 and the user request reception unit 21 as compared with the blood flow rate detection device 1 of the first embodiment. It differs in that it is equipped.
  • Other components included in the blood flow detection device 2 of the second embodiment may be the same as other components included in the blood flow detection device 1 of the first embodiment.
  • the user request accepting unit 21 accepts an instruction from a user (that is, a user of the blood flow detection device 2) instructing the timing at which the sampling circuit 19 outputs the average blood flow Qh.
  • the user request receiving unit 21 generates a timing at which the sampling circuit 19 outputs the average blood flow Qh in synchronization with the received user instruction.
  • FIG. 6 is a flowchart showing the operation flow of the blood flow rate detection device 2 of the second embodiment.
  • the blood flow rate detection device 2 of the second embodiment performs the operations from step S11 to step S17 in the same manner as the blood flow rate detection device 1 of the first embodiment.
  • the sampling circuit 19 generates the average blood flow rate Qh calculated by the averaging circuit 17 at the timing when the timing generator 18 generates (that is, the timing according to the user's instruction). In synchronization with (), the blood is output to the outside of the blood flow detection device 1 (step S21).
  • Such a blood flow volume detection device 2 of the second embodiment can preferably enjoy the same effects as the various effects that the blood flow volume detection device 1 of the first embodiment can enjoy.
  • the average blood flow rate Qh is output in synchronization with the timing instructed by the user. Therefore, the output rate (in other words, transfer rate) of the average blood flow Qh output from the blood flow detection device 1 can be set to a user's favorite rate. Therefore, the average blood flow Qh is output in a user-friendly manner.
  • FIG. 7 is a block diagram illustrating a configuration of the blood flow rate detection device 3 according to the third embodiment.
  • the blood flow rate detection device 3 according to the third embodiment has a pulse rate calculation circuit 31, an averaging circuit 32, and a sampling circuit 33 as compared with the blood flow rate detection device 1 according to the first embodiment. Is different in that it is further equipped.
  • Other components included in the blood flow detection device 3 of the third embodiment may be the same as other components included in the blood flow detection device 1 of the first embodiment.
  • the pulse rate calculation circuit 31 calculates the pulse rate R (i (where i is a variable assigned in time series according to the timing of calculating the pulse rate)) from the sample number n calculated by the pulse wave period detection circuit 16. Is calculated.
  • the averaging circuit 32 calculates an average value of the predetermined number w of pulse rates R (i) (hereinafter, referred to as “average pulse rate Rh” as appropriate).
  • the sampling circuit 33 synchronizes the average pulse rate Rh calculated by the averaging circuit 32 to the outside of the blood flow detection device 1 (or the blood flow detection device 1 inside) in synchronization with the timing generated by the timing generator 18. Output to a processing block (not shown).
  • the blood flow rate detection apparatus 1 may output the average pulse rate Rh calculated by the averaging circuit 32 as it is. In this case, the blood flow rate detection device 1 may not include the sampling circuit 33.
  • FIG. 8 is a flowchart showing the flow of the operation of the blood flow rate detection device 3 of the third embodiment.
  • the blood flow rate detection device 3 of the third embodiment performs the operations from step S11 to step S17 in the same manner as the blood flow rate detection device 1 of the first embodiment.
  • the pulse rate calculation circuit 31 further calculates the pulse rate R (i) from the sample number n calculated by the pulse wave cycle detection circuit 16 (step S31).
  • the pulse rate R (i) may be calculated using a mathematical formula of sampling period Ta) ⁇ 60 of the converter 14.
  • the pulse rate calculation circuit 31 outputs the calculated pulse rate R (i) to the averaging circuit 32.
  • an average value that is, average pulse rate Rh
  • the sampling circuit 33 synchronizes the average pulse rate Rh calculated by the averaging circuit 32 to the outside of the blood flow detection device 1 (or the blood flow detection device 1 is synchronized with the timing generated by the timing generator 18).
  • the data is output to a processing block (not shown) provided inside (step S33).
  • Such a blood flow volume detection device 3 of the third embodiment can preferably enjoy the same effects as the various effects that the blood flow volume detection device 1 of the first embodiment can enjoy.
  • the blood flow rate detection device 2 of the second embodiment not only the average blood flow rate Qh but also the average pulse rate Rh is output.
  • the blood flow detection device 4 of the fourth example will be described with reference to FIGS. 9 to 11.
  • the same configuration and operation as those of the blood flow rate detection device 1 of the first embodiment are denoted by the same reference numerals and step numbers, and detailed description thereof is omitted.
  • FIG. 9 is a block diagram showing the configuration of the blood flow rate detection device 4 of the fourth embodiment.
  • the blood flow volume detection device 4 of the fourth embodiment is compared with the blood flow volume detection device 1 of the first embodiment, and the pulse amplitude detection circuit 41, the level decrease detection circuit 42, and the averaging The difference is that a circuit 43 and a selection circuit 44 are further provided.
  • Other components included in the blood flow detection device 4 of the fourth embodiment may be the same as other components included in the blood flow detection device 1 of the first embodiment.
  • the pulse amplitude detection circuit 41 constitutes a specific example of “third detection means”, and analyzes the time change of the blood flow volume Q (k) calculated by the arithmetic circuit 15 to thereby detect the pulse amplitude (that is, blood (Amplitude of flow rate Q (k)) A is calculated.
  • the level drop detection circuit 42 determines whether or not the pulse amplitude A detected by the pulse amplitude detection circuit 41 is less than a predetermined threshold value.
  • the averaging circuit 43 calculates a predetermined number p of sample numbers n (j (where j is a variable assigned in time series according to the calculation timing of the number of samples)) (hereinafter referred to as “average samples” as appropriate). A number Nh ′′).
  • the selection circuit 44 selects one of the number of samples n (j) detected by the pulse wave period detection circuit 16 and the average sample value Nh calculated by the averaging circuit 43 based on the determination result in the level decrease detection circuit 42. Output to the averaging circuit 17.
  • FIG. 10 is a flowchart showing the flow of the operation of the blood flow rate detection device 4 of the fourth embodiment.
  • FIG. 11 is a timing chart showing waveforms of various signals observed during the operation of the blood flow rate detection device 4 of the fourth embodiment.
  • the blood flow rate detection device 4 of the fourth embodiment performs the operations from step S11 to step S16 in the same manner as the blood flow rate detection device 1 of the first embodiment.
  • the pulse wave amplitude detection circuit 41 analyzes the time change of the blood flow volume Q (k) calculated by the arithmetic circuit 15 to obtain the pulse amplitude (that is, the amplitude of the blood flow volume Q (k), which is shown in FIG. (Refer to the second-stage waveform) A is calculated (step S42).
  • the pulse wave amplitude detection circuit 41 may handle the amplitude of the blood flow rate Q (k) as the pulse amplitude A as it is.
  • the pulse wave amplitude detection circuit 41 may handle a value obtained by performing a predetermined calculation process on the amplitude of the blood flow volume Q (k) as the pulse amplitude A.
  • the level decrease detection circuit 42 determines whether or not the pulse amplitude A detected by the pulse wave amplitude detection circuit 41 is less than a predetermined threshold (see the second waveform in FIG. 11) (step S43).
  • the predetermined threshold value may be set to a lower limit value of the pulse amplitude A that can suitably calculate the pulse wave period Tr.
  • the predetermined threshold may be set to a lower limit value of the pulse amplitude A that can calculate a reliable pulse wave period Tr.
  • step S43 when it is determined that the pulse amplitude A is not less than the predetermined threshold value (step S43: No), it is determined that the pulse amplitude A is not relatively decreased. That is, it is determined that the pulse amplitude A has a magnitude that allows the reliable pulse wave period Tr to be calculated. Therefore, the sample value n (j) currently calculated by the pulse wave period detection circuit 16 is estimated to be highly reliable. Therefore, in this case, the selection circuit 44 averages the number of samples n (j) detected by the pulse wave period detection circuit 16 as shown on the left side of the waveform of the fifth stage in FIG. It outputs to the circuit 17 (step S44). After that, the averaging circuit 17 calculates an average value (that is, average blood flow Qh) of the blood flow volume Q (k) of the sample number n (j) (step S17).
  • an average value that is, average blood flow Qh
  • step S43 when it is determined that the pulse amplitude A is less than the predetermined threshold (step S43: Yes), it is determined that the pulse amplitude A is relatively decreased. That is, it is determined that the pulse amplitude A does not have a magnitude that allows the reliable pulse wave cycle Tr to be calculated. Therefore, the sample value n (j) currently calculated by the pulse wave period detection circuit 16 is estimated to be low in reliability. For this reason, in this case, the selection circuit 44 outputs the average number of samples Nh calculated by the averaging circuit 43 to the averaging circuit 17 as shown on the right side of the fifth waveform in FIG. (Step S45). Thereafter, the averaging circuit 17 calculates an average value (that is, average blood flow Qh) of the blood flow Q (k) of the average number of samples Nh (step S17).
  • an average value that is, average blood flow Qh
  • the sampling circuit 19 synchronizes the average blood flow Qh calculated by the averaging circuit 17 with the outside of the blood flow detection device 1 (or the blood flow detection device 1 is synchronized with the timing generated by the timing generator 18).
  • the data is output to a processing block (not shown) provided inside (step S18).
  • Such a blood flow volume detection device 4 of the fourth embodiment can preferably enjoy the same effects as the various effects that the blood flow volume detection device 1 of the first embodiment can enjoy.
  • the blood flow rate detection device 4 of the fourth embodiment it is difficult to calculate the reliable sample number n due to the relative decrease in the pulse amplitude A (or reliable pulse wave period).
  • Tr an average value of the sample number n calculated in the past (that is, the average sample number Nh) can be used instead of the sample number n calculated in real time. Therefore, the blood flow rate detection device 4 of the fourth embodiment can calculate the average blood flow rate Qh that is correspondingly reliable even when the pulse amplitude A is relatively decreased.
  • the blood flow rate detection device 4 performs the pulse calculated in the past in addition to or instead of calculating the average value of the number of samples n calculated in the past (that is, the average number of samples Nh).
  • the average value of the wave period Tr may be calculated, and the number of samples n may be calculated from the average value of the pulse wave period Tr.
  • the blood flow rate detection device 4 according to the fourth embodiment uses a statistically calculated average in addition to or instead of using the average value of the number of samples n calculated in the past (that is, the average number of samples Nh).
  • the sample number n may be calculated from a reference value of a typical pulse wave period Tr.
  • the present invention can be appropriately changed without departing from the gist or the idea of the invention that can be read from the claims and the entire specification, and a blood flow detection device and method involving such a change are also included in the present invention. Included in technical thought.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Physiology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Hematology (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Description

流体評価装置及び方法Fluid evaluation apparatus and method

 本発明は、例えば血流等の流体の流量を検出することで当該流体を評価する流体評価装置及び方法の技術分野に関する。 The present invention relates to a technical field of a fluid evaluation apparatus and method for evaluating a fluid by detecting a flow rate of the fluid such as a blood flow.

 この種の流体評価装置として、例えば、レーザ血流計を用いて血流量を検出する装置が存在する(例えば、特許文献1参照)。特許文献1に開示されたレーザ血流計は、レーザ光を生体に照射し、その反射又は散乱の際に生ずるドップラシフトに起因した波長の変化に基づいて、生体の血流量を検出する。より具体的には、レーザ光が生体に照射されると、生体の血管内の血液の流れ(即ち、散乱体である赤血球の移動)に起因した散乱光が発生する。この散乱光の周波数は、元のレーザ光の周波数と比較して、散乱体である赤血球の移動速度に対応したレーザドップラ作用によって変化している。この散乱光と生体からのレーザ光の反射光との相互干渉により、いわゆる周波数差分信号に相当する信号成分が得られる。この信号成分が解析されることで、生体の血流量が検出される。 As this type of fluid evaluation device, for example, there is a device that detects a blood flow using a laser blood flow meter (see, for example, Patent Document 1). The laser blood flow meter disclosed in Patent Document 1 irradiates a living body with laser light and detects the blood flow rate of the living body based on a change in wavelength caused by Doppler shift that occurs during reflection or scattering. More specifically, when the living body is irradiated with laser light, scattered light is generated due to the flow of blood in the blood vessels of the living body (that is, movement of red blood cells that are scatterers). The frequency of the scattered light is changed by the laser Doppler action corresponding to the moving speed of the red blood cells that are the scatterers, as compared with the frequency of the original laser light. A signal component corresponding to a so-called frequency difference signal is obtained by mutual interference between the scattered light and the reflected light of the laser light from the living body. By analyzing this signal component, the blood flow of the living body is detected.

 尚、血流量を検出する装置を開示する文献ではないものの、本発明に関連するその他の先行技術文献として、特許文献2及び3がある。 In addition, although it is not the literature which discloses the apparatus which detects a blood flow rate, there exists patent document 2 and 3 as another prior art literature relevant to this invention.

特許第4685705号公報Japanese Patent No. 4658705 特許第3350521号公報Japanese Patent No. 3350521 特開2001-224593号公報JP 2001-224593 A

 血流量は、一般的には、心臓の拍動に合わせて大きく変動する(例えば、サインカーブ状に変化する)。このため、検出した血流量をそのまま出力する(例えば、レーザドップラ血流計の外部に接続されているディスプレイ等に出力する)と、時間軸に沿って刻一刻と変化する血流量が出力されることになる。場合によっては、このように出力される血流量(つまり、検出結果)が活用しにくいものとなってしまうおそれがある。従って、検出した血流量の時間平均値を算出すると共に、当該算出した時間平均値を出力することで、出力される血流量(つまり、検出結果)を活用しやすくする試みが考えられる。 The blood flow generally varies greatly according to the heart beat (for example, changes in a sine curve). For this reason, when the detected blood flow is output as it is (for example, output to a display connected to the outside of the laser Doppler blood flow meter), the blood flow that changes every moment along the time axis is output. It will be. In some cases, the blood flow volume (that is, the detection result) output in this way may be difficult to utilize. Therefore, it is conceivable to try to make it easier to utilize the output blood flow volume (that is, the detection result) by calculating the time average value of the detected blood flow volume and outputting the calculated time average value.

 ここで、血流量の時間平均値の算出方法によっては、算出された血流量の時間平均値が、本来算出されるべき(つまり、理想的な)血流量の時間平均値からずれてしまいかねないということが、本願発明者等の実験によって明らかになっている。例えば、所定の固定期間毎の血流量の時間平均値を算出するとする。一方で、血流の脈動周期は、心臓の拍動に応じて常に変動し得る。このため、所定の固定期間毎の血流量の時間平均値を算出すると、血流の脈動周期と血流量の時間平均値を算出する固定周期との間のずれに起因して、算出された血流量の時間平均値が、本来算出されるべき(つまり、理想的な)血流量の時間平均値よりも大きくなってしまう又は小さくなってしまいかねないということが、本願発明者等の実験によって明らかになっている。このような血流量の時間平均値の変動(つまり、誤差)は、血流量の高精度な検出という観点から好ましいとは言い難い。 Here, depending on the method of calculating the time average value of blood flow, the calculated time average value of blood flow may deviate from the time average value of blood flow that should be calculated (that is, ideal). This is clarified by experiments by the inventors of the present application. For example, suppose that the time average value of the blood flow volume for every predetermined fixed period is calculated. On the other hand, the pulsation cycle of blood flow can always vary according to the pulsation of the heart. For this reason, when the time average value of the blood flow rate for each predetermined fixed period is calculated, the calculated blood is caused by the difference between the pulsation cycle of the blood flow and the fixed cycle for calculating the time average value of the blood flow rate. It is clear from experiments by the present inventors that the time average value of the flow rate may become larger or smaller than the time average value of blood flow that should be originally calculated (that is, ideal). It has become. Such a change (that is, error) in the time average value of the blood flow rate is not preferable from the viewpoint of highly accurate detection of the blood flow rate.

 尚、上述した技術的問題は、血液の流量(血流量)の時間平均値が算出される場合に限らず、任意の流体の流量の時間平均値が算出される場合にも同様に発生し得る。 The above-described technical problem is not limited to the case where the time average value of the blood flow rate (blood flow rate) is calculated, but may similarly occur when the time average value of the flow rate of any fluid is calculated. .

 本発明は、例えば上記問題点に鑑みてなされたものであり、流体の流量の平均値を好適に算出することが可能な流体評価装置及び方法を提供することを課題とする。 The present invention has been made in view of the above-mentioned problems, for example, and an object of the present invention is to provide a fluid evaluation apparatus and method capable of suitably calculating an average value of fluid flow rates.

 上記課題を解決するための流体評価装置は、内部に流体が流れている被測定対象に対してレーザ光を照射する照射手段と、前記被測定対象に照射された前記レーザ光を検出することで得られる、前記レーザ光のドップラシフトに起因する周波数変化に基づいて、前記流体の流量を検出する第1検出手段と、前記流量の変化に基づいて、前記流量の脈動周期を検出する第2検出手段と、前記脈動周期に同期した期間毎の前記流量の平均値を算出する算出手段とを備える。 A fluid evaluation apparatus for solving the above-described problem is to detect an irradiation unit that irradiates a laser beam to a measurement target in which a fluid flows, and the laser beam irradiated to the measurement target. First detection means for detecting the flow rate of the fluid based on a change in frequency resulting from a Doppler shift of the laser light, and a second detection for detecting a pulsation cycle of the flow rate based on the change in the flow rate. And means for calculating an average value of the flow rate for each period synchronized with the pulsation cycle.

 上記課題を解決するための流体評価方法は、内部に流体が流れている被測定対象に対してレーザ光を照射する照射工程と、前記被測定対象に照射された前記レーザ光を検出することで得られる、前記レーザ光のドップラシフトに起因する周波数変化に基づいて、前記流体の流量を検出する第1検出工程と、前記流量の変化に基づいて、前記流体の脈動周期を検出する第2検出工程と、前記脈動周期に同期した期間毎の前記流量の平均値を算出する算出工程とを備える。 A fluid evaluation method for solving the above-described problem includes an irradiation step of irradiating a measurement target with a fluid flowing therein to a laser beam, and detecting the laser beam irradiated to the measurement target. A first detection step for detecting the flow rate of the fluid based on a frequency change resulting from a Doppler shift of the laser light, and a second detection for detecting a pulsation cycle of the fluid based on the change in the flow rate. And a calculation step of calculating an average value of the flow rate for each period synchronized with the pulsation cycle.

第1実施例の血流量検出装置の構成を示すブロック図である。It is a block diagram which shows the structure of the blood flow rate detection apparatus of 1st Example. 第1実施例の血流量検出装置の動作の流れを示すフローチャートである。It is a flowchart which shows the flow of operation | movement of the blood flow rate detection apparatus of 1st Example. 第1実施例の血流量検出装置の動作の際に観測される各種信号の波形を示すタイミングチャートである。It is a timing chart which shows the waveform of the various signals observed in the case of operation | movement of the blood flow rate detection apparatus of 1st Example. 脈波周期に対応するサンプル数を算出することなく、固定サンプル数個の血流量の平均値(平均血流量)を算出する比較例の血流量検出装置の動作の際に観測される各種信号の波形を示すタイミングチャートである。Without calculating the number of samples corresponding to the pulse wave period, various signals observed during the operation of the blood flow detection device of the comparative example that calculates the average value (average blood flow) of the blood flow of a fixed number of samples It is a timing chart which shows a waveform. 第2実施例の血流量検出装置の構成を示すブロック図である。It is a block diagram which shows the structure of the blood flow rate detection apparatus of 2nd Example. 第2実施例の血流量検出装置の動作の流れを示すフローチャートである。It is a flowchart which shows the flow of operation | movement of the blood flow rate detection apparatus of 2nd Example. 第3実施例の血流量検出装置の構成を示すブロック図である。It is a block diagram which shows the structure of the blood flow rate detection apparatus of 3rd Example. 第3実施例の血流量検出装置の動作の流れを示すフローチャートである。It is a flowchart which shows the flow of operation | movement of the blood flow rate detection apparatus of 3rd Example. 第4実施例の血流量検出装置の構成を示すブロック図である。It is a block diagram which shows the structure of the blood flow rate detection apparatus of 4th Example. 第4実施例の血流量検出装置の動作の流れを示すフローチャートである。It is a flowchart which shows the flow of operation | movement of the blood flow rate detection apparatus of 4th Example. 第4実施例の血流量検出装置の動作の際に観測される各種信号の波形を示すタイミングチャートである。It is a timing chart which shows the waveform of various signals observed in the case of operation of the blood flow rate detection device of the 4th example.

 以下、発明を実施するための形態として、流体評価装置及び方法の実施形態について順に説明する。 Hereinafter, embodiments of a fluid evaluation apparatus and method will be described in order as modes for carrying out the invention.

 (流体評価装置の実施形態)
 <1>
 本実施形態の流体評価装置は、内部に流体が流れている被測定対象に対してレーザ光を照射する照射手段と、前記被測定対象に照射された前記レーザ光を検出することで得られる、前記レーザ光のドップラシフトに起因する周波数変化に基づいて、前記流体の流量を検出する第1検出手段と、前記流量の変化に基づいて、前記流量の脈動周期を検出する第2検出手段と、前記脈動周期に同期した期間毎の前記流量の平均値を算出する算出手段とを備える。
(Embodiment of fluid evaluation apparatus)
<1>
The fluid evaluation apparatus according to the present embodiment is obtained by detecting an irradiation unit that irradiates laser light to a measurement target in which a fluid flows, and the laser light irradiated to the measurement target. First detection means for detecting the flow rate of the fluid based on a frequency change resulting from the Doppler shift of the laser light; and second detection means for detecting a pulsation cycle of the flow rate based on the change in the flow rate; Calculating means for calculating an average value of the flow rate for each period synchronized with the pulsation cycle.

 本実施形態の流体評価装置によれば、照射手段は、被測定対象に対してレーザ光を照射する。被測定対象の内部には、流体が流れている。レーザ光は、当該流体に対して照射されることが好ましい。尚、被測定対象及び流体の一例として、人間や動物等の生体及び血液があげられる。或いは、被測定対象及び流体の他の一例として、水道管及び水があげられる。 According to the fluid evaluation apparatus of this embodiment, the irradiation means irradiates the measurement target with laser light. A fluid flows inside the object to be measured. The laser light is preferably applied to the fluid. Examples of the measurement target and the fluid include living bodies such as humans and animals and blood. Or a water pipe and water are mention | raise | lifted as another example of a to-be-measured object and a fluid.

 第1検出手段は、レーザ光のドップラシフトに起因する周波数変化に基づいて、流体の流量を検出する。具体的には、例えば、レーザ光が被測定対象に照射されると、被測定対象の内部の流体の流れに起因した散乱光が発生する。この散乱光の周波数は、元のレーザ光の周波数と比較して、流体の移動速度に対応したレーザドップラ作用によって変化している。この散乱光と被測定対象からのレーザ光の反射光との相互干渉により、いわゆる周波数差分信号に相当する信号成分が得られる。この信号成分が解析されることで、流体の流量が検出される。 The first detection means detects the flow rate of the fluid based on the frequency change caused by the Doppler shift of the laser light. Specifically, for example, when laser light is irradiated on the measurement target, scattered light is generated due to the flow of fluid inside the measurement target. The frequency of the scattered light is changed by the laser Doppler action corresponding to the moving speed of the fluid as compared with the frequency of the original laser light. A signal component corresponding to a so-called frequency difference signal is obtained by the mutual interference between the scattered light and the reflected light of the laser light from the object to be measured. By analyzing this signal component, the flow rate of the fluid is detected.

 第2検出手段は、第1検出手段が検出した流量の変化に基づいて、流量の脈動周期を検出する。尚、「脈動周期」とは、例えば脈打つように変化する流量の変化周期(例えば、流量のピーク値から次のピーク値までの期間等)を意味する。尚、流体が血液である場合には、脈動周期は、実質的には、脈拍の周期(典型的には、1回の脈の周期)に一致する。 The second detection means detects the pulsation cycle of the flow rate based on the change in the flow rate detected by the first detection means. Note that the “pulsation cycle” means, for example, a flow rate change cycle that changes so as to pulsate (for example, a period from the peak value of the flow rate to the next peak value). When the fluid is blood, the pulsation cycle substantially coincides with the pulse cycle (typically, the cycle of one pulse).

 算出手段は、第1検出手段が検出した流量の平均値(例えば、所定時間毎の流量のサンプル値の和の平均値である時間平均値)を算出する。本実施形態では特に、算出手段は、脈動周期に同期した期間毎の流量の平均値(例えば、脈動周期に同期した期間毎の流量のサンプル値の和の平均値である時間平均値)を算出する。言い換えれば、算出手段は、脈動周期に同期した期間毎に、流量の平均値を算出する(例えば、脈動周期に同期した期間毎に、当該期間の間に検出された流量のサンプル値の和の平均値である時間平均値を算出する)。更に言い換えれば、算出手段は、脈動周期に同期して、流量の平均値を算出する(例えば、脈動周期に同期して、流量のサンプル値の和の平均値である時間平均値を算出する)。 The calculating means calculates an average value of the flow rate detected by the first detecting means (for example, a time average value that is an average value of the sum of the sample values of the flow rate every predetermined time). Particularly in the present embodiment, the calculation means calculates an average value of the flow rate for each period synchronized with the pulsation cycle (for example, a time average value that is an average value of the sum of the sample values of the flow rate for each period synchronized with the pulsation cycle). To do. In other words, the calculation means calculates the average value of the flow rate for each period synchronized with the pulsation cycle (for example, for each period synchronized with the pulsation cycle, the sum of the sample values of the flow rate detected during the period) Calculate the average time value). In other words, the calculation means calculates an average value of the flow rate in synchronization with the pulsation cycle (for example, calculates a time average value that is an average value of the sum of the sample values of the flow rate in synchronization with the pulsation cycle). .

 具体的な例をあげる。例えば、第2検出手段が検出した脈動周期が1秒であるとする。また、第1検出手段は、0.1秒毎に1回、流量を検出するとする。この場合には、算出手段は、脈動周期である1秒毎の流量の平均値を算出する。つまり、脈動周期である1秒の間に第1検出手段が10回流量を検出することを考慮すれば、算出手段は、脈動周期である1秒の間に検出される10個の流量のサンプル値の和の平均値を算出する。一方で、第2検出手段が検出した脈動周期が1.2秒に変化したとすると、算出手段は、脈動周期である1.2秒の間に検出される12個の流量のサンプル値の和の平均値を算出する。 I will give a specific example. For example, it is assumed that the pulsation cycle detected by the second detection means is 1 second. The first detection means detects the flow rate once every 0.1 seconds. In this case, the calculation means calculates an average value of the flow rate per second, which is a pulsation cycle. In other words, considering that the first detection means detects the flow rate 10 times during one second, which is the pulsation cycle, the calculation means can obtain 10 flow rate samples detected during the 1 second, which is the pulsation cycle. The average value of the sum of values is calculated. On the other hand, if the pulsation cycle detected by the second detection unit is changed to 1.2 seconds, the calculation unit calculates the sum of the sample values of 12 flow rates detected during the pulsation cycle of 1.2 seconds. The average value of is calculated.

 以上説明した本実施形態の流体評価装置によれば、何らかの要因によって変化し得る流体の流量の脈動周期に合わせて、流体の流量の平均値が算出される。つまり、流体の流量の平均値を算出する際の算出周期は、流体の流量の脈動周期に合わせて適宜調整される。従って、本実施形態の流体評価装置によれば、流体の流量の平均値を算出する際の算出周期が固定されている比較例の流体評価装置と比較して、流体の流量の脈動周期と流体の流量の平均値を算出する算出周期との間にずれが生ずる可能性が低減される。或いは、本実施形態の流体評価装置によれば、流体の流量の脈動周期と流体の流量の平均値を算出する算出周期との間にずれが生ずることは殆ど又は全くなくなる。従って、本実施形態の流体評価装置によれば、流体の流量の平均値を算出する際の算出周期が固定されている比較例の流体評価装置と比較して、算出手段が算出した流体の流量の平均値が、本来算出されるべき(つまり、理想的な)流体の流量の平均値からずれてしまう可能性が低減される。或いは、本実施形態の流体評価装置によれば、算出手段が算出した流体の流量の平均値が、本来算出されるべき(つまり、理想的な)流体の流量の平均値からずれてしまうことが殆ど又は全くなくなる。このため、本実施形態の流体評価装置は、流体の流量の平均値を好適に算出することができる。 According to the fluid evaluation apparatus of the present embodiment described above, the average value of the fluid flow rate is calculated in accordance with the pulsation cycle of the fluid flow rate that may change due to some factor. That is, the calculation cycle for calculating the average value of the fluid flow rate is appropriately adjusted according to the pulsation cycle of the fluid flow rate. Therefore, according to the fluid evaluation device of the present embodiment, the pulsation cycle of the fluid flow rate and the fluid are compared with the fluid evaluation device of the comparative example in which the calculation cycle when calculating the average value of the fluid flow rate is fixed. The possibility of a deviation from the calculation cycle for calculating the average value of the flow rates of is reduced. Alternatively, according to the fluid evaluation apparatus of the present embodiment, there is little or no deviation between the pulsation cycle of the fluid flow rate and the calculation cycle for calculating the average value of the fluid flow rate. Therefore, according to the fluid evaluation apparatus of the present embodiment, the flow rate of the fluid calculated by the calculation unit as compared with the fluid evaluation apparatus of the comparative example in which the calculation cycle when calculating the average value of the flow rate of the fluid is fixed. Is less likely to deviate from the average value of the flow rate of the fluid to be originally calculated (that is, ideal). Alternatively, according to the fluid evaluation apparatus of the present embodiment, the average value of the flow rate of the fluid calculated by the calculation unit may deviate from the average value of the flow rate of the fluid that should be originally calculated (that is, ideal). Little or no. For this reason, the fluid evaluation apparatus of this embodiment can calculate the average value of the flow volume of a fluid suitably.

 <2>
 本実施形態の流体評価装置の他の態様では、前記算出手段は、前記脈動周期と同一の期間毎の又は前記脈動周期のN(但し、Nは、2以上の整数)倍の期間毎の前記流量の平均値を算出する。
<2>
In another aspect of the fluid evaluation apparatus according to the present embodiment, the calculation unit is configured to perform the calculation every the same period as the pulsation period or every N times (where N is an integer of 2 or more) the pulsation period. Calculate the average flow rate.

 この態様によれば、算出手段は、脈動周期に同期した期間毎の流量の平均値(例えば、脈動周期に同期した期間毎の流量のサンプル値の和の平均値である時間平均値)を好適に算出することができる。 According to this aspect, the calculation means preferably uses an average value of the flow rate for each period synchronized with the pulsation cycle (for example, a time average value that is an average value of the sum of the sample values of the flow rate for each period synchronized with the pulsation cycle). Can be calculated.

 <3>
 本実施形態の流体評価装置の他の態様では、所望のタイミングで、前記算出手段が算出した前記平均値を出力する出力手段を更に備える。
<3>
In another aspect of the fluid evaluation apparatus of the present embodiment, the fluid evaluation apparatus further includes an output unit that outputs the average value calculated by the calculation unit at a desired timing.

 この態様によれば、算出手段が算出した流体の流量の平均値が、流体評価装置が備える又は流体評価装置の外部に取り付けられた他の装置(例えば、ディスプレイや解析回路等)に対して出力される。 According to this aspect, the average value of the flow rate of the fluid calculated by the calculation unit is output to another device (for example, a display or an analysis circuit) provided in the fluid evaluation device or attached to the outside of the fluid evaluation device. Is done.

 <4>
 上述の如く出力手段を備える流体評価装置の他の態様では、前記所望のタイミングは、周期的に現れる固定されたタイミングである。
<4>
In another aspect of the fluid evaluation apparatus including the output unit as described above, the desired timing is a fixed timing that appears periodically.

 この態様によれば、算出手段が算出した流体の流量の平均値が、一定周期で出力される。このため、算出手段が平均値を算出する期間(つまり、脈動周期に同期した期間)が変動したとしても、算出手段が算出した流体の流量の平均値が出力手段から出力される周期が一定に保たれる。従って、算出手段が算出した流体の流量の平均値がランダムな周期で出力される場合と比較して、出力手段から平均値を取得する他の装置にとっては、当該平均値を利用しやすくなる。 According to this aspect, the average value of the flow rate of the fluid calculated by the calculation means is output at a constant cycle. For this reason, even if the period during which the calculation unit calculates the average value (that is, the period synchronized with the pulsation cycle) fluctuates, the cycle in which the average value of the fluid flow rate calculated by the calculation unit is output from the output unit is constant. Kept. Therefore, compared with the case where the average value of the flow rate of the fluid calculated by the calculation unit is output in a random cycle, the average value can be easily used for other apparatuses that acquire the average value from the output unit.

 <5>
 上述の如く出力手段を備える流体評価装置の他の態様では、前記所望のタイミングは、ユーザによって指定されるタイミングである。
<5>
In another aspect of the fluid evaluation apparatus including the output unit as described above, the desired timing is a timing designated by the user.

 この態様によれば、算出手段が算出した流体の流量の平均値が、ユーザによって指定されるタイミングで出力される。このため、算出手段が平均値を算出する期間(つまり、脈動周期に同期した期間)が変動したとしても、算出手段が算出した流体の流量の平均値が出力手段から出力される周期は、ユーザの指示によって適宜調整される。従って、算出手段が算出した流体の流量の平均値がユーザの指示とは無関係な周期で出力される場合と比較して、出力手段から平均値を取得する他の装置にとっては、当該平均値を利用しやすくなる。 According to this aspect, the average value of the flow rate of the fluid calculated by the calculation means is output at a timing designated by the user. For this reason, even if the period during which the calculation unit calculates the average value (that is, the period synchronized with the pulsation cycle) fluctuates, the cycle in which the average value of the fluid flow rate calculated by the calculation unit is output from the output unit is It is adjusted appropriately according to the instructions. Therefore, compared with the case where the average value of the flow rate of the fluid calculated by the calculation means is output in a cycle unrelated to the user's instruction, for other devices that acquire the average value from the output means, the average value is It becomes easy to use.

 <6>
 本実施形態の流体評価装置の他の態様では、前記流量の振幅を検出する第3検出手段を更に備え、前記算出手段は、(i)前記振幅が所定量以上となる場合に、前記第1検出手段が検出した前記脈動周期に同期した期間毎の前記流量の平均値を算出し、(ii)前記振幅が前記所定量以上とならない場合に、前記第1検出手段が過去に検出した前記脈動周期に同期した期間毎の前記流量の平均値を算出する。
<6>
In another aspect of the fluid evaluation apparatus of this embodiment, the fluid evaluation device further includes third detection means for detecting the amplitude of the flow rate, and the calculation means (i) when the amplitude is equal to or greater than a predetermined amount, An average value of the flow rate for each period synchronized with the pulsation cycle detected by the detection unit is calculated, and (ii) the pulsation detected by the first detection unit in the past when the amplitude does not exceed the predetermined amount. The average value of the flow rate for each period synchronized with the cycle is calculated.

 この態様によれば、流量の振幅が相対的に小さい(つまり、所定量以上とならない)場合には、流量の振幅が相対的に大きい(つまり、所定量以上となる)場合と比較して、第2検出手段が検出する脈動周期の信頼性が相対的に低下していると推測される。従って、算出手段は、流量の振幅が相対的に小さい(つまり、所定量以上とならない)場合には、第2検出手段がリアルタイムで検出している信頼性が相対的に低い脈動周期に代えて、信頼性が相対的に高かった第2検出手段が過去に検出した脈動周期に同期した期間毎の流量の平均値を算出する。このため、第2検出手段がリアルタイムで検出している脈動周期に同期した期間毎の流量の平均値を常に算出する場合と比較して、算出手段が算出する平均値の信頼性を相対的に高めることができる。 According to this aspect, when the amplitude of the flow rate is relatively small (that is, not larger than the predetermined amount), compared to the case where the amplitude of the flow rate is relatively large (that is, larger than the predetermined amount), It is presumed that the reliability of the pulsation cycle detected by the second detection means is relatively lowered. Therefore, when the amplitude of the flow rate is relatively small (that is, not more than a predetermined amount), the calculating means replaces the pulsation cycle with relatively low reliability that the second detecting means detects in real time. The second detection means having relatively high reliability calculates the average value of the flow rate for each period synchronized with the pulsation cycle detected in the past. For this reason, the reliability of the average value calculated by the calculation means is relatively compared with the case where the average value of the flow rate for each period synchronized with the pulsation cycle detected in real time by the second detection means is always calculated. Can be increased.

 尚、算出手段は、第2検出手段が過去に検出した流量の脈動周期そのものに同期した期間(つまり、脈動周期と同一の期間又は脈動周期のN(但し、Nは、2以上の整数)倍若しくは1/N倍の期間)毎の流量の平均値を算出してもよい。或いは、算出手段は、第2検出手段が過去に検出した流量の脈動周期の時間平均値に同期した期間毎の流量の平均値を算出してもよい。 The calculation means is a period synchronized with the pulsation period itself of the flow rate detected by the second detection means in the past (that is, the same period as the pulsation period or N times the pulsation period (where N is an integer of 2 or more). Alternatively, the average value of the flow rate for each 1 / N times) may be calculated. Alternatively, the calculation unit may calculate an average value of the flow rate for each period synchronized with the time average value of the pulsation cycle of the flow rate detected in the past by the second detection unit.

 <7>
 前記流体は、前記被測定対象としての生体の内部を流れる血流であり、前記脈動周期は、前記血流の脈波周期である。
<7>
The fluid is a blood flow flowing inside the living body as the measurement target, and the pulsation cycle is a pulse wave cycle of the blood flow.

 この態様によれば、流体評価装置は、算出手段は、血液の脈波周期(つまり、脈拍の周期であり、典型的には、1回の脈の周期)に同期した期間毎の血液の流量(つまり、血流量)の平均値を算出することができる。 According to this aspect, in the fluid evaluation device, the calculation means has a blood flow rate for each period synchronized with the pulse wave cycle of blood (that is, the pulse cycle, typically one pulse cycle). An average value of (that is, blood flow volume) can be calculated.

 (流体評価方法の実施形態)
 <8>
 本実施形態の流体評価方法は、内部に流体が流れている被測定対象に対してレーザ光を照射する照射工程と、前記被測定対象に照射された前記レーザ光を検出することで得られる、前記レーザ光のドップラシフトに起因する周波数変化に基づいて、前記流体の流量を検出する第1検出工程と、前記流量の変化に基づいて、前記流体の脈動周期を検出する第2検出工程と、前記脈動周期に同期した期間毎の前記流量の平均値を算出する算出工程とを備える。
(Embodiment of fluid evaluation method)
<8>
The fluid evaluation method of the present embodiment is obtained by irradiating a laser beam to a measurement target in which a fluid flows inside, and detecting the laser beam irradiated to the measurement target. A first detection step for detecting a flow rate of the fluid based on a frequency change caused by a Doppler shift of the laser light; a second detection step for detecting a pulsation cycle of the fluid based on the change in the flow rate; A calculation step of calculating an average value of the flow rate for each period synchronized with the pulsation cycle.

 本実施形態の流体評価方法によれば、上述した本実施形態の流体評価装置が享受する各種効果を好適に享受することができる。 According to the fluid evaluation method of the present embodiment, various effects enjoyed by the fluid evaluation apparatus of the present embodiment described above can be suitably enjoyed.

 尚、本実施形態の流体評価装置が採用する各種態様に対応して、本実施形態の流体評価方法も、各種態様を採用してもよい。 Incidentally, in response to various aspects adopted by the fluid evaluation apparatus of the present embodiment, the fluid evaluation method of the present embodiment may adopt various aspects.

 本実施形態のこのような作用及び他の利得は次に説明する実施例から明らかにされる。 Such an operation and other advantages of the present embodiment will be clarified from examples described below.

 以上説明したように、本実施形態の流体評価装置は、照射手段と、第1検出手段と、第2検出手段と、算出手段とを備える。本実施形態の流体評価方法は、照射工程と、第1検出工程と、第2検出工程と、算出工程とを備える。従って、流体の流量の平均値を好適に算出することができる。 As described above, the fluid evaluation apparatus according to the present embodiment includes an irradiation unit, a first detection unit, a second detection unit, and a calculation unit. The fluid evaluation method of this embodiment includes an irradiation process, a first detection process, a second detection process, and a calculation process. Therefore, the average value of the flow rate of the fluid can be suitably calculated.

 以下、図面を参照しながら、流体評価装置の実施例について説明する。尚、以下では、流体評価装置を、生体の血管内を流れる血液の流量(つまり、血流量)を検出する血流量検出装置に適用した例について説明を進める。但し、流体評価装置は、血液以外の任意の流体(例えば、水道管の中を流れる水や、人工透析装置の血流回路を構成するチューブの中を流れる血液等)の流量を検出する任意の流量検出装置に適用されてもよい。 Hereinafter, embodiments of the fluid evaluation apparatus will be described with reference to the drawings. Hereinafter, an example in which the fluid evaluation device is applied to a blood flow rate detection device that detects a flow rate (that is, blood flow rate) of blood flowing in a blood vessel of a living body will be described. However, the fluid evaluation device can detect any flow rate of any fluid other than blood (for example, water flowing in a water pipe or blood flowing in a tube constituting a blood flow circuit of an artificial dialysis device). The present invention may be applied to a flow rate detection device.

 (1)第1実施例
 はじめに、図1から図4を参照しながら、第1実施例の血流量検出装置1について説明を進める。
(1) First Embodiment First , the blood flow detection device 1 of the first embodiment will be described with reference to FIGS.

 (1-1)血流量検出装置の構成
 はじめに、図1を参照しながら、第1実施例の血流量検出装置1の構成について説明する。図1は、第1実施例の血流量検出装置1の構成を示すブロック図である。
(1-1) Configuration of Blood Flow Detection Device First, the configuration of the blood flow detection device 1 of the first embodiment will be described with reference to FIG. FIG. 1 is a block diagram illustrating a configuration of a blood flow rate detection device 1 according to the first embodiment.

 図1に示すように、第1実施例の血流量検出装置1は、レーザ素子11と、受光素子12と、増幅器13と、A/D(Analogue to Digital)コンバータ14と、演算回路15と、脈波周期検出回路16と、平均化回路17と、タイミング生成器18と、サンプリング回路19とを備えている。 As shown in FIG. 1, the blood flow rate detection device 1 according to the first embodiment includes a laser element 11, a light receiving element 12, an amplifier 13, an analog / digital (A / D) converter 14, an arithmetic circuit 15, A pulse wave period detection circuit 16, an averaging circuit 17, a timing generator 18, and a sampling circuit 19 are provided.

 レーザ素子11は、「照射手段」の一具体例を構成しており、人間や動物等の生体100に対してレーザ光LBを照射する。このとき、レーザ素子11は、生体100内の血管に対してレーザ光LBを照射することが好ましい。 The laser element 11 constitutes a specific example of “irradiation means” and irradiates the living body 100 such as a human being or an animal with the laser beam LB. At this time, it is preferable that the laser element 11 irradiates the blood vessel in the living body 100 with the laser beam LB.

 受光素子12は、「第1検出手段」の一具体例を構成しており、生体100からのレーザ光LBの反射光と生体100からのレーザ光LBの散乱光との相互干渉によって生ずるビート信号光を受光する。受光素子12は、受光したビート信号光を電気信号に変換することで得られる検出電流を生成する。 The light receiving element 12 constitutes a specific example of “first detection means”, and beat signals generated by mutual interference between the reflected light of the laser light LB from the living body 100 and the scattered light of the laser light LB from the living body 100. Receives light. The light receiving element 12 generates a detection current obtained by converting the received beat signal light into an electrical signal.

 増幅器13は、「第1検出手段」の一具体例を構成しており、受光素子12から出力される検出電流を、電圧信号に変換した上で増幅する。 The amplifier 13 constitutes a specific example of “first detection means”, which amplifies the detection current output from the light receiving element 12 after converting it into a voltage signal.

 A/Dコンバータ14は、「第1検出手段」の一具体例を構成しており、増幅器13の出力(つまり、受光素子12が受光したビート信号光に応じた電圧信号)に対してA/D変換処理(つまり、量子化処理)を行う。その結果、A/Dコンバータ14は、受光素子12が受光したビート信号光に応じた電圧信号のサンプル値(つまり、量子化された電圧信号)を、演算回路15に出力する。 The A / D converter 14 constitutes a specific example of “first detection means”, and an A / D converter 14 outputs A / D with respect to the output of the amplifier 13 (that is, the voltage signal corresponding to the beat signal light received by the light receiving element 12). D conversion processing (that is, quantization processing) is performed. As a result, the A / D converter 14 outputs the sample value of the voltage signal (that is, the quantized voltage signal) corresponding to the beat signal light received by the light receiving element 12 to the arithmetic circuit 15.

 演算回路15は、「第1検出手段」の一具体例を構成しており、A/Dコンバータ14の出力(つまり、受光素子12が受光したビート信号光に応じた電圧信号のサンプル値)に対して、FFT(Fast Fourier Transform)を用いた周波数解析を行う。その結果、演算回路15は、血流量Q(k(尚、kは、血流量の演算のタイミングに応じて時系列的に割り当てられる変数))を算出する。 The arithmetic circuit 15 constitutes a specific example of “first detection means”, and outputs the output of the A / D converter 14 (that is, the sample value of the voltage signal corresponding to the beat signal light received by the light receiving element 12). On the other hand, frequency analysis using FFT (Fast Fourier Transform) is performed. As a result, the arithmetic circuit 15 calculates a blood flow rate Q (k (where k is a variable assigned in time series in accordance with the blood flow calculation timing)).

 脈波周期検出回路16は、「第2検出手段」の一具体例を構成しており、演算回路15が算出した血流量Q(k)の時間変化を解析することで、脈波周期(つまり、血液の脈動の周期であり、実質的には、脈拍の周期)Trを算出する。その後、脈波周期検出回路16は、脈波周期Trに対応する血流量Q(k)のサンプル数n(つまり、脈波周期と同一の期間に検出される血流量Q(k)のサンプル数n)を算出する。 The pulse wave cycle detection circuit 16 constitutes a specific example of “second detection means”, and by analyzing the temporal change of the blood flow rate Q (k) calculated by the calculation circuit 15, the pulse wave cycle (that is, , The period of blood pulsation, and substantially the period of pulse) Tr. Thereafter, the pulse wave cycle detection circuit 16 counts the number of samples n of the blood flow rate Q (k) corresponding to the pulse wave cycle Tr (that is, the number of samples of the blood flow rate Q (k) detected in the same period as the pulse wave cycle. n) is calculated.

 平均化回路17は、「算出手段」の一具体例を構成しており、n個の血流量Q(k)の平均値(以下、適宜“平均血流量Qh”と称する)を算出する。 The averaging circuit 17 constitutes a specific example of “calculation means”, and calculates an average value of n blood flow volumes Q (k) (hereinafter referred to as “average blood flow volume Qh” as appropriate).

 タイミング生成器18は、サンプリング回路19が平均血流量Qhを出力するタイミングを生成する。 The timing generator 18 generates a timing at which the sampling circuit 19 outputs the average blood flow rate Qh.

 サンプリング回路19は、「出力手段」の一具体例を構成しており、平均化回路17が算出した平均血流量Qhを、タイミング生成器18が生成したタイミングに同期して、血流量検出装置1の外部に(或いは、血流量検出装置1が内部に備える不図示の処理ブロックに)出力する。 The sampling circuit 19 constitutes a specific example of “output means”, and the blood flow detection device 1 synchronizes the average blood flow Qh calculated by the averaging circuit 17 with the timing generated by the timing generator 18. (Or to a processing block (not shown) included in the blood flow rate detection device 1).

 尚、血流量検出装置1は、平均化回路17が算出した平均血流量Qhをそのまま出力してもよい。この場合、血流量検出装置1は、タイミング生成器18及びサンプリング回路19を備えていなくともよい。 Note that the blood flow detection device 1 may output the average blood flow Qh calculated by the averaging circuit 17 as it is. In this case, the blood flow rate detection device 1 may not include the timing generator 18 and the sampling circuit 19.

 (1-2)血流量検出装置の動作
 続いて、図2及び図3を参照して、第1実施例の血流量検出装置1の動作の流れについて説明する。図2は、第1実施例の血流量検出装置1の動作の流れを示すフローチャートである。図3は、第1実施例の血流量検出装置1の動作の際に観測される各種信号の波形を示すタイミングチャートである。
(1-2) Operation of Blood Flow Detection Device Next, the flow of operation of the blood flow detection device 1 of the first embodiment will be described with reference to FIGS. FIG. 2 is a flowchart showing the operation flow of the blood flow rate detection device 1 of the first embodiment. FIG. 3 is a timing chart showing waveforms of various signals observed during the operation of the blood flow rate detection device 1 of the first embodiment.

 図2に示すように、レーザ素子11は、人間や動物等の生体100に対してレーザ光LBを照射する(ステップS11)。このとき、レーザ素子11は、生体100内の血管に対してレーザ光LBを照射することが好ましい。 As shown in FIG. 2, the laser element 11 irradiates a living body 100 such as a human being or an animal with the laser beam LB (step S11). At this time, it is preferable that the laser element 11 irradiates the blood vessel in the living body 100 with the laser beam LB.

 その後、受光素子12は、生体100からのレーザ光LBの反射光と生体100からのレーザ光LBの散乱光との相互干渉によって生ずるビート信号光を受光する(ステップS12)。具体的には、レーザ光LBが生体100に照射されると、生体100内の血管の内部の血液の流れ(即ち、散乱体である赤血球の移動)に起因した散乱光が発生する。この散乱光の周波数は、元のレーザ光の周波数と比較して、血液の移動速度に対応したレーザドップラ作用によって変化している。受光素子12は、このような散乱光と生体100からのレーザ光LBの反射光との相互干渉により生ずるビート信号光(いわゆる、周波数差分信号)を受光する。尚、ビート信号光を生じさせる散乱光として、生体100に照射されたレーザ光LBの透過光に相当する前方散乱光が用いられてもよいし、生体100に照射されたレーザ光LBの反射光に相当する後方散乱光が用いられてもよい。 Thereafter, the light receiving element 12 receives the beat signal light generated by the mutual interference between the reflected light of the laser light LB from the living body 100 and the scattered light of the laser light LB from the living body 100 (step S12). Specifically, when the living body 100 is irradiated with the laser light LB, scattered light is generated due to blood flow inside the blood vessel in the living body 100 (that is, movement of red blood cells that are scatterers). The frequency of the scattered light is changed by the laser Doppler action corresponding to the moving speed of the blood as compared with the frequency of the original laser light. The light receiving element 12 receives beat signal light (so-called frequency difference signal) generated by mutual interference between the scattered light and the reflected light of the laser light LB from the living body 100. Note that as the scattered light that generates the beat signal light, forward scattered light corresponding to the transmitted light of the laser light LB irradiated on the living body 100 may be used, or the reflected light of the laser light LB irradiated on the living body 100. Backscattered light corresponding to may be used.

 その後、受光素子12は、受光したビート信号光を電気信号に変換することで得られる検出電流を生成する。受光素子12は、生成した検出電流を、増幅器13に出力する。増幅器13は、受光素子12から出力される検出電流(つまり、受光素子12が受光したビート信号光に応じた検出電流)を、電圧信号に変換した上で増幅する。増幅器は、電圧信号をA/Dコンバータ14に出力する。 Thereafter, the light receiving element 12 generates a detection current obtained by converting the received beat signal light into an electric signal. The light receiving element 12 outputs the generated detection current to the amplifier 13. The amplifier 13 converts the detection current output from the light receiving element 12 (that is, the detection current corresponding to the beat signal light received by the light receiving element 12) into a voltage signal and then amplifies it. The amplifier outputs a voltage signal to the A / D converter 14.

 その後、A/Dコンバータ14は、増幅器13の出力(つまり、受光素子12が受光したビート信号光に応じた電圧信号)に対してA/D変換処理(つまり、量子化処理)を行う(ステップS13)。その結果、A/Dコンバータ14は、受光素子12が受光したビート信号光に応じた電圧信号のサンプル値(つまり、量子化された電圧信号)を、演算回路15に出力する。具体的には、例えば、A/Dコンバータ14は、A/Dコンバータ14のサンプリング周期をTaとすると、周期Ta毎に、受光素子12が受光したビート信号光に応じた電圧信号のサンプル値(つまり、量子化された電圧信号)を出力する。 Thereafter, the A / D converter 14 performs an A / D conversion process (that is, a quantization process) on the output of the amplifier 13 (that is, a voltage signal corresponding to the beat signal light received by the light receiving element 12) (steps). S13). As a result, the A / D converter 14 outputs the sample value of the voltage signal (that is, the quantized voltage signal) corresponding to the beat signal light received by the light receiving element 12 to the arithmetic circuit 15. Specifically, for example, if the sampling period of the A / D converter 14 is Ta, the A / D converter 14 samples the voltage signal sample value (corresponding to the beat signal light received by the light receiving element 12 for each period Ta ( That is, a quantized voltage signal) is output.

 演算回路15は、A/Dコンバータ14の出力(つまり、受光素子12が受光したビート信号光に応じた電圧信号のサンプル値)に対して、FFT(Fast Fourier Transform)を用いた周波数解析を行う。その結果、演算回路15は、血流量Q(k)を算出する(ステップS14:図3の1段目の波形参照)。具体的には、例えば、演算回路15は、ビート信号光に応じた電圧信号のサンプル値であって且つNf個のサンプル値に対してFFTを行う。演算回路15は、当該FFTを行うことで得られるパワースペクトルと周波数ベクトルとの乗算結果である1次モーメントを用いて、血流量Q(k)を算出する。FFTを用いた周波数解析による血流量Q(k)の算出方法については、公知の方法(例えば、特許第3313841号公報に開示された方法等)が用いられてもよいため、詳細な説明を省略する。演算回路15は、算出した血流量Q(k)を、脈波周期検出回路16及び平均化回路17の夫々に出力する。 The arithmetic circuit 15 performs frequency analysis using FFT (Fast Fourier Transform) on the output of the A / D converter 14 (that is, the sample value of the voltage signal corresponding to the beat signal light received by the light receiving element 12). . As a result, the arithmetic circuit 15 calculates a blood flow rate Q (k) (step S14: refer to the first waveform in FIG. 3). Specifically, for example, the arithmetic circuit 15 performs FFT on the sample value of the voltage signal corresponding to the beat signal light and the Nf sample values. The arithmetic circuit 15 calculates the blood flow rate Q (k) using a first moment that is a result of multiplication of the power spectrum and the frequency vector obtained by performing the FFT. Since a known method (for example, the method disclosed in Japanese Patent No. 3313841) may be used as a method for calculating the blood flow rate Q (k) by frequency analysis using FFT, detailed description thereof is omitted. To do. The arithmetic circuit 15 outputs the calculated blood flow volume Q (k) to each of the pulse wave period detection circuit 16 and the averaging circuit 17.

 尚、血流量Q(k)が算出される周期(つまり、1個の(言い換えれば、1サンプルの)血流量Q(k)が算出される周期)Tqは、図3の1段目及び2段目の夫々の波形に示すように、Nf×Ta(つまり、FFTのポイント数Nf×A/Dコンバータ14のサンプリング周期Ta)となる。 Note that the period in which the blood flow volume Q (k) is calculated (that is, the period in which one (in other words, one sample) blood flow volume Q (k) is calculated) Tq is the first and second stages in FIG. As shown in each waveform of the stage, Nf × Ta (that is, the number of FFT points Nf × the sampling period Ta of the A / D converter 14).

 その後、脈波周期検出回路16は、演算回路15が算出した血流量Q(k)の時間変化(つまり、複数サンプルの血流量Q(k)の時間変化)を解析することで、脈波周期Trを算出する(ステップS15)。例えば、脈波周期検出回路16は、図3の1段目及び3段目の夫々の波形に示すように、血流量Q(k)の振幅中心を基準として血流量Q(k)を2値化すると共に、当該2値化した血流量Q(k)の立ち上がりから立ち下がりまで(つまり、振幅中心以上のパルスと振幅中心未満のパルスとの組み合わせ)の期間を、脈波周期Trを算出してもよい。言い換えれば、例えば、脈波周期検出回路16は、図3の1段目及び3段目の夫々の波形に示すように、血流量Q(k)が当該血流量Q(k)の振幅中心と交わる点(いわゆる、ゼロクロス点)を検出することで、脈波周期Trを算出してもよい。或いは、例えば、脈波周期検出回路16は、血流量Q(k)のピーク値を検出すると共に、時間軸上で連続する2つのピーク値の間の間隔を、脈波周期Trとして算出してもよい。或いは、例えば、脈波周期検出回路16は、血流量Q(k)のボトム値を検出すると共に、時間軸上で連続する2つのボトム値の間の間隔を、脈波周期Trとして算出してもよい。 Thereafter, the pulse wave cycle detection circuit 16 analyzes the time change of the blood flow rate Q (k) calculated by the arithmetic circuit 15 (that is, the time change of the blood flow rate Q (k) of a plurality of samples) to thereby calculate the pulse wave cycle. Tr is calculated (step S15). For example, the pulse wave period detection circuit 16 binarizes the blood flow rate Q (k) with reference to the amplitude center of the blood flow rate Q (k) as shown in the first and third waveforms in FIG. And the pulse wave period Tr is calculated for a period from the rise to the fall of the binarized blood flow rate Q (k) (that is, a combination of a pulse with an amplitude greater than or equal to a pulse less than the amplitude center). May be. In other words, for example, the pulse wave period detection circuit 16 determines that the blood flow rate Q (k) is the amplitude center of the blood flow rate Q (k), as shown in the first and third waveforms in FIG. The pulse wave period Tr may be calculated by detecting a crossing point (so-called zero cross point). Alternatively, for example, the pulse wave period detection circuit 16 detects the peak value of the blood flow rate Q (k) and calculates the interval between two peak values that are continuous on the time axis as the pulse wave period Tr. Also good. Alternatively, for example, the pulse wave cycle detection circuit 16 detects the bottom value of the blood flow rate Q (k) and calculates the interval between two bottom values that are continuous on the time axis as the pulse wave cycle Tr. Also good.

 その後、脈波周期検出回路16は、図3の3段目の波形に示すように、脈波周期Trに対応する血流量Q(k)のサンプル数nを算出する(ステップS16)。具体的には、脈波周期検出回路16は、脈波周期Trと同一の期間に検出される血流量Q(k)のサンプル数nを算出する。言い換えれば、脈波周期検出回路16は、脈波周期Trと同一の期間に含まれる血流量Q(k)のサンプル数nを算出する。ここで、脈波周期検出回路16は、n=Tr/Tq=Tr/(Nf×Ta)という数式を用いて、サンプル数nを算出してもよい。言い換えれば、脈波周期検出回路16は、サンプル数n=脈波周期Tr/血流量Q(k)が算出される周期Tq=脈波周期Tr/(FFTのポイント数Nf×A/Dコンバータ14のサンプリング周期Ta)という数式を用いて、サンプル数nを算出してもよい。 Thereafter, the pulse wave period detection circuit 16 calculates the sample number n of the blood flow rate Q (k) corresponding to the pulse wave period Tr as shown in the third waveform in FIG. 3 (step S16). Specifically, the pulse wave cycle detection circuit 16 calculates the number of samples n of the blood flow rate Q (k) detected in the same period as the pulse wave cycle Tr. In other words, the pulse wave period detection circuit 16 calculates the number of samples n of the blood flow rate Q (k) included in the same period as the pulse wave period Tr. Here, the pulse wave period detection circuit 16 may calculate the number of samples n using a mathematical formula of n = Tr / Tq = Tr / (Nf × Ta). In other words, the pulse wave period detection circuit 16 is configured such that the number of samples n = pulse wave period Tr / period of blood flow Q (k) is calculated Tq = pulse wave period Tr / (number of FFT points Nf × A / D converter 14 The number of samples n may be calculated using a mathematical formula of sampling period Ta).

 具体的な数値を例にあげて説明する。A/Dコンバータ14のサンプリング周期Taが55マイクロ秒であり、FFTのポイント数Nfが1024個であり、脈波周期Trが0.9秒であるとする。この場合、サンプル数nは、n=0.9/(55×10-6×1024)≒16個となる。従って、脈波周期Tr=0.9秒に対応するサンプル数nは、16個となる。或いは、脈波周期Trが1.03秒に変わったとすると、サンプル数nは、n=1.03/(55×10-6×1024)≒18個となる。 A specific numerical value will be described as an example. It is assumed that the sampling period Ta of the A / D converter 14 is 55 microseconds, the FFT point number Nf is 1024, and the pulse wave period Tr is 0.9 seconds. In this case, the number of samples n is n = 0.9 / (55 × 10 −6 × 1024) ≈16. Therefore, the number of samples n corresponding to the pulse wave period Tr = 0.9 seconds is 16. Alternatively, if the pulse wave period Tr is changed to 1.03 seconds, the number of samples n is n = 1.03 / (55 × 10 −6 × 1024) ≈18.

 尚、脈波周期検出回路16が算出するサンプル数nは、整数であることが好ましい。従って、脈波周期検出回路16は、算出したサンプル数nが小数点以下の桁を有している場合には、小数点以下の桁が繰り上げられてもよい。或いは、小数点以下の桁が繰り下げられてもよい。或いは、小数点以下の桁が四捨五入されてもよい。 Note that the number of samples n calculated by the pulse wave period detection circuit 16 is preferably an integer. Therefore, the pulse wave period detection circuit 16 may carry up the digits after the decimal point when the calculated sample number n has digits after the decimal point. Alternatively, digits after the decimal point may be carried down. Alternatively, digits after the decimal point may be rounded off.

 図3の3段目の波形を見て分かるように、脈波周期検出回路16が算出するサンプル数nは、血流量Q(n)の時間変動に応じて適宜変わる。つまり、脈波周期検出回路16が算出するサンプル数nは、脈波周期Trに応じて適宜変わる。図3の3段目の波形には、ある期間においてはサンプル数n1が算出される一方で、その後の期間においてはサンプル数n2が算出される例が示されている。 As can be seen from the waveform at the third stage in FIG. 3, the number of samples n calculated by the pulse wave period detection circuit 16 changes as appropriate according to the time variation of the blood flow rate Q (n). That is, the number of samples n calculated by the pulse wave cycle detection circuit 16 changes appropriately according to the pulse wave cycle Tr. The third waveform in FIG. 3 shows an example in which the number of samples n1 is calculated in a certain period while the number of samples n2 is calculated in the subsequent period.

 尚、脈波周期検出回路16は、脈波周期TrのN(但し、Nは、2以上の整数)倍の期間に対応する血流量Q(k)のサンプル数nを算出してもよい。この場合、脈波周期検出回路16は、n=Tr×N/Tq=Tr×N/(Nf×Ta)という数式を用いて、サンプル数nを算出してもよい。 Note that the pulse wave period detection circuit 16 may calculate the number n of samples of the blood flow rate Q (k) corresponding to a period N (where N is an integer of 2 or more) times the pulse wave period Tr. In this case, the pulse wave period detection circuit 16 may calculate the number of samples n using a mathematical formula of n = Tr × N / Tq = Tr × N / (Nf × Ta).

 脈波周期検出回路16は、算出したサンプル数nを、平均化回路17に出力する。 The pulse wave period detection circuit 16 outputs the calculated sample number n to the averaging circuit 17.

 その後、平均化回路17は、ステップS16で算出されたサンプル数n個の血流量Q(k)の平均値(つまり、平均血流量Qh)を算出する(ステップS17)。例えば、平均化回路17は、(Q(0)+Q(1)+・・・+(Q(n-1))/nを、平均血流量Qhとして算出する。例えば、サンプル数n=16個である場合には、平均化回路17は、時間軸に沿って連続的に順次算出された16個の血流量Q(k)のサンプル値を加算すると共に当該加算結果をサンプル数n=16で除算することで、平均血流量Qhを算出する。 Thereafter, the averaging circuit 17 calculates an average value (that is, average blood flow Qh) of the blood flow volumes Q (k) of n samples calculated in step S16 (step S17). For example, the averaging circuit 17 calculates (Q (0) + Q (1) +... + (Q (n−1)) / n as the average blood flow rate Qh, for example, the number of samples n = 16 In this case, the averaging circuit 17 adds the sample values of the 16 blood flow rates Q (k) sequentially and sequentially calculated along the time axis, and outputs the addition result with the number of samples n = 16. By dividing, the average blood flow rate Qh is calculated.

 より具体的には、図3の波形を例にあげれば、平均化回路17は、ある期間(サンプル数n1が算出されている期間)において、(Q(0)+Q(1)+・・・+(Q(n1-1))/n1を、平均血流量Qhとして算出する。続いて、平均化回路17は、別の期間(サンプル数n2が算出されている期間)において、(Q(n1)+Q(n1+1)+・・・+(Q(n1+n2-1))/n2を、平均血流量Qhとして算出する。平均化回路17は、算出した平均血流量Qhを、サンプリング回路19に出力する。 More specifically, taking the waveform shown in FIG. 3 as an example, the averaging circuit 17 determines (Q (0) + Q (1) +... In a certain period (a period in which the number of samples n1 is calculated). Then, + (Q (n1-1)) / n1 is calculated as the average blood flow rate Qh, and then the averaging circuit 17 calculates (Q (n11) in another period (period in which the number of samples n2 is calculated). ) + Q (n1 + 1) +... + (Q (n1 + n2-1)) / n2 is calculated as the average blood flow Qh, and the averaging circuit 17 outputs the calculated average blood flow Qh to the sampling circuit 19. .

 その後、サンプリング回路19は、平均化回路17が算出した平均血流量Qhを、タイミング生成器18が生成したタイミング(図3の5段目の波形に示す再サンプリングタイミングに相当)に同期して、血流量検出装置1の外部に(或いは、血流量検出装置1が内部に備える不図示の処理ブロックに)出力する(ステップS18)。その結果、血流量検出装置1からは、図3の6段目(最終段)の波形に示す平均血流量Qhが出力される。 Thereafter, the sampling circuit 19 synchronizes the average blood flow Qh calculated by the averaging circuit 17 with the timing generated by the timing generator 18 (corresponding to the resampling timing shown in the fifth waveform in FIG. 3), The data is output to the outside of the blood flow detection device 1 (or to a processing block (not shown) included in the blood flow detection device 1) (step S18). As a result, the average blood flow Qh shown in the waveform of the sixth stage (final stage) in FIG.

 尚、タイミング生成器18が生成したタイミングは、固定タイミング(つまり、固定的な周期で現れるタイミング)であることが好ましい。この場合、血流量検出装置1から出力される平均血流量Qhの出力レート(言い換えれば、転送レート)が一定になる。 It should be noted that the timing generated by the timing generator 18 is preferably a fixed timing (that is, a timing that appears in a fixed cycle). In this case, the output rate (in other words, transfer rate) of the average blood flow Qh output from the blood flow detection device 1 is constant.

 また、サンプリング回路19は、平均血流量Qhを所定時間毎に更に平均化することで算出される時間平均値Qh’を出力してもよい。この場合、血流量検出装置1から出力する値(この場合には、時間平均値Qh’)の変動を少なくすることができる。 Further, the sampling circuit 19 may output a time average value Qh ′ calculated by further averaging the average blood flow rate Qh every predetermined time. In this case, fluctuations in the value output from the blood flow rate detection device 1 (in this case, the time average value Qh ′) can be reduced.

 ここで、比較例として、図4を参照しながら、脈波周期Trに対応するサンプル数nを算出することなく、固定サンプル数m個の血流量Q(k)の平均値(平均血流量Qh)を算出する血流量検出装置について説明する。図4は、脈波周期Trに対応するサンプル数nを算出することなく、固定サンプル数m個の血流量Q(k)の平均値(平均血流量Qh)を算出する比較例の血流量検出装置の動作の際に観測される各種信号の波形を示すタイミングチャートである。 Here, as a comparative example, with reference to FIG. 4, without calculating the sample number n corresponding to the pulse wave period Tr, the average value (average blood flow rate Qh) of the blood flow rate Q (k) of m fixed samples. ) Will be described. FIG. 4 shows a blood flow detection of a comparative example in which an average value (average blood flow Qh) of blood flow Q (k) of m fixed samples is calculated without calculating the number n of samples corresponding to the pulse wave period Tr. It is a timing chart which shows the waveform of various signals observed at the time of operation of an apparatus.

 図4に示すように、比較例の血流量検出装置は、(Q(0)+Q(1)+・・・+(Q(m-1))/mを、平均血流量Qhとして算出する。例えば、固定サンプル数m=21個である場合には、比較例の血流量検出装置は、時間軸に沿って連続的に順次算出された21個の血流量Q(k)のサンプル値を加算すると共に当該加算結果を固定サンプル数m=21で除算することで、平均血流量Qhを算出する。 As shown in FIG. 4, the blood flow rate detection device of the comparative example calculates (Q (0) + Q (1) +... + (Q (m−1)) / m as the average blood flow rate Qh. For example, when the fixed sample number m = 21, the blood flow rate detection device of the comparative example adds 21 sample values of the blood flow rate Q (k) sequentially and sequentially calculated along the time axis. At the same time, the average blood flow Qh is calculated by dividing the addition result by the fixed sample number m = 21.

 しかしながら、平均血流量Qhを算出する際のサンプル数mが固定されている一方で、脈波周期Trは変動する。というのも血流量は、心拍に伴って脈動しているがゆえに、当該心拍の変化に伴って脈波周期Trも当然に変動するからである。従って、脈波周期Trと平均血流量Qhを算出する算出周期(つまり、固定サンプル数mに対応する周期)との間にずれが生ずる。このため、平均血流量Qhを算出するタイミングによっては、図4の1段目の波形の相対的に左側に示すように、平均血流量Qhを算出するためのm個の血流量Q(k)のうち振幅中心よりも大きい側に分布する血流量Q(k)の個数が、振幅中心よりも小さい側に分布する血流量Q(k)の個数よりも多くなることがある。このようなm個の血流量Q(k)から算出される平均血流量Qhは、図4の4段目の波形の中央付近に示すように、本来算出されるべき理想的な平均血流量Qhよりも大きくなってしまう可能性が高い。或いは、平均血流量Qhを算出するタイミングによっては、図4の1段目の波形の相対的に右側に示すように、平均血流量Qhを算出するためのm個の血流量Q(k)のうち振幅中心よりも大きい側に分布する血流量Q(k)の個数が、振幅中心よりも小さい側に分布する血流量Q(k)の個数よりも少なくなることがある。このようなm個の血流量Q(k)から算出される平均血流量Qhは、図4の4段目の波形の相対的に右側に示すように、本来算出されるべき理想的な平均血流量Qhよりも小さくなってしまう可能性が高い。従って、平均血流量Qhを算出する際のサンプル数mが固定されている場合には、平均血流量Qhが誤差を含んでしまいかねない。 However, while the number m of samples for calculating the average blood flow rate Qh is fixed, the pulse wave cycle Tr varies. This is because the blood flow pulsates with the heartbeat, and the pulse wave period Tr naturally varies with the change of the heartbeat. Therefore, a deviation occurs between the pulse wave period Tr and the calculation period for calculating the average blood flow rate Qh (that is, the period corresponding to the fixed sample number m). Therefore, depending on the timing at which the average blood flow Qh is calculated, m blood flows Q (k) for calculating the average blood flow Qh, as shown on the relatively left side of the waveform in the first stage of FIG. Of these, the number of blood flows Q (k) distributed on the side larger than the amplitude center may be larger than the number of blood flows Q (k) distributed on the side smaller than the amplitude center. The average blood flow Qh calculated from the m blood flows Q (k) is ideally calculated as shown in the vicinity of the center of the fourth waveform in FIG. Is likely to become bigger than. Alternatively, depending on the timing of calculating the average blood flow Qh, as shown on the relatively right side of the waveform in the first stage of FIG. 4, m blood flows Q (k) for calculating the average blood flow Qh Of these, the number of blood flows Q (k) distributed on the side larger than the amplitude center may be smaller than the number of blood flows Q (k) distributed on the side smaller than the amplitude center. The average blood flow Qh calculated from such m blood flow volumes Q (k) is an ideal average blood to be originally calculated, as shown on the right side of the waveform in the fourth row in FIG. The possibility of becoming smaller than the flow rate Qh is high. Therefore, if the number of samples m when calculating the average blood flow Qh is fixed, the average blood flow Qh may include an error.

 加えて、脈波周期Trよりも低い周波数の平均血流量Qhの揺らぎとしては、生体100の自律神経系の活動が関与する揺らぎが存在することが本願発明者等の実験等によって確かめられている。この揺らぎは、生体100の反応として本来の揺らぎであり、検出されることが好ましい。ところで、比較例の血流量検出装置から出力される平均血流量Qhには、図4の4段目の波形に示すように、揺らぎ(つまり、平均血流量Qhの大小のゆらぎ)が生じている。しかしながら、この揺らぎの周期は、上述した平均血流量Qhの誤差(つまり、平均血流量Qhを算出する際に用いられる血流量Q(k)の振幅中心に対する偏り)に起因した誤った揺らぎであり、上述した生体100の反応として本来の揺らぎとは異なる。しかしながら、このような誤った揺らぎもまた、血流量Q(k)の脈波周期Trよりも低い周波数(つまり、長い周期)を有している。従って、比較例の血流量検出装置では、生体100の反応として本来の揺らぎと、上述した平均血流量Qhの誤差(つまり、平均血流量Qhを算出する際に用いられる血流量Q(k)の振幅中心に対する偏り)に起因した誤った揺らぎとが混在してしまう。このため、生体100の反応として本来の揺らぎが好適に検出されないおそれがある。 In addition, it has been confirmed by experiments of the inventors of the present application and the like that the fluctuation of the average blood flow Qh having a frequency lower than the pulse wave period Tr is a fluctuation involving the activity of the autonomic nervous system of the living body 100. . This fluctuation is an original fluctuation as a reaction of the living body 100 and is preferably detected. By the way, fluctuations (that is, fluctuations in the average blood flow Qh) are generated in the average blood flow Qh output from the blood flow detection device of the comparative example, as shown in the waveform at the fourth stage in FIG. . However, this fluctuation cycle is an erroneous fluctuation due to the above-described error in the average blood flow Qh (that is, the deviation of the blood flow Q (k) used for calculating the average blood flow Qh from the amplitude center). The reaction of the living body 100 is different from the original fluctuation. However, such an erroneous fluctuation also has a frequency (that is, a long period) lower than the pulse wave period Tr of the blood flow rate Q (k). Therefore, in the blood flow rate detection device of the comparative example, the original fluctuation as the reaction of the living body 100 and the error of the above average blood flow rate Qh (that is, the blood flow rate Q (k) used when calculating the average blood flow rate Qh). Incorrect fluctuations due to (bias with respect to the amplitude center) are mixed. For this reason, there is a possibility that the original fluctuation is not suitably detected as a reaction of the living body 100.

 しかるに、本実施例の血流量検出装置1によれば、脈波周期Trに合わせて算出されるサンプル数n個の血流量Q(k)の平均値(平均血流量Qh)が算出される。つまり、平均血流量Qhを算出する際の算出周期(つまり、サンプル数nに対応する周期)は、血流量Q(k)の脈動周期Trに合わせて適宜調整される。従って、本実施例の血流量検出装置1によれば、比較例の血流量検出装置と比較して、脈波周期Trと平均血流量Qhを算出する算出周期との間にずれが生ずる可能性が低減される。或いは、本実施例の血流量検出装置1によれば、脈波周期Trと平均血流量Qhを算出する算出周期との間にずれが生ずることは殆ど又は全くなくなる。本実施例の血流量検出装置1によれば、比較例の血流量検出装置と比較して、平均血流量Qhを算出するためのn個の血流量Q(k)のうち振幅中心よりも大きい側に分布する血流量Q(k)の個数が、振幅中心よりも小さい側に分布する血流量Q(k)の個数よりも多くなる可能性が低減される。或いは、本実施例の血流量検出装置1によれば、平均血流量Qhを算出するためのn個の血流量Q(k)のうち振幅中心よりも大きい側に分布する血流量Q(k)の個数が、振幅中心よりも小さい側に分布する血流量Q(k)の個数よりも多くなることが殆ど又は全くなくなる。同様に、本実施例の血流量検出装置1によれば、比較例の血流量検出装置と比較して、平均血流量Qhを算出するためのn個の血流量Q(k)のうち振幅中心よりも大きい側に分布する血流量Q(k)の個数が、振幅中心よりも小さい側に分布する血流量Q(k)の個数よりも少なくなる可能性が低減される。或いは、本実施例の血流量検出装置1によれば、平均血流量Qhを算出するためのn個の血流量Q(k)のうち振幅中心よりも大きい側に分布する血流量Q(k)の個数が、振幅中心よりも小さい側に分布する血流量Q(k)の個数よりも少なくなることが殆ど又は全くなくなる。つまり、本実施例の血流量検出装置1によれば、比較例の血流量検出装置と比較して、平均血流量Qhを算出するためのn個の血流量Q(k)のうち振幅中心よりも大きい側に分布する血流量Q(k)の個数と振幅中心よりも小さい側に分布する血流量Q(k)の個数とが概ね一致する可能性が高くなる。このため、本実施例の血流量検出装置1によれば、比較例の血流量検出装置と比較して、平均血流量Qhが、本来算出されるべき(つまり、理想的な)平均血流量Qhからずれてしまう可能性が低減される。或いは、本実施例の血流量検出装置1によれば、平均血流量Qhが、本来算出されるべき(つまり、理想的な)平均血流量Qhからずれてしまうことが殆ど又は全くなくなる。このため、本実施例の血流量検出装置1は、平均血流量Qhを好適に算出することができる。 However, according to the blood flow volume detection device 1 of the present embodiment, the average value (average blood flow volume Qh) of the blood flow volume Q (k) of n samples calculated in accordance with the pulse wave cycle Tr is calculated. That is, the calculation cycle (that is, the cycle corresponding to the number of samples n) when calculating the average blood flow Qh is appropriately adjusted according to the pulsation cycle Tr of the blood flow Q (k). Therefore, according to the blood flow rate detection device 1 of the present embodiment, there is a possibility that a deviation occurs between the pulse wave cycle Tr and the calculation cycle for calculating the average blood flow rate Qh, as compared with the blood flow rate detection device of the comparative example. Is reduced. Alternatively, according to the blood flow detection device 1 of the present embodiment, there is little or no deviation between the pulse wave cycle Tr and the calculation cycle for calculating the average blood flow Qh. According to the blood flow detection device 1 of the present embodiment, compared to the blood flow detection device of the comparative example, the n blood flows Q (k) for calculating the average blood flow Qh is larger than the amplitude center. The possibility that the number of blood flows Q (k) distributed on the side is larger than the number of blood flows Q (k) distributed on the side smaller than the amplitude center is reduced. Alternatively, according to the blood flow volume detection device 1 of the present embodiment, the blood flow volume Q (k) distributed on the side larger than the amplitude center among the n blood flow volumes Q (k) for calculating the average blood flow volume Qh. Of the blood flow volume Q (k) distributed on the side smaller than the amplitude center is hardly or not at all. Similarly, according to the blood flow volume detection device 1 of the present embodiment, the amplitude center of the n blood flow volumes Q (k) for calculating the average blood flow volume Qh as compared with the blood flow volume detection device of the comparative example. The possibility that the number of blood flow volumes Q (k) distributed on the larger side is smaller than the number of blood flow volumes Q (k) distributed on the smaller side than the amplitude center is reduced. Alternatively, according to the blood flow volume detection device 1 of the present embodiment, the blood flow volume Q (k) distributed on the side larger than the amplitude center among the n blood flow volumes Q (k) for calculating the average blood flow volume Qh. Is less than the number of blood flows Q (k) distributed on the side smaller than the amplitude center. That is, according to the blood flow volume detection device 1 of the present embodiment, compared to the blood flow volume detection device of the comparative example, from the center of amplitude among the n blood flow volumes Q (k) for calculating the average blood flow volume Qh. There is a high possibility that the number of blood flows Q (k) distributed on the larger side and the number of blood flows Q (k) distributed on the smaller side of the amplitude center substantially coincide with each other. Therefore, according to the blood flow detection device 1 of the present embodiment, the average blood flow Qh that should be originally calculated (that is, the ideal) average blood flow Qh is compared with the blood flow detection device of the comparative example. The possibility of shifting from is reduced. Alternatively, according to the blood flow rate detection device 1 of the present embodiment, the average blood flow rate Qh is hardly or completely deviated from the originally calculated (that is, ideal) average blood flow rate Qh. For this reason, the blood flow rate detection device 1 of the present embodiment can suitably calculate the average blood flow rate Qh.

 加えて、本実施例の血流量検出装置1によれば、比較例の血流量検出装置と比較して、上述した平均血流量Qhの誤差(つまり、平均血流量Qhを算出する際に用いられる血流量Q(k)の振幅中心に対する偏り)に起因した誤った揺らぎが平均血流量Qhに生ずる可能性が低減される。或いは、本実施例の血流量検出装置1によれば、上述した平均血流量Qhの誤差(つまり、平均血流量Qhを算出する際に用いられる血流量Q(k)の振幅中心に対する偏り)に起因した誤った揺らぎが平均血流量Qhに生ずることが殆ど又は全くなくなる。このため、本実施例の血流量検出装置1によれば、生体100の反応として本来の揺らぎが好適に検出される。 In addition, according to the blood flow detection device 1 of the present embodiment, compared to the blood flow detection device of the comparative example, it is used when calculating the error of the average blood flow Qh (that is, the average blood flow Qh). The possibility that an erroneous fluctuation due to the deviation of the blood flow volume Q (k) from the amplitude center) occurs in the average blood flow volume Qh is reduced. Alternatively, according to the blood flow detection device 1 of the present embodiment, the above-described error in the average blood flow Qh (that is, the deviation of the blood flow Q (k) used for calculating the average blood flow Qh from the amplitude center). The resulting false fluctuations will occur little or no in the mean blood flow Qh. For this reason, according to the blood flow rate detection device 1 of the present embodiment, the original fluctuation is suitably detected as the reaction of the living body 100.

 尚、上述した平均血流量Qhの誤差(つまり、平均血流量Qhを算出する際に用いられる血流量Q(k)の振幅中心に対する偏り)に起因した誤った揺らぎが発生することを防止するために、固定サンプル数mを相対的に大きくすることが考えられる。しかしながら、固定サンプル数mを相対的に大きくすると、生体100の反応として本来の揺らぎもあまた平均化されて抑圧されてしまう。このため、結果的に、生体100の反応として本来の揺らぎが好適に検出されないおそれがある。従って、サンプル数nを脈波周期Trに応じて算出する本実施例の血流量検出装置1は、実践上大変有利である。 In order to prevent an erroneous fluctuation due to the above-described error in the average blood flow Qh (that is, the deviation of the blood flow Q (k) with respect to the amplitude center used when calculating the average blood flow Qh). In addition, it is conceivable to relatively increase the fixed sample number m. However, if the number m of fixed samples is relatively large, the original fluctuation is also averaged and suppressed as a reaction of the living body 100. For this reason, as a result, the original fluctuation may not be suitably detected as the reaction of the living body 100. Therefore, the blood flow rate detection device 1 of the present embodiment that calculates the number of samples n according to the pulse wave period Tr is very advantageous in practice.

 加えて、本実施例の血流量検出装置1によれば、平均血流量Qhは、タイミング生成器18が生成したタイミングに同期して出力される。従って、血流量検出装置1から出力される平均血流量Qhの出力レート(言い換えれば、転送レート)を一定にすることができる。従って、血流量検出装置1と外部機器との間のインタフェース構築が比較的容易になる。 In addition, according to the blood flow detection device 1 of the present embodiment, the average blood flow Qh is output in synchronization with the timing generated by the timing generator 18. Therefore, the output rate (in other words, transfer rate) of the average blood flow Qh output from the blood flow detection device 1 can be made constant. Therefore, it is relatively easy to construct an interface between the blood flow rate detection device 1 and the external device.

 加えて、本実施例の血流量検出装置1は、血流量Q(k)が算出される周期Tqよりも長い周期で、平均化血流量Qhを出力することができる。このため、血流量Q(k)をそのまま出力する場合と比較して、単位時間当たりの出力データ量を低減することができる。更には、本実施例の血流量検出装置1は、血流量Q(k)が平均化された後に平均化した結果である平均血流量Qhが出力することができる。このため、血流量Q(k)をそのまま出力する場合と比較して、単位時間当たりの出力データ量を低減することができる。 In addition, the blood flow rate detection device 1 according to the present embodiment can output the averaged blood flow rate Qh in a cycle longer than the cycle Tq in which the blood flow rate Q (k) is calculated. For this reason, compared with the case where the blood flow rate Q (k) is output as it is, the output data amount per unit time can be reduced. Furthermore, the blood flow rate detection device 1 of the present embodiment can output an average blood flow rate Qh that is the result of averaging after the blood flow rate Q (k) is averaged. For this reason, compared with the case where the blood flow rate Q (k) is output as it is, the output data amount per unit time can be reduced.

 (2)第2実施例
 続いて、図5及び図6を参照しながら、第2実施例の血流量検出装置2について説明を進める。尚、以下の説明では、第1実施例の血流量検出装置1と同一の構成及び動作については、同一の参照符号及びステップ番号を付してその詳細な説明を省略する。
(2) Second Example Next, the blood flow detection device 2 of the second example will be described with reference to FIGS. 5 and 6. In the following description, the same configuration and operation as those of the blood flow rate detection device 1 of the first embodiment are denoted by the same reference numerals and step numbers, and detailed description thereof is omitted.

 (2-1)血流量検出装置の構成
 はじめに、図5を参照しながら、第2実施例の血流量検出装置2の構成について説明する。図5は、第2実施例の血流量検出装置2の構成を示すブロック図である。
(2-1) Configuration of Blood Flow Detection Device First, the configuration of the blood flow detection device 2 of the second embodiment will be described with reference to FIG. FIG. 5 is a block diagram showing the configuration of the blood flow rate detection device 2 of the second embodiment.

 図5に示すように、第2実施例の血流量検出装置2は、第1実施例の血流量検出装置1と比較して、タイミング生成器18を備えていない一方でユーザ要求受付部21を備えているという点で異なっている。第2実施例の血流量検出装置2が備えるその他の構成要素は、第1実施例の血流量検出装置1が備えるその他の構成要素と同一であってもよい。 As shown in FIG. 5, the blood flow rate detection device 2 of the second embodiment does not include the timing generator 18 and the user request reception unit 21 as compared with the blood flow rate detection device 1 of the first embodiment. It differs in that it is equipped. Other components included in the blood flow detection device 2 of the second embodiment may be the same as other components included in the blood flow detection device 1 of the first embodiment.

 ユーザ要求受付部21は、サンプリング回路19が平均血流量Qhを出力するタイミングを指示するユーザ(つまり、血流量検出装置2のユーザ)の指示を受け付ける。ユーザ要求受付部21は、受け付けたユーザの指示に同期して、サンプリング回路19が平均血流量Qhを出力するタイミングを生成する。 The user request accepting unit 21 accepts an instruction from a user (that is, a user of the blood flow detection device 2) instructing the timing at which the sampling circuit 19 outputs the average blood flow Qh. The user request receiving unit 21 generates a timing at which the sampling circuit 19 outputs the average blood flow Qh in synchronization with the received user instruction.

 (2-2)血流量検出装置の動作
 続いて、図6を参照して、第2実施例の血流量検出装置2の動作の流れについて説明する。図6は、第2実施例の血流量検出装置2の動作の流れを示すフローチャートである。
(2-2) Operation of Blood Flow Detection Device Next, the flow of operation of the blood flow detection device 2 of the second embodiment will be described with reference to FIG. FIG. 6 is a flowchart showing the operation flow of the blood flow rate detection device 2 of the second embodiment.

 図6に示すように、第2実施例の血流量検出装置2は、第1実施例の血流量検出装置1と同様に、ステップS11からステップS17までの動作を行う。 As shown in FIG. 6, the blood flow rate detection device 2 of the second embodiment performs the operations from step S11 to step S17 in the same manner as the blood flow rate detection device 1 of the first embodiment.

 その後、第2実施例の血流量検出装置2では、サンプリング回路19は、平均化回路17が算出した平均血流量Qhを、タイミング生成器18が生成したタイミング(つまり、ユーザの指示に応じたタイミング)に同期して、血流量検出装置1の外部に出力する(ステップS21)。 Thereafter, in the blood flow rate detection device 2 of the second embodiment, the sampling circuit 19 generates the average blood flow rate Qh calculated by the averaging circuit 17 at the timing when the timing generator 18 generates (that is, the timing according to the user's instruction). In synchronization with (), the blood is output to the outside of the blood flow detection device 1 (step S21).

 このような第2実施例の血流量検出装置2は、第1実施例の血流量検出装置1が享受することができる各種効果と同様の効果を、好適に享受することができる。 Such a blood flow volume detection device 2 of the second embodiment can preferably enjoy the same effects as the various effects that the blood flow volume detection device 1 of the first embodiment can enjoy.

 加えて、第2実施例の血流量検出装置2によれば、平均血流量Qhは、ユーザが指示したタイミングに同期して出力される。従って、血流量検出装置1から出力される平均血流量Qhの出力レート(言い換えれば、転送レート)を、ユーザの好みのレートに設定することができる。従って、ユーザフレンドリな態様で、平均血流量Qhが出力される。 In addition, according to the blood flow rate detection device 2 of the second embodiment, the average blood flow rate Qh is output in synchronization with the timing instructed by the user. Therefore, the output rate (in other words, transfer rate) of the average blood flow Qh output from the blood flow detection device 1 can be set to a user's favorite rate. Therefore, the average blood flow Qh is output in a user-friendly manner.

 (3)第3実施例
 続いて、図7及び図8を参照しながら、第3実施例の血流量検出装置3について説明を進める。尚、以下の説明では、第1実施例の血流量検出装置1と同一の構成及び動作については、同一の参照符号及びステップ番号を付してその詳細な説明を省略する。
(3) Third Example Next, the blood flow detection device 3 of the third example will be described with reference to FIGS. In the following description, the same configuration and operation as those of the blood flow rate detection device 1 of the first embodiment are denoted by the same reference numerals and step numbers, and detailed description thereof is omitted.

 (3-1)血流量検出装置の構成
 はじめに、図7を参照しながら、第3実施例の血流量検出装置3の構成について説明する。図7は、第3実施例の血流量検出装置3の構成を示すブロック図である。
(3-1) Configuration of Blood Flow Detection Device First, the configuration of the blood flow detection device 3 of the third embodiment will be described with reference to FIG. FIG. 7 is a block diagram illustrating a configuration of the blood flow rate detection device 3 according to the third embodiment.

 図7に示すように、第3実施例の血流量検出装置3は、第1実施例の血流量検出装置1と比較して、脈拍数算出回路31と、平均化回路32と、サンプリング回路33とを更に備えているという点で異なっている。第3実施例の血流量検出装置3が備えるその他の構成要素は、第1実施例の血流量検出装置1が備えるその他の構成要素と同一であってもよい。 As shown in FIG. 7, the blood flow rate detection device 3 according to the third embodiment has a pulse rate calculation circuit 31, an averaging circuit 32, and a sampling circuit 33 as compared with the blood flow rate detection device 1 according to the first embodiment. Is different in that it is further equipped. Other components included in the blood flow detection device 3 of the third embodiment may be the same as other components included in the blood flow detection device 1 of the first embodiment.

 脈拍数算出回路31は、脈波周期検出回路16が算出したサンプル数nから、脈拍数R(i(尚、iは、脈拍数の算出のタイミングに応じて時系列的に割り当てられる変数))を算出する。 The pulse rate calculation circuit 31 calculates the pulse rate R (i (where i is a variable assigned in time series according to the timing of calculating the pulse rate)) from the sample number n calculated by the pulse wave period detection circuit 16. Is calculated.

 平均化回路32は、所定数w個の脈拍数R(i)の平均値(以下、適宜“平均脈拍数Rh”と称する)を算出する。 The averaging circuit 32 calculates an average value of the predetermined number w of pulse rates R (i) (hereinafter, referred to as “average pulse rate Rh” as appropriate).

 サンプリング回路33は、平均化回路32が算出した平均脈拍数Rhを、タイミング生成器18が生成したタイミングに同期して、血流量検出装置1の外部に(或いは、血流量検出装置1が内部に備える不図示の処理ブロックに)出力する。 The sampling circuit 33 synchronizes the average pulse rate Rh calculated by the averaging circuit 32 to the outside of the blood flow detection device 1 (or the blood flow detection device 1 inside) in synchronization with the timing generated by the timing generator 18. Output to a processing block (not shown).

 尚、血流量検出装置1は、平均化回路32が算出した平均脈拍数Rhをそのまま出力してもよい。この場合、血流量検出装置1は、サンプリング回路33を備えていなくともよい。 In addition, the blood flow rate detection apparatus 1 may output the average pulse rate Rh calculated by the averaging circuit 32 as it is. In this case, the blood flow rate detection device 1 may not include the sampling circuit 33.

 (3-2)血流量検出装置の動作
 続いて、図8を参照して、第3実施例の血流量検出装置3の動作の流れについて説明する。図8は、第3実施例の血流量検出装置3の動作の流れを示すフローチャートである。
(3-2) Operation of Blood Flow Detection Device Next, the flow of operation of the blood flow detection device 3 of the third embodiment will be described with reference to FIG. FIG. 8 is a flowchart showing the flow of the operation of the blood flow rate detection device 3 of the third embodiment.

 図8に示すように、第3実施例の血流量検出装置3は、第1実施例の血流量検出装置1と同様に、ステップS11からステップS17までの動作を行う。 As shown in FIG. 8, the blood flow rate detection device 3 of the third embodiment performs the operations from step S11 to step S17 in the same manner as the blood flow rate detection device 1 of the first embodiment.

 第3実施例の血流量検出装置3では、更に、脈拍数算出回路31は、脈波周期検出回路16が算出したサンプル数nから、脈拍数R(i)を算出する(ステップS31)。ここで、脈拍数算出回路31は、R(i)=n×Tq×60=n×(Nf×Ta)×60という数式を用いて、脈拍数R(i)を算出してもよい。言い換えれば、脈拍数算出回路31は、脈拍数R(i)=サンプル数n×血流量Q(k)が算出される周期Tq×60=サンプル数n×(FFTのポイント数Nf×A/Dコンバータ14のサンプリング周期Ta)×60という数式を用いて、脈拍数R(i)を算出してもよい。 In the blood flow rate detection device 3 of the third embodiment, the pulse rate calculation circuit 31 further calculates the pulse rate R (i) from the sample number n calculated by the pulse wave cycle detection circuit 16 (step S31). Here, the pulse rate calculation circuit 31 may calculate the pulse rate R (i) using a mathematical formula R (i) = n × Tq × 60 = n × (Nf × Ta) × 60. In other words, the pulse rate calculation circuit 31 calculates the pulse rate R (i) = the number of samples n × the cycle in which the blood flow rate Q (k) is calculated Tq × 60 = the number of samples n × (the number of FFT points Nf × A / D). The pulse rate R (i) may be calculated using a mathematical formula of sampling period Ta) × 60 of the converter 14.

 具体的な数値を例にあげて説明する。A/Dコンバータ14のサンプリング周期Taが55マイクロ秒であり、FFTのポイント数Nfが1024個であり、サンプル数nが16個であるとする。この場合、脈拍数R(i)は、R(i)=16×(1024×55×10-6)×60≒54回/分となる。 A specific numerical value will be described as an example. It is assumed that the sampling period Ta of the A / D converter 14 is 55 microseconds, the FFT point number Nf is 1024, and the sample number n is 16. In this case, the pulse rate R (i) is R (i) = 16 × (1024 × 55 × 10 −6 ) × 60≈54 times / minute.

 脈拍数算出回路31は、算出した脈拍数R(i)を、平均化回路32に出力する。 The pulse rate calculation circuit 31 outputs the calculated pulse rate R (i) to the averaging circuit 32.

 その後、平均化回路32は、所定数w個の脈拍数R(i)の平均値(つまり、平均脈拍数Rh)を算出する(ステップS32)。例えば、平均化回路32は、(R(0)+R(1)+・・・+(R(w-1))/wを、平均脈拍数Rhとして算出する。例えば、所定数w=16である場合には、平均化回路32は、時間軸に沿って連続的に順次算出された6個の脈拍数R(i)のサンプル値を加算すると共に当該加算結果を所定数w=16で除算することで、平均脈拍数Rhを算出する。平均化回路32は、算出した平均脈拍数Rhを、サンプリング回路33に出力する。 Thereafter, the averaging circuit 32 calculates an average value (that is, average pulse rate Rh) of a predetermined number w of pulse rates R (i) (step S32). For example, the averaging circuit 32 calculates (R (0) + R (1) +... + (R (w−1)) / w as the average pulse rate Rh, for example, when the predetermined number w = 16. In some cases, the averaging circuit 32 adds the sample values of the six pulse rates R (i) sequentially and sequentially calculated along the time axis, and divides the addition result by the predetermined number w = 16. Then, the average pulse rate Rh is calculated, and the averaging circuit 32 outputs the calculated average pulse rate Rh to the sampling circuit 33.

 その後、サンプリング回路33は、平均化回路32が算出した平均脈拍数Rhを、タイミング生成器18が生成したタイミングに同期して、血流量検出装置1の外部に(或いは、血流量検出装置1が内部に備える不図示の処理ブロックに)出力する(ステップS33)。 After that, the sampling circuit 33 synchronizes the average pulse rate Rh calculated by the averaging circuit 32 to the outside of the blood flow detection device 1 (or the blood flow detection device 1 is synchronized with the timing generated by the timing generator 18). The data is output to a processing block (not shown) provided inside (step S33).

 このような第3実施例の血流量検出装置3は、第1実施例の血流量検出装置1が享受することができる各種効果と同様の効果を、好適に享受することができる。加えて、第2実施例の血流量検出装置2によれば、平均血流量Qhのみならず、平均脈拍数Rhもまた出力される。 Such a blood flow volume detection device 3 of the third embodiment can preferably enjoy the same effects as the various effects that the blood flow volume detection device 1 of the first embodiment can enjoy. In addition, according to the blood flow rate detection device 2 of the second embodiment, not only the average blood flow rate Qh but also the average pulse rate Rh is output.

 (4)第4実施例
 続いて、図9から図11を参照しながら、第4実施例の血流量検出装置4について説明を進める。尚、以下の説明では、第1実施例の血流量検出装置1と同一の構成及び動作については、同一の参照符号及びステップ番号を付してその詳細な説明を省略する。
(4) Fourth Example Next, the blood flow detection device 4 of the fourth example will be described with reference to FIGS. 9 to 11. In the following description, the same configuration and operation as those of the blood flow rate detection device 1 of the first embodiment are denoted by the same reference numerals and step numbers, and detailed description thereof is omitted.

 (4-1)血流量検出装置の構成
 はじめに、図9を参照しながら、第4実施例の血流量検出装置4の構成について説明する。図9は、第4実施例の血流量検出装置4の構成を示すブロック図である。
(4-1) Configuration of Blood Flow Detection Device First, the configuration of the blood flow detection device 4 of the fourth embodiment will be described with reference to FIG. FIG. 9 is a block diagram showing the configuration of the blood flow rate detection device 4 of the fourth embodiment.

 図9に示すように、第4実施例の血流量検出装置4は、第1実施例の血流量検出装置1と比較して、脈振幅検出回路41と、レベル低下検出回路42と、平均化回路43と、選択回路44とを更に備えているという点で異なっている。第4実施例の血流量検出装置4が備えるその他の構成要素は、第1実施例の血流量検出装置1が備えるその他の構成要素と同一であってもよい。 As shown in FIG. 9, the blood flow volume detection device 4 of the fourth embodiment is compared with the blood flow volume detection device 1 of the first embodiment, and the pulse amplitude detection circuit 41, the level decrease detection circuit 42, and the averaging The difference is that a circuit 43 and a selection circuit 44 are further provided. Other components included in the blood flow detection device 4 of the fourth embodiment may be the same as other components included in the blood flow detection device 1 of the first embodiment.

 脈振幅検出回路41は、「第3検出手段」の一具体例を構成しており、演算回路15が算出した血流量Q(k)の時間変化を解析することで、脈振幅(つまり、血流量Q(k)の振幅)Aを算出する。 The pulse amplitude detection circuit 41 constitutes a specific example of “third detection means”, and analyzes the time change of the blood flow volume Q (k) calculated by the arithmetic circuit 15 to thereby detect the pulse amplitude (that is, blood (Amplitude of flow rate Q (k)) A is calculated.

 レベル低下検出回路42は、脈振幅検出回路41が検出した脈振幅Aが所定閾値未満であるか否かを判定する。 The level drop detection circuit 42 determines whether or not the pulse amplitude A detected by the pulse amplitude detection circuit 41 is less than a predetermined threshold value.

 平均化回路43は、所定数p個のサンプル数n(j(尚、jは、サンプル数の算出のタイミングに応じて時系列的に割り当てられる変数))の平均値(以下、適宜“平均サンプル数Nh”と称する)を算出する。 The averaging circuit 43 calculates a predetermined number p of sample numbers n (j (where j is a variable assigned in time series according to the calculation timing of the number of samples)) (hereinafter referred to as “average samples” as appropriate). A number Nh ″).

 選択回路44は、レベル低下検出回路42での判定結果に基づいて、脈波周期検出回路16が検出するサンプル数n(j)及び平均化回路43が算出する平均サンプル値Nhのいずれかを、平均化回路17に出力する。 The selection circuit 44 selects one of the number of samples n (j) detected by the pulse wave period detection circuit 16 and the average sample value Nh calculated by the averaging circuit 43 based on the determination result in the level decrease detection circuit 42. Output to the averaging circuit 17.

 (4-2)血流量検出装置の動作
 続いて、図10及び図11を参照して、第4実施例の血流量検出装置4の動作の流れについて説明する。図10は、第4実施例の血流量検出装置4の動作の流れを示すフローチャートである。図11は、第4実施例の血流量検出装置4の動作の際に観測される各種信号の波形を示すタイミングチャートである。
(4-2) Operation of Blood Flow Rate Detection Device Next, the flow of operation of the blood flow rate detection device 4 of the fourth embodiment will be described with reference to FIGS. FIG. 10 is a flowchart showing the flow of the operation of the blood flow rate detection device 4 of the fourth embodiment. FIG. 11 is a timing chart showing waveforms of various signals observed during the operation of the blood flow rate detection device 4 of the fourth embodiment.

 図10に示すように、第4実施例の血流量検出装置4は、第1実施例の血流量検出装置1と同様に、ステップS11からステップS16までの動作を行う。 As shown in FIG. 10, the blood flow rate detection device 4 of the fourth embodiment performs the operations from step S11 to step S16 in the same manner as the blood flow rate detection device 1 of the first embodiment.

 その後、第4実施例の血流量検出装置4では、平均化回路43は、所定数p個のサンプル数n(j)の平均値(つまり、平均サンプル数Nh)を算出する。例えば、平均化回路43は、(n(0)+n(1)+・・・+(n(p-1))/pを、平均サンプル数Nhとして算出する。例えば、所定数p=16である場合には、平均化回路43は、時間軸に沿って連続的に順次算出された16個のサンプル数n(i)のサンプル値を加算すると共に当該加算結果を所定数p=16で除算することで、平均サンプル数Nhを算出する。平均化回路43は、算出した平均サンプル数Nhを、選択回路44に出力する。 Thereafter, in the blood flow rate detection device 4 of the fourth embodiment, the averaging circuit 43 calculates the average value (that is, the average sample number Nh) of the predetermined number p of sample numbers n (j). For example, the averaging circuit 43 calculates (n (0) + n (1) +... + (N (p−1)) / p as the average sample number Nh, for example, when the predetermined number p = 16. In some cases, the averaging circuit 43 adds the sample values of the 16 sample numbers n (i) sequentially and sequentially calculated along the time axis, and divides the addition result by the predetermined number p = 16. The averaging circuit 43 outputs the calculated average sample number Nh to the selection circuit 44.

 その後、脈波振幅検出回路41は、演算回路15が算出した血流量Q(k)の時間変化を解析することで、脈振幅(つまり、血流量Q(k)の振幅であり、図11の2段目の波形参照)Aを算出する(ステップS42)。尚、脈波振幅検出回路41は、血流量Q(k)の振幅を、そのまま脈振幅Aとして取り扱ってもよい。或いは、脈波振幅検出回路41は、血流量Q(k)の振幅に対して所定の演算処理を行うことで得られる値を、脈振幅Aとして取り扱ってもよい。 After that, the pulse wave amplitude detection circuit 41 analyzes the time change of the blood flow volume Q (k) calculated by the arithmetic circuit 15 to obtain the pulse amplitude (that is, the amplitude of the blood flow volume Q (k), which is shown in FIG. (Refer to the second-stage waveform) A is calculated (step S42). The pulse wave amplitude detection circuit 41 may handle the amplitude of the blood flow rate Q (k) as the pulse amplitude A as it is. Alternatively, the pulse wave amplitude detection circuit 41 may handle a value obtained by performing a predetermined calculation process on the amplitude of the blood flow volume Q (k) as the pulse amplitude A.

 その後、レベル低下検出回路42は、脈波振幅検出回路41が検出した脈振幅Aが所定閾値(図11の2段目の波形参照)未満であるか否かを判定する(ステップS43)。尚、所定閾値は、脈波周期Trを好適に算出することが可能な脈振幅Aの下限値が設定されていてもよい。或いは、所定閾値は、信頼性のある脈波周期Trを算出することが可能な脈振幅Aの下限値が設定されていてもよい。 Thereafter, the level decrease detection circuit 42 determines whether or not the pulse amplitude A detected by the pulse wave amplitude detection circuit 41 is less than a predetermined threshold (see the second waveform in FIG. 11) (step S43). The predetermined threshold value may be set to a lower limit value of the pulse amplitude A that can suitably calculate the pulse wave period Tr. Alternatively, the predetermined threshold may be set to a lower limit value of the pulse amplitude A that can calculate a reliable pulse wave period Tr.

 ステップS43の判定の結果、脈振幅Aが所定閾値未満でないと判定される場合には(ステップS43:No)、脈振幅Aが相対的に低下していないと判定される。つまり、脈振幅Aは、信頼性のある脈波周期Trを算出することが可能な大きさを有していると判定される。従って、現時点で脈波周期検出回路16が算出しているサンプル値n(j)は信頼性が高いと推測される。このため、この場合には、選択回路44は、図11の5段目の波形の相対的に左側に示すように、脈波周期検出回路16が検出するサンプル数n(j)を、平均化回路17に出力する(ステップS44)。その後、平均化回路17は、サンプル数n(j)個の血流量Q(k)の平均値(つまり、平均血流量Qh)を算出する(ステップS17)。 As a result of the determination in step S43, when it is determined that the pulse amplitude A is not less than the predetermined threshold value (step S43: No), it is determined that the pulse amplitude A is not relatively decreased. That is, it is determined that the pulse amplitude A has a magnitude that allows the reliable pulse wave period Tr to be calculated. Therefore, the sample value n (j) currently calculated by the pulse wave period detection circuit 16 is estimated to be highly reliable. Therefore, in this case, the selection circuit 44 averages the number of samples n (j) detected by the pulse wave period detection circuit 16 as shown on the left side of the waveform of the fifth stage in FIG. It outputs to the circuit 17 (step S44). After that, the averaging circuit 17 calculates an average value (that is, average blood flow Qh) of the blood flow volume Q (k) of the sample number n (j) (step S17).

 他方で、ステップS43の判定の結果、脈振幅Aが所定閾値未満であると判定される場合には(ステップS43:Yes)、脈振幅Aが相対的に低下していると判定される。つまり、脈振幅Aは、信頼性のある脈波周期Trを算出することが可能な大きさを有していないと判定される。従って、現時点で脈波周期検出回路16が算出しているサンプル値n(j)は信頼性が低いと推測される。このため、この場合には、選択回路44は、図11の5段目の波形の相対的に右側に示すように、平均化回路43が算出する平均サンプル数Nhを、平均化回路17に出力する(ステップS45)。その後、平均化回路17は、平均サンプル数Nh個の血流量Q(k)の平均値(つまり、平均血流量Qh)を算出する(ステップS17)。 On the other hand, as a result of the determination in step S43, when it is determined that the pulse amplitude A is less than the predetermined threshold (step S43: Yes), it is determined that the pulse amplitude A is relatively decreased. That is, it is determined that the pulse amplitude A does not have a magnitude that allows the reliable pulse wave cycle Tr to be calculated. Therefore, the sample value n (j) currently calculated by the pulse wave period detection circuit 16 is estimated to be low in reliability. For this reason, in this case, the selection circuit 44 outputs the average number of samples Nh calculated by the averaging circuit 43 to the averaging circuit 17 as shown on the right side of the fifth waveform in FIG. (Step S45). Thereafter, the averaging circuit 17 calculates an average value (that is, average blood flow Qh) of the blood flow Q (k) of the average number of samples Nh (step S17).

 その後、サンプリング回路19は、平均化回路17が算出した平均血流量Qhを、タイミング生成器18が生成したタイミングに同期して、血流量検出装置1の外部に(或いは、血流量検出装置1が内部に備える不図示の処理ブロックに)出力する(ステップS18)。 Thereafter, the sampling circuit 19 synchronizes the average blood flow Qh calculated by the averaging circuit 17 with the outside of the blood flow detection device 1 (or the blood flow detection device 1 is synchronized with the timing generated by the timing generator 18). The data is output to a processing block (not shown) provided inside (step S18).

 このような第4実施例の血流量検出装置4は、第1実施例の血流量検出装置1が享受することができる各種効果と同様の効果を、好適に享受することができる。 Such a blood flow volume detection device 4 of the fourth embodiment can preferably enjoy the same effects as the various effects that the blood flow volume detection device 1 of the first embodiment can enjoy.

 加えて、第4実施例の血流量検出装置4によれば、脈振幅Aの相対的な低下に起因して信頼性のあるサンプル数nの算出が困難(或いは、信頼性のある脈波周期Trの算出が困難)である場合には、リアルタイムで算出されているサンプル数nに代えて、過去に算出されたサンプル数nの平均値(つまり、平均サンプル数Nh)を用いることができる。従って、第4実施例の血流量検出装置4は、脈振幅Aが相対的に低下した場合であっても、相応に信頼性のある平均血流量Qhを算出することができる。 In addition, according to the blood flow rate detection device 4 of the fourth embodiment, it is difficult to calculate the reliable sample number n due to the relative decrease in the pulse amplitude A (or reliable pulse wave period). In the case where it is difficult to calculate Tr, an average value of the sample number n calculated in the past (that is, the average sample number Nh) can be used instead of the sample number n calculated in real time. Therefore, the blood flow rate detection device 4 of the fourth embodiment can calculate the average blood flow rate Qh that is correspondingly reliable even when the pulse amplitude A is relatively decreased.

 尚、サンプル数nの信頼性は、実質的には、脈波周期Trの信頼性と等価である。従って、第4実施例の血流量検出装置4は、過去に算出されたサンプル数nの平均値(つまり、平均サンプル数Nh)を算出することに加えて又は代えて、過去に算出された脈波周期Trの平均値を算出すると共に当該脈波周期Trの平均値からサンプル数nを算出してもよい。或いは、第4実施例の血流量検出装置4は、過去に算出されたサンプル数nの平均値(つまり、平均サンプル数Nh)を用いることに加えて又は代えて、統計的に算出される平均的な脈波周期Trの基準値からサンプル数nを算出してもよい。 Note that the reliability of the number of samples n is substantially equivalent to the reliability of the pulse wave period Tr. Therefore, the blood flow rate detection device 4 according to the fourth embodiment performs the pulse calculated in the past in addition to or instead of calculating the average value of the number of samples n calculated in the past (that is, the average number of samples Nh). The average value of the wave period Tr may be calculated, and the number of samples n may be calculated from the average value of the pulse wave period Tr. Alternatively, the blood flow rate detection device 4 according to the fourth embodiment uses a statistically calculated average in addition to or instead of using the average value of the number of samples n calculated in the past (that is, the average number of samples Nh). The sample number n may be calculated from a reference value of a typical pulse wave period Tr.

 尚、第1実施例から第4実施例で説明した各構成の一部を適宜組み合わせてもよい。この場合であっても、第1実施例から第4実施例で説明した各構成の一部を適宜組み合わせることで得られる血流量検出装置は、上述した各種効果を好適に享受することができる。 In addition, you may combine suitably a part of each structure demonstrated in 1st Example-4th Example. Even in this case, the blood flow rate detection device obtained by appropriately combining a part of the configurations described in the first to fourth embodiments can suitably enjoy the various effects described above.

 また、本発明は、請求の範囲及び明細書全体から読み取るこのできる発明の要旨又は思想に反しない範囲で適宜変更可能であり、そのような変更を伴う血流量検出装置及び方法もまた本発明の技術思想に含まれる。 In addition, the present invention can be appropriately changed without departing from the gist or the idea of the invention that can be read from the claims and the entire specification, and a blood flow detection device and method involving such a change are also included in the present invention. Included in technical thought.

 1、2、3、4 血流量検出装置
 11 レーザ素子
 12 受光素子
 13 増幅器
 14 A/Dコンバータ
 15 演算回路
 16 脈波周期検出回路
 17 平均化回路
 18 タイミング生成器
 19 サンプリング回路
 21 ユーザ要求受付部
 31 脈拍数算出回路
 32 平均化回路
 33 サンプリング回路
 41 脈波振幅検出回路
 42 レベル低下検出回路
 43 平均化回路
 44 選択回路
1, 2, 3, 4 Blood flow rate detection device 11 Laser element 12 Light receiving element 13 Amplifier 14 A / D converter 15 Arithmetic circuit 16 Pulse wave period detection circuit 17 Averaging circuit 18 Timing generator 19 Sampling circuit 21 User request receiving unit 31 Pulse rate calculation circuit 32 Averaging circuit 33 Sampling circuit 41 Pulse wave amplitude detection circuit 42 Level drop detection circuit 43 Averaging circuit 44 Selection circuit

Claims (8)

 内部に流体が流れている被測定対象に対してレーザ光を照射する照射手段と、
 前記被測定対象に照射された前記レーザ光を検出することで得られる、前記レーザ光のドップラシフトに起因する周波数変化に基づいて、前記流体の流量を検出する第1検出手段と、
 前記流量の変化に基づいて、前記流量の脈動周期を検出する第2検出手段と、
 前記脈動周期に同期した期間毎の前記流量の平均値を算出する算出手段と
 を備えることを特徴とする流体評価装置。
An irradiating means for irradiating a laser beam to a measurement target in which a fluid flows;
First detection means for detecting a flow rate of the fluid based on a frequency change caused by Doppler shift of the laser light, which is obtained by detecting the laser light irradiated on the measurement target;
Second detection means for detecting a pulsation cycle of the flow rate based on the change in the flow rate;
A fluid evaluation apparatus comprising: a calculating unit that calculates an average value of the flow rate for each period synchronized with the pulsation cycle.
 前記算出手段は、前記脈動周期と同一の期間毎の又は前記脈動周期のN(但し、Nは、2以上の整数)倍の期間毎の前記流量の平均値を算出することを特徴とする請求項1に記載の流体評価装置。 The said calculating means calculates the average value of the said flow volume for every period of the same period as the said pulsation period, or every period (N is an integer greater than or equal to 2) times the said pulsation period. Item 2. The fluid evaluation apparatus according to Item 1.  所望のタイミングで、前記算出手段が算出した前記平均値を出力する出力手段を更に備えることを特徴とする請求項1に記載の流体評価装置。 2. The fluid evaluation apparatus according to claim 1, further comprising output means for outputting the average value calculated by the calculation means at a desired timing.  前記所望のタイミングは、周期的に現れる固定されたタイミングであることを特徴とする請求項3に記載の流体評価装置。 4. The fluid evaluation apparatus according to claim 3, wherein the desired timing is a fixed timing that appears periodically.  前記所望のタイミングは、ユーザによって指定されるタイミングであることを特徴とする請求項3に記載の流体評価装置。 4. The fluid evaluation apparatus according to claim 3, wherein the desired timing is a timing designated by a user.  前記流量の振幅を検出する第3検出手段を更に備え、
 前記算出手段は、(i)前記振幅が所定量以上となる場合に、前記第1検出手段が検出した前記脈動周期に同期した期間毎の前記流量の平均値を算出し、(ii)前記振幅が前記所定量以上とならない場合に、前記第1検出手段が過去に検出した前記脈動周期に同期した期間毎の前記流量の平均値を算出することを特徴とする請求項1に記載の流体評価装置。
Further comprising third detecting means for detecting the amplitude of the flow rate;
The calculating means calculates (i) an average value of the flow rate for each period synchronized with the pulsation cycle detected by the first detecting means when the amplitude is a predetermined amount or more, and (ii) the amplitude 2. The fluid evaluation according to claim 1, wherein when the value does not exceed the predetermined amount, the first detection unit calculates an average value of the flow rate for each period synchronized with the pulsation cycle detected in the past. apparatus.
 前記流体は、前記被測定対象としての生体の内部を流れる血流であり、
 前記脈動周期は、前記血流の脈波周期であることを特徴とする請求項1に記載の流体評価装置。
The fluid is a blood flow flowing inside a living body as the measurement target,
The fluid evaluation apparatus according to claim 1, wherein the pulsation cycle is a pulse wave cycle of the blood flow.
 内部に流体が流れている被測定対象に対してレーザ光を照射する照射工程と、
 前記被測定対象に照射された前記レーザ光を検出することで得られる、前記レーザ光のドップラシフトに起因する周波数変化に基づいて、前記流体の流量を検出する第1検出工程と、
 前記流量の変化に基づいて、前記流体の脈動周期を検出する第2検出工程と、
 前記脈動周期に同期した期間毎の前記流量の平均値を算出する算出工程と
 を備えることを特徴とする流体評価方法。
An irradiation step of irradiating a measurement target in which a fluid flows inside with a laser beam;
A first detection step for detecting a flow rate of the fluid based on a frequency change caused by Doppler shift of the laser light, which is obtained by detecting the laser light irradiated on the measurement target;
A second detection step of detecting a pulsation cycle of the fluid based on the change in the flow rate;
A fluid evaluation method comprising: a calculation step of calculating an average value of the flow rate for each period synchronized with the pulsation cycle.
PCT/JP2012/057379 2012-03-22 2012-03-22 Fluid evaluation device and method Ceased WO2013140583A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PCT/JP2012/057379 WO2013140583A1 (en) 2012-03-22 2012-03-22 Fluid evaluation device and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2012/057379 WO2013140583A1 (en) 2012-03-22 2012-03-22 Fluid evaluation device and method

Publications (1)

Publication Number Publication Date
WO2013140583A1 true WO2013140583A1 (en) 2013-09-26

Family

ID=49222077

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/057379 Ceased WO2013140583A1 (en) 2012-03-22 2012-03-22 Fluid evaluation device and method

Country Status (1)

Country Link
WO (1) WO2013140583A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017113320A (en) * 2015-12-24 2017-06-29 パイオニア株式会社 Fluid evaluation apparatus and method, computer program, and recording medium
EP3534122A4 (en) * 2016-10-25 2020-06-17 Pioneer Corporation FLUID MEASURING DEVICE

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07124162A (en) * 1993-11-05 1995-05-16 Toshiba Corp Cardiac function test system
JP2003019119A (en) * 2001-07-10 2003-01-21 Canon Inc Fundus blood flow meter
JP2008506472A (en) * 2004-07-14 2008-03-06 エドワーズ ライフサイエンシーズ コーポレイション Real-time measurement of ventricular stroke volume variance by continuous arterial pulse contour analysis
JP2008080106A (en) * 2006-08-28 2008-04-10 Toshiba Corp Ultrasonic diagnostic apparatus and data processing program for ultrasonic diagnostic apparatus
JP2010523286A (en) * 2007-04-10 2010-07-15 ユニヴァーシティー オブ サザン カリフォルニア Method and system for blood flow measurement using Doppler optical coherence tomography

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07124162A (en) * 1993-11-05 1995-05-16 Toshiba Corp Cardiac function test system
JP2003019119A (en) * 2001-07-10 2003-01-21 Canon Inc Fundus blood flow meter
JP2008506472A (en) * 2004-07-14 2008-03-06 エドワーズ ライフサイエンシーズ コーポレイション Real-time measurement of ventricular stroke volume variance by continuous arterial pulse contour analysis
JP2008080106A (en) * 2006-08-28 2008-04-10 Toshiba Corp Ultrasonic diagnostic apparatus and data processing program for ultrasonic diagnostic apparatus
JP2010523286A (en) * 2007-04-10 2010-07-15 ユニヴァーシティー オブ サザン カリフォルニア Method and system for blood flow measurement using Doppler optical coherence tomography

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017113320A (en) * 2015-12-24 2017-06-29 パイオニア株式会社 Fluid evaluation apparatus and method, computer program, and recording medium
EP3534122A4 (en) * 2016-10-25 2020-06-17 Pioneer Corporation FLUID MEASURING DEVICE
US11486892B2 (en) 2016-10-25 2022-11-01 Air Water Biodesign Inc. Fluid measuring device

Similar Documents

Publication Publication Date Title
JP5806390B2 (en) Fluid evaluation apparatus and method
EP2730302B1 (en) Hypotension-predicting device
US9936888B2 (en) Pulse period calculation device and biosensor equipped with the same
JP6090424B2 (en) Pulse wave propagation time measurement device
US20160113530A1 (en) Method and apparatus for acquiring biological information, and wrist watch-type terminal using the same
JP2003153881A (en) Blood component measuring instrument
KR101571500B1 (en) Apparatus for measuring bio signal using bandwidth of pulse signal, and system for monitoring user comprising the same
RU2012101554A (en) ECG CONTROL WITH REDUCED FALSE ALARMS IN CONNECTION WITH ASYSTOLIA
JP5562805B2 (en) Pulse rate measuring method and blood oxygen saturation measuring method
JP2011115459A (en) Device and method for detecting biological information
KR20090096803A (en) Method and apparatus for removing noise from pulse signal and recording medium
JP2018007894A (en) Measuring apparatus, measuring method and measuring program
JP2003265422A (en) Sphygmograph device
CN101940476A (en) Blood oxygen saturation detection method and system
US10568581B2 (en) Pulsimeter, frequency analysis device, and pulse measurement method
WO2013140583A1 (en) Fluid evaluation device and method
JP6194031B2 (en) Fluid evaluation apparatus and method
JP2008132012A (en) Pulse wave detector
JP2003153882A (en) Blood component measuring instrument
JPWO2013140583A1 (en) Fluid evaluation apparatus and method
JP2017192629A (en) Measurement device, measurement method, and electronic apparatus
JP2005160944A (en) Apparatus and method for measuring blood flow rate
JP4322302B1 (en) Biological signal analyzer
JPWO2015129025A1 (en) Measuring device, pulse oximeter, measuring method, computer program, and recording medium
JP2017049266A (en) Fluid evaluation apparatus and method

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12871955

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2014505910

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 12871955

Country of ref document: EP

Kind code of ref document: A1