Finite Difference Elastic Wave Modeling Including Surface Topography

Surface topography and the weathered zone (i.e., heterogeneity near the earth’s surface) have great effects on elastic wave propagation. Both surface waves and body waves are contaminated by scattering and conversion by the irregular surface topographic features. In this paper, we present a 2D numerical solver for the elastic wave equation that combines a 4th-order ADER scheme (Arbitrary high-order accuracy using DERivatives) with the characteristic variable method at the free surface boundary. The method is based on the velocity-stress formulation. We demonstrate the method by calculating synthetic seismograms for simple features

Predicting stress-induced anisotropy around a borehole

Formation elastic properties near a borehole may be altered from their original state due to the stress concentration around the borehole. This could result in a biased estimation of formation properties but could provide a means to estimate in situ stress from sonic logging data. In order to properly account for the formation property alteration, we propose an iterative numerical approach to calculate the stress-induced anisotropy around a borehole by combining Mavko’s rock physics model and a finite-element method. We show the validity and accuracy of our approach by comparing numerical results to laboratory measurements of the stress-strain relation of a sample of Berea sandstone, which contains a borehole and is subjected to uniaxial stress loading. Our iterative approach converges very fast and can be applied to calculate the spatially varying stiffness tensor of the formation around a borehole for any given stress state.Massachusetts Institute of Technology. Earth Resources Laboratory (Founding Member Postdoctoral Fellowship

Lithospheric processes

This is the final report of a three-year, Laboratory Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL). The main objective was to improve understanding of the origin and evolution of the Earth`s lithosphere by studying selected processes, such as deformation and magmatic intrusion during crustal extension, formation and extraction of mantle melts, fluid transport of heat and mass, and surface processes that respond to deep-seated events. Additional objectives were to promote and develop innovative techniques and to support relevant educational endeavors. Seismic studies suggest that underplating of crust by mantle melts is an important crustal-growth mechanism, that low-angle faults can be seismogenic, and that shear deformation creates mantle anisotropy near plate boundaries. Results of geochemical work determined that magmas from oceanic intraplate islands are derived from a uniform depth in the upper mantle, whereas melts erupted at mid-ocean ridges are mixed from a range of depths. The authors have determined the extent and style of fluid infiltration and trace-element distribution in natural magmatic systems, and, finally, investigated {sup 21}Ne as a tool for dating of surficial materials

Computational 3-D inversion for seismic exploration

This is the final report of a four-month, Laboratory Directed Research and Development (LDRD) project carried out at the Los Alamos National Laboratory (LANL). There is a great need for a new and effective technology with a wide scope of industrial applications to investigate media internal properties of which can be explored only from the backscattered data. The project was dedicated to the development of a three-dimensional computational inversion tool for seismic exploration. The new computational concept of the inversion algorithm was suggested. The goal of the project was to prove the concept and the practical validity of the algorithm for petroleum exploration

Bayesian Neural Networks for Geothermal Resource Assessment: Prediction with Uncertainty

We consider the application of machine learning to the evaluation of geothermal resource potential. A supervised learning problem is defined where maps of 10 geological and geophysical features within the state of Nevada, USA are used to define geothermal potential across a broad region. We have available a relatively small set of positive training sites (known resources or active power plants) and negative training sites (known drill sites with unsuitable geothermal conditions) and use these to constrain and optimize artificial neural networks for this classification task. The main objective is to predict the geothermal resource potential at unknown sites within a large geographic area where the defining features are known. These predictions could be used to target promising areas for further detailed investigations. We describe the evolution of our work from defining a specific neural network architecture to training and optimization trials. Upon analysis we expose the inevitable problems of model variability and resulting prediction uncertainty. Finally, to address these problems we apply the concept of Bayesian neural networks, a heuristic approach to regularization in network training, and make use of the practical interpretation of the formal uncertainty measures they provide.Comment: 27 pages, 12 figure

Antibodies to Chlamydia trachomatis in patients presenting with ectopic pregnancy at Groote Schuur hospital

Objectives. To detennine the prevalence of antibodies to Chlamydia trachomatis in women presenting with ectopic pregnancies to Groote Schuur Hospital.Methods. C. trachomatis antibody titres were measured using a modified micro-immunofluorescence test in women presenting with ectopic pregnancy. Control subjects were drawn from women with term pregnancies and an uneventful reproductive history.Results. Seventy-four patients and controls were studied. Demographic variables were controlled for at time of entry into the study. A significant association between the number of lifetime sexual partners and exposure to C. trachomatis was noted (P = 0.001). Patients with ectopic pregnancies had significantly higher antibody titres than control subjects (P = 0.001), and in both groups the prevalence of background antichlamydial antibody was high (ectopic pregnancies 59%, pregnant controls 32%).Conclusions. While the role of C. trachomatis infection in women who develop ectopic pregnancies needs to be explored further, it seems wise to treat them all with empirical antibiotics at the time of presentation

Prestack depth migration for complex 2D structure using phase-screen propagators

We present results for the phase-screen propagator method applied to prestack depth migration of the Marmousi synthetic data set. The data were migrated as individual common-shot records and the resulting partial images were superposed to obtain the final complete Image. Tests were performed to determine the minimum number of frequency components required to achieve the best quality image and this in turn provided estimates of the minimum computing time. Running on a single processor SUN SPARC Ultra I, high quality images were obtained in as little as 8.7 CPU hours and adequate images were obtained in as little as 4.4 CPU hours. Different methods were tested for choosing the reference velocity used for the background phase-shift operation and for defining the slowness perturbation screens. Although the depths of some of the steeply dipping, high-contrast features were shifted slightly the overall image quality was fairly insensitive to the choice of the reference velocity. Our jests show the phase-screen method to be a reliable and fast algorithm for imaging complex geologic structures, at least for complex 2D synthetic data where the velocity model is known

In vivo imaging of microenvironmental and anti-PD-L1-mediated dynamics in cancer using S100A8/S100A9 as an imaging biomarker

Purpose: As a promotor of tumor invasion and tumor microenvironment (TME) formation, the protein complex S100A8/S100A9 is associated with poor prognosis. Our aim was to further evaluate its origin and regulatory effects, and to establish an imaging biomarker for TME activity. Methods: S100A9−/−cells (ko) were created from syngeneic murine breast cancer 4T1 (high malignancy) and 67NR (low malignancy) wildtype (wt) cell lines and implanted into either female BALB/c wildtype or S100A9−/− mice (n = 10 each). Anti-S100A9-Cy5.5-targeted fluorescence reflectance imaging was performed at 0 h and 24 h after injection. Potential early changes of S100A9-presence under immune checkpoint inhibition (anti-PD-L1, n = 7 vs. rat IgG2b as isotype control, n = 3) were evaluated. Results: In S100A9−/−mice contrast-to-noise-ratios were significantly reduced for wt and S100A9−/−tumors. No significant differences were detected for 4T1 ko and 67NR ko cells as compared to wildtype cells. Under anti-PD-L1 treatment S100A9 presence significantly decreased compared with the control group. Conclusion: Our results confirm a secretion of S100A8/S100A9 by the TME, while tumor cells do not apparently release the protein. Under immune checkpoint inhibition S100A9-imaging reports an early decrease of TME activity. Therefore, S100A9-specific imaging may serve as an imaging biomarker for TME formation and activity

Hot dry rock fracture propagation and reservoir characterization

North America's largest hydraulic fracturing opeations have been conducted at Fenton hill, New mexico to creae hot dry rock geothermal reservoirs. Microearthquakes induced by these fracturing operations were measured with geophones. The large volume of rock over which the microearthquakes were distributed indicates a mechanism of hydraulic stimulation which is at odds with conventional fracturing theory, which predicts failure along a plane which is perpendicular to the least compressive earth stress. Shear slippage along pre-existing joints in the rock is more easily induced than conventional tensile failure, particularly when the difference between minimum and maximum earth stresses is large and the pre-existing joints are oriented at angles between 30 and 60)degree) to the principal earth stresses, and a low viscosity fluid like water is injected. Shear slippage results in local redistribution of stresses, which allows a branching, or dendritic, stimulation pattern to evolve, in agreement with the patterns of microearthquake locations. Field testing of HDR reservoirs at the Fenton Hill site shows that significant reservoir growth occurred as energy was extracted. Tracer, microseismic, and geochemical measurements provided the primary quantitative evidence for the increases in accessible reservoir volume and fractured rock surface area. These temporal increases indicate that augmentation of reservoir heat production capacity in hot dry rock system occurred. For future reservoir testing, Los Alamos is developing tracer techniques using reactive chemicals to track thermal fronts. Recent studies have focused on the kinetics of hydrolysis of derivatives of bromobenzene, which can be used in reservoirs as hot as 275)degree)C

Coda-wave interferometry analysis of time-lapse VSP data for monitoring geological carbon sequestration

Injection and movement/saturation of carbon dioxide (CO2) in a geological formation will cause changes in seismic velocities. We investigate the capability of coda-wave interferometry technique for estimating CO2-induced seismic velocity changes using time-lapse synthetic vertical seismic profiling (VSP) data and the field VSP datasets acquired for monitoring injected CO2 in a brine aquifer in Texas, USA. Synthetic VSP data are calculated using a finite-difference elastic-wave equation scheme and a layered model based on the elastic Marmousi model. A possible leakage scenario is simulated by introducing seismic velocity changes in a layer above the CO2 injection layer. We find that the leakage can be detected by the detection of a difference in seismograms recorded after the injection compared to those recorded before the injection at an earlier time in the seismogram than would be expected if there was no leakage. The absolute values of estimated mean velocity changes, from both synthetic and field VSP data, increase significantly for receiver positions approaching the top of a CO2 reservoir. Our results from field data suggest that the velocity changes caused by CO2 injection could be more than 10% and are consistent with results from a crosswell tomogram study. This study demonstrates that time-lapse VSP with coda-wave interferometry analysis can reliably and effectively monitor geological carbon sequestration