Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Based Procedures”

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.

1,198 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

TU-E-201-03: Eye Lens Dosimetry in Radiotherapy Using Contact Lens-Shaped Applicator

Madan M. Rehani, Massachusetts General Hospital and Harvard Medical School, Boston Methods for Eye Lens Dosimetry and Studies On Lens Opacities with Interventionalists Radiation induced cataract is a major threat among staff working in interventional suites. Nearly 16 million interventional procedures are performed annually in USA. Recent studies by the principal investigator’s group, primarily among interventional cardiologists, on behalf of the International Atomic Energy Agency, show posterior subcapsular (PSC) changes in the eye lens in 38–53% of main operators and 21–45% of support staff. These changes have potential to lead to cataract in future years, as per information from A-Bomb survivors. The International Commission on Radiological Protection has reduced dose limit for staff by a factor of 7.5 (from 150 mSv/y to 20 mSv/y). With increasing emphasis on radiation induced cataracts and reduction in threshold dose for eye lens, there is a need to implement strategies for estimating eye lens dose. Unfortunately eye lens dosimetry is at infancy when it comes to routine application. Various approaches are being tried namely direct measurement using active or passive dosimeters kept close to eyes, retrospective estimations and lastly correlating patient dose in interventional procedures with staff eye dose. The talk will review all approaches available and ongoing active research in this area, as well as data from surveys done in Europe on status of eye dose monitoring in interventional radiology and nuclear medicine. The talk will provide update on how good is Hp(10) against Hp(3), estimations from CTDI values, Monte Carlo based simulations and current status of eye lens dosimetry in USA and Europe. The cataract risk among patients is in CT examinations of the head. Since radiation induced cataract predominantly occurs in posterior sub-capsular (PSC) region and is thus distinguishable from age or drug related cataracts and is also preventable, actions on awareness can lead to avoidance or even prevention. Learning Objectives: To understand recent changes in eye lens dose limits and thresholds for tissue reactions To understand different approaches to dose estimation for eye lens To learn about challenges in eye lens opacities among staff in interventional fluoroscopy Di Zhang, Toshiba America Medical Systems, Tustin, CA, USA Eye lens radiation dose from brain perfusion CT exams CT perfusion imaging requires repeatedly exposing one location of the head to monitor the uptake and washout of iodinated contrast. The accumulated radiation dose to the eye lens can be high, leading to concerns about potential radiation injury from these scans. CTDIvol assumes continuous z coverage and can overestimate eye lens dose in CT perfusion scans where the table do not increment. The radiation dose to the eye lens from clinical CT brain perfusion studies can be estimated using Monte Carlo simulation methods on voxelized patient models. MDCT scanners from four major manufacturers were simulated and the eye lens doses were estimated using the AAPM posted clinical protocols. They were also compared to CTDIvol values to evaluate the overestimation from CTDIvol. The efficacy of eye lens dose reduction techniques such as tilting the gantry and moving the scan location away from the eyelens were also investigated. Eye lens dose ranged from 81 mGy to 279 mGy, depending on the scanner and protocol used. It is between 59% and 63% of the CTDIvol values reported by the scanners. The eye lens dose is significantly reduced when the eye lenses were not directly irradiated. CTDIvol should not be interpreted as patient dose; this study has shown it to overestimate dose to the eye lens. These results may be used to provide more accurate estimates of actual dose to ensure that protocols are operated safely below thresholds. Tilting the gantry or moving the scanning region further away from the eyes are effective for reducing lens dose in clinical practice. These actions should be considered when they are consistent with the clinical task and patient anatomy. Learning Objectives: To become familiar with method of eye dose estimation for patient in specific situation of brain perfusion CT To become familiar with level of eye lens radiation doses in patients undergoing brain perfusion MDCT To understand methods for reducing eye lens dose to patient Jong Min Park, Department of Radiation Oncology, Seoul National University Hospital, Seoul, Korea Eye lens dosimetry in radiotherapy using contact lens-shaped applicator Dose calculation accuracy of commercial treatment planning systems is relatively low at shallow depths. Therefore, in-vivo measurements are often performed in the clinic to verify delivered doses to eye lens which are located at shallow depth. Current in-vivo dosimetry for eye lens during radiotherapy is generally performed with small in-vivo dosimeters on the surface of patient eyelid. Since this procedure potentially contains considerable uncertainty, a contact lens-shaped applicator made of acrylic (lens applicator) was developed for in-vivo measurements of eye lens dose during radiotherapy to reduce uncertainty. The lens applicator allows the insertion of commercially available metal oxide semiconductor field effect transistor (MOSFET) dosimeters. Computed tomography (CT) images of an anthropomorphic phantom with and without the lens applicator were acquired. A total of 20 VMAT plans were delivered to an anthropomorphic phantom and the doses with the lens applicator and the doses at the surface of the eyelid were measured using both micro and standard MOSFET dosimeters. The differences in measured dose at the surface of the eyelid from the calculated lens dose were acquired. The differences between the measured and the calculated doses at the lens applicator, as well as the differences between the measured and the calculated doses at the surface of the eyelid were acquired. The statistical significance of the differences was analyzed. The average difference between the measured and the calculated dose with the lens applicator was 16.8 % ± 10.4 % with a micro MOSFET dosimeter and 16.6 % ± 10.9% with a standard MOSFET dosimeter. The average difference without the lens applicator was 35.9% ± 41.5% with micro MOSFET dosimeter and 42.9% ± 52.2% with standard MOSFET dosimeter. The maximum difference with micro MOSFET dosimeter was 46% with the applicator and 188.4% without the applicator. For the standard MOSFET dosimeter, the maximum difference was 44.4% with the applicator and 246.4% without the applicator. The lens applicator allowed reduction of the differences between the calculated and the measured dose during in-vivo measurement for the eye lens as compared to in-vivo measurement at the surface of the eyelid. Learning Objectives: To understand limitations of dose calculation with commercial treatment planning system for eye lens during radiotherapy To learn about current in-vivo dosimetry methods for eye lens in the clinic To understand limitations of in-vivo dosimetry for eye lens during radiotherapy Di Zhang is an employee of Toshiba America Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-E-201-00: Eye Lens Dosimetry for Patients and Staff

Madan M. Rehani, Massachusetts General Hospital and Harvard Medical School, Boston Methods for Eye Lens Dosimetry and Studies On Lens Opacities with Interventionalists Radiation induced cataract is a major threat among staff working in interventional suites. Nearly 16 million interventional procedures are performed annually in USA. Recent studies by the principal investigator’s group, primarily among interventional cardiologists, on behalf of the International Atomic Energy Agency, show posterior subcapsular (PSC) changes in the eye lens in 38–53% of main operators and 21–45% of support staff. These changes have potential to lead to cataract in future years, as per information from A-Bomb survivors. The International Commission on Radiological Protection has reduced dose limit for staff by a factor of 7.5 (from 150 mSv/y to 20 mSv/y). With increasing emphasis on radiation induced cataracts and reduction in threshold dose for eye lens, there is a need to implement strategies for estimating eye lens dose. Unfortunately eye lens dosimetry is at infancy when it comes to routine application. Various approaches are being tried namely direct measurement using active or passive dosimeters kept close to eyes, retrospective estimations and lastly correlating patient dose in interventional procedures with staff eye dose. The talk will review all approaches available and ongoing active research in this area, as well as data from surveys done in Europe on status of eye dose monitoring in interventional radiology and nuclear medicine. The talk will provide update on how good is Hp(10) against Hp(3), estimations from CTDI values, Monte Carlo based simulations and current status of eye lens dosimetry in USA and Europe. The cataract risk among patients is in CT examinations of the head. Since radiation induced cataract predominantly occurs in posterior sub-capsular (PSC) region and is thus distinguishable from age or drug related cataracts and is also preventable, actions on awareness can lead to avoidance or even prevention. Learning Objectives: To understand recent changes in eye lens dose limits and thresholds for tissue reactions To understand different approaches to dose estimation for eye lens To learn about challenges in eye lens opacities among staff in interventional fluoroscopy Di Zhang, Toshiba America Medical Systems, Tustin, CA, USA Eye lens radiation dose from brain perfusion CT exams CT perfusion imaging requires repeatedly exposing one location of the head to monitor the uptake and washout of iodinated contrast. The accumulated radiation dose to the eye lens can be high, leading to concerns about potential radiation injury from these scans. CTDIvol assumes continuous z coverage and can overestimate eye lens dose in CT perfusion scans where the table do not increment. The radiation dose to the eye lens from clinical CT brain perfusion studies can be estimated using Monte Carlo simulation methods on voxelized patient models. MDCT scanners from four major manufacturers were simulated and the eye lens doses were estimated using the AAPM posted clinical protocols. They were also compared to CTDIvol values to evaluate the overestimation from CTDIvol. The efficacy of eye lens dose reduction techniques such as tilting the gantry and moving the scan location away from the eyelens were also investigated. Eye lens dose ranged from 81 mGy to 279 mGy, depending on the scanner and protocol used. It is between 59% and 63% of the CTDIvol values reported by the scanners. The eye lens dose is significantly reduced when the eye lenses were not directly irradiated. CTDIvol should not be interpreted as patient dose; this study has shown it to overestimate dose to the eye lens. These results may be used to provide more accurate estimates of actual dose to ensure that protocols are operated safely below thresholds. Tilting the gantry or moving the scanning region further away from the eyes are effective for reducing lens dose in clinical practice. These actions should be considered when they are consistent with the clinical task and patient anatomy. Learning Objectives: To become familiar with method of eye dose estimation for patient in specific situation of brain perfusion CT To become familiar with level of eye lens radiation doses in patients undergoing brain perfusion MDCT To understand methods for reducing eye lens dose to patient Jong Min Park, Department of Radiation Oncology, Seoul National University Hospital, Seoul, Korea Eye lens dosimetry in radiotherapy using contact lens-shaped applicator Dose calculation accuracy of commercial treatment planning systems is relatively low at shallow depths. Therefore, in-vivo measurements are often performed in the clinic to verify delivered doses to eye lens which are located at shallow depth. Current in-vivo dosimetry for eye lens during radiotherapy is generally performed with small in-vivo dosimeters on the surface of patient eyelid. Since this procedure potentially contains considerable uncertainty, a contact lens-shaped applicator made of acrylic (lens applicator) was developed for in-vivo measurements of eye lens dose during radiotherapy to reduce uncertainty. The lens applicator allows the insertion of commercially available metal oxide semiconductor field effect transistor (MOSFET) dosimeters. Computed tomography (CT) images of an anthropomorphic phantom with and without the lens applicator were acquired. A total of 20 VMAT plans were delivered to an anthropomorphic phantom and the doses with the lens applicator and the doses at the surface of the eyelid were measured using both micro and standard MOSFET dosimeters. The differences in measured dose at the surface of the eyelid from the calculated lens dose were acquired. The differences between the measured and the calculated doses at the lens applicator, as well as the differences between the measured and the calculated doses at the surface of the eyelid were acquired. The statistical significance of the differences was analyzed. The average difference between the measured and the calculated dose with the lens applicator was 16.8 % ± 10.4 % with a micro MOSFET dosimeter and 16.6 % ± 10.9% with a standard MOSFET dosimeter. The average difference without the lens applicator was 35.9% ± 41.5% with micro MOSFET dosimeter and 42.9% ± 52.2% with standard MOSFET dosimeter. The maximum difference with micro MOSFET dosimeter was 46% with the applicator and 188.4% without the applicator. For the standard MOSFET dosimeter, the maximum difference was 44.4% with the applicator and 246.4% without the applicator. The lens applicator allowed reduction of the differences between the calculated and the measured dose during in-vivo measurement for the eye lens as compared to in-vivo measurement at the surface of the eyelid. Learning Objectives: To understand limitations of dose calculation with commercial treatment planning system for eye lens during radiotherapy To learn about current in-vivo dosimetry methods for eye lens in the clinic To understand limitations of in-vivo dosimetry for eye lens during radiotherapy Di Zhang is an employee of Toshiba America Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-E-201-01: Methods for Eye Lens Dosimetry and Studies On Lens Opacities with Interventionists

Madan M. Rehani, Massachusetts General Hospital and Harvard Medical School, Boston Methods for Eye Lens Dosimetry and Studies On Lens Opacities with Interventionalists Radiation induced cataract is a major threat among staff working in interventional suites. Nearly 16 million interventional procedures are performed annually in USA. Recent studies by the principal investigator’s group, primarily among interventional cardiologists, on behalf of the International Atomic Energy Agency, show posterior subcapsular (PSC) changes in the eye lens in 38–53% of main operators and 21–45% of support staff. These changes have potential to lead to cataract in future years, as per information from A-Bomb survivors. The International Commission on Radiological Protection has reduced dose limit for staff by a factor of 7.5 (from 150 mSv/y to 20 mSv/y). With increasing emphasis on radiation induced cataracts and reduction in threshold dose for eye lens, there is a need to implement strategies for estimating eye lens dose. Unfortunately eye lens dosimetry is at infancy when it comes to routine application. Various approaches are being tried namely direct measurement using active or passive dosimeters kept close to eyes, retrospective estimations and lastly correlating patient dose in interventional procedures with staff eye dose. The talk will review all approaches available and ongoing active research in this area, as well as data from surveys done in Europe on status of eye dose monitoring in interventional radiology and nuclear medicine. The talk will provide update on how good is Hp(10) against Hp(3), estimations from CTDI values, Monte Carlo based simulations and current status of eye lens dosimetry in USA and Europe. The cataract risk among patients is in CT examinations of the head. Since radiation induced cataract predominantly occurs in posterior sub-capsular (PSC) region and is thus distinguishable from age or drug related cataracts and is also preventable, actions on awareness can lead to avoidance or even prevention. Learning Objectives: To understand recent changes in eye lens dose limits and thresholds for tissue reactions To understand different approaches to dose estimation for eye lens To learn about challenges in eye lens opacities among staff in interventional fluoroscopy Di Zhang, Toshiba America Medical Systems, Tustin, CA, USA Eye lens radiation dose from brain perfusion CT exams CT perfusion imaging requires repeatedly exposing one location of the head to monitor the uptake and washout of iodinated contrast. The accumulated radiation dose to the eye lens can be high, leading to concerns about potential radiation injury from these scans. CTDIvol assumes continuous z coverage and can overestimate eye lens dose in CT perfusion scans where the table do not increment. The radiation dose to the eye lens from clinical CT brain perfusion studies can be estimated using Monte Carlo simulation methods on voxelized patient models. MDCT scanners from four major manufacturers were simulated and the eye lens doses were estimated using the AAPM posted clinical protocols. They were also compared to CTDIvol values to evaluate the overestimation from CTDIvol. The efficacy of eye lens dose reduction techniques such as tilting the gantry and moving the scan location away from the eyelens were also investigated. Eye lens dose ranged from 81 mGy to 279 mGy, depending on the scanner and protocol used. It is between 59% and 63% of the CTDIvol values reported by the scanners. The eye lens dose is significantly reduced when the eye lenses were not directly irradiated. CTDIvol should not be interpreted as patient dose; this study has shown it to overestimate dose to the eye lens. These results may be used to provide more accurate estimates of actual dose to ensure that protocols are operated safely below thresholds. Tilting the gantry or moving the scanning region further away from the eyes are effective for reducing lens dose in clinical practice. These actions should be considered when they are consistent with the clinical task and patient anatomy. Learning Objectives: To become familiar with method of eye dose estimation for patient in specific situation of brain perfusion CT To become familiar with level of eye lens radiation doses in patients undergoing brain perfusion MDCT To understand methods for reducing eye lens dose to patient Jong Min Park, Department of Radiation Oncology, Seoul National University Hospital, Seoul, Korea Eye lens dosimetry in radiotherapy using contact lens-shaped applicator Dose calculation accuracy of commercial treatment planning systems is relatively low at shallow depths. Therefore, in-vivo measurements are often performed in the clinic to verify delivered doses to eye lens which are located at shallow depth. Current in-vivo dosimetry for eye lens during radiotherapy is generally performed with small in-vivo dosimeters on the surface of patient eyelid. Since this procedure potentially contains considerable uncertainty, a contact lens-shaped applicator made of acrylic (lens applicator) was developed for in-vivo measurements of eye lens dose during radiotherapy to reduce uncertainty. The lens applicator allows the insertion of commercially available metal oxide semiconductor field effect transistor (MOSFET) dosimeters. Computed tomography (CT) images of an anthropomorphic phantom with and without the lens applicator were acquired. A total of 20 VMAT plans were delivered to an anthropomorphic phantom and the doses with the lens applicator and the doses at the surface of the eyelid were measured using both micro and standard MOSFET dosimeters. The differences in measured dose at the surface of the eyelid from the calculated lens dose were acquired. The differences between the measured and the calculated doses at the lens applicator, as well as the differences between the measured and the calculated doses at the surface of the eyelid were acquired. The statistical significance of the differences was analyzed. The average difference between the measured and the calculated dose with the lens applicator was 16.8 % ± 10.4 % with a micro MOSFET dosimeter and 16.6 % ± 10.9% with a standard MOSFET dosimeter. The average difference without the lens applicator was 35.9% ± 41.5% with micro MOSFET dosimeter and 42.9% ± 52.2% with standard MOSFET dosimeter. The maximum difference with micro MOSFET dosimeter was 46% with the applicator and 188.4% without the applicator. For the standard MOSFET dosimeter, the maximum difference was 44.4% with the applicator and 246.4% without the applicator. The lens applicator allowed reduction of the differences between the calculated and the measured dose during in-vivo measurement for the eye lens as compared to in-vivo measurement at the surface of the eyelid. Learning Objectives: To understand limitations of dose calculation with commercial treatment planning system for eye lens during radiotherapy To learn about current in-vivo dosimetry methods for eye lens in the clinic To understand limitations of in-vivo dosimetry for eye lens during radiotherapy Di Zhang is an employee of Toshiba America Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-E-201-02: Eye Lens Dosimetry From CT Perfusion Studies

Madan M. Rehani, Massachusetts General Hospital and Harvard Medical School, Boston Methods for Eye Lens Dosimetry and Studies On Lens Opacities with Interventionalists Radiation induced cataract is a major threat among staff working in interventional suites. Nearly 16 million interventional procedures are performed annually in USA. Recent studies by the principal investigator’s group, primarily among interventional cardiologists, on behalf of the International Atomic Energy Agency, show posterior subcapsular (PSC) changes in the eye lens in 38–53% of main operators and 21–45% of support staff. These changes have potential to lead to cataract in future years, as per information from A-Bomb survivors. The International Commission on Radiological Protection has reduced dose limit for staff by a factor of 7.5 (from 150 mSv/y to 20 mSv/y). With increasing emphasis on radiation induced cataracts and reduction in threshold dose for eye lens, there is a need to implement strategies for estimating eye lens dose. Unfortunately eye lens dosimetry is at infancy when it comes to routine application. Various approaches are being tried namely direct measurement using active or passive dosimeters kept close to eyes, retrospective estimations and lastly correlating patient dose in interventional procedures with staff eye dose. The talk will review all approaches available and ongoing active research in this area, as well as data from surveys done in Europe on status of eye dose monitoring in interventional radiology and nuclear medicine. The talk will provide update on how good is Hp(10) against Hp(3), estimations from CTDI values, Monte Carlo based simulations and current status of eye lens dosimetry in USA and Europe. The cataract risk among patients is in CT examinations of the head. Since radiation induced cataract predominantly occurs in posterior sub-capsular (PSC) region and is thus distinguishable from age or drug related cataracts and is also preventable, actions on awareness can lead to avoidance or even prevention. Learning Objectives: To understand recent changes in eye lens dose limits and thresholds for tissue reactions To understand different approaches to dose estimation for eye lens To learn about challenges in eye lens opacities among staff in interventional fluoroscopy Di Zhang, Toshiba America Medical Systems, Tustin, CA, USA Eye lens radiation dose from brain perfusion CT exams CT perfusion imaging requires repeatedly exposing one location of the head to monitor the uptake and washout of iodinated contrast. The accumulated radiation dose to the eye lens can be high, leading to concerns about potential radiation injury from these scans. CTDIvol assumes continuous z coverage and can overestimate eye lens dose in CT perfusion scans where the table do not increment. The radiation dose to the eye lens from clinical CT brain perfusion studies can be estimated using Monte Carlo simulation methods on voxelized patient models. MDCT scanners from four major manufacturers were simulated and the eye lens doses were estimated using the AAPM posted clinical protocols. They were also compared to CTDIvol values to evaluate the overestimation from CTDIvol. The efficacy of eye lens dose reduction techniques such as tilting the gantry and moving the scan location away from the eyelens were also investigated. Eye lens dose ranged from 81 mGy to 279 mGy, depending on the scanner and protocol used. It is between 59% and 63% of the CTDIvol values reported by the scanners. The eye lens dose is significantly reduced when the eye lenses were not directly irradiated. CTDIvol should not be interpreted as patient dose; this study has shown it to overestimate dose to the eye lens. These results may be used to provide more accurate estimates of actual dose to ensure that protocols are operated safely below thresholds. Tilting the gantry or moving the scanning region further away from the eyes are effective for reducing lens dose in clinical practice. These actions should be considered when they are consistent with the clinical task and patient anatomy. Learning Objectives: To become familiar with method of eye dose estimation for patient in specific situation of brain perfusion CT To become familiar with level of eye lens radiation doses in patients undergoing brain perfusion MDCT To understand methods for reducing eye lens dose to patient Jong Min Park, Department of Radiation Oncology, Seoul National University Hospital, Seoul, Korea Eye lens dosimetry in radiotherapy using contact lens-shaped applicator Dose calculation accuracy of commercial treatment planning systems is relatively low at shallow depths. Therefore, in-vivo measurements are often performed in the clinic to verify delivered doses to eye lens which are located at shallow depth. Current in-vivo dosimetry for eye lens during radiotherapy is generally performed with small in-vivo dosimeters on the surface of patient eyelid. Since this procedure potentially contains considerable uncertainty, a contact lens-shaped applicator made of acrylic (lens applicator) was developed for in-vivo measurements of eye lens dose during radiotherapy to reduce uncertainty. The lens applicator allows the insertion of commercially available metal oxide semiconductor field effect transistor (MOSFET) dosimeters. Computed tomography (CT) images of an anthropomorphic phantom with and without the lens applicator were acquired. A total of 20 VMAT plans were delivered to an anthropomorphic phantom and the doses with the lens applicator and the doses at the surface of the eyelid were measured using both micro and standard MOSFET dosimeters. The differences in measured dose at the surface of the eyelid from the calculated lens dose were acquired. The differences between the measured and the calculated doses at the lens applicator, as well as the differences between the measured and the calculated doses at the surface of the eyelid were acquired. The statistical significance of the differences was analyzed. The average difference between the measured and the calculated dose with the lens applicator was 16.8 % ± 10.4 % with a micro MOSFET dosimeter and 16.6 % ± 10.9% with a standard MOSFET dosimeter. The average difference without the lens applicator was 35.9% ± 41.5% with micro MOSFET dosimeter and 42.9% ± 52.2% with standard MOSFET dosimeter. The maximum difference with micro MOSFET dosimeter was 46% with the applicator and 188.4% without the applicator. For the standard MOSFET dosimeter, the maximum difference was 44.4% with the applicator and 246.4% without the applicator. The lens applicator allowed reduction of the differences between the calculated and the measured dose during in-vivo measurement for the eye lens as compared to in-vivo measurement at the surface of the eyelid. Learning Objectives: To understand limitations of dose calculation with commercial treatment planning system for eye lens during radiotherapy To learn about current in-vivo dosimetry methods for eye lens in the clinic To understand limitations of in-vivo dosimetry for eye lens during radiotherapy Di Zhang is an employee of Toshiba America Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Subtask 1.3 – Integrated Carbon Capture and Storage for North Dakota Ethanol Production

The Energy & Environmental Research Center (EERC), in partnership with Red Trail Energy, LLC (RTE), a North Dakota ethanol producer; the North Dakota Industrial Commission (NDIC) Renewable Energy Program (REP); and the U.S. Department of Energy (DOE), conducted a feasibility and implementation study for a commercial carbon capture and storage (CCS) effort. This subtask provided technical support and developed recommended practices of how small-scale industrial CO2 emitters (<1,000,000 tonnes of CO2 emitted annually) may economically deploy CCS. The 64-million-gallon dry mill RTE ethanol facility, emitting an average 180,000 tonnes of CO2 annually, was the subject of the case study. The positive outcome of this research, which shows technical and economic potential for ethanol–CCS in North Dakota, has resulted in RTE acquiring an approved Permit to Drill on December 2, 2019, for a stratigraphic test well in early 2020, a necessary step toward a North Dakota CO2 Storage Facility Permit (SFP) for the RTE CCS effort. Outcomes include 1) validating the Broom Creek Formation as a regional target for CCS, 2) determining the full carbon life cycle of an industrial fuel production facility with CCS, 3) developing a field implementation plan (FIP) for small-scale CCS, and 4) determining the validity and pathway for using CCS to meet low-carbon fuel (LCF) standards. The technical team included the EERC, RTE, Trimeric Corporation, Schlumberger Carbon Services, and Computer Modelling Group (CMG). Findings from activities conducted November 2016 – May 2020 are summarized as follows. Several key steps have been accomplished toward validating the Broom Creek Formation as a CO2 injection and storage target at the RTE CCS site, including verification of the presence and structure of sandstone layers that may comprise the potential CO2 storage reservoir and the several thousand feet of overlying confining zone. Work included site characterization using existing data, interpretation of seismic data within the study area, and geologic modeling simulation of CO2 injection. Interpreted results estimate 3000 feet of confining zone between the Broom Creek Formation (storage target) and the lowermost underground source of drinking water (e.g., the Fox Hills Formation). The thickness of the Broom Creek injection target varies 230–420 ft within the survey area. Results were used to inform the location of a stratigraphic test well and associated characterization test program. No impediments were identified within the targeted CO2 storage complex that would prevent the project from moving forward. A stratigraphic test well is the next step to validate these results and to acquire remaining data necessary to develop a North Dakota CO2 SFP application. The EERC generated full carbon life cycle estimates for CCS integration with the RTE ethanol facility. Results indicated that an average 40% reduction in CO2 emissions is possible through CCS implementation. Approximately half of the carbon in the overall life cycle for a dry mill ethanol plant is generated through the fermentation process and emitted to the atmosphere; this is the CO2 stream targeted for CCS. The remaining carbon is attributed to corn feedstock farming (diesel, fertilizer), energy for fuel processing (natural gas, electricity), and transportation (diesel) of the fuel product. Life cycle carbon estimates are also affected by the anticipated energy consumption of a potential capture facility, which depend on the type of CO2 product generated. For example, a 30%–40% net CO2 emission reduction is estimated if a liquefied CO2 facility were incorporated compared to a 40%–50% net CO2 reduction if a supercritical “injection-grade” CO2 product is generated; i.e., more energy is required to further refine the CO2 stream, affecting the full carbon life cycle estimates. A CCS FIP was developed and initiated at the RTE CCS site, resulting in the development of several guidance documents: a CO2 Capture Process Design Package, a North Dakota CO2 geologic Storage Permits Template, and a Public Outreach Package for CCS in North Dakota. General FIP components include CO2 capture system, pipeline and well designs; monitoring, verification, and accounting (MVA) plans; geologic characterization and testing programs; and permitting and outreach plans. Vendor bids were also acquired for the CO2 liquefaction facility. Near-surface characterization (groundwater and soil gas sampling) and geologic characterization (seismic survey) were initiated to inform development of a UIC Class VI-compliant MVA plan compliant with a North Dakota CO2 SFP. Designs (well and geologic testing) were completed for a stratigraphic test well compliant with a North Dakota CO2 SFP. In addition, the outreach plan was executed, including community open houses, meetings with city/county/state officials, and development of public materials. Although other entities continued to mature incentive programs in 2019–2020, California and the Internal Revenue Service (IRS) currently provide the most advanced economic opportunities for CCS integrated with fuel production. The California Low-Carbon Fuel Standard (LCFS) adopted a CCS Protocol in January 2019, allowing submittal of a design-based pathway (DBP) application for an approved temporary (not certified) carbon intensity value for a fully engineered facility. An ethanol–CCS DBP application to the California LCFS Program (officially approved February 28, 2020) was developed to show that the RTE CCS effort meets LCFS requirements. The approved DBP provides confidence to advance the project and supports potential investments. Other entities continue to mature incentive programs. The IRS issued guidance in February 2020 that addresses the definition of beginning of construction and revenue procedure on partnerships for the Enhancement of Carbon Dioxide Sequestration Credit (a.k.a. Section 45Q) CCS tax credit program; the IRS anticipates issuing further guidance on issues such as secure geologic storage, utilization qualifications, and recapture of claimed credits. Maturing incentive programs coupled with workable permitting regulations provide confidence to advance CCS projects in North Dakota and support financial investment to proceed with designing, constructing, and implementing CCS projects at small-scale fuel production facilities. The largest hurdles for CCS implementation at small-scale industrial systems are often business/economic-related (i.e., not technical). Market uncertainty already exists for agriculture-based alternative fuels such as corn ethanol, for which production has increased by ~33%, and prices have correspondingly lowered since 2015. Passing the Section 45Q tax credit program improves economic feasibility for CCS but may require external investors for a small business to achieve maximum benefits. Public–private partnerships with NDIC and DOE have resulted in foundational technical and regulatory knowledge, growing stakeholder confidence, and a pathway to implementation that enables similar industrial CCS projects in the region to advance. This subtask was funded through the EERC–DOE Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by NDIC and RTE. The authors would also like to thank CMG, ESRI, IHS, Neuralog, and Schlumberger for allowing the use of their software packages in support of this work.

42 ENGINEERING↗
DOE OSTI2022

Precursor Analysis Report: Cyber Attack on Thyssenkrupp Blast Furnace 2014

The Cyber Attack on Thyssenkrupp Blast Furnace 2014 Precursor Analysis Report leverages publicly available information about the Thyssenkrupp Steel Mill cyber attack and catalogs anomalous observables for each technique employed in the attack. This analysis is based upon the methodology of the Cybersecurity for the Operational Technology Environment (CyOTE) program. In December 2014, the German Government’s Federal Office for Information Security (BSI) released a report detailing a cyber attack on a German steel mill that occurred earlier that year, though exact dates and details of the attack were not revealed. While the report did not specify the name of the company, multiple sources identified the victim as one of Europe’s largest steel manufacturers, Thyssenkrupp AG. Further, Thyssenkrupp announced on 16 May of that year that Europe’s largest blast furnace, “Schwelgern 2,” located at its facility in Duisburg, Germany, would be offline for several weeks for repairs and upgrades, suggesting Schwelgern 2 was likely the target of the attack. The attack began in early 2014, when adversaries infiltrated the victim steel mill’s Information Technology (IT) network via a spearphishing campaign, then worked their way into the Operational Technology (OT) environment, where they executed software that caused denial of service, denial of control, and eventually a loss of control. This led to the blast furnace shutting down without proper safety procedures, resulting in catastrophic physical damage. No lives were lost in the incident, but ThyssenKrupp suffered $4 million in damage to the blast furnace and an additional $6 million in lost revenue. The adversaries required specialized knowledge and expertise in steel production, which enabled them to compromise a variety of internal systems and components across both IT and OT networks. The attack also demonstrated detailed knowledge of the industrial control systems (ICS) and production processes being used. This combination resulted in one of the earliest known publicly reported cybersecurity incidents resulting in physical damage to ICS equipment. Researchers and analysts identified 19 unique techniques (used in a sequence of 20 steps) utilized during the attack with a total of 454 observables using MITRE ATT&CK® for Industrial Control Systems. The CyOTE program assesses observables accompanying techniques used prior to the triggering event to identify opportunities to detect malicious activity. If observables accompanying the attack techniques are perceived and investigated prior to the triggering event, earlier comprehension of malicious activity can take place. Fifteen of the identified techniques used during the Thyssenkrupp cyber attack were precursors to the triggering event. Analysis identified 369 observables associated with these precursor techniques, 316 of which were assessed to have an increased likelihood of being perceived in the 120 days preceding the triggering event. The response and comprehension time could have been reduced if the observables had been identified earlier. The information gathered in this report contributes to a library of observables tied to a repository of artifacts, data sources, and technique detection references for practitioners and developers to support the comprehension of indicators of attack. Asset owners and operators can use these products if they experience similar observables or to prepare for comparable scenarios.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2019

Support for Reactor Operators in Case of Cyber-Security Threats (NEUP Final Report)

Increasing use of digital technology in nuclear power plants (NPPs) poses cyber-security as a crucial threat to public safety and to continuous energy production. Cyber-security risks are comprised of complex known and unknown interactions between various entities, system vulnerabilities, network protocols, human users and malicious attacks. There is little understanding or research geared towards plant operators’ response under cyber-security threats and operational procedures to cope with such threats. This is particularly critical when a cyber-security event masquerades as a safety incident or an evolving accident. Instead of leading operators to the remediation of the accident, the masked cyber event may lead them to circumvent the safety mechanisms of the plant. Given these challenges, the objective of this project is to develop a tool to characterize abnormal NPP events as “cyber” or “safety” incidents, and to further develop this tool as a real-time operator aid to assist in response to the incident. This NEUP Final Report summarizes the research efforts towards the above goals and the results in the following three aspects. First, a methodology for classifying abnormal NPP events as cyber-attacks or safety incidents was developed. This provides the basis for responding to potential cyber-attacks. In this methodology, the causal relationships between various variables (i.e., physical variables, network variables, and control variables) for a nuclear system under consideration are modeled using dynamic Bayesian networks. Based on observations of certain variables (i.e., the evidence), inference is performed to identify the state of hidden variables, which helps indicate the occurrence of cyber-attacks or safety incidents. The proposed methodology was applied to the full scope NPP simulator system at The Ohio State University as a case study. The result shows that the proposed methodology is capable of distinguishing between cyber-attacks and safety incidents. Second, a game-theoretic approach for responding to cyber-attacks on nuclear power plants was developed. This is the second major goal of the project, and the method is used to recommend optimal responses to reactor operators in the case of cyber-attacks. We consider two cases in this research. In the first case, the interaction between the defender (i.e., the reactor operator) and the attacker is modeled as a two-player, nonzero-sum, Markov game with an infinite time horizon. In the second case, we consider a semi-Markov game between the defender and the attacker within a finite time horizon. The model in the second case enables us to consider the time sensitive response of the defender. In each case, the method for deriving the Nash equilibrium is developed. The Nash Equilibrium of the game provides the valid prediction of both players’ actions because no single player can benefit from unilaterally deviating from the equilibrium policy if the other player adheres to his/hers, hence it provides the best response of the defender to cyber-attacks. The proposed method in each case is demonstrated through a simplified benchmark digital feedwater control system. Third, a prototype display is developed to integrate the above information. The prototype display depicts the communication information to the operator and incorporates the cyber-security tool developed by the team into the interface. It consists of the regular display in current control rooms, as well as a panel used to indicate the occurrence of cyber-attacks or safety events and a panel used to guide the operator’s response to detected cyber-attacks.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

WIPP PA Sensitivity Studies in Support of CRA-2014 Completeness Determination and Compliance Decision - 18209

The Waste Isolation Pilot Plant (WIPP) has been developed by the U.S. Department of Energy (DOE) for the geologic (deep underground) disposal of transuranic (TRU) waste. Containment of TRU waste at the WIPP is regulated by the U.S. Environmental Protection Agency (EPA). The DOE demonstrates compliance with the containment requirements by means of performance assessment (PA) calculations. WIPP PA calculations estimate the probability and consequence of potential radionuclide releases from the repository to the accessible environment for a regulatory period of 10,000 years after facility closure. Since TRU waste was first received at the WIPP in 1999, compliance recertification applications (CRAs) have been required by the EPA on five-year cycle. The latest CRA (CRA-2014) was submitted by the DOE to the EPA in 2014. As part of the application, a formal set of PA calculations (CRA-2014 PA) was performed and results were documented in detail. The CRA-2014 PA incorporated several parameter updates based on new data, as well as the implementation of a different model for panel closures, barriers planned for the opening of each waste panel. The results presented as part of the CRA-2014 showed via the required complementary cumulative distribution function (CCDF) release curves that the estimated releases were within compliance limits. As part of their process for reviewing the completeness of the CRA-2014 [1], the EPA formally requested that the DOE perform sensitivity studies based on EPA-requested changes to model parameters. The DOE performed and documented the results of three sensitivity studies described below: CRA14-SEN2 [2], CRA14-SEN3 [3], and CRA14-SEN4 [4]. In each case, the new analysis results were compared to CRA- 2014 PA (CRA14) results to determine the impact of the EPA-requested changes. CRA14-SEN2: the application of EPA-requested modified parameters in the operations and experimental (non-waste) areas of the repository to simulate an accelerated (instantaneous) creep closure, the inclusion of capillary pressure effects on relative permeability, and an increase in initial/residual brine saturation and residual gas saturation. The modifications to the repository model resulted in increased pressures and decreased brine saturations in waste areas and increased pressures and brine saturations in the operations and experimental areas. However, the overall effects on total high-probability mean releases from the repository were entirely insignificant. It was concluded that the modeling assumptions associated with the operations and experimental areas of the repository have an insignificant effect on the prediction of total releases from the repository and/or adequacy of the current (CRA14) model to demonstrate compliance with the regulatory limits. CRA14-SEN3: the application of EPA-requested modified parameters in the panel closure areas and operations and experimental (non-waste) areas of the repository to simulate an accelerated (instantaneous) creep closure, the inclusion of capillary pressure effects on relative permeability, and an increase in initial/residual brine saturation and residual gas saturation. The modifications to the repository model predominantly result in increased pressures and decreased brine saturations in waste areas, increased pressures and brine saturations in the operations and experimental areas, and decreased brine and gas flows across all panel closures. The pressure increases in repository waste regions yielded slightly decreased brine saturations (on average) in those areas. However, the overall effects on total high-probability mean releases from the repository are entirely insignificant. It was concluded that the modeling assumptions made for this sensitivity study which were associated with the panel closures and north end areas of the repository have only a small effect on the prediction of total releases from the repository and/or adequacy of the current (CRA14) model to demonstrate compliance with the regulatory limits. CRA14-SEN4: the application of several changes to the current WIPP PA model (CRA14), including: parameter changes, use of an updated code, and a revised computational grid for the BRAGFLO code. The modifications to the repository model resulted in increased releases in all primary release mechanisms. The impacts of each EPA-requested change to CRA14 were analyzed with respect to each release mechanism in the CRA14-SEN4 study. Overall, total high-probability predicted mean releases from the repository were increased by about 15%. Total low-probability predicted mean releases were increased by about 107%. It was concluded that the EPA-requested changes to the CRA14 result in increases to the predicted total releases from the repository. However, releases calculated in the CRA14-SEN4 analysis remained below regulatory limits, demonstrating continued compliance of the WIPP. In January 2017, the EPA determined that CRA-2014 was complete and then began a technical evaluation of CRA-2014 to determine whether the WIPP continued to comply with the radiation protection standards for disposal. During their evaluation, the EPA utilized the sensitivity study results to inform their ultimate determination in July 2017 that the WIPP continues to comply with regulatory standards. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Comparative Assessment of Experimental Testing of Instrument with an Embedded Digital Device Using Model-Based and Conventional Methods

The use of digital equipment has not been widespread in U.S. nuclear facilities. A major contributing factor to this situation arises from regulatory concerns about common cause failures (CCFs) stemming from the use of identical software in redundant systems, questions as to how to quantify the reliability of software-based products, validity of tools that are used for software verification and validation (V&V), and the effectiveness of software QA procedures. Compounding the concerns are sporadic incidents of digital equipment malfunctions in nuclear power plants, including lockups resulting from software memory leaks, feedwater control malfunctions due to excessive traffic on Ethernet networks, software logic errors, and other problems leading to plant trips. Additionally, the limited amount of demonstrated research and a lack of industry consensus on quantitative software reliability methodologies has hindered the widespread adoption of digital I&C in the nuclear industry. In a recent Regulatory Issue Summary (RIS), the U.S. Nuclear Regulatory Commission (NRC) identified concerns about the impact of greater use of instrumentation and control (I&C) equipment with an embedded digital device (EDD) [1]. Specifically, the NRC staff states that increased use of such devices “may increase a facility’s vulnerability to a CCF” (i.e., common-cause failure). The prevailing NRC guidance regarding software CCF (SCCF) identifies two design attributes that are acceptable for eliminating CCF concerns: (1) diversity or (2) testability (specifically, 100% testability). Either solution can result in high costs and remaining licensing uncertainty (i.e., how much diversity is enough? how to ensure test coverage of every possible sequence of device states?). Consequently, many utilities and reactor designers have limited or avoided more extensive use of digital technology to minimize licensing, scheduling, and financial risk. Without development of cost effective qualification methods to satisfy regulatory requirements and address the potential for CCF vulnerability associated with EDDs, the nuclear power industry may not be able to realize the benefits of digital technology achieved by other industries. In order to address these challenges, the U.S. Department of Energy (DOE) is sponsoring this research to assess the prospect of CCF in digital I&C equipment. The purpose of the current research is to develop an effective approach employing science-based methods to resolve concerns about CCF vulnerability that serve to inhibit deployment of advanced instrumentation (e.g., sensors, actuators, microcontrollers) with EDDs in nuclear power applications. The research objectives address the challenge of establishing high levels of safety and reliability assurance needed for the qualification of EDDs (e.g., microprocessors, programmable logic devices) that are subject to software design faults, complex failure modes, and CCF vulnerability. Specific objectives are: (1) assess the regulatory context for treatment of CCF vulnerability a graded approach to their qualification, (3) develop and extend model-based testing methods to enable effective demonstration of whether devices are subject to CCF, which may arise from vulnerabilities introduced at any stage of the design lifecycle, (4) establish a cost-effective testing framework that incorporates automation and test scenario prioritization, and (5) demonstrate the qualification approach through selection and testing of candidate digital device(s). This report documents an assessment of the results from experimental testing of a representative instrument with an embedded digital device. The experiment serves as a proof-of-concept demonstration of a model-based testing (MBT) framework. The objective of MBT is to ensure that an EDD undergoes suitably comprehensive testing to establish that it is not subject to SCCF. The approach by which MBT accomplishes this objective is through development and application of an effective test suite that provides the appropriate coverage to detect the full range of postulated faults that may arise at the requirements, design, and code (implementation) level of the EDD software. The findings of the experiment involve the application of the MBT framework by Ohio State University (OSU) through testing of the selected demonstration instrument (i.e., the prototype smart pressure sensor developed by Virginia Commonwealth University, VCU). In addition, the experiment findings also involve application of a more conventional black-box testing approach that serves as a baseline against which the MBT methodology can be compared. The primary outcome documented in this report is the comparative analysis of the capabilities demonstrated for each method. The results of the analysis demonstrate clear benefits that arise from the MBT methodology. The comparison indicates greater effectiveness from the more comprehensive evidence available through MBT and enhanced efficiency achieved through automation within the MBT framework. The structure of the report addresses the approach and application of each test method as well as the comparative analysis of the experimental results. Chapter 2 gives a summary overview of the demonstration approach. Chapter 3 covers the baseline testing and Chapter 4 documents the model-based testing. Chapter 5 describes the findings of the comparative analysis and presents the conclusions that can be drawn about the effectiveness and efficiency of the MBT methodology.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Quantitative Characterization of Impacts of Coupled Geomechanics and Flow on Safe and Permanent Geological Storage of CO 2 in Fractured Aquifers

This is the final scientific report for the award DE-FE0023305, entitled “Quantitative Characterization of Impacts of Coupled Geomechanics and Flow on Safe and Permanent Geological Storage of CO 2 in Fractured Aquifers.” The work has been divided into six tasks. In Task 2, we characterized rock properties, which are important when developing a quantitative approach for understanding and predicting geomechanical effects on large-scale CO 2 injection and long-term storage in the subsurface. Rock properties of interest for this characterization include porosity, permeability, elastic constants, strength, and heat capacity. We measured rock properties for three different rock types: concrete, sandstone and shale. These properties were obtained from the acoustic measurement, permeability and porosity measurement, Brazilian test, the uniaxial compressive test, and heat capacity measurement. Then, using the Brazilian technique to fracture samples, we measured their permeability under brine injection for varying effective stresses, followed by similar measurements of permeability under sc-CO 2 injection. Permeability to brine and sc-CO 2 decreased as the effective stress increased. The apparent permeability to sc-CO 2 was an order of magnitude less than that for brine, a consequence of CO 2 being the non-wetting fluid. In Task 3, we developed understandings and correlations for CO 2 injection pressure induced fracturing. We initially fractured four eight inch cubes of concrete with brine to establish a standard test procedure. Then, additional concrete samples were fractured using brine. Confining stresses were exerted on the samples and brine was injected at a constant rate into a borehole that was previously drilled into the concrete cube. The borehole pressure was measured and fracture initiation was identified as a peak in the pressure profile. That was followed by fracture propagation, which was identified as a plateau in the pressure profile, and finally the fracture reached the sample surface. The sample surfaces were photographed before and after fracturing, which enabled determination of where the fracture reached the sample surface. Afterwards, aqueous food dye solution was injected into the fractured concrete samples to color the fracture planes. The borehole was filled with food dye solution in advance, a gas pulse was transmitted to the borehole, and the borehole was then shut in to allow the food dye to be driven into existing fractures. Then, the concrete samples were broken down by high-pressure gas injection to reveal the geometry and morphology of the hydraulic fractures. In addition, acoustic measurements were conducted on multiple locations on each face of the concrete samples before and after injection for comparison. This comparison was used to confirm changes in the internal structure of the concrete sample, namely the formation of fractures. There were six concrete samples fractured with brine. Afterwards, concrete samples were fracture using sc-CO 2 . The main change in experimental procedure was the inclusion of temperature control, since in the field, we mostly encounter temperature and pressure conditions above the supercritical point of CO 2 . The concrete samples were pre-heated in an oven to elevate their temperature, and thermal tape was used to heat the CO 2 to a temperature above its supercritical temperature. Two concrete samples were fractured using liquid CO 2 , four concrete samples were fractured using gaseous CO 2 , and twenty two concrete samples were fractured using sc-CO 2 . For the concrete samples fractured using sc-CO 2 , three had pre-existing hydraulic fractures on the surface near the wellbore, and two were water saturated (as opposed to containing only air). Two concrete samples were composite samples, with a high permeability ball in the center, in order to represent CO 2 injection into a high permeability zone surrounded by a low permeability sealing formation. Later, we fractured five shale samples, obtained from the Niobrara shale outcrop, with sc-CO 2 . The injection rate was constant, confining stresses were applied, the shale samples were pre-heated in an oven to the desired temperature, and the sample was then fractured. An analysis of the failure type was done for many of the concrete samples, based on the predicted break down pressure for tensile and shear failure, and it was found that most samples fractured due to shear failure. Also, the presence of induced fractures has little to no effect on the fracture initiation pressure, but significantly changes the fracture morphology. The CO 2 fracturing in water saturated samples behaved much similar to brine induced fracturing. There was a significant difference of breakdown pressure for injection of CO 2 and brine, where the breakdown pressure for CO 2 induced fracturing is generally around the minimum horizontal stress and that for brine is much higher. Finally, the fracture orientation for most samples with large confining stress differential was perpendicular to the minimum horizontal stress. When the stress differential, especially the difference between the two directions with the smallest stresses was small, the fracture orientation was determined not only by the minimum stress direction, but also the magnitude of the confining stress, the tensile strength of the rock, and the breakdown pressure. In Task 4, we modified our coupled flow-geomechanical models to model fracture growth and propagation in storage formations and caprocks. These flow-geomechanical models are TOUGH2-CSM and TOUGH2-FLAC. The TOUGH2-CSM fluid and heat flow formulation is based on the TOUGH2 formulation of mass and energy conservation equations that govern fluid and heat flow in general multiphase, multicomponent, multi-porosity systems. The TOUGH2-CSM simplified geomechanical formulation is based on the linear theory of elasticity applied to multi-porosity non-isothermal (thermo-multi-poroelastic) media. We previously derived, from the fundamental relations of the linear theory of elasticity, an equation relating mean stress, pore pressures, temperatures, and body force (the Mean Stress equation) that we added to the TOUGH2-CSM fluid and heat flow equations. We extended the TOUGH2-CSM simplified geomechanical formulation to calculate the entire stress tensor by deriving equations for stress tensor components from derivatives of the Cartesian thermo-multi-poroelastic Navier equation components. In addition, we derived equations for stress tensor components in rz-coordinates. We solved these geomechanical equations using the integral finite difference method. This method utilizes momentum fluxes obtained by “factoring” a divergence operator out of each geomechanical equation. In Task 5, we did a literature survey to determine suitable correlations for fracture initiation, growth, and propagation. After surveying numerous papers, we modified TOUGH2-CSM to model stress dependent fracture initiation and growth using the Mohr-Coulomb criterion for shear failure in faults and caprock and the condition of negative minimum effective stress for tensile failure. Tensile failure results in a fracture whose width depends on the difference between fracture pressure and minimum stress and whose propagation is based on the difference between the stress intensity factor at the fracture tip and the rock toughness. The TOUGH-FLAC model is based on the FLAC3D ubiquitous joint model, extended for the modeling of permeability changes induced by changes in effective normal stress across fractures as well permeability changes caused by shearing along fractures. In Task 6, we validated our coupled flow-geomechanical models using analytical solutions and problems from the literature. We validated the TOUGH2-CSM modifications using analytical solutions to the displacement from a uniform load on a semi-infinite elastic medium and the two-dimensional Mandel-Cryer effect. Those were followed by a sample problem to test our stress tensor calculations, which was for injection into a single-phase reservoir with constant properties. Comparisons of our simulator to published results were done on for the depletion of a single-phase reservoir with stress dependent porosity, two examples from the In Salah gas project, an axisymmetric baserock–reservoir–caprock system with a normal faulting stress regime that had thermally induced shear stresses, simulations of CO 2 leakage through caprock in a two-dimensional reservoir in Cartesian coordinates, and predictions of caprock failure from CO 2 injection into an axisymmetric reservoir. Finally, we simulated CO 2 injection pressure-induced fracturing from one of our laboratory studies, Sample 39, and obtained a fracture profile. We also validated the TOUGH2-FLAC model for fluid driven fracture growth against solutions based on the KGD model and a case with an inclined fracture that was loaded from the boundaries to achieve wing-crack propagation. The strain softening tensile behavior and softening of modulus considering a simple damage approach was verified by a simulation representing hydraulic fracturing stress measurement around a vertical well. Finally, the model was validated against deep fracture zone opening and surface uplift at In Salah with inverses analysis used to improve the match of simulation with field data. In Task 7, we developed a scheme based on inverse modeling that can be used to identify caprock leakage. Injection of fluid into a reservoir results in a time varying pressure profile that depends on the properties of the rock and the fluid. The presence of significant caprock leakage would effect this pressure profile and we used inverse modeling, namely the Levenberg-Marquardt method, to determine the leakage location from the pressure profile. A two-dimensional and a three-dimensional example based on published data were presented. The published simulation was run and its results were considered as “field” data. In the two-dimensional example, a simulation of CO 2 injection into a two-dimensional aquifer-caprock system, the caprock had a vertical fault through which leakage could occur. The location of this leakage was obtained from our inverse modeling scheme. In the three-dimensional example, the simulation domain had four geological layers, one of which was the injection zone that contained a horizontal injection well. We simulated fractures as high permeability gridblocks, introduced a fracture that spanned the caprock depth, and ran a simulation to obtain “field” data. Then, using our inverse modeling scheme, we obtained the location of that fracture.

58 GEOSCIENCES↗