Attenuation Coefficients for Layered Ceiling and Floor Shields in PET/CT Clinics

Similar documents
A Monte Carlo shielding model for PET/CT clinics

IMPLEMENTATION OF SALIVARY GLANDS IN THE BODYBUILDER ANTHROPOMORPHIC PHANTOMS

Automated ADVANTG Variance Reduction in a Proton Driven System. Kenneth A. Van Riper1 and Robert L. Metzger2

I. INTRODUCTION. Figure 1. Radiation room model at Dongnai General Hospital

Validation of GEANT4 for Accurate Modeling of 111 In SPECT Acquisition

ELECTRON DOSE KERNELS TO ACCOUNT FOR SECONDARY PARTICLE TRANSPORT IN DETERMINISTIC SIMULATIONS

Dosimetry Simulations with the UF-B Series Phantoms using the PENTRAN-MP Code System

Absorbed fractions in a voxel-based phantom calculated with the MCNP-4B code

Application of MCNP Code in Shielding Design for Radioactive Sources

Mesh Human Phantoms with MCNP

Electron Dose Kernels (EDK) for Secondary Particle Transport in Deterministic Simulations

Outline. Monte Carlo Radiation Transport Modeling Overview (MCNP5/6) Monte Carlo technique: Example. Monte Carlo technique: Introduction

MCNP4C3-BASED SIMULATION OF A MEDICAL LINEAR ACCELERATOR

Michael Speiser, Ph.D.

Image-based Monte Carlo calculations for dosimetry

Evaluation of RayXpert for shielding design of medical facilities

UNIVERSITY OF SOUTHAMPTON

Improved Detector Response Characterization Method in ISOCS and LabSOCS

Three Dimensional Dosimetry Analyses In Radionuclide Therapy Using IDL And MCNP-based Software Tools

Effects of the difference in tube voltage of the CT scanner on. dose calculation

Av. Professor Lineu Prestes São Paulo, SP ABSTRACT

Breaking Through the Barriers to GPU Accelerated Monte Carlo Particle Transport

Introduction to Positron Emission Tomography

Radiation Oncology treatment facility design. Simulators, CT scanners, HDR Brachytherapy

Modeling Radiation Transport Using MCNP6 and Abaqus/CAE Chelsea A. D Angelo, Steven S. McCready, Karen C. Kelley Los Alamos National Laboratory

Shielding factors for traditional safety glasses

Suitability Study of MCNP Monte Carlo Program for Use in Medical Physics

SESAME: A NEW TOOL FOR THE NUMERICAL DOSIMETRIC RECONSTRUCTION OF A RADIOLOGICAL ACCIDENT CAUSED BY AN

CALCULATION OF DOSE DISTRIBUTION AND DOSE RATE AT EACH DWELL POSITION IN 60 CO IRRADIATOR SVST-CO-60/B IN VIETNAM

Ch. 4 Physical Principles of CT

Artifact Mitigation in High Energy CT via Monte Carlo Simulation

Limitations in the PHOTON Monte Carlo gamma transport code

Proton dose calculation algorithms and configuration data

THE SIMULATION OF THE 4 MV VARIAN LINAC WITH EXPERIMENTAL VALIDATION

Medical Physics Research Center, Mashhad University of Medical Sciences, Mashhad, Iran.

Comparison of internal and external dose conversion factors using ICRP adult male and MEET Man voxel model phantoms.

CALCULATION OF THE ACTIVITY INVENTORY FOR THE TRIGA REACTOR AT THE MEDICAL UNIVERSITY OF HANNOVER (MHH) IN PREPARATION FOR DISMANTLING THE FACILITY

Radiology. Marta Anguiano Millán. Departamento de Física Atómica, Molecular y Nuclear Facultad de Ciencias. Universidad de Granada

OBJECT ORIENTED DESIGN OF ANTHROPOMORPHIC PHANTOMS AND GEANT4-BASED IMPLEMENTATIONS ABSTRACT

SHIELDING DEPTH DETERMINATION OF COBALT PHOTON SHOWER THROUGH LEAD, ALUMINUM AND AIR USING MONTE CARLO SIMULATION

Corso di laurea in Fisica A.A Fisica Medica 5 SPECT, PET

UNCOMPROMISING QUALITY

ANALYSIS OF CT AND PET/SPECT IMAGES FOR DOSIMETRY CALCULATION

MCNP Variance Reduction technique application for the Development of the Citrusdal Irradiation Facility

Modeling the ORTEC EX-100 Detector using MCNP

ISOCS Characterization of Sodium Iodide Detectors for Gamma-Ray Spectrometry

Basics of treatment planning II

OPTIMIZATION OF MONTE CARLO TRANSPORT SIMULATIONS IN STOCHASTIC MEDIA

Construction of Voxel-type Phantom Based on Computed Tomographic Data of RANDO Phantom for the Monte Carlo Simulations

Investigation of tilted dose kernels for portal dose prediction in a-si electronic portal imagers

Monte Carlo Simulation for Neptun 10 PC Medical Linear Accelerator and Calculations of Electron Beam Parameters

Combining Analytical and Monte Carlo Modelling for Industrial Radiology

Production of neutrons in laminated barriers of radiotherapy rooms: comparison between the analytical methodology and Monte Carlo simulations

Spectral analysis of non-stationary CT noise

REMOVAL OF THE EFFECT OF COMPTON SCATTERING IN 3-D WHOLE BODY POSITRON EMISSION TOMOGRAPHY BY MONTE CARLO

Towards patient specific dosimetry in nuclear medicine associating Monte Carlo and 3D voxel based approaches

Reconstruction in CT and relation to other imaging modalities

Neutronics Analysis of TRIGA Mark II Research Reactor. R. Khan, S. Karimzadeh, H. Böck Vienna University of Technology Atominstitute

DUAL energy X-ray radiography [1] can be used to separate

Monte Carlo methods in proton beam radiation therapy. Harald Paganetti

DETERMINISTIC 3D RADIATION TRANSPORT SIMULATION FOR DOSE DISTRIBUTION AND ORGAN DOSE EVALUATION IN DIAGNOSTIC CT

Impact of X-ray Scatter When Using CT-based Attenuation Correction in PET: A Monte Carlo Investigation

Dose Distributions. Purpose. Isodose distributions. To familiarize the resident with dose distributions and the factors that affect them

CT vs. VolumeScope: image quality and dose comparison

Monte Carlo Method for Medical & Health Physics

Investigation of Scattered Radiation Dose at the Door of a Radiotherapy Vault When the Maze Intersects the Primary Beam

A PRACTICAL LOOK AT MONTE CARLO VARIANCE REDUCTION METHODS IN RADIATION SHIELDING

Graphical User Interface for Simplified Neutron Transport Calculations

A dedicated tool for PET scanner simulations using FLUKA

Detailed analysis of scatter contribution from different simulated geometries of X-ray detectors

Basic Radiation Oncology Physics

Comparison of absorbed dose distribution 10 MV photon beam on water phantom using Monte Carlo method and Analytical Anisotropic Algorithm

Room scatter effects in Total Skin Electron Irradiation: Monte Carlo simulation study

3/27/2012 WHY SPECT / CT? SPECT / CT Basic Principles. Advantages of SPECT. Advantages of CT. Dr John C. Dickson, Principal Physicist UCLH

An Investigation into the Limitations of Tomographic Scanning for Quantitative Nondestructive Gamma-ray Measurements of High Density Radioactive Waste

Muon imaging for innovative tomography of large volume and heterogeneous cemented waste packages

Shielding Design Methods for Radiation Oncology Departments

Mathematical methods and simulations tools useful in medical radiation physics

Disclosure 7/24/2014. Validation of Monte Carlo Simulations For Medical Imaging Experimental validation and the AAPM Task Group 195 Report

VISIPLAN 3D ALARA planning tool Version 3.0

OPTIMIZATION OF MONTE CARLO TRANSPORT SIMULATIONS IN STOCHASTIC MEDIA

C a t p h a n / T h e P h a n t o m L a b o r a t o r y

IMPROVEMENTS TO MONK & MCBEND ENABLING COUPLING & THE USE OF MONK CALCULATED ISOTOPIC COMPOSITIONS IN SHIELDING & CRITICALITY

GPU implementation for rapid iterative image reconstruction algorithm

3-D PET Scatter Correction

Indrin Chetty Henry Ford Hospital Detroit, MI. AAPM Annual Meeting Houston 7:30-8:25 Mon 08/07/28 1/30

Introduction to Biomedical Imaging

Diagnostic imaging techniques. Krasznai Zoltán. University of Debrecen Medical and Health Science Centre Department of Biophysics and Cell Biology

VIEWING MORSE-CG RADIATION TRANSPORT WITH. 3-D COLOR GRAPHICS* Abstract

Total Measurement Uncertainty Estimation For Tomographic Gamma Scanner R.Venkataraman 1, S.Croft 1, M.Villani 1, R.D. McElroy 1, and R.J.

MapCHECK 2 & 3DVH. The Gold Standard for 2D Arrays

CT Protocol Review: Practical Tips for the Imaging Physicist Physicist

Identification of Shielding Material Configurations Using NMIS Imaging

SPECT QA and QC. Bruce McBride St. Vincent s Hospital Sydney.

Large Plastic Scintillation Detectors for the Nuclear Materials Identification System

Monte Carlo modelling and applications to imaging

Simulation of Radiographic Testing for POD Assessment

Assesing multileaf collimator effect on the build-up region using Monte Carlo method

Design and performance characteristics of a Cone Beam CT system for Leksell Gamma Knife Icon

Basics of treatment planning II

Transcription:

Attenuation Coefficients for Layered Ceiling and Floor Shields in PET/CT Clinics Robert L. Metzger and Kenneth A. Van Riper Radiation Safety Engineering, Inc 3245 North Washington Street, Chandler, AZ 85225 rlmetzger@radsafe.com; kvr@rt66.com INTRODUCTION Positron Emission Tomography (PET) is enjoying explosive growth due to its ability to accurately stage many types of cancer and follow the progress of treatments. The facilities present unique challenges in shielding design due to the nature of the study, and the desire to place the facilities near other imaging and treatment clinics. A typical PET patient will receive 555 MBq (15 mci) of 18 F labeled 2-Fluoro-2- Deoxy-D-Glucose ( 18 FDG) and will rest in a darkened room (called a quiet room) for 45 minutes to allow the drug to localize in the lesions of interest. A normal clinic will have one to three quiet rooms. While in the quiet room, the patient will normally recline in a comfortable lounge chair. The distribution of the 18 FDG in the patient will vary from patient to patient, but much of the isotope will be located in the brain and the bladder. After resting for approximately 45 minutes, the patient is asked to empty his or her bladder, and is placed on the scanner bed. Images are collected for 30 to 45 minutes on the PET or combined PET/CT unit. The patient is then released after the scanning is complete. Shielding design for PET/CT clinics is frequently difficult as many of the facilities are being retrofitted into existing imaging centers and hospitals where space is cramped. The patient quiet rooms, the hot laboratory, and the scanners are commonly positioned in close proximity to uncontrolled areas where non-occupational dose limits apply. Of particular concern are the ceiling and floor shields for the quiet and the scanner rooms because these are constructed by adding sheets of lead to the existing concrete ceiling and floor decks. Much of the cost of shielding in a PET/CT clinic is driven by the installation costs of these shields. The attenuation provided by these layered shields is difficult to calculate by point kernel methods. An underestimation of the shielding effectiveness will result in additional lead and higher cost. CLINIC ROOM MODELS Wall shielding is constructed of commercial sheet lead and can easily be calculated using point kernel techniques [1, 2] with appropriate buildup factors. When occupied areas exist above or below the PET clinic, however, the shielding is constructed of the existing concrete floor deck and commercial sheet lead applied to the top of the floor, or, more commonly, suspended from the floor joists. Concrete floor decks in medical office buildings typically range from 5 to 15 cm in thickness, while hospitals designed for heavy floor loading can exceed 20 cm in thickness. The concrete floor decks in existing Page 1

structures may be constructed of normal density concrete (2.35 g/cm3) or lightweight concrete, and the poured thickness may vary from the specified amount by 1 cm or more. The layered lead and concrete shield is difficult to estimate by point kernel techniques as the energy of the photons incident on each layer is difficult to determine, as is the pass off of scattered photons from layer to layer. Nonetheless, the accuracy of the floor and ceiling shielding estimate is important as much of the cost of the total clinic shielding is driven by the construction costs related to the suspended lead. Consequently, three-dimensional Monte Carlo models for MCNP5 [3] were developed for the quiet and scanner rooms of a generic PET clinic. The models are defined in the Moritz geometry editing code [4] using parameters [5] representing dimensions and shield thicknesses. Moritz propagates a change in one of these values throughout the model before it prepares an input file for MCNP. The models may be used to provide facility-specific shielding designs by simply modifying the dimensions of the rooms and concrete deck thickness to match those found in individual facility design plans. In this work we use the models to calculate the attenuation provided by common thicknesses of concrete and layers of lead used for ceiling and floor shields. Both models assume a floor to ceiling distance of 4.27 m (14 feet). Fig. 1. The quiet room model. Fig. 2. Detail of wall and ceiling showing wallboard (blue), lead (black), steel decking (green), and concrete (grey). Page 2

Quiet Room The quiet room model, shown in Fig. 1, consists of a male MIRD [6] anthropomorphic model modified to a sitting posture [7] inside a room. The inside room dimensions, taken from design plans for one clinic we evaluated, are 2.43 m (8 feet) along the human model axis, 2.13 m (7 feet) across, and 3.67 m (12 feet) high. The bottom of the human trunk is centered in the room approximately 0.67 m above the floor. A transformation applied to the room results in a reclining position with the vertical axis of the human model at 30º relative to the room vertical. (The human model was not transformed to facilitate source positioning.). The walls consist of a 0.635 cm (1/4 inch) lead layer sandwiched between 1.59 cm (5/8 inch) layers of gypsum wallboard (blue in Fig. 2). Fig. 3. The scanner room model. Page 3

The source consists of two 0.511 MeV photons per 18 F disintegration with half of the 555 MBq (15 mci) activity uniformly distributed throughout each of the brain and bladder. Self-attenuation of photons in the patient is correctly modeled in this manner. The Moritz code was used to plot the positions of a few hundred source points to verify that the source was confined to the brain and bladder. Scanner Room The scanner room model, shown in Fig. 3, is derived from common plans of scanner rooms in PET/CT facilities and uses a MIRD phantom placed supine in a combination PET/CT scanner. The model uses the same floor to ceiling height as the quiet room model and interior dimensions of 2.43 m (8 feet) x 3.05 m (10 feet). The PET/CT scanner was modeled based on figures provided by GE for their current PET and Multi-Slice CT unit, but may be extrapolated to all modern PET/CT units with multi-slice CTs. The model, shown in Fig. 3, consists of iron PET and CT rings with a thickness to match the attenuation of the actual rings. The base plates for the scanner and lead isolation rings between the PET and CT gantries were also modeled. The information for the PET/CT unit was provided by GE Healthcare [8] based on their current designs. This model provides a realistic estimate of the attenuation provided by the PET/CT unit. The source term assumes a patient load of 277 MBq (7.5 mci) of 18 FDG remains after the patient has rested in the quiet room for 45 minutes and then urinated. Thirty percent of the activity is located in the brain and the rest is distributed throughout the body. The correct source distribution was verified by plotting several hundred source points with the Moritz code. MCNP Calculations Since the source terms and the geometry for these models are complex, grids of tallies were used in air above and below the rooms to determine the position and value of maximum dose rate for the floors above and below. A dose response function was applied to the volume flux (F4) tallies. The tally volumes were ¼ m thick in the vertical direction centered 0.5 m in the upper rooms and 1.7 m above the floor of the room below. Additional tallies just above the floor and below the ceiling of the quiet and scanner rooms were taken to measure the variation of dose rates within the rooms. Numerous runs were made to provide a matrix of attenuation factors for common concrete deck (5.08, 10.16, 15.24, and 20.32 cm [2, 4, 6, and 8 inches]) and lead (0, 0.15875, 0.33175, and 0.635 cm [0, 1/16, 1/8, and ¼ inches]) thicknesses used to construct these shields. Importance splitting was used through the layers of the shields. Sufficient histories were run to ensure all relative errors were < 8% and the MCNP statistical tests were passed. RESULTS The tallies show hat the largest dose rates are directly above and below the patient's head in the quiet room. For the scanner room, the high point on the floor above was also above the patient's head, but is shifted to the center of the phantom in the floor below the scanner due to added attenuation provided under the patient's head by the scanner base plate. The variation among the tallies the upper and lower rooms is small no more than a factor of 4 for both the quiet and scanner room models and less than a factor of 2 in ½ of the cases. Fig. 4 shows a 3D view of the scanner room and the tallies in the upper and Page 4

lower rooms with the tallies color coded by dose rate in arbitrary units. Fig. 5 shows 2D views of the same tallies. Fig. 4. 3D view of the tallies in the rooms above and below the scanner room colored coded by the dose rate in arbitrary units. Attenuation Factors We define the attenuation as the ratio of the dose rate from one region in a model with no floor or ceiling shielding to the rate from the same region in a shielded model. We consider two locations the tally volumes in the upper and lower rooms. The attenuations were calculated for each of these regions. Because of the variation in dose rates in these regions, the attenuation depends on the specific tally or tallies chosen for the ratio. Two methods were considered. One uses the ratio of the largest dose rates in Page 5

each location. The other takes the average dose rates in each location. Neither method gives attenuation factors that were consistently higher or lower than the other method. Tables I through IV show the attenuation factors based on the average method. The columns are labeled by the lead thickness (in cm) and the rows by the concrete thickness (in cm). Any application of the factors should take note of the variation of the dose rate with position, the height at which the tallies were taken, and the floor to ceiling height assumed in the models. Fig. 5. 2D views of the tallies in the rooms above (left) and below (right) the scanner room colored by the dose rate in arbitrary units. TABLE I. Attenuation Factors for Room Above Quiet Room. Concrete/Lead 0 0.15875 ( 1 / 16 ) 0.33175 ( 1 / 8 ) 0.635 ( 1 / 4 ) 5.08 (2 ) 1.86 2.77 3.74 6.53 10.16 (4 ) 3.89 5.58 7.48 13.01 15.24 (6 ) 8.52 11.91 15.89 27.54 20.32 (8 ) 19.23 26.39 34.99 60.49 TABLE II. Attenuation Factors for Room Below Quiet Room. Concrete/Lead 0 0.15875 ( 1 / 16 ) 0.33175 ( 1 / 8 ) 0.635 ( 1 / 4 ) 5.08 (2 ) 1.95 2.99 4.10 7.28 10.16 (4 ) 4.10 6.08 8.31 14.93 15.24 (6 ) 9.04 13.04 17.79 31.78 20.32 (8 ) 20.61 29.00 39.41 70.62 Page 6

TABLE III. Attenuation Factors for Room Above Scanner Room. Concrete/Lead 0 0.15875 ( 1 / 16 ) 0.33175 ( 1 / 8 ) 0.635 ( 1 / 4 ) 5.08 (2 ) 1.82 2.73 3.70 6.44 10.16 (4 ) 3.80 5.45 7.34 12.77 15.24 (6 ) 8.28 11.55 15.48 26.81 20.32 (8 ) 18.64 25.57 34.05 58.81 TABLE IV. Attenuation Factors for Room Below Scanner Room. Concrete/Lead 0 0.15875 ( 1 / 16 ) 0.33175 ( 1 / 8 ) 0.635 ( 1 / 4 ) 5.08 (2 ) 1.75 2.57 3.47 5.99 10.16 (4 ) 3.52 5.05 6.81 11.89 15.24 (6 ) 7.62 10.71 14.39 25.15 20.32 (8 ) 17.21 23.82 31.97 55.75 Comparison with AAPM Values The American Association of Physicists in Medicine (AAPM) has issued a guide for shielding in PET clinics [9}. The guide includes a table of attenuation factors for various thicknesses of lead, concrete, and iron. Table V compares our results for concrete only shields to the AAPM results. Our results show more attenuation by factors between 1.5 and 1.7. The AAPM results are based on Monte Carlo calculations by one of the guide s authors (D. Simpkin) who assumed a 511 kev broad beam incident on the shields. In our models, the average energy of photons entering the shields is between 300 and 250 kev because of down scattering in the patient. The lower photon energies result in more attenuation by the shields. TABLE V. Comparison of Attenuation Factors without Lead Layer. Concrete Thickness Table I AAPM Tab. I/AAPM 5.08 (2 ) 1.75 1.27 1.47 10.16 (4 ) 3.52 2.50 1.56 15.24 (6 ) 7.62 5.00 1.70 20.32 (8 ) 17.21 11.11 1.73 Scanner Attenuation Some shielding is provided to the areas above and below the scanner room by the double gantry of the scanner and the base plate. To investigate how much it contributes to the attenuation, we compared the tallies from the model of the scanner room with no floor and ceiling shielding and the scanner model present to the results from a model where the scanner gantries and base plate are replaced by air (but with the patient present). For the upper and lower rooms, the additional average attenuation factors are 1.6 and 2.3, respectively. The scanner model includes a large steel plate on the floor that contributes to the greater attenuation for the lower room. Page 7

CONCLUSIONS We have run a number of Monte Carlo calculations for various combinations of lead and concrete shielding in the floor and ceiling of typical quiet and scanner rooms in a PET clinic. We have developed models in which the room dimensions and shielding thicknesses can easily be adjusted to match a specific configuration. Down scattering of the photon energy within the patient results in more attenuation through the shields than when assuming a 511 kev photon energy incident on the shield. The PET and CT rings contribute a factor of 2 to the total attenuation. REFERENCES 1. A. CHILTON, R. FAW, J. SHULTIS, Principles of Radiation Shielding, Prentice Hall, Englewood Cliffs, USA (1984). 2. D. K. TRUBY, et al., Gamma-ray attenuation coefficients and buildup factors for engineering materials, American Nuclear Society, LaGrange Park, IL, ANSI/ANS-6.4.3 1991 (1992). 3. J. F. BRIESMEISTER, Editor, MCNP A General Monte Carlo N-Particle Transport Code, Los Alamos National Laboratory Report LA-13709-M, (2000). 4. K. A. VAN RIPER, Interactive 3D Display of MCNP Geometry Models, Proceedings of the ANS International Meeting on Mathematical Methods for Nuclear Applications, Salt Lake, UT, Sept. 10 13, (2001). 5. K. A. VAN RIPER, Parameterized Combinatorial Geometry Modeling In Moritz, Proceedings of Mathematics and Computation, Supercomputing, Reactor Physics and Nuclear and Biological Applications, Palais des Papes, Avignon, France, September 12-15, 2005, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2005). 6. M. CRISTY, K. F. ECKERMAN, Specific Absorbed Fractions of Energy at Various Ages from Internal Photon Sources. I. Methods, Oak Ridge National Laboratory Report ORNL/TM 8381/VI (1987). 7. R. H. OLSHER, K. A. VAN RIPER, Application of a Sitting MIRD Phantom for Effective Dose Calculations, Proceedings of the ICRS 10/RPS 2004 Conference, 9 14 May 2004, Funchal, Portugal (Radiation Protection Dosimetry, in press) (2005). 8. J. WILLIAMS, Ph.D., GE Healthcare, private communication (2005). 9. M. T. MADSEN, J. A. ANDERSON, J. R. HALAMA, J. KLECK, D. J. SIMPKIN, J. R. VOTAW, R. E. WENDT III, L. E. WILLIAMS, and M. V. YESTER, Draft AAPM Guide on PET and PET/CT Shielding Requirements, AAPM (2004). Page 8