Search

Deep underground geologic formations are being evaluated for long-term storage of CO 2 , including large continental flood basalt formations. At the GHGT-11 and GHGT-12 conferences, progress was reported on the initial phases for Wallula Basalt Pilot demonstration test (located in Eastern Washington state), where nearly 1,000 metric tons of CO 2 was injected over a 3-week period during July/August 2013. The target CO 2 injection intervals were two permeable basalt interflow reservoir zones with a combined thickness of ∼20 m that occur within a layered basalt sequence between a depth of 830-890 m below ground surface. During the two-year post-injection monitoring period, downhole fluid samples were periodically collected, coupled with limited wireline borehole logging surveys that provided indirect evidence of on-going chemical geochemical reactions/alterations and CO 2 distribution. A detailed pre-closure field characterization program included downhole fluid sampling, and performance of hydrologic tests and wireline geophysical surveys. Side-wall cores also were retrieved from within the targeted injection zones. Visual observations of the core material identified small globular nodules, translucent to yellow in color, residing within vugs and small cavities, which were not evident in pre-injection side-wall cores obtained from the native basalt formation. Characterization by x-ray diffraction identified these nodular precipitates as ankerite, the identical iron and calcium rich carbonate observed to form in laboratory tests with Columbia River basalts. Isotopic characterization (δ 13 C, δ 18 O) conducted on the ankerite nodules indicate a distinct isotopic signature that is closely aligned with that of the injected CO 2 . Final post-injection wireline geophysical logging results also indicate the presence of free-phase CO 2 at the top of the two injection interflow zones, with no vertical migration of CO 2 above the injection horizons. These findings support previous assumptions regarding storage feasibility and rapid mineralization of CO 2 injected into a suitable basalt formation.

Carbon capture and storage (CCS) is one of a suite of technology options that might play a significant role in reducing greenhouse gas emissions. However, outside of traditional enhanced oil and gas recovery operations with a well established business model, CCS project deployments are struggling with adoption of a federal or international climate policy driver appearing unlikely for the foreseeable future. As part of a feasibility study for an industrial CCS project at the Boise White Paper mill in Washington State, a business model and CCS system design was developed that provided financial surety for the project developers while recognizing that uncertainty in revenue forecasts for CO2 storage would persist for the foreseeable future. Key to the business model was installation of a new 37 MWe biomass-fueled power island at the plant that would replace antiquated boilers and that could still supply the necessary steam to run pulp and paper production operations when capturing CO2 emissions under favorable market conditions for monetizing CO2 storage credits. Under unfavorable market conditions, CO2 capture would be suspended and excess power generated from the plant would reduce external electrical energy purchases required for paper mill operations. The net CO2 reduction for the project versus current operations is 1.0 MMT of CO2 per year with the CCS system online, and 139 ktons per year with the CCS system offline due to reduced natural gas co-firing. So, both operational modes offered a significant net reduction in CO2 emissions. Based on assumptions regarding electricity and natural gas pricing, and CO2 storage revenue generated through an open market mechanism priced at $15 per metric ton CO2, the payback period for the project capturing 0.5 MMT/yr (62% of total CO2 production) was estimated at 6.7 years versus 7.4 years with the CCS system offline. Geologic storage was evaluated in the deep flood basalt formations that dominate the storage capacity in Eastern Washington. Sub-basalt sediments were also examined at the site to provide alternative or supplemental storage capacity in the event that storage in the basalt sequences did not prove technically or economically feasible.

The U.S. Department of Energy Big Sky Regional Carbon Sequestration Partnership completed drilling the world’s first continental flood basalt sequestration pilot borehole to a total depth (TD) of 1253 m at a paper mill site near the town of Wallula located in Southeastern Washington State. Site suitability was assessed prior to drilling by acquisition, processing and analysis of a four-mile, five-line, three component seismic swath, which was processed as a single data-dense line. Analysis of the seismic survey data indicated absence of major geologic structures that would preclude CO2 injection at the site. Drilling of Wallula pilot borehole was initiated on January 13, 2009 and reached TD on April 6, 2009. Hydrogeologic information was obtained primarily during borehole drilling/advancement utilizing a progressive drill-and-test characterization strategy. A general decreasing transmissivity trend with depth pattern was observed, which is consistent with results exhibited for Columbia River basalt interflow zones at a number of other deep wells in the region. Based on the comparative results from 10 test intervals, a candidate injection test zone was identified between the general depth interval of ∼828 and 875 m bgs. Over this interval, three brecciated interflow zones were intersected and isolated for CO2 injection. The flow tops have moderate permeability (75 to 150 millidarcies) and are bounded by thick flow interiors that have extremely low (microdarcy) permeability. The borehole configuration established at the Wallula pilot site provides a unique opportunity to scientifically study the reservoir behaviour of three connected reservoir intervals confined between primary and secondary caprock zones. The permitting process for the CO2 injection has proceeded in accordance with formal rules for geologic sequestration projects enacted in June 2008 into the underground injection control program administered by the Washington State Department of Ecology. The permitting process is expected to conclude in October 2010 and injection would begin soon thereafter. Post-injection monitoring includes long-term sampling of water retrieved from the injection zone, shallow groundwater and soil gas monitoring, and PSInSAR.

Mafic continental flood basalts form a globally important, but under-characterized CO 2 sequestration target. The Columbia River Basalt Group (CRBG) in the northwestern U.S. is up to 5 km thick and covers over 168,000 km 2 . In India, flood basalts are 3 km thick and cover greater than 500,000 km 2 . Laboratory experiments demonstrate that the CRBG and other basalts react with formation water and super critical (sc) CO 2 to precipitate carbonates, thus adding a potential mineral trapping mechanism to the standard trapping mechanisms of most other types of CO 2 sequestration reservoirs. Brecciated tops of individual basalt flows in the CRBG form regional aquifers that locally have greater than 30% porosity and three Darcies of permeability. Porous flow tops are potential sites for sequestration of gigatons of sc CO 2 in areas where the basalts contain unpotable water and are at depths greater than 800 m. In this paper we report on the U.S. DOE Big Sky Regional Carbon Sequestration Partnership surface seismic and borehole geophysical characterization that supports a field test of capacity, integrity, and geochemical reactivity of CRBG reservoirs in eastern Washington, U.S.A. Traditional surface seismic methods have had little success in imaging basalt features in on-shore areas where the basalt is thinly covered by sediment. Processing of the experimental 6.5 km, 5 line 3C seismic swath included constructing an elastic wavefield model, identifying and separating seismic wave modes, and processing the swath as a single 2D line. Important findings include: (1) a wide variety of shear wave energy modes swamp the P-wave seismic records; (2) except at very short geophone offsets, ground roll overprints P-wave signal; and (3) because of extreme velocity contrasts, P-wave events are refracted at incidence angles greater than 7–15 degrees. Subsequent removal of S-wave and other noise during processing resulted in tremendous improvement in image quality. The application of wireline logging to onshore basalts is underexploited. Full waveform sonic logs and resistivity-based image logs acquired in the 1250 m basalt pilot borehole provide powerful tools for evaluating geomechanics and lithofacies. The azimuth of the fast shear wave is parallel to S H and records the changes through geologic time in basalt flow and tectonic stress tensors. Combined with image log data, azimuthal S-wave data provide a borehole technique for assessing basalt emplacement and cooling history that is related to the development of reservoirs and seals, as well as the orientation of tectonic stresses and fracture systems that could affect CO 2 transport or containment. Reservoir and seal properties are controlled by basalt lithofacies, and rescaled P- and S-wave slowness curves, integrated with image logs, provide a tool for improved recognition of subsurface lithofacies.

Numerical representation of the geologic framework and its hydrologic and geochemical properties is an integral part of all vadose zone flow and transport modeling. Historically, the geologic framework and the physical distribution of flow and transport properties have been represented by simple homogeneous and horizontally stratified hydrogeologic units. To capture more of the heterogeneity, small-scale variability, and uncertainty within a model, the physical and geochemical parameters may be represented by probability density functions loosely correlated to individual hydrogeologic units. As computer-processing capabilities have become more advanced, there has been more emphasis on improving spatial resolution and quantifying uncertainty in key model parameters. One of the more popular approaches has focused on geostatistical simulation of the flow and transport properties themselves, with little regard to the geologic strata and sedimentary sequences. However, newer approaches are focusing more on geostatistical simulation of the sequence-stratigraphic relations of lithofacies and the geostatistical distributions of flow and transport properties within those facies. These approaches require more rigorous quantitative treatment of geologic data than are normally supported by the mostly qualitative nature of borehole geologic information. Thus at the Hanford Site, efforts are being made to systematize borehole geologic data to be used in a moremore » quantitative way to define the spatial distribution of flow and transport properties in support of vadose flow and transport simulations. New detailed procedures translate qualitative descriptive information into categorical data, and translate inconsistent quantitative and semi-quantitative data into common parametric data sets. A geologic data management system is being developed to manage and integrate these new interpreted categorical data sets with other existing databases to support synergistic analysis and to improve numerical representation of the geologic architecture.« less

The spherically symmetric model (SSM) for electron--hydrogen scattering is one in which the Coulomb repulsion between the electrons, 2/r12, is replaced by 2/r>, where r> is the greater of r1 or r2. The partial wave decomposition of the wave function for this model is carried out with emphasis on L > 0 partial waves, in which case it is shown that singlet and triplet results are identical. The model, nevertheless, retains all inelastic and ionization channels, and therefore provides a nontrivial testing ground for methods of calculation. The decomposition results in two-dimensional partial differential equations (PDEs). Here numerical calculations are carried out in two ways using (i) a propagation method, and (ii) a noniterative PDE method. The latter is augmented by an integral formula, which greatly enhances the practicality of the method. Elastic phase shifts for L = 1,2,3 partial waves in the elastic region are presented. PACS No. 34.80

Experimental research work has been conducted and is undergoing at Pacific Northwest National Laboratory (PNNL) to address a variety of scientific issues related with the potential leaks of the carbon dioxide (CO2) gas from deep storage reservoirs. The main objectives of this work are as follows: • Develop a systematic understanding of how CO2 leakage is likely to influence pertinent geochemical processes (e.g., dissolution/precipitation, sorption/desorption and redox reactions) in the aquifer sediments. • Identify prevailing environmental conditions that would dictate one geochemical outcome over another. • Gather useful information to support site selection, risk assessment, policy-making, and public education efforts associated with geological carbon sequestration. In this report, we present results from experiments conducted at PNNL to address research issues related to the main objectives of this effort. A series of batch and column experiments and solid phase characterization studies (quantitative x-ray diffraction and wet chemical extractions with a concentrated acid) were conducted with representative rocks and sediments from an unconfined, oxidizing carbonate aquifer, i.e., Edwards aquifer in Texas, and a confined aquifer, i.e., the High Plains aquifer in Kansas. These materials were exposed to a CO2 gas stream simulating CO2 gas leaking scenarios, and changes in aqueous phase pHmore » and chemical composition were measured in liquid and effluent samples collected at pre-determined experimental times. Additional research to be conducted during the current fiscal year will further validate these results and will address other important remaining issues. Results from these experimental efforts will provide valuable insights for the development of site-specific, generation III reduced order models. In addition, results will initially serve as input parameters during model calibration runs and, ultimately, will be used to test model predictive capability and competency. The results from these investigations will provide useful information to support site selection, risk assessment, and public education efforts associated with geological, deep subsurface CO2 storage and sequestration.« less

This white paper accompanies a demonstration model that implements methods for the risk-informed design of monitoring, verification and accounting (RI-MVA) systems in geologic carbon sequestration projects. The intent is that this model will ultimately be integrated with, or interfaced with, the National Risk Assessment Partnership (NRAP) integrated assessment model (IAM). The RI-MVA methods described here apply optimization techniques in the analytical environment of NRAP risk profiles to allow systematic identification and comparison of the risk and cost attributes of MVA design options.

This report documents the efforts taken to develop a preliminary design for the first commercial-scale CO2 capture and sequestration (CCS) project associated with biomass power integrated into a pulp and paper operation. The Boise Wallula paper mill is located near the township of Wallula in Southeastern Washington State. Infrastructure at the paper mill will be upgraded such that current steam needs and a significant portion of the current mill electric power are supplied from a 100% biomass power source. A new biomass power system will be constructed with an integrated amine-based CO2 capture plant to capture approximately 550,000 tons of CO2 per year for geologic sequestration. A customized version of Fluor Corporation’s Econamine Plus™ carbon capture technology will be designed to accommodate the specific chemical composition of exhaust gases from the biomass boiler. Due to the use of biomass for fuel, employing CCS technology represents a unique opportunity to generate a net negative carbon emissions footprint, which on an equivalent emissions reduction basis is 1.8X greater than from equivalent fossil fuel sources (SPATH and MANN, 2004). Furthermore, the proposed project will offset a significant amount of current natural gas use at the mill, equating to an additional 200,000 tons of avoided CO2 emissions. Hence, the total net emissions avoided through this project equates to 1,100,000 tons of CO2 per year. Successful execution of this project will provide a clear path forward for similar kinds of emissions reduction that can be replicated at other energy-intensive industrial facilities where the geology is suitable for sequestration. This project also represents a first opportunity for commercial development of geologic storage of CO2 in deep flood basalt formations. The Boise paper mill site is host to a Phase II pilot study being carried out under DOE’s Regional Carbon Partnership Program. Lessons learned from this pilot study and other separately funded projects studying CO2 sequestration in basalts will be heavily leveraged in developing a suitable site characterization program and system design for permanent sequestration of captured CO2. The areal extent, very large thickness, high permeability in portions of the flows, and presence of multiple very low permeability flow interior seals combine to produce a robust sequestration target. Moreover, basalt formations are quite reactive with water-rich supercritical CO2 and formation water that contains dissolved CO2 to generate carbonate minerals, providing for long-term assurance of permanent sequestration. Sub-basalt sediments also exist at the site providing alternative or supplemental storage capacity.

Pacific Northwest National Laboratory (PNNL) has embarked on an initiative to develop world-class capabilities for performing experimental and computational analyses associated with geologic sequestration of carbon dioxide. The ultimate goal of this initiative is to provide science-based solutions for helping to mitigate the adverse effects of greenhouse gas emissions. This Laboratory-Directed Research and Development (LDRD) initiative currently has two primary focus areas—advanced experimental methods and computational analysis. The experimental methods focus area involves the development of new experimental capabilities, supported in part by the U.S. Department of Energy’s (DOE) Environmental Molecular Science Laboratory (EMSL) housed at PNNL, for quantifying mineral reaction kinetics with CO2 under high temperature and pressure (supercritical) conditions. The computational analysis focus area involves numerical simulation of coupled, multi-scale processes associated with CO2 sequestration in geologic media, and the development of software to facilitate building and parameterizing conceptual and numerical models of subsurface reservoirs that represent geologic repositories for injected CO2. This report describes work in support of the computational analysis focus area. The computational analysis focus area currently consists of several collaborative research projects. These are all geared towards the development and application of conceptual and numerical models for geologic sequestration of CO2. The software beingmore » developed for this focus area is referred to as the Geologic Sequestration Software Suite or GS3. A wiki-based software framework is being developed to support GS3. This report summarizes work performed in FY09 on one of the LDRD projects in the computational analysis focus area. The title of this project is Data Assimilation Tools for CO2 Reservoir Model Development. Some key objectives of this project in FY09 were to assess the current state-of-the-art in reservoir model development, the data types and analyses that need to be performed in order to develop and parameterize credible and robust reservoir simulation models, and to review existing software that is applicable to these analyses. This report describes this effort and highlights areas in which additional software development, wiki application extensions, or related GS3 infrastructure development may be warranted.« less

Ice coverage of the Arctic Ocean is predicted to become thinner and to cover less area with time. The combination of more ice-free waters for exploration and navigation, along with increasing demand for hydrocarbons and improvements in technologies for the discovery and exploitation of new hydrocarbon resources have focused attention on the hydrocarbon potential of the Arctic Basin and its margins. The purpose of this document is to 1) summarize results of a review of published hydrocarbon resources in the Arctic, including both conventional oil and gas and methane hydrates and 2) develop a set of digital maps of the hydrocarbon potential of the Arctic Ocean. These maps can be combined with predictions of ice-free areas to enable estimates of the likely regions and sequence of hydrocarbon production development in the Arctic. In this report, conventional oil and gas resources are explicitly linked with potential gas hydrate resources. This has not been attempted previously and is particularly powerful as the likelihood of gas production from marine gas hydrates increases. Available or planned infrastructure, such as pipelines, combined with the geospatial distribution of hydrocarbons is a very strong determinant of the temporal-spatial development of Arctic hydrocarbon resources. Significant unknowns decrease themore » certainty of predictions for development of hydrocarbon resources. These include: 1) Areas in the Russian Arctic that are poorly mapped, 2) Disputed ownership: primarily the Lomonosov Ridge, 3) Lack of detailed information on gas hydrate distribution, and 4) Technical risk associated with the ability to extract methane gas from gas hydrates. Logistics may control areas of exploration more than hydrocarbon potential. Accessibility, established ownership, and leasing of exploration blocks may trump quality of source rock, reservoir, and size of target. With this in mind, the main areas that are likely to be explored first are the Bering Strait and Chukchi Sea, in spite of the fact that these areas do not have highest potential for future hydrocarbon reserves. Opportunities for improving the mapping and assessment of Arctic hydrocarbon resources include: 1) Refining hydrocarbon potential on a basin-by-basin basis, 2) Developing more realistic and detailed distribution of gas hydrate, and 3) Assessing the likely future scenarios for development of infrastructure and their interaction with hydrocarbon potential. It would also be useful to develop a more sophisticated approach to merging conventional and gas hydrate resource potential that considers the technical uncertainty associated with exploitation of gas hydrate resources. Taken together, additional work in these areas could significantly improve our understanding of the exploitation of Arctic hydrocarbons as ice-free areas increase in the future.« less

We are using viral oncogene probes to study the pathways by which osteoblast-specific gene expression is induced in ascorbic acid-treated MC3T3-E1 cells. The 12S product of the adenovirus E1A gene binds directly to key cellular regulators and, as a result, represses tissue specific gene expression and blocks differentiation in a wide variety of cell types. The main cellular targets of the E1A 12S product are the pRB family and p300/CBP family. The p300 family appears to be the primary target for E1A-mediated repression of tissue-specific gene expression in a variety of cell types. We have generated MC3T3-E1 cell lines that stably express either the wild-type 12S product or a mutant that targets p300/CBP, but not the pRB family. Using these constructs to dissect osteoblast differentiation, we found that targeting of p300/CBP appears to be sufficient to repress alkaline phosphatase expression, although a low but functional level of expression can be maintained if the pRB family is not targeted as well. Induction of alkaline phosphatase expression and activity can be dissociated from expression of late-stage markers such as osteocalcin and osteopontin. Surprisingly, cell lines exhibiting severe repression of alkaline phosphatase activity differentiate to a mineral-secreting phenotype much like normal MC3T3-E1 cells. Osteopontin induction is dependent on at least a minimal level of alkaline phosphatase activity, although it is not dependent on induction of alkaline phosphatase at the RNA level. If alkaline phosphatase is supplied exogenously, osteopontin expression can be induced in conditions in which endogenous alkaline phosphatase is severely repressed.

The noniterative method of solving elliptic partial differential equations (pde) is generalized from two dimensions (2d) to (3d). The generalization is done in the context of electron (or positron)‐molecule scattering. Application is being made to e−N2, where, beyond the (r, θ) coordinates of the scattered electron, the third dimension is the internuclear dependence (R) of the target molecule. A brief derivation of the 3d‐pde is given together with its boundary conditions. In a separate section the numerical methodology, whereby the corresponding difference equation is solved, is outlined; this represents the core of the paper. However, more generally the chief significance of the 3d methodology is likely to be in e±‐scattering from polyatomic molecules in the fixed‐nuclei approximation. The essential features of this 3d‐pde and its solution are sketched.

Simplified multicomponent diffusion calculations are made for the ionosphere of Venus. Large differences in temperature between electrons and ions and appreciable temperature gradients that are near those of recent measurements are used. Compositions for which binary thermal diffusion coefficients for ions are the same as multi-ion ones are examined as well as those that are quite different. An attempt to combine binary coefficients to give multi-ions ones has not been particularly successful.

The adiabatic-nuclei approximation of vibrational-rotational excitation of homonuclear diatomic molecules can be simply augmented to describe the vibrational-rotational coupling by including the dependence of the vibrational wave function on j. Appropriate formulas are given, and the theory, is applied to e-H2 excitation, whereby it is shown that deviations from the simple Born-Oppenheimer approximation measured by Wong and Schultz can be explained. More important, it can be seen that the inclusion of the j-dependent centrifugal term is essential for transitions involving high-rotational quantum numbers.

The method of polarized orbitals is modified to treat low-energy scattering of electrons from highly polarizable systems, specifically alkali-metal atoms. The modification is carried out in the particular context of the $e$-Li system, but the procedure is general; it consists of modifying the polarized orbital, so that when used in the otherwise orthodox form of the method, it gives (i) the correct electron affinity of the negative ion (in this case ${\mathrm{Li}}^{\ensuremath{-}}$), (ii) the proper (i.e., Levinson-Swan) number of nodes of the associated zero-energy scattering orbital, and (iii) the correct polarizability. A procedure is devised whereby the scattering length can be calculated from the (known) electron affinity without solving the bound-state equation. Using this procedure we adduce a $^{1}S$ scattering length of $8.69{a}_{0}$. (The $^{3}S$ scattering length is $\ensuremath{-}9.22{a}_{0}$.) The above modifications can also be carried out in the (lesser) exchange adiabatic approximation. However, they lead to qualitatively incorrect $^{3}S$ phase shifts. The modified polarized-orbital phase shifts are qualitatively similar to close-coupling and elaborate variational calculations. Quantitative differences from the latter calculations, however, remain; they are manifested most noticeably in the very-low-energy total and differential spin-flip cross sections.

Using plane waves to describe the incident and scattered particles, and screened hydrogenic and Coulomb functions to describe the atomic electrons before and after ejections, we have calculated the differential and total ionization cross sections of 11 atoms and one ion by electron impact, and ionization of helium by proton impact. The effective charges of the screened hydrogenic functions are fixed by the Hartree-Fock calculations. Calculations have been carried out for the atomic s, p, and d electrons. For low atomic numbers, we find reasonable agreement with the experimental data. For intermediate atomic numbers, we expect our results to overestimate the actual cross sections, since our choice of a unit charge for the Coulomb function of the ejected electrons will overestimate the atomic dipole potential strength, and in turn the high-energy cross sections. The advantage of the method presented here is that the ionization amplitude is given in analytic form. This may allow further analysis on this amplitude, and facilitates extension of the numerical integration for the cross section to high impact energies.

n/a

n/a