New method for quantification of vuggy porosity from digital optical borehole images as applied to the karstic Pleistocene limestone of the Biscayne aquifer, southeastern Florida

K.J. Cunningham; J.I. Carlson; N.F. Hurley

doi:10.1016/j.jappgeo.2003.06.006

Home/Archive/USGS

USGS · 70027766

New method for quantification of vuggy porosity from digital optical borehole images as applied to the karstic Pleistocene limestone of the Biscayne aquifer, southeastern Florida

K.J. Cunningham·J.I. Carlson·N.F. Hurley

View original source ↗Export metadata ↓

Abstract

Vuggy porosity is gas- or fluid-filled openings in rock matrix that are large enough to be seen with the unaided eye. Well-connected vugs can form major conduits for flow of ground water, especially in carbonate rocks. This paper presents a new method for quantification of vuggy porosity calculated from digital borehole images collected from 47 test coreholes that penetrate the karstic Pleistocene limestone of the Biscayne aquifer, southeastern Florida. Basically, the method interprets vugs and background based on the grayscale color of each in digital borehole images and calculates a percentage of vuggy porosity. Development of the method was complicated because environmental conditions created an uneven grayscale contrast in the borehole images that makes it difficult to distinguish vugs from background. The irregular contrast was produced by unbalanced illumination of the borehole wall, which was a result of eccentering of the borehole-image logging tool. Experimentation showed that a simple, single grayscale threshold would not realistically differentiate between the grayscale contrast of vugs and background. Therefore, an equation was developed for an effective subtraction of the changing grayscale contrast, due to uneven illumination, to produce a grayscale threshold that successfully identifies vugs. In the equation, a moving average calculated around the circumference of the borehole and expressed as the background grayscale intensity is defined as a baseline from which to identify a grayscale threshold for vugs. A constant was derived empirically by calibration with vuggy porosity values derived from digital images of slabbed-core samples and used to make the subtraction from the background baseline to derive the vug grayscale threshold as a function of azimuth. The method should be effective in estimating vuggy porosity in any carbonate aquifer. ?? 2003 Published by Elsevier B.V.

Explore related subjects

Journal of Applied Geophysics Article

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

K.J. Cunningham, J.I. Carlson, N.F. Hurley. 2004. New method for quantification of vuggy porosity from digital optical borehole images as applied to the karstic Pleistocene limestone of the Biscayne aquifer, southeastern Florida. https://doi.org/10.1016/j.jappgeo.2003.06.006

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2024

Imaging freshwater and saline aquifers beneath Bradford County, Pennsylvania, USA, using Audio-Magnetotelluric (AMT) data

We investigate the utility of audio-magnetotelluric (AMT) data for detecting high electrical conductivity zones indicative of fresh and saline groundwater lenses at depths of < ∼1000 m at seven survey sites in the southwestern corner of Bradford County, PA. Dimensionality analysis indicates the AMT data can be interpreted using two-dimensional (2-D) models. 2-D electrical resistivity models obtained from inverting the AMT data reveal conductive zones at depths of <30 m at all sites, between depths of 50–150 m at six sites, and between depths of 400–600 m at one site. The models show slightly lower resistivities (∼15–30 Ω m) in the deeper zones compared to the shallower zones (>30 Ω m). Consistent with well log data from two boreholes in the study area and salinity estimates obtained from the imaged resistivities, we interpret the <30 m deep conductive zone as the freshwater aquifer, and the two deeper conductive zones as salty/briny groundwater lenses. Our salinity estimates for the deeper zones vary from 1 to 380 (unitless Practical Salinity Scale), which agrees well with previous estimates of Appalachian brine salinity of 10–343. We attribute the high-salinity groundwater lenses to Appalachian Basin brine migrating upwards from deeper in the basin via bedrock fracture and fault systems. Our findings indicate AMT surveying can be useful for imaging fresh and saline aquifers at shallow depths (< ∼1000 m), at least within the northeastern section of the Appalachian Basin in Pennsylvania.

Pennsylvania↗

USGS2023

Determining seasonal recharge, storage changes, and specific yield using repeat microgravity and water-level measurements in the Mesilla Basin alluvial aquifer, New Mexico, 2016–2018

Increasing water demand and multi-year drought conditions within the Mesilla/Conejos-Médanos Basin near the New Mexico-Texas- Chihuahua border have resulted in diminished surface-water supplies and increased groundwater withdrawals. To better understand recharge to the shallow aquifer, the spatial and temporal groundwater storage changes, and the variability of specific yield (Sy) in the aquifer, seasonal groundwater elevation and repeat microgravity measurements were made during the irrigation release and non-release seasons of 2016, 2017, and 2018 at a network of locations near Las Cruces, New Mexico. The data collected during this investigation were able to capture seasonal change in groundwater elevations and storage from various sources of recharge at multiple sites in the shallow aquifer. Seasonal recharge in the study area was attributed to streamflow, the application and conveyance of irrigation water, and large or sustained precipitation events. However, increasing groundwater gradients in recent decades between piezometers close to the river and those more than a kilometer from the river suggests that recharge from river seepage has become localized at the seasonal scale. Overall, there was a net increase in storage of almost 8.4 cubic hectometers in the study reach between the start and end of the study, largely following the increased surface-water availability and above average precipitation in 2017. Specific yield, estimated by comparing the groundwater-level changes and storage changes at six sites in the study area, ranged from 0.14 (+/− 0.05) to 0.30 (+/− 0.06), which is slightly greater than previously reported estimates (0.10 to 0.25), but still within the error of the estimates. Most of the variability in the estimated storage change, that was not well-correlated with groundwater elevation change, is thought to be from soil moisture in the unsaturated zone. This investigation demonstrates the value of adding repeat microgravity measurements to conventional groundwater monitoring to better understand the sources and extent of recharge as well as the variability of Sy in the aquifer. Continued monitoring, under a variety of available surface water and meteorological conditions, could provide a more comprehensive understanding of the water budget and reduce the specific yield estimation uncertainty. Evaluating water-levels and storage conditions prior to, and following, local recharge events may help managers identify threshold conditions for aquifer storage depletions and recoveries.

New Mexico↗

USGS2023

Using DC resistivity ring array surveys to resolve conductive structures around tunnels or mine-workings

In underground environments, conventional direct current (DC) resistivity surveys with a single linear array of electrodes produce fundamentally non-unique inversions. These non-uniqueness and model resolution issues stem from limitations placed on the location of transmitters (TXs) and receivers (RXs) by the geometry of existing tunnels and boreholes. Poor excitation and/or sampling of the region of interest (ROI) can create artifacts and reduce the resolution of the recovered model. To address these problems we propose the use of an ensemble of ring arrays, which are created by placing one or more electrodes in each face (sidewalls, floor, and ceiling) of the tunnel to form a ring of electrodes at each along-tunnel location. Using a series of increasingly complex synthetic models, we assess the benefits of ring arrays and show that they can be used to better constrain the location and shape of anomalous bodies around the tunnel. Although ring arrays significantly improve the resolution of the recovered model, the size of the comprehensive ring array survey increases rapidly with the number of electrodes used. To balance model resolution and survey size, we developed a physics-based survey design methodology. In this methodology, TXs are selected based upon secondary charge accumulations on a test block that is moved throughout the ROI. Although this survey design methodology does not produce a strictly optimal survey, it balances model resolution and survey size in a practical and computationally efficient manner. Since the ring array more accurately estimates the around-tunnel location of targets and ensures that targets on all sides of the tunnel are detected, it is ideally suited to tunnel-based environments. Our results show that only about 6 % "> 6% of the possible TXs and 0.5 % "> 0.5% of the RXs in the comprehensive ring array survey are needed to retain the improvements in resolution. Therefore, economical ring array surveys can be designed for both reconnaissance and target characterization. Following the inversion of the reconnaissance dataset, additional rings can be added to reduce the inter-ring spacing or off-tunnel boreholes can be added to the region around identified anomalies to increase resolution as required.

Journal of Applied Geophysics↗