Next Article in Journal
TrackInk: An IoT-Enabled Real-Time Object Tracking System in Space
Next Article in Special Issue
Smart Bio-Impedance-Based Sensor for Guiding Standard Needle Insertion
Previous Article in Journal
Hash-Chain Fog/Edge: A Mode-Based Hash-Chain for Secured Mutual Authentication Protocol Using Zero-Knowledge Proofs in Fog/Edge
Previous Article in Special Issue
A Way of Bionic Control Based on EI, EMG, and FMG Signals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Patient Specific Numerical Modeling for Renal Blood Monitoring Using Electrical Bio-Impedance

Department of Medical and Technical Information Technology, Bauman Moscow State Technical University, 105005 Moscow, Russia
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(2), 606; https://doi.org/10.3390/s22020606
Submission received: 29 November 2021 / Revised: 30 December 2021 / Accepted: 10 January 2022 / Published: 13 January 2022
(This article belongs to the Special Issue Bioimpedance Sensors: Instrumentation, Models, and Applications)

Abstract

:
Knowledge of renal blood circulation is considered as an important physiological value, particularly for fast detection of acute allograft rejection as well as the management of critically ill patients with acute renal failure. The electrical impedance signal obtained from kidney with an appropriate electrode system and optimal electrode system position regarding to the kidney projection on skin surface reflects the nature of renal blood circulation and tone of renal blood vessels. This paper proposes a specific numerical modelling based on prior information from MRI-data. The numerical modelling was conducted for electrical impedance change estimation due to renal blood distribution. The proposed model takes into the account the geometrical and electrophysiological parameters of tissues around the kidney as well as the actual blood distribution within the kidney. The numerical modelling had shown that it is possible to register the electrical impedance signal caused by renal blood circulation with an electrode system commensurate with the size of kidney, which makes it possible to reduce the influence of surrounding tissues and organs. Experimental studies were obtained to prove the numerical modelling and the effectiveness of developed electrode systems based on the obtained simulation results. The obtained electrical impedance signal with the appropriate electrode system shows very good agreement with the renal blood change estimated using Doppler ultrasound. For the measured electrical impedance signal, it is possible to obtain the amplitude-time parameters, which reflect the hemodynamic characteristics of the kidneys and used in diagnostics, which is the subject of further research.

1. Introduction

Renal function is intimately dependent on renal blood circulation [1]. However, the measurement of renal blood circulation is time consuming and not widely available, thus the current methods are either invasive, rely on the use of radioactive tracers which are diffusible, have kinetics which are complicated by tubular transport or glomerular filtration, or cannot allow the assessment of small variations and rapid changes in renal blood flow and cannot be used for bedside monitoring [2]. The real time measurement of renal blood flow using contrast-enhanced ultrasound was considered in [2,3], although, this method is cost-effective and widely available, it still has considerable limitations in sensitivity and specificity for the diagnosis of allograft rejection [3]. Renal near-infrared spectroscopy is used to evaluate regional oximetry in a non-invasive continuous real-time fashion, and reflects tissue perfusion [4], however the near-infrared method has a depth limitation, which makes this approach applicable only to special patients like newborns.
Several studies on animals and subjects was conducted to prove that the electrical impedance signal obtained from kidney correlate with the change of blood circulation and reflect the tone of blood vessels [5,6]. Nevertheless, due to the anatomical location of kidney especially with obese patients as well as the adjacent organs such as liver, spleen even the surrounds muscles and fat tissues made the authors, who studied this method previously to conclude that the electrical impedance signal obtained using skin attached electrodes gives a characteristic of the complex blood flow change of all structures located in the inter-electrode space. For precise assessment, implantable electrical impedance device, which can provide measurements on wide range of frequencies was considered to diagnose the rejection in kidney transplantation in previous studies [5,7].
However, in order to achieve high accuracy of determining the amount of blood per 100 G of kidney tissues using attached non-invasive electrodes, it is required to develop a detailed specific mathematical model of the region of interest as well as a method to eliminate the effect of surrounding organs and tissues. Improvement of this technique can allow its application for monitoring patients undergo different type of cardiovascular surgeries, who need a monitoring to prevent the acute kidney injury, which is a common complication after cardiac surgery and is associated with increased morbidity and mortality and increased length of stay in the intensive care unit [8]. Furthermore, the develop of such system with precise measurements can be used for acute kidney allograft rejection diagnosis. In related previous work [9] belongs to the authors of this paper, different biophysical models were considered to describe the kidney and take in to the account the electro physiological parameters, and the depth of kidney. Nevertheless, these models are more reliable for non-obese patients and children because as the kidneys are located deeper the impact of the surrounding tissues become greater and hence leads to an error of determining the amount of blood per 100 G of kidney tissues. However, the desirable area to acquire the electrical impedance signal can be presented as a multi-layered model; these layers specify the biological tissues of the region of interest, which are: the skin, subcutaneous fat layer, the soft muscle layer, the fat layer around the kidney, and the kidney layers. According to the analysis of the electrical properties of biological tissues at several excitations frequencies that has been considered in related works [9,10], it was concluded that the most adequate frequency range is 80–100 kHz where active resistance is the main component of impedance. Thus, at lower frequencies the impact of capacitive part increases, which is determined by the heterogeneity of tissue structures in the region of interest, furthermore the low frequencies can cause seizures and polarizing electronic phenomena. Thus, an electrical current with 100 kHz has been used during this study. Figure 1 shows the electrical resistivity of biological tissues in the study area in different frequencies including 100 kHz [11]. Four electrode system configuration was proposed in this study to provide the desired accuracy for bio-impedance measurement, thus two outer electrodes are used to inject a small alternating current into the human body the voltage drop is measured to calculate the complex impedance [12].
However, the aim of this study is to conduct a numerical modeling based on prior information obtained from MRI-images. This modelling should take into the account the geometrical and electrophysiological parameters of the tissues above the kidney as well as the actual kidney dimensions and blood distribution in order to estimate the contribution of the surrounding tissues around the kidney located in the inter-electrode space. Additionally, the task of verifying the developed model with the help of an experimental study was set.

2. Materials and Methods

Numerical modeling has been conducted using a simulation platform COMSOL Multiphysics 5 [13,14,15]. MRI images have been used as an input for 3-D model creation. The preliminary processing was carried out using MicroDicom [16]. In order to reduce the calculation time, the right kidney was considered during the numerical modelling. During this study, the actual kidney size, the dimensions of the patient’s torso, the thickness of subcutaneous fat and muscle layers, as well as the distance of the kidney relative to the center of the trunk cross-section were outlined as shown in Figure 2.
The 3-D model was implemented in COMSOL Multiphysics, while the electrical resistivity of tissues in the region of interest was set according to Gabriel et al. [11]. Figure 3 shows the obtained 3-D model of perineal space with an electrode system located at the center of right kidney.
In this model, the kidney is represented by a two-layer ellipse. These layers mimic the renal cortex and medulla. The kidney is located in the fibrous layer, presented in the form of a cylinder with major and minor axes (2a1 − 2hm) and (2b1 − 2hm), respectively. The electrodes are modeled as cylinders equidistant from each other with a base diameter of 3 mm. Table 1 shows the geometrical dimensions of the designed model obtained from MRI-images.
The change of kidney electrical resistivity due to the blood filling during cardiac cycle has been estimated using the parallel conductor model, illustrated in Figure 4 [17].
As the kidneys receive a total blood flow of approximately 1000 mL per minute (20% of the cardiac output). This equates to 300–400 mL per minute per 100 g of tissue [18]. The Equation (1), which deduced from the parallel conductor model, has been used to calculate the electrical resistivity change of healthy kidney which receives a total blood flow of approximately 400 mL per minute per 100 g of tissue.
Δ ρ = ρ k     ρ total = ( ρ k     mb   ( ρ k     ρ b ) ) / ( mk     ρ b   +   mb     ρ k )
where ρk is the electrical resistivity of kidney tissue, ρb is the electrical resistivity of blood; ρtotal is the electrical resistivity of a blood-filled kidney.
In previous work [10] it has been concluded that the estimation using the parallel conductor model showed that the dynamic electrical resistivity change of the kidney due to nature of renal blood circulation in normal cases is approximately 1 Ω·m and the cortical as well as the medulla layers perceive such a change in different percentages according to the data illustrated in Table 2. The table shows the static electrical resistivity of kidney layers before and after blood filling, two different cases were considered in this study to estimate the sensitivity of electrical impedance measurement to different blood flow distribution in kidney layers. The first case: the cortex receives a 90% of the blood flow to the kidney which corresponds to normal perfusion while the second case suggests less perfusion of cortex.
During the numerical simulation, in line four electrodes with equal distance between the adjacent electrodes was used, thus the distance between the measuring electrodes (b) was varied in the range from 7 to 40 mm with a step of 3 mm. Accordingly, the distance between the current electrodes, which is 3 × b has been changed. The static electrical impedance values as well as the bio-impedance change have been calculated for every electrode system iteration change. The different cases of layered distribution of renal blood flow were considered during the modeling.
The experimental studies on subject were carried out in the medical and technological center of Bauman Moscow state technical university. The multichannel system (REO-32) has been used for electrical impedance measurements. The system provides 4-electrode technique by applying alternating current of frequency 100 kHz and constant amplitude to a current electrodes 3 mA [19]. The technical specifications of REO32-system are shown in Table 3.
The experimental studies were carried out on the same subject, whose MRI images were used for numerical modeling. In the study, 3 different electrode systems were designed as shown in Figure 5. The distances between the current (2a) and measuring (2b) electrodes for each electrode system are illustrated in Table 4.
Figure 6 shows the electrode system positioning during the experiment. The studies were carried out in the area of the right kidney and electrical impedance mapping were performed with 20 mm shifts to the right due to the design dimensions of the electrode system used. In order to reduce the influence of kidney displacement during breathing, the measurements were taken with breath holding. The measurements were carried out several times to confirm the reproducibility of the obtained signals; uniform conditions as well as the degree of pressing of the electrode system were provided.

3. Results

3.1. Simulation and Experimental Results

The presented plots in Figure 7 show the results of numerical modeling of two cases of layered renal blood distribution. In Figure 7a we show the static electrical impedance change in relation to the electrode spacing change, thus the impedance change is a function which depends on the apparent resistivity of tissues in the reign of interest as well as the geometrical factor of the electrode system. In the case of using a symmetric tetrapolar measurement system, the impedance between the measuring electrodes can be represented with Equation (2) [20]:
Z   = 2 ρ b π ( a 2 b 2 )
The significant contribution on static electrical bio-impedance in the area of interest is the fat layer as it has the largest electrical resistivity. The involved subject in this study is 35 year man and weight 55 kg, thus as the electrode spacing increases the impedance decreases for every iteration change which mean that more conductive layers are captured. The corresponding plots of the dynamic electrical impedance change due to blood filling is shown in Figure 7b, during this study it was suggested the origin of the recorded electrical impedance change is the renal blood circulation, thus the influence of the layers above the kidney was neglected during the numerical modelling. In contract with static impedance, the dynamic electrical impedance increases due the increase of the electrode spacing, this indicates the increase of sensitivity to kidney blood filling.
The electrical impedance change for the considered two cases in this study is illustrated in Table 5 with details. The presented values show a high sensitivity of electrode systems with distances of 20 mm and more between the measuring electrodes to the changes in electrical impedance of kidney due to its blood filling, since the probing depth also increases with an increase in the inter-electrode distance [21]. The results analysis shows an impedance value in a range of 140 mΩ can be obtained, which enough for further inverse solution to estimate the hemodynamic parameters of kidney.
The big issue of renal blood flow estimation is the kidney displacement during respiration, thus this displacement can reach few centimeter depending on the breathing depth [22]. The kidney displacement in one direction was considered during this study and modelled by the increase of kidney depth, which can be occurred during inhale. Figure 7c shows the change of electrical impedance due to different kidney depth, hence for real time monitoring this influence should be considered. However, in normal cases when the kidney is located far from the body surface the small electrode system cannot provide the desire sensitivity to kidney blood circulation as shown in Figure 7d, thus the desired impedance signal can be achieved just with increasing the electrode spacing, this increase leads to an increase of the surrounding tissues, however the multi-channel measurement with different electrode spacing can be a solution to increase the accuracy of determining the amount of blood per 100 G of kidney tissues.
The result of electrical bio-impedance recording is shown in Figure 8, before measurement the location of kidney was identified by ultrasound system, the skin was preliminarily treated with isopropyl alcohol (2-propanol, isopropanolm, IPA) 70% (w/v) to remove the stratum corneum and increase the conductivity. The same electrode system location was achieved as proposed on numerical modeling, the subject was asked to hold breathing for 10 s which is enough to obtain a series of cardiac cycles to make sure of signal reliability. The first measurement was obtained using the smallest electrode system; the static electrical impedance as well as the dynamic was acquired. The same measurement was repeated with the second and the third electrode systems. As shown from the obtained data, the signal becomes clear and more specific to renal blood circulation with appropriate electrode spacing selection, which confirms the requirement of patient specific modelling before measurements to select the appropriate electrode system as well as the right position for every patient individually. The results of experimental studies on subject and illustrated in Table 6 shows the same behavior of numerical modelling, thus the amplitude of pulsatile impedance becomes greater and specifically clear with electrode spacing changes, however the observed difference between the experimental and the modelling values is because of the electrical resistivity of biological tissues used during the numerical simulation differ from the real values as the first were taken from data obtained in vitro. However, this difference should be adjusted during the further inverse problem solution.

3.2. Verification with Doppler Ultrasound

To assess the contribution of pulse changes in the blood filling of the kidney vessels to the electrical impedance signals, Doppler ultrasound of the kidney vessels was performed. The General Electric models LOGIC S8 with an ML6-15-D ultrasonic transducer at a frequency of 6.3 MHz. was used during this study to obtain prior information such as renal artery and renal vein velocity as well as the corresponding diameters. Breath-holding was employed during ultrasound acquisitions. Figure 9 shows the recoded curves of blood velocities (cm/s) in the renal artery Va and renal vein Vv, the corresponding cross sectional area of the renal artery and renal vein are Sa = 0.30 cm2 Sv = 0.38 cm2 respectively. However, the kidney blood volume change per one cardio cycle Q(t) (cm3) can be calculated based on the difference in the volumetric inflow and outflow according to the Equation (3).
Q ( t ) = 0 t V a ( t ) S a dt 0 t V v ( t ) S v dt ,
Using the obtained curves, it is easy to achieve to the dependences of the volumetric flow rate. Thus, the values of the blood flow velocities at each time point were multiplied by the cross sectional area of the renal artery and renal vein to get the corresponding volumetric flow rate for renal artery and renal vein as shown Figure 10. The obtained volumetric blood flow curves were used for the integral calculation of the change in blood volume over one cardio-cycle.
Figure 11 show the result of the bio-impedance pulse behavior comparing with the blood volume change within the kidney obtained using ultrasound.
The correlation coefficient between the change in kidney blood volume, calculated from the results of ultrasound measurements, and the change in the electrical impedance during one cardio cycle is 0.89. This results show that the electrical impedance signal obtained with attached electrode reflects the nature of blood circulation in the kidney and showed high correlation with the renal blood behavior of kidney.

4. Discussion

Electrical Impedance is a simple and non-invasive method to acquire data on hemodynic parameter of kidney. However, the accuracy of determining this parameter can be affected by various dynamic sources. The origin of these sources can be the adjacent organs and tissues around the kidney as well as the kidney displacement during respiration. Thus, the analysis of the geometrical and electrophysiological parameters of the kidney and tissues around shows that for the study of signal mechanisms formation and the transition from measurable electrical quantities to hemodynamic parameters it is required to develop a justified mathematical model which characterizes the region of interest with full details. The simple analytical models do not fully describe the area of study, which leads to a significant error when solving the inverse function, thus the proposed model in this study takes into account the geometric and electrophysiological parameters of the tissues located around the kidney, as well as the actual distribution of blood between the layers of the kidney. The goal of this study was to develop a patient specific model of perirenal space based on MRI data, which can be used to investigate the mechanisms of electrical bio-impedance signal obtained from kidney using skin attached electrodes. Numerical modeling has shown that it is possible to record the electrical impedance signal caused by renal blood circulation using a specially selected electrode system, commensurate with the size of the kidney and located relative to the longitudinal axis of the kidney, which makes it possible to reduce the influence of surrounding tissues and organs. The kidneys can move during respiration in a tilted coronal and sagittal plane [22], however in order to reduce the influence of kidney replacement during breathing, the measurements were taken with breath holding.
The experimental study on subject with different electrode system shows that electrical impedance signal quality can be influenced by the electrode spacing and some specific points can be disappear in case of measurement with no appropriate electrode system, thus previously identified electrode spacing based on the forward problem solution obtained by the numerical problem can avoid this problem. An important point in this work is the correlation between the registered electrical impedance signals with the arising pulse blood flow fluctuations from kidney obtained using Doppler ultrasound.

5. Conclusions

The paper proposed a new approach which allows the monitoring of renal blood circulation based on electrical impedance measurement from specific location regarding to the kidney. Specific numerical modeling has been conducted in order to choose the appropriate electrode system parameters, thus the electrical impedance signal which reflects the nature of renal blood circulation for every patient depends on such parameters. Taking in to the consideration the impact of the surrounding tissues by developing such specific model allows increasing the accuracy of measurement and the hemodynamic parameters, which is the subject of further study.

Author Contributions

Methodology, M.A.-h., E.C. and S.S.; data organization, M.A.-h. and E.C.; investigation, M.A.-h. and E.C.; data visualization, M.A.-h. and E.C.; writing—review and editing, M.A.-h., E.C. and S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This research was approved by the medical center of Bauman Moscow State technical University and approved by the University of Bauman Institutional Review Board.

Informed Consent Statement

Informed consent was obtained from the subject involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Regan, M.C.; Young, L.S.; Geraghty, J.; Fitzpatrick, J.M. Regional renal blood flow in normal and disease states. Urol. Res. 1995, 23, 1–10. [Google Scholar] [CrossRef] [PubMed]
  2. Kalantarinia, K.; Belcik, J.T.; Patrie, J.T.; Wei, K. Real-time measurement of renal blood flow in healthy subjects using contrast-enhanced ultrasound. Am. J. Physiol. Renal Physiol. 2009, 297, 1129–1134. [Google Scholar] [CrossRef] [PubMed]
  3. Thölking, G.; Schuette-Nuetgen, K.; Kentrup, D.; Pawelski, H.; Reuter, S. Imaging-based diagnosis of acute renal allograft rejection. World J. Transplant. 2016, 6, 174–182. [Google Scholar] [CrossRef] [Green Version]
  4. Ruf, B.; Bonelli, V.; Balling, G.; Hörer, J.; Nagdyman, N.; Braun, S.L.; Ewert, P.; Reiter, K. Intraoperative renal near-infrared spectroscopy indicates developing acute kidney injury in infants undergoing cardiac surgery with cardiopulmonary bypass: A case–control study. Crit. Care 2015, 19, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Bogonez, P.; Riu, P.J. Implantable bioimpedance system for measuring the impedance of kidney. In Proceedings of the 13th International Conference on Electrical Bioimpedance and the 8th Conference on Electrical Impedance Tomography, Graz, Austria, 29 August–2 September 2007; Springer: Graz, Austria, 2007; pp. 256–259. [Google Scholar]
  6. Hanssen, O.; Erpicium, P.; Lovinfosse, P.; Meunier, P.; Weekers, L.; Tshibanda, L.; Krzesinski, J.M.; Hustinx, R.; Jouret, F. Non-invasive approaches in the diagnosis of acute rejection in kidney transplant recipients. Part I. In vivo imaging methods. Clin. Kidney J. 2017, 10, 97–105. [Google Scholar] [PubMed] [Green Version]
  7. Rodriguez, S.; Ollmar, S.; Waqar, M.; Rusu, A. A Batteryless Sensor ASIC for Implantable Bio-Impedance Applications. IEEE Trans. Biomed. Circuits Syst. 2016, 10, 533–544. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Lopez-Gomez, J.M. Evolution and applications of bioimpedance in managing chronic kidney disease. Nefrologia 2011, 31, 630–634. [Google Scholar] [PubMed]
  9. Al-harosh, M.; Yangirov, M.; Kolesnikov, D.; Shchukin, S. Bio-Impedance Sensor for Real-Time Artery Diameter Waveform Assessment. Sensors 2021, 21, 8438. [Google Scholar] [CrossRef] [PubMed]
  10. Al-Harosh, M.B.; Chernikov, E.S.; Shchukin, S.I.; Gries, T.; Leonhardt, S. Renal Blood Monitoring System Using Bio-impedance Measurement: Pilot Study. In Proceedings of the 2021 IEEE Ural Symposium on Biomedical Engineering, Radioelectronics and Information Technology (USBEREIT), Yekaterinburg, Russia, 13–14 May 2020; pp. 199–202. [Google Scholar]
  11. Gabriel, C.; Gabriel, S.; Corthout, E. The dielectric properties of biological tissues: I. Literature survey. Phys. Med. Biol. 1996, 41, 2231–2249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Al-harosh, M.B.; Shchukin, S. Peripheral vein detection using electrical impedance method. J. Electr. Bioimpedance 2017, 8, 79–83. [Google Scholar] [CrossRef] [Green Version]
  13. Dokos, S. Modelling Organs, Tissues, Cells and Devices: Using Matlab and Comsol Multiphysics; Springer: Berlin/Heidelberg, Germany, 2017; pp. 355–379. [Google Scholar]
  14. Tikhomirov, A.N.; Briko, A.N.; Seleznev, N.V.; Shchukin, S.I.; Levando, A.M.; Murashko, M.A. Development of a Geometric Model of the Heart and Chest for Multichannel Electrical Impedance Computer Cardiography Technology. In Proceedings of the 2020 IEEE Ural Symposium on Biomedical Engineering, Radioelectronics and Information Technology (USBEREIT), Yekaterinburg, Russia, 13–14 May 2020; pp. 32–35. [Google Scholar]
  15. Tikhomirov, A.N.; Shchukin, S.I.; Leonhardt, S.; Volkov, A.K.; Murashko, M.A. Multichannel Electrical Impedance Methods for Monitoring Cardiac Activity Indicators. Biomed. Eng. 2019, 52, 365–370. [Google Scholar] [CrossRef]
  16. MicroDicom—Free DICOM Viewer and Software. Available online: http://www.microdicom.com/ (accessed on 28 October 2021).
  17. Tagawa, T.; Tamura, T.; Oberg, P.A. Biomedical Sensors and Instruments; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
  18. Schmidt-Nielsen, B.; Schmidt, D. Renal function of Sphenodon punctatum. Comp. Biochem. Physiol. Part A Physiol. 1973, 44, 121–129. [Google Scholar] [CrossRef]
  19. Al-Harosh, M.B.; Shchukin, S.I. The Venous Occlusion Effect to Increase the Accuracy of Electrical Impedance Peripheral Veins Detection. In Proceedings of the 2017 Joint Conference of the European Medical and Biological Engineering Conference (EMBEC) and the Nordic-Baltic Conference on Biomedical Engineering and Medical Physics (NBC), Tampere, Finland, 11–15 June 2017; Springer: Berlin/Heidelberg, Germany, 2017; pp. 538–541. [Google Scholar]
  20. Al-Harosh, M.B.; Shchukin, S.I. Numerical modeling of the electrical impedance method of peripheral veins localization. In Proceedings of the World Congress on Medical Physics and Biomedical Engineering, Toronto, ON, Canada, 7–12 June 2015; Springer: Cham, Switzerland, 2015; pp. 1683–1686. [Google Scholar]
  21. Goidina, T.; Kobelev, A.; Gulyaev, Y. Precision Electrode System for Electrical Impedance Myography. In Proceedings of the 2021 IEEE Ural Symposium on Biomedical Engineering, Radioelectronics and Information Technology (USBEREIT), Yekaterinburg, Russia, 13–14 May 2020; pp. 143–146. [Google Scholar]
  22. Siva, S.; Pham, D.; Gill, S.; Bressel, M.; Dang, K.; Devereux, T. An analysis of respiratory induced kidney motion on four-dimensional computed tomography and its implications for stereotactic kidney radiotherapy. Radiat. Oncol. 2013, 8, 1–8. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The electrical resistivity of the biological tissues in the study area.
Figure 1. The electrical resistivity of the biological tissues in the study area.
Sensors 22 00606 g001
Figure 2. Assessment of the kidney location, its geometric dimensions and the thickness of the perirenal layers according to the MRI image.
Figure 2. Assessment of the kidney location, its geometric dimensions and the thickness of the perirenal layers according to the MRI image.
Sensors 22 00606 g002
Figure 3. (a) Model of the perirenal space in the COMSOL Multiphysics 5.5 software; (b) Geometric parameters of the model (cross section), where 2a1—is the major axis of the ellipse (cross section of the trunk), 2b1—is the minor axis of the ellipse (cross section of the trunk), hsubc—is the thickness of the subcutaneous fat layer, hm—thickness of the muscle layer, hc—thickness of the renal cortex.
Figure 3. (a) Model of the perirenal space in the COMSOL Multiphysics 5.5 software; (b) Geometric parameters of the model (cross section), where 2a1—is the major axis of the ellipse (cross section of the trunk), 2b1—is the minor axis of the ellipse (cross section of the trunk), hsubc—is the thickness of the subcutaneous fat layer, hm—thickness of the muscle layer, hc—thickness of the renal cortex.
Sensors 22 00606 g003
Figure 4. The parallel conductor model and the equivalent circuit, Where: V0, Vb The volume of kidney and blood respectively, ρ0, ρb—the electrical resistivity of kidney tissue and blood respectively.
Figure 4. The parallel conductor model and the equivalent circuit, Where: V0, Vb The volume of kidney and blood respectively, ρ0, ρb—the electrical resistivity of kidney tissue and blood respectively.
Sensors 22 00606 g004
Figure 5. The designed electrode systems: (a) the first, (b) the second, (c) the third.
Figure 5. The designed electrode systems: (a) the first, (b) the second, (c) the third.
Sensors 22 00606 g005
Figure 6. The electrode system position corresponding to the projection of kidney on skin surface.
Figure 6. The electrode system position corresponding to the projection of kidney on skin surface.
Sensors 22 00606 g006
Figure 7. The change of electrical impedance obtained by numerical modelling: (a) the variation of static impedance with increasing distance between measuring electrodes, (b) the variation of pulsatile impedance with increasing distance between measuring electrodes, (c) the variation of pulsatile impedance with increasing kidney depth, (d) the variation of pulsatile impedance with increasing distance between measuring electrodes at the maximum kidney depth.
Figure 7. The change of electrical impedance obtained by numerical modelling: (a) the variation of static impedance with increasing distance between measuring electrodes, (b) the variation of pulsatile impedance with increasing distance between measuring electrodes, (c) the variation of pulsatile impedance with increasing kidney depth, (d) the variation of pulsatile impedance with increasing distance between measuring electrodes at the maximum kidney depth.
Sensors 22 00606 g007
Figure 8. The electrical impedance signals obtained by (a) the first electrode system, (b) the second electrode system, (c) the third electrode system.
Figure 8. The electrical impedance signals obtained by (a) the first electrode system, (b) the second electrode system, (c) the third electrode system.
Sensors 22 00606 g008
Figure 9. This velocity curves obtained using ultrasound: (a) The renal artery blood velocity; (b) The renal vein blood velocity; (c) The averaged renal artery blood velocity; (d) The averaged renal vein blood velocity.
Figure 9. This velocity curves obtained using ultrasound: (a) The renal artery blood velocity; (b) The renal vein blood velocity; (c) The averaged renal artery blood velocity; (d) The averaged renal vein blood velocity.
Sensors 22 00606 g009
Figure 10. This volumetric flow rate in: (a) The renal artery; (b) The renal vein.
Figure 10. This volumetric flow rate in: (a) The renal artery; (b) The renal vein.
Sensors 22 00606 g010
Figure 11. The correlation between bio-impedance signal and the blood volume change of kidney using ultrasound.
Figure 11. The correlation between bio-impedance signal and the blood volume change of kidney using ultrasound.
Sensors 22 00606 g011
Table 1. The geometrical parameters of the designed model.
Table 1. The geometrical parameters of the designed model.
Model ParametersValue [mm]
2a1253.7
2b1166.71
hsubc5.54
hm10.74
hc5
Table 2. The electrical resistivity of kidney tissues a according to blood distribution.
Table 2. The electrical resistivity of kidney tissues a according to blood distribution.
CasesElectrical Resistivity
ρ [Ω·m]
Before blood-filling (cortex)5.84
Before blood-filling (medulla)5.84
After blood-filling (90–10%, cortex)4.99
After blood-filling (90–10%, medulla)5.75
After blood-filling (70–30%, cortex)5.18
After blood-filling (70–30%, medulla)5.56
Table 3. The technical specification of the bio-impedance system.
Table 3. The technical specification of the bio-impedance system.
SpecificationRange
Channels number32
the current3 mA
the sampling frequency500 Hz
frequency band 0.01–137 Hz
the measuring range value1–250 Ω
ADC12 bit
Calibrationintegrated
The accuracy of pulsatile impedance1.0 mΩ
The accuracy of static impedance50 mΩ
Table 4. The electrode system parameters.
Table 4. The electrode system parameters.
The Electrode System2a [mm]2b [mm]
1st3010
2nd4515
3rd6020
Table 5. The obtained electrical impedance values during the numerical modeling.
Table 5. The obtained electrical impedance values during the numerical modeling.
Distance between Measuring Electrodes [mm]Before [Ω]After (90–10%) [Ω]After (70–30%) [Ω]
7394.707394.665394.669
10199.349199.284199.289
13111.898111.801111.808
1670.64270.51370.521
1955.92455.76855.777
2246.87946.69952.492
2543.22443.02643.034
2840.48140.27140.278
3138.89138.67138.678
3437.80237.57637.581
3736.58936.36136.366
4035.72935.49935.503
Table 6. The results of bio-impedance measurement on subject.
Table 6. The results of bio-impedance measurement on subject.
Distance between the Current Electrodes [mm]Base Impedance [Ω]dZ [mΩ]
3095.00232.94
4569.44843.18
6047.25753.42
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Al-harosh, M.; Chernikov, E.; Shchukin, S. Patient Specific Numerical Modeling for Renal Blood Monitoring Using Electrical Bio-Impedance. Sensors 2022, 22, 606. https://doi.org/10.3390/s22020606

AMA Style

Al-harosh M, Chernikov E, Shchukin S. Patient Specific Numerical Modeling for Renal Blood Monitoring Using Electrical Bio-Impedance. Sensors. 2022; 22(2):606. https://doi.org/10.3390/s22020606

Chicago/Turabian Style

Al-harosh, Mugeb, Egor Chernikov, and Sergey Shchukin. 2022. "Patient Specific Numerical Modeling for Renal Blood Monitoring Using Electrical Bio-Impedance" Sensors 22, no. 2: 606. https://doi.org/10.3390/s22020606

APA Style

Al-harosh, M., Chernikov, E., & Shchukin, S. (2022). Patient Specific Numerical Modeling for Renal Blood Monitoring Using Electrical Bio-Impedance. Sensors, 22(2), 606. https://doi.org/10.3390/s22020606

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop