Geatches D.,Energy Resources Engineering |
McCarty D.,Chevron |
Wilcox J.,Energy Resources Engineering
American Mineralogist | Year: 2014
The focus of the study presented here is the illite component of the clay mineral components found within gas shale, specifically the creation of robust, atomistic models of illite particles bearing characteristics of the 2M1 polytype. The template for the illite models is derived from crystallographic data, which is used together with the general formula, K1.67(Al3.29Fe0.38Mg0.32)[(Si6.68Al1.32)O20(OH)4] to create different particles of illite, each containing potassium plus either ammonium (NH4 +) or hydronium (H3O+) ions, or water molecules in two different proportions. These atomistic models are optimized from first principles using plane waves and pseudo-potentials within the formalism of density functional theory (DFT). The resulting lattice lengths, lattice dynamics in the form of infrared frequencies, positions of interlayer molecules and simulated Xâray powder diffraction (XRPD) patterns are analyzed with reference to the available experimental data. We conclude that all four illite models display properties that fall within the range of available experimental data, and that therefore these are state-of-the-art atomistic illites, ready for use in further studies. Furthermore, the analytical data from these models will enable the characterization of physical samples of illite with varying interlayer constituents. © 2014 by Walter de Gruyter Berlin/Boston.
Having an accurate picture of the planet's subsurface is crucial for rational decision-making across a wide variety of activities, including environmental cleanup, oil and natural gas well-drilling, and the underground sequestration of greenhouse gases that contribute to climate change. But oftentimes, it's impossible to know all of the details about a subsurface system. It's therefore important to be able to systematically calculate the myriad unknown configurations the system might take – in other words, the uncertainty associated with that system – and then collect data to rule out unlikely scenarios. "If you're using very simple models that represent the subsurface as simply layered, or homogenous, then you could be in for a rude shock when you start drilling for oil or water, or injecting cleanup chemicals into a contaminated aquifer," said geostatistician Jef Caers, a professor of geological sciences at Stanford School of Earth, Energy & Environmental Sciences. "The wells you drill could be dry and the cleanup project that you thought would take a few months could take years. These things have actually happened, and it was due to faulty assumptions about uncertainties," Caers said. The solution is to work from more realistic uncertainty models of the subsurface, but studying the deep Earth is difficult. The steady buildup and erosion of sediments that leads to the creation of complex structures and layers in the subsurface takes place over thousands or hundred of thousands of years and is not easily observable. "Ideally, we would slice a riverbed to see its vertical profile, but in reality you can't do that," said Céline Scheidt, a senior research engineer in the Department of Energy Resources Engineering. Instead, scientists use a combination of real-world and computer-generated models to replicate the physical process or system that is of interest – for example, the deposition of sediments at the mouth of a river or the complex accumulation of rock layers. They then use geostatistical tools such as multiple-point statistics (MPS) to quantify uncertainty about the subsurface. "We want to fundamentally understand uncertainty in naturally variable systems," said Caers, who recently published a book on MPS. "For example, in a deltaic environment, you want to know all the possible deltaic configurations you can possibly imagine. Then, when you get subsurface data about the site, you can rule out those configurations that are not consistent with the data." MPS algorithms have gained popularity in the past decade due to their ability to produce geologically realistic representations of actual variability. MPS works by using "training images" to narrow down the possible permutations that a system can take. Training images contain numerical patterns that are deemed representative of a particular geologic system. Typically, only one training image is used, but for situations where large uncertainty is present, this is not sufficient. In a new study, Scheidt, Caers and their colleagues devised a new method for selecting multiple training images that can be plugged into MPS algorithms to help simulate natural variability. The first step involved creating a miniature river basin in a lab and running an experiment that was designed to replicate the dynamic deposition of sediments in a delta. The experiment, conducted by Anjali Fernandes and Chris Paola, a postdoctoral researcher at Tulane University and a sedimentologist at the University of Minnesota, respectively, was designed to simulate tens of thousands of years' worth of river deposition in a few days. An overhead camera took snapshots of the experiment at one-minute intervals, and the team selected 136 images – representing about two hours – for analysis. Once the team had a pool of images from which to select training images, the questions then became, how should the training images be selected? And how many images are required to capture the variations in water flow and sediment deposition in the experiment? "By examining the images, it became clear that more than one training image was needed," Scheidt said, "but using all of the images as training images would be impractical and unnecessary, as most likely not all of the images provide new information." The solution to their dilemma came from an unusual source. Lewis Li, a graduate student in Caers' lab, mentioned the problem to his roommate, a neuroscientist, who told Li about the "demon algorithm," which doctors use to gauge the efficacy of cancer treatments. "The demon algorithm calculates the strain acting on a moving image over time," Caers said. "Doctors use it to determine whether a tumor's shape change is due to radiation therapy by seeing where major morphological changes to the tumor are occurring and how those changes correlate with times of treatment." While the demon algorithm is primarily used to spot major changes in a system, Scheidt realized it could also be used to analyze a series of images and select those taken during periods when changes to a system are small and similar in nature. When the team applied the demon algorithm to their 136 snapshots, it chose six images that it deemed to be representative of the overall types of depositional changes occurring throughout the experiment. When the team used those six images as training images for the MPS algorithm, they found that the variability generated by their geostatistical method matched well with the "natural variability" represented by the overhead snapshots of the experimental basin. "This research is a first step toward bridging the fundamental gap that exists between those fields of science focusing on the physical understanding of natural systems and the statistical representation of uncertainty when modeling complex systems," Caers said. Explore further: A picture is worth 1000 words, but how many emotions?
Geatches D.L.,Energy Resources Engineering |
Wilcox J.,Energy Resources Engineering
European Journal of Mineralogy | Year: 2014
In this investigation we create a variety of dioctahedral, interlayer (cation)-deficient mica clay mineral models of the 1Millite series, using the ab initio method of density functional theory (DFT), employing plane waves, pseudo potentials, and super cells. The structures we create are modelled on an average formula of illite such as (Ca0.0590.059,K0.655) (Si3.597, Al0.403)(Fe0.628, Al 0.696,Mg0.428)O10(OH)2, as well as existing crystal structures and experimental data. This study investigates the relative positions of both octahedral and tetrahedral cation substitutions (e.g. Fe3+ for Al3+ and Al3+ for Si4+) in both trans-and cis-vacant structures. The purpose of this work is to create robust, representative models of single layers (I) and double layers (I-I) of illite found within gas shale, that can be used in future investigations. We examine lattice parameters, interatomic distances of tetrahedral and octahedral sheets, K-O distances, and K positions, and find reasonable agreement with experimental data, and excellent agreement with simulated data of similar structures. Finally, to enable experimental characterization of samples of illite with different compositions, we provide simulated X-ray powder diffraction (XRPD) patterns including up to the first ten (h k l) reflections, and find characteristic identifiers of trans-and cis-vacant micas in agreement with previous theoretical and experimental studies. We conclude that the models we create are feasible models of I and I-I layers of illite found within shale and are ready for use in future work. © 2013 E. Schweizerbart'sche Verlagsbuchhandlung, D-70176 Stuttgart.