Three-Dimensional Modeling of the Casing Effect in Onshore Controlled-Source Electromagnetic Surveys

The presence of steel-cased wells and other infrastructure causes a significant change in the electromagnetic fields that has to be taken into consideration in modeling and interpretation of field data. A realistic and accurate simulation requires the borehole casing to be incorporated into the modeling scheme, which is numerically challenging. Due to the huge conductivity contrast between the casing and surrounding media, a spatial discretization that provides accurate results at different spatial scales ranging from millimeters to hundreds of meters is required. In this paper, we present a full 3D frequency-domain electromagnetic modeling based on a parallel finite-difference algorithm considering the casing effect and investigate its applicability on the borehole-to-surface configuration of the Hontomín CO₂ storage site. To guarantee a robust solution of linear systems with highly ill-conditioned matrices caused by huge conductivity contrasts and multiple spatial scales in the model, we employ direct sparse solvers. Different scenarios are simulated in order to study the influence of the source position, conductivity model, and the effect of the steel casing on the measured data. Several approximations of the real hollow casing that allow for a large reduction in the number of elements in the resulting meshes are studied. A good agreement between the modeled responses and the real field data demonstrates the feasibility of simulating casing effects in complex geological areas. The steel casing of the well greatly increases the amplitude of the surface electromagnetic fields and thus improves the signal-to-noise ratio and the sensitivity to deep targets.     

Optical water type discrimination and tuning remote sensing band-ratio algorithms: Application to retrieval of chlorophyll and Kd(490) in the Irish and Celtic Seas

This paper assesses the feasibility of applying remote sensing algorithms based on blue/green reflectance ratios to Case 2 waters. Two algorithms from the SeaDAS (NASA) image processing package, OC4v4 for surface chlorophyll concentration, Chl, and K(490) for the attenuation coefficient for downward irradiance at 490 nm, K_d(490), were investigated using an extensive set of observations from the Irish and Celtic Seas. In situ data from a profiling radiometer were used as inputs for the algorithms to avoid uncertainties in atmospheric correction procedures, and direct measurements of Chl and K_d490 and were used for validation purposes. The standard versions of the algorithms performed poorly: OC4v4 generally overestimated Chl (with a very low coefficient of determination), and K(490) progressively underestimated K_d490 for values greater than 0.3 m⁻¹. A two-step procedure for level 2 product generation was therefore devised in which the numerical coefficients of OC4v4 and K(490) were tuned for the two optical water types known to occur most frequently in this region (McKee, D., Cunningham, A., 2006. Identification and characterisation of two optical water types in the Irish Sea from in situ inherent optical properties and seawater constituents. Estuarine, Coastal and Shelf Science 68, 305–316) by statistical regression on a data set of 102 stations from the Irish and Celtic Seas. The water types were distinguished by the magnitude of their normalised water leaving radiance signals at 665 nm, nL_w(665), and appropriate versions of the tuned algorithms applied to each water type. When this procedure was tested on an independent data set of 19 stations from the Bristol Channel, Chl values were recovered with an RMS error of 0.36 mg m⁻³ and K_d(490) values with an RMS error of 0.095 m⁻¹. The identification of water types from water-leaving radiance signals, and the application of band-ratio algorithms tuned for specific water types, may therefore provide a simple means of improving the quality of remote sensing products in optically complex shelf seas.     

Woodleigh impact structure,Australia: Shock petrography and geochemical studies

Abstract— The large, complex Woodleigh structure in the Carnarvon basin of Western Australia has recently been added to the terrestrial impact crater record. Many aspects of this structure are, however, still uncertain. This work provides a detailed petrographic assessment of a suite of representative drill core samples from the borehole Woodleigh 1 that penetrated uplifted basement rocks of the central part of this structure. Fundamental rock and mineral deformation data and high‐precision chemical data, including results of PGE and oxygen isotopic analysis, are presented. The sampled interval displays likely impact‐produced macrodeformation in the form of fracturing and breccia veining at the microscopic scale. Contrary to earlier reports that these breccias represent pseudotachylite (friction melt) or even shock/shear‐produced pseudotachylitic melt breccia cannot be confirmed due to pervasive post‐impact alteration. Abundant planar deformation features (PDFs) in quartz, in addition to diaplectic glass and partial isotropization, are the main shock deformation effects observed, confirming that Woodleigh is of impact origin. Over the investigated depth interval, the statistics of quartz grains with a variable number of sets of PDFs does not change significantly, and the patterns of crystallographic orientations of PDFs in randomly selected quartz grains does not indicate a change in absolute shock pressure with depth either. The value of oxygen isotopes for the recognition of meteoritic contamination, as proposed by earlier Woodleigh workers, is critically assessed. Neither INA nor PGE analyses of our samples support the presence of a meteoritic component within this basement section, as had been claimed in earlier work.     

DUN AND BRADSTREET'S MARKET IDENTIFIERS AND ECONOMIC GEOGRAPHY*

An increasing number of studies of business start-ups, closings, and relocations rely on a single data base, DUNS Market Identifiers. This data file to U.S. business establishments provides researchers with essential information about individual firms and establishments. It also exhibits serveral idiosyncracies and requires careful consideration of the meaning of such basic terms as “location”and “firm.”Drawing from a recent survey of 500 New Jersey firms. I consider some of the methodological amd conceptual issues Dun and Bradstreet data pose for industrial location analysts. The Dun and Bradstreet listings exaggerate the exaggerate the incidence and magnitude of establishment openings and closing.     

Late Quaternary glaciation and equilibrium-line altitudes of the Mayan Ice Cap, Guatemala, Central America

The Sierra los Cuchumatanes (3837 m), Guatemala, supported a plateau ice cap and valley glaciers around Montaña San Juan (3784 m) that totaled ∼ 43 km² in area during the last local glacial maximum. Former ice limits are defined by sharp-crested lateral and terminal moraines that extend to elevations of ∼ 3450 m along the ice cap margin, and to ca. 3000-3300 m for the valley glaciers. Equilibrium-line altitudes (ELAs) estimated using the area-altitude balance ratio method for the maximum late Quaternary glaciation reached as low as 3470 m for the valley glaciers and 3670 m for the Mayan Ice Cap. Relative to the modern altitude of the 0°C isotherm of ∼ 4840 m, we determined ELA depressions of 1110-1436 m. If interpreted in terms of a depression of the freezing level during maximal glaciation along the modern lapse rate of − 5.3°C km^− 1, this ΔELA indicates tropical highland cooling of ∼ 5.9 to 7.6 ± 1.2°C. Our data support greater glacial highland cooling than at sea level, implying a high tropical sensitivity to global climate changes. The large magnitude of ELA depression in Guatemala may have been partially forced by enhanced wetness associated with southward excursions of the boreal winter polar air mass.     

Modelling particle residence times in agricultural river basins using a sediment budget model and fallout radionuclide tracers

Contemporary patterns in river basin sediment dynamics have been widely investigated but the timescales associated with current sediment delivery processes have received much less attention. Furthermore, no studies have quantified the effect of recent land use change on the residence or travel times of sediment transported through river basins. Such information is crucial for understanding contemporary river basin function and responses to natural and anthropogenic disturbances or management interventions. To address this need, we adopt a process‐based modelling approach to quantify changes in spatial patterns and residence times of suspended sediment in response to recent agricultural land cover change. The sediment budget model SedNet was coupled with a mass balance model of particle residence times based on atmospheric and fluvial fluxes of three fallout radionuclide tracers (⁷Be, excess ²¹⁰Pb and ¹³⁷Cs). Mean annual fluxes of suspended sediment were simulated in seven river basins (38–920 km²) in south‐west England for three land cover surveys (1990, 2000 and 2007). Suspended sediment flux increased across the basins from 0.5–15 to 1.4–37 kt y^‐1 in response to increasing arable land area between consecutive surveys. The residence time model divided basins into slow (upper surface soil) and rapid (river channel and connected hillslope sediment source area) transport compartments. Estimated theoretical residence times in the slow compartment decreased from 13–48 to 5.6–14 ky with the increase in basin sediment exports. In contrast, the short residence times for the rapid compartment increased from 185–256 to 260–368 d as the modelled connected source area expanded with increasing sediment supply from more arable land. The increase in sediment residence time was considered to correspond to longer sediment travel distances linked to larger connected source areas. This novel coupled modelling approach provides unique insight into river basin responses to recent environmental change not otherwise available from conventional measurement techniques. © 2014 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Laplace-transform analytic element solution of transient flow in porous media

A Laplace-transform analytic element method (LT-AEM) is described for the solution of transient flow problems in porous media. Following Laplace transformation of the original flow problem, the analytic element method (AEM) is used to solve the resultant time-independent modified Helmholtz equation, and the solution is inverted numerically back into the time domain. The solution is entirely general, retaining the mathematical elegance and computational efficiency of the AEM while being amenable to parallel computation. It is especially well suited for problems in which a solution is required at a limited number of points in space–time, and for problems involving materials with sharply contrasting hydraulic properties. We illustrate the LT-AEM on transient flow through a uniform confined aquifer with a circular inclusion of contrasting hydraulic conductivity and specific storage. Our results compare well with published analytical solutions in the special case of radial flow.     

Experimental validation of the wavefield transform of electromagnetic fields

The wavefield transform is a mathematical technique for transforming low-frequency electromagnetic (EM) signals to a non-diffusive wave domain. The ray approximation is valid in the transform space and this makes traveltime tomography for 3D mapping of the electrical conductivity distribution in the subsurface possible. The transform, however, imposes stringent frequency bandwidth and signal-to-noise ratio requirements on the data. Here we discuss a laboratory scale experiment designed to collect transform quality EM data, and to demonstrate the practical feasibility of transforming these data to the wavefield domain.
We have used the scalable nature of EM fields to design a time-domain experiment using graphite blocks to simulate realistic field conditions while leaving the time scale undisturbed. The spatial dimensions have been scaled down by a factor of a thousand by scaling conductivity up by a factor of a million. The graphite blocks have two holes drilled into them to carry out cross-well and borehole-to-surface experiments. Steel sheets have been inserted between the blocks to simulate a conductive layer.
Our experiments show that accurate EM data can be recorded on a laboratory scale model even when the scaling of some features, such as drill-hole diameters, is not maintained. More importantly, the time-domain EM data recorded in cross-well and surface-to-borehole modes can be usefully and accurately transformed to the wavefield domain. The observed wavefield propagation delay is proportional to the direct distance between the transmitter and receiver in a homogeneous medium. In a layered medium, data accuracy is reduced and, hence, our results are not so conclusive. On the basis of the experimental results we conclude that the wavefield transform could constitute a valid approach to the interpretation of accurate, undistorted time-domain data if further improvement in the transform can be realized.