Search

This paper describes the latest developments at the Institute for Energy Technology (IFE) in Norway, in the field of real-time 3D (three-dimensional) radiation risk assessment for the support of work simulation in nuclear environments. 3D computer simulation can greatly facilitate efficient work planning, briefing, and training of workers. It can also support communication within and between work teams, and with advisors, regulators, the media and public, at all the stages of a nuclear installation's lifecycle. Furthermore, it is also a beneficial tool for reviewing current work practices in order to identify possible gaps in procedures, as well as to support the updating of international recommendations, dissemination of experience, and education of the current and future generation of workers.IFE has been involved in research and development into the application of 3D computer simulation and virtual reality (VR) technology to support work in radiological environments in the nuclear sector since the mid 1990s. During this process, two significant software tools have been developed, the VRdose system and the Halden Planner, and a number of publications have been produced to contribute to improving the safety culture in the nuclear industry.This paper describes the radiation risk assessment techniques applied in earlier versions of the VRdose system and the Halden Planner, for visualising radiation fields and calculating dose, and presents new developments towards implementing a flexible and up-to-date dosimetric package in these 3D software tools, based on new developments in the field of radiation protection. The latest versions of these 3D tools are capable of more accurate risk estimation, permit more flexibility via a range of user choices, and are applicable to a wider range of irradiation situations than their predecessors.

Of the biological effects of human exposure to radiofrequency and microwave radiation, the best-established are those due to elevation of tissue temperature. To prevent harmful levels of heating, restrictions have been proposed on the specific absorption rate (SAR). However, the relationship between SAR and temperature rise is not an invariant, since not only the heat capacity but also the efficiency of heat dissipation varies between different tissues and exposure scenarios. For small enough SAR, the relationship is linear and may be characterized by a 'heating factor' deltaT/SAR. Under whole-body irradiation the SAR may be particularly high in the ankles due to the concentration of current flowing through a relatively small cross-sectional area. In a previous paper, the author has presented calculations of the SAR distribution in a human leg in the high frequency (HF) band. In this paper, the heating factor for this situation is derived using a finite element approximation of the Pennes bioheat equation. The sensitivity of the results to different blood perfusion rates is investigated, and a simple local thermoregulatory model is applied. Both time-dependent and steady-state solutions are considered. Results confirm the appropriateness of the ICNIRP reference level of 100 mA on current through the leg, but suggest that at higher currents significant thermoregulatory adjustments to muscle blood flow will occur.

This article assesses the radiological environment at the nuclear site of the Mayak PA, Russian Federation, during and after the accident in 1957, the so-called 'Kyshtym Accident', and the radiation doses to those who participated in the eradication of its consequences. Based on numerous archival documents, this paper presents the radiation data for 1957-1960, including individual dosimetry monitoring data and estimated doses to the Mayak workers, as well as the to personnel in supporting organisations and the military involved in the remediation operations. From 1957-1959 some 38 500 individuals took part in the clean-up actions and remediation of contaminated areas of the Mayak PA industrial site after the accident, including individuals exposed at the time of the accident: Mayak PA employees, militarypersonnel, and civil construction workers. External equivalent doses to 10 500 individuals were estimated in the range of 220 to 265 mSv, while there were cases of doses up to 950 mSv and higher. The collective dose received during the accident and its aftermath was in the region of 7300 person-Sv. By October-December 1957, the collective dose was about 4500 person-Sv. Collective doses recorded in 1958 and 1959 amounted to 2250 person-Sv and 480 person-Sv, respectively.

In order to develop the basis of an approach, through which we can comprehend the relationships of radiation exposures and dose, and dose and radiological effects, large numbers of dosimetric data are required. It is complex to estimate the dosimetry involved due to enormous variability of organism shapes, masses, and their natural habitats. This paper aims to develop an empirical equation, which correlates dose conversion coefficients (DCCs) with organism volume V, density ρ, source energy E and curve fitting coefficients. Compared with the values calculated by a Monte Carlo method, the mean relative deviation of the empirical-fitted DCCs are reasonably low for gamma external exposure and beta internal exposure. The empirical-fitted values are very close to those calculated using a simplified anatomical phantom, the ERICA Tool or a voxel phantom.

At its 58th session in May 2011, the United Nations Committee on the Effects of Atomic Radiation (UNSCEAR) decided to carry out, once sufficient information was available, a full assessment of the levels of exposure and radiation risks attributable to the Fukushima accident. It envisages a preliminary document for consideration at its 59th session in May of 2012 and a more complete report for the 60th session of the Committee in 2013. This paper summarises the aims and objectives of the project, the scope, the working arrangements as well as the relation of the work to other activities.

In this study, a comprehensive electron dosimetry model of newborn skeletal tissues is presented. The model is constructed using the University of Florida newborn hybrid phantom of Lee et al (2007 Phys. Med. Biol. 52 3309-33), the newborn skeletal tissue model of Pafundi et al (2009 Phys. Med. Biol. 54 4497-531) and the EGSnrc-based Paired Image Radiation Transport code of Shah et al (2005 J. Nucl. Med. 46 344-53). Target tissues include the active bone marrow (surrogate tissue for hematopoietic stem cells), shallow marrow (surrogate tissue for osteoprogenitor cells) and unossified cartilage (surrogate tissue for chondrocytes). Monoenergetic electron emissions are considered over the energy range 1 keV to 10 MeV for the following source tissues: active marrow, trabecular bone (surfaces and volumes), cortical bone (surfaces and volumes) and cartilage. Transport results are reported as specific absorbed fractions according to the MIRD schema and are given as skeletal-averaged values in the paper with bone-specific values reported in both tabular and graphic format as electronic annexes (supplementary data). The method utilized in this work uniquely includes (1) explicit accounting for the finite size and shape of newborn ossification centers (spongiosa regions), (2) explicit accounting for active and shallow marrow dose from electron emissions in cortical bone as well as sites of unossified cartilage, (3) proper accounting of the distribution of trabecular and cortical volumes and surfaces in the newborn skeleton when considering mineral bone sources and (4) explicit consideration of the marrow cellularity changes for active marrow self-irradiation as applicable to radionuclide therapy of diseased marrow in the newborn child.

The recent development in linear accelerator control systems, named volumetric-modulated arc therapy (VMAT), has generated significant interest in arc-based intensity-modulated radiation therapy (IMRT). The VMAT delivery technique features simultaneous changes in dose rate, gantry angle and gantry rotation speed as well as multi-leaf collimator (MLC) leaf positions while radiation is on. In this paper, we describe a generalized VMAT planning tool that is designed to take full advantage of the capabilities of the new linac control systems. The algorithm incorporates all of the MLC delivery constraints such as restrictions on MLC leaf interdigitation and the MLC leaf velocity constraints. A key feature of the algorithm is that it is able to plan for both single- and multiple-arc deliveries. Compared to conventional step-and-shoot IMRT plans, our VMAT plans created using this tool can achieve similar or better plan quality with less MU and better delivery efficiency. The accuracy of the obtained VMAT plans is also demonstrated through plan verifications performed on an Elekta Synergy linear accelerator equipped with a conventional MLC of 1 cm leaf width using a PreciseBeam? VMAT linac control system.

Radiotherapy using kilovoltage x-rays in conjunction with contrast agents incorporated into the tumor, gold nanoparticles in particular, could represent a potential alternative to current techniques based on high-energy linear accelerators. In this paper, using the voxelized Zubal phantom in conjunction with the Monte Carlo code PENELOPE to model a prostate cancer treatment, it is shown that in combination with a 360 degrees arc delivery technique, tumoricidal doses of radiation can be delivered to deep-seated tumors while still providing acceptable doses to the skin and other organs at risk for gold concentrations in the tumor within the range of 7-10 mg-Au per gram of tissue. Under these conditions and using a x-ray beam with 90% of the fluence within the range of 80-200 keV, a 72 Gy physical absorbed dose to the prostate can be delivered, while keeping the rectal wall, bladder, skin and femoral heads below 65 Gy, 55 Gy, 40 Gy and 30 Gy, respectively. However, it is also shown that non-uniformities in the contrast agent concentration lead to a severe degradation of the dose distribution and that, therefore, techniques to locally quantify the presence of the contrast agent would be necessary in order to determine the incident x-ray fluence that best reproduces the dosimetry obtained under conditions of uniform contrast agent distribution.

Diagnostic radiological imaging such as conventional radiography, fluoroscopy and computed tomography (CT) examinations will continue to provide tremendous benefits in modern healthcare. The benefit derived by the patient should far outweigh the risk associated with a properly conducted imaging examination. Nonetheless, it is very important to be able to quantify the risk associated with any radiological examination of patients, and effective dose has been considered a useful indicator of patient exposure. Quantification of the risks associated with radiological imaging is very important as such information will be helpful to physicians and their patients for comparing risks from various imaging examinations and for making informed decisions whenever there is a need for any radiological imaging. The determination of equivalent and effective doses in diagnostic radiology is of interest as a basis for estimates of risk from medical exposures. In this paper we describe a simple computer program OrgDose, which calculates the doses to 27 organs in the body and then calculates the organ equivalent and effective doses and the risk from various procedures in the radiology department including conventional radiography, fluoroscopy and computed tomography examinations. The program will be a useful tool for the medical and paramedical personnel who are involved with assessing organ and effective doses and risks from diagnostic radiology procedures.

Treatment plan optimization is a multi-criteria process. Optimizing solely on one objective or on a sum of a priori weighted objectives may result in inferior treatment plans. Manually adjusting weights or constraints in a trial and error procedure is time consuming. In this paper we introduce a novel multi-criteria optimization approach to automatically optimize treatment constraints (dose-volume and maximum-dose). The algorithm tries to meet these constraints as well as possible, but in the case of conflicts it relaxes lower priority constraints so that higher priority constraints can be met. Afterwards, all constraints are tightened, starting with the highest priority constraints. Applied constraint priority lists can be used as class solutions for patients with similar tumour types. The presented algorithm does iteratively apply an underlying algorithm for beam profile optimization, based on a quadratic objective function with voxel-dependent importance factors. These voxel-dependent importance factors are automatically adjusted to reduce dose-volume and maximum-dose constraint violations.

Intensity-modulated radiation therapy (IMRT) uses radiation beams of varying intensities to deliver varying doses of radiation to different areas of the tissue. The use of IMRT has allowed the delivery of higher doses of radiation to the tumor and lower doses to the surrounding healthy tissue. It is not uncommon for head and neck tumors, for example, to have large treatment widths that are not deliverable using a single field. In such cases, the intensity matrix generated by the optimizer needs to be split into two or three matrices, each of which may be delivered using a single field. Existing field-splitting algorithms used the pre-specified arbitrary split line or region where the intensity matrix is split along a column, i.e., all rows of the matrix are split along the same column (with or without the overlapping of split fields, i.e., feathering). If three fields result, then the two splits are along the same two columns for all rows. In this paper we study the problem of splitting a large field into two or three subfields with the field width as the only constraint, allowing for an arbitrary overlap of the split fields, so that the total MU efficiency of delivering the split fields is maximized. Proof of optimality is provided for the proposed algorithm. An average decrease of 18.8% is found in the total MUs when compared to the split generated by a commercial treatment planning system and that of 10% is found in the total MUs when compared to the split generated by our previously published algorithm.

Currently, standards for the compliance testing of wireless devices are being extended to cover a wider frequency band and different usage patterns of mobile phones as well as of novel body-worn and handheld devices. As a consequence, not only the head but also strongly varying tissue distributions of the body are exposed to electromagnetic radiation. Several authors have reported changes in the SAR absorption of body tissue due to the presence of a low permittivity fat layer. This paper identifies two different effects which can lead to increased SAR in layered tissue in comparison to the SAR assessed using homogeneous tissue simulating liquid: (1) for larger distances between the tissue and the antenna, standing wave effects occur depending on the frequency and fat layer thickness. (2) In the very close near-field (distances approximately lambda/40), reactive E-field components lead to high local absorption in the skin. The latter effect occurs at lower frequencies and depends on the antenna type. Modification of the parameters of the homogeneous liquids cannot compensate for these effects. However, a conservative exposure estimate can be obtained by applying a multiplication factor between 1 and 3 to the values assessed using current experimental dosimetric techniques.

Margins about a target volume subject to external beam radiation therapy are designed to assure that the target volume of tissue to be sterilized by treatment is adequately covered by a lethal dose. Thus, margins are meant to guarantee that all potential variation in tumour position relative to beams allows the tumour to stay within the margin. Variation in tumour position can be broken into two types of dislocations, reducible and irreducible. Reducible variations in tumour position are those that can be accommodated with the use of modern image-guided techniques that derive parameters for compensating motions of patient bodies and/or motions of beams relative to patient bodies. Irreducible variations in tumour position are those random dislocations of a target that are related to errors intrinsic in the design and performance limitations of the software and hardware, as well as limitations of human perception and decision making. Thus, margins in the era of image-guided treatments will need to accommodate only random errors residual in patient setup accuracy (after image-guided setup corrections) and in the accuracy of systems designed to track moving and deforming tissues of the targeted regions of the patient's body. Therefore, construction of these margins will have to be based on purely statistical data. The characteristics of these data have to be determined through the central limit theorem and Gaussian properties of limiting error distributions. In this paper, we show how statistically determined margins are to be designed in the general case of correlated distributions of position errors in three-dimensional space. In particular, we show how the minimal margins for a given level of statistical confidence are found. Then, how they are to be used to determine geometrically minimal PTV that provides coverage of GTV at the assumed level of statistical confidence. Our results generalize earlier recommendations for statistical, central limit theorem-based recommendations for margin construction that were derived for uncorrelated distributions of errors (van Herk, Remeijer, Rasch and Lebesque 2000 Int. J. Radiat. Oncol. Biol. Phys. 47 1121-35; Stroom, De Boer, Huizenga and Visser 1999 Int. J. Radiat. Oncol. Biol. Phys. 43 905-19).

Many materials in everyday use such as window glass in homes and offices, glass in sunrooms and greenhouses, vehicle glass and some brands of sunscreens act as a barrier to the shorter UVB wavelengths while transmitting some of the longer UVA wavelengths. This paper reports on the erythemal exposures due to the UVA waveband encountered over a 12-month period for a solar zenith angle (SZA) range of 4 degrees to 80 degrees and the resulting times required for an erythemal exposure of one standard erythemal dose (SED) due to the erythemal exposures to the UVA wavelengths. The minimum time for an exposure of one SED due to the UVA wavelengths in winter is approximately double that what it is in summer. The time period of 40 to 60 min was the most frequent length of time for an exposure of one SED with 60 to 80 min the next frequent length of time required for a one SED exposure.

Because proton beams activate positron emitters in patients, positron emission tomography (PET) has the potential to play a unique role in the in vivo verification of proton radiotherapy. Unfortunately, the PET image is not directly proportional to the delivered radiation dose distribution. Current treatment verification strategies using PET therefore compare the actual PET image with full-blown Monte Carlo simulations of the PET signal. In this paper, we describe a simpler and more direct way to reconstruct the expected PET signal from the local radiation dose distribution near the distal fall-off region, which is calculated by the treatment planning programme. Under reasonable assumptions, the PET image can be described as a convolution of the dose distribution with a filter function. We develop a formalism to derive the filter function analytically. The main concept is the introduction of 'Q' functions defined as the convolution of a Gaussian with a powerlaw function. Special Q functions are the Gaussian itself and the error function. The convolution of two Q functions is another Q function. By fitting elementary dose distributions and their corresponding PET signals with Q functions, we derive the Q function approximation of the filter. The new filtering method has been validated through comparisons with Monte Carlo calculations and, in one case, with measured data. While the basic concept is developed under idealized conditions assuming that the absorbing medium is homogeneous near the distal fall-off region, a generalization to inhomogeneous situations is also described. As a result, the method can determine the distal fall-off region of the PET signal, and consequently the range of the proton beam, with millimetre accuracy. Quantification of the produced activity is possible. In conclusion, the PET activity resulting from a proton beam treatment can be determined by locally filtering the dose distribution as obtained from the treatment planning system. The filter function can be calculated analytically using convolutions of Gaussians and powerlaw functions.

For the last two decades, the organ and tissue equivalent dose as well as effective dose conversion coefficients recommended by the International Commission on Radiological Protection (ICRP) have been determined with exposure models based on stylized MIRD5-type phantoms representing the human body with its radiosensitive organs and tissues according to the ICRP Reference Man released in Publication No. 23, on Monte Carlo codes sometimes simulating rather simplified radiation physics and on tissue compositions from different sources. Meanwhile the International Commission on Radiation Units and Measurements (ICRU) has published reference data for human tissue compositions in Publication No. 44, and the ICRP has released a new report on anatomical reference data in Publication No. 89. As a consequence many of the components of the traditional stylized exposure models used to determine the effective dose in the past have to be replaced: Monte Carlo codes, human phantoms and tissue compositions. This paper presents results of comprehensive investigations on the dosimetric consequences to be expected from the replacement of the traditional stylized exposure models by the voxel-based exposure models. Calculations have been performed with the EGS4 Monte Carlo code for external and internal exposures to photons and electrons with the stylized, gender-specific MIRD5-type phantoms ADAM and EVA on the one hand and with the recently developed tomographic or voxel-based phantoms MAX and FAX on the other hand for a variety of exposure conditions. Ratios of effective doses for the voxel-based and the stylized exposure models will be presented for external and internal exposures to photons and electrons as a function of the energy and the geometry of the radiation field. The data indicate that for the exposure conditions considered in these investigations the effective dose may change between +60% and -50% after the replacement of the traditional exposure models by the voxel-based exposure models.

Automatic beam angle selection is an important but challenging problem for intensity-modulated radiation therapy (IMRT) planning. Though many efforts have been made, it is still not very satisfactory in clinical IMRT practice because of overextensive computation of the inverse problem. In this paper, a new technique named BASPSO (Beam Angle Selection with a Particle Swarm Optimization algorithm) is presented to improve the efficiency of the beam angle optimization problem. Originally developed as a tool for simulating social behaviour, the particle swarm optimization (PSO) algorithm is a relatively new population-based evolutionary optimization technique first introduced by Kennedy and Eberhart in 1995. In the proposed BASPSO, the beam angles are optimized using PSO by treating each beam configuration as a particle (individual), and the beam intensity maps for each beam configuration are optimized using the conjugate gradient (CG) algorithm. These two optimization processes are implemented iteratively. The performance of each individual is evaluated by a fitness value calculated with a physical objective function. A population of these individuals is evolved by cooperation and competition among the individuals themselves through generations. The optimization results of a simulated case with known optimal beam angles and two clinical cases (a prostate case and a head-and-neck case) show that PSO is valid and efficient and can speed up the beam angle optimization process. Furthermore, the performance comparisons based on the preliminary results indicate that, as a whole, the PSO-based algorithm seems to outperform, or at least compete with, the GA-based algorithm in computation time and robustness. In conclusion, the reported work suggested that the introduced PSO algorithm could act as a new promising solution to the beam angle optimization problem and potentially other optimization problems in IMRT, though further studies need to be investigated.

This paper outlines the development of a 2 mm resolution voxel model, NAOMI (aNAtOMIcal model), designed to be representative of the average adult female. The primary medical imaging data were derived from a high-resolution MRI scan of a 1.65 m tall, 23 year old female subject with a mass of 58 kg. The model was rescaled to a height of 1.63 m and a mass of 60 kg, the dimensions of the International Commission on Radiological Protection reference adult female. There are 41 tissue types in the model. The application of NAOMI to the calculations of induced current densities and electric fields from applied low frequency magnetic and electric fields is described. Comparisons are made with values from the male voxel model, NORMAN. The calculations were extended from 50 Hz up to 10 MHz. External field reference levels are compared with the ICNIRP guidelines.

The paper describes dosimetric models that allow the estimation of average radiation exposures to terrestrial biota due to environmental sources in the soil as well as internal uniform distributions of radionuclides. Simple three-dimensional phantoms for 13 faunal reference organisms are specified. The calculation of absorbed dose per unit source strength for these targets is based on photon and electron transport simulations using the Monte Carlo method. The presented absorbed dose rate conversion coefficients are derived for terrestrial reference species. This allows the assessment of internal exposure as well as external photon exposure depending on the nuclide, habitat, target size and environmental contamination. To enable the application of specific radiation weighting factors for α-, low energy β- (E 0 < 10 keV), β- and γ-radiations, their partial contributions to the total absorbed dose are provided separately. The coefficients for external exposure are listed for organisms living above the ground for an infinite plane source 3 mm deep in soil, as well as for a horizontally infinite volume source uniformly distributed to a depth of 10 cm. Furthermore, the coefficients are also presented for organisms living in a contaminated 50 cm thick soil layer. A multi-layer canopy model for plants is also described. The conversion coefficients are given for 3 H, 1 4 C, 4 0 K, 3 6 Cl, 5 9 , 6 3 Ni, 8 9 , 9 0 Sr, 9 4 Nb, 9 9 TC 1 0 6 Ru, 1 2 9 , 1 3 1 I, 1 3 4 , 1 3 5 , 1 3 7 Cs, 2 1 0 Po, 2 1 0 Pb, 2 2 6 Ra, 2 2 7 , 2 2 8 , 2 3 0 , 2 3 1 , 2 3 2 , 2 3 4 Th, 2 3 4 , 2 1 5 , 2 3 8 U, 2 3 8 , 2 3 9 , 2 4 0 , 2 4 1 Pu, 2 4 1 Am, 2 3 7 Np and 2 4 2 , 2 4 3 , 2 4 4 Cm, together with their progeny.

An epidemiological study of mortality and cancer rates in plutonium workers at the Sellafield nuclear installation has been carried out by researchers at the London School of Hygiene and Tropical Medicine. The study required the assessment of organ-specific doses from plutonium for more than 5000 workers over a period of 40 years. This was a major undertaking as it involved the reconstruction of annual received doses from the results of some 223,000 urine samples that had been provided by the workers in the study. This paper outlines the techniques and strategies adopted in order to generate best estimates of dose from the available data, and presents summaries of these doses.