Estimating Travel Time in Bank Filtration Systems from a Numerical Model Based on DTS Measurements

An approach is presented to determine the seasonal variations in travel time in a bank filtration system using a passive heat tracer test. The temperature in the aquifer varies seasonally because of temperature variations of the infiltrating surface water and at the soil surface. Temperature was measured with distributed temperature sensing along fiber optic cables that were inserted vertically into the aquifer with direct push equipment. The approach was applied to a bank filtration system consisting of a sequence of alternating, elongated recharge basins and rows of recovery wells. A SEAWAT model was developed to simulate coupled flow and heat transport. The model of a two‐dimensional vertical cross section is able to simulate the temperature of the water at the well and the measured vertical temperature profiles reasonably well. MODPATH was used to compute flowpaths and the travel time distribution. At the study site, temporal variation of the pumping discharge was the dominant factor influencing the travel time distribution. For an equivalent system with a constant pumping rate, variations in the travel time distribution are caused by variations in the temperature‐dependent viscosity. As a result, travel times increase in the winter, when a larger fraction of the water travels through the warmer, lower part of the aquifer, and decrease in the summer, when the upper part of the aquifer is warmer.     

Effects of discharge,wind, and tide on sedimentation in a recently restored tidal freshwater wetland

Sediment deposition is one of the key mechanisms to counteract the impact of sea level rise in tidal freshwater wetlands (TFWs). However, information about sediment deposition rates in TFWs is limited, especially for those located in the transition zone between the fluvially dominated and tidally dominated sections of a river delta where sedimentation rates are affected by the combined impact of river discharge, wind, and tides. Using a combined hydrodynamic–morphological model, we examined how hydrometeorological boundary conditions control sedimentation rates and patterns in a TFW located in the Rhine–Meuse estuary in the Netherlands. The modelling results show that net sedimentation rate increases with the magnitude of the river discharge, whereas stronger wind increasingly prevents sedimentation. Sediment trapping efficiency decreases for both increasing river discharge and wind magnitude. The impact of wind storms on the trapping efficiency becomes smaller for higher water discharge. The spatial sedimentation patterns are affected by all controls. Our study illustrates the importance of evaluating both the separate and the joint impact of discharge, wind, and tides when estimating sedimentation rates in a TFW affected by these controls. Such insights are relevant to design measures to reactivate the sedimentation process in these areas.     

Three-dimensional Marchenko internal multiple attenuation on narrow azimuth streamer data of the Santos Basin,Brazil

In recent years, a variety of Marchenko methods for the attenuation of internal multiples has been developed. These methods have been extensively tested on two-dimensional synthetic data and applied to two-dimensional field data, but only little is known about their behaviour on three-dimensional synthetic data and three-dimensional field data. Particularly, it is not known whether Marchenko methods are sufficiently robust for sparse acquisition geometries that are found in practice. Therefore, we start by performing a series of synthetic tests to identify the key acquisition parameters and limitations that affect the result of three-dimensional Marchenko internal multiple prediction and subtraction using an adaptive double-focusing method. Based on these tests, we define an interpolation strategy and use it for the field data application. Starting from a wide azimuth dense grid of sources and receivers, a series of decimation tests are performed until a narrow azimuth streamer geometry remains. We evaluate the effect of the removal of sail lines, near offsets, far offsets and outer cables on the result of the adaptive double-focusing method. These tests show that our method is most sensitive to the limited aperture in the crossline direction and the sail line spacing when applying it to synthetic narrow azimuth streamer data. The sail line spacing can be interpolated, but the aperture in the crossline direction is a limitation of the acquisition. Next, we apply the adaptive Marchenko double-focusing method to the narrow azimuth streamer field data from the Santos Basin, Brazil. Internal multiples are predicted and adaptively subtracted, thereby improving the geological interpretation of the target area. These results imply that our adaptive double-focusing method is sufficiently robust for the application to three-dimensional field data, although the key acquisition parameters and limitations will naturally differ in other geological settings and for other types of acquisition.     

Declines in oil-rates of stranded birds in the North Sea highlight spatial patterns in reductions of chronic oil pollution

Strandings of oiled seabirds are used to signal the problem of chronic oil pollution. Species-specific oil rates reflect the risk for marine birds to become oiled at sea. High oil rates were characteristic for seabirds common in areas with frequent oil spills; low oil rates for birds wintering away from the busiest shipping lanes. Declining trends in oil-rates were found, reflecting a reduction in the amount of oil intentionally discharged over the past 50 years. Spatial patterns in the risk to become oiled could be identified, when the winter distribution patterns of the affected birds were incorporated in the analysis. Declines in oil rates were most pronounced in coastal birds. These trends were consistent with tendencies to police nearshore waters more effectively than offshore waters. While levels of chronic oil pollution are much reduced, future emphasis should be to reduce chronic oiling more effectively in offshore waters.     

Virtual signals combination – phase analysis for 2D and 3D data representation (acoustic medium)

We formulate the Kirchhoff‐Helmholtz representation theory for the combination of seismic interferometry signals synthesized by cross‐correlation and by cross‐convolution in acoustic media. The approach estimates the phase of the virtual reflections from the boundary encompassing a volume of interest and subtracts these virtual reflections from the total seismic‐interferometry wavefield. The reliability of the combination result, relevant for seismic exploration, depends on the stationary‐phase and local completeness in partial coverage regions. The analysis shows the differences in the phase of the corresponding seismic interferometry (by cross‐correlation) and virtual reflector (by cross‐convolution) signals obtained by 2D and 3D formulations, with synthetic examples performed to remove water layer multiples in ocean bottom seismic (OBS) acoustic data.     

Retrieving virtual reflection responses at drill‐bit positions using seismic interferometry with drill‐bit noise

In the field of seismic interferometry, researchers have retrieved surface waves and body waves by cross‐correlating recordings of uncorrelated noise sources to extract useful subsurface information. The retrieved wavefields in most applications are between receivers. When the positions of the noise sources are known, inter‐source interferometry can be applied to retrieve the wavefields between sources, thus turning sources into virtual receivers. Previous applications of this form of interferometry assume impulsive point sources or transient sources with similar signatures. We investigate the requirements of applying inter‐source seismic interferometry using non‐transient noise sources with known positions to retrieve reflection responses at those positions and show the results using synthetic drilling noise as source. We show that, if pilot signals (estimates of the drill‐bit signals) are not available, it is required that the drill‐bit signals are the same and that the phases of the virtual reflections at drill‐bit positions can be retrieved by deconvolution interferometry or by cross‐coherence interferometry. Further, for this case, classic interferometry by cross‐correlation can be used if the source power spectrum can be estimated. If pilot signals are available, virtual reflection responses can be obtained by first using standard seismic‐while‐drilling processing techniques such as pilot cross‐correlation and pilot deconvolution to remove the drill‐bit signatures in the data and then applying cross‐correlation interferometry. Therefore, provided that pilot signals are reliable, drill‐bit data can be redatumed from surface to borehole depths using this inter‐source interferometry approach without any velocity information of the medium, and we show that a well‐positioned image below the borehole can be obtained using interferometrically redatumed reflection responses with just a simple velocity model. We discuss some of the practical hurdles that restrict the application of the proposed method offshore.     

Annama H chondrite—Mineralogy,physical properties,cosmic ray exposure,and parent body history

The fall of the Annama meteorite occurred early morning (local time) on April 19, 2014 on the Kola Peninsula (Russia). Based on mineralogy and physical properties, Annama is a typical H chondrite. It has a high Ar‐Ar age of 4.4 Ga. Its cosmic ray exposure history is atypical as it is not part of the large group of H chondrites with a prominent 7–8 Ma peak in the exposure age histograms. Instead, its exposure age is within uncertainty of a smaller peak at 30 ± 4 Ma. The results from short‐lived radionuclides are compatible with an atmospheric pre‐entry radius of 30–40 cm. However, based on noble gas and cosmogenic radionuclide data, Annama must have been part of a larger body (radius >65 cm) for a large part of its cosmic ray exposure history. The ¹⁰Be concentration indicates a recent (3–5 Ma) breakup which may be responsible for the Annama parent body size reduction to 30–35 cm pre‐entry radius.     

Tutorial: unified 1D inversion of the acoustic reflection response

Acoustic inversion in one-dimension gives impedance as a function of travel time. Inverting the reflection response is a linear problem. Recursive methods, from top to bottom or vice versa, are known and use a fundamental wave field that is computed from the reflection response. An integral over the solution to the Marchenko equation, on the other hand, retrieves the impedance at any vertical travel time instant. It is a non-recursive method, but requires the zero-frequency value of the reflection response. These methods use the same fundamental wave field in different ways. Combining the two methods leads to a non-recursive scheme that works with finite-frequency bandwidth. This can be used for target-oriented inversion. When a reflection response is available along a line over a horizontally layered medium, the thickness and wave velocity of any layer can be obtained together with the velocity of an adjacent layer and the density ratio of the two layers. Statistical analysis over 1000 noise realizations shows that the forward recursive method and the Marchenko-type method perform well on computed noisy data.