Pre‐drill Groundwater Geochemistry in the Karoo Basin,South Africa

Enhanced production of unconventional hydrocarbons in the United States has driven interest in natural gas development globally, but simultaneously raised concerns regarding water quantity and quality impacts associated with hydrocarbon extraction. We conducted a pre‐development assessment of groundwater geochemistry in the critically water‐restricted Karoo Basin, South Africa. Twenty‐two springs and groundwater samples were analyzed for major dissolved ions, trace elements, water stable isotopes, strontium and boron isotopes, hydrocarbons and helium composition. The data revealed three end‐members: a deep, saline groundwater with a sodium‐chloride composition, an old, deep freshwater with a sodium‐bicarbonate‐chloride composition and a shallow, calcium‐bicarbonate freshwater. In a few cases, we identified direct mixing of the deep saline water and shallow groundwater. Stable water isotopes indicate that the shallow groundwater was controlled by evaporation in arid conditions, while the saline waters were diluted by apparently fossil meteoric water originated under wetter climatic conditions. These geochemical and isotopic data, in combination with elevated helium levels, suggest that exogenous fluids are the source of the saline groundwater and originated from remnant seawater prior to dilution by old meteoric water combined with further modification by water‐rock interactions. Samples with elevated methane concentrations (>14 ccSTP/kg) were strongly associated with the sodium‐chloride water located near dolerite intrusions, which likely provide a preferential pathway for vertical migration of deeply sourced hydrocarbon‐rich saline waters to the surface. This pre‐drill evaluation indicates that the natural migration of methane‐ and salt‐rich waters provides a source of geogenic contamination to shallow aquifers prior to shale gas development in the Karoo Basin.     

Evaluating new fault-controlled hydrothermal dolomitization models: Insights from the Cambrian Dolomite,Western Canadian Sedimentary Basin

Fault-controlled hydrothermal dolomitization in tectonically complex basins can occur at any depth and from different fluid compositions, including ‘deep-seated’, ‘crustal’ or ‘basinal’ brines. Nevertheless, many studies have failed to identify the actual source of these fluids, resulting in a gap in our knowledge on the likely source of magnesium of hydrothermal dolomitization. With development of new concepts in hydrothermal dolomitization, the study aims in particular to test the hypothesis that dolomitizing fluids were sourced from either seawater, ultramafic carbonation or a mixture between the two by utilizing the Cambrian Mount Whyte Formation as an example. Here, the large-scale dolostone bodies are fabric-destructive with a range of crystal fabrics, including euhedral replacement (RD1) and anhedral replacement (RD2). Since dolomite is cross-cut by low amplitude stylolites, dolomitization is interpreted to have occurred shortly after deposition, at a very shallow depth (<1 km). At this time, there would have been sufficient porosity in the mudstones for extensive dolomitization to occur, and the necessary high heat flows and faulting associated with Cambrian rifting to transfer hot brines into the near surface. While the δ¹⁸O_water and ⁸⁷Sr/⁸⁶Sr ratios values of RD1 are comparable with Cambrian seawater, RD2 shows higher values in both parameters. Therefore, although aspects of the fluid geochemistry are consistent with dolomitization from seawater, very high fluid temperature and salinity could be suggestive of mixing with another, hydrothermal fluid. The very hot temperature, positive Eu anomaly, enriched metal concentrations, and cogenetic relation with quartz could indicate that hot brines were at least partially sourced from ultramafic rocks, potentially as a result of interaction between the underlying Proterozoic serpentinites and CO₂-rich fluids. This study highlights that large-scale hydrothermal dolostone bodies can form at shallow burial depths via mixing during fluid pulses, providing a potential explanation for the mass balance problem often associated with their genesis.     

勘查植物地球化学在我国不同地球化学景观区的应用现状及展望

勘查植物地球化学是通过分析植物器官来获取成矿元素信息的地球化学方法。我国学者在干旱荒漠区、森林沼泽区、半干旱和湿润半湿润中低山景观区开展了系统的勘查植物地球化学试验,厘定了一批适用于勘查植物地球化学调查的植物,为该方法的推广及应用奠定了良好基础。今后应进一步加强对勘查植物地球化学异常产生的机理、勘查植物地球化学数据库及指南建立、勘查植物地球化学异常遥感应用等方面的研究,更好地利用勘查植物地球化学方法服务矿产勘查工作。     

Periglacial trimlines,former nunataks and the altitude of the last ice sheet in Wester Ross,northwest Scotland

High-level weathering limits separating ice-scoured topography from frost-weathered detritus were identified on 28 mountains in Wester Ross at altitudes of 700–960 m, and a further 22 peaks support evidence of ice scouring to summit level. Weathering limits are defined most clearly on sandstone and gneiss, which have resisted frost shattering during the Late Devensian Lateglacial, but can also be distinguished on schists and quartzite. Schmidt hammer measurements and analyses of clay mineral assemblages indicate significantly more advanced rock and soil weathering above the weathering limits. The persistence of gibbsite above weathering limits indicates that they represent the upper limit of Late Devensian glacial erosion. The regular decline of weathering-limit altitudes along former flowlines eliminates the possibility that the weathering limits represent former thermal boundaries between protective cold-based and erosive warm-based ice. The weathering limits are therefore interpreted as periglacial trimlines that define the maximum surface altitude of the last ice sheet. Calculated basal shear stresses of 50–95 kPa are consistent with this interpretation. Reconstruction of ice-sheet configuration indicates that the former ice-shed lay above 900 m along the present watershed, and that the ice surface descended northwestwards, with broad depressions along major troughs and localised domes around independent centres of ice dispersal. Extrapolation of the ice surface gradient and altitude suggests that the ice sheet did not overrun the Outer Hebrides, but was confluent with the independent Outer Hebrides ice-cap in the North Minch basin. Erratics located up to 140 m above the reconstructed ice surface are inferred to have been emplaced by a pre-Late Devensian ice sheet (or ice sheets) of unknown age. © 1997 John Wiley & Sons, Ltd.     

Deep neural network application for 4D seismic inversion to changes in pressure and saturation: Optimizing the use of synthetic training datasets

In this work, we tackle the challenge of quantitative estimation of reservoir dynamic property variations during a period of production, directly from four-dimensional seismic data in the amplitude domain. We employ a deep neural network to invert four-dimensional seismic amplitude maps to the simultaneous changes in pressure, water and gas saturations. The method is applied to a real field data case, where, as is common in such applications, the data measured at the wells are insufficient for properly training deep neural networks, thus, the network is trained on synthetic data. Training on synthetic data offers much freedom in designing a training dataset, therefore, it is important to understand the impact of the data distribution on the inversion results. To define the best way to construct a synthetic training dataset, we perform a study on four different approaches to populating the training set making remarks on data sizes, network generality and the impact of physics-based constraints. Using the results of a reservoir simulation model to populate our training datasets, we demonstrate the benefits of restricting training samples to fluid flow consistent combinations in the dynamic reservoir property domain. With this the network learns the physical correlations present in the training set, incorporating this information into the inference process, which allows it to make inferences on properties to which the seismic data are most uncertain. Additionally, we demonstrate the importance of applying regularization techniques such as adding noise to the synthetic data for training and show a possibility of estimating uncertainties in the inversion results by training multiple networks.     

Fluvial dynamics and 14C-10Be disequilibrium on the Bolivian Altiplano

Determining sediment transfer times is key to understanding source-to-sink dynamics and the transmission of environmental signals through the fluvial system. Previous work on the Bolivian Altiplano applied the in situ cosmogenic ¹⁴C-¹⁰Be-chronometer to river sands and proposed sediment storage times of ~10–20 kyr in four catchments southeast of Lake Titicaca. However, the fidelity of those results hinges upon isotopic steady-state within sediment supplied from the source area. With the aim of independently quantifying sediment storage times and testing the ¹⁴C-¹⁰Be steady-state assumption, we dated sediment storage units within one of the previously investigated catchments using radiocarbon dating, cosmogenic ¹⁰Be-²⁶Al isochron burial dating, and ¹⁰Be-²⁶Al depth-profile dating. Palaeosurfaces appear to preserve remnants of a former fluvial system, which has undergone drainage reversal, reduction in catchment area, and local isostatic uplift since ~2.8 Ma. From alluvium mantling the palaeosurfaces we gained a deposition age of ~580 ka, while lower down fluvial terraces yielded ≤34 ka, and floodplains ~3–1 ka. Owing to restricted channel connectivity with the terraces and palaeosurfaces, the main source of channel sediment is via reworking of the late Holocene floodplain. Yet modelling a set of feasible scenarios reveals that floodplain storage and burial depth are incompatible with the ¹⁴C-¹⁰Be disequilibrium measured in the channel. Instead we propose that the ¹⁴C-¹⁰Be offset results from: (i) non-uniform erosion whereby deep gullies supply hillslope-derived debris; and/or (ii) holocene landscape transience associated with climate or human impact. The reliability of the ¹⁴C-¹⁰Be chronometer vitally depends upon careful evaluation of sources of isotopic disequilibrium in a wide range of depositional and erosional landforms in the landscape. © 2018 John Wiley & Sons, Ltd.     

Intra-survey reservoir fluctuations – implications for quantitative 4D seismic analysis

During the time taken for seismic data to be acquired, reservoir pressure may fluctuate as a consequence of field production and operational procedures and fluid fronts may move significantly. These variations prevent accurate quantitative measurement of the reservoir change using 4D seismic data. Modelling studies on the Norne field simulation model using acquisition data from ocean-bottom seismometer and towed streamer systems indicate that the pre-stack intra-survey reservoir fluctuations are important and cannot be neglected. Similarly, the time-lapse seismic image in the post-stack domain does not represent a difference between two states of the reservoir at a unique base and monitor time, but is a mixed version of reality that depends on the sequence and timing of seismic shooting. The outcome is a lack of accuracy in the measurement of reservoir changes using the resulting processed and stacked 4D seismic data. Even for perfect spatial repeatability between surveys, a spatially variant noise floor is still anticipated to remain. For our particular North Sea acquisition data, we find that towed streamer data are more affected than the ocean-bottom seismometer data. We think that this may be typical for towed streamers due to their restricted aperture compared to ocean-bottom seismometer acquisitions, even for a favourable time sequence of shooting and spatial repeatability. Importantly, the pressure signals on the near and far offset stacks commonly used in quantitative 4D seismic inversion are found to be inconsistent due to the acquisition timestamp. Saturation changes at the boundaries of fluid fronts appear to show a similar inconsistency across sub-stacks. We recommend that 4D data are shot in a consistent manner to optimize aerial time coverage, and that additionally, the timestamp of the acquisition should be used to optimize pre-stack quantitative reservoir analysis.     

Late Holocene floodplain incision and alluvial fan formation in the central Grampian Highlands,Scotland: chronology,environment and implications

A section cut across an alluvial fan and the underlying floodplain terrace in the central Grampian Highlands provides an unusually complete record of late Holocene events. At ca. 2.7–2.4 cal kyr BP floodplain aggradation was replaced by net floodplain incision. Pollen evidence and charcoal counts provide no evidence for contemporaneous anthropogenic landscape change, and the timing of the transition suggests that it reflects an increase in high-magnitude erosive flood events following overall climatic deterioration. The overlying fan was deposited by torrential hyperconcentrated flows during three brief storm-generated depositional events at ca. 2.2–2.1, 1.9–1.8 and 0.9–0.7 cal kyr BP, separated and succeeded by prolonged periods of stability and peat accumulation. During these three events, a cumulative total of ca. 6750 m³ of sediment was deposited, probably in no more than a few hours over a timescale of two millennia. These findings imply that proposed links between human activity and the development of alluvial fans or debris cones require reassessment, and that different elements of the Holocene alluvial landscape have responded in different ways to the same climatic inputs. Aggregation of dating evidence relating to aggradation or incision of alluvial landforms at different scales therefore may produce misleading results. Copyright © 1999 John Wiley & Sons, Ltd.