Search NASA⌕ Search

SEARCH · Search NASA

Results for “combined risk segment map”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

Low Temperature Geothermal Play Fairway Analysis For The Appalachian Basin: Phase 1 Revised Report November 18, 2016

Geothermal energy is an attractive sustainable energy source. Yet project developers need confirmation of the resource base to warrant their time and financial resources. The Geothermal Play Fairway Analysis of the Appalachian Basin evaluated risk metrics that communicate the favorability of potential low-temperature geothermal energy resources in reservoirs more than 1000 m below the surface. This analysis is focused on the direct use of the heat, rather than on electricity production. Four risk factors of concern for direct-use geothermal plays in the Appalachian Basin portions of New York, Pennsylvania, and West Virginia are examined individually, and then in combination: 1) thermal resource quality, 2) natural reservoir quality, 3) induced seismicity, and 4) utilization opportunities. Uncertainty in the risk estimation is quantified. Based on these metrics, geothermal plays in the Appalachian Basin were identified as potentially viable for a variety of direct-use-heat applications. The methodologies developed in this project may be applied in other sedimentary basins as a foundation for low temperature (50-150 °C), direct use geothermal resource, risk, and uncertainty assessment. Three methods with which to combine the four risk factors were used. Among these, the averaging of the individual risk factors indicates the most favorable counties within the study area are the West Virginia counties of Monongalia, Harrison, Lewis (dubbed the Morgantown–Clarksburg play fairway), Putnam, and Kanawha (Charleston play fairway), the New York counties of Chemung and Steuben plus adjacent Bradford county in Pennsylvania (Corning–Ithaca play fairway), and the Pennsylvania counties of Mercer, Crawford, Erie, and Warren, and adjacent Chautauqua county in New York (together, the Meadville–Jamestown play fairway). These higher priority regions are surrounded by broader medium priority zones. Also worthy of additional exploration is a broad region near Pittsburgh Pennsylvania, for which the available geological data are insufficient to fully analyze the geological risks but yet the population is high. First, to assess the spatial variation in the depth to which one would need to drill to obtain geothermal temperatures that are useful to a future project, the project used bottom-hole temperature data from Appalachian Basin oil and gas exploration. These bottom hole temperature data are abundant but of low quality. Second, the project examined the potential for sufficient water flow rates through rocks to harvest heat from a geothermal well field, considering only natural reservoirs. This analysis provides a very incomplete picture of spatial variability of natural reservoirs because the oil and gas reservoir data lack key properties and are spatially biased toward those locations with profitable amounts of hydrocarbons in the rock pore spaces. Third, in light of the fact that earthquake activity has been induced in several states by subsurface work related to the oil and gas industry, this project examined the potential for similar activity in the Appalachian Basin. Acknowledging that data for such a task are insufficient, we utilized what was available: records of seismic activity, regional estimates of the orientations of stress in the rocks, and locations and orientations of zones of lateral change in rock properties at depths down to several kilometers below Earth’s surface. With these data, we created a first approximation of spatially variable risks for induced earthquakes. Because no data existed with which to test the reliability of these methods, the results have a high degree of uncertainty. Fourth, we examined the spatial variability of the above-the-ground factors that contribute to the economical viability of projects to tap low-temperature geothermal resources for direct-use. We worked principally with population density as a regionally known variable that would impact the cost of district heating. The resulting maps omit the costs of producing the hot water from the ground, because the below-ground costs are directly coupled to the thermal resource risk factor and natural reservoir risk factor – later analyses of those costs will be needed. The result of the district heating analysis is highly skewed: few census locations yielded a low estimated surface cost. The team also identified more than 165 prospects for high value direct-use geothermal energy opportunities such as industrial sites, university campuses, and federal facilities, among others. At the closure of this regional analysis, the most significant technical uncertainties are 1) reservoir distribution and capacities; 2) validity of thermal resource maps, and 3) the holistic estimation of Levelized Cost of Heat for favorable geological situations.

15 GEOTHERMAL ENERGY↗

Appalachian Basin Play Fairway Analysis: Revised 2016 Combined Risk Factor Analysis

This submission contains information used to compute the combined risk factors for deep geothermal energy opportunities in the Appalachian Basin, in the context of a the Play Fairway Analysis project. The risk factors are sedimentary rock reservoir quality, thermal resource quality, potential for induced seismicity, and utilization for direct-use heating of neighborhoods. The methods used to combine the risk factors included taking the average, the geometric mean, and the minimum of the four risk factors. Combined risk maps are provided for three different sedimentary rock reservoir metrics. Combined risk maps are also provided for the three geologic risk factors alone (thermal, reservoir, and seismic), and for the three risk factors that exclude reservoir quality (utilization, seismicity, and thermal qualities). The 2015 data submission should be visited to obtain associated shapefiles, which include: 1) definition of the High and Medium priority play fairways (Inner_Fairway, and Outer_Fairway), 2) definition of the US Census Places (usCensusPlaces), 3) places (cities) of interest in the region (Places_of_Interest) identified as geothermal play fairways, 4) the point centers of the raster cells (Raster_Center_Locations), and 5) locations of industries and special-use communities (e.g., colleges and military bases) identified as low temperature heat users (Industries). The 2015 submission also includes: 1) a methodology memo that explains how the risk factors were combined (GPFA-AB_combining_risk_factors.pdf), 2) the earthquake-based seismic risk map, and 3) supporting information with details of the calculations or processing used in generating these data files. More details on each file are given in the spreadsheet "list_of_contents.xlsx" in the folder "Supporting_Information". Code used to calculate values is available at https://github.com/calvinwhealton/geothermal_pfa under the folder "combining_metrics". Note that the 2016 code is currently under the branch named "combining_metrics_2016" in the folder called "combining_metrics". This branch may be merged with the master branch in the future. Many files contained within this submission update and replace the indicated files contained in: Cornell University. (2015). Risk Factor Analysis in Low-Temperature Geothermal Play Fairway Analysis for the Appalachian Basin (GPFA-AB) [data set]. Retrieved from https://gdr.openei.org/submissions/622. doi:10.15121/1261942

15 GEOTHERMAL ENERGY↗

GPFA-AB_Phase1RiskAnalysisTask5DataUpload

This submission contains information used to compute the risk factors for the GPFA-AB project (DE-EE0006726). The risk factors are natural reservoir quality, thermal resource quality, potential for induced seismicity, and utilization. The methods used to combine the risk factors included taking the product, sum, and minimum of the four risk factors. The files are divided into images, rasters, shapefiles, and supporting information. The image files show what the raster and shapefiles should look like. The raster files contain the input risk factors, calculation of the scaled risk factors, and calculation of the combined risk factors. The shapefiles include definition of the fairways, definition of the US Census Places, the center of the raster cells, and locations of industries. Supporting information contains details of the calculations or processing used in generating the files. An image of the raster will have the same name except *.png as the file ending instead of *.tif. Images with “fairways” or “industries” added to the name are composed of a raster with the relevant shapefile added. The file About_GPFA-AB_Phase1RiskAnalysisTask5DataUpload.pdf contains information the citation, special use considerations, authorship, etc. More details on each file are given in the spreadsheet “list_of_contents.csv” in the folder “SupportingInfo”. Code used to calculate values is available at https://github.com/calvinwhealton/geothermal_pfa under the folder “combining_metrics”.

15 GEOTHERMAL ENERGY↗

Risk Factor Analysis in Low-Temperature Geothermal Play Fairway Analysis for the Appalachian Basin (GPFA-AB)

This submission contains information used to compute the risk factors for the GPFA-AB project. The risk factors are natural reservoir quality, thermal resource quality, potential for induced seismicity, and utilization. The methods used to combine the risk factors included taking the product, sum, and minimum of the four risk factors. The files are divided into images, rasters, shapefiles, and supporting information. The image files show what the raster and shapefiles should look like. The raster files contain the input risk factors, calculation of the scaled risk factors, and calculation of the combined risk factors. The shapefiles include definition of the fairways, definition of the US Census Places, the center of the raster cells, and locations of industries. Supporting information contains details of the calculations or processing used in generating the files. An image of the raster will have the same name except *.png as the file ending instead of *.tif. Images with 'fairways' or 'industries' added to the name are composed of a raster with the relevant shapefile added. The file About_GPFA-AB_Phase1RiskAnalysisTask5DataUpload.pdf contains information the citation, special use considerations, authorship, etc. ***See 'GPFA-AB.zip' at bottom for compressed and organized version of the files associated with this submission*** **More details (including location) on each file are given in the spreadsheet 'list_of_contents.csv' in the folder 'SupportingInfo'** Code used to calculate values is available: https://github.com/calvinwhealton/geothermal_pfa under the folder 'combining_metrics' - *See link below*

15 GEOTHERMAL ENERGY↗

Geothermal play fairway analysis, part 1: Example from the Snake River Plain, Idaho

The Snake River Plain (SRP) volcanic province overlies the track of the Yellowstone hotspot, a thermal anomaly that extends deep into the mantle. Most of the area is underlain by a basaltic volcanic province that overlies a mid-crustal intrusive complex, which in turn provides the long-term heat flux needed to sustain geothermal systems. Previous studies have identified several known geothermal resource areas within the SRP. For the geothermal study presented herein, our goals were to: (1) adapt the methodology of Play Fairway Analysis (PFA) for geothermal exploration to create a formal basis for its application to geothermal systems, (2) assemble relevant data for the SRP from publicly available and private sources, and (3) build a geothermal PFA model for the SRP and identify the most promising plays, using GIS-based software tools that are standard in the petroleum industry. The study focused on identifying three critical resource parameters for exploitable hydrothermal systems in the SRP: heat source, reservoir and recharge permeability, and cap or seal. Data included in the compilation for heat source were heat flow, distribution and ages of volcanic vents, groundwater temperatures, thermal springs and wells, helium isotope anomalies, and reservoir temperatures estimated using geothermometry. Reservoir and recharge permeability was inferred from the analysis of stress orientations and magnitudes, post-Miocene faults, and subsurface structural lineaments based on magnetics and gravity data. Data for cap or seal included the distribution of impermeable lake sediments and clay-seal associated with hydrothermal alteration below the regional aquifer. These data were used to compile Common Risk Segment maps for heat, permeability, and seal, which were combined to create a Composite Common Risk Segment map for all southern Idaho that reflects the risk associated with geothermal resource exploration and identifies favorable resource tracks. Our regional assessment indicated that undiscovered geothermal resources may be located in several areas of the SRP. Two of these areas, the western SRP and Camas Prairie, were selected for more detailed assessment, during which heat, permeability, and seal were evaluated using newly collected field data and smaller grid parameters to refine the location of potential resources. These higher resolution assessments illustrate the flexibility of our approach over a range of scales.

54 ENVIRONMENTAL SCIENCES↗

Geothermal Play Fairway Analysis, Part 2: GIS methodology

Play Fairway Analysis (PFA) in geothermal exploration originates from a systematic methodology developed within the petroleum industry and is based on a geologic, geophysical, and hydrologic framework of identified geothermal systems. We tailored this methodology to study the geothermal resource potential of the Snake River Plain and surrounding region, but it can be adapted to other geothermal resource settings. We adapted the PFA approach to geothermal resource exploration by cataloging the critical elements controlling exploitable hydrothermal systems, establishing risk matrices that evaluate these elements in terms of both probability of success and level of knowledge, and building a code-based ‘processing model’ to process results. A geographic information system was used to compile a range of different data types, which we refer to as elements (e.g., faults, vents, heat flow, etc.), with distinct characteristics and measures of confidence. Discontinuous discrete data (points, lines, or polygons) for each element were transformed into continuous interpretive 2D grid surfaces called evidence layers. Because different data types have varying uncertainties, most evidence layers have an accompanying confidence layer which reflects spatial variations in these uncertainties. Risk layers, as defined here, are the product of evidence and confidence layers, and are the building blocks used to construct Common Risk Segment (CRS) maps for heat, permeability, and seal, using a weighted sum for permeability and heat, but a different approach with seal. CRS maps quantify the variable risk associated with each of these critical components. In a final step, the three CRS maps were combined into a Composite Common Risk Segment (CCRS) map, using a modified weighted sum, for results that reveal favorable areas for geothermal exploration. Additional maps are also presented that do not mix contributions from evidence and confidence (to allow an isolated view of evidence and confidence), as well as maps that calculate favorability using the product of components instead of a weighted sum (to highlight where all components are present). Our approach helped to identify areas of high geothermal favorability in the western and central Snake River Plain during the first phase of study and helped identify more precise local drilling targets during the second phase of work. By identifying favorable areas, this methodology can help to reduce uncertainty in geothermal energy exploration and development.

15 GEOTHERMAL ENERGY↗

MO-D-213-02: Quality Improvement Through a Failure Mode and Effects Analysis of Pediatric External Beam Radiotherapy

Purpose: To conduct a failure mode and effects analysis (FMEA) as per AAPM Task Group 100 on clinical processes associated with teletherapy, and the development of mitigations for processes with identified high risk. Methods: A FMEA was conducted on clinical processes relating to teletherapy treatment plan development and delivery. Nine major processes were identified for analysis. These steps included CT simulation, data transfer, image registration and segmentation, treatment planning, plan approval and preparation, and initial and subsequent treatments. Process tree mapping was utilized to identify the steps contained within each process. Failure modes (FM) were identified and evaluated with a scale of 1–10 based upon three metrics: the severity of the effect, the probability of occurrence, and the detectability of the cause. The analyzed metrics were scored as follows: severity – no harm = 1, lethal = 10; probability – not likely = 1, certainty = 10; detectability – always detected = 1, undetectable = 10. The three metrics were combined multiplicatively to determine the risk priority number (RPN) which defined the overall score for each FM and the order in which process modifications should be deployed. Results: Eighty-nine procedural steps were identified with 186 FM accompanied by 193 failure effects with 213 potential causes. Eighty-one of the FM were scored with a RPN > 10, and mitigations were developed for FM with RPN values exceeding ten. The initial treatment had the most FM (16) requiring mitigation development followed closely by treatment planning, segmentation, and plan preparation with fourteen each. The maximum RPN was 400 and involved target delineation. Conclusion: The FMEA process proved extremely useful in identifying previously unforeseen risks. New methods were developed and implemented for risk mitigation and error prevention. Similar to findings reported for adult patients, the process leading to the initial treatment has an associated high risk.

60 APPLIED LIFE SCIENCES↗

Common risk segment mapping: Streamlining exploration for carbon storage sites, with application to coastal Texas and Louisiana

Large-scale deployment of Carbon Capture and Storage (CCS) will require a commensurately large number of sites. Efficient screening methods are needed to create investment assurance and focus efforts on the most promising sites. The problem is similar to petroleum exploration, for which there are well-developed (though seldom published) workflows, including Common Risk Segment (CRS) mapping. In brief, the process requires 1) defining the key play elements; 2) identifying candidate geologic intervals for each; 3) creating fact-based maps for those intervals; 4) determining minimum criteria for the success of each element; 5) reinterpreting the fact-based maps in terms of chance of success; and 6) combining the individual maps to form a composite, basin-scale view of prospectivity. We adapt the CRS process to screening for CO 2 storage sites. Critically, we redefine the process in terms of cost of characterization and development, rather than chance of success. For illustration, we apply the process to the example of the Lower Miocene on the Texas and Louisiana Gulf Coast. We show that the predictions are consistent with historic hydrocarbon production volumes and rates. The power of the CRS method is that it creates a systematic approach to geologic evaluation and translates complex, multidimensional analysis into clear, graphical and easily comprehended business inputs. The results highlight sweet spots and identifies critical risks, suggesting a focus for further data collection and analysis. Furthermore, the method developed here can be applied to both surface and subsurface factors anywhere that there is interest in geologic storage of CO 2 .

54 ENVIRONMENTAL SCIENCES↗

NASA Langley Research Center Systems Analysis & Concepts Directorate Participation in the Exploration Systems Architecture Study

The NASA Langley Research Center (LaRC) Systems Analysis & Concepts Directorate (SACD) began studying human exploration missions beyond low Earth orbit (LEO) in the year 1999. This included participation in NASA s Decadal Planning Team (DPT), the NASA Exploration Team (NExT), Space Architect studies and Revolutionary Aerospace Systems Concepts (RASC) architecture studies that were used in formulating the new Vision for Space Exploration. In May of 2005, NASA initiated the Exploration Systems Architecture Study (ESAS). The primary outputs of the ESAS activity were concepts and functional requirements for the Crewed Exploration Vehicle (CEV), its supporting launch vehicle infrastructure and identification of supporting technology requirements and investments. An exploration systems analysis capability has evolved to support these functions in the past and continues to evolve to support anticipated future needs. SACD had significant roles in supporting the ESAS study team. SACD personnel performed the liaison function between the ESAS team and the Shuttle/Station Configuration Options Team (S/SCOT), an agency-wide team charged with using the Space Shuttle to complete the International Space Station (ISS) by the end of Fiscal Year (FY) 2010. The most significant of the identified issues involved the ability of the Space Shuttle system to achieve the desired number of flights in the proposed time frame. SACD with support from the Kennedy Space Center performed analysis showing that, without significant investments in improving the shuttle processing flow, that there was almost no possibility of completing the 28-flight sequence by the end of 2010. SACD performed numerous Lunar Surface Access Module (LSAM) trades to define top level element requirements and establish architecture propellant needs. Configuration trades were conducted to determine the impact of varying degrees of segmentation of the living capabilities of the combined descent stage, ascent stage, and other elements. The technology assessment process was developed and implemented by SACD as the ESAS architecture was refined. SACD implemented a rigorous and objective process which included (a) establishing architectural functional needs, (b) collection, synthesis and mapping of technology data, and (c) performing an objective decision analysis resulting in technology development investment recommendations. The investment recommendation provided budget, schedule, and center/program allocations to develop required technologies for the exploration architecture, as well as the identification of other investment opportunities to maximize performance and flexibility while minimizing cost and risk. A summary of the trades performed and methods utilized by SACD for the Exploration Systems Mission Directorate (ESAS) activity is presented along with how SACD is currently supporting the implementation of the Vision for Space Exploration.

Keyes, Jennifer↗

SU-E-J-100: The Combination of Deformable Image Registration and Regions-Of-Interest Mapping Technique to Accomplish Accurate Dose Calculation On Cone Beam Computed Tomography for Esophageal Cancer

Purpose: We introduce a new method combined with the deformable image registration (DIR) and regions-of-interest mapping (ROIM) technique to accurately calculate dose on daily CBCT for esophageal cancer. Methods: Patients suffered from esophageal cancer were enrolled in the study. Prescription was set to 66 Gy/30 F and 54 Gy/30 F to the primary tumor (PTV66) and subclinical disease (PTV54) . Planning CT (pCT) were segmented into 8 substructures in terms of their differences in physical density, such as gross target volume (GTV), venae cava superior (SVC), aorta, heart, spinal cord, lung, muscle and bones. The pCT and its substructures were transferred to the MIM software to readout their mean HU values. Afterwards, a deformable planning CT to daily KV-CBCT image registration method was then utilized to acquire a new structure set on CBCT. The newly generated structures on CBCT were then transferred back to the treatment planning system (TPS) and its HU information were overridden manually with mean HU values obtained from pCT. Finally, the treatment plan was projected onto the CBCT images with the same beam arrangements and monitor units (MUs) to accomplish dose calculation. Planning target volume (PTV) and organs at risk (OARs) from both of the pCT and CBCT were compared to evaluate the dose calculation accuracy. Results: It was found that the dose distribution in the CBCT showed little differences compared to the pCT, regardless of whether PTV or OARs were concerned. Specifically, dose variation in GTV, PTV54, PTV66, SVC, lung and heart were within 0.1%. The maximum dose variation was presented in the spinal cord, which was up to 2.7% dose difference. Conclusion: The proposed method combined with DIR and ROIM technique to accurately calculate dose distribution on CBCT for esophageal cancer is feasible.

60 APPLIED LIFE SCIENCES↗

Merging high resolution geophysical and geochemical surveys to reduce exploration risk at glass buttes, Oregon

The engineering and studies phase of the Glass Buttes project was aimed at reducing risk during the early stages of geothermal project development. The project’s inclusion of high resolution geophysical and geochemical surveys allowed Ormat to evaluate the value of these surveys both independently and in combination to quantify the most valuable course of action for exploration in an area where structure, permeability, and temperature are the most pressing questions. The sizes of the thermal anomalies at Glass Buttes are unusually large. Over the course of Phase I Ormat acquired high resolution LIDAR data to accurately map fault manifestations at the surface and collected detailed gravity and aeromagnetic surveys to map subsurface structural features. In addition, Ormat collected airborne hyperspectral data to assist with mapping the rock petrology and mineral alteration assemblages along Glass Buttes faults and magnetotelluric (MT) survey to try to better constrain the structures at depth. Direct and indirect identification of alteration assemblages reveal not only the geochemical character and temperature of the causative hydrothermal fluids but can also constrain areas of upflow along specific fault segments. All five datasets were merged along with subsurface lithologies and temperatures to predict the most likely locations for high permeability and hot fluids. The Glass Buttes temperature anomalies include 2 areas, totaling 60 km2 (23 mi2) of measured temperature gradients over 165° C/km (10° F/100ft). The Midnight Point temperature anomaly includes the Strat-1 well with 90°C (194 °F) at 603 m (1981 ft) with a 164 °C/km (10°F/100ft) temperature gradient at bottom hole and the GB-18 well with 71°C (160 °F) at 396 m (1300 ft) with a 182°C/km (11°F/100ft) gradient. The primary area of alteration and elevated temperature occurs near major fault intersections associated with Brothers Fault Zone and Basin and Range systems. Evidence for faulting is observed in each data set as follows. Field observations include fault plane orientations, complicated fault intersections, and hydrothermal alteration apparently pre-dating basalt flows. Geophysical anomalies include large, linear gradients in gravity and aeromagnetic data with magnetic lows possibly associated with alteration. Resistivity low anomalies also appear to have offsets associated with faulting. Hyperspectral and XRF identified alteration and individual volcanic flow units, respectively. When incorporated into a 3D geologic model, the fault intersections near the highest proven temperature and geophysical anomalies provide the first priority targets at Midnight Point. Ormat geologists selected the Midnight Point 52-33 drilling target based on a combination of pre-existing drilling data, geologic field work, geophysical interpretation, and geochemical analysis. Deep temperatures of well 52-33 was lower than anticipated. Temperature gradients in the well mirrored those found in historical drilling, but they decreased below 1500 ft and were isothermal below 2000 ft.

15 GEOTHERMAL ENERGY↗

Advancements in 3D Structural Analysis of Geothermal Systems

Robust geothermal activity in the Great Basin, USA is a product of both anomalously high regional heat flow and active fault-controlled extension. Elevated permeability associated with some fault systems provides pathways for circulation of geothermal fluids. Constraining the local-scale 3D geometry of these structures and their roles as fluid flow conduits is crucial in order to mitigate both the costs and risks of geothermal exploration and to identify blind (no surface expression) geothermal resources. Ongoing studies have indicated that much of the robust geothermal activity in the Great Basin is associated with high density faulting at structurally complex fault intersection/interaction areas, such as accommodation/transfer zones between discrete fault systems, step-overs or relay ramps in fault systems, intersection zones between faults with different strikes or different senses of slip, and horse-tailing fault terminations. These conceptualized models are crucial for locating and characterizing geothermal systems in a regional context. At the local scale, however, pinpointing drilling targets and characterizing resource potential within known or probable geothermal areas requires precise 3D characterization of the system. Employing a variety of surface and subsurface data sets, we have conducted detailed 3D geologic analyses of two Great Basin geothermal systems. Using EarthVision (Dynamic Graphics Inc., Alameda, CA) we constructed 3D geologic models of both the actively producing Brady’s geothermal system and a ‘greenfield’ geothermal prospect at Astor Pass, NV. These 3D models allow spatial comparison of disparate data sets in 3D and are the basis for quantitative structural analyses that can aid geothermal resource assessment and be used to pinpoint discrete drilling targets. The relatively abundant data set at Brady’s, ~80 km NE of Reno, NV, includes 24 wells with lithologies interpreted from careful analysis of cuttings and core, a 1:24,000 scale detailed geologic map and cross-sections, 2D seismic reflection profiles and other geophysical data, and downhole temperature data. The 3D geologic model based on these data consists of 61 fault planes, 25 distinct stratigraphic units, and 2 intrusive bodies. Geothermal fluids are produced from a left step-over/relay ramp within the Brady’s Fault Zone (BFZ). Under local stress conditions, fault segments that strike NNE-to-NE are most likely to slip and/or dilate, and therefore transmit geothermal fluids. The 3D model defines the locations of discrete fault intersections within the BFZ and indicates that the densest zones of structurally controlled fracture permeability are ~10-to-10s of meters in diameter and plunge ~55° NW-NNW beneath the heart of the BFZ step over. The locations of high intersection density, high fault slip and dilation tendency, high subsurface temperature, and lithologies known to support high fracture permeability are combined to produce 3D ‘fairway’ maps useful in both assessments of geothermal resource potential and for defining drilling targets. Astor Pass is located on the Pyramid Lake Paiute Reservation, ~80 km north of Reno, NV. It is a prospective ‘greenfield’ geothermal area, and thus subsurface data are relatively sparse. Available data include: two relatively deep wells (~1400 m) and one shallower well (~500 m) with lithologies interpreted from drill cuttings, several 2D seismic reflection profiles, a 1:24,000 scale geologic map and cross-section, a shallow temperature survey, and downhole temperature data. 3D modeling based on these data has defined 19 distinct fault planes and 16 stratigraphic units. Based on the stress field calculated from borehole breakouts, drilling induced tensile cracks and petal-centerline cracks in the two relatively deep wells, 3D slip and dilation tendency analysis indicates that northerly striking fault segments are most likely to slip and/or dilate, and therefore transmit geothermal fluids. Analysis of fault intersection density indicates that the highest density of structurally controlled permeability within the field lies in a narrow (10-to-10s of m) zone plunging moderately (~35°) to the NNW beneath Pleistocene tufa deposits. This zone of increased fracture density, which we interpret as the primary upflow zone, is controlled by the intersection of N-to-NNW striking normal faults and a WNW striking dextral fault zone and represents the most promising target for future drilling. Construction of a 3D geologic model involves integration of a variety of data into an internally consistent framework. A robust model allows for spatial comparison between the various types of data (structural, stratigraphic, geophysical, temperature, etc.) that are commonly used independently to site geothermal wells. Furthermore, highly detailed 3D geologic models provide the basis for additional quantitative analysis, including 3D fault slip and dilation tendency analysis and the precise location of structurally controlled permeability pathways. These analyses provide detailed information relating to the internal dynamics of geothermal systems and can mitigate the costs and risks of geothermal exploration and development by contributing to better well targeting and more accurate evaluations of resource potential.

15 GEOTHERMAL ENERGY↗

SU-E-J-141: Comparison of Dose Calculation On Automatically Generated MRBased ED Maps and Corresponding Patient CT for Clinical Prostate EBRT Plans

Purpose: To analyze the effect of computing radiation dose on automatically generated MR-based simulated CT images compared to true patient CTs. Methods: Six prostate cancer patients received a regular planning CT for RT planning as well as a conventional 3D fast-field dual-echo scan on a Philips 3.0T Achieva, adding approximately 2 min of scan time to the clinical protocol. Simulated CTs (simCT) where synthesized by assigning known average CT values to the tissue classes air, water, fat, cortical and cancellous bone. For this, Dixon reconstruction of the nearly out-of-phase (echo 1) and in-phase images (echo 2) allowed for water and fat classification. Model based bone segmentation was performed on a combination of the DIXON images. A subsequent automatic threshold divides into cortical and cancellous bone. For validation, the simCT was registered to the true CT and clinical treatment plans were re-computed on the simCT in pinnacle{sup 3}. To differentiate effects related to the 5 tissue classes and changes in the patient anatomy not compensated by rigid registration, we also calculate the dose on a stratified CT, where HU values are sorted in to the same 5 tissue classes as the simCT. Results: Dose and volume parameters on PTV and risk organs as used for the clinical approval were compared. All deviations are below 1.1%, except the anal sphincter mean dose, which is at most 2.2%, but well below clinical acceptance threshold. Average deviations are below 0.4% for PTV and risk organs and 1.3% for the anal sphincter. The deviations of the stratifiedCT are in the same range as for the simCT. All plans would have passed clinical acceptance thresholds on the simulated CT images. Conclusion: This study demonstrated the clinical usability of MR based dose calculation with the presented Dixon acquisition and subsequent fully automatic image processing. N. Schadewaldt, H. Schulz, M. Helle and S. Renisch are employed by Phlips Technologie Innovative Techonologies, a subsidiary of Royal Philips NV.

60 APPLIED LIFE SCIENCES↗

MCTP system model based on linear programming optimization of apertures obtained from sequencing patient image data maps

Purpose: The authors present a hybrid direct multileaf collimator (MLC) aperture optimization model exclusively based on sequencing of patient imaging data to be implemented on a Monte Carlo treatment planning system (MC-TPS) to allow the explicit radiation transport simulation of advanced radiotherapy treatments with optimal results in efficient times for clinical practice. Methods: The planning system (called CARMEN) is a full MC-TPS, controlled through aMATLAB interface, which is based on the sequencing of a novel map, called “biophysical” map, which is generated from enhanced image data of patients to achieve a set of segments actually deliverable. In order to reduce the required computation time, the conventional fluence map has been replaced by the biophysical map which is sequenced to provide direct apertures that will later be weighted by means of an optimization algorithm based on linear programming. A ray-casting algorithm throughout the patient CT assembles information about the found structures, the mass thickness crossed, as well as PET values. Data are recorded to generate a biophysical map for each gantry angle. These maps are the input files for a home-made sequencer developed to take into account the interactions of photons and electrons with the MLC. For each linac (Axesse of Elekta and Primus of Siemens) and energy beam studied (6, 9, 12, 15 MeV and 6 MV), phase space files were simulated with the EGSnrc/BEAMnrc code. The dose calculation in patient was carried out with the BEAMDOSE code. This code is a modified version of EGSnrc/DOSXYZnrc able to calculate the beamlet dose in order to combine them with different weights during the optimization process. Results: Three complex radiotherapy treatments were selected to check the reliability of CARMEN in situations where the MC calculation can offer an added value: A head-and-neck case (Case I) with three targets delineated on PET/CT images and a demanding dose-escalation; a partial breast irradiation case (Case II) solved with photon and electron modulated beams (IMRT + MERT); and a prostatic bed case (Case III) with a pronounced concave-shaped PTV by using volumetric modulated arc therapy. In the three cases, the required target prescription doses and constraints on organs at risk were fulfilled in a short enough time to allow routine clinical implementation. The quality assurance protocol followed to check CARMEN system showed a high agreement with the experimental measurements. Conclusions: A Monte Carlo treatment planning model exclusively based on maps performed from patient imaging data has been presented. The sequencing of these maps allows obtaining deliverable apertures which are weighted for modulation under a linear programming formulation. The model is able to solve complex radiotherapy treatments with high accuracy in an efficient computation time.

60 APPLIED LIFE SCIENCES↗

Evaluation of Functional Marrow Irradiation Based on Skeletal Marrow Composition Obtained Using Dual-Energy Computed Tomography

Purpose: To develop an imaging method to characterize and map marrow composition in the entire skeletal system, and to simulate differential targeted marrow irradiation based on marrow composition. Methods and Materials: Whole-body dual energy computed tomography (DECT) images of cadavers and leukemia patients were acquired, segmented to separate bone marrow components, namely, bone, red marrow (RM), and yellow marrow (YM). DECT-derived marrow fat fraction was validated using histology of lumbar vertebrae obtained from cadavers. The fractions of RM (RMF = RM/total marrow) and YMF were calculated in each skeletal region to assess the correlation of marrow composition with sites and ages. Treatment planning was simulated to target irradiation differentially at a higher dose (18 Gy) to either RM or YM and a lower dose (12 Gy) to the rest of the skeleton. Results: A significant correlation between fat fractions obtained from DECT and cadaver histology samples was observed (r=0.861, P<.0001, Pearson). The RMF decreased in the head, neck, and chest was significantly inversely correlated with age but did not show any significant age-related changes in the abdomen and pelvis regions. Conformity of radiation to targets (RM, YM) was significantly dependent on skeletal sites. The radiation exposure was significantly reduced (P<.05, t test) to organs at risk (OARs) in RM and YM irradiation compared with standard total marrow irradiation (TMI). Conclusions: Whole-body DECT offers a new imaging technique to visualize and measure skeletal-wide marrow composition. The DECT-based treatment planning offers volumetric and site-specific precise radiation dosimetry of RM and YM, which varies with aging. Our proposed method could be used as a functional compartment of TMI for further targeted radiation to specific bone marrow environment, dose escalation, reduction of doses to OARs, or a combination of these factors.

62 RADIOLOGY AND NUCLEAR MEDICINE↗