Induced Electromagnetic Fields Estimation in Spine Examinations of MRI Patients: A Re-Evaluation of Existing Clinical Protocols at a Hospital in Accra, Ghana

Show more

1. Introduction

Magnetic Resonance Imaging (MRI) of patients undergoing investigations, are exposed to strong, non-uniform static magnetic fields outside the central imaging region, in which the movement of the body may induce electric currents in tissues which could possibly be harmful [1] . According to Faraday’s Law, MRI scans do expose patients to time varying magnetic fields [2] , which can induce electrical currents that can cause muscle contraction, nerve stimulation, magnetophosphenes, vertigo, nausea, and metallic taste in the mouth [3] .

In Ghana, most patients undergoing MRI investigations are as a result of spinal related injuries. The human spine is one of the delicate organs of the human body because it contains the spinal cord that controls coordination and movement of the human body. As a result, any effect on the spine will have an adverse effect on the body. Hospital data indicate that these injuries may be due to their nature of work―the building industry (carpentry, masonry, laborers), motor mechanics, poor sitting postures, road accidents, carrying of loads on the head, a bad fall, and other such activities. All these activities in one way or the other affect the spine, with the lumbar spine being the most affected because the majority of patients who undergo MRI of the spine complain of lower back pain.

The number of MRI centers in Ghana increased from one in 2008 to ten in 2013. With the increasing number of MRI examinations in Ghana, it is therefore, a quality control measure to ensure that EM field exposure to workers is well below the acceptable maximum limit of 3.8 V/m according to EU Directive 2004/40/EC [4] , which has since been repealed and replaced with the Directive 2013/35/EU of the European Parliament and of the Council of 26 June 2013 [5] . MR patient exposure is covered by the IEC 60601-2-33 directive [6] and limits are maintained by MR systems and these limits are not to be exceeded. Further, there is the requirement that exposure of workers to radiation does not exceed those of the patients. Therefore, there is the need to further understand and evaluate the effects of EM field exposure to patients undergoing MRI examinations in the local Ghanaian setting. This study was undertaken to re-evaluate established clinical protocols, by estimating the induced EM-fields in patients who underwent MRI spine examination. The study was also to serve as an auditing and additional quality control system for training of MRI technicians and medical physicists in Ghana.

Patient participation was voluntary. The procedures involved in the study were verbally explained to patients using an acceptable language of the patient and their written consent sort. Data were anonymously collected and the procedures met the Ethical and Protocol Review Committees of the 37 Military Hospital, Accra, and Ministry of Health, Accra, Ghana which was approved on 12th September, 2012.

2. Materials and Methods

The MRI instrument used was a Philips MR Systems Achieva, (Model No. 43000-574; PHILIPS, Eindhoven, The Netherlands), with a 1.50 T superconducting magnet. The range of magnetic field strength used for medical imaging is typically from 0.15T to 3T [7] . In this study, a magnetic field of 1.50 T was used. The MRI scanner is made up of a long cylindrical hole which is approximately 60 cm in diameter and open on both ends, a sliding table (couch) which rests in the opening and a sense spine coil 1.50 T which is placed under the area of the body being examined.

A simulation test was first performed using the mathematical software MATLAB, Mathworks [8] taking into consideration the following factors: mean length, h, of the spine (41.0 cm was assumed), the mean distance of movement, d, (the displacement) and time, t, represents the duration of movement of the patient on the couch to the central axis of the MRI scanner. The actual evaluation then extended to cover 120 patients who had been referred for MRI of the spine from October 2012 to January 2013, at the 37 Military Hospital, Accra, Ghana. Patient motion through the MR scanner was recorded using a measuring tape and a stop watch. A combined weight, kg, and height, m, scale was also determined. The results were outputted as means ± standard deviations (stddev).

2.1. Model Description and Formulation

2.1.1. Field Distribution of Time Varying Magnetic Fields

A numerical model based on Faraday’s equation and experimental data from MRI spine examinations were used to compute the induced EM fields of the human spine. The computation was carried out using the finite difference method (FDM). The FDM is ideally suited for modeling of electromagnetic waves scattering from complex non-canonical objects [9] and MATLAB software.

The computation method was developed as follows:

The expression for a patient’s motion, v is

$v=\frac{d}{t}$ (1)

where, v is the velocity of the motion, d is distance and t is time.

Static magnetic fields effects are likely to be caused by the magnetic field

vector, B, the gradient or the magnetic field,
$\frac{\partial B}{\partial z}$ , the “force product”, P_{F} may

be expressed as:

${P}_{F}=B\frac{\partial B}{\partial z}$ (2)

where, z is taken as the direction of B.

Flow induced electric field, E_{F}, depends on the product of the flow velocity, v, the magnetic field vector, B and the angle
$\theta $ between them, given as:

${E}_{F}=vB\mathrm{sin}\theta $ (3)

The motion-induced electric field depends on the geometry and the time rate of change of the magnetic field.

$\oint \underset{\_}{E}\cdot \text{d}\underset{\_}{l}}=-{\displaystyle \underset{s}{\int}\frac{\text{d}\underset{\_}{B}}{\text{d}t}\cdot \text{d}\underset{\_}{s}$ (4)

The length vector of the body part enclosing the changing magnetic field is dl and ds the normal vector to the incremental area dA.

2.1.2. Transformation of Faraday’s Equation

Transforming the Faraday‘s equation from Maxwell’s time-dependent equations into difference equations with respect to specific field positions on an elementary lattice, the FD technique provides a simple and effective technique for the modeling of a field distribution within the spine of an MRI patient. The system is then solved in a time marching sequence by alternately calculating the electric and magnetic fields in an interlaced spatial grid cell. Given the location of an electromagnetic source and a complete description of the environment in terms of its dielectric parameters; permittivity (ε), permeability (μ), conductivity (σ), this method provides the ability to assess the electric and magnetic fields.

Therefore, the total magnetic flux passing through a loop is,

${\varphi}_{m}\left(t\right)={\displaystyle \underset{s}{\int}B\times \text{d}s}$ (5)

Electric potential can also be expressed as;

$v\left(t\right)={\displaystyle \oint E\times \text{d}l}$ (6)

But from Faraday’s law;

$v\left(t\right)=-\frac{\text{d}{\varphi}_{m}\left(t\right)}{\text{d}t}$ (7)

where, ${\varphi}_{m}\left(t\right)$ is the total time varying magnetic flux passing through the surface.

Considering a close cylindrical loop, substituting Equation (5) into Equation (7),

$v\left(t\right)=-{\displaystyle \int \frac{B}{\text{d}t}\times \text{d}s}$ (8)

$v\left(t\right)=-\frac{\text{d}}{\text{d}t}{\displaystyle \int B\times \text{d}s}$ (9)

where, ds = Area (A).

The differential equation of

$\frac{{\partial}^{2}A}{\partial {r}^{2}}=2\text{\pi}{r}^{2}h$

Therefore, Equation (9) becomes

$v\left(t\right)=-\frac{\text{d}}{\text{d}t}{\displaystyle \int B\times A}$ (10)

$=-\frac{\partial B}{\partial t}\times \frac{\partial A}{\partial t}$ (11)

From Equations (6) and (9)

$\oint E\times \text{d}l}=-\frac{\text{d}}{\text{d}t}{\displaystyle \int B\times \text{d}s$ (12)

Equation (11) is the integral form of the Faraday’s Law.

Using Stokes’ theorem [10] to transform the line integral into a surface integral;

$\oint \underset{\_}{E}\times \text{d}l}={\displaystyle \underset{s}{\int}\left(\underset{\_}{\nabla}\times \underset{\_}{E}\right)\times \text{d}a$ (13)

Therefore, the differential form of the Faraday’s law of induction is;

$\underset{\_}{\nabla}\times \underset{\_}{E}=-\frac{\partial \underset{\_}{B}}{\partial t}$ (14a)

$\begin{array}{l}\underset{\_}{\nabla}\times \underset{\_}{E}=\left(\frac{1}{r}\frac{\partial {E}_{z}}{\partial \theta}-\frac{\partial {E}_{\theta}}{\partial z}\right)\stackrel{^}{r}+\left(\frac{\partial {E}_{r}}{\partial z}-\frac{\partial {E}_{z}}{\partial r}\right)\stackrel{^}{\theta}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{1}{r}\left(\frac{\partial}{\partial r}\left(r{E}_{\theta}\right)-\frac{\partial {E}_{r}}{\partial \theta}\right)\stackrel{^}{z}\end{array}$ (14b)

Note that the determinant;

$\underset{\_}{\nabla}\times \underset{\_}{E}=\left[\begin{array}{ccc}\stackrel{^}{r}& r\stackrel{^}{\theta}& \stackrel{^}{z}\\ \frac{\partial}{\partial r}& \frac{\partial}{\partial \theta}& \frac{\partial}{\partial z}\\ {E}_{r}& rE\theta & {E}_{z}\end{array}\right]$

Formulating EM wave equation of the spine from Maxwell’s equation;

$\nabla \times B=\varphi \nabla \times E=0$ (15)

$\underset{\_}{\nabla}\times \underset{\_}{E}=-\frac{\partial \underset{\_}{B}}{\partial t}\underset{\_}{\nabla}\times \underset{\_}{B}={\epsilon}_{\circ}{\mu}_{\circ}\frac{\partial \underset{\_}{E}}{\partial t}$ (16)

$\underset{\_}{\nabla}\times \underset{\_}{\nabla}\times \underset{\_}{E}=-\underset{\_}{\nabla}\times \frac{\partial \underset{\_}{B}}{\partial t}$ (17)

$\underset{\_}{\nabla}\times \underset{\_}{\nabla}\times \underset{\_}{E}=\underset{\_}{\nabla}\left(\underset{\_}{\nabla}\cdot E\right)-{\underset{\_}{\nabla}}^{2}\underset{\_}{E}$ (18)

$-\frac{\partial}{\partial t}\underset{\_}{\nabla}\times \underset{\_}{\nabla}=\underset{\_}{\nabla}\left(\underset{\_}{\nabla}.\underset{\_}{\nabla}\right)-{\underset{\_}{\nabla}}^{2}\underset{\_}{E}$ (19)

$-\frac{\partial \underset{\_}{B}}{\partial t}\underset{\_}{\nabla}\times \underset{\_}{B}={\epsilon}_{\circ}{\mu}_{\circ}\frac{{\partial}^{2}\underset{\_}{E}}{\partial {t}^{2}}$ (20)

Combing Equations (19) and (20);

$\underset{\_}{\nabla}\left(\underset{\_}{\nabla}\times \underset{\_}{\nabla}\right)-{\underset{\_}{\nabla}}^{2}\underset{\_}{E}=-{\epsilon}_{o}{\mu}_{o}\frac{{\partial}^{2}\underset{\_}{E}}{\partial {t}^{2}}$ (21)

where, $\nabla \left(\nabla \times E\right)=0$

$-{\underset{\_}{\nabla}}^{2}\underset{\_}{E}=-{\epsilon}_{o}{\mu}_{o}\frac{{\partial}^{2}\underset{\_}{E}}{\partial {t}^{2}}$ (22)

Equation (22) indicates the modified Maxwell’s equation, where, ${\epsilon}_{\circ}{\mu}_{\circ}=\frac{1}{{v}^{2}}.$

Therefore, Equation (22) becomes

$-\underset{\_}{\nabla}\underset{\_}{E}=-\frac{1}{{v}^{2}}\frac{{\partial}^{2}E}{\partial {t}^{2}}$ (23)

Since the EM field is polarized in the y direction, the induced EM field in the 3D plane within the spine becomes;

From Equation (14) ${\varphi}_{m}\left(t\right)={\displaystyle \underset{s}{\int}B}\times \text{d}s$

$\frac{{\partial}^{2}{\varphi}_{m}}{\partial {r}^{2}}+\frac{{\partial}^{2}{\varphi}_{m}}{\partial {z}^{2}}=\frac{1}{2\text{\pi}rh}\frac{{\partial}^{2}{\varphi}_{m}}{\partial {t}^{2}}$ (24)

The discretization of Equation (25) yields the EM field within various positions of the spine.

$2\text{\pi}h\left(\frac{{\partial}^{2}{\varphi}_{m}}{\partial {r}^{2}}+\frac{{\partial}^{2}{\varphi}_{m}}{\partial {z}^{2}}\right)=\frac{{\partial}^{2}{\varphi}_{m}}{\partial {t}^{2}}$ (25)

where, h is the length of the spine.

2.1.3. Discretization of Spine Model

The discretization of the partial differential equation, Equation (25) is illustrated in Figure 1, giving rise to the solution, Equation (26), where, i, j and k represent variables temperature, length and time, respectively.

${\varphi}_{i,j,k}=\frac{{\displaystyle \underset{m=0}{\overset{1}{\sum}}\left({\varphi}_{i-m,j,k}-2{\varphi}_{i,j,k}+{\varphi}_{i+m,j,k}\right)}}{\Delta {r}^{2}}+\frac{{\displaystyle \underset{m=0}{\overset{1}{\sum}}\left({\varphi}_{i,j-m,k}-2{\varphi}_{i,j,k}+{\varphi}_{i,j+1,k}\right)}}{\Delta {z}^{2}}$ (26)

Therefore, the final discretized EM-field equation from the Faraday’s equation becomes;

Figure 1. The discretization of the partial differential Equation (25) into discrete counterparts.

$\Rightarrow \frac{{\varphi}_{i,j}^{n+1}-2{\varphi}_{i,j}^{n}+{\varphi}_{i,j}^{n-1}}{\left(\Delta {t}^{2}\right)}={v}^{2}\omega \left[\frac{{\varphi}_{i+1,j}^{n}-2{\varphi}_{i,j}^{n}+{\varphi}_{i-1,j}^{n}}{{\left(\Delta h\right)}^{2}}+\frac{{\varphi}_{i,j+1}^{n}-2{\varphi}_{i,j}^{n}+{\varphi}_{i,j-1}^{n}}{{\left(\Delta h\right)}^{2}}\right]$ (27)

where, ϕ^{n}, ϕ^{n}^{+1} represent the iterated EM fields within the spatial grids of the spine and
$\Delta h$ is the step size.

3. Results and Discussions

3.1. Simulation Results

The induced EM fields were computed for the three sections of the spine as the patients underwent the MRI examinations using FDM and MATLAB. The simulated results of the induced electric field (0.39 V/m) and current density (0.039 A/m^{2}) were lower than the measured results (0.44 V/m and 0.043 A/m^{2}, respectively) and this was due to the fact that during the simulation an assumed length of the spine (41.0 cm) was assumed for all patients.

3.2. Experimental Results

A total of 120 patients were MRI examined of which 75 patients ranging in ages from 8 to 90 years underwent lumbar spine examinations (Table 1). The study revealed that most of the patients undergoing the lumbar spine examination were mainly drivers, civil servants, artisans, and traders, whose activities involve several hours of sitting or standing for extended periods (Table 1). There were a few cases in which the patients were involved in motor accidents. The lumbar spine had the largest displacement with a mean (±stddev) value of 1.46 ± 0.06 m compared with the thoracic spine, 1.33 ± 0.36 m and the cervical spine, 1.10 ± 0.10 m (Table 2).

The lumbar spine had the highest duration with a mean of 5.98 ± 0.48 s compared with the cervical spine, 5.92 ± 0.40 s and the thoracic spine, 5.88 ± 0.67 s.

Table 1. Gender, age range and number of patients for the various spine examinations.

Table 2. Results of distance and time measurements.

The results showed that, a maximum velocity of 0.29 m/s, with peak induced electric field and current density of 0.44 V/m and of 0.043 A/m^{2} were recorded, respectively. The results of the induced electric field were compared to the limits given in the directive IEC 60601-2-33 standard [6] . The comparison shows that the derived peak electric field of 0.44 V/m was much lower than the IEC standard of 3.8 V/m in the normal operating mode and 4.7 V/m in the first level operating mode. These values are lower than the onset thresholds of 6 - 7 V/m for peripheral nerve stimulation (PNS). According to Directive 2013/35/EU [5] for EMF, exposure limit of current density for the head and torso is 3 A/m^{2}. The result of the induced current densities in the spinal cord was 0.043 A/m^{2} which as much lower than the limit of 3 A/m^{2}. Similarly, the induced current densities in the skin and spine were lower than the Directive 2013/35/EU exposure limit [5] .

Solutions to the discretized partial differential equations obtained with the finite difference method showing the variations in induced EM fields for male (Figure 2 and Figure 3) and female (Figure 4 and Figure 5) patients, who underwent the MRI examination of the lumbar spine are displayed in two-dimension (2D) and three-dimension (3D). For the 2D profiles for males (Figure 2), the mean induced EM field of 0.24 V/m (range 0.0 to 0.38 V/m) was recorded while the female profile (Figure 4), recorded a mean value of 0.22 V/m and ranges from 0.0 to 0.36 V/m. For the 3D profiles, the induced EM fields were static with time, t and varied with position, r. For both the 3D profiles, concentration of the induced EM fields was evenly distributed which depicts a hyperbolic shape. The color maps represent the concentrations of the EM fields; the blue color represents low EM field region whiles the red color represents high EM field region. Comparing the 3D profiles for both genders, the female lumbar (Figure 5) was more highly concentrated than the male profile (Figure 3).

Figure 2. Variations in induced EM-Field for male lumbar spine in 2D.

Figure 3. Variations in induced EM-Field for male lumbar spine in 3D profiles.

Figure 4. Variations in induced EM-Field for female lumbar spine in 2D profiles.

Figure 5. Variations in induced EM-Field for female lumbar spine in 3D profiles.

From this study, the simulated results of the induced electric field (0.39 V/m) and current density (0.039 A/m^{2}) the experimental results; induced electric field (0.44 V/m) and current density (0.043 A/m^{2}) were much lower compared with the report by Liu et al. [1] , where the simulated results of the induced electric field and current density (1.8 V/m and 0.22 A/m^{2} respectively) and the experimental results (1.8 V/m and 0.21 A/m^{2} respectively). Here, a lower magnetic field strength of 1.5 T was used as compared with the 4.0 T used for simulation, with the maximum current of about 220 mA/m^{2} in the report by Liu et al., [1] . In modern MRI scanners, Faraday currents can reach a value of 386 mA/m^{2} [11] . Further, Bowtell and Bowley [12] , indicated that, the electric field peak value of about 1.8 V/m (for the 4 T MRI scanner, with a corresponding duration of movement of 0.5 m/s movement), is much lower than 6.2 V/m, which has been suggested to be the threshold for nerve stimulation at low frequencies [5] .

4. Conclusion

The results of the level of induced EM-fields were such that they might not pose any potential health hazards to the patients as these values were well below the recommended guidance levels set by Directives 2013/35/EU [5] and IEC 60601-2-33 standards [6] exposure limits for workers and patients, respectively.

We also recommend periodic reviews of such protocols of newer NMR instruments to help shape safety policies for both patients and workers.

Conflict of Interest

The authors declare that there are no conflicts of interests regarding the study and publication of this paper.

Acknowledgements

We would like to express our sincere thanks to Prof. AWK Kyere, Head of Medical Physics Department, Graduate School of Nuclear and Allied Sciences, Legon, Ghana, the staff of the MRI Unit at 37 Military Hospital and Ghana Standards Authority.

References

[1] Liu, F., Zhao, H. and Crozier, S. (2003) Calculation of the Electric Fields Induced by Body and Head in High-Field MRI. Journal of Magnetic Resonance, 161, 99-107.

http://pfeifer.phas.ubc.ca/refbase/files/Liu-JMagRes-2003-161-99.pdf

https://doi.org/10.1016/S1090-7807(02)00180-5

[2] Mcrobbie, D.W. (2012) Occupational Exposure in MRI. British Journal of Radiology, 85, 293-312.

https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3486652/

https://doi.org/10.1259/bjr/30146162

[3] ICNIRP (2014) Guidelines for Limiting Exposure to Electric Fields Induced by Movement of the Human Body in a Static Magnetic Field and by Time-Varying Magnetic Fields below 1 Hz. Health Physics, 106, 418-425.

https://www.icnirp.org/cms/upload/publications/ICNIRPmvtgdl_2014.pdf

https://doi.org/10.1097/HP.0b013e31829e5580

[4] Directive 2004/40/EC of the European Parliament and of the Council (2004). Official Journal of the European Union, L184/1, L184-9.

https://archive.senseaboutscience.org/data/files/PhysicalAgentsDirective.pdf

[5] Directive 2013/35/EU-Electromagnetic Fields.

https://osha.europa.eu/en/legislation/directives/directive-2013-35-eu-electromagnetic-fields

[6] IEC 60601-2-33:2010+AMD1:2013+AMD2:2015 CSV (2010) Medical Electrical Equipment—Part 2-33: Particular Requirements for the Basic Safety and Essential Performance of Magnetic Resonance Equipment for Medical Diagnosis. Consolidated Version, 2010, 2013, 2015.

https://webstore.iec.ch/preview/info_iec60601-2-33%7Bed3.0%7Db.pdf

[7] MR-TIP (2012) MRI.

http://www.mr-tip.org/mri

[8] MATLAB (2013) Mathworks.

http://mathworks.com/products/matlab

[9] Taflove, A. and Hagness, S.C. (2005) Computational Electrodynamics: The Finite-Difference Time-Domain Method. 3rd Edition, Artech House, Norwood, MA.

[10] Green, P. and Rosenberg, J. (2010) Line Integrals around Closed Curves, and the Theorems of Green and Stokes.

https://qmplus.qmul.ac.uk/pluginfile.php/283380/mod_resource/content/3/Chapter%204%20-%20Vector%20intergration%2C%202010.pdf

[11] Schmitt, F., Stehling, M.K. and Turner, R. (1998) Echo-Planar Imaging Theory, Technique and Application. Springer, New York, 9-14.

http://radioramanuja.com/file/download/echo-planar-imaging-theory-technique-and-application

https://doi.org/10.1007/978-3-642-80443-4

[12] Bowtell, R. and Bowley, R.M. (2000) Analytic Calculations of the E-Fields Induced by Time-Varying Magnetic Fields Generated by Cylindrical Gradient Coils. Magnetic Resonance in Medicine, 44, 782-790.

https://doi.org/10.1002/1522-2594(200011)44:5<782::AID-MRM16>3.0.CO;2-7