BJR
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

British Journal of Radiology (2003) 76, 638-647
© 2003 British Institute of Radiology
doi: 10.1259/bjr/21647806

This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Schultz, F W
Right arrow Articles by Zoetelief, J
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Schultz, F W
Right arrow Articles by Zoetelief, J

Full Paper

Monte Carlo calculations for assessment of radiation dose to patients with congenital heart defects and to staff during cardiac catheterizations

F W Schultz, PhD 1 J Geleijns, PhD 2 F M Spoelstra, PhD 2 and J Zoetelief, PhD 1

1 Interfaculty Reactor Institute, Delft University of Technology, Mekelweg 15, 2629 JB Delft and 2 Department of Radiology, Leiden University Medical Center, PO Box 9600, 2300 RC Leiden, The Netherlands

Correspondence: F W Schultz, Interfaculty Reactor Institute (IRI/ST-MSF), Mekelweg 15, 2629 JB Delft, The Netherlands


    Abstract
 Top
 Abstract
 Introduction
 Methods
 Results and discussion
 Conclusions
 References
 
Effective dose is an important quantity in relation to assessment of radiation risk. Organ and effective doses to paediatric patients undergoing diagnostic and therapeutic heart catheterization procedures can be assessed by combining relatively simple measurements, e.g. of dose–area product (DAP), and calculated dose conversion factors (DCF). This also holds for the radiation dose to the hospital staff, e.g. the cardiologist. Monte Carlo (MC) simulation of radiation transport in mathematical anthropomorphic phantoms is used to obtain the DCFs, which strongly depend on beam quality and geometrical parameters. The performance of a dedicated fast MC code (PCXMC) for patient dosimetry is compared with that of a more elaborate general purpose MC code (MCNP). Resulting organ doses sometimes may differ considerably, partly due to phantom differences. While MCNP uses separate male and female mathematical phantoms, PCXMC uses a hermaphrodite. However, both codes yield effective doses that agree rather well, so PCXMC can be used for convenience. The MCNP code is used to calculate the effective dose to the cardiologist exposed to radiation scattered from the patient. Without protective clothing, effective dose per procedure to the cardiologist is at least two orders of magnitude lower than that to the patient. The effectiveness of various types and thickness of protective clothing has been evaluated for one view of one cardiac catheterization. The results of the calculations do not contradict experimental studies from the literature. MC simulation may serve as a useful tool to improve the accuracy of estimating occupational effective dose from personal dose monitors.


    Introduction
 Top
 Abstract
 Introduction
 Methods
 Results and discussion
 Conclusions
 References
 
Interventional radiology and cardiac catheterizations are developing into preferred alternatives to surgical procedures that tend to lead to hospitalization of patients for longer time periods. Inflicting less physical damage due to minimal invasive procedures not only means a decreased burden to the patient, but society in general also profits when the patient returns sooner to normal life. The medical application of radiation, however, requires proper justification and optimization. The consumption of fluoroscopy time for guiding a procedure may be considerable. This has led to radiation injuries in the past [1]. Accurate estimates of radiation doses are necessary to assess the risk to the patient. Also medical staff are at risk, as continuous presence near the patient is inevitable and results in prolonged exposure to (scattered) radiation.

X-ray guided cardiac catheterizations are a form of radiology that is applied to examine and treat paediatric patients with congenital heart defects relatively often. High levels of radiation exposure may occur due to using biplane X-ray equipment and prolonged fluoroscopy time and cine imaging to monitor and record access of the flexible wire to the diseased region after entering via a peripheral blood vessel. As children are relatively more radiosensitive than adults [2] and have longer mean lifetime expectancy, it is very important to optimize procedures and keep the dose to the paediatric patient as low as possible. On the other hand, the occupational dose to the medical staff tends to be high in comparison with other radiological procedures. Furthermore, it may build up quickly because the same personnel will be involved regularly in similar catheterizations. The risk to the cardiologist is most probably highest, as this person remains longest with and nearest to the patient during the procedure.

Effective dose, a weighted sum of organ doses [2], is considered the best dosimetric quantity for estimating the risk of exposure to ionizing radiation. To compute effective dose it is necessary to know the organ doses with good accuracy. Direct measurements are very difficult if possible at all. The usual approach is to take a dosimetric quantity for which the value is simpler to establish, e.g. air kerma free in air (Ka), entrance skin dose (ESD) with or without backscatter, or dose–area product (DAP). The organ doses are obtained subsequently by multiplying by an appropriate dose conversion factor (DCF). DCFs depend on the exposure conditions and are generally determined through Monte Carlo (MC) simulation of radiation transport.

As a result of elaborate and careful development of software and the advance of cheaper and faster computers, MC simulation has become a well-accepted technique for solving radiation transport problems in a multitude of applications [3] including medical physics. One of the general-purpose codes is MCNP, which originates from the Los Alamos National Laboratory [4] and has been made widely available. Owing to its great flexibility with respect to the type of application, a relatively large effort from the user is required to construct the proper input and retrieve the desired output information.

In contrast, a commercially available MC code like PCXMC [5] is designed to serve just one type of application, i.e. the calculation of patient doses (effective and organ) in diagnostic radiology. Therefore, the variation in the configuration of input and output variables is limited, which enables a user-friendly interface with an interactive menu structure. Furthermore, by allowing a simplified treatment of the physics, the required processing time for the actual MC simulations can be kept very short.

At the Leiden University Medical Center (LUMC) radiation exposure to patients and staff during paediatric heart catheterizations has been investigated. Patient doses due to various procedures were determined by measuring DAP and using DCFs obtained with PCXMC. PCXMC was chosen because a considerable number of DCF sets had to be evaluated, which would have been too time-consuming with MCNP.

The main results of the investigation are presented elsewhere [6]. The present paper is meant to demonstrate that the choice of PCXMC for calculating the DCFs is justified. To that purpose the more extensive code, MCNP is applied to a small sample of cases. The results are compared with those from PCXMC.

Furthermore, the present paper looks deeper into the radiation burden of the cardiologist, being the staff member who collects the largest dose. In [6], doses to staff are expressed as readings from an electronic personal dosimeter (EPD). For one procedure, DCFs for the scattered radiation are calculated with MCNP, as PCXMC is not suitable for this type of problem. The influence of several shielding measures (apron, collar) is examined. Effective dose is compared with the EPD readings and with calculated values of the personal dose equivalent for penetrating radiation, HP(10) [7]. It is important to know to what extent these quantities agree, because the measured value is assumed to represent the personal dose equivalent, which, as an operational quantity, is in turn supposed to reflect effective dose.


    Methods
 Top
 Abstract
 Introduction
 Methods
 Results and discussion
 Conclusions
 References
 
The Monte Carlo codes
MCNP is a well-known general purpose Monte Carlo code for transport of neutrons, photons and electrons developed at the Los Alamos National Laboratory [4]. The user can apply up to second order surfaces (boxes, ellipsoids, cones, etc.) and fourth order torii to build a three-dimensional (3D) geometry, which can be filled with materials of arbitrary composition and density. Point, surface or volume sources of radiation can be defined, from which the mentioned particles are emitted with user-specified probability distributions for energy and direction. The code then simulates the particle tracks and interactions with the materials, according to probability density distributions implied by particle and material properties. Taking comprehensive account of the underlying physics of radiation-matter interaction, it creates secondary particles (which are also transported) and keeps a record of quantities like particle fluence, energy deposition and dose. The user indicates at what points, surfaces or volumes these quantities are reported.

For this paper MCNP version 4C is used, which has been implemented on a Compaq XP900 Alpha workstation (Compaq Computer B.V., Gouda, The Netherlands. Unix operating system). Without attempting optimization, i.e. application of additional variance reduction techniques, it typically takes some 6 h of computer processing time (20 million starting particles) to yield less than 0.5% relative statistical uncertainty in the calculated effective dose conversion factor for patients. To calculate the effective DCF for the cardiologist exposed to radiation scattered from the patient, about four times as much processing time is required to yield approximately 4% uncertainty.

PCXMC, developed at the Medical Radiation Laboratory of the Finnish Radiation and Nuclear Safety Authority [5, 8], is a special purpose code for calculating patient dose in diagnostic radiology from either measured dose–area product or entrance air kerma (without backscatter). Operating under Windows (95/98/NT) on a personal computer, the code is fast (minutes) but transports only photons in a limited energy range (up to 150 keV). Focus-to-skin distance (FSD), degrees rotation about the patient's longitudinal axis and in craniocaudal direction, skin entry point of the central beam axis and field size at this position determine shape, position and orientation of the beam. A graphical interface facilitates the beam positioning. A hermaphrodite mathematical phantom represents the patient. In addition to adults, one of five age classes can be chosen for paediatric examinations. To some extent a standard phantom can be adapted to the actual length and weight of the patient. Doses to the eye lenses are not computed. This may be a disadvantage when evaluating examinations with beams in the head region.

Description of the geometry for MCNP
All components of the geometry are placed within an air-filled box. In the practical situation the patient is placed in supine position on the centreline of a table (not actually modelled for the radiation transport simulations) of adjustable height (Figure 1Go). The width of the table is 55 cm, including the mattress. The cardiologist leans against the table at the right hand side of the patient, at the level of the groin (or more toward the feet for shorter patients). Thus, the patient's heart can be imaged either from left to right (LLAT) using the lateral tube – to the left of the cardiologist – or from back to front (PA) with the undercouch tube of the biplane X-ray system (Philips Diagnost, Philips Medical Systems, Best, The Netherlands). Imaging is not always performed using exact PA or LLAT projections (as defined in [9]) but deviations of some 10–20 degrees may occur. The undercouch tube can be rotated (if to the right side of the patient this yields a LAO, left anterior oblique view) and/or tilted toward the head (positive angle) or the feet. The lateral tube can also be rotated toward a LAO view. Only the focus of the X-ray tube is modelled, i.e. as a standard MCNP point source emitting a cone of photons (of which the cone axis and top angle can be adjusted). Photon stopping material is used as a collimator model to reduce the circular cross section of the cone to a rectangle. Actual FSD, entrance point of the central beam axis on the patient and field size in the entrance plane determine the pyramid shape of the simulated X-ray beam. In reality a calibrated Diamentor M4 DAP meter (PTW, Freiburg, Germany) mounted on the X-ray machine measures the characteristic physical quantity for normalization in the dose calculations. In the simulations the DAP is calculated as the dose in air averaged over the beam area at, arbitrarily, 20 cm from the focus, multiplied by the beam area at this distance. The mathematical anthropomorphic phantom ADAM [10] represents the cardiologist (male adult). It contains all organs that are required for the calculation of effective dose. One of the paediatric MIRD phantoms [11] is chosen to represent the patient. Based on gender, body weight and age of the patient, a male or female phantom from the age class newborn, 1, 5, 10 or 15 years is selected. Composition and density of the tissues are the same as for ADAM. For lateral views the arms are removed. The phantoms, largely described by first and second order equations, are adapted to the input format of the MCNP code.



View larger version (16K):
[in this window]
[in a new window]
 
Figure 1. Geometry of the exposure conditions, (a) side and (b) top view. A and A': focus of lateral and undercouch tubes, respectively, the latter being rotated in this case from pure posteroanterior to 21 degrees LAO. The former presents a left lateral (LLAT) view. B and B' dose–area product meter positions at 20 cm from the focus. C: lead equivalent glass (50 x 40 x 0.1 cm3) with curtain attached below (50 x 35 x 0.05 cm3), at 8 cm from the lateral focus. D: cardiologist (Adam phantom, 70.6 kg, length 170 cm). The arrow indicates the position of the personal dosimeter (EPD). E: patient (here male, 5-year-old MIRD phantom, 18.8 kg, length 109 cm) with circumference of the heart shown. For lateral exposure the arms are removed (broken lines). The whole is enclosed in an air-filled box.

 
In practice a combination of lead equivalent glass and curtain is placed against the lateral tube to shield the cardiologist. The top of the sheet of glass (50 x 40 x 0.1 cm3) is at the level of the cardiologist's diaphragm. The curtain (50 x 35 x 0.05 cm3) is attached to the bottom of the glass sheet and reaches approximately the cardiologist's knees. The distance to the focus is 8 cm. This shield was added to the simulated geometry. In practice the cardiologist wears a long wrap-around apron (front: 0.5 mm Pb equivalent, back: 0.25 mm Pb equivalent) and a thyroid collar (0.5 mm Pb equivalent). He is equipped with an electronic personal dosimeter (EPD Mk2, Siemens Environmental Systems Ltd, Poole, UK) sensitive to X-rays in the range of 0.02–6 MeV. It is attached to the apron, centred on and just under the right clavicle. For the simulations, either no protective clothing or one of three types of apron is added to ADAM. The first type is a frontal apron having the shape of the surface of half an elliptical cylinder. It covers the forward half of ADAM from the shoulders to the knees, including the arms. The second type is similar to the first one, but material is removed from left and right to expose the arms. Finally, the third type is a wrap-around apron, covering front, back and sides of the body but leaving the arms free. To achieve this a duplicate of the second apron is placed at the back, connecting lead sheets are inserted in the trunk between arm bones and ribs up to the armpits, and the hips and sides of the upper legs are covered. The thyroid collar shields the front of the neck, the lower jaw and the forward part of the trunk top between the shoulders. Instead of modelling the actual device, the dose at the EPD position is calculated as an HP(10) value, i.e. dose in muscle tissue at a depth of 10 mm below the surface. For this purpose, in the simulation a hemisphere (1 cm radius, soft tissue material) is placed at the EPD location on either the skin or the apron. An internal HP(10) value is also calculated, i.e. in soft tissue inside the phantom at 1 cm below the skin surface under the EPD position. As for the organ doses, MCNP yields HP(10) in the form of a dose conversion factor, i.e. normalized to DAP. It is derived from the MCNP output as "dose per photon at the origin of the hemisphere", i.e. dose at 10 mm depth in tissue, divided by the calculated "DAP per photon".

X-ray spectra
The relative intensities of the photons emitted from the point source, as a function of the energy, must be known for MC simulation. The X-ray spectra to be used in the MCNP simulations are generated with spectrum generator software [12], taking the X-ray voltage and Al filtration applied in the clinic as input parameters. The PCXMC code automatically generates its own spectra, based on the same input parameters.

A spectrum of scattered X-rays was required for the assessment of the dose to the cardiologist. It was obtained with MCNP by recording the energies of photons passing a surface between the patient and the cardiologist, when the former – in this case the newborn male phantom – is being exposed (PA) to the primary beam.

Patients and staff
The procedures to assess and treat congenital heart defects comprised diagnostic catheterizations (n=24), atrial septic defect (ASD) closure (n=14), radio-frequency (RF) ablation (n=10), balloon dilatation (n=11) and patent ductus arteriosus (DA) occlusion (n=10). The mean age of the patients is 5 years and weight 19 kg with ranges of 0–18 years and 2–68 kg, respectively.

Four types of imaging are distinguished for each procedure, i.e. posteroanterior (PA) and lateral (LLAT) continuous fluoroscopy, and PA and LLAT cine exposures for angiography (at 12.5 or 25 frames per second). The DAP during each view was recorded, as were the parameters needed for MC input.

Although four categories of specialists can be distinguished, i.e. cardiologist, echo-cardiographer, anaesthetist and nurse, only the first is considered here as they are assumed to be the most exposed. Various types of protective clothing are worn (lead equivalent apron: short or long, frontal or wrap-around; lead equivalent collar, sometimes). In contrast to the practice in the UK, the radiation exposure is measured with an EPD outside the apron.


    Results and discussion
 Top
 Abstract
 Introduction
 Methods
 Results and discussion
 Conclusions
 References
 
Doses to paediatric patients
To validate some of the PCXMC results used in [6], DCFs are recalculated with MCNP for three of the five catheterization procedures, each of which is applied to a patient in a different age class. The cases were selected because of the best match, in terms of age and weight, of the available phantoms for MCNP and the patients, who happened to be males. Patient and phantom data are shown for the three procedures in Table 1Go. The technical (tube voltage, filtration) and geometric (direction, FSD, field size) parameters of the beam are also presented, for the four views per procedure (LLAT or PA in cine or fluoroscopy mode). Next in Table 1Go the DCFs per view for effective dose are shown with their relative (statistical) uncertainty. The results from simulations with the MCNP and PCXMC codes are compared in terms of a ratio (which ideally should tend to 1) and as percent relative difference. MCNP tends to yield slightly lower DCFs. The relative difference has a median value of -2.5% (range 0.6% to -11%), and (in absolute sense) is larger than 10% in only two cases. As DAP per view was actually measured, the DCFs can be used to calculate the effective dose per view and per procedure, according to MC code. These results are also shown in Table 1Go. For total effective dose per procedure, the relative difference between MCNP and PCXMC is not larger (absolutely) than approximately 5%. From the viewpoint of radiation protection, disagreement at such level is not unacceptable, hence it can be concluded a posteriori that the choice of PCXMC for faster results is justified.


View this table:
[in this window]
[in a new window]
 
Table 1. Comparison of three paediatric cardiac catheterizations with respect to effective dose (conversion factor) to the patient according to MCNP and PCXMC

 
It might be interesting to have a closer look at any differences between PCXMC and MCNP at organ level. Two views are selected from Table 1Go, i.e. those for which the codes yield (absolutely) the smallest and largest difference in effective dose. These are, respectively, ASD closure/PA fluoroscopy and RF ablation/LA cine view. In Table 2Go the corresponding DCFs for separate organs and for effective dose are listed with their statistical uncertainties. The latter always are larger for PCXMC owing to fewer particles being simulated in this code. Again, for comparing the two codes the ratio and the relative difference are used. Ratios often deviate from the ideal value of 1 (the codes yielding equal results), even for the view for which both codes agree best on the effective dose. Several relative differences in organ DCF between the two MC codes may also seem rather large. This holds both for organs inside (e.g. heart) and outside (e.g. small intestines or kidneys) the primary beam. Except for the lowest doses, (organs far from the beam, like testes and urinary bladder, and therefore the least reliable results), corresponding organ DCFs remain of comparable magnitude. Thus, it can be concluded that both codes yield similar dose distributions over the body. Furthermore, in the present cases only a few organs (lungs, oesophagus, bone marrow) contribute substantially to the effective dose, which further limits the influence of many of the larger differences.


View this table:
[in this window]
[in a new window]
 
Table 2. Paediatric heart catheterizations; comparison of organ and effective dose conversion factors by MCNP and PCXMC for the two views with (absolutely) smallest or largest relative difference in effective dose conversion factor

 
A conclusion regarding calculation of individual organ doses is that the two MC codes may sometimes disagree. The main reason for deviations has not been further investigated, but probably the phantom geometry (hermaphrodite versus gender-specific) and perhaps material composition are factors of importance.

Dose to the cardiologist
Using two large-size phantoms in an MC simulation, i.e. one patient and one cardiologist, with all necessary internal organs to be defined, puts a substantial load on the computer resources. By its design purpose, PCXMC does not even allow such a configuration. MCNP does, but its performance would improve if a less complicated set-up of structures in the simulation were allowed. Therefore, tests were made to see if the cardiologist's organ doses (complete phantom) could be calculated with a simplified patient phantom. Only the external shape of the 5-year-old male phantom was used. The body was filled homogeneously with polymethylmethacrylate (PMMA). The PA fluoroscopy exposure of the simplified phantom was compared with the complete patient phantom. No protective clothing was used. The difference in calculated HP(10), just below the centre of the cardiologist's right clavicle, was only 4% (normalized to DAP: 3.07 and 2.96 µSV (Gy.cm2)-1, respectively, for simplified and full patient). However, effective dose to the cardiologist differed by 49% (4.17 and 2.80 µSv (Gy.cm2)-1 for simplified and full patient, respectively). Scatter from the cardiologist himself causes about 22% of the calculated HP(10) value, coming mainly from muscle (9%), skin (5%) and bone, i.e. ribs: 1.5%. The remaining 78% of the dose is due to scatter from the patient. In the detailed phantom the principal contributors are muscle (46.5%), lungs (11%), heart (7%), ribs (4%), liver (4%) and skin (2%). It was concluded that the patient phantom's composition has a strong impact on the scatter properties. Hence, the complete patient phantom had to be used.

The radiation burden to the cardiologist due to two heart catheterization procedures is shown in Table 3Go. If the cardiologist does not wear protective clothing he receives an effective dose that is either less than 1% (ASD closure) or less than 0.2% (diagnostic procedure) of the effective dose to the patient. The main contributions to the effective dose of the cardiologist come from testes, stomach, colon, bladder and lung doses. The percentages are view-dependent. The calculated HP(10) is only a rough approximation of the effective dose. In the case of ASD closure there is an overestimate (relative difference 19%). For the diagnostic views there is an underestimate (relative difference -43%), which from the viewpoint of radiation protection is worse than overestimation. It should be noted, however, that in this case the internal HP(10) values are listed. As for the other type of procedure, calculated HP(10) just outside the skin overestimates E (relative difference 79%). The reading of the EPD underestimates the effective dose in the absence of protective clothing. This may be caused by non-uniform directional sensitivity and/or insensitivity to lower X-ray energies.


View this table:
[in this window]
[in a new window]
 
Table 3. Paediatric heart catheterizations: comparison of two procedures and effect of shielding with respect to the radiation burden of the cardiologist

 
Kicken and Bos [13] have investigated the effectiveness of a 0.5 mm Pb equivalent apron in vascular radiology. Supplying staff with thermoluminescent dosimeters (TLD) of standard sensitivity, to be worn above and under the apron, they measured an overall attenuation factor of 120. With high sensitivity TLDs the attenuation ranged from 22 to 375, with mean value 77, for varying exposure conditions. The present attenuation of over 100 – (internal) HP(10) being reduced from 6.84 µSv to 0.06 µSv owing to the use of a 0.5 mm Pb equivalent apron in the diagnostic procedure – agrees well with these literature data.

As expected, by applying a 0.5 mm Pb equivalent frontal apron (type 1) in the case of diagnostic catheterization the effective dose to the cardiologist is dramatically reduced. The distribution of organ doses also changes, with principal contributions now coming from bone marrow, skeleton, skin, thyroid and brain. Thus, the EPD's above-mentioned underestimation of effective dose to the cardiologist changes to a 100 fold overestimation.

It is obvious that a correction must be applied to the reading of an EPD worn on the apron to obtain a proper indication of the magnitude of the effective dose. Huyskens et al [14] introduced the divider, multiplier and protection factor as characteristic quantities and calculated values for various exposure conditions with regard to shielding and beam orientation. In their model the assumption of a broad unidirectional X-ray beam was made. The divider is the number by which the individual dose equivalent measured on the apron has to be divided to obtain effective dose. The multiplier is a number by which the individual dose equivalent measured under the apron has to be multiplied to obtain the effective dose. Protection factor is the ratio of effective doses without and with the apron being used, which indicates the effectiveness of the shielding.

These quantities are listed in Table 4Go for various types of aprons and lead equivalent thickness. The quantities from which they are derived, i.e. HP(10) on and under the apron and effective dose (all normalized to DAP) are also shown. The PA fluoroscopy view of the diagnostic procedure performed with the male baby patient was chosen because it yields the highest contribution to the effective dose of the cardiologist.


View this table:
[in this window]
[in a new window]
 
Table 4. Paediatric heart catheterizations: characterization of effectiveness of protective measures with respect to the radiation burden of the cardiologist (in all cases DAP=3.52 Gy.cm2 and the effective dose to the patient is 327.3 µSv (Gy.cm2)-1)

 
In Table 4Go two beams are considered. Beam A consists of radiation scattered on and in the patient when the latter is exposed to the primary beam from the X-ray unit. This represents the actual situation of exposing the cardiologist. Beam B is rather "artificial" in this respect, being a broad unidirectional beam like the one used by Huyskens et al Beam B is aimed to be incident on the cardiologist at a 45° angle coming from front left (RPO, right posterior oblique view). The photon energies are according to the spectrum of scattered X-rays. The simulations with Beam B were conducted to establish that discrepancies in results observed for the present Beam A and Huyskens' results can be explained by the differences in the beam properties. As can be seen in Table 4Go the values of the characteristic quantities for Beam B closely match the ranges stated by Huyskens et al.

If no apron is worn, the effective dose calculated with the unidirectional Beam B is lower than for exposure to actual scatter as calculated with Beam A (1.52 versus 2.48 µSv (Gy cm2)-1). When an apron is worn the opposite is true (e.g. 0.21 versus 0.06 µSv (Gy cm2)-1 for the 0.25 mm lead-equivalent wrap-around apron). Thus, considering unidirectional exposure results in underestimating the protection that the apron actually offers. Depending on the apron properties, the protection factor is 2.5 to 7 times higher according to the more realistic simulation than when using Beam B. For instance the lower regions from which the scatter beam originates reduces the contribution of the thyroid to the effective dose. Consequently, the actual value of the quantity divider is higher. In other words, the overestimation of E by HP(10) measured on the apron is increased by a factor of 2 to 5 when exposure is approximated with a unidirectional beam.

The characteristic quantity that applies to the situation in the UK, where staff wear dosimeters underneath protective clothing, is the multiplier. According to the present study the dosimeter reading should be multiplied by 6 and 2.5 to yield appropriate estimates of effective dose for frontal and wrap-around aprons, respectively. The present values differ to some extent from the ones stated by Huyskens et al (Table 4Go), depending on apron and exposure conditions.

For a number of fluoroscopy procedures (70–110 kV) Faulkner et al [15] have measured doses with TLDs in and on a Rando-Alderson phantom exposed to radiation scattered from another anthropomorphic phantom in the primary beam. For commonly used aprons they found divider values of 2–60 and multiplier values of less than 7. The present study, under different exposure conditions (e.g. 59 kV), is in agreement with their results.

If no apron is used, the on-apron value of HP(10) in Table 4Go is lower than the under-apron value. The difference is due to the fact that the under-apron position is on the skin of the phantom and the on-apron position is 2 cm away, i.e. on the spot where the apron would be if it were present. The closer to the phantom, the larger is the backscatter component to the dose. Hence the under-apron value is higher.

That backscatter from the phantom is effectively stopped by the aprons can be deduced from the constant, lower value of HP(10) on the apron with the aprons present compared with the absence of an apron. (Note that the layer of lead itself will not generate much scattered radiation at 60 kV.)

It might be concluded, wrongly, from looking at the HP(10) under-apron values that only thickness and not the type of apron (frontal or wrap-around) matters. However, the influence of the type of apron is evident from the effective dose. Obviously, wrap-around and thicker aprons offer more protection. As expected, addition of a thyroid collar causes a further, though relatively small, improvement.

Padovani et al [16] have evaluated a number of empirical rules that are meant to estimate occupational effective dose from measured operational quantities.

The first holds for a single dosimeter positioned at the neck, measuring HP(10) on the apron. In that case, the divider value of 21 is to be used [17]. The shape of the Adam phantom is not well suited to calculate the individual dose equivalent at the neck in the situation that the radiation source is close to the phantom and has a low frontal position. The protruding trunk then forms an unnatural shield, hiding the "monitor" spot. Assuming that, in reality, the position – on-apron, either high on the chest or on the neck – does not affect the reading of the dosimeter too much, in the present study the rule is about correct for the 0.25 mm Pb equivalent frontal apron. Applying the rule to the other aprons would overestimate E by a factor of up to almost 4.

The second rule applies to twin dosimeter situations. The NCRP [17] recommends the Rosenstein-Webster [18] algorithm: Go


where E is the effective dose, Hwaist equals HP(10) under-apron at the waist and Hneck equals HP(10) on-apron at the neck. Since the dosimeters are positioned differently, this formula cannot be compared directly with the results of the present study. After comparison with literature data, Padovani et al state that, in general, both rules lead to an underestimation of effective dose. They find better agreement when applying a third rule, by Niklason et al [19]: Go


where k=0.02 or 0.06 in the presence or absence of a thyroid collar, respectively, and Hos is the HP(0.07) on the collar and Hu is the HP(10) under the apron at the waist. Niklason's rule is mentioned here only for the sake of completeness. Again due to the unsuitable shape of the anthropomorphic phantom in this situation, it is not possible to check the validity of the rule for the present source position, from which the neck is, incorrectly, "out of view". Substituting the chest value of HP(10) for Hos, the (normalized) effective doses derived with Niklason's rule become 0.07 and 0.06 µSv (Gy.cm2)-1 for, respectively, the frontal and wrap-around apron plus collar (all 0.5 mm Pb). These values underestimate or overestimate the corresponding ones calculated with MCNP, which amount to 0.12 and 0.04 µSv (Gy.cm2)-1.

A study by Mateya et al [20] reveals that two-dosimeter algorithms underestimate experimentally determined effective dose values from low-energy X-rays by a factor of 1.4 to 3.3. They conclude that accurate estimation of effective dose from personal dosimeters under conditions as described here remains problematic.

With a slight adaptation of the Adam phantom to smooth the transition from chest to neck, which Ulanovsky et al [21] have shown to be feasible, MC calculation can be of advantage for solving such problems, i.e. by helping to clarify the relationship between effective dose and personal dosimeter readings for cases like those described above.


    Conclusions
 Top
 Abstract
 Introduction
 Methods
 Results and discussion
 Conclusions
 References
 
The Monte Carlo code PCXMC has proved to be a useful tool for fast evaluation of the radiation burden to paediatric patients from catheterizations in cases of (suspected) congenital heart defects. Although for individual organs the calculated dose may sometimes differ considerably from the value obtained with the well-established general purpose MC code MCNP, with respect to effective dose, both codes show good agreement.

To generate appropriate scattered radiation for exposure of the cardiologist, the complete patient phantom is required. Using a simplified PMMA phantom of equal external shape is not sufficient.

Using a broad unidirectional beam to expose the cardiologist rather than radiation actually scattering off the patient yields overestimation of occupational effective dose and underestimation of the shielding offered by protective clothing.

The effectiveness of various types and thickness of protective clothing has been evaluated for one view of one catheterization procedure. The results of the calculations are supported reasonably well by experimental studies from the literature, although the conditions are not always exactly comparable.

The derivation of occupational effective dose from single or multiple dosimeter readings after interventional radiology and cardiac catheterizations requires further optimization, e.g. through MC simulation studies, to improve the accuracy of the estimates.


    Footnotes
 
This study was supported by a grant from The Netherlands Organisation for Health Research and Development (ZonMw, project number 97-18-003). Back

Received for publication October 10, 2002. Revision received May 6, 2003. Accepted for publication May 30, 2003.


    References
 Top
 Abstract
 Introduction
 Methods
 Results and discussion
 Conclusions
 References
 

  1. Vañó E, Arranz L, Sastre JM, Moro C, Ledo A, Garate MT, et al. Dosimetric and radiation protection considerations based on some cases of patient skin injuries in interventional cardiology. Br J Radiol 1998;71:510–6.[Abstract]
  2. ICRP. 1990 Recommendations of the International Commission on Radiological Protection. ICRP Publication 60. Ann ICRP 1991;21(1–3). New York NY, USA: Pergamon Press, 1991.
  3. Kling A, Barão F, Nakagawa M, Távora L, Vaz P, editors. Advanced Monte Carlo for radiation physics, particle transport simulation and applications. Proceedings of the Monte Carlo 2000 Conference; 2000 Oct 23–26; Lisbon, Portugal. Berlin, Germany: Springer Verlag, 2001.
  4. Briesmeister J, editor. MCNP: A general Monte Carlo N-particle transport code, version 4C. Manual LA-13709-M. Los Alamos NM, USA: Los Alamos National Laboratory, 2000.
  5. Tapiovaara M, Lakkisto M, Servomaa A. PCXMC-A PC-based Monte Carlo program for calculating patient doses in medical x-ray examinations. Helsinki, Finland: Finnish Centre for Radiation and Nuclear Safety (STUK), 1997.
  6. Spoelstra FM, Geleijns J, Helbing WA, Blom NA, Ottenkamp J, Schultz FW, et al. Radiation exposure to patients and staff during paediatric heart catheterisations. Submitted to Circulation, April 2002.
  7. ICRU. Determination of dose equivalents resulting from external radiation sources. ICRU Report 39. Bethesda MD, USA: International Commission on Radiation Units and Measurements, 1985.
  8. Servomaa A, Tapiovaara M. Organ dose calculation in medical x-ray examinations by the program PCXMC. Radiat Prot Dosim 1998;80:213–9.[Abstract]
  9. Hart D, Jones DG, Wall BF. Estimation of effective dose in diagnostic radiology from entrance surface dose and dose-area product measurements. Report NRPB-R262. Chilton, Didcot, UK: National Radiological Protection Board, 1994.
  10. Kramer R, Zankl M, Williams G, Drexler G. The calculation of dose from external photon exposure using reference human phantoms and Monte Carlo methods; Part I: The male (ADAM) and female (EVA) adult mathematical phantoms. Report S-885. Munich, Germany: GSF Forschungszentrum für Umwelt und Gesundheit, 1986.
  11. Eckerman KF, Cristy M, Ryman JC. The ORNL mathematical phantom series. http://homer.ornl.gov/vlab/mird2.pdf, Dec 1996.
  12. Cranley K, Gilmore BJ, Fogarty GWA, Desponds L. Catalogue of diagnostic x-ray spectra and other data. Report No 78. York, UK: Institute of Physics and Engineering in Medicine, 1997.
  13. Kicken PJH, Bos AJJ. Effectiveness of lead aprons in vascular radiology: results of clinical measurements. Radiology 1995;197:473–8.[Abstract/Free Full Text]
  14. Huyskens CJ, Franken Y, Hummel WA. Guidance on personal dosimetry for occupational exposure in interventional radiology. J Radiat Prot 1994;14:229–34.[CrossRef]
  15. Faulkner K, Marshall NW. The relationship of effective dose to personnel and monitor reading for simulated fluoroscopic irradiation conditions. Health Phys 1993;64:502–8.[Medline]
  16. Padovani R, Foti C, Malisan MR. Staff dosimetry protocols in interventional radiology. Radiat Prot Dosim 2001;94:193–6.[Abstract]
  17. NCRP. Use of personal monitors to estimate effective dose equivalent and effective dose to workers for external exposure to low-LET radiation. NCRP report 122. Bethesda MD, USA: National Council on Radiation Protection and Measurements, 1995.
  18. Rosenstein M, Webster EW. Effective dose to personnel wearing protective aprons during fluoroscopy and interventional radiology. Health Phys 1994;67:88–9.[Medline]
  19. Niklason LT, Marx MV, Chan HP. The estimation of occupational effective dose in diagnostic radiology with two dosimeters. Health Phys 1994;67:611–5.[Medline]
  20. Mateya CF, Claycamp HG. Phantom-derived estimation of effective dose equivalent from x rays with and without a lead apron. Health Phys 1998;74:489–91.[Medline]
  21. Ulanovsky AV, Eckerman KF. Modifications to the ORNL phantom series in simulation of the responses of thyroid detectors. Radiat Prot Dosim 1998;79:429–31.[Abstract]



This article has been cited by other articles:


Home page
Br. J. Radiol.Home page
A KARAMBATSAKIDOU, B SAHLGREN, B HANSSON, M LIDEGRAN, and A FRANSSON
Effective dose conversion factors in paediatric interventional cardiology
Br. J. Radiol., September 1, 2009; 82(981): 748 - 755.
[Abstract] [Full Text] [PDF]


Home page
Radiat Prot DosimetryHome page
V. Tsapaki, S. Kottou, S. Korniotis, N. Nikolaki, S. Rammos, and S. C. Apostolopoulou
RADIATION DOSES IN PAEDIATRIC INTERVENTIONAL CARDIOLOGY PROCEDURES
Radiat Prot Dosimetry, January 2, 2009; (2009) ncn313v1.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Neuroradiol.Home page
C.B. Nauer, F. Kellner-Weldon, G. Von Allmen, D. Schaller, and J. Gralla
Effective Doses from Scan Projection Radiographs of the Head: Impact of Different Scanning Practices and Comparison with Conventional Radiography
AJNR Am. J. Neuroradiol., January 1, 2009; 30(1): 155 - 159.
[Abstract] [Full Text] [PDF]


Home page
Radiat Prot DosimetryHome page
B. J. Irving, G. J. Maree, E. R. Hering, and T. S. Douglas
Radiation dose from a linear slit scanning X-ray machine with full-body imaging capabilities
Radiat Prot Dosimetry, July 1, 2008; 130(4): 482 - 489.
[Abstract] [Full Text] [PDF]


Home page
Radiat Prot DosimetryHome page
F. W. Schultz and J. Zoetelief
DOSEMETER READINGS AND EFFECTIVE DOSE TO THE CARDIOLOGIST WITH PROTECTIVE CLOTHING IN A SIMULATED INTERVENTIONAL PROCEDURE
Radiat Prot Dosimetry, April 3, 2008; (2008) ncn085v1.
[Abstract] [Full Text] [PDF]


Home page
Radiat Prot DosimetryHome page
K. Smans, E. Vano, R. Sanchez, F. W. Schultz, J. Zoetelief, T. Kiljunen, C. Maccia, H. Jarvinen, R. Bly, A. Kosunen, et al.
Results of a European survey on patient doses in paediatric radiology
Radiat Prot Dosimetry, March 1, 2008; 129(1-3): 204 - 210.
[Abstract] [Full Text] [PDF]


Home page
Radiat Prot DosimetryHome page
S. H. Park, J. K. Lee, and C. Lee
Dose conversion coefficients calculated using tomographic phantom, KTMAN-2, for X-ray examination of cardiac catheterisation
Radiat Prot Dosimetry, February 1, 2008; 128(3): 351 - 358.
[Abstract] [Full Text] [PDF]


Home page
Br. J. Radiol.Home page
T Siiskonen, M Tapiovaara, A Kosunen, M Lehtinen, and E Vartiainen
Monte Carlo simulations of occupational radiation doses in interventional radiology
Br. J. Radiol., June 1, 2007; 80(954): 460 - 468.
[Abstract] [Full Text] [PDF]


Home page
Br. J. Radiol.Home page
R R Bridcut, E Murphy, A Workman, P Flynn, and R J Winder
Patient dose from 3D rotational neurovascular studies
Br. J. Radiol., May 1, 2007; 80(953): 362 - 366.
[Abstract] [Full Text] [PDF]


Home page
Br. J. Radiol.Home page
D G W Onnasch, F K Schroder, G Fischer, and H-H Kramer
Diagnostic reference levels and effective dose in paediatric cardiac catheterization
Br. J. Radiol., March 1, 2007; 80(951): 177 - 185.
[Abstract] [Full Text] [PDF]


Home page
Radiat Prot DosimetryHome page
J. Th. M. Jansen and J. Zoetelief
Dose constraints and guidance for exposure of individuals knowingly and willingly helping in the support and comfort of individuals undergoing medical exposure
Radiat Prot Dosimetry, June 1, 2006; 118(3): 315 - 324.
[Abstract] [Full Text] [PDF]


Home page
Radiat Prot DosimetryHome page
F. W. Schultz and J. Zoetelief
Estimating effective dose for a cardiac catheterisation procedure with single or double personal dosemeters
Radiat Prot Dosimetry, May 1, 2006; 118(2): 196 - 204.
[Abstract] [Full Text] [PDF]


Home page
Radiat Prot DosimetryHome page
F. W. Schultz and J. Zoetelief
Dose conversion coefficients for interventional procedures
Radiat Prot Dosimetry, December 1, 2005; 117(1-3): 225 - 230.
[Abstract] [Full Text] [PDF]


Home page
Radiat Prot DosimetryHome page
J. N. Kroon
3-Dimensional rotational X-ray imaging, 3D-RX: image quality and patient dose simulation for optimisation studies
Radiat Prot Dosimetry, May 17, 2005; 114(1-3): 341 - 349.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Schultz, F W
Right arrow Articles by Zoetelief, J
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Schultz, F W
Right arrow Articles by Zoetelief, J


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
BJR DMFR IMAGING  ALL BIR JOURNALS