Numerical assessment of potential impacts of hydraulically fractured Bowland Shale on overlying aquifers

Natural gas extracted from hydraulically fractured shale formations potentially has a big impact on the global energy landscape. However, there are concerns of potential environmental impacts of hydraulic fracturing of the shale formations, particularly those related to water quality. To evaluate the potential impact of hydraulically fractured shale on overlying aquifers, we conduct realizations of numerical modeling simulations to assess fluid flow and chloride transport from a synthetic Bowland Shale over a period of 11,000 years. The synthetic fractured shale was represented by a three-dimensional discrete fracture model that was developed by using the data from a Bowland Shale gas exploration in Lancashire, UK. Chloride mass exchange between fractures and the rock matrix was fully accounted for in the model. The assessment was carried out to investigate fluid and chloride mass fluxes before, during, and after hydraulic fracturing of the Bowland Shale. Impacts of the upward fracture height and aperture, as well as hydraulic conductivity of the multilayered bedrock system, are also included this assessment. This modeling revealed that the hydraulically fractured Bowland Shale is unlikely to pose a risk to its overlying groundwater quality when the induced fracture aperture is ≤200 µm. With the fracture aperture ≥1000 µm, the upward chloride flux becomes very sensitive to the upward fracture height growth and hydraulic conductivity of the multilayered bedrock system. In the extremely unlikely event of the upward fracture growth directly connecting the shale formation to the overlying Sherwood Sandstone aquifer with the fracture aperture ≥1000 µm, the upward chloride mass flux could potentially pose risks to the overlying aquifer in 100 years. The model study also revealed that the upward mass flux is significantly intercepted by the horizontal mass flux within a high permeable layer between the Bowland Shale and its overlying aquifers, reducing further upward flux toward the overlying aquifers.

Integrating geophysical data and upscaling techniques to support regional groundwater flow modelling: A practical example of the Lagan Valley, Northern Ireland

When studying heterogeneous aquifer systems, especially at regional scale, a degree of generalization is anticipated. This can be due to sparse sampling regimes, complex depositional environments or lack of accessibility to measure the subsurface. This can lead to an inaccurate conceptualization which can be detrimental when applied to groundwater flow models. It is important that numerical models are based on observed and accurate geological information and do not rely on the distribution of artificial aquifer properties. This can still be problematic as data will be modelled at a different scale to which it was collected. It is proposed here that integrating geophysics and upscaling techniques can assist in a more realistic and deterministic groundwater flow model. In this study, the sedimentary aquifer of the Lagan Valley in Northern Ireland is chosen due to intruding sub-vertical dolerite dykes. These dykes are of a lower permeability than the sandstone aquifer. The use of airborne magnetics allows the delineation of heterogeneities, confirmed by field analysis. Permeability measured at the field scale is then upscaled to different levels using a correlation with the geophysical data, creating equivalent parameters that can be directly imported into numerical groundwater flow models. These parameters include directional equivalent permeabilities and anisotropy. Several stages of upscaling are modelled in finite element. Initial modelling is providing promising results, especially at the intermediate scale, suggesting an accurate distribution of aquifer properties. This deterministic based methodology is being expanded to include stochastic methods of obtaining heterogeneity location based on airborne geophysical data. This is through the Direct Sample method of Multiple-Point Statistics (MPS). This method uses the magnetics as a training image to computationally determine a probabilistic occurrence of heterogeneity. There is also a need to apply the method to alternate geological contexts where the heterogeneity is of a higher permeability than the host rock.

Integrating Multiple Point Statistics with Aerial Geophysical Data to assist Groundwater Flow Models

The process of accounting for heterogeneity has made significant advances in statistical research, primarily in the framework of stochastic analysis and the development of multiple-point statistics (MPS). Among MPS techniques, the direct sampling (DS) method is tested to determine its ability to delineate heterogeneity from aerial magnetics data in a regional sandstone aquifer intruded by low-permeability volcanic dykes in Northern Ireland, UK. The use of two two-dimensional bivariate training images aids in creating spatial probability distributions of heterogeneities of hydrogeological interest, despite relatively ‘noisy’ magnetics data (i.e. including hydrogeologically irrelevant urban noise and regional geologic effects). These distributions are incorporated into a hierarchy system where previously published density function and upscaling methods are applied to derive regional distributions of equivalent hydraulic conductivity tensor K. Several K models, as determined by several stochastic realisations of MPS dyke locations, are computed within groundwater flow models and evaluated by comparing modelled heads with field observations. Results show a significant improvement in model calibration when compared to a simplistic homogeneous and isotropic aquifer model that does not account for the dyke occurrence evidenced by airborne magnetic data. The best model is obtained when normal and reverse polarity dykes are computed separately within MPS simulations and when a probability threshold of 0.7 is applied. The presented stochastic approach also provides improvement when compared to a previously published deterministic anisotropic model based on the unprocessed (i.e. noisy) airborne magnetics. This demonstrates the potential of coupling MPS to airborne geophysical data for regional groundwater modelling.

Diagenetic Alterations and Reservoir Quality Evolution of Lower Cretaceous Fluvial Sandstones: Nubian Formation, Sirt Basin, North-Central Libya

Lower Cretaceous meandering and braided fluvial sandstones of the Nubian Formation form some of the most important subsurface reservoir rocks in the Sirt Basin, north-central Libya. Mineralogical, petrographical and geochemical analyses of sandstone samples from well BB6-59, Sarir oilfield, indicate that the meandering fluvial sandstones are fine- to very fine-grained subarkosic arenites (av. Q91F5L4), and that braided fluvial sandstones are medium- to very coarse-grained quartz arenites (av. Q96F3L1). The reservoir qualities of these sandstones were modified during both eodiagenesis (ca. <70oC; <2 km) and mesodiagenesis (ca. >70oC; >2km). Reservoir quality evolution was controlled primarily by the dissolution and kaolinitization of feldspars, micas and mud intraclasts during eodiagenesis, and by the amount and thicknessof grain-coating clays, chemical compaction and quartz overgrowths during mesodiagenesis. However, dissolution and kaolinitization of feldspars, micas and mud intraclasts resulted in the creation of intercrystalline micro- and mouldic macro-porosity and permeability during eodiagenesis, which were more widespread in braided fluvial than in meandering fluvial sandstones. This was because of the greater depositional porosity and permeability in the braided fluvial sandstones which enhanced percolation of meteoric waters. The development of only limited quartz overgrowths in the braided fluvial sandstones, in which quartz grains are coated by thick illite layers, retained high porosity and permeability (12-23 % and 30- 600 mD). By contrast, meandering fluvial sandstones underwent porosity loss as a result of quartz overgrowth development on quartz grains which lack or have thin and incomplete grain-coating illite (2-15 % and 0-0.1mD). Further loss of porosity in the meandering fluvial sandstones occurred as a result of chemical compaction (pressuredissolution) induced by the occurrence of micas along grains contacts. Otherdiagenetic alterations, such as the growth of pyrite, siderite, dolomite/ankerite and albitization, had little impact on reservoir quality. The albitization of feldspars may have had minor positive influence on reservoir quality throughthe creation of intercrystalline micro-porosity between albite crystals.The results of this study show that diagenetic modifications of the braided and meandering fluvial sandstones in the Nubian Formation, and resulting changes in reservoir quality, are closely linked to depositional porosity and permeability. They are also linked to the thickness of grain-coating infiltrated clays, and to variations in detrital composition, particularly the amounts of mud intraclasts, feldspars and mica grains as well as climatic conditions.

Development of a biocidal treatment regime to inhibit biological growths on cultural heritage: BIODAM

Existing chemical treatments to prevent biological damage to monuments often involve considerable amounts of potentially dangerous and even poisonous biocides. The scientific approach described in this paper aims at a drastic reduction in the concentration of biocide applications by a polyphasic approach of biocides combined with cell permeabilisers, polysaccharide and pigment inhibitors and a photodynamic treatment. A variety of potential agents were screened to determine the most effective combination. Promising compounds were tested under laboratory conditions with cultures of rock deteriorating bacteria, algae, cyanobacteria and fungi. A subsequent field trial involved two sandstone types with natural biofilms. These were treated with multiple combinations of chemicals and exposed to three different climatic conditions. Although treatments proved successful in the laboratory, field trials were inconclusive and further testing will be required to determine the most effective treatment regime. While the most effective combination of chemicals and their application methodology is still being optimised, results to date indicate that this is a promising and effective treatment for the control of a wide variety of potentially damaging organisms colonising stone substrates

The influence of substrate texture on early biological colonization

Biological colonization of stone is a major concern in the preservation and presentation of cultural heritage. Colonization is typically associated with unpleasant soiling, and varying degrees of biodeterioration. A better understanding of why organisms grow where they do, will aid in
developing preventative, and treatment methods for biosoiling of cultural heritage. Sandstone exposure trials were set up at nine different locations across Northern Ireland to investigate the influences of local climate, local environmental,and micro-climatic factors on the early stages (up to 21 months) of biological colonization.
Results showed that, green and yellow soiling occurred on tooled stone surfaces, whereas darkening occurred preferentially on smooth surfaces. It is likely that different populations of organisms occur on these surfaces with green algae occurring on tooled surfaces due to slower drying rates (i.e. prolonged moisture retention), and cyanobacteria and fungi thriving on smooth surfaces due to their ability to withstand moisture fluctuation.

The inventorying of building stone performance-in-use to assess susceptibility to decay

The Natural Stone Database for Northern Ireland was constructed to address the paucity of information available to stone conservation practitioners within the region. Almost 2000 listed buildings and monuments were surveyed over three years to produce an interactive GIS database. This contained information on stone sources, together with details of stone condition and decay processes. This paper uses elements of this GIS to investigate stone decay patterns across Northern Ireland. The results demonstrate that as the level of stone decay increases, so does the proportion of buildings with sandstone as the primary stone type. It appears that a
higher open porosity level combined with Northern Ireland’s wetter climate and maritime location leads to rapid wetting and drying cycles within sandstones. This is coupled with the ingress and crystallisation of marine and other salts within stone pores leading to considerably
higher rates of decay than for any other stone type.

Geology of the Gander Group, Gander River Ultramafic Belt and Davidsville Group in the Jonathan's Pond - Weir's Pond Area, northeast Newfoundland /

The study area is situated in NE Newfoundland between Gander Lake and the north coast and on the boundary between the Gander and Botwood tectonostratigraphic zones (Williams et al., 1974). The area is underlain by three NE trending units; the Gander Group, the Gander River Ultramafic Belt (the GRUB) and the Davidsville Group. The easternmost Gander Group consists of a thick, psammitic unit composed predominantly of psammitic schist and a thinner, mixed unit of semipelitic and pelitic schist with minor psammite. The mixed unit may stratigraphically overlie the psammitic unit or be a lateral facies equivalent of the latter. No fossils have been recovered from the Gander Group. The GRUB is a terrain of mafic and ultramafic plutonic rocks with minor pillow lava and plagiogranite. It is interpreted to be a dismembered ophiolite in thrust contact with the Gander Group. The westernmost Davidsville Group consists of a basal conglomerate, believed deposited unconformably upon the GRUB from which it was derived, and an upper unit of greywacke and slate, mostly of turbidite origin, with minor limestone and calcareous sandstone. The limestone, which lies near the base of the unit, contains Upper Llanvirn to Lower Llandeilo fossils. The Gander and Davidsville Groups display distinctly different sedimentological , structural and metamorphic histories. The Gander Group consists of quartz-rich, relatively mature sediment. It has suffered three pre-Llanvirn deformations, of which the main deformation, Dp produced a major, NE-N-facing recumbent anticline in the southern part of the study area. Middle greenschist conditions existed from D^ to D- with growth of metamorphic minerals during each dynamic and static phase. In contrast, the mineralogically immature Davidsville Group sediment contains abundant mafic and ultramafic detritus which is absent from the Gander Group. The Davidsville Group displays the effects of a single penetrative deformation with localized D_ and D_ features, all of which can be shown to postdate D_ in the Gander Group. Rotation of the flat Gander S- into a subvertical orientation near the contact with the GRUB and the Davidsville Group is believed to be a Davidsville D^ feature. Regional metamorphism in the Davidsville Group is lower greenschist with a single growth phase, MS . These sedimentological, structural and metamorphic differences between the Gander and Davidsville Groups persist even where the GRUB is absent and the two units are in contact, indicating that the tectonic histories of the Gander and Davidsville Groups are distinctly different. Structural features in the GRUB, locally the result of multiple deformations, may be the result of Gander and/or Davidsville deformations. Metamorphism is in the greenschist facies. Geochemical analyses of the pillow lava suggest that these rocks were formed in a back-arc basin. Mafic intrusives in the Gander Group appear to be the result of magraatism separate from that producing the pillow lava. The Gander Group is interpreted to be a continental rise prism deposited on the eastern margin of the Late Precambrian-Lower Paleozoic lapetus Ocean. The GRUB, oceanic crust possibly formed in a marginal basin to the west, is believed to have been thrust eastward over the Gander Group, deforming the latter, during the pre-Llanvirnian, possibly Precambrian, Ganderian Orogeny. The Middle Ordovician and younger Davidsville Group was derived from, and deposited unconformably on, this deformed terrain. Deformation of the Davidsville Group occurred during the Middle Devonian Acadian Orogeny.

Regional variations of deformation mechanisms in the Lorrain Quartzite near Whitefish Falls, Ontario, Canada /

The quartzite microfabric found in the Lorrain Formation was studied across the La Cloche syncline, along a regional north-south transect along highway 6, near Whitefish Falls, Ontario. The complete stratigraphic sequence across the syncline is preserved, and is present on each fold limb. The lithostratigraphic units with the smallest grains size and lowest mica content are located close to the core of the fold, while coarser grained mica and feldspar rich units are situated at the northern and southern most extent of the transect. Deformation mechanisms vary with lithology and with position across the fold. Pressure solution appears to be the dominant deformation mechanism in the feldspathic, micaceous and ferruginous units. In the finer grained, mica poor white medium grained and cherty sandstone units, grain boundary migration (GBM) characteristics show dominance over those of pressure solution and show high amounts of fracturing which cut migrated boundaries and therefore post date GBM. All samples across the fold display a preferred orientation of quartz c-axes. The senses of asymmetry of fabrics are found to be similar across the syncline, with the exception of the ferruginous sandstone unit. Formation of these similar fabrics synmietries can not be the result of strain related to first order folding. The mica content appears to be related to the percentage of quartz lost due to pressure solution as a result of strain; the more mica present, the less quartz was lost. Calculations based on the shape of initial grains suggest that conservatively 30% of the quartz volume has been dissolved out of the Lorrain quartzite, and potentially migrated hundreds of meters to other members of the Huronian Supergroup as there was no meso or macroscopic evidence observed in outcrop.

Textual characteristics of coarse sediments in selected streams of the Niagara Peninsula, Ontario

The streams flowing through the Niagara Escarpment are paved by coarse carbonate and sandstone sediments which have originated from the escarpment units and can be traced downstream from their source. Fifty-nine sediment samples were taken from five streams, over distances of 3,000 to 10,000 feet (915 to 3050 m), to determine downstream changes in sediment composition, textural characteristics and sorting. In addition, fluorometric velocity measurements were used in conjunction with measured -discharge and flow records to estimate the frequency of sediment movement. The frequency of sediments of a given lithology changes downstream in direct response to the outcrop position of the formations in the channels. Clasts derived from a single stratigraphic unit usually reach a maximum frequency within the first 1,000 feet (305 m) of transport. Sediments derived from formations at the top of waterfalls reach a modal frequency farther downstream than material originating at the base of waterfalls. Downstream variations in sediment size over the lengths of the study reaches reflect the changes in channel morphology and lithologic composition of the sediment samples. Linear regression analyses indicate that there is a decrease in the axial lengths between the intial and final samples and that the long axis decreases in length more rapidly than the intermediate, while the short axis remains almost constant. Carbonate sediments from coarse-grained, fossiliferous units - iii - are more variable in size than fine-grained dolostones and sandstones. The average sphericity for carbonates and sandstones increases from 0.65 to 0.67, while maximum projection sphericity remains nearly constant with an average value of 0.52. Pebble roundness increases more rapidly than either of the sphericity parameters and the sediments change from subrounded to rounded. The Hjulstrom diagram indicates that the velocities required to initiate transport of sediments with an average intermediate diameter of 10 cm range from 200 cm/s to 300 cm/s (6.6 ft./sec. to 9.8 ft./sec.). From the modal velocitydischarge relations, the flows corresponding to these velocities are greater than 3,500 cfs (99 m3s). These discharges occur less than 0.01 p~r cent (0.4 days) of the time and correspond to a discharge occurring during the spring flood.

Interlake Tissue Mills Fire Slides

The Interlake Tissue Mills storage location was at the intersection of Glendale and Merritt Street in Merritton. In the 1850s, the Beaver Cotton Mills (later Merritton Cotton Mills) was built. It was a frame building that burned down in 1881. A sandstone structure replaced the previous building. In 1906, the cotton mill closed and in 1912, the Independent Rubber Company took over the building. The mill was doubled in size. The rubber company remained there for ten years and then the building remained vacant until 1935. Interlake Tissue Mills acquired the building for storage purposes. A fire razed the larger of the plant’s buildings on Memorial Day weekend in 1961. After the fire, Domtar used the smaller building for storage for a short time. This is currently the site of the Keg Restaurant.

Interpretación sedimentológica y paleogeográfica del nivel de areniscas eocénicas de Sarrià de Ter (Prov. de Gerona)

There is in this work a description of sedimentological characters of the medium Eocene sandstone level from Sarriá de Ter (Gerona province Spain) and an interpretation of their genesis and situation in the regional paleogeography. The conclusion is that these sandstones form a deltaic channel(delta front). To wards the E and NE the become fluvial deposits, to wards the W, Plana de Vic, they become clay and silstones considered as a prodelta

Spatangoid-produced ichnofabrics (Bateig Limestone, Miocene, Spain) and the preservation of spatangoid trace fossils

A spatangoid-produced ichnofabric is described from the Miocene Bateig Limestone, SE Spain. This ichnofabric is characterized by the dominant presence of large meniscate burrows (Bichordites) produced by irregular echinoids. This constitutes an unusual mode of occurrence for spatangoid bioturbation, as their traces are most typically preserved in bases and tops of sandstone event beds. In fact, despite their important role as burrowers in modern settings (that can be extended back to the Early Cretaceous based on their body fossil record), spatangoid trace fossils (Scolicia and Bichordites) are comparatively rare. Several factors play an important role in their preservation: mechanism of burrowing, sediment characteristics, early diagenesis and presence/absence of deep-tier burrowers. Spatangoid-produced ichnofabrics, such as those from the Bateig Limestone, characterize depositional settings with intermittent deposition of event beds where there is an absence of deeper-tier bioturbation. (C) 2008 Elsevier B.V. All rights reserved

Seismic activity triggered by water wells in the Parana Basin, Brazil

Triggered seismicity is commonly associated with deep water reservoirs or injection wells where water is injected at high pressure into the reservoir rock. However, earth tremors related solely to the opening of groundwater wells are extremely rare. Here we present a clear case of seismicity induced by pore-pressure changes following the drilling of water wells that exploit a confined aquifer in the intracratonic Parana Basin of southeastern Brazil. Since 2004, shallow seismic activity, with magnitudes up to 2.9 and intensities V MM, has been observed near deep wells (120-200 m) that were drilled in early 2003 near the town of Bebedouro. The wells were drilled for irrigation purposes, cross a sandstone layer about 60-80 m thick and extract water from a confined aquifer in fractured zones between basalt flow layers. Seismic activity, mainly event swarms, has occurred yearly since 2004, mostly during the rainy season when the wells are not pumped. During the dry season when the wells are pumped almost continuously, the activity is very low. A seismographic network, installed in March 2005, has located more than 2000 microearthquakes. The events are less than 1 km deep (mostly within the 0.5 km thick basalt layer) and cover an area roughly 1.5 km x 5 km across. The seismicity generally starts in a small area and expands to larger distances with an equivalent hydraulic diffusivity ranging from 0.06 to 0.6 m(2)/s. Geophysical and geothermal logging of several wells in the area showed that water from the shallow sandstone aquifer enters the well at the top and usually forms waterfalls. The waterfalls flow down the sides of the wells and feed the confined, fractured aquifer in the basalt layer at the bottom. Two seismic areas are observed: the main area surrounds several wells that are pumped continuously during the dry season, and a second area near another well (about 10 km from the first area) that is not used for irrigation and not pumped regularly. The main area displays cyclic annual activity, but the second area does not. We explain the earthquake swarms as being triggered by pore pressure diffusion in the fractured basalt layer due to additional pressure from the newly connected surface aquifer. This reaches critically prestressed areas up to a few kilometers away from the wells. During periods of continuous pumping, the reduction of pore pressure in the confined aquifer stops the seismic activity. Our study suggests that this kind of activity may be more common than previously thought and implies that many other cases of small tremors associated with the drilling of water wells may have gone unnoticed.

Spectral analysis and modeling of microcyclostratigraphy in late Paleozoic glaciogenic rhythmites, Parana Basin, Brazil

We investigate the depositional time scale of lithological couplets (fine sandstone/siltstone-siltstone/mudstone) from two distinctive outcrops of Permo-Carboniferous glacial rhythmites in the Itarare Group (Parana Basin, Brazil). Resolving the fundamental issue of time scale for these rhythmites is important in light of recent evidence for paleosecular variation measured in these sequences. Spectral analysis and tuning of high-resolution gray scale scans of sediment core microstratigraphy, which comprises pervasive laminations, reveal a comparable spectral content at both localities, with a frequency suite interpreted as that of short-term climate variability of Recent and modern times. This evidence for decadal- to centennial-scale deposition of these lithological couplets is discussed in light of the `varvic` character, i.e., annual time scale that was previously assumed for the rhythmites.