Abstract
In this study, the Visual Monte Carlo radiation transport code and the female voxel phantom FAX were used to calculate organ and effective doses delivered by target–source irradiation geometries associated with radioiodine therapy treatments. Methods: Specific situations were considered: when a patient was accompanied during hospitalization, when a patient was accompanied on return to his or her residence, and when a patient received daily care at home. Results: This simulation study showed that, in the 3 situations considered, the total effective dose to an individual in normal contact with the patient was less than 0.85 mSv for up to 11.1 GBq (300 mCi) of administered activity. Conclusion: The results of this study suggest that for these patients receiving radioiodine therapy, radiation protection procedures after hospital discharge are unnecessary.
When 131I is administered to a patient who has undergone a total thyroidectomy, the patient is typically hospitalized for at least 2 d after the administration to allow the dose rate to fall to acceptable levels. The patient is instructed to maintain a certain distance from family members for several days after hospital discharge (1). A total thyroidectomy removes almost all of the thyroid gland tissue. To achieve the total elimination of thyroid tissue, 131I activity is administered in the range of 3.7 GBq (100 mCi) to 11.1 GBq (300 mCi) (2). This activity is concentrated in the remaining thyroid tissue; the high localized dose causes the death of the remaining thyroid tissue cells. Current radiation protection procedures establish that nonradiation worker medical staff and caregivers should not receive an effective dose of greater than 5 mSv for each treatment and that members of the general public should not receive an effective dose of greater than 1 mSv when exposed to a patient after treatment (3).
At present, the effective dose to a caregiver is calculated by considering the 131I to be concentrated in a point source in a patient, without shielding. The effective dose to the caregiver and members of the public is also calculated from a point source near the patient. This approach does not take into consideration the attenuation and scattering of the radiation by the human body. Other studies have used the line-source approximation instead of the point-source model for estimating the dose to an individual exposed to a patient who has received radioactivity (4), and anthropomorphic mathematic phantoms with Monte Carlo transport techniques have been proposed for a more accurate dose calculation (5,6). The present study proposes a more realistic method for calculating radiation transport between patients and exposed individuals: the Monte Carlo method and anthropomorphic voxel phantoms.
To evaluate the doses received by caregivers and the general public, 3 exposure situations were considered. In the first, the caregiver spends one night close to the patient while hospitalized; in the second, the patient is accompanied by the caregiver during the journey from hospital to home; and in the third, the caregiver accompanies the patient in his or her daily routine. The calculations were done with mean biokinetic parameters obtained experimentally (7,8); with the Visual Monte Carlo code (9), which simulates radiation transport by use of the Monte Carlo method; and with the female adult voxel phantom (FAX) (10), which is based on CT images and represents more realistic anatomic characteristics. The female adult voxel phantom FAX was used as the patient (source) and the target in the simulations because thyroid cancer is more common in women than in men.
MATERIALS AND METHODS
Voxel Phantom FAX
This phantom presents anatomic structures comparable to those described for the “reference woman” in ICRP Publication 89 (11). The phantom consists of 453 slices obtained from a CT scan, resulting in a voxel phantom of an adult woman with a height of 1.63 m and a weight of 59.76 kg. FAX includes 22 segmented organs and tissues and a total of 1,272,504 voxels of sides measuring 0.36 cm × 0.36 cm × 0.36 cm.
Visual Monte Carlo Code
Various codes that simulate photon transport through nonhomogeneous structures have been written; among these codes are GEANT4 (12), MCNP (13), and FLUKA (14). The Visual Monte Carlo program was chosen for this study because it is written in the programming language Visual Basic version 6 and has visualization resources through its graphic interface. This software has already been applied to areas of radiation protection, including the calibration of in vivo measurement systems, calculations of doses from external sources of radiation, and calculations of specific effective energies (15).
Calculation of Cumulated Activity
Mean biokinetic parameters that describe the processes of uptake, retention, and clearance of 131I were obtained experimentally. Among these were coefficients of the biokinetic equations that describe the behavior of 131I in a patient's body. These coefficients were obtained from measurements of 131I activities in thyroid remnants in various patients by the SPECT technique (7,8). Activities calculated for each patient were plotted as a function of time, and the uptake and elimination phases were fitted with nonlinear regression to a 2-component exponential function.
The activity of 131I in the thyroid and the time that it stays in this region depend on the properties and biologic characteristics of the radiopharmaceutical and of the individual's metabolism as well as on the physical half-life of the radionuclide. The total number of all nuclear transformations in this organ over a given time interval is defined here as Ã:(Eq. 1)In Equation 1, A(t) is the activity in the thyroid as a function of time. The cumulated activity in any particular source organ is calculated by integrating the instantaneous activity over time. A(t) can be described as the subtraction of 2 exponential terms:(Eq. 2)In Equation 2, t is the time since the ingestion of 131I, a0 is the initial maximum activity in the thyroid, and a1 and a2 are obtained experimentally. a0 is calculated as follows:(Eq. 3)In Equation 3, f is the maximum fraction captured by residues of the thyroid, A0 is the initial activity administered to the patient, and constant C is calculated as follows:Substituting Equation 2 into Equation 1 yields:(Eq. 4)
In this work, a1 and a2 values of 0.0129 and 0.577 h−1 were used; these values were experimentally derived from patients who had undergone a total thyroidectomy and had right and left thyroid residues with a total mass of 3 g. The average percentage uptake in the remaining thyroid tissue was 3.8% (f) of 131I. A total activity intake of 11.1 GBq (300 mCi) (A0) of 131I via ingestion was assumed; this value represented the most conservative dose. The calculations were made on the basis of the assumptions that all of the 131I was concentrated in the thyroid remnants and that the rest of the 131I was quickly eliminated through the urine. These assumptions were based on the biokinetics of iodine in humans; that is, after 48 h, all of the iodine present in the body has been taken up by the remaining thyroid tissue (1).
Irradiation Geometries
Three irradiation geometries representing conservative situations that a patient and a caregiver would be likely to experience during and after hospitalization were considered. The first exposure geometry simulated a situation in which the caregiver stayed in the patient's room in the hospital overnight. Before sleeping, the patient ingested 131I at 07:00 am and then began sleeping at 21:00 pm. This time interval represented 14 h of biologic and physical elimination of 131I from thyroid residues. From 21:00 pm to 05:00 am, the patient and the caregiver slept in the same room and were assumed to be stationary during this 8-h sleeping period; this is the period that was simulated. The value of the cumulated activity in the thyroid residues in this case was 1.07 × 1013 Bq·s. The exposure geometry was assumed to be left lateral–right lateral for the 2 voxel phantoms. The patient phantom (source) and the caregiver phantom (target) were assumed to be side by side at a distance of 1.0 m from inside arm to inside arm (Fig. 1).
In the second exposure geometry, which simulated a return home, the patient returned home in the company of a caregiver or the general public, such as a taxi driver. In this situation, 2 scenarios were considered: the patient beside the caregiver or taxi driver at a distance of 0.3 m and the patient behind the caregiver or taxi driver at a distance of 0.5 m (Fig. 2). The exposure time for this simulation was the time required to travel from the hospital to the home. It was assumed that the patient was discharged from the hospital after 2 d of treatment. Cumulated activities were calculated for various travel times (Table 1).
The third exposure geometry simulated a situation in which the patient was in close contact (a distance of 0.3 m) with the caregiver for 16 h/d on days 2–7 after the intake of 131I. A left lateral–right lateral exposure geometry was assumed for the 2 voxel phantoms. After current procedures, the patient is normally instructed to resume “normal contact” only 7 d after intake. Cumulated activities (number of nuclear transformations) were calculated for various days (Table 2).
For all of the exposure geometries, 108 photon histories were used for the Visual Monte Carlo code runs to achieve an error in organ dose assessment of less than 5%.
RESULTS
All effective dose calculations were made with the tissue weighting factors given in ICRP Publication 103 (16). Calculated equivalent organ doses to the caregiver during the overnight stay in the hospital are shown in Figure 3. An effective dose of 0.05 mSv was calculated for the caregiver.
Table 3 shows the calculated effective doses for an individual accompanying a patient immediately after discharge from the hospital. The highest dose was observed for the 5-h journey with the patient behind the driver at a distance of 0.5 m.
Table 4 shows the calculated effective doses for a caregiver in close contact with a patient.
DISCUSSION
The present study was the first step in elaborating an alternative method for estimating doses from patients receiving radioiodine therapy by use of the Monte Carlo method and anthropomorphic voxel phantoms. One limitation of our study was the use of only the thyroid and not the whole body as the source of radioactivity. This choice was made because only experimental data (biokinetic parameters) for thyroid uptake from a specific group of patients were available to us. In the future, we intend to develop new irradiation geometries with more experimental data (from the thyroid and the rest of the body), other therapeutic radiopharmaceuticals, and combinations of adult male and female phantoms.
CONCLUSION
The results of this work show that the calculated effectives doses were less than 0.85 mSv when doses attributable to all of the simulated exposure geometries were considered. This finding implies that even for the overnight stay of the caregiver, the presence of the caregiver on a 5-h trip home, and the close contact of the caregiver with the patient at home for 16 h/d on days 2–7, the total effective doses were considerably lower than the 5-mSv limit established in U.S. Nuclear Regulatory Commission Regulatory Guide 8.39 and lower than the 1-mSv limit above which release instructions are required (3). This finding suggests that for these patients receiving radioiodine therapy, radiation protection procedures after hospital discharge are unnecessary. However, other radiation protection measures, such as control of urine, the primary via of excretion of iodine, must be implemented to minimize the spread of contamination when such patients are released. The relaxing of conservative radiation protection procedures will result in psychologic benefits for patients and their families.
Footnotes
-
COPYRIGHT © 2009 by the Society of Nuclear Medicine, Inc.
References
- Received for publication September 15, 2008.
- Accepted for publication November 12, 2008.