1. Introduction
Non-invasive blood glucose monitoring plays a crucial role in the management of diabetes and patient health. Near-infrared spectroscopy emerges as a promising technique for non-invasive glucose measurement, primarily due to its strong tissue penetration capability [
1]. This enables the acquisition of glucose information from both tissue fluid and blood within the dermis layer. The methodology and associated hardware systems have relatively matured. However, similar to other optical measurement techniques, in vivo glucose signals captured through near-infrared light are susceptible to interference from physiological variations and other components present in tissues [
2,
3]. The human body, being a dynamic biological system, presents intricate physiological processes and regulatory mechanisms that pose significant challenges in accurately extracting glucose information from optical signals.
An understanding of the magnitude and characteristics of the sensitivity caused by a unit change in blood glucose concentration is necessary. Developers of instruments are required to design the hardware resolution of the measurement system based on this sensitivity, as well as select the wavelengths of light and the source–detector separations (SDSs) to be used and develop appropriate signal processing approaches.
As directly obtaining the signal of glucose for living tissues is typically challenging, most research groups combine near-infrared diffuse transmittance/reflection spectroscopy with chemometric methods to extract in vivo blood glucose information [
4,
5,
6]. These approaches largely rely on data correlation or projection [
7], extracting glucose information from a statistical standpoint without quantifying glucose sensitivity based on physical or chemical principles.
Only a few groups focus on the direct sensing of glucose under discrete single wavelengths. Han et al. of our research team has made significant progress in in vivo signal acquisition strategy both in sensors and differential data processing, and they have succeeded in observing glucose signals directly [
8]. They measured the glucose in human arm skin at six wavelengths (within the 1000–1700 nm range) using differentiating processing on the diffuse reflectance of two SDSs to eliminate common-mode interference. However, they still cannot observe the glucose signals through the diffuse reflectance at the singular SDS. In this paper, we draw on Han’s method and utilize differential measurement techniques to estimate the changes in skin scattering coefficients caused by glucose. Then, we indirectly calculate the glucose sensitivity under a different single SDS. This differential measurement approach aids in obtaining reliable in vivo measurement results. The glucose signals derived from differential measurements exhibit distinct magnitudes and spectral characteristics compared to those obtained at a single SDS. Only by understanding the sensitivity characteristics at different SDSs can we flexibly choose the method and corresponding SDSs according to the measurement needs.
Many studies have employed simulation techniques, particularly the MC method, to identify the magnitude and characteristics of the sensitivity of glucose signals in human tissues. This approach has seen rapid development in recent years and is now capable of simulating the absorption and scattering processes of photons in highly complex human tissue structures [
9,
10]. However, this modeling and simulation approach requires the determination of the change of both absorption and scattering coefficients in living tissues in response to blood glucose variations. Changes in glucose lead to corresponding changes in the tissue’s absorption and scattering characteristics. Among these effects, the absorption change is relatively more straightforward, mainly due to the distinct absorption properties of glucose molecules at specific NIR wavelengths [
11], while the scattering is more complex, being closely associated with the state of the scattering media and susceptible to the tissue structure and composition [
12]. Therefore, figuring out the change in the scattering of living tissues caused by glucose is essential for acquiring an accurate evaluation of in vivo glucose signal sensitivity.
Attempts have been made to use biological phantoms or tissue slices as substitutes for live animal tissues to study the scattering characteristics caused by glucose changes [
13]. Experiments have shown that the scattering changes induced by glucose in phantoms or slice experiments are noticeably different from those in live animal experiments. In living tissues, the variability of tissue scattering coefficients caused by glucose may be significantly higher than those observed in similar experiments with slice tissues or phantom solutions [
14,
15,
16]. Matti Kinnunen et al. conducted a study comparing the scattering changes caused by glucose in 5% intralipid solution, tissue slices, and live mice. Using 910 nm light, a glucose concentration change of 1 mmol/L resulted in only a −0.07% change in the scattering coefficient in intralipid solution compared to −1.37% and −2.7% changes in sliced skin and live mice, respectively [
15]. This suggests that the changes in scattering coefficients induced by glucose in living tissues compared to slice tissues can differ by several times, even up to tenfold.
In obtaining accurate data through in vivo animal tests, Larin et al. utilized Optical Coherence Tomography (OCT) on live animals. They found that with a 1 mmol/L increase in blood glucose concentration and using 1300 nm light, the animal’s living tissue scattering coefficient decreased by 0.22% [
12]. However, due to the limited depth of OCT, these results were obtained using wavelengths where tissue absorption is weaker. To date, no studies have verified whether the scattering coefficient changes measured with these wavelengths are applicable to wavelengths with stronger absorption. Therefore, obtaining scattering coefficient variation data in the 1000–1700 nm band is crucial for fully understanding the impact of glucose fluctuations on scattering coefficients. These data will help determine the sensitivity within using a wavelength band, selecting the measuring wavelength to achieve the optimal possible glucose response.
This paper begins by defining the sensitivity of absorbance to glucose changes in in vivo measurements and analyzes the parameters required to estimate this absorbance sensitivity. In light of the challenges in acquiring data on scattering coefficient variations caused by glucose, we introduce a method that integrates MC simulations with in vivo experiments. This method is employed to estimate the alterations in scattering coefficients due to glucose concentration changes when using monochromatic light in in vivo measurements. Furthermore, by applying the Mie scattering theory [
17], we extend the observed changes in scattering coefficients with monochromatic light to the 1000–1700 nm wavelength range. Additionally, we conducted MC simulations for a three-layer skin model to acquire the absorbance sensitivity of glucose in this wavelength range. The human OGTTs are conducted with five typical wavelengths in the 1000–1700 nm range to validate the calculated sensitivity, while the method of extending the scattering properties of 1314 nm light to the 1000–1700 nm range is also validated as feasible.
We utilized the high-precision monitoring equipment and in vivo measurement approach described in Han’s paper, ensuring the accuracy and reliability of the human experimental data in our study. Building upon Han’s foundation, we conducted rigorous theoretical calculations to determine the sensitivity of blood glucose at any given wavelength and source–detector separation within the 1000–1700 nm range and then validated these calculations through extensive human experimentation. Our primary focus was to measure the change in the scattering coefficient caused by variations in glucose concentration in human skin, a detail not addressed in Han’s findings. Therefore, our work provides significant theoretical insights relevant to understanding the components of the glucose signal. It enables us to assess sensitivity regardless of whether the signal arises from glucose absorption or scattering, offering valuable data that are crucial for designing spectroscopic measurement instruments.
3. Calculated Results of Sensitivity
Glucose absorbance sensitivity, denoted as
, depends on four factors:
and
, and the resultant absorbance changes due to the absorption and scattering coefficient shifts. In this equation,
was established in
Section 2.4 Human Subject Experiments, while
can use the reported data [
27], as shown in
Figure 4.
and
are obtained through MC simulations. To calculate
, the scattering coefficient at each wavelength is kept constant. The MC simulation yields the absorbance when the absorption coefficient rises by 1 mm
−1, from which
is determined, as shown in
Figure 5a. Similarly, to calculate
, the absorption coefficient at each wavelength is kept constant. The absorbance when the scattering coefficient rises by 1 mm
−1 is calculated to obtain
, as shown in
Figure 5b.
Based on
and
depicted in
Figure 4 and
Figure 5a, we can calculate the absorbance sensitivity due to glucose’s absorption,
, as illustrated in
Figure 5c. Similarly, based on
and
presented in
Figure 3c and
Figure 5b, we can calculate the absorbance sensitivity due to glucose’s scattering,
, as demonstrated in
Figure 5d.
From
Figure 5c, it can be observed that
is relatively small in the 1000–1400 nm range, with values approaching 0. However, in the 1400–1700 nm range, it is significantly larger, varying between 0.0004 a.u. and 0.0012 a.u., and is dependent on the SDSs. Within the used SDS range (1.7–2.9 mm), the larger the SDS, the greater the
.
Figure 5d reveals that the
varies noticeably in the 1350–1600 nm range with SDS (1.7–2.9 mm), reaching its maximum absolute value at SDS = 2.9 mm. Additionally, for SDS > 2.3 mm in the wavelength range where tissue absorption is stronger (1350–1700 nm),
and
have opposite signs, implying that these two effects partially offset each other’s sensitivity. Overall, by comparing
Figure 5c,d, it is evident that the impact of the scattering effect,
, is greater than that of the absorption effect. Therefore, changes in tissue scattering due to glucose are the primary contributors to glucose sensitivity in the human body.
Figure 5e presents the sensitivity curves of the combined effects of glucose absorption and scattering on the human body at various SDSs across different wavelengths. It is observable that in the 1000–1350 nm range, there is little variation in sensitivity across wavelengths, with higher sensitivity achievable at the closer SDS of 1.7 mm. The sensitivity in this range is primarily determined by the scattering effect,
. In the 1350–1700 nm range, the sensitivity is jointly determined by the effects of
and
. At the closer SDS of 1.7 mm, higher ‘positive’ sensitivity can be obtained, while at the farther SDS of 2.9 mm, higher absolute ‘negative’ sensitivity can be achieved.
In the practical application of human measurement, multiple wavelengths and multivariate analysis can be employed, considering the influence of fluctuations in other tissue components. The selected wavelengths should be not only sensitive to changes in glucose but also to hemoglobin, water, albumin, urea, and lactic acid, enabling multivariate analysis to effectively extract the blood glucose signal. According to
Figure 5e, glucose sensitivity is high within the range of 1350–1700 nm. However, given the strong absorption of water near 1450 nm, the range of 1500 to 1700 nm is deemed more suitable for glucose sensing. The choice of measurement wavelength for other major blood components, such as albumin, hemoglobin, water, urea, and lactic acid, can refer to reference [
27,
28]. For instance, Ren et al. examined the absorption spectra of glucose, lactic acid, and urea in the near-infrared band, while Deng et al. obtained the absorption spectra of water. Similarly, the strong absorption peak of water around 1450 nm should be circumvented when sensing these components.
4. Results of Human Validation Experiment
To validate the estimated of the 1000–1700 nm range and the simulated absorbance sensitivity , 25 human OGTTs were conducted. The used wavelengths are 1050 nm, 1219 nm, 1314 nm, 1380 nm, 1550 nm, and 1609 nm. The 25 human OGTT cases demonstrated a good correlation between and and were considered reliable experimental results for the validation of and sensitivity.
The differential absorbance
at the SDSs of 1.7 mm and 2.0 mm was recorded;
Figure 6a–d shows the human experiment results of 4 volunteers among 25 human OGTT cases. The sensitivity results for all 25 subjects are shown in
Figure 6e. The average sensitivity
of 25 cases is shown in
Figure 6f and compared with the MC results.
As shown in
Figure 6e, for the sensitivity results of each of the 25 subjects, the 1550 nm light exhibited the maximal sensitivity among the six wavelengths, with an average sensitivity of −0.00107 a.u.·(mmol/L)⁻
1. As shown in
Figure 6f, it can be observed that at the six wavelengths, the calculated and experimental sensitivity of differential absorbance
AD are close to each other. The relative errors between the average value of experimental sensitivity and the calculated sensitivity at 1050 nm, 1219 nm, 1314 nm, 1380 nm, 1550 nm, and 1609 nm are, respectively, 13%, 26%, 0%, 21%, 3%, and 9%. This indicates that the glucose signal sensitivity calculated in this study shows good consistency with the sensitivity acquired in the human experiment. Moreover, the sensitivity measured at 1550 nm closely matches the calculated sensitivity, indicating that the
obtained at 1314 nm can be extended to wavelengths where human tissue absorption is stronger.
The measured sensitivity in the trial at 1219 nm shows a greater difference from the calculated sensitivity, primarily due to the weaker response of glucose changes at this wavelength, which is more susceptible to the influence of fatty substances in tissue. The 1380 nm light has strong absorption due to water, and in measurements, it is easily affected by changes in body tissue water content, which is a major reason for the discrepancy between the actual values and the calculated sensitivity.
5. Discussion
In the human experiment, it was found that using light at 1314 nm, the relative change in the scattering coefficient due to 1 mmol/L blood glucose concentration change was −0.29%. Based on these data and the Mie scattering theory, we assume that the relative change in the scattering coefficient due to glucose variation is −0.29% for all lights within the 1000–1700 nm range, which are used for the calculation of the sensitivity of the absorbance. The human OGTTs were conducted to validate the calculated sensitivity results of glucose signals.
In our measurements on 25 subjects using 1314 nm, the relative change rate of scattering coefficient with blood glucose concentration was estimated to be −0.29% ± 0.098%. The variation among subjects is likely attributable to differences in skin structure, tissue components, and the physiological state of the skin.
The
obtained in this paper is compared with the published results, as shown in
Table 1.
It can be seen from
Table 1 that the
measured in vivo in the reported studies is much greater than that in the phantom (such as intralipid), which is about 11–306 times that in 2% intralipid. The result of this paper is a factor of 14.5, which is also in this range and is on the small side.
The results measured in vivo in
Table 1 were mostly obtained by the OCT method, and there was a significant difference between their values. The reason for this significant difference may be due to the different measurement objects, measurement depths, and wavelengths used in various studies. In addition, the research of GAO et al. provided an explanation, that is, they found that the scattering coefficient change with glucose in blood and interstitial fluid is opposite [
29]. The
in blood vessels with a diameter of more than 20 µm is +2.7%, while the
is −0.75% in the interstitial fluid of the dermis. Therefore, different blood and interstitial fluid distributions in tissues will lead to different measurement results. In
Table 1, the other studies using OCT measurements, except GAO, only give the results of
in a specific depth range without distinguishing between blood and interstitial fluid. Therefore, the different tissue objects with different probing depths can significantly affect the magnitude of
.
In this paper, we used near-infrared diffuse reflectance spectroscopy. The measurement principle is different from the OCT method. It does not locate at a specific depth in detection. The measured diffuse light passes through all layers of skin tissue, including the epidermis, dermis, and a part of subcutaneous tissue. The average tissue scattering change caused by glucose in the whole light path will affect the results, so there is a certain difference with the OCT results.
5.1. Comparison of Calculated Sensitivity with Experimental Sensitivity
The MC method demonstrates the sensitivity of absorbance caused by glucose for the wavelengths within the 1000–1700 nm range. Then, the MC results were well validated by the human experiment at six typical wavelengths within this range. The experimental results of the differential absorbance of SDSs of 1.7 mm and 2.0 mm show good consistency with the theoretical sensitivity calculated by the MC method. Additionally, we observed that among these six wavelengths, the sensitivity of differential absorbance is slightly higher for the three longer wavelengths (1380 nm, 1550 nm, and 1609 nm) and slightly lower for the three shorter wavelengths (1050 nm, 1219 nm, and 1314 nm). Therefore, for the differential measurement method, longer wavelengths should be preferred while avoiding wavelengths with strong water absorption, such as the 1500–1600 nm range.
5.2. Reasonably Choosing SDSs Based on Sensitivity Results
Our method can determine the absorbance sensitivity caused by glucose,
, at different SDSs for the wavelengths within the 1000–1700 nm range, serving as a reference for selecting reasonable SDS. Based on the calculated result, we can identify the “floating reference point” on the human body [
30], i.e., the SDS where absorbance sensitivity is close to zero, as shown in
Figure 5e. When using light within the 1000–1300 nm range, the sensitivity is close to zero at SDS = 2.6 mm. The detection signal at this SDS only varies with the human body background and can serve as a reference for eliminating physiological background interference. When using light within the 1300–1560 nm range, the “floating reference point” is about 1.7–2.0 mm; for light within the 1560–1700 nm range, it is within the 2.6–2.9 mm range.
The choice of SDS according to the calculated sensitivity also dictates the requirements for the instrument noise level. To achieve higher sensitivity levels, a larger SDS should be chosen, but the larger the SDS, the weaker the detected light energy. In this case, it is necessary to minimize instrument noise. For instance, when using light within the 1500–1560 nm range, the sensitivity at SDS = 2.9 mm is −0.00083 a.u. To discern a 1 mmol/L glucose concentration change, the random noise of the absorbance at this SDS needs to be controlled to 0.00061 a.u. or less. If the instrument’s performance cannot meet this requirement, then the measurements should use a smaller SDS to ensure accuracy.
5.3. Comparison of Glucose Sensitivity in Human Body and Intralipid Solution
The methodology presented in this paper can also be applied to estimate the sensitivity of glucose measurement in intralipid phantoms.
Figure 7a shows the absorbance sensitivity for detecting glucose in intralipid solutions, while
Figure 7b shows the absorbance sensitivity due to glucose’s absorption change, and
Figure 7c shows that due to glucose’s scattering change. Comparing
Figure 5e and
Figure 7a, it is evident that the glucose detection sensitivity in the human body is greater than that in intralipid. The primary reason is that in the human body, the change in the scattering coefficient caused by glucose is greater than in intralipid, as shown in
Figure 3c. Moreover, by comparing
Figure 5c,d, we can understand that the glucose detection sensitivity in the human body is mainly contributed by a glucose-caused scattering change. However, the situation is different in the intralipid solutions; as observed in
Figure 7b,c, the effects of absorption and scattering changes on sensitivity in the medium are similar in magnitude but opposite in direction, which counterbalance each other. Consequently, the glucose sensitivity in the intralipid solution is lower than in the human body. In short, despite the absorption changes induced by glucose contributing similarly to sensitivity in both intralipid and human tissues, the scattering variations caused by glucose are markedly different in these two media. This distinction is the primary factor leading to the differences in glucose sensitivity observed between intralipid and human tissues.
5.4. Preliminary Research on Diversity in Sensitivity among Individuals
Through experimentation, we learned that glucose detection sensitivity is influenced by the optical properties of human tissue. Differences in skin’s optical properties among individuals directly determine the level of sensitivity. In our study, we measured the sensitivity of differential absorbance for the forearm skin in 25 subjects and found their difference due to skin characteristics. It needs to be clarified that the discussion on the diversity in sensitivity among individuals is preliminary, as it is only for the tests on the forearm (the extensor side) of these Asian subjects.
Skin’s optical properties vary with age and gender. In the experiment, we divided the 25 cases into four groups based on the subject gender and age: elderly male group (6 individuals), elderly female group (7 individuals), young male group (6 individuals), and young female group (6 individuals), and we calculated the average sensitivity for each group, as shown in
Figure 8.
The results indicate that at wavelengths of 1550 nm and 1609 nm, the amplitude of glucose detection sensitivity is larger with less measurement interference, so particular attention should be paid to results at these wavelengths. In
Figure 8a, we can see that the glucose sensitivity in females is slightly higher than in males. However,
Figure 8b shows no significant difference in sensitivity between elderly subjects (over 65 years) and younger subjects (under 35 years).
According to research by Kono et al. [
31], gender does not significantly affect the skin’s absorption coefficient, but there is a notable gender difference in the scattering coefficient when using light in the 1000–1600 nm range, with the average scattering coefficient in female skin being about 18.8% higher than in male skin. Our study found that due to the higher scattering coefficient, their
is larger, and the change in the scattering coefficient has a greater impact on the sensitivity of glucose signals. This might explain the slightly higher glucose detection sensitivity in females.
We also noted that skin differences among people of different ages are mainly manifested in water content. According to research by Nakagawa et al. [
32], the water content in the dermis of the elderly (73.3 ± 1.6%) is higher, while in younger people, it is slightly lower at 69.3 ± 1.7%. Although elderly people might have a higher absorption coefficient due to increased skin water content, our experiments did not observe a significant impact of age-related differences in the absorption coefficient on the sensitivity, as shown in
Figure 8b. However, elderly females, compared to the other three groups, have higher absorption and scattering coefficients. These combined effects could lead to slightly higher glucose signals in this group, as shown in
Figure 8c.