SEARCH · USGS Science

Results for “Computers & Geosciences”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

178 records · Page 10Linked to original sources

USGS2018

Geometric versus anemometric surface roughness for a shallow accumulating snowpack

When applied to a snow-covered surface, aerodynamic roughness length, z 0 , is typically considered as a static parameter within energy balance equations. However, field observations show that z 0 changes spatially and temporally, and thus z 0 incorporated as a dynamic parameter may greatly improve models. To evaluate methods for characterizing snow surface roughness, we compared concurrent estimates of z 0 based on (1) terrestrial light detection and ranging derived surface geometry of the snowpack surface (geometric, z 0 G ) and (2) vertical wind profile measurements (anemometric, z 0 A ). The value of z 0 G was computed from Lettau’s equation and underestimated z 0 A but compared well when scaled by a factor of 2.34. The Counihan method for computing z 0 G was found to be unsuitable for estimating z 0 on a snow surface. During snowpack accumulation in early winter, z 0 varied as a function of the snow-covered area (SCA). Our results show that as the SCA increases, z 0 decreases, indicating there is a topographic influence on this relation.

Geosciences↗

USGS2003

Rockfall hazard and risk assessment in the Yosemite Valley, California, USA

Rock slides and rock falls are the most frequent types of slope movements in Yosemite National Park, California. In historical time (1857-2002) 392 rock falls and rock slides have been documented in the valley, and some of them have been mapped in detail. We present the results of an attempt to assess rock fall hazards in the Yosemite Valley. Spatial and temporal aspects of rock falls hazard are considered. A detailed inventory of slope movements covering the 145-year period from 1857 to 2002 is used to determine the frequency-volume statistics of rock falls and to estimate the annual frequency of rock falls, providing the temporal component of rock fall hazard. The extent of the areas potentially subject to rock fall hazards in the Yosemite Valley were obtained using STONE, a physically-based rock fall simulation computer program. The software computes 3-dimensional rock fall trajectories starting from a digital elevation model (DEM), the location of rock fall release points, and maps of the dynamic rolling friction coefficient and of the coefficients of normal and tangential energy restitution. For each DEM cell the software calculates the number of rock falls passing through the cell, the maximum rock fall velocity and the maximum flying height. For the Yosemite Valley, a DEM with a ground resolution of 10 ?? 10 m was prepared using topographic contour lines from the U.S. Geological Survey 1:24 000-scale maps. Rock fall release points were identified as DEM cells having a slope steeper than 60??, an assumption based on the location of historical rock falls. Maps of the normal and tangential energy restitution coefficients and of the rolling friction coefficient were produced from a surficial geologic map. The availability of historical rock falls mapped in detail allowed us to check the computer program performance and to calibrate the model parameters. Visual and statistical comparison of the model results with the mapped rock falls confirmed the accuracy of the model. The model results are compared with a previous map of rockfall talus and with a geomorphic assessment of rock fall hazard based on potential energy referred to as a shadow angle approach, recently completed for the Yosemite Valley. The model results are then used to identify the roads and trails more subject to rock fall hazard. Of the 166.5 km of roads and trails in the Yosemite Valley 31.2% were found to be potentially subject to rock fall hazard, of which 14% are subject to very high hazard. ?? European Geosciences Union 2003.

Natural Hazards and Earth System Sciences↗

USGS2016

Dynamic reusable workflows for ocean science

Digital catalogs of ocean data have been available for decades, but advances in standardized services and software for catalog search and data access make it now possible to create catalog-driven workflows that automate — end-to-end — data search, analysis and visualization of data from multiple distributed sources. Further, these workflows may be shared, reused and adapted with ease. Here we describe a workflow developed within the US Integrated Ocean Observing System (IOOS) which automates the skill-assessment of water temperature forecasts from multiple ocean forecast models, allowing improved forecast products to be delivered for an open water swim event. A series of Jupyter Notebooks are used to capture and document the end-to-end workflow using a collection of Python tools that facilitate working with standardized catalog and data services. The workflow first searches a catalog of metadata using the Open Geospatial Consortium (OGC) Catalog Service for the Web (CSW), then accesses data service endpoints found in the metadata records using the OGC Sensor Observation Service (SOS) for in situ sensor data and OPeNDAP services for remotely-sensed and model data. Skill metrics are computed and time series comparisons of forecast model and observed data are displayed interactively, leveraging the capabilities of modern web browsers. The resulting workflow not only solves a challenging specific problem, but highlights the benefits of dynamic, reusable workflows in general. These workflows adapt as new data enters the data system, facilitate reproducible science, provide templates from which new scientific workflows can be developed, and encourage data providers to use standardized services. As applied to the ocean swim event, the workflow exposed problems with two of the ocean forecast products which led to improved regional forecasts once errors were corrected. While the example is specific, the approach is general, and we hope to see increased use of dynamic notebooks across the geoscience domains.

Journal of Marine Science and Engineering↗

USGS1987

Microprocessor-based data-acquisition system for a borehole radar

An efficient microprocessor-based system has been implemented that permits real-time acquisition, stacking, and digital recording of data generated by a borehole radar system. Although the system digitizes, stacks, and records independently of a computer, it is interfaced to a desktop computer for program control over system parameters such as sampling interval, number of samples, number of times the data are stacked prior to recording on nine-track tape, and for graphics display of the digitized data. The data can be transferred to the desktop computer during recording, or played back from a tape at a later time. Using the desktop computer, the operator can observe results while recording data and generate hard-copy graphics in the field. Thus, the radar operator can immediately evaluate the quality of data being obtained, modify system parameters, study the radar logs before leaving the field, and re-run borehole logs if necessary. The system has proven to be reliable in the field and has increased productivity both in the field and in the laboratory.

IEEE Transactions on Geoscience and Remote Sensing↗

USGS2009

A new numerical method for calculating extrema of received power for polarimetric SAR

A numerical method called cross-step iteration is proposed to calculate the maximal/minimal received power for polarized imagery based on a target's Kennaugh matrix. This method is much more efficient than the systematic method, which searches for the extrema of received power by varying the polarization ellipse angles of receiving and transmitting polarizations. It is also more advantageous than the Schuler method, which has been adopted by the PolSARPro package, because the cross-step iteration method requires less computation time and can derive both the maximal and minimal received powers, whereas the Schuler method is designed to work out only the maximal received power. The analytical model of received-power optimization indicates that the first eigenvalue of the Kennaugh matrix is the supremum of the maximal received power. The difference between these two parameters reflects the depolarization effect of the target's backscattering, which might be useful for target discrimination. ?? 2009 IEEE.

IEEE Geoscience and Remote Sensing Letters↗

USGS2008

Numerical modeling of higher order magnetic moments in UXO discrimination

The surface magnetic anomaly observed in unexploded ordnance (UXO) clearance is mainly dipolar, and consequently, the dipole is the only magnetic moment regularly recovered in UXO discrimination. The dipole moment contains information about the intensity of magnetization but lacks information about the shape of the target. In contrast, higher order moments, such as quadrupole and octupole, encode asymmetry properties of the magnetization distribution within the buried targets. In order to improve our understanding of magnetization distribution within UXO and non-UXO objects and to show its potential utility in UXO clearance, we present a numerical modeling study of UXO and related metallic objects. The tool for the modeling is a nonlinear integral equation describing magnetization within isolated compact objects of high susceptibility. A solution for magnetization distribution then allows us to compute the magnetic multipole moments of the object, analyze their relationships, and provide a depiction of the anomaly produced by different moments within the object. Our modeling results show the presence of significant higher order moments for more asymmetric objects, and the fields of these higher order moments are well above the noise level of magnetic gradient data. The contribution from higher order moments may provide a practical tool for improved UXO discrimination. ?? 2008 IEEE.

IEEE Transactions on Geoscience and Remote Sensing↗

USGS2026

The path to FAIR research models: Lessons learned

Numerical modeling of Earth surface processes emerged as an important scientific tool in the late 1960s to mid-1970s, driven by the development of finite element methods in computer science. These advancements, initially applied in civil engineering, enabled scientists to simulate complex geological phenomena. At that time, models were often only described in publications, access was limited to researchers with direct connections to the developers, and the code was rarely documented for reuse, limiting their application beyond the original research context. The FAIR principles (Findability, Accessibility, Interoperability, and Reusability) as applied to data began to take shape in the 21st century with the rise of open science, digital repositories, and standardized data sharing frameworks. In the late 2010s, grassroots movements began to apply some of the FAIRness goals to numerical models. Subsequently, more formalized FAIR model principles were developed that addressed the specific needs of the scientific modeling community, resulting in the formulation of the FAIR principles for research software (FAIR4RS). In this study, we examine the development and implementation of strategies by two geoscience research infrastructures – the CSDMS (Community Surface Dynamics Modeling System) Model Repository and the U.S. Geological Survey Model Catalog – to enhance the FAIRness of models guided by FAIR4RS. Some of the development and implementation efforts described predate the formalization of FAIR and FAIR4RS principles, making this an ongoing and adaptive process. We evaluate the temporal progression towards increased FAIR4RS alignment across three phases of research infrastructure development: prototype, refinement, and growth and iteration. Although certain principles were more straightforward to implement early in prototypes of the catalog infrastructures, others required broader community collaboration during refinement, and some continue to pose practical challenges in the growth and iteration phase. By tracing these dynamics, our aim is to provide insights that can guide other modeling initiatives in effectively adopting FAIR4RS principles within their communities.

Geoscientific Model Development↗

USGS2023

Differentiable modelling to unify machine learning and physical models for geosciences

Process-based modelling offers interpretability and physical consistency in many domains of geosciences but struggles to leverage large datasets efficiently. Machine-learning methods, especially deep networks, have strong predictive skills yet are unable to answer specific scientific questions. In this Perspective, we explore differentiable modelling as a pathway to dissolve the perceived barrier between process-based modelling and machine learning in the geosciences and demonstrate its potential with examples from hydrological modelling. ‘Differentiable’ refers to accurately and efficiently calculating gradients with respect to model variables or parameters, enabling the discovery of high-dimensional unknown relationships. Differentiable modelling involves connecting (flexible amounts of) prior physical knowledge to neural networks, pushing the boundary of physics-informed machine learning. It offers better interpretability, generalizability, and extrapolation capabilities than purely data-driven machine learning, achieving a similar level of accuracy while requiring less training data. Additionally, the performance and efficiency of differentiable models scale well with increasing data volumes. Under data-scarce scenarios, differentiable models have outperformed machine-learning models in producing short-term dynamics and decadal-scale trends owing to the imposed physical constraints. Differentiable modelling approaches are primed to enable geoscientists to ask questions, test hypotheses, and discover unrecognized physical relationships. Future work should address computational challenges, reduce uncertainty, and verify the physical significance of outputs.

Nature Reviews Earth & Environment↗

USGS2018

Karachi tides during the 1945 Makran tsunami

This paper extends and detides a Karachi tide-gauge record as an observational basis for assessing Indian Ocean tsunami risk. The extended marigram encompasses the time of the great 1945 Makran earthquake of early November 28, local time, and of the ensuing tsunami, which continued into November 29. The marigram was published previously as a 9-h excerpt that begins 1 h after the earthquake. The full marigram presented here covers most of 17 days from November 15 to December 1. Gaps include a tsunami-induced outage that may help explain why the highest water level gauged is 1 m below the maximum water level reported nearby. The detiding method computes a reference tidal curve that disregards all observations from November 28 and 29. For those 2 days, the reference tide is guided by Admiralty tide tables and, secondarily, by high waters and low waters gauged before and after. As in previous estimates, the detided tsunami crests about 0.5 m above ambient tide, but now with the possibility that the gauge failed to record a higher wave. Anomalies described for the first time include an early one that likely resulted from a recognized problem with the Karachi tide station, but which might instead represent an earthquake precursor.

Karachi↗

USGS2021

Implications of aggregating and smoothing daily production data on estimates of the transition time between flow regimes in horizontal hydraulically fractured Bakken oil wells

The level to which data are aggregated or smoothed can impact analytical and predictive modeling results. This paper discusses findings regarding such impacts on estimating change points in production flow regimes of horizontal hydraulically fractured shale oil wells producing from the middle member of the Bakken Formation. Change points that signal transitions in flow regimes are important because they subsequently affect estimates of ultimate recovery from wells producing from shale plays. Extending our earlier work, we employ two different statistical approaches, Bacon–Watts Bayesian regression and nonlinear constrained least squares regression, and a designed computational experiment to estimate the time of transition from the transient to the boundary-dominated flow regime for 14 different wells using daily production data rather than aggregated monthly data, as previously considered. The daily data were also smoothed to reduce noise. Computational experiments suggest that both statistical approaches can lead to plausible estimates of the transition point under different data aggregation or smoothing regimes, but that daily data are likely too granular to produce credible estimates. Although the expected value of transition points using smoothed daily data and monthly disaggregated data are generally comparable, the confidence intervals bounding the estimates based on smoothed daily data are generally wider. Our results not only inform the operational practices of oil producers engaged in economic evaluation of their shale resources and additional play development activities, but also the activities of petroleum research groups, government agencies, and financial organizations seeking to improve the trustworthiness of resource projections.

Mathematical Geosciences↗

USGS2011

Liquefaction probability curves for surficial geologic deposits

Liquefaction probability curves that predict the probability of surface manifestations of earthquake-induced liquefaction are developed for 14 different types of surficial geologic units. The units consist of alluvial fan, beach ridge, river delta topset and foreset beds, eolian dune, point bar, flood basin, natural river and alluvial fan levees, abandoned river channel, deep-water lake, lagoonal, sandy artificial fill, and valley train deposits. Probability is conditioned on earthquake magnitude and peak ground acceleration. Curves are developed for water table depths of 1.5 and 5.0 m. Probabilities are derived from complementary cumulative frequency distributions of the liquefaction potential index (LPI) that were computed from 927 cone penetration tests. For natural deposits with a water table at 1.5 m and subjected to a M7.5 earthquake with peak ground acceleration (PGA) = 0.25 g , probabilities range from <0.03 for alluvial fan and lacustrine deposits to >0.5 for beach ridge, point bar, and deltaic deposits. The curves also were used to assign ranges of liquefaction probabilities to the susceptibility categories proposed previously for different geologic deposits. For the earthquake described here, probabilities for susceptibility categories have ranges of 0–0.08 for low, 0.09–0.30 for moderate, 0.31–0.62 for high, and 0.63–1.00 for very high. Retrospective predictions of liquefaction during historical earthquakes based on the curves compare favorably to observations.

Environmental & Engineering Geoscience↗

USGS2015

Interactive access to LP DAAC satellite data archives through a combination of open-source and custom middleware web services

Current methods of searching for and retrieving data from satellite land remote sensing archives do not allow for interactive information extraction. Instead, Earth science data users are required to download files over low-bandwidth networks to local workstations and process data before science questions can be addressed. New methods of extracting information from data archives need to become more interactive to meet user demands for deriving increasingly complex information from rapidly expanding archives. Moving the tools required for processing data to computer systems of data providers, and away from systems of the data consumer, can improve turnaround times for data processing workflows. The implementation of middleware services was used to provide interactive access to archive data. The goal of this middleware services development is to enable Earth science data users to access remote sensing archives for immediate answers to science questions instead of links to large volumes of data to download and process. Exposing data and metadata to web-based services enables machine-driven queries and data interaction. Also, product quality information can be integrated to enable additional filtering and sub-setting. Only the reduced content required to complete an analysis is then transferred to the user.

IEEE Geoscience and Remote Sensing Letters↗

USGS2009

Atmospheric correction at AERONET locations: A new science and validation data set

This paper describes an Aerosol Robotic Network (AERONET)-based Surface Reflectance Validation Network (ASRVN) and its data set of spectral surface bidirectional reflectance and albedo based on Moderate Resolution Imaging Spectroradiometer (MODIS) TERRA and AQUA data. The ASRVN is an operational data collection and processing system. It receives 50 ?? 50 km 2 ; subsets of MODIS level 1B (L1B) data from MODIS adaptive processing system and AERONET aerosol and water-vapor information. Then, it performs an atmospheric correction (AC) for about 100 AERONET sites based on accurate radiative-transfer theory with complex quality control of the input data. The ASRVN processing software consists of an L1B data gridding algorithm, a new cloud-mask (CM) algorithm based on a time-series analysis, and an AC algorithm using ancillary AERONET aerosol and water-vapor data. The AC is achieved by fitting the MODIS top-of-atmosphere measurements, accumulated for a 16-day interval, with theoretical reflectance parameterized in terms of the coefficients of the Li SparseRoss Thick (LSRT) model of the bidirectional reflectance factor (BRF). The ASRVN takes several steps to ensure high quality of results: 1) the filtering of opaque clouds by a CM algorithm; 2) the development of an aerosol filter to filter residual semitransparent and subpixel clouds, as well as cases with high inhomogeneity of aerosols in the processing area; 3) imposing the requirement of the consistency of the new solution with previously retrieved BRF and albedo; 4) rapid adjustment of the 16-day retrieval to the surface changes using the last day of measurements; and 5) development of a seasonal backup spectral BRF database to increase data coverage. The ASRVN provides a gapless or near-gapless coverage for the processing area. The gaps, caused by clouds, are filled most naturally with the latest solution for a given pixel. The ASRVN products include three parameters of the LSRT model (kL, kG, and kV), surface albedo, normalized BRF (computed for a standard viewing geometry, VZA = 0, SZA = 45??), and instantaneous BRF (or one-angle BRF value derived from the last day of MODIS measurement for specific viewing geometry) for the MODIS 500-m bands 17. The results are produced daily at a resolution of 1 km in gridded format. We also provide a cloud mask, a quality flag, and a browse bitmap image. The ASRVN data set, including 6 years of MODIS TERRA and 1.5 years of MODIS AQUA data, is available now as a standard MODIS product (MODASRVN) which can be accessed through the Level 1 and Atmosphere Archive and Distribution System website ( http://ladsweb.nascom.nasa.gov/data/search.html). It can be used for a wide range of applications including validation analysis and science research. ?? 2006 IEEE.

IEEE Transactions on Geoscience and Remote Sensing↗

USGS2012

Archive of digital chirp subbottom profile data collected during USGS cruise 10BIM04 offshore Cat Island, Mississippi, September 2010

In September of 2010, the U.S. Geological Survey (USGS), in cooperation with the U.S. Army Corps of Engineers (USACE), conducted a geophysical survey to investigate the geologic controls on barrier island framework of Cat Island, Miss., as part of a broader USGS study on Barrier Island Mapping (BIM). These surveys were funded through the Mississippi Coastal Improvements Program (MsCIP) and the Northern Gulf of Mexico (NGOM) Ecosystem Change and Hazard Susceptibility Project as part of the Holocene Coastal Evolution of the Mississippi-Alabama Region Subtask. This report serves as an archive of unprocessed digital chirp subbottom data, trackline maps, navigation files, GIS files, Field Activity Collection System (FACS) logs, and formal FGDC metadata. Gained (showing a relative increase in signal amplitude) digital images of the seismic profiles are also provided. Refer to the Acronyms page for expansions of acronyms and abbreviations used in this report. The USGS Saint Petersburg Coastal and Marine Science Center (SPCMSC) assigns a unique identifier to each cruise or field activity. For example, 10BIM04 tells us the data were collected in 2010 during the fourth field activity for that project in that calendar year. Refer to http://walrus.wr.usgs.gov/infobank/programs/html/definition/activity.html for a detailed description of the method used to assign the field activity identification (ID). All chirp systems use a signal of continuously varying frequency; the EdgeTech SB-512i system used during this survey produces high-resolution, shallow-penetration (typically less than 50 milliseconds (ms)) profile images of sub-seafloor stratigraphy. The towfish contains a transducer that transmits and receives acoustic energy; it was housed within a float system (built at the SPCMSC), which allows the towfish to be towed at a constant depth of 1.07 meters (m) below the sea surface. As transmitted acoustic energy intersects density boundaries, such as the seafloor or sub-surface sediment layers, some energy is reflected back toward the transducer, received, and recorded by a Personal Computer (PC)-based seismic acquisition system. This process is repeated at regular time intervals (for example, 0.125 seconds (s)), and returned energy is recorded for a specific duration (for example, 50 ms). In this way, a two-dimensional (2-D) vertical image of the shallow geologic structure beneath the ship track is produced. Figure 1 displays the acquisition geometry. Refer to table 1 for a summary of acquisition parameters and table 2 for trackline statistics. The archived trace data are in standard Society of Exploration Geophysicists (SEG) SEG Y rev. 0 format (Barry and others, 1975); the first 3,200 bytes of the card image header are in American Standard Code for Information Interchange (ASCII) format instead of Extended Binary Coded Decimal Interchange Code (EBCDIC) format. The SEG Y files may be downloaded and processed with commercial or public domain software such as Seismic Unix (SU) (Cohen and Stockwell, 2010). See the How To Download SEG Y Data page for download instructions. The printable profiles provided here are GIF images that were processed and gained using SU software, and they can be viewed from the Profiles page or from links located on the trackline maps; refer to the Software page for links to example SU processing scripts. The SEG Y files are available on the DVD version of this report or on the Web, downloadable via the USGS Coastal and Marine Geoscience Data System (http://cmgds.marine.usgs.gov). The data are also available for viewing using GeoMapApp (http://www.geomapapp.org) and Virtual Ocean (http://www.virtualocean.org) multi-platform open source software.

Mississippi↗

USGS2013

Multitemporal cross-calibration of the Terra MODIS and Landsat 7 ETM+ reflective solar bands

In recent years, there has been a significant increase in the use of remotely sensed data to address global issues. With the open data policy, the data from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Enhanced Thematic Mapper Plus (ETM+) sensors have become a critical component of numerous applications. These two sensors have been operational for more than a decade, providing a rich archive of multispectral imagery for analysis of mutitemporal remote sensing data. This paper focuses on evaluating the radiometric calibration agreement between MODIS and ETM+ using the near-simultaneous and cloud-free image pairs over an African pseudo-invariant calibration site, Libya 4. To account for the combined uncertainties in the top-of-atmosphere (TOA) reflectance due to surface and atmospheric bidirectional reflectance distribution function (BRDF), a semiempirical BRDF model was adopted to normalize the TOA reflectance to the same illumination and viewing geometry. In addition, the spectra from the Earth Observing-1 (EO-1) Hyperion were used to compute spectral corrections between the corresponding MODIS and ETM+ spectral bands. As EO-1 Hyperion scenes were not available for all MODIS and ETM+ data pairs, MODerate resolution atmospheric TRANsmission (MODTRAN) 5.0 simulations were also used to adjust for differences due to the presence or lack of absorption features in some of the bands. A MODIS split-window algorithm provides the atmospheric water vapor column abundance during the overpasses for the MODTRAN simulations. Additionally, the column atmospheric water vapor content during the overpass was retrieved using the MODIS precipitable water vapor product. After performing these adjustments, the radiometric cross-calibration of the two sensors was consistent to within 7%. Some drifts in the response of the bands are evident, with MODIS band 3 being the largest of about 6% over 10 years, a change that will be corrected in Collection 6 MODIS processing.

IEEE Transactions on Geoscience and Remote Sensing↗

USGS2013

Time-lapse analysis of methane quantity in Mary Lee group of coal seams using filter-based multiple-point geostatistical simulation

Coal seam degasification and its success are important for controlling methane, and thus for the health and safety of coal miners. During the course of degasification, properties of coal seams change. Thus, the changes in coal reservoir conditions and in-place gas content as well as methane emission potential into mines should be evaluated by examining time-dependent changes and the presence of major heterogeneities and geological discontinuities in the field. In this work, time-lapsed reservoir and fluid storage properties of the New Castle coal seam, Mary Lee/Blue Creek seam, and Jagger seam of Black Warrior Basin, Alabama, were determined from gas and water production history matching and production forecasting of vertical degasification wellbores. These properties were combined with isotherm and other important data to compute gas-in-place (GIP) and its change with time at borehole locations. Time-lapsed training images (TIs) of GIP and GIP difference corresponding to each coal and date were generated by using these point-wise data and Voronoi decomposition on the TI grid, which included faults as discontinuities for expansion of Voronoi regions. Filter-based multiple-point geostatistical simulations, which were preferred in this study due to anisotropies and discontinuities in the area, were used to predict time-lapsed GIP distributions within the study area. Performed simulations were used for mapping spatial time-lapsed methane quantities as well as their uncertainties within the study area. The systematic approach presented in this paper is the first time in literature that history matching, TIs of GIPs and filter simulations are used for degasification performance evaluation and for assessing GIP for mining safety. Results from this study showed that using production history matching of coalbed methane wells to determine time-lapsed reservoir data could be used to compute spatial GIP and representative GIP TIs generated through Voronoi decomposition. Furthermore, performing filter simulations using point-wise data and TIs could be used to predict methane quantity in coal seams subjected to degasification. During the course of the study, it was shown that the material balance of gas produced by wellbores and the GIP reductions in coal seams predicted using filter simulations compared very well, showing the success of filter simulations for continuous variables in this case study. Quantitative results from filter simulations of GIP within the studied area briefly showed that GIP was reduced from an initial ∼73 Bcf (median) to ∼46 Bcf (2011), representing a 37 % decrease and varying spatially through degasification. It is forecasted that there will be an additional ∼2 Bcf reduction in methane quantity between 2011 and 2015. This study and presented results showed that the applied methodology and utilized techniques can be used to map GIP and its change within coal seams after degasification, which can further be used for ventilation design for methane control in coal mines.

Alabama↗