A Critical Comparison of Three Methods for Time-Lapse Time-Shift Calculation

Time-shift, one of the most popular time-lapse seismic attributes, has been widely used in dynamic reservoir characterization by linking it with pressure and geomechanical changes. Therefore, it is important to select appropriate calculation methods according to different time-lapse seismic data quality and time-shift magnitude. To date, there have been various published works comparing different time-shift calculation methods and discussing their advantages and disadvantages. However, most of these comparisons are based only on synthetic tests or single field applications. As the quality of time-lapse seismic data and time-shift magnitude can vary in different fields, one method may not work consistently well for each case. In this paper, a critical comparison of three different time-shift calculation techniques (Hale’s fast cross-correlation, Rickett’s non-linear inversion, and Whitcombe’s correlated leakage method) is provided. The three methods are applied to a set of synthetic data sets that are designed to account for various seismic noise and time-shift magnitudes. They are also applied to four real time-lapse seismic data sets from three North Sea fields. The calculated time-shift results are compared with the input (in synthetic tests) or the real observations from information such as seabed subsidence and compaction (in field applications). Both qualitative and quantitative comparisons are performed. At the end, each of the time-shift methods is evaluated based on different aspects, and the most appropriate method is suggested for each data scenario. All three time-shift methods are found to successfully measure time-shifts. However, Rickett’s non-linear inversion is the most outstanding method, as it gives smooth time-shifts with relatively good accuracy, and the derived time strains are more stable and interpretable.

     

Simulating the effects of spatial configurations of agricultural ditch drainage networks on surface runoff from agricultural catchments

The study of runoff is a crucial issue because it is closely related to flooding, water quality and erosion. In cultivated catchments, agricultural ditch drainage networks are known to influence runoff. As anthropogenic elements, agricultural ditch drainage networks can therefore be altered to better manage surface runoff in cultivated catchments. However, the relationship between the spatial configuration, i.e. the density and the topology, of agricultural ditch drainage networks and surface runoff in cultivated catchments is not understood. We studied this relationship by using a random network simulator that was coupled to a distributed hydrological model. The simulations explored a large variety of spatial configurations corresponding to a thousand stochastic agricultural ditch drainage networks on a 6.4 km² Mediterranean cultivated catchment. Next, several distributed hydrological functions were used to compute water flow paths and runoff for each simulation. The results showed that (i) denser networks increased the drained volume and the peak discharge and decreased hillslopes runoff, (ii) greater network density did not affect the surface runoff any further above a given network density, (iii) the correlation between network density and runoff was weaker for small subcatchments (< 2 km²) where the variability in the drained area that resulted from changes in agricultural ditch drainage networks increased the variability of runoff and (iv) the actual agricultural ditch drainage network appeared to be well optimized for managing runoff as compared with the simulated networks. Finally, our results highlighted the role of agricultural ditch drainage networks in intercepting and decreasing overland flow on hillslopes and increasing runoff in drainage networks. Copyright © 2011 John Wiley & Sons, Ltd.     

Chronology of glaciation and deglaciation during the Loch Lomond (Younger Dryas) Stade in the Scottish Highlands: implications of recalibrated 10Be exposure ages

Considerable uncertainty surrounds the timing of glacier advance and retreat during the Younger Dryas or Loch Lomond Stade (LLS) in the Scottish Highlands. Some studies favour ice advance until near the end of the stade (c. 11.7 ka), whereas others support the culmination of glacier advance in mid‐stade (c. 12.6–12.4 ka). Most published ¹⁰ Be exposure ages reported for boulders on moraines or deglacial sites post‐date the end of the LLS, and thus appear to favour the former view, but recalibration of 33 ¹⁰ Be ages using a locally derived ¹⁰ Be production rate and assuming rock surface erosion rates of zero to 1 mm ka^?1 produces exposure ages 130–980 years older than those originally reported. The recalibrated ages are filtered to exclude anomalous data, and then employed to generate aggregate probability density distributions for the timing of moraine deposition and deglaciation. The results suggest that the most probable age for the timing of the deposition of the sampled outermost moraines lies in the interval 12.4–12.1 ka or earlier. Deglacial ages obtained for sites inside Loch Lomond Stadial glacier limits imply that glaciers at some or all of the sampled sites were retreating prior to 12.1 ka. Use of aggregated data does not exclude the possibility of asynchronous glacier behaviour at different sites, but confirms that some glaciers reached their maximum limits and began to retreat several centuries before the rapid warming that terminated the LLS at 11.7–11.6 ka, consistent with the retrodictions of recent numerical modelling experiments and with geomorphological evidence for gradual oscillatory ice‐margin retreat under stadial conditions.     

Gold dispersion in a tropical rainforest weathering profile at Dondo Mobi, Gabon

At Dondo Mobi in the gold district of Eteke, south Gabon, gold mineralization was studied in a tropical forest setting. The mineralization occurs in quartz veins within amphibolites of an Archaean gneiss-amphibolite series. Gold distribution patterns were studied in the different horizons of the weathering profile and in the different grain-size fractions of the materials sampled from three pits sunk in the weathered zone: upslope (pit P1), directly overlying (pit P2) and downslope (pit P3) of the mineralization.The weathering profile consists of an upper, thick, loose sandy argillaceous horizon (H1), an irregular nodular horizon with laterite nodules or blocks (H2) and a saprolite layer (H3) up to 70 m thick. In each pit, the specific geochemical signature of the bedrock is recognized in the three horizons of the weathering profile. Some groups of elements, e.g., Cr-Ni or Ba-V-P₂O₃, characterize amphibolite or black shale, respectively. Signal attenuations in the upper horizon can be explained by a homogenization effect. That can be related to a mushroom supergene dispersion, simultaneously inducing vertical decrease and lateral enrichment in element contents. Gold also mushrooms extensively in the different horizons of the weathering profile. The dispersion pattern is somewhat anisotropic, a strong enlargement being observed in the upper (H1) and intermediate nodular (H2) horizons alike. However, the evolution of Au distribution is not the same for all grain-size fractions: (a) in the finest fraction, Au is regularly distributed in the weathering profile at the anomalous top and tends to preserve an equivalent level of concentration laterally; (b) in the coarsest fraction, the maximum Au content is found in the upper horizon just above the mineralized structure; it is rapidly decreases laterally, in the surface halo.The Au geochemical dispersion halo consists mainly of fine-grained gold developed from rather coarse-grained mineralization. Its characteristics suggest that Au dispersion was produced as a result of Au dissolution together with collapse of the weathering profile causing gold migration over short distances.The minimal volumetric reduction of saprolite in the upper horizon can be estimated to be in the range of 60%, on the basis of the content of the different grain-size fractions and the density variations in the various horizons.In exploration, geochemical techniques, even using wide sample spacing, are proving quite adequate to prospect for Au in tropical rainforest weathering environment. Indeed, significant, well-contrasted dispersion haloes facilitate detection of rather confined mineralization.     

Inversion for seismic anisotropy using genetic algorithms1

A general inversion scheme based on a genetic algorithm is developed to invert seismic observations for anisotropic parameters. The technique is applied to the inversion of shear-wave observations from two azimuthal VSP data sets from the Conoco test site in Oklahoma. Horizontal polarizations and time-delays are inverted for hexagonal and orthorhombic symmetries. The model solutions are consistent with previous studies using trial and error matching of full waveform synthetics. The shear-wave splitting observations suggest the presence of a shear-wave line singularity and are consistent with a dipping fracture system which is known to exist at the test site. Application of the inversion scheme prior to full waveform modelling demonstrates that a considerable saving in time is possible whilst retaining the same degree of accuracy.     

Prediction of tension crack location and riverbank erosion hazards along destabilized channels

Riverbank erosion, associated sedimentation and land loss hazards are a land management problem of global significance and many attempts to predict the onset of riverbank instability have been made. Recently, Osman and Thorne (1988) have presented a Culmann-type analysis of the stability of steep, cohesive riverbanks; this has the potential to be a considerable improvement over previous bank stability theories, which do not account for bank geometry changes due to toe scour and lateral erosion. However, in this paper it is shown that the existing Osman-Thorne model does not properly incorporate the influence of tension cracking on bank stability since the location of the tension crack on the floodplain is indirectly determined via calculation or arbitrary specification of the tension crack depth. Furthermore, accurate determination of tension crack location is essential to the calculation of the geometry of riverbank failure blocks and hence prediction of land loss and bank sediment yield associated with riverbank instability and channel widening. In this paper, a rational, physically based method to predict the location of tension cracks on the floodplain behind the eroding bank face is presented and tested. A case study is used to illustrate the computational procedure required to apply the model. Improved estimates of failure block geometry using the new method may potentially be applied to improve predictions of bank retreat and floodplain land loss along river channels destabilized as a result of environmental change.     

Bias in estimating fractal dimension with the rescaled-range (R/S) technique

Fractal geostatistics are being applied to subsurface geological data as a way of predicting the spatial distribution of hydrocarbon reservoir properties. The fractal dimension is the controlling parameter in stochastic methods to produce random fields of porosity and permeability. Rescaled range (R/S)analysis has become a popular way of estimating the fractal dimension, via determination of the Hurst exponent (H). A systematic investigation has been undertaken of the bias to be expected due to a range of factors commonly inherent in borehole data, particularly downhole wireline logs. The results are integrated with a review of previous work in this area. Small datasets. overlapping samples, drift and nonstationariry of means can produce a very large bias, and convergence of estimates of H around 0.85–0.90 regardless of original fractal dimension. Nonstationarity can also account for H>1, which has been reported in the literature but which is theoretically impossible for fractal time series. These results call into question the validity of fractal stochastic models built using fractal dimensions estimated with the R/Smethod.     

On the kinetics of textural equilibration in forsterite marbles

In the aureole of the Beinn an Dubhaich granite, Skye, the minimum observed median forsterite-calcite-calcite dihedral angle varies from 110° at the olivinein isograd to about 165° (the equilibrium value) at the granite-limestone contact. Laboratory experiments were performed to investigate the kinetics of this textural change. It was found that the rate of change of the forsterite-calcite-calcite dihedral angle followed approximately first-order kinetics with an activation energy of 48±4 kJ mol^-1 for fluid-present conditions, and 90 ±4 kJ mol^-1 for fluid-absent conditions. Scanning ion imaging demonstrated that, at least in the early stages of textural change, solution-reprecipitation of the calcite was the rate determining step in the fluid-present runs. The latter result and the value of the activation energy are both consistent with the activation energy found by previous authors for (albeit zeroth order) silicate-aqueous solution solution-reprecipitation reactions. The value of activation energy for the dry data does not correspond to those for either grain boundary or volume diffusion of oxygen or magnesium in forsterite. The mechanism for textural equilibration in the fluid-absent runs is uncertain. Application of the experimentally-derived rate equation to the Beinn an Dubhaich marbles gave activation energies higher than those obtained experimentally. It is concluded from consideration of grain growth effects that activation energies derived from the Beinn an Dubhaich marbles probably reflect textural equilibration under predominantly fluid-absent conditions.