Peter J. Biggs Ph.D. Massachusetts General Hospital Harvard Medical School Boston, MA 02114

Similar documents
Peter J. Biggs Ph.D., Massachusetts General lhospital, Harvard Medical School, Boston, MA 02114

Radiation Shielding for Megavoltage Photon Therapy Machines

Units S H I E L D I N H = D Q. H: Dose equivalent (Sv) D: Dose (Gy) Q: Quality Factor. 1Sv = 1 J/Kg. 1Gy = 1 J/Kg

DESIGN OF A TREATMENT ROOM FOR AN 18-MV LINAC

Progress in Nuclear Science and Technology, Volume 6,

A study on the cost of concrete shielding in a standard radiotherapy facility room

Chapter 4. QUANTIFYING THE HAZARD II: DATA & ANALYSIS. The dose equivalents for spheres in air with 10 cm radius centred at a point in the

Calculations of Photoneutrons from Varian Clinac Accelerators and Their Transmissions in Materials*

Neutron source strength measurements for Varian, Siemens, Elekta, and General Electric linear accelerators

Radiation shielding for gamma stereotactic radiosurgery units

Proposed Room Requirements for CT System

Radiation Shielding. PTCOG 57 Cincinnati, USA (2018)

Proposed Room Requirements for CT System

Radiation protection issues in proton therapy

ASSESSING LEAKAGE WORKLOADS OF MEDICAL LINEAR ACCELERATORS FOR IMRT AND TBI TECHNIQUES

The neutron dose equivalent evaluation and shielding at the maze entrance of a Varian Clinac 23EX treatment room

Secondary Radiation and Shielding Design for Particle Therapy Facilities

Shielding Design Considerations for Proton Therapy Facilities

COMPARISON OF COMPUTER CODES APPLICABILITY IN SHIELDING DESIGN FOR HADRON THERAPY FACILITIES *

Radiation safety of the Danish Center for Proton Therapy (DCPT) Lars Hjorth Præstegaard Dept. of Medical Physics, Aarhus University Hospital

Shielding Calculation Techniques

Estimation of neutron and gamma radiation doses inside the concrete shield wall for 10 and 15 MV medical linear accelerators

Applications of MCBEND

Radiation Protection Dosimetry (2007), Vol. 126, No. 1 4, pp Advance Access publication 11 May 2007

Estimate of Photonuclear Reaction in a Medical Linear Accelerator Using a Water-Equivalent Phantom

Neutron Fluence and Energy Spectra Around the Varian Clinac 2lOOC/23OOC Medical Accelerator

Shielding verification and neutron dose evaluation of the Mevion S250 proton therapy unit

CHARACTERISTICS OF DEGRADED ELECTRON BEAMS PRODUCED BY NOVAC7 IORT ACCELERATOR

Michael Dunn Nuclear Data Group Leader Nuclear Science & Technology Division Medical Physics Working Group Meeting October 26, 2005

Accelerator Facility Shielding Design

Radiation Shielding of a 230 MeV Proton Cyclotron For Cancer Therapy

I. INTRODUCTION EXPERIMENTAL

Secondary Particles Produced by Hadron Therapy

NEUTRON RADIATION FROM MEDICAL ELECTRON ACCELERATORS*

Safety Reports Series No.47. Radiation Protection in the Design of Radiotherapy Facilities

Monte Carlo Calculations Using MCNP4B for an Optimal Shielding Design. of a 14-MeV Neutron Source * James C. Liu and Tony T. Ng

Hospital Cyclotrons: Radiation Safety Aspects. Matthew Griffiths

Radiation Protection At Synchrotron Radiation Facilities

Derivation of factors for estimating the scatter of diagnostic x-rays from walls and ceiling slabs

Shielding Design for the Imaging and Medical Beamline at the Australian Synchrotron

QUESTION 8 GIVEN. Electron beam kinetic energy = 20 MeV Peak current = 1 A Beam pulse length = 1 microsecond Beam pulse frequency = 10 Hz

Beam diagnostics: Alignment of the beam to prevent for activation. Accelerator physics: using these sensitive particle detectors.

Photon-beams monitor-unit calculations

Neutron Shielding Properties Of Concrete With Boron And Boron Containing Mineral

Work supported by Department of Energy contract DE-AC03-76SF00515

X- AND GAMMA RAY EXPOSURE CALCULATION OF SHIELDING

Physics of Radiotherapy. Lecture II: Interaction of Ionizing Radiation With Matter

Neutronic Design on a Small Accelerator based 7 Li (p, n) Neutron Source for Neutron Scattering Experiments

International Journal of Scientific & Engineering Research, Volume 7, Issue 2, February-2016 ISSN

Physics of particles. H. Paganetti PhD Massachusetts General Hospital & Harvard Medical School

ESTIMATION OF 90 SCATTERING COEFFICIENT IN THE SHIELDING CALCULATION OF DIAGNOSTIC X-RAY EQUIPMENT

Secondary Neutron Dose Measurement for Proton Line Scanning Therapy

Radiation Shielding of Extraction Absorbers for a Fermilab Photoinjector

Measurement of induced radioactivity in air and water for medical accelerators

Volume 1 No. 4, October 2011 ISSN International Journal of Science and Technology IJST Journal. All rights reserved

CALCULATION OF SHIELDING AND RADIATION DOSES FOR PET/CT NUCLEAR MEDICINE FACILITY

Radiation exposure of personnel during IORT: radiation protection aspects.

NEUTRON EXPOSURE FROM ELECTROM LINEAR ACCELERATORS AND A PROTON ACCELERATOR: MEASUREMENTS AND SIMULATIONS

PECULIARITIES OF FORMING THE RADIATION SITUATION AT AN AREA OF NSC KIPT ACCELERATORS LOCATION

Laboratoire National Henri Becquerel (CEA/LIST/LNHB), France (2) ENEA-Radiation Protection Institute, Bologna, Italy (3)

Radiation Quantities and Units

In today s lecture, we want to see what happens when we hit the target.

Neutron Skyshine Calculations with the Integral Line-Beam Method

City University of Hong Kong

Doing the right thing, doing it right

A Measuring System with Recombination Chamber for Photoneutron Dosimetry at Medical Linear Accelerators

Particle Beam Technology and Delivery

Comparison with simulations to experimental data for photoneutron reactions using SPring-8 Injector

Introduction to neutron sources

(Tandem Collimators for the Tangential GammaRay Spectrometer - KM6T-TC)

NEUTRON SPECTROMETRY WITH BUBBLE DETECTORS

Prompt Radiation Fields at Accelerators

Calibration of the GNU and HSREM neutron survey instruments

Radiation Safety Considerations for the TPS Accelerators

(10%) (c) What other peaks can appear in the pulse-height spectrum if the detector were not small? Give a sketch and explain briefly.

Shielding design for a laser-accelerated proton therapy system

1. What would be the dose rate of two curies of 60Co with combined energies of 2500 kev given off 100% of the time?

Monte Carlo simulation of the photoneutron field in linac radiotherapy treatments with different collimation systems

Prompt gamma measurements for the verification of dose deposition in proton therapy. Contents. Two Proton Beam Facilities for Therapy and Research

Objectives. Tolerance Limits and Action Levels for IMRT QA. The Overall Process of IMRT Planning and Delivery. Chain of IMRT Process

1.1.4 What are the differences between the Varian 600C and the Siemens KD2?

Reactor radiation skyshine calculations with TRIPOLI-4 code for Baikal-1 experiments

Determination of Ambient Dose Equivalent at INFLPR 7 MeV Linear Accelerator

Measurements with the new PHE Neutron Survey Instrument

Bulk shielding design for the MAX IV facility

Overview and Status of the Austrian Particle Therapy Facility MedAustron. Peter Urschütz

Low energy linear acceleratordriven subcritical assembly

Journal of Radiation Protection and Research

Implementation of the IAEA-AAPM Code of Practice for the dosimetry of small static fields used in external beam radiotherapy

Effect of Cosmic-ray Shielding in Passive Neutron Coincidence Counting

Evaluate Shielding Design of the Brachytherapy Unit by Using Monte Carlo Simulation Code

USE OF DLP FOR ESTABLISHING THE SHIELDING OF MULTI- DETECTOR COMPUTED TOMOGRAPHY ROOMS

Georgia Institute of Technology. Radiation Detection & Protection (Day 3)

Outline. Absorbed Dose in Radioactive Media. Introduction. Radiation equilibrium. Charged-particle equilibrium

BULK-I. Radiation Shielding Tool for Proton Accelerator Facilities

Optimization studies of photo-neutron production in high-z metallic targets using high energy electron beam for ADS and transmutation

IRPA 2 nd European Congress, Paris, May Session: RP in Medicine - Workers. Activation Products In Medical Linear Accelerators

Neutron collimator design of neutron radiography based on the BNCT facility arxiv: v1 [physics.ins-det] 3 May 2013

Calculations of Neutron Yield and Gamma Rays Intensity by GEANT4

NEUTRONS FROM MEDICAL EJ,ECTRON ACCELERATORS '

Transcription:

Basic Design Principles and Latest Recommended Data, Neutrons and Structural Concerns Winter on Cape Cod Peter J. Biggs Ph.D. Massachusetts General Hospital Harvard Medical School Boston, MA 02114

Outline 1 Basic principles and definitions 2 Calculations 3 Workload, use and occupancy factors 4 Shielding materials 5 Tenth value layers (TVL) 6 Room layout & features, construction details 7 Neutrons & laminated barriers 8 Mazes and doors low energies 9 Ducts 10 Skyshine 11 Examples of calculations

Definitions P: Weekly design dose limit (Sv/wk) d: Distance from target to measurement point W: Workload (Gy/wk) U: Use Factor T: Occupancy factor a: Scatter fraction; (θ, E) d sec : Distance from scatterer to measurement point d sca : Distance from target to scatterer d l: Distance from the target to measurement point F: Area of the beam in the plane of the scatterer (cm 2 ) B: Barrier transmission factor

Basic Principles The purpose of radiation shielding is to reduce the effective equivalent dose from a linear accelerator to a point outside the room to a level that is determined by individual states. In the US: Public or uncontrolled area - 0.02 msv per week - 0.02 msv in any one hour Controlled area - 1 msv/wk k(in practice 0.1 msv/wk)

Basic Shielding Equations: NCRP #151 Methodology Primary: Scatter: Leakage: B p B B s = l = 2 Pd WUT = I in I out P awt d 2 sec 1000* Pd WT d 2 l B x is the barrier transmission factor 2 sca 400 F B x = I I out in The factor 1000 is due to the 0.1% leakage requirement

Number of Tenth Value Layers (Primary & Leakage) The number, n, of tenth value layers (TVLs) required to reduce the dose to the value P is given by: n = log10 ( B) and the thickness for primary pimay and leakage radiation is given by: T = TVL + ( n 11 ) TVL p, l 1 p, l e

Number of Tenth Value Layers for Scattered dr Radiation Scatter TVLs are also dependent on scattering angle and there is no shoulder to the attenuation curve, so: T = n s s TVL s ( ) ϑ

Tenth Value Layers Primary TVLs (cm)* Energy (MV) 4 6 10 15 18 20 Concrete 35(30) 37(33) 41(37) 44(41) 45(53) 46(44) Steel 9.9 10 11 11 11 11 Lead 5.7 5.7 5.7 5.7 5.7 5.7 Leakage TVLs (90 )* Energy (MV) 4 6 10 15 18 20 Concrete 33(28) 34(29) 35(31) 36(33) 36(34) 36(34) Steel** 8 8.5 8.7 8.7 Lead** 4.5 4.6 4.7 4.7 * First term is 1st TVL and term in brackets is for all other TVLs. Data from NCRP #151 ** Data from McGinley

Obliquity Factor -I P

Obliquity yf factor - II Rule of thumb: - for θ >45 - and B<0.001 - Add ~2 HVL for low-energy photons - And ~1 HVL for high-energy photons

Obliquity factor - III - Experimental data by Kirn et. al. (1954) - verified by Biggs (1996) using the Monte Carlo for a large range of clinical energies - Monte Carlo results are in good agreement with previous recommendations at low energy, but more detail is given - Caveat: beware of applying obliquity factor corrections for large angles!!

Conventional Workload - Number of patients treated per week multiplied by the average dose (MU)/patient at isocenter - Include calibration, service - Recommended dvalues: 500 or 1000 Gywk -1 - Survey shows values from 250 to 450 Gy wk -1 - Balance between highh energy and low energy use - Electron-onlymachines? l

Workload TBI - Workload for TBI>>workload for conventional therapy due to the extended distances: W = D TBI 2 TBI - Leakage workload is also higher, but patient- and wall-scattered workload is not. TBI - Radiation is usually directed at one barrier - Treatment time is much longer than for conventional R x d

Workload IMRT (1) - Workload for IMRT is complicated by the fact that the #MU is much larger than for conventional therapy. The IMRT factor is defined as: C = I MU MU IMRT conv and varies bt between 5 and 10 or more - Thisonly affects leakage radiation, not primary or scattered radiation.

Workload IMRT (2) - Serial tomotherapy has the highest C values, >10. Helical tomotherapy is much less. - Linac based IMRT at MGH gives an average C value of about 6 - Anecdote: For a 6/18 MV machine the energy use prior to IMRT was 20%/80% (MU). With 28% IMRT patient load, the use was 70%/30%

Copnventional Use Factors 90 interval: - 0 (31%); 90 and 270 (21.3 % each); 180 (26.3 %) - previously, all barriers were assigned 0.25 45 interval: - 0 (25.6 %); 45 and 315 (5.8 % each); 90 and 270 (15.9 % each); 135 and 225 (4 % each); 180 (23 %)

Use Factor for Conventional Therapy

Use Factors Special Cases - A significant TBI load will require one wall to have an increased use factor - IMRT may also require a change in the values assigned to the use factor for the primary barrier

Use Factor Dedicated Machines - One should pay py attention to the types of procedures to be carried out in new room - A room dedicated to SRS and SRT requires different considerations from a general therapy room - A room dedicated to breast treatments will use specific c ranges of primary ybeam angles ges

Example: Use Factor for Stereotactic Radiosurgery

Typical Angles Used in Breast Therapy

Workload - Summation (1) Primary: WU p = ( W U W U + W U + W U +...) conv = WU wall scat conv Patient scatter: W + TBI TBI IMRT IMRT QA QA scatter, iso = conv IMRT QA + = ( W + W + W +...) (note that the TBI contribution is calculated separately source not at isocenter)

Workload - Summation (2) Leakage: W L = conv TBI I IMRT QA QA ( W + W + C W + C W +...)

Occupancy Factors - Occupancy factors down to 1/40 are now allowed (previously 1/16) - Lowest occupancy allows 0.02 msv in any one hour and 0,02 msv/wk - Of note, the occupancy outside a treatment room door is 1/8, instead of the conventional unity for a controlled area

If the Dose Rate is High, Should I Freak ko Out? Note that at the weekly limit of 0.0202 msv/wk, the dose rate outside a barrier could be quite high. e.g. for a primary barrier, an occupancy of 1/40 would imply a weekly dose of 3.2 msv and for a beam on time of 2.8 hr, h,the instantaneous dose rate is 1.1 msv/hr!! Moral: Choose factors wisely!!

TADR - TADR is the dose rate averaged over a week, a day or an hour: Daily beam on time TADR = IDR * Length of working week - The TADR (1 hour) or dose limit it in any one hour is: WU R * h h = IDR DR 0

TADR - NCRP recommends use of 1 hr as the minimum period for measuring IDR - Not an issue for occupancy of f1 - Only an issue when treatments t t are not given uniformly throughout the 40 hr week

Structural Details

Shielding Materials Materials Density Comments (g/cm 2 ) Concrete 235;385 2.35; 3.85 High density concrete is very expensive Concrete blocks 2.35; 3.85; 4.62 Lack structural integrity of concrete Lead 11.35 Great for photons; bad for neutrons. Needs structural support Steel 7.87 Not as efficient as lead for photons Earth 1.5 Cheap! Therefore, build underground Brick 1.65-2.05 Polyethylene; l ~1.0 10 Used to shield against neutrons in doors, ducts etc. borated polyethylene ** Note: For h.d. concrete blocks, mortar of equiv. density should be used.

Possible Beam Orientations for Room with Maze

Width of Primary Barrier -I Allow 0.3 m (1 foot) either side of the primary beam w = 0.4 2 d + 0.6 - allow for fixed primary collimator diameter (0.5 m) w = 0.5 d + 0.6 - Allow for greater width at the top of the wall w = 0.5d cos( ϑ) + 0.6 where d is the distance between the source and a point in the barrier and θ is the angle subtended by the central axis with the horizontal

B A C D

Contouring the Primary Shielding in the Ceiling

Physics Conduit

Lintel over Maze Wall

Lead-Only Room: Groundshine

Construction Techniques for Upgrading Linear Accelerator Rooms

A B C

Neutrons

Neutron Photoproduction Cross-Sections on Pb Cross: (γ,xn) Plus: (γ,1n)

Neutron Yields vs Energy and Z kw -1 ) ns sec -1 k Y (neutron Electron energy (MeV)

Photoneutron Energy Spectrum Re elative nu umber of neutron ns per ½ MeV Neutron energy, E n (MeV)

Shielding Thickness Required to Halve the Average Neutron Energy Half Ene ergy Lay yer (cm) most probable energy E(MeV)

Neutron Spectra from Medical Linacs - Neutrons produced d in the head of the linac are first moderated by the x-ray shielding - Neutrons are further moderated by scattering off the concrete walls of the therapy room

Neutron Spectra from Medical Linacs - The total neutron fluence therefore consists of direct ( fast ) fast) neutrons, scattered neutrons and thermal neutrons Φtotal = Φdir + Φsca + Φth - Fast neutrons obey the inverse law, but scattered and thermal neutrons are isotropically distributed; hence the neutron fluence drops off less fast than inverse square

Facure et. al.: 10 MV

Production of Neutrons by Primary Beam in a Laminated ated Barrier Incoming photon t t m t2 1 neutron Capture γ ray For high energy x-ray beams, the x-ray dose is enhanced by a factor of 2.7 to be conservative

Production of Neutrons by Primary Beam in a Laminated Barrier H DRF 0 max 1 x 2 = *10 *10 tm + t2 + 0.3 2 t / TVL t / TVL n where e H is the neutron dose equiv.(µsv wk -1 ) D 0 is the x-ray dose at isocenter (cgy wk -1 ) R is the neutron prod n. rate (µsv cgy -1 m -2 ) and F max is the max. beam area at isocenter (m 2 )

Total Dose Behind a Laminated Barrier (1) H tot = H n + H ph = H n1 +H n2 +2.7*H tr The 2.7 factor accounts for the production of capture gamma rays by the photo-produced neutrons. This is an conservative, empirical figure givenbymcginley.h n1 and H n2 represent direct neutrons and photo-produced neutrons respectively..

Total Dose Behind a Laminated Barrier (2) - Photoneutron production data is given for 18 MV for lead and steel and 15 MV for lead - The capture gamma ray enhancement factor was derived empirically for 18 MV on lead. No data for lower energies. However, stated to be good for all energies.

Mazes & Doors Low Energies

Contributions to the Dose at the End of a Maze 1. Primary scatter 2. Patient scatter 3. Leakage along maze 4. Leakage through maze wall

Secondary Radiation at Door due to Wall Scatter of Primary Beam S s = W U B α Aα A 1 1 2 ( d d d ) 2 i r1 r 2 2

Secondary Radiation at Door due to Patient Scatter S p = aw U B ( F / 400) ( d d d ) 2 sca sec s α A 1 1

Secondary Radiation at Door due to Leakage L = L 0 W U B α1 ( d d ) 2 d s d 1 A 1

Secondary Radiation at Door due to Leakage Through Maze Wall L d = L0 W U B10 ( d ) 2 dl d ( t / TVL) Note that this radiation is more energetic than radiation along the maze

Mazes The total dose when the beam is directed towards the wall is the sum of the four components: D t = f*s s + S p + L + L d Note that the S s component is reduced by a factor f to account for patient attenuation. Assuming use factors of ¼ for the 4 gantry angles, McGinley gives an average total dose of: D T = 264D 2.64 t

Comments on NCRP Formalism Several authors have noted that the foregoing formalism underestimates the true dose at the door. This problem has been partly rectified by redefining the geometry of the scatter and assuming lower energies for the scattered radiation which, in turn, increase the scatter coefficients Numark &Kase;Al-Affan; McGinley & James

Reflection Coefficients for Monoenergetic X- Rays on Concrete NCRP #51

Maze Special Case

Ducts

Principle of Minimizing Loss of Shielding due to Presence of Duct Outside room Shielding wall Inside room Duct Leakage or scatter radiation Plane of gantry rotation Isocenter

A Ideal arrangement Additional shielding Ceiling Duct B C Ceiling duct McGinley s recommendation Isocenter

Skyshine Photons Neutrons

Skyshine: Measurements etc. - McGinley has shown that, for x-rays, the calculations are an underestimate at short distances (factor of 4) and an overestimate at long distances (factor of 5). - McGinley has shown that, for neutrons, the calculations are a gross underestimate at all distances. - Zavgorodni has presented a formalism for side scattered radiation for radiation penetrating thin barriers.

Practical Problems with NCRP #151 - Use of dual energy workloads, dual TVLs, laminated barriers creates confusion in the calculations - Use of the obliquity factor for scattered radiation (not leakage) is unclear two citations in the text give opposing views. Experts are divided on this issue, but it can be used, but with caution!

Handling Dual Energies in NCRP #151 - A dual energy machine having say 6 and 18 MV will have different workloads for (P,S) and L radiation for 3D CRT and IMRT (6 MV assumed). In general, a conservative approach assumes: - use of the maximum energy in a full 3D CRT environment for the primary and - the IMRT workload for the secondary barrier orthogonal to the plane of gantry rotation

Handling Dual Energies in NCRP #151 - in the secondary barrier region adjacent to the primary, the dominance of 3D CRT vs 3D CRT+IMRT depends on the explicit workloads - one should be careful to consider 20 scatter radiation for laminated barriers since the primary beam angle is only 14 (exclusive of extra 2 barrier width recommended by NCRP)

Handling Dual Energies in NCRP #151 - In that same angular zone in the secondary barrier, there will be be TVLs for two energies for scattered and leakage radiation - Rule is that: - scatter should be added to scatter to determine if an extra HVL is required - leakage should be added to leakage to determine if an extra HVL is required - Scatter is then compared to leakage to determine if an additional HVL is required

Handling Dual Energies in NCRP #151 - To recap: (i) S1 + S2 S (ii) L1 + L2 L (iii) S +L Total secondary thickness - Question remains as to which of the TVLs to use; higher or lower for L (ii); S or L for (iii) - However, because of these multiple additions, this may overestimate t the required shielding

Example: Barrier Adjacent to Primary 3D CRT + IMRT Energy Radiation Workload (Gy) Thickness Comment 15 MV Leakage 225 36.1 6 MV Leakage 1775 42.9 Combined 6 +15 MV Leakage 2000 47.2 Use TVL 1 for 15 MV 15 MV Scatter 225 45.6 6MV Scatter 525 42.0 Combined 6 +15 MV Scatter 750 49.3 Use TVL sca for 15 MV Combined 6 + 15 MV Leakage + scatter 53.3 Use TVL 1 for 15 MV N o t e Note: and only one shielding material is used!

Example: Barrier Orthogonal to Gantry Rotation Plane 3D CRT + IMRT Energy Radiation Workload (Gy) Thickness Comment 15 MV Leakage 225 37.3 6 MV Leakage 1775 43.9 Combined 6 +15 MV Leakage 2000 47.8 Use TVL 1 for 15 MV 15 MV Scatter 225 19.8 6 MV Scatter 525 22.6 Combined 6 +15 MV Scatter 750 24.7 Use TVL sca for 15 MV Combined 6 + 15 MV Leakage + scatter 47.8 Use TVL 1 for 15 MV

TVLs for Laminated Secondary Barriers

Summary - The basic formalism of the calculations is well explained - NCRP #151 addresses the shielding issues raised by modern radiotherapy accelerators and techniques hi and allows the designer to perform the necessary calculations (but see second talk) - Data is up-to date - A complete set of examples are given for every type of barrier plus robotic arm (~30% of book

Thank Thank you foryou for your attention! your attention! Muito obrigado pela sua atenção