ESEM study of authigenic chlorite acid sensitivity in sandstone reservoirs

The effect of HCl on authigenic chlorite in three different sandstones has been examined uisng an Environmental Scanning Electron Microscope (ESEM), together with conventional analytical techniques. The ESEM enabled chlorites to be directly observed in situ at high magnifications during HCl treatment, and was particularly effective in allowing the same chlorite areas to be closely compared before and after acid treatment. Chlorites were reacted with 1M to 10M HCl at temperatures up to 80°C and for periods up to five months. After all treatments, chlorites show extensive leaching of iron, magnesium and aluminum, and their crystalline structure is destroyed. However, despite these major compositional and structural changes, chlorites show little or no visible evidence of acid attack, with precise morphological detail of individual plates preserved in all samples following acid treatments. Chlorite dissolution, sensu stricto, did not occur as a result of acidization of the host sandstones. Acid-treated chlorides are likely to exits in a structurally weakened state that may make them susceptible to physical disintegration during fluid flow. Accordingly, fines migration may be a significant engineering problem associated with the acidization of chlorite-bearing sandstones. © 1993.

ESEM study of illite/smectite freshwater sensitivity in sandstone reservoirs

The water sensitivity of authigenic smectite- and illite-rich illite/smectites in sandstone reservoirs has been investigated using an Environmental Scanning Electron Microscope (ESEM). The ESEM enabled the illite/smectites to be directly observed in situ at high magnification during freshwater immersion, and was also particularly effective in allowing the same selected illite/smectite areas to be closely compared before and after freshwater treatments. The tendency of authigenic smectite-rich illite/smectite to swell on contact with fresh water varies greatly. Smectite-rich illite/smectite may osmotically swell to many times its original volume to form a gel which greatly reduces porosity and permeability, or may undergo only a subtle morphological change which has little or no adverse effect on reservoir quality. Authigenic illite-rich illite/smectite in sandstones does not swell when immersed in fresh water. Even after prolonged soaking in fresh water, illite-rich illite/smectite particles retain their original morphology. Accordingly, illite-rich illite/smectite in sandstones is unlikely to cause formation damage if exposed to freshwater-based fluids. © 1993.

Pore structure characterization of North American shale gas reservoirs using USANS/SANS, gas adsorption, and mercury intrusion

Small-angle and ultra-small-angle neutron scattering (SANS and USANS), low-pressure adsorption (N2 and CO2), and high-pressure mercury intrusion measurements were performed on a suite of North American shale reservoir samples providing the first ever comparison of all these techniques for characterizing the complex pore structure of shales. The techniques were used to gain insight into the nature of the pore structure including pore geometry, pore size distribution and accessible versus inaccessible porosity. Reservoir samples for analysis were taken from currently-active shale gas plays including the Barnett, Marcellus, Haynesville, Eagle Ford, Woodford, Muskwa, and Duvernay shales. Low-pressure adsorption revealed strong differences in BET surface area and pore volumes for the sample suite, consistent with variability in composition of the samples. The combination of CO2 and N2 adsorption data allowed pore size distributions to be created for micro–meso–macroporosity up to a limit of �1000 Å. Pore size distributions are either uni- or multi-modal. The adsorption-derived pore size distributions for some samples are inconsistent with mercury intrusion data, likely owing to a combination of grain compression during high-pressure intrusion, and the fact that mercury intrusion yields information about pore throat rather than pore body distributions. SANS/USANS scattering data indicate a fractal geometry (power-law scattering) for a wide range of pore sizes and provide evidence that nanometer-scale spatial ordering occurs in lower mesopore–micropore range for some samples, which may be associated with inter-layer spacing in clay minerals. SANS/USANS pore radius distributions were converted to pore volume distributions for direct comparison with adsorption data. For the overlap region between the two methods, the agreement is quite good. Accessible porosity in the pore size (radius) range 5 nm–10 lm was determined for a Barnett shale sample using the contrast matching method with pressurized deuterated methane fluid. The results demonstrate that accessible porosity is pore-size dependent.

Modelling surface and shallow groundwater interactions in an ungauged subtropical coastal catchment using the SWAT model, Elimbah Creek, southeast Queensland, Australia

The study presented here applies the highly parameterised semi-distributed U.S. Department of Agriculture Soil and Water Assessment Tool (SWAT) to an Australian subtropical catchment. SWAT has been applied to numerous catchments worldwide and is considered to be a useful tool that is under ongoing development with contributions coming from different research groups in different parts of the world. In a preliminary run the SWAT model application for the Elimbah Creek catchment has estimated water yield for the catchment and has quantified the different sources. For the modelling period of April 1999 to September 2009 the results show that the main sources of water in Elimbah Creek are total surface runoff and lateral flow (65%). Base-flow contributes 36% to the total runoff. On a seasonal basis modelling results show a shift in the source of water contributing to Elimbah Creek from surface runoff and lateral flow during intense summer storms to base-flow conditions during dry months. Further calibration and validation of these results will confirm that SWAT provides an alternative to Australian water balance models.

Image-based flow visualisation (IBFV) to enhance interpretation of complex flow patterns within a shallow tidal barrier estuary

We applied a texture-based flow visualisation technique to a numerical hydrodynamic model of the Pumicestone Passage in southeast Queensland, Australia. The quality of the visualisations using our flow visualisation tool, are compared with animations generated using more traditional drogue release plot and velocity contour and vector techniques. The texture-based method is found to be far more effective in visualising advective flow within the model domain. In some instances, it also makes it easier for the researcher to identify specific hydrodynamic features within the complex flow regimes of this shallow tidal barrier estuary as compared with the direct and geometric based methods.

Carbonate cementation in sandstones of the Plover Formation, North West Shelf, W.A

Ancient sandstones include important reservoirs for hydrocarbons (oil and gas), but, in many cases, their ability to serve as reservoirs is heavily constrained by the effects of carbonate cements on porosity and permeability. This study investigated the controls on distribution and abundance of carbonate cements within the Jurassic Plover Formation, Browse Basin, North West Shelf, Australia. Samples were analysed petrographically with point counting of 59 thin sections and mineralogically with x-ray diffraction from two wells within the Torosa Gas Field. Selected samples were also analysed for stable isotopes of O and C. Sandstones are classified into eleven groups. Most abundant are quartzarenites and then calcareous quartzarenites. Lithology ranged between sandstones consisting of mostly quartz with scant or no carbonate in the form of cement or allochems, to sandstones with as much as 40% carbonate. The major sources of carbonate cement in Torosa 1 and Torosa 4 sandstones were found to be early, shallow marine diagenetic processes (including cementation), followed by calcite cementation and recrystallisation of cements and allochems during redistribution by meteoric waters. Blocky and sparry calcite cements, indicative of meteoric environments on the basis of stable isotope values and palaeotemperature assessment, overprinted the initial shallow marine cement phase in all cases and meteoric cements are dominant. Torosa 4 was influenced more by marine settings than Torosa 1, and thus has the greater potential for calcite cement. The relatively low compaction of calcite-cemented sandstones and the stable isotope data suggest deep burial cementation was not a major factor. Insufficient volcanic rock fragments or authigenic clay content infers alteration of feldspars was not a major source of calcite. Very little feldspar is present, altered or otherwise. Hence, increased alkalinity from feldspar dissolution is not a contributing factor in cement formation. Increased alkalinity from bacterial sulphate reduction in organic–rich fine sediments may have driven limited cementation in some samples. The main definable and significant source of diagenetic marine calcite cement originated from original marine cements and the nearby dissolution of biogenic sources (allochems) at relatively shallow depths. Later diagenetic fluids emplaced minor dolomite, but this cement did not greatly affect the reservoir quality in the samples studied.

Electrical imaging and fluid modeling of convective fingering in a shallow water-table aquifer

Unstable density-driven flow can lead to enhanced solute transport in groundwater. Only recently has the complex fingering pattern associated with free convection been documented in field settings. Electrical resistivity (ER) tomography has been used to capture a snapshot of convective instabilities at a single point in time, but a thorough transient analysis is still lacking in the literature. We present the results of a 2 year experimental study at a shallow aquifer in the United Arab Emirates that was designed to specifically explore the transient nature of free convection. ER tomography data documented the presence of convective fingers following a significant rainfall event. We demonstrate that the complex fingering pattern had completely disappeared a year after the rainfall event. The observation is supported by an analysis of the aquifer halite budget and hydrodynamic modeling of the transient character of the fingering instabilities. Modeling results show that the transient dynamics of the gravitational instabilities (their initial development, infiltration into the underlying lower-density groundwater, and subsequent decay) are in agreement with the timing observed in the time-lapse ER measurements. All experimental observations and modeling results are consistent with the hypothesis that a dense brine that infiltrated into the aquifer from a surficial source was the cause of free convection at this site, and that the finite nature of the dense brine source and dispersive mixing led to the decay of instabilities with time. This study highlights the importance of the transience of free convection phenomena and suggests that these processes are more rapid than was previously understood.

Turbulence measurements in a shallow tidal estuary : analysis based on triple decomposition

In an estuary, mixing and dispersion are the result of the combination of large scale advection and small scale turbulence which are both complex to estimate. A field study was conducted in a small sub-tropical estuary in which high frequency (50 Hz) turbulent data were recorded continuously for about 48 hours. A triple decomposition technique was introduced to isolate the contributions of tides, resonance and turbulence in the flow field. A striking feature of the data set was the slow fluctuations which exhibited large amplitudes up to 50% the tidal amplitude under neap tide conditions. The triple decomposition technique allowed a characterisation of broader temporal scales of high frequency fluctuation data sampled during a number of full tidal cycles.

High-resolution GNSS-tracked drifter for studying surface dispersion in shallow water

The use of GNSS tracked Lagrangian drifters allows more realistic quantification of fluid motion and dispersion coefficients than Eulerian techniques because such drifters are analogues of particles that are relevant to flow field characterisation and pollutant dispersion. Using the fast growing Real Time Kinematic (RTK) positioning technique derived from Global Satellite Navigation Systems (GNSS), drifters are developed for high frequency (10 Hz) sampling with position estimates to centimetre accuracy. The drifters are designed with small size and less direct wind drag to follow the sub-surface flow which characterizes dispersion in shallow waters. An analysis of position error from stationary observation indicates that the drifter can efficiently resolve motion up to 1 Hz. The result of the field deployments of the drifter in conjunction with acoustic Eulerian devices shows higher estimate of the drifter streamwise velocities. Single particle statistical analysis of field deployments in a shallow estuarine zone yielded dispersion coefficients estimate comparable to those of dye tracer studies. The drifters capture the tidal elevation during field studies in a tidal estuary.

Integrating field survey data with satellite image data to improve shallow water seagrass maps : the role of AUV and snorkeller surveys

Repeatable and accurate seagrass mapping is required for understanding seagrass ecology and supporting management decisions. For shallow (< 5 m) seagrass habitats, these maps can be created by integrating high spatial resolution imagery with field survey data. Field survey data for seagrass is often collected via snorkelling or diving. However, these methods are limited by environmental and safety considerations. Autonomous Underwater Vehicles (AUVs) are used increasingly to collect field data for habitat mapping, albeit mostly in deeper waters (>20 m). Here we demonstrate and evaluate the use and potential advantages of AUV field data collection for calibration and validation of seagrass habitat mapping of shallow waters (< 5 m), from multispectral satellite imagery. The study was conducted in the seagrass habitats of the Eastern Banks (142 km2), Moreton Bay, Australia. In the field, georeferenced photos of the seagrass were collected along transects via snorkelling or an AUV. Photos from both collection methods were analysed manually for seagrass species composition and then used as calibration and validation data to map seagrass using an established semi-automated object based mapping routine. A comparison of the relative advantages and disadvantages of AUV and snorkeller collected field data sets and their influence on the mapping routine was conducted. AUV data collection was more consistent, repeatable and safer in comparison to snorkeller transects. Inclusion of deeper water AUV data resulted in mapping of a larger extent of seagrass (~7 km2, 5 % of study area) in the deeper waters of the site. Although overall map accuracies did not differ considerably, inclusion of the AUV data from deeper water transects corrected errors in seagrass mapped at depths to 5 m, but where the bottom is visible on satellite imagery. Our results demonstrate that further development of AUV technology is justified for the monitoring of seagrass habitats in ongoing management programs.

The origin of a large (>3km) maar volcano by coalescence of multiple shallow craters: Lake Purrumbete maar, southeastern Australia

Lake Purrumbete maar is located in the intraplate, monogenetic Newer Volcanics Province in southeastern Australia. The extremely large crater of 3000. m in diameter formed on an intersection of two fault lines and comprises at least three coalesced vents. The evolution of these vents is controlled by the interaction of the tectonic setting and the properties of both hard and soft rock aquifers. Lithics in the maar deposits originate from country rock formations less than 300. m deep, indicating that the large size of the crater cannot only be the result of the downwards migration of the explosion foci in a single vent. Vertical crater walls and primary inward dipping beds evidence that the original size of the crater has been largely preserved. Detailed mapping of the facies distributions, the direction of transport of base surges and pyroclastic flows, and the distribution of ballistic block fields, form the basis for the reconstruction of the complex eruption history,which is characterised by alternations of the eruption style between relatively dry and wet phreatomagmatic conditions, and migration of the vent location along tectonic structures. Three temporally separated eruption phases are recognised, each starting at the same crater located directly at the intersection of two local fault lines. Activity then moved quickly to different locations. A significant volcanic hiatus between two of the three phases shows that the magmatic system was reactivated. The enlargement of especially the main crater by both lateral and vertical growth led to the interception of the individual craters and the formation of the large circular crater. Lake Purrumbete maar is an excellent example of how complicated the evolution of large, seemingly simple, circular maar volcanoes can be, and raises the question if these systems are actually monogenetic.

Estimate of Lagrangian integral scales in shallow tidal water using high resolution GPS-tracked drifters

In estuaries and natural water channels, the estimate of velocity and dispersion coefficients is critical to the knowledge of scalar transport and mixing. This estimate is rarely available experimentally at sub-tidal time scale in shallow water channels where high frequency is required to capture its spatio-temporal variation. This study estimates Lagrangian integral scales and autocorrelation curves, which are key parameters for obtaining velocity fluctuations and dispersion coefficients, and their spatio-temporal variability from deployments of Lagrangian drifters sampled at 10 Hz for a 4-hour period. The power spectral densities of the velocities between 0.0001 and 0.8 Hz were well fitted with a slope of 5/3 predicted by Kolmogorov’s similarity hypothesis within the inertial subrange, and were similar to the Eulerian power spectral previously observed within the estuary. The result showed that large velocity fluctuations determine the magnitude of the integral time scale, TL. Overlapping of short segments improved the stability of the estimate of TL by taking advantage of the redundant data included in the autocorrelation function. The integral time scales were about 20 s and varied by up to a factor of 8. These results are essential inputs for spatial binning of velocities, Lagrangian stochastic modelling and single particle analysis of the tidal estuary.

Postglacial Early Permian (late Sakmarian-early Artinskian) shallow-marine carbonate deposition along a 2000 km transect from Timor to west Australia

Late Sakmarian to early Artinskian (Early Permian) carbonate deposition was widespread in the marine intracratonic rift basins that extended into the interior of Eastern Gondwana from Timor in the north to the northern Perth Basin in the south. These basins spanned about 20° of paleolatitude (approximately 35°S to 55°S). This study describes the type section of the Maubisse Limestone in Timor-Leste, and compares this unit with carbonate sections in the Canning Basin (Nura Nura Member of the Poole Sandstone), the Southern Carnarvon Basin (Callytharra Formation) and the northern Perth Basin (Fossil Cliff Member of the Holmwood Shale). The carbonate units have no glacial influence and formed part of a major depositional cycle that, in the southern basins, overlies glacially influenced strata and lies a short distance below mudstone containing marine fossils and scattered dropstones (perhaps indicative of sea ice). In the south marine conditions became more restricted and were replaced by coal measures at the top of the depositional sequence. In the north, the carbonate deposits are possibly bryozoan–crinoidal mounds; whereas in the southern basins they form laterally continuous relatively thin beds, deposited on a very low-gradient seafloor, at the tops of shale–limestone parasequences that thicken upward in parasequence sets. All marine deposition within the sequence took place under very shallow (inner neritic) conditions, and the limestones have similar grain composition. Bryozoan and crinoidal debris dominate the grain assemblages and brachiopod shell fragments, foraminifera and ostracod valves are usually common. Tubiphytes ranged as far south as the Southern Carnarvon Basin, albeit rarely, but is more common to the north. Gastropod and bivalve shell debris, echinoid spines, solitary rugose corals and trilobite carapace elements are rare. The uniformity of the grain assemblage and the lack of tropical elements such as larger fusulinid foraminifera, colonial corals or dasycladacean algae indicate temperate marine conditions with only a small increase in temperature to the north. The depositional cycle containing the studied carbonate deposits represents a warmer phase than the preceding glacially influenced Asselian to early Sakmarian interval and the subsequent cool phase of the “mid” Artinskian that is followed by significant warming during the late Artinskian–early Kungurian. The timing of cooler and warmer intervals in the west Australian basins seems out-of-phase with the eastern Australian succession, but this may be a problem of chronostratigraphic miscorrelation due to endemic faunas and palynofloras.

Turbulent mixed convection in a shallow enclosure with a series of heat generating components

Turbulent mixed convection flow and heat transfer in a shallow enclosure with and without partitions and with a series of block-like heat generating components is studied numerically for a range of Reynolds and Grashof numbers with a time-dependent formulation. The flow and temperature distributions are taken to be two-dimensional. Regions with the same velocity and temperature distributions can be identified assuming repeated placement of the blocks and fluid entry and exit openings at regular distances, neglecting the end wall effects. One half of such module is chosen as the computational domain taking into account the symmetry about the vertical centreline. The mixed convection inlet velocity is treated as the sum of forced and natural convection components, with the individual components delineated based on pressure drop across the enclosure. The Reynolds number is based on forced convection velocity. Turbulence computations are performed using the standard k– model and the Launder–Sharma low-Reynolds number k– model. The results show that higher Reynolds numbers tend to create a recirculation region of increasing strength in the core region and that the effect of buoyancy becomes insignificant beyond a Reynolds number of typically 5×105. The Euler number in turbulent flows is higher by about 30 per cent than that in the laminar regime. The dimensionless inlet velocity in pure natural convection varies as Gr1/3. Results are also presented for a number of quantities of interest such as the flow and temperature distributions, Nusselt number, pressure drop and the maximum dimensionless temperature in the block, along with correlations.

Simplified dispersion curves for circular cylindrical shells using shallow shell theory.

An alternative derivation of the dispersion relation for the transverse vibration of a circular cylindrical shell is presented. The use of the shallow shell theory model leads to a simpler derivation of the same result. Further, the applicability of the dispersion relation is extended to the axisymmetric mode and the high frequency beam mode.