Multi‐dimensional analyses of the SEAM controlled source electromagnetic data—the story of a blind test of interpretation workflows

Using a subset of the SEG Advanced Modeling Program Phase I controlled‐source electromagnetic data, we apply our standard controlled‐source electromagnetic interpretation workflows to delineate a simulated hydrocarbon reservoir. Experience learned from characterizing such a complicated model offers us an opportunity to refine our workflows to achieve better interpretation quality. The exercise proceeded in a blind test style, where the interpreting geophysicists did not know the true resistivity model until the end of the project. Rather, the interpreters were provided a traditional controlled‐source electromagnetic data package, including electric field measurements, interpreted seismic horizons, and well log data. Based on petrophysical analysis, a background resistivity model was established first. Then, the interpreters started with feasibility studies to establish the recoverability of the prospect and carefully stepped through 1D, 2.5D, and 3D inversions with seismic and well log data integrated at each stage. A high‐resistivity zone is identified with 1D analysis and further characterized with 2.5D inversions. Its lateral distribution is confirmed with a 3D anisotropic inversion. The importance of integrating all available geophysical and petrophysical data to derive more accurate interpretation is demonstrated.     

Mineralogy and thermal properties of V-type Asteroid 956 Elisa: Evidence for diogenitic material from the Spitzer IRS (5-35 μm) spectrum

We present the thermal infrared (5-35 μm) spectrum of 956 Elisa as measured by the Spitzer Infrared Spectrograph (“IRS”; Houck, J.R. et al. [2004]. Astrophys. J. Suppl. 154, 18-24) together with new groundbased lightcurve data and near-IR spectra. From the visible lightcurve photometry, we determine a rotation period of 16.494 ± 0.001 h, identify the rotational phase of the Spitzer observations, and estimate the visible absolute magnitude (H_V) at that rotational phase to be 12.58 ± 0.04. From radiometric analysis of the thermal flux spectrum, we find that at the time of observation 956 Elisa had a projected radius of 5.3 ± 0.4 km with a visible albedo p_V = 0.142 ± 0.022, significantly lower than that of the prototype V-type asteroid, 4 Vesta. (This corresponds to a radius of 5.2 ± 0.4 km at lightcurve mean.) Analysis with the standard thermal model (STM) results in a sub-solar temperature of 292.3 ±  2.8 K and beaming parameter η = 1.16 ± 0.05. Thermophysical modeling places a lower limit of on the thermal inertia of the asteroid’s surface layer (if the surface is very smooth) but more likely values fall between 30 and depending on the sense of rotation.The emissivity spectrum, calculated by dividing the measured thermal flux spectrum by the modeled thermal continuum, exhibits mineralogically interpretable spectral features within the 9-12 μm reststrahlen band, the 15-16.5 μm Si-O-Si stretching region, and the 16-25 μm reststrahlen region that are consistent with pyroxene of diogenitic composition: extant diogenitic pyroxenes fall within the narrow compositional range Wo_2±1En_74±2Fs_24±1. Spectral deconvolution of the 9-12 μm reststrahlen features indicates that up to ≈20% olivine may also be present, suggesting an olivine-diogenite-like mineralogy. The mid-IR spectrum is inconsistent with non-cumulate eucrite as the major component on the surface of 956 Elisa, although cumulate eucrite material may be present at abundances lower than that of the diogenite component.Analysis of new near-IR spectra of 956 Elisa with the Modified Gaussian Model (MGM; Sunshine, J.M., Pieters, C.M., Pratt, S.F. [1990]. J. Geophys. Res. 95 (May), 6955-6966) results in two pyroxene compositions: 75% magnesian low-Ca pyroxene and 25% high-Ca pyroxene. High-Ca pyroxene is not evident in the mid-IR data, but may belong to a component that is underrepresented in the mid-IR spectrum either because of its spatial distribution on the asteroid or because of its particle size. High-Ca pyroxenes that occur as exsolution lamellae may also be more evident spectrally in the NIR than in the mid-IR. In any case, we find that the mid-IR spectrum of 956 Elisa is dominated by emission from material of diogenite-like composition, which has very rarely been observed among asteroids.     

Ultraviolet spectroscopy of Asteroid (4) Vesta

We report a comprehensive review of the UV–visible spectrum and rotational lightcurve of Vesta combining new observations by Hubble Space Telescope and Swift Gamma-ray Burst Observatory with archival International Ultraviolet Explorer observations. The geometric albedos of Vesta from 220 nm to 953 nm are derived by carefully comparing these observations from various instruments at different times and observing geometries. Vesta has a rotationally averaged geometric albedo of 0.09 at 250 nm, 0.14 at 300 nm, 0.26 at 373 nm, 0.38 at 673 nm, and 0.30 at 950 nm. The linear spectral slope as measured between 240 and 320 nm in the ultraviolet displays a sharp minimum near a sub-Earth longitude of 20°, and maximum in the eastern hemisphere. This is consistent with the longitudinal distribution of the spectral slope in the visible wavelength. The photometric uncertainty in the ultraviolet is ∼20%, and in the visible wavelengths it is better than 10%. The amplitude of Vesta’s rotational lightcurves is ∼10% throughout the range of wavelengths we observed, but is smaller at 950 nm (∼6%) near the 1-μm band center. Contrary to earlier reports, we found no evidence for any difference between the phasing of the ultraviolet and visible/near-infrared lightcurves with respect to sub-Earth longitude. Vesta’s average spectrum between 220 and 950 nm can well be described by measured reflectance spectra of fine particle howardite-like materials of basaltic achondrite meteorites. Combining this with the in-phase behavior of the ultraviolet, visible, and near-infrared lightcurves, and the spectral slopes with respect to the rotational phase, we conclude that there is no global ultraviolet/visible reversal on Vesta. Consequently, this implies a lack of global space weathering on Vesta, as previously inferred from visible–near-infrared data.     

The influence of volcanological and sedimentological processes on diamond grade distribution in kimberlites: examples from the EKATI Diamond Mine,NWT, Canada

The distribution of diamonds within individual kimberlite pipes is poorly documented in the public domain due to the proprietary nature of the data. The study of the diamond distribution within two pipes, Fox and Koala, from the EKATI Diamond Mine, NWT, Canada, in conjunction with detailed facies models has shown several distinct relationships of deposit type and grade distribution. In both pipes, the lithological facies represent grade units which can be distinguished from each other in terms of relative size and abundance of diamonds. A positive relationship between olivine grain size and abundance with diamond size and abundance is observed, indicating that sorting of fragmental kimberlites influences diamond distribution. Though surface geological processes do not control the diamond potential of the erupting magma, they can be responsible for concentrating diamonds into economically significant proportions. A good understanding of the eruption, transport and depositional processes responsible for the individual lithological units and the diamond distribution within them is important for successful resource estimation. This may lead to recognition of areas suitable for selective mining, making a marginal deposit economic.     

The use of reference models from a priori data to guide 2D inversion of electrical resistivity tomography data

Clay-rich till plains cover much of the UK. Such sites are attractive locations for landfills, since clay aquitards lower the risk of landfill leachate entering groundwater. However, such tills often contain sand and gravel bodies that can act as leachate flow routes. Such bodies may not be detected by conventional site investigation techniques such as drilling boreholes and trial pitting. A method of guided inversion, where a priori data are used to construct structural reference models for use in inverting electrical resistivity tomography data, was proposed as a tool to improve the detection of sand and gravel bodies within clay-rich till sequences.
Following a successful 2D guided inversion synthetic modelling study, a field study was undertaken. Wenner 2D electrical resistivity tomography lines, resistivity cone penetrometry bores and electromagnetic induction ground resistivity data were collected over a site on the East Yorkshire coast, England, where sand and gravel lenses were known to exist from cliff exposures. A number of equally valid geoelectrical models were constructed using the electromagnetic and resistivity cone data. These were used as structural reference models in the inversion of the resistivity tomography data. Blind inversion using an homogenous reference model was also carried out for comparison.
It was shown for the first time that the best solution model produced by 2D inversion of one data set with a range of structural reference models could be determined by using the l ₂ model misfit between the solution models and associated reference models (reference misfit) as a proxy for the l ₂ misfit between the solution models and the synthetic model or 'best-guess' geoelectrical model (true misfit). The 2D methodology developed here is applicable in clay-rich till plains containing sand and gravel bodies throughout the UK.     

Modeling the natural degradation of earthworks

Stratigraphic, topographic, and ground‐penetrating radar data obtained from a ca. 1800‐year‐old embankment and adjacent ditch at the Hopewell Mound Group (Chillicothe, Ohio) are used to validate the archaeological application of a simple finite‐differences diffusion model employed frequently to assess geomorphic change in natural landscapes. Although diffusion models have been used to describe the topographic degradation of landforms in a variety of geomorphic terrains, the approach has not been applied to ancient earthworks in an archaeological context. The results of this study indicate that a variety of initial earthwork forms can result in the sinusoidal profile apparent on the current landscape. Using the model results to interpret the field data, we suggest that the initial embankment form was steeper and the adjacent ditch was deeper. As a result of natural degradation processes, the earthwork widened and flattened over time. These results have broad implications for any study aimed at: (1) assessing the function of original earthwork forms, (2) determining the formation processes of complicated stratigraphies or artifact assemblages, (3) estimating the time and labor investment required for construction, or (4) identifying the socio‐political structures necessary to build earthworks. © 2005 Wiley Periodicals, Inc.     

利用大地电磁资料对玄武岩成像的可行性

The magnetotelluric （MT） method has been among the favorite supporting tools for seismic imaging of sub-salt and sub-basalt targets. In this paper we present an example from Kachchh, India （where basaltic rocks overlie Mesozoic sedimentary rocks）, and discuss the feasibility of using MT method as an exploration tool in this geological setting. Our results highlight the difference in magnetotelluric response caused by the thin intrabasalt layering. The key issue addressed in this paper is what MT can and cannot provide in such geological settings. First, we compute apparent resistivity and phase response curves using representative resistivity-depth models and borehole data from the study area. Later, we compare these results to assess the plausibility of using MT to image the sub-volcanic sediments at Kachchh. Finally, we substantiate our discussion through one-dimensional inversion of the field observed MT data from this region that exhibits poor sensitivity of MT for thin basalt layers.     

Measuring the ocean wave directional spectrum ‘First Five’ with HF radar

Ocean Dynamics - An ability to reliably measure the first five Fourier coefficients of the directional distribution of ocean wave energy is becoming an international requirement for any directional...     

The ability of detainment bunds to decrease sediments transported from pastoral catchments in surface runoff

Erosion leading to sedimentation in surface water may disrupt aquatic habitats and deliver sediment-bound nutrients that contribute to eutrophication. Land use changes causing loss of native vegetation have accelerated already naturally high erosion rates in New Zealand and increased sedimentation in streams and lakes. Sediment-bound phosphorus (P) makes up 71–79% of the 17–19 t P y⁻¹ delivered from anthropogenic sources to Lake Rotorua in New Zealand. Detainment bunds (DBs) were first implemented in the Lake Rotorua catchment in 2010 as a strategy to address P losses from pastoral agriculture. The bunds are 1.5–2 m high earthen stormwater retention structures constructed across the flow path of targeted low-order ephemeral streams with the purpose of temporarily ponding runoff on productive pastures. The current DB design protocol recommends a minimum pond volume of 120 m³ ha⁻¹ of contributing catchment with a maximum pond storage capacity of 10 000 m³. No previous study has investigated the ability of DBs to decrease annual suspended sediment (SS) loads leaving pastoral catchments. Annual SS yields delivered to two DBs with 20 ha and 55 ha catchments were 109 and 28 kg SS ha⁻¹, respectively, during this 12-month study. The DBs retained 1280 kg (59%) and 789 kg (51%) of annual SS loads delivered from the catchments as a result of the bunds' ability to impede stormflow and facilitate soil infiltration and sediment deposition. The results of this study highlight the ability of DBs to decrease SS loads transported from pastures in surface runoff, even during large storm events, and suggests DBs are able to reduce P loading in Lake Rotorua.