WO2013140583A1 - Fluid evaluation device and method - Google Patents
Fluid evaluation device and method Download PDFInfo
- 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
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
- A61B5/024—Measuring pulse rate or heart rate
- A61B5/02416—Measuring pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
- A61B5/026—Measuring blood flow
- A61B5/0261—Measuring 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
本発明は、例えば血流等の流体の流量を検出することで当該流体を評価する流体評価装置及び方法の技術分野に関する。 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
尚、血流量を検出する装置を開示する文献ではないものの、本発明に関連するその他の先行技術文献として、特許文献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.
血流量は、一般的には、心臓の拍動に合わせて大きく変動する(例えば、サインカーブ状に変化する)。このため、検出した血流量をそのまま出力する(例えば、レーザドップラ血流計の外部に接続されているディスプレイ等に出力する)と、時間軸に沿って刻一刻と変化する血流量が出力されることになる。場合によっては、このように出力される血流量(つまり、検出結果)が活用しにくいものとなってしまうおそれがある。従って、検出した血流量の時間平均値を算出すると共に、当該算出した時間平均値を出力することで、出力される血流量(つまり、検出結果)を活用しやすくする試みが考えられる。 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.
以下、発明を実施するための形態として、流体評価装置及び方法の実施形態について順に説明する。 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
(1-1)血流量検出装置の構成
はじめに、図1を参照しながら、第1実施例の血流量検出装置1の構成について説明する。図1は、第1実施例の血流量検出装置1の構成を示すブロック図である。
(1-1) Configuration of Blood Flow Detection Device First, the configuration of the blood
図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
レーザ素子11は、「照射手段」の一具体例を構成しており、人間や動物等の生体100に対してレーザ光LBを照射する。このとき、レーザ素子11は、生体100内の血管に対してレーザ光LBを照射することが好ましい。
The
受光素子12は、「第1検出手段」の一具体例を構成しており、生体100からのレーザ光LBの反射光と生体100からのレーザ光LBの散乱光との相互干渉によって生ずるビート信号光を受光する。受光素子12は、受光したビート信号光を電気信号に変換することで得られる検出電流を生成する。
The
増幅器13は、「第1検出手段」の一具体例を構成しており、受光素子12から出力される検出電流を、電圧信号に変換した上で増幅する。
The
A/Dコンバータ14は、「第1検出手段」の一具体例を構成しており、増幅器13の出力(つまり、受光素子12が受光したビート信号光に応じた電圧信号)に対してA/D変換処理(つまり、量子化処理)を行う。その結果、A/Dコンバータ14は、受光素子12が受光したビート信号光に応じた電圧信号のサンプル値(つまり、量子化された電圧信号)を、演算回路15に出力する。
The A /
演算回路15は、「第1検出手段」の一具体例を構成しており、A/Dコンバータ14の出力(つまり、受光素子12が受光したビート信号光に応じた電圧信号のサンプル値)に対して、FFT(Fast Fourier Transform)を用いた周波数解析を行う。その結果、演算回路15は、血流量Q(k(尚、kは、血流量の演算のタイミングに応じて時系列的に割り当てられる変数))を算出する。
The
脈波周期検出回路16は、「第2検出手段」の一具体例を構成しており、演算回路15が算出した血流量Q(k)の時間変化を解析することで、脈波周期(つまり、血液の脈動の周期であり、実質的には、脈拍の周期)Trを算出する。その後、脈波周期検出回路16は、脈波周期Trに対応する血流量Q(k)のサンプル数n(つまり、脈波周期と同一の期間に検出される血流量Q(k)のサンプル数n)を算出する。
The pulse wave
平均化回路17は、「算出手段」の一具体例を構成しており、n個の血流量Q(k)の平均値(以下、適宜“平均血流量Qh”と称する)を算出する。
The averaging
タイミング生成器18は、サンプリング回路19が平均血流量Qhを出力するタイミングを生成する。
The
サンプリング回路19は、「出力手段」の一具体例を構成しており、平均化回路17が算出した平均血流量Qhを、タイミング生成器18が生成したタイミングに同期して、血流量検出装置1の外部に(或いは、血流量検出装置1が内部に備える不図示の処理ブロックに)出力する。
The
尚、血流量検出装置1は、平均化回路17が算出した平均血流量Qhをそのまま出力してもよい。この場合、血流量検出装置1は、タイミング生成器18及びサンプリング回路19を備えていなくともよい。
Note that the blood
(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
図2に示すように、レーザ素子11は、人間や動物等の生体100に対してレーザ光LBを照射する(ステップS11)。このとき、レーザ素子11は、生体100内の血管に対してレーザ光LBを照射することが好ましい。
As shown in FIG. 2, the
その後、受光素子12は、生体100からのレーザ光LBの反射光と生体100からのレーザ光LBの散乱光との相互干渉によって生ずるビート信号光を受光する(ステップS12)。具体的には、レーザ光LBが生体100に照射されると、生体100内の血管の内部の血液の流れ(即ち、散乱体である赤血球の移動)に起因した散乱光が発生する。この散乱光の周波数は、元のレーザ光の周波数と比較して、血液の移動速度に対応したレーザドップラ作用によって変化している。受光素子12は、このような散乱光と生体100からのレーザ光LBの反射光との相互干渉により生ずるビート信号光(いわゆる、周波数差分信号)を受光する。尚、ビート信号光を生じさせる散乱光として、生体100に照射されたレーザ光LBの透過光に相当する前方散乱光が用いられてもよいし、生体100に照射されたレーザ光LBの反射光に相当する後方散乱光が用いられてもよい。
Thereafter, the
その後、受光素子12は、受光したビート信号光を電気信号に変換することで得られる検出電流を生成する。受光素子12は、生成した検出電流を、増幅器13に出力する。増幅器13は、受光素子12から出力される検出電流(つまり、受光素子12が受光したビート信号光に応じた検出電流)を、電圧信号に変換した上で増幅する。増幅器は、電圧信号をA/Dコンバータ14に出力する。
Thereafter, the
その後、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 /
演算回路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
尚、血流量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
その後、脈波周期検出回路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
具体的な数値を例にあげて説明する。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 /
尚、脈波周期検出回路16が算出するサンプル数nは、整数であることが好ましい。従って、脈波周期検出回路16は、算出したサンプル数nが小数点以下の桁を有している場合には、小数点以下の桁が繰り上げられてもよい。或いは、小数点以下の桁が繰り下げられてもよい。或いは、小数点以下の桁が四捨五入されてもよい。
Note that the number of samples n calculated by the pulse wave
図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
尚、脈波周期検出回路16は、脈波周期TrのN(但し、Nは、2以上の整数)倍の期間に対応する血流量Q(k)のサンプル数nを算出してもよい。この場合、脈波周期検出回路16は、n=Tr×N/Tq=Tr×N/(Nf×Ta)という数式を用いて、サンプル数nを算出してもよい。
Note that the pulse wave
脈波周期検出回路16は、算出したサンプル数nを、平均化回路17に出力する。
The pulse wave
その後、平均化回路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
より具体的には、図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
その後、サンプリング回路19は、平均化回路17が算出した平均血流量Qhを、タイミング生成器18が生成したタイミング(図3の5段目の波形に示す再サンプリングタイミングに相当)に同期して、血流量検出装置1の外部に(或いは、血流量検出装置1が内部に備える不図示の処理ブロックに)出力する(ステップS18)。その結果、血流量検出装置1からは、図3の6段目(最終段)の波形に示す平均血流量Qhが出力される。
Thereafter, the
尚、タイミング生成器18が生成したタイミングは、固定タイミング(つまり、固定的な周期で現れるタイミング)であることが好ましい。この場合、血流量検出装置1から出力される平均血流量Qhの出力レート(言い換えれば、転送レート)が一定になる。
It should be noted that the timing generated by the
また、サンプリング回路19は、平均血流量Qhを所定時間毎に更に平均化することで算出される時間平均値Qh’を出力してもよい。この場合、血流量検出装置1から出力する値(この場合には、時間平均値Qh’)の変動を少なくすることができる。
Further, the
ここで、比較例として、図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
しかるに、本実施例の血流量検出装置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
加えて、本実施例の血流量検出装置1によれば、比較例の血流量検出装置と比較して、上述した平均血流量Qhの誤差(つまり、平均血流量Qhを算出する際に用いられる血流量Q(k)の振幅中心に対する偏り)に起因した誤った揺らぎが平均血流量Qhに生ずる可能性が低減される。或いは、本実施例の血流量検出装置1によれば、上述した平均血流量Qhの誤差(つまり、平均血流量Qhを算出する際に用いられる血流量Q(k)の振幅中心に対する偏り)に起因した誤った揺らぎが平均血流量Qhに生ずることが殆ど又は全くなくなる。このため、本実施例の血流量検出装置1によれば、生体100の反応として本来の揺らぎが好適に検出される。
In addition, according to the blood
尚、上述した平均血流量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
加えて、本実施例の血流量検出装置1によれば、平均血流量Qhは、タイミング生成器18が生成したタイミングに同期して出力される。従って、血流量検出装置1から出力される平均血流量Qhの出力レート(言い換えれば、転送レート)を一定にすることができる。従って、血流量検出装置1と外部機器との間のインタフェース構築が比較的容易になる。
In addition, according to the blood
加えて、本実施例の血流量検出装置1は、血流量Q(k)が算出される周期Tqよりも長い周期で、平均化血流量Qhを出力することができる。このため、血流量Q(k)をそのまま出力する場合と比較して、単位時間当たりの出力データ量を低減することができる。更には、本実施例の血流量検出装置1は、血流量Q(k)が平均化された後に平均化した結果である平均血流量Qhが出力することができる。このため、血流量Q(k)をそのまま出力する場合と比較して、単位時間当たりの出力データ量を低減することができる。
In addition, the blood flow
(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
(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
ユーザ要求受付部21は、サンプリング回路19が平均血流量Qhを出力するタイミングを指示するユーザ(つまり、血流量検出装置2のユーザ)の指示を受け付ける。ユーザ要求受付部21は、受け付けたユーザの指示に同期して、サンプリング回路19が平均血流量Qhを出力するタイミングを生成する。
The user
(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
その後、第2実施例の血流量検出装置2では、サンプリング回路19は、平均化回路17が算出した平均血流量Qhを、タイミング生成器18が生成したタイミング(つまり、ユーザの指示に応じたタイミング)に同期して、血流量検出装置1の外部に出力する(ステップS21)。
Thereafter, in the blood flow rate detection device 2 of the second embodiment, the
このような第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
加えて、第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
(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
(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
脈拍数算出回路31は、脈波周期検出回路16が算出したサンプル数nから、脈拍数R(i(尚、iは、脈拍数の算出のタイミングに応じて時系列的に割り当てられる変数))を算出する。
The pulse
平均化回路32は、所定数w個の脈拍数R(i)の平均値(以下、適宜“平均脈拍数Rh”と称する)を算出する。
The averaging
サンプリング回路33は、平均化回路32が算出した平均脈拍数Rhを、タイミング生成器18が生成したタイミングに同期して、血流量検出装置1の外部に(或いは、血流量検出装置1が内部に備える不図示の処理ブロックに)出力する。
The
尚、血流量検出装置1は、平均化回路32が算出した平均脈拍数Rhをそのまま出力してもよい。この場合、血流量検出装置1は、サンプリング回路33を備えていなくともよい。
In addition, the blood flow
(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
第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
具体的な数値を例にあげて説明する。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 /
脈拍数算出回路31は、算出した脈拍数R(i)を、平均化回路32に出力する。
The pulse
その後、平均化回路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
その後、サンプリング回路33は、平均化回路32が算出した平均脈拍数Rhを、タイミング生成器18が生成したタイミングに同期して、血流量検出装置1の外部に(或いは、血流量検出装置1が内部に備える不図示の処理ブロックに)出力する(ステップS33)。
After that, the
このような第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
(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
(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
脈振幅検出回路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
レベル低下検出回路42は、脈振幅検出回路41が検出した脈振幅Aが所定閾値未満であるか否かを判定する。
The level
平均化回路43は、所定数p個のサンプル数n(j(尚、jは、サンプル数の算出のタイミングに応じて時系列的に割り当てられる変数))の平均値(以下、適宜“平均サンプル数Nh”と称する)を算出する。
The averaging
選択回路44は、レベル低下検出回路42での判定結果に基づいて、脈波周期検出回路16が検出するサンプル数n(j)及び平均化回路43が算出する平均サンプル値Nhのいずれかを、平均化回路17に出力する。
The
(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
その後、第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
その後、脈波振幅検出回路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
その後、レベル低下検出回路42は、脈波振幅検出回路41が検出した脈振幅Aが所定閾値(図11の2段目の波形参照)未満であるか否かを判定する(ステップS43)。尚、所定閾値は、脈波周期Trを好適に算出することが可能な脈振幅Aの下限値が設定されていてもよい。或いは、所定閾値は、信頼性のある脈波周期Trを算出することが可能な脈振幅Aの下限値が設定されていてもよい。
Thereafter, the level
ステップ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
他方で、ステップ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
その後、サンプリング回路19は、平均化回路17が算出した平均血流量Qhを、タイミング生成器18が生成したタイミングに同期して、血流量検出装置1の外部に(或いは、血流量検出装置1が内部に備える不図示の処理ブロックに)出力する(ステップS18)。
Thereafter, the
このような第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
加えて、第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
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.
前記算出手段は、(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.
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)
| 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)
| 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 |
-
2012
- 2012-03-22 WO PCT/JP2012/057379 patent/WO2013140583A1/en not_active Ceased
Patent Citations (5)
| 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)
| 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 |