This Week's Lecture | Current Lecture Schedule | Fall 2001 Abstracts | Spring 2001 Abstracts | Fall 2000 Abstracts
FALL 2001 HEILAND LECTURE SERIES
Abstracts
Ernie Hull
M.S. Candidate
Colorado School of Mines

Friday, November 30, 2001 • Brown Building, Room 204A • 4:00 p.m.
3D Seismic Interpretation of Tectonics and Sedimentation within the Upper Cretaceous Lewis Shale and Fox Hills Sandstone Baggs, Wyoming
The purpose of this research is to sue 3-D seismic data to examine sedimentation within the Lewis Shale and Fox Hills Sandstone of the southern Washakie and northerm Sand Wash basins near Baggs, Wyoming. The study area is dominated by left lateral movement on a wrench fault interpreted to be the western extension of the Cheyenne Belt, a Precambrian suture. The eastern portion fo the study area contains the location of the South Baggs and West Side Cnal gas fields, a significant gas field in the eastern Green River Basin. The two fields are separated by a northwest trending fault within the student area.

Subsurface mapping indicates that tectonics played an active role during deposition of the Lower and Middle Lewis (Dad member) sediments. During deposition of the Upper Lewis and Fox Hills Sandstone, tectonics ceased to affect sedimentation and a pattern of differential compaction exists with respect to the underlying Dad member. Clinoforms within the Dad member prograded from the southwest to the northeast. Faults strike in three primary directions, northwest, northeast and east west. These strike directions fit well with left lateral wrench fault patterns predicted by published outcrop, subsurface and clay model studies. Folds and dip directions occur in orientations consistent with such models for left lateral wrench faults.

Ernie received his B.S. degree in Applied Geophysics from the University of California, Los Angeles, in June 1995. He worked for Western Geophysical from 1995-1999 and held various positions in Europe and the former Soviet Union. ernie is a M.S. candidate and successfully defended his thesis earlier this month. He will study in russia during the Spring 2002 and will join British Petroleum in Anchorage, Alaska, shortly thereafter.

Return to Top
Carlos Moita
M.S. Candidate
Colorado School of Mines
Friday, November 2, 2001 • Brown Building, Room 204A • 4:00 p.m.
Formation Invasion from a Horizontal Borehole
When a borehole is drilled, the drilling fluids used to remove drill cuttings and to counter balance reservoir pressure, invade permeable formations. Mud filtrate invasion changes the properties of the formation surrounding the borehole and thus, the measurements of borehole logging tools. Since, most borehole logging tool responses will be affected by the flushed zone, it is important to know the volume and shape of the mud filtrate invaded region. In order to investigate the feasibility of using a Ground-Penetrating Radar (GPR) tool to image the invaded volume, laboratory experiments were carried out by injecting fluids in a highly permeable sand body. For this purpose a transparent Plexiglas? box was built and filled with 100% fresh-water saturated sand, which allowed us to visually follow the path of dyed mineral oil, injected through a small horizontal borehole. Two GPR systems were used. A 900 MHz pulse EKKO 1000? radar system to record from the outside of the Plexiglas box, and a 1500 MHz borehole directional radar tool (DBOR) to record invasion front reflections from the inside of the borehole. After processing the GPR data, we found that the interpretation of the invasion fronts in the radar reflection data, despite the intense horizontal fingering, agreed very well with the same fronts traced on the front and back of the sandbox. Furthermore, the estimate for the irreducible water saturation and relative permeability's, using the radar image of the invasion shape matched the estimates based on of the invasion front traced on the physical model. It is concluded that GPR is an effective tool to assess the shape and extent of the invaded zone, provided that no highly conductive fluids surround the directional radar tool.

Carlos Moita, Has an MS in Geology from the University of Lisbon in 1996. I worked for the "Instituto Geologico e Mineiro" in Portugal since 1994 to present. And I'm a candidate to a Master degree in Geophysics at CSM.

Return to Top
Daniel T. Georgi
Director of Science Baker Hughes / Inteq
Friday, October 26, 2001• Metals Hall, Green Center • 4:00 p.m.
Petrophysical Aspects of Resistivity and Permeability Anisotropy
Cost effective, efficient production of hydrocarbons depends on accurate reserves estimates, reservoir architecture and hydrocarbon distribution. In addition complete reservoir, description including both horizontal and vertical permeability is critical to the development strategy. Hydrocarbon recovery efficiency depends on many factors but one key parameter is the ratio of the vertical to horizontal permeability. There are numerous ways to estimate permeability. It is now even possible to estimate permeability from continuous formation evaluation measurements. However, it remains difficult to determine directional permeability. Nuclear Magnetic Resonance (NMR) logs are used routinely to estimate permeability, but the NMR-derived permeability is based on scalar properties (e.g. ??porosity, T2 decay time distributions, bulk volume irreducible water) and is inherently a scalar property. Few Formation Evaluation measurements provide directional information. Dip and image logs provide bed thickness and layer dip information. However, only the new multi-component induction (3DEXO) and crossed-dipole shear-wave acoustic tools provide direct measurements of macroscopic formation anisotropy.ium parameters.

In this presentation, I will discuss the new 3-D Induction logging tool (3DEX) and then address both theoretically and with real data the computation of permeability anisotropy. I use both laminated sand models as well as macroscopic models based on the resistivity anisotropy measurements to estimate vertical over horizontal (kv:kh) permeability ratios. In addition to multi-component induction data, I will explore both stand alone and joint interpretation of NMR property variations to predict the macroscopic reservoir permeability anisotropy.

Dr. Daniel T. Georgi (dan.georgi@bakerhughes.com), Director of Science, BakerAtlas/INTEQ, holds an M.A. in Geophysics and a Ph.D. in Earth Sciences from Columbia University. He has been involved with fundamental research, tool development, and field studies throughout his career. He has dealt with many aspects of formation evaluation of conventional, fractured, and heavy-oil reservoirs as well as open and cased hole log analysis.

Return to Top
Petr Jilek
Ph.D. Candidate
CSM Department of Geophysics
Friday, October 19, 2001• Metals Hall, Green Center • 4:00 p.m.
Using PP- and PS-wave Reflection Coefficients for Quantitative AVO Analysis in Anisotropic Media

Over the last two decades, amplitude variation with offset (AVO) analysis has been used for direct hydrocarbon detection, lithology identification and, more recently, for reservoir monitoring. The improved acquisition design and higher quality of seismic data have allowed the AVO methodology to become more quantitative and sometimes provide us with estimates of important reservoir characteristics including porosity and crack density. Quantitative AVO inversion, however, is difficult to carry out due to a number of practical and theoretical problems, such as the instability of the conventional AVO inversion of P-wave data. Consequently, AVO analysis is often incapable of constraining the individual physical parameters needed for reservoir characterization, particularly in anisotropic media. We suggest to stabilize the AVO inversion by using PP-and PS-waves collected in multicomponent surveys for a wide range of source-receiver azimuths (i.e., using 3D-3C data). Converted PS-waves are generated simultaneously with PP-waves and, therefore, represent a relatively inexpensive source of additional information. Joint inversion of PP- and PS-wave reflection coefficients for the parameters of azimuthally anisotropic formations is the main focus of this work. The developed inversion algorithms are designed not only to estimate the medium parameters, but also to provide simple analytic insight into the parameter combinations constrained by the available data. The analytic aspect of the inversion is especially important in anisotropic media where "blind" numerical algorithms often suffer from non-uniqueness and trade-offs between medium parameters.

Received MS in geophysics from the Charles University, Prague, Czech Republic. Currently pursuing Ph.D. degree in geophysics, CSM.

Return to Top
Kristen W. Sneddon
M.S. Candidate
CSM Department of Geophysics
Friday, October 19, 2001• Metals Hall, Green Center • 4:00 p.m.
Modeling Porosity, DNAPL Saturation, and
Layer Thickness from Non-Invasive GPR
Dense non-aqueous phase liquids (DNAPL) are a pervasive and persistent category of groundwater contamination. In an effort to better understand their unique subsurface behavior, a controlled and carefully monitored injection of PCE (perchloroethylene), a typical DNAPL, was performed in conjunction with the University of Waterloo at Canadian Forces Base Borden in 1991. Of the various geophysical methods used to monitor the migration of injected PCE, the U.S. Geological Survey collected 500-MHz ground penetrating radar (GPR) data. GPR data were acquired over time on a fixed two-dimensional surficial grid as the DNAPL was injected into the subsurface. Using these data, interactive full-waveform GPR modeling of regularized field traces resolves relative dielectric permittivity versus depth profiles for pre-injection and later-time data. Modeled values are end members in recursive calculations of the Bruggeman-Hanai-Sen (BHS) mixing formula, yielding interpreted pre-injection porosity, post-injection DNAPL saturation values, and layer thickness.

Kristen received her B.S. degree in Applied Geophysics from the University of California at Los Angeles (UCLA) in 1996. She currently works at the U.S. Geological Survey in Denver with the Crustal Imaging and Characterization Team (CICT) and hopes to complete her M.S. degree in Geophysics from CSM in the spring. Her thesis advisor is Gary Olhoeft.
Return to Top
Ryan E. North
M.S. Candidate
CSM Department of Geophysics
Friday, October 12, 2001• Metals Hall, Green Center • 4:00 p.m.
Numerical modeling of acoustic wave propagation through a shock tube
In recent years experimental groups at Delft University and Colorado School of Mines have experimentally investigated the acoustic properties of porous and permeable media using a shock tube geometry. The water- saturated sample is placed within one end of a presumably rigid cylinder; an impulsive plane wave insonifies the sample from the other end. Pressure gauges imbedded in the wall monitor the fluid pressures in the small gap between the sample and the wall. These pressures have, to date, been assumed to be the same as in the adjacent fluid in the sample, and been interpreted within the context of one-dimensional wave propagation according to the Biot theory, which did not take the gap into account.

Although the 1-D code successfully predicted travel times of various waveforms it failed to predict amplitudes correctly. For this reason we used a code developed by the Keldysh Institute of Applied Mathematics (KIAM) to model the radial and axial variations inherent in the problem. The finite difference time-stepping code utilizes the full Biot theory, including the frequency dependent tortuosity and permeability. The results of the present study indicate that a one-dimensional theory may suffice, but the character of the normal modes may differ significantly from those of a plane wave in the bulk, due to the presence of the small gap. The code allows us to determine whether the Biot slow wave exists for differing rock types.

B.S. Engineering Physics (1997) West Virginia Wesleyan College Future Degree (2001)

Return to Top
Philip J. Brown II
M.S. Candidate
CSM Department of Geophysics
Friday, October 12, 2001• Metals Hall, Green Center • 4:00 p.m.
The Observation of the Biot Slow-P Wave in Natural Rocks and the Determination of Petrophysical Parameters Using a Shock Tube
Modeling of acoustic wave propagation in rocks is critical in seismic interpretation, acoustic log analysis, and rock and fluid property evaluations. Maurice Biot's relations are among the most commonly applied equations to predict fluid and frequency effects.

Biot's analysis of wave propagation predicts that a second compressional, or 'slow' P-wave can occur in fluid saturated porous media. This wave has been detected directly in natural rocks using a shock tube. Since the observation of this wave mode in sedimentary rocks, many questions have arisen. These include, when does the slow P-wave need to be taking into account when interpreting acoustic data and can the behavior of the slow P-wave be used to ascertain petrophysical parameters of a rock?

In order to help answer these questions, shock tube experiments were performed. Acoustic travel times and amplitudes of fast and slow compressional waves, as well as reflection coefficients from naturally occurring sandstone samples were determined. Sandstones studied include cores from the Lyons, Fox Hills, Berea, and Bentheimer formations. These samples had a permeability range of 0.3 milliDarcies to 5 Darcies, and a porosity range of 9 to 30 percent. Empirical data collected on wave velocities, attenuation and pore pressure were compared to a one-dimensional Biot model which demonstrated that the model accurately predicts wave velocity but overestimates amplitudes. Current observations show that the slow-compressional wave cannot be detected using the present experimental configuration in samples with permeabilities below 200 mD. Shock Tube data has shown that the reflection coefficients at the water to sample interface decreased sharply when the slow P-wave was present.

Philip Brown received a B.S. in Geophysical Engineering from the Colorado School of Mines in 1994. After graduation, Phil was employed by Sperry-Sun Drilling Services working as a Formation Evaluation while Drilling Engineer offshore in the Gulf of Mexico until 1996. Since 1992, Phil has also participated in a number of projects involving near-surface geophysical survey design, data collection, processing and interpretation for hydrological, environmental, engineering, mining, and forensic applications. Such projects included work for The Colorado School of Mines, Necro Search International, Pegasus Gold Corporation, MicroGeophysics Corporation, Blackhawk Geosciences, and Infraseis Incorporated. In 1998, Phil joined Energy Data Services in Lakewood, CO working on well log data clean up, well log normalization, formation evaluation and research. Since August of 1999, Phil has been a full time student at the Colorado School of Mines with the Center for Petrophysics where he worked as the Geophysics Department Webmaster and held the Treasurer position with the Society of Geophysics Graduate Students (SGGS). Phil is currently working for the Crustal Imaging and Characterization Team at the United States Geologic Survey in Lakewood CO as a Research Geophysicist and is Vice President of the SGGS. Further information about Phil can be found on his web site, http://www.mines.edu/students/p/pjbrown/.

Return to Top
Dr. Biondo Biondi
Associate Professor (Research)
Co-director Stanford Exploration Project
Friday, October 5, 2001• Metals Hall, Green Center • 4:00 p.m.
Seismic Imaging beyond Kirchhoff: Challenges and Opportunities
Kirchhoff migration has theoretical and practical limitations when multipathing and other complex wave phenomena occur, such as in subsalt exploration. Wavefield extrapolation methods have the potential to overcome these limitations, but we encounter several new challenges when applying wavefields methods to 3-D prestack imaging.

At the Stanford Exploration Project we have made substantial progress to make wavefield methods competitive with Kirchhoff methods, but we also found new and stimulating problems. I will present our results and discuss some of the unsolved problems related to: computational cost, velocity estimation, amplitude-preserved migration, data regularization, and poor reflector illumination. New insights in seismic and inversion theory are required to harness the foreseeable increase in computer power for imaging reservoir in complex areas.

Biondo Biondi is Associate Professor (Research) of Geophysics at Stanford University and leads the Stanford Exploration Project efforts in 3-D imaging. Biondo graduated from Politecnico di Milano in 1984 and received a M.S. (1988) and a Ph.D. (1990) in geophysics from Stanford. He has made contributions on several aspects of seismic imaging, ranging from velocity estimation to parallel algorithms for seismic migration. Since the early nineties he has been at the forefront of the development of wave-equation 3-D prestack migration methods. Biondo is in close contact with the practical applications of seismic imaging by his involvement with 3DGeo Development, that he founded in 1994. 3DGeo Development is a growing seismic services start-up that brings innovative technologies to the exploration Industry such as wave-equation imaging and Internet based seismic processing. Biondo is a member of SEG and EAGE and is a member of the SEG Research Committee and SEG Translation Committee.

Return to Top
Ed Jenner
Ph.D. Candidate
Friday, September 7, 2001• Metals Hall, Green Center • 4:00 p.m.
A method for azimuthal velocity analysis and application to 3D data
The azimuthal variation of normal-moveout velocity (ANMO) for pure modes is described by an ellipse, even if the medium is arbitrarily anisotropic and heterogeneous. The simplest anisotropic model that produces azimuthal velocity variation is that of an isotropic, heterogeneous medium that has vertical, aligned fractures (HTI). It should be noted, however, that heterogeneity (e.g. reflector dip or lateral velocity variation over a horizontal reflector) can also cause ANMO, even in an isotropic medium. The HTI model is particularly important since fractures are commonly observed in reservoir rocks where, if connected, may form paths of higher permeability compared to the remaining pores. Analyses of fracture intensity and direction are thus particularly important as a reservoir characterization tool since the fractures can dominate fluid flow during production. In this talk I will discuss a method for azimuthal velocity analysis which does not rely on picking semblance maxima on azimuth sectored data. Instead the method uses a flattening procedure that is akin to a time-variant trim static. A least squares fit to the traveltimes gives the fast and slow velocity magnitudes and direction. I will show examples of this method applied to several 3D land datasets and discuss the implications for future processing of both land and marine data. Ed received his Bachelors degree in Physics with Astrophysics from the University of Birmingham, U.K. and his Masters in Geophysics from the University of Leeds, U.K. He has worked for BP on multiple removal, CONOCO on prestack migration and AXIS Geophysics on azimuthal processing as well as interpretation of North Slope data with ARCO Alaska. In 1999 Ed spent a year with CGG in Paris working on developing a methodology for converted wave prestack depth migration. Ed is currently with the Reservoir Characterization Project and expects to complete his Ph.D. by the end of the year.
Return to Top