1. Field of the Invention
The present invention relates to non-invasive measurement of blood sugar levels for measuring glucose concentration in a living body without blood sampling, and an optical measurement apparatus suitable therefor.
2. Description of Related Art
Hilson et al. report facial and sublingual temperature changes in diabetics following intravenous glucose injection (Diabete & Metabolisme, “Facial and sublingual temperature changes following intravenous glucose injection in diabetics” by R. M. Hilson and T. D. R. Hockaday, 1982, 8, 15-19). Scott et al. discuss the issue of diabetics and thermoregulation (Can. J. Physiol. Pharmacol., “Diabetes mellitus and thermoregulation”, by A. R. Scott, T. Bennett, I. A. MacDonald, 1987, 65, 1365-1376). Based on the knowledge gained from such researches, Cho et al. suggest a method and apparatus for determining blood glucose concentration by temperature measurement without requiring the collection of a blood sample (U.S. Pat. Nos. 5,924,996, and 5,795,305).
Various other attempts have been made to determine glucose concentration without blood sampling. For example, a method has been suggested (JP Patent Publication (Kokai) No. 2000-258343 A) whereby a measurement site is irradiated with near-infrared light of three wavelengths, and the intensity of transmitted light as well as the temperature of the living body is detected. A representative value of the second-order differentiated value of absorbance is then calculated, and the representative value is corrected in accordance with the difference between the living body temperature and a predetermined reference temperature. The blood sugar concentration corresponding to the thus corrected representative value is then determined. An apparatus is also provided (JP Patent Publication (Kokai) No. 10-33512 A (1998)) whereby a measurement site is heated or cooled while monitoring the living body temperature. The degree of attenuation of light based on light irradiation is measured at the moment of temperature change so that the glucose concentration responsible for the temperature-dependency of the degree of light attenuation can be measured. Further, an apparatus is reported (JP Patent Publication (Kokai) No. 10-108857 A (1998)) whereby an output ratio between reference light and transmitted light following the irradiation of the sample is taken, and then a glucose concentration is calculated in accordance with a linear expression of the logarithm of the output ratio and the living body temperature. Another apparatus for measuring glucose concentration is provided (U.S. Pat. No. 5,601,079) whereby the result of irradiation using two light sources is detected by three infrared light detectors and also temperature is detected.
Glucose (blood sugar) in blood is used for glucose oxidation reaction in cells to produce necessary energy for the maintenance of living bodies. In the basal metabolism state, in particular, most of the produced energy is converted into heat energy for the maintenance of body temperature. Thus, it can be expected that there is some relationship between blood glucose concentration and body temperature. However, as is evident from the way sicknesses cause fever, the body temperature also fluctuates due to factors other than blood glucose concentration. While methods have been proposed to determine blood glucose concentration by temperature measurement without blood sampling, they could hardly be considered sufficiently accurate.
Further, a method has also been proposed that detects the result of irradiation of light from two light sources using three infrared light detectors, and that also detects temperature for determining glucose concentration. The method, which only detects two kinds of optical intensity, is unable to provide sufficient accuracy.
It is an object of the invention to provide a method and apparatus for determining blood glucose concentration with high accuracy based on temperature data and optical data of a test subject without blood sampling.
Blood sugar is delivered to the cells throughout the human body via blood vessel systems, particularly the capillary blood vessels. In the human body, complex metabolic pathways exist. Glucose oxidation is a reaction in which, fundamentally, blood sugar reacts with oxygen to produce water, carbon dioxide, and energy. Oxygen herein refers to the oxygen delivered to the cells via blood. The volume of oxygen supply is determined by the blood hemoglobin concentration, the hemoglobin oxygen saturation, and the volume of blood flow. On the other hand, the heat produced in the body by glucose oxidation is dissipated from the body by convection, heat radiation, conduction, and so on. On the assumption that the body temperature is determined by the balance between the amount of energy produced in the body by glucose burning, namely heat production, and heat dissipation such as mentioned above, the inventors set up the following model:
(1) The amount of heat production and the amount of heat dissipation are considered equal.
(2) The amount of heat production is a function of the blood glucose concentration and the volume of oxygen supply.
(3) The volume of oxygen supply is determined by the blood hemoglobin concentration, the blood hemoglobin oxygen saturation, and the volume of blood flow in the capillary blood vessels.
(4) The amount of heat dissipation is mainly determined by heat convection and heat radiation.
According to this model, we achieved the present invention after realizing that blood sugar levels can be accurately determined on the basis of the results of measuring the temperature of the body surface and parameters relating to the blood oxygen concentration and the blood flow volume. The parameters can be measured, e.g., from a part of the human body, such as the fingertip. The parameters relating to convection and radiation can be determined by measuring the temperature on the fingertip. The parameters relating to the blood hemoglobin concentration and the blood hemoglobin oxygen saturation can be determined by spectroscopically measuring blood hemoglobin and then finding the ratio between hemoglobin bound with oxygen and hemoglobin not bound with oxygen. The parameter relating to the volume of blood flow can be determined by measuring the amount of heat transfer from the skin.
The invention provides an optical measurement apparatus comprising: a first light source for producing light of a first wavelength; a first optical fiber for irradiating a light incident point on the surface of a subject with the light from said first light source; a second light source for producing light of a second wavelength; a second optical fiber for irradiating said light incident point on the surface of the subject with the light from said second light source in a direction different from that of the light of said first wavelength; a first photodetector on which reflected light of the light of said first wavelength reflected by said light incident point and scattered light of the light of said second wavelength are incident without via fiber; a second photodetector on which reflected light of the light of said second wavelength reflected at said light incident point and scattered light of the light of said first wavelength are incident without via fiber; a third photodetector; and a third optical fiber having an incident end thereof disposed at such a position as to be in contact with the surface of said subject, said third optical fiber being adapted to receive, on an incident end thereof, light exiting from an area spaced apart from said light incident point on the surface of said subject and then transmit the light to said third detector.
Preferably, the plane of incidence of the light of the first wavelength and the plane of incidence of the light of the second wavelength are substantially perpendicular to each other with respect to the light incident point on the subject surface. The plane of incidence herein refers to a plane that includes the incident ray and a normal at the incident point on the subject surface. Further, in the present specification, the ray that enters the incident plane after having been irradiated onto the incident point on the subject surface will be referred to as reflected light. The light that leaves in directions other than that of the incident plane from near the incident point will be referred to as scattered light. The scattered light that leaves out of a position on the subject surface that is spaced apart from the incident point will be referred to as traveled photon.
Preferably, the outgoing light from each light source is irradiated onto the light incident point on the subject surface via an optical fiber. The reflected light and scattered light from the examined subject are directly incident on the photodetector, and the traveled photon is incident on the photodetector via an optical fiber. An exiting end of the light-irradiating optical fiber and an incident end of the optical fiber for detecting reflected or scattered are preferably disposed near the plane of a cone whose apex corresponds to the light incident point on the subject surface. The first wavelength may be a wavelength at which the molar absorption coefficient of oxyhemoglobin is equal to that of reduced hemoglobin, and the second wavelength may be a wavelength for detecting the difference in absorbance between the oxyhemoglobin and reduced hemoglobin.
The invention further provides a blood sugar level measuring apparatus comprising: (1) a heat-amount measuring portion for measuring a plurality of temperatures deriving from a body surface and acquiring information that is used for calculating a convective heat transfer amount and radiation heat transfer amount that are related to the dissipation of heat from the body surface; (2) a blood flow volume measuring portion for acquiring information about the volume of blood flow; (3) an optical measurement portion including a light source for producing light of at least two different wavelengths, an optical system for irradiating the body surface with the light emitted by said light source, and at least three different photodetectors for detecting the light that has been irradiated onto the body surface, said optical measurement portion providing hemoglobin concentration and hemoglobin oxygen saturation in blood; (4) a memory portion in which relationships between parameters respectively corresponding to said plurality of temperatures, blood flow volume, hemoglobin concentration and hemoglobin oxygen saturation in blood, and blood sugar levels are stored; (5) a calculation portion for converting a plurality of measurement values inputted from said heat amount measuring portion, said blood flow volume measuring portion, and said optical measurement portion respectively into said parameters, and then calculating a blood sugar value by applying said parameters to said relationships stored in said memory portion; and (6) a display portion for displaying the blood sugar level calculated by said calculation portion, wherein said optical measurement portion includes a first light source for producing light of a first wavelength, a second light source for producing light of a second wavelength, a first optical fiber, a second optical fiber, a third optical fiber, a first photodetector, a second photodetector, and a third photodetector, a light incident point on the surface of a subject is irradiated with light emitted by said first light source via said first optical fiber, said light incident point on the surface of said subject is irradiated with light emitted by said second light source via said second optical fiber in a direction different from that of the light of said first wavelength, reflected light of the light of said first wavelength reflected at said light incident point and scattered light of the light of said second wavelength are incident on said first photodetector without via fiber, reflected light of the light of said second wavelength reflected at said light incident point and scattered light of the light of said first wavelength are incident on said second photodetector without via fiber, said third optical fiber has an incident end thereof disposed at such a position as to be in contact with the surface of said subject in an area spaced apart from said light incident point on the subject surface, and said third photodetector is adapted to receive light exiting from an area spaced apart from said light incident point on the surface of the subject via said third optical fiber.
The invention will now be described by way of preferred embodiments thereof with reference made to the drawings.
Initially, the above-mentioned model will be described in more specific terms. Regarding the amount of heat dissipation, convective heat transfer, which is one of the main causes of heat dissipation, is related to temperature difference between the ambient (room) temperature and the body-surface temperature. The amount of heat dissipation due to radiation, which is another main cause of dissipation, is proportional to the fourth power of the body-surface temperature according to the Stefan-Boltzmann law. Thus, it can be seen that the amount of heat dissipation from the human body is related to the room temperature and the body-surface temperature. On the other hand, the amount of oxygen supply, which is a major factor related to the amount of heat production, is expressed as the product of hemoglobin concentration, hemoglobin oxygen saturation, and blood flow volume.
The hemoglobin concentration can be measured from the absorbance at the wavelength (equal-absorbance wavelength) at which the molar absorbance coefficient of the oxyhemoglobin is equal to that of the reduced (deoxy-) hemoglobin. The hemoglobin oxygen saturation can be measured by measuring the absorbance at the equal-absorbance wavelength and the absorbance at at least one different wavelength at which the ratio between the molar absorbance coefficient of the oxyhemoglobin and that of the reduced (deoxy-) hemoglobin is known, and then solving simultaneous equations. Namely, the hemoglobin concentration and hemoglobin oxygen saturation can be obtained by conducting the measurement of absorbance at at least two wavelengths. However, in order to accurately determine the hemoglobin concentration and hemoglobin oxidation saturation from absorbance, the influence of interfering components must be corrected. The interfering components affecting the absorbance include the thickness of the skin (epidermis), for example. These interfering components can be measured in various manners, of which one example will be described below.
The thickness of the skin can be measured by measuring the intensity of only that light that has traveled in the skin by a distance d from where light was shone on the skin.
The measurements are carried out using at least three detectors, namely a reflected-light detector for detecting mainly reflected light, a scattered-light detector for detecting mainly scattered light, and a traveled-photon detector for detecting traveled photon.
The reflected-light detector can detect part of the scattered light produced by the light traveling inside the body and then exiting from the body surface, as well as mainly the reflected light reflected by the body surface. The scattered-light detector can detect part of the scattered light scattered on the body surface, as well as mainly the scattered light produced by the light passing inside the body and then exiting through the body surface. The path of the traveled photon up to the traveled-photon detector is optically blocked in order to prevent the detection of light other than the traveled photon, namely the light deriving from reflected light and scattered light. The traveled-photon detector is thus adapted to detect only traveled photon, so that the skin thickness can be estimated. During detection, a total of at least three detectors, namely at least one each of the reflected-light detector, scattered-light detector, and traveled-photon detector, are used. Preferably, additional detectors with similar functions and high detection sensitivities adapted for particular kinds of wavelength may be used. Further, a transmitted-light detector may be added for detecting light that has passed through the detection area, as necessary.
The wavelength values described herein are most appropriate values for obtaining absorbance for various intended purposes, such as for obtaining the absorbance at the equal molar absorbance coefficient wavelengths, or for obtaining the peak of absorbance. Thus, in addition to the wavelengths described herein, other wavelengths in the vicinities thereof may be used and still similar measurements can be performed.
The rest is the blood flow volume, which can be measured by various methods. One example will be described below.
Before the block comes into contact with the body surface, the temperatures T1 and T2 at the two points of the block are equal to the room temperature Tr. When a body-surface temperature Ts is higher than the room temperature Tr as the block comes into contact with the body surface, the temperature T1 swiftly rises due to the transfer of heat from the skin, and it approaches the body-surface temperature Ts. On the other hand, the temperature T2 is lowered from the temperature T1 as the heat conducted through the block is dissipated from the block surface, and it rises more gradually. The temporal variation of the temperatures T1 and T2 depends on the amount of heat transferred from the body surface to the block, which in turn depends on the blood flow volume in the capillary blood vessels under the skin. If the capillary blood vessels are regarded as a heat exchanger, the coefficient of transfer of heat from the capillary blood vessels to the surrounding cell tissues is given as a function of the blood flow volume. Thus, by measuring the amount of heat transfer from the body surface to the block by monitoring the temporal variation of the temperatures T1 and T2, the amount of heat transferred from the capillary blood vessels to the cell tissues can be estimated. Based on this estimation, the blood flow volume can then be estimated.
Then, the T1 measured value between tstart and tend is approximated by an S curve, such as a logistic curve. A logistic curve is expressed by the following equation:
where T is temperature, and t is time.
The measured value can be approximated by determining coefficients a, b, c, and d using the non-linear least-squares method. For the resultant approximate expression, T is integrated between time tstart and time tend to obtain a value S1.
Similarly, an integrated value S2 is calculated from the T2 measured value. The smaller (S1−S2) is, the larger the amount of transfer of heat is from the body surface to the position of T2. (S1−S2) becomes larger with increasing body-surface contact time tCONT (=tend−tstart). Thus, a5/(tCONT×(S1−S2)) is designated as a parameter X5 indicating the volume of blood flow, using a5 as a proportionality coefficient.
Thus, it will be seen that the measured amounts necessary for the determination of blood glucose concentration by the above-described model are the room temperature (ambient temperature), body surface temperature, temperature changes in the block brought into contact with the body surface, the temperature due to radiation from the body surface, the absorbance of reflected light or scattered light at at least two wavelengths, and the intensity of traveled photon.
Hereafter, an example of an apparatus for non-invasively measuring blood sugar levels according to the principle of the invention will be described.
On the top surface of the apparatus are provided an operating portion 1, a measuring portion 12 where the finger to be measured is to be placed, and a display portion 13 for displaying measurement results, the state of the apparatus, measured values, for example. The operating portion 11 includes four push buttons 11a to 11d for operating the apparatus. The measuring portion 12 has a cover 14 which, when opened (as shown), reveals a finger rest portion 15 with an oval periphery. The finger rest portion 15 accommodates an opening end 16 of a radiation-temperature sensor portion, a contact-temperature sensor portion 17, and an optical sensor portion 18.
First, the process of measuring temperatures by the non-invasive blood sugar level measuring apparatus according to the invention will be described. In a portion of the measuring portion with which the examined portion (ball of the finger) is to come into contact, a thin plate 21 of a highly heat-conductive material, such as gold, is placed. A bar-shaped heat-conductive member 22, which is made of a material with a heat conductivity lower than that of the plate 21, such as polyvinylchloride, is thermally connected to the plate 21 and extends into the apparatus. The temperature sensors include a thermistor 23 that is an adjacent-temperature detector with respect to the examined portion for measuring the temperature of the plate 21, and a thermistor 24 that is an indirect-temperature detector with respect to the examined portion for measuring the temperature of a portion of the heat-conducting member which is spaced apart from the plate 21 by a certain distance. An infrared lens 25 is disposed inside the apparatus at such a position that the examined portion (ball of the finger) placed on the finger rest portion 15 can be seen through the lens. Below the infrared lens 25 is disposed a pyroelectric detector 27 via an infrared radiation-transmitting window 26. Another thermistor 28 is disposed in close proximity to the pyroelectric detector 27.
Thus, the temperature sensor portion of the measuring portion has four temperature sensors, and they measure four kinds of temperatures as follows:
(1) Temperature on the finger surface (thermistor 23): T1
(2) Temperature of the heat-conducting member (thermistor 24): T2
(3) Temperature of radiation from the finger (pyroelectric detector 27): T3
(4) Room temperature (thermistor 28): T4
The optical sensor portion 18 is described hereafter. The optical sensor portion 18 measures the hemoglobin concentration and the hemoglobin oxygen saturation necessary for the determination of the oxygen supply volume. In order to measure the hemoglobin concentration and the hemoglobin oxygen saturation, it is necessary to measure the absorbance of scattered light at at least two wavelengths, the absorbance of reflected light at at least one wavelength, and the intensity of traveled photon at at least one wavelength. The accuracy of the absorbance of reflected light can be improved by measuring at a plurality of wavelengths, if possible, and then using a mean value. Thus, in the present embodiment, the absorbance of reflected light is measured at two different wavelengths. The accuracy of the measurement of the intensity of traveled photon can also be improved by measuring at a plurality of wavelengths, if possible, and then using a mean value.
The ends of three optical fibers 31 to 33 are located in the optical sensor portion 18. The optical fibers 31 and 32 are for optical irradiation, while the optical fiber 33 is for receiving light. As shown in
The two light-emitting diodes 36 and 37 emit light in a time-sharing manner. The finger of an examined subject is irradiated with the light emitted by the light-emitting diodes 36 and 37 via the irradiating optical fibers 31 and 32. The light shone on the finger from the light-irradiating optical fiber 31 is reflected by the skin, and the reflected light is detected by the photodiode 38, the scattered light is detected by the photodiode 39, and the traveled photon is incident on the light-receiving optical fiber 33 and is then detected by the photodiode 40. The traveled photon-receiving optical fiber 33 is adapted to be in close contact with the finger surface such that it can avoid the direct entry of reflected and/or scattered light. The light with which the finger is irradiated via the light-irradiating optical fiber 32 is reflected by the skin of the finger, and the reflected light is detected by the photodiode 39, the scattered light is detected by the photodiode 38, and the traveled photon is incident on the light-receiving optical fiber 35 and is then detected by the photodiode 40. Thus, by causing the two light-emitting diodes 36 and 37 to emit light in a time-divided manner, the photodiodes 38 and 39 can detect different light depending on the irradiating position of the light-irradiating optical fibers. By this structure, the size of the optical sensor portion 18 can be reduced. With regard to the light shone on the finger via the light-irradiating fiber 32, the light-receiving optical fiber 33 may be adapted not to detect traveled photon.
By thus employing the structure such that the light reflected by the skin of the finger and the scattered light are received directly by photodiodes, the received amount of light detected by each photodiode can be increased. Regarding the light-receiving optical fiber 33, it is similarly possible to increase the amount of received light by disposing the photodiode 40 directly at the position corresponding to the end of the light-receiving optical fiber 33. However, putting the photodiode 40 directly at the end of the light-receiving optical fiber 33 would result in an increased size of the optical sensor portion 18. Accordingly, it is desirable to use the light-receiving optical fiber 33 if the size of the optical sensor portion 18 is to be further reduced.
(1) Sends a control signal 1 to a controller 2 in synchronism with the clock of a clock generator for a certain duration of time in order to select a control signal 2. As a result, a switching circuit 51 is turned on, thereby turning power on and causing the light-emitting diode 36 to emit light.
(2) After a certain duration of time has elapsed, sends a control signal 1 to the controller 2 in synchronism with the clock of the clock generator for a certain duration of time in order to select a control signal 3. As a result, a switching circuit 52 is turned on, thereby turning power on and causing the light-emitting diode 37 to emit light.
It is also possible to cause the two light-emitting diodes 36 and 37 to emit nearly simultaneously, rather than in a time-divided manner. In an example of a method of separately detecting light from a plurality of light sources, the individual light sources are driven by modulating them with different modulation frequencies. In this method, light from each light source can be separately detected by focusing on the frequency components contained in a detection signal from photodetectors.
The arrangement of the light-irradiating optical fibers, photodiodes, and the light-receiving optical fibers in the optical sensor portion 18 is determined based on the following theories (1) to (3).
(1) Regarding the position of the reflected-light receiving photodiode with respect to the light-irradiating optical fiber, it is most appropriate to position the light-receiving plane of the photodetector at a position where the reflected light is theoretically received, namely at a position within the plane of incidence of light on the subject where light reflected in a direction with an outgoing angle that is equal to the angle of incidence on a light-incident point of the subject is received. By locating the light-receiving plane of the reflected-light receiving photodiode at such a position, the ratio of reflected light in the amount of received light can be maximized.
(2) The scattered light-receiving photodetector has the light-receiving plane thereof disposed in a plane that forms an angle of approximately 90° with respect to the plane of incidence of light on the subject. The scattered-light receiving photodetector is disposed at approximately 90° relative to the reflected-light receiving photodetector because the source of light detected as scattered light is desired to be narrowed to the scattering phenomena as much as possible, as opposed to theory (1), or because the range of the phenomena as the object of detection of scattering is desired to be increased by the provision of the large angle of approximately 90°.
(3) The light-receiving end of the light-receiving optical fiber for traveled photon is disposed at a position in the plane of incidence of light on the subject that is farther than the light-receiving plane of the reflected-light receiving photodetector with respect to the light-irradiating optical fiber. The light-receiving end of the traveled-photon receiving optical fiber is thus disposed in the plane of incidence of light on the subject for the following reason. During the process in which light enters the skin and is scattered inside, the distribution of light spreads, and yet the distribution is greatest in the direction of incidence. As a result, the amount of light exiting from the skin is also greatest in this direction, so that the traveled photon can be most efficiently detected. Further, the light-receiving end of the traveled-photon receiving optical fiber is disposed farther than the light-receiving end of the reflected-light receiving photodetector with respect to the light-irradiating optical fiber. By so doing, a large amount of information can be detected, such as information relating to the absorption of light by hemoglobin in blood flowing in the capillary blood tubes during the process of light penetrating the skin and being scattered inside, or information relating to the thickness of skin, for example. It is also possible, however, to dispose the light-receiving optical fiber for traveled photon at a position other than that in the plane of incidence of light on the subject, though in that case the amount of traveled photon that is detected would be reduced.
In accordance with those theories (1) to (3), the exiting ends of the light-irradiating optical fibers, the photodetectors, and the receiving end of the light-receiving optical fiber are disposed in the optical sensor portion 18 as shown in the plan view of
Regarding the angles of irradiation and detection of light by the light-irradiating optical fiber 31 and photodiode 38, the reflected-light receiving photodiode 38 is positioned such that it can receive a beam of light reflected at a point (light incident point on the subject) y, namely the point y in
By thus disposing the light-irradiating optical fiber 31, the photodiode 38 and the traveled-photon receiving optical fiber 33 along the same line as shown in
Alternatively, the exiting ends of the light-irradiating optical fibers, the photodiodes, and the receiving end of the light-receiving optical fiber may be disposed in the optical sensor portion 18 as shown in a plan view of
By thus disposing the traveled-photon receiving optical fiber 33 on the line ZZ as shown in
Further alternatively, the exiting ends of the light-irradiating optical fibers, the photodiodes, and the receiving end of the light-receiving optical fiber in the optical sensor portion 18 may be disposed as shown in a plan view of
In this arrangement, the traveled-photon receiving optical fiber 33 is positioned in a direction opposite to that in which the light-irradiating optical fibers 31 and 32 radiate, so that, although the amount of light received by the light-receiving optical fiber 33 is fairly small, the received light contains hardly any reflected light or scattered light and consists mostly of traveled photon.
Regarding the arrangement of the light-irradiating optical fibers, photodiodes, and light-receiving optical fiber in the optical sensor portion 18 shown in
Regarding the optical sensor portion 18, the exiting ends and the receiving planes of the light-irradiating optical fibers 31 and 32 and the photodiodes 38 and 39, respectively, may be displaced a little in their optical axial directions as long as they are aimed at the light incident point γ on the subject (see
Further regarding the optical sensor portion 18 shown in
The photodiodes 38 and 39 provide reflectance R as measurement data, and absorbance can be approximately calculated from log(1/R). Light of wavelengths 810 nm and 950 nm is irradiated, and R is measured for each and log(1/R) is obtained for each, so that absorbance AD11 and AD21 at wavelength 810 nm and absorbance AD12 and AD22 at wavelength 950 nm can be measured. Part of the light penetrates into the skin and travels a certain distance d while being scattered therein repeatedly. The intensity ID3i of traveled photon is measured by a photodiode 40. (The absorbance of reflected light of wavelength λi detected by the photodiode for detecting reflected light is referenced by AD1i, the absorbance of scattered light of wavelength λi detected by the photodiode for detecting scattered light is referenced by AD2i, and the intensity of traveled photon of wavelength λi detected by the photodiode 40 is referenced by ID3i.)
When the reduced hemoglobin concentration is [Hb] and the oxyhemoglobin concentration is [HbO2], scattered-light absorbance AD2i at wavelength λi is expressed by the following equations:
where AHb(λi) and AHb02(λi) are the molar absorbance coefficients of the reduced hemoglobin and the oxyhemoglobin, respectively, and are known at the respective wavelengths. Terms a, b, and c are proportionality coefficients. AD1i is the reflected-light absorbance at wavelength λi, and ID3i is the traveled photon intensity at wavelength λi. From the above equations, the parameter aR, which is determined by the relationship between reflected light and scattered light, and the parameter D of the skin thickness can be determined as constants, and can be substituted in the equation of AD2i. The parameter determined by the relationship between reflected light and scattered light is a parameter relating to the roughness of the skin surface, for example, and the influence of the roughness of the skin surface, for example, can be corrected using that parameter. The parameter relating to the thickness of the skin can be determined from the measurement value obtained by the traveled-photon detector, and the influence of the thickness of the skin can be corrected using that parameter. Since i=2 wavelengths, two equations of AD2i are produced. By solving these simultaneous equations, the two variables to be obtained, namely [Hb] and [HbO2], can be obtained. The hemoglobin concentration [Hb]+[HbO2], and the hemoglobin oxygen saturation [HbO2]/([Hb]+[HbO2]) can be determined from the above-obtained [Hb] and [HbO2].
Although the present example has been described with regard to the measurement of the hemoglobin concentration and the hemoglobin oxygen saturation based on the measurement of absorbance at two wavelengths, absorbance may be measured by adding one or more wavelengths at which the difference in molar absorbance coefficient between the oxyhemoglobin and the deoxyhemoglobin is large so as to increase the measurement accuracy.
For example, when six wavelengths are used for measurement, any of the configurations shown in
The light-emitting diode 36a emits light of 810 nm, light-emitting diode 36b light of 880 nm, light-emitting diode 36c light of 950 nm, light-emitting diode 37a light of 450 nm, light-emitting diode 37b light of 520 nm, and light-emitting diode 37c light of 660 nm, for example. Using the result of detection of irradiated light having these six wavelengths, corrections can be made for the influences of interfering components on the determination of hemoglobin concentration and hemoglobin oxygen saturation from absorbance, the interfering components including melanin pigment, bilirubin and the turbidity of blood, for example. Thus, the accuracy of measurement can be improved.
As mentioned above, by providing the light-irradiating optical fiber with three branches, an ideal configuration can be obtained in which the three light-emitting diodes share the same point of light irradiation. However, since the actual fiber-irradiated light has certain spread, the same function can be provided by employing a structure in which a light-irradiating fiber is provided to each light-emitting diode and the tips of the fibers are bundled. This structure, which can employ conventional fibers and can therefore be made inexpensively, is shown in
In practice, the light-emitting diodes 36 and 37 and the photodiode 40 are equipped with a light-blocking cap 42 for preventing the leakage of light to the outside and the reception of light from the outside. The light-blocking cap 42 is formed by a soft material, such as a silicon resin, so that it can be easily mounted on the light-emitting diode or optical fiber assembly. With regard to the photodiodes 38 and 39, which are disposed inside the sensor portion, there might be no need to provide such a measure because they are disposed inside the sensor portion and are therefore already blocked against external light.
The seven kinds of analog signals are supplied via individual amplifiers A1 to A7 to analog/digital converters AD1 to AD7, where they are converted into digital signals. Based on the digitally converted values, parameters xi (i=1, 2, 3, 4, 5) are calculated. The following are specific descriptions of xi (where e1 to e5 are proportionality coefficients):
Parameter proportional to heat radiation
x1=e1×(T3)4
Parameter proportional to heat convection
x2=e2×(T4−T3)
Parameter proportional to hemoglobin concentration
x3=e3×([Hb]+[Hb2])
Parameter proportional to hemoglobin saturation
Parameter proportional to blood flow volume
Then, normalized parameters are calculated from mean values and standard deviations of parameter xi obtained from actual data pertaining to large numbers of able-bodied people and diabetic patients. A normalized parameter Xi (where i=1, 2, 3, 4, 5) is calculated from each parameter xi according to the following equation:
where
xi: parameter
{overscore (x)}i: mean value of the parameter
SD(xi): standard deviation of the parameter
Using the above five normalized parameters, calculations are conducted for conversion into a glucose concentration to be eventually displayed. A program necessary for the processing calculations is stored in a ROM in the microprocessor built inside the apparatus. The memory area required for the processing calculations is secured in a RAM similarly built inside the apparatus. The results of calculation are displayed on the LCD.
The ROM stores, as a constituent element of the program necessary for the processing calculations, a function for determining glucose concentration C in particular. The function is defined as follows. C is expressed by the below-indicated equation (1), where ai (i=0, 1, 2, 3, 4, 5) is determined from a plurality of pieces of measurement data in advance according to the following procedure:
(1) A multiple regression equation is created that indicates the relationship between the normalized parameters and the glucose concentration C.
(2) Normalized equations (simultaneous equations) relating to the normalized parameters are obtained from equations obtained by the least-squares method.
(3) Values of coefficient ai (i=0, 1, 2, 3, 4, 5) are determined from the normalized equations and then substituted into the multiple regression equation.
Initially, the regression equation (1) indicating the relationship between the glucose concentration C and the normalized parameters X1, X2, X3, X4, and X5 is formulated.
Then, the least-squares method is employed to obtain a multiple regression equation that would minimize the error with respect to a measured value Ci of glucose concentration according to an enzyme electrode method. When the sum of squares of the residual is E, E is expressed by the following equation (2):
The sum E of squares of the residual becomes minimum when partial differentiation of equation (2) with respect to a0, a2, . . . , a5 gives zero. Thus, we have the following equations:
When the mean values of C and X1 to X5 are Cmean and X1mean to X5mean, respectively, since Ximean=0 (i=1 to 5), equation (1) yields:
The variation and covariation between the normalized parameters are expressed by equation (5). Covariation between the normalized parameter Xi (i=1 to 5) and C is expressed by equation (6).
Substituting equations (4), (5), and (6) into equation (3) and rearranging yields a set of simultaneous equations (normalized equations) (7). Solving the set of equations (7) yields a1 to a5.
a1S11+a2S12+a3S13+a4S14+a5S15=S1C
a1S21+a2S22+a3S23+a4S24+a5S25=S2C
a1S31+a2S32+a3S33+a4S34+a5S35=S3C
a1S41+a2S42+a3S43+a4S44+a5S45=S4C
a1S51+a2S52+a3S53+a4S54+a5S55=S5C (7)
Constant term a0 is obtained by means of equation (4). The thus obtained ai (i=0, 1, 2, 3, 4, 5) is stored in ROM at the time of manufacture of the apparatus. In actual measurement using the apparatus, the normalized parameters X1 to X5 obtained from the measured values are substituted into regression equation (1) to calculate the glucose concentration C.
Hereafter, an example of the process of calculating parameter Xi will be described. The example concerns measurement values obtained from able-bodied persons. Coefficients for the parameter calculation equations are determined by temperature data and optical measurement data that have been measured in advance. The ROM in the microprocessor stores the following formula for the calculation of the parameter:
When T3=36.5° C. is substituted in the above equations as a measurement value, for example, x1=1.74×103. When T4=19.7° C. is substituted in the above equations, x2=2.08×10. Then, before finding X3, it is necessary to find [Hb] and [HbO2]. The coefficients for a concentration calculation formula are determined by the scattered-light absorbance coefficient of each substance that has been measured in advance. Using that formula, [Hb] and [HbO2] can be determined by solving the following set of simultaneous equations in the case of measurement using two wavelengths:
Solving this set of simultaneous equations gives [Hb]=0.17 mmol/L and [HbO2]=2.17 mmol/L. Thus we have x3=3.18 and x4=2.48. Then, substituting S1=1.76×102, S2=1.89×10, and tCONT=22 seconds gives x5=4.40×102.
The hemoglobin concentration ([Hb]+[HbO2]) was calculated to be 2.34 mmol/L. When the hemoglobin concentration was measured at the same time by an invasive method, i.e. by blood sampling, the value was 2.28 mmol/L.
When the traveled photon is not similarly detected by the light-receiving optical fiber 33 at the same time, the information about the parameter of the thickness of the skin would not be obtained. In that case, the below-indicated simultaneous equations would be obtained, and solving them would yield [Hb]=0.18 mmol/L and [HbO2]=2.26 mmol/L. Thus, the hemoglobin concentration ([Hb]+[HbO2]) would be 2.44 mmol/L.
Thus, it has been confirmed that the result of calculation in the case where traveled photon is detected by the light-receiving optical fiber 33 is closer to the value of hemoglobin concentration measured by blood sampling than the calculation result in the case where traveled photon is not detected by the light-receiving optical fiber 33. Thus, it has been shown that the measurement accuracy can be improved by providing the optical sensor portion 18 with the light-receiving optical fiber 33.
Next, X1 to X5 are obtained. X1 to X5 are the results of normalization of the above-obtained parameters X1 to X5. Assuming the distribution of a parameter is normal, 95% of a normalized parameter takes on values between −2 and +2. The normalized parameters can be determined by the following equations:
From the above equations, we have normalized parameters X1=−0.06, X2=+0.04, X3=+0.05, X4=−0.12, and X5=+0.10.
Hereafter, an example of the process of calculating the glucose concentration will be described. The coefficients for regression equation (1) are determined in advance based on many items of data obtained from able-bodied persons and diabetics, and the ROM in the microprocessor stores the following formula for calculating the glucose concentration:
C=99.1+18.3×X1−20.2×X2−24.4×X3−21.8×X4−25.9×X5
Substituting X1 to X5 into the above equation gives C=96 mg/dl. Substituting normalized parameters X1=+1.15, X2=−1.02, X3=−0.83, X4=−0.91, and X5=−1.24, which can be obtained as an example of the measured values for a diabetic patient, in the equation yields C=213 mg/dl.
The following describes the results of measurement by the conventional enzymatic electrode method in which a blood sample is reacted with a reagent and the amount of resultant electrons is measured to determine glucose concentration, and the results of measurement by an embodiment of the invention. When the glucose concentration for an able-bodied person was 89 mg/dl according to the enzymatic electrode method in one example, substituting the normalized parameters X1=−0.06, X2=+0.04, X3=+0.07, X4=−0.10, and X5=+0.10, which were obtained by measurement at the same time according to the invention, into the above equation yields C=95 mg/dl. In another example, when the measured value of glucose concentration for a diabetic patient was 238 mg/dl according to the enzymatic electrode method, substituting the normalized parameters X1=+1.15, X2=−1.02, X3=−0.86, X4=1.02, and X5=−1.24, which were obtained by measurement at the same time according to the invention, into the above equation yields C=216 mg/dl. The results thus indicated that the method according to the invention can provide highly accurate glucose concentration values.
Thus, the invention makes it possible to determine blood sugar levels in a non-invasive measurement with similar levels of accuracy to the conventional invasive method.
All publication, patents, and patent applications cited herein are incorporated herein by reference to their entirety.
U.S. patent application Ser. No. 10/620,689 is a co-pending application of this application. The disclosures of the co-pending application are incorporated herein by cross-reference.