SEARCH · Search USGS

Results for “Computers and 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.

143 records · Page 8Linked to original sources

USGS2005

Porphyry copper deposits of the world: database, map, and grade and tonnage models

Mineral deposit models are important in exploration planning and quantitative resource assessments for two reasons: (1) grades and tonnages among deposit types are significantly different, and (2) many types occur in different geologic settings that can be identified from geologic maps. Mineral deposit models are the keystone in combining the diverse geoscience information on geology, mineral occurrences, geophysics, and geochemistry used in resource assessments and mineral exploration. Too few thoroughly explored mineral deposits are available in most local areas for reliable identification of the important geoscience variables or for robust estimation of undiscovered deposits-thus we need mineral-deposit models. Globally based deposit models allow recognition of important features because the global models demonstrate how common different features are. Well-designed and -constructed deposit models allow geologists to know from observed geologic environments the possible mineral deposit types that might exist, and allow economists to determine the possible economic viability of these resources in the region. Thus, mineral deposit models play the central role in transforming geoscience information to a form useful to policy makers. The foundation of mineral deposit models is information about known deposits-the purpose of this publication is to make this kind of information available in digital form for porphyry copper deposits. This report is an update of an earlier publication about porphyry copper deposits. In this report we have added 84 new porphyry copper deposits and removed 12 deposits. In addition, some errors have been corrected and a number of deposits have had some information, such as grades, tonnages, locations, or ages revised. This publication contains a computer file of information on porphyry copper deposits from around the world. It also presents new grade and tonnage models for porphyry copper deposits and for three subtypes of porphyry copper deposits and a map showing the location of all deposits. The value of this information and any derived analyses depends critically on the consistent manner of data gathering. For this reason, we first discuss the rules used in this compilation. Next, the fields of the data file are considered. Finally, we provide new grade and tonnage models.

Open-File Report↗

USGS2000

The 1999 Southern California Seismic Network bulletin

The Pasadena office of the U.S. Geological Survey (USGS), together with the Caltech Seismological Laboratory, operates a network of more than 350 remote seismometers in southern California called the S outhern C alifornia S eismic N etwork (SCSN). SCSN is part of TriNet, a cooperative project between the USGS, Caltech, and the California Division of Mines and Geology (CDMG). The TriNet project is halfway completed and is upgrading the existing network to digital, adding new stations, and developing real-time and earthquake-alert capabilities. Signals from the SCSN sites are telemetered to a central processing location at the Caltech Seismological Lab in Pasadena. Computers that detect and record thousands of earthquakes each year continuously monitor these signals. Phase arrival times for these events are picked by analysts and are archived along with digital seismograms. Data acquisition, processing, and archiving are achieved using the Caltech/USGS Seismic Processing (CUSP) system ( Dollar, 1989 ). These data have been compiled into the SCSN Catalog of Earthquakes, a list beginning in 1932 that currently contains more than 356,300 events. Waveform, phase, and catalog data are archived by the Southern California Earthquake Center Data Center (SCEC_DC). This data set is critical to the evaluation of earthquake hazards in California and to the advancement of geoscience as a whole.

California↗

USGS1992

Toward digital geologic map standards: a progress report

Establishing modern scientific and technical standards for geologic maps and their derivative map products is vital to both producers and users of such maps as we move into an age of digital cartography. Application of earth-science data in complex geographic information systems, acceleration of geologic map production, and reduction of population costs require that national standards be developed for digital geologic cartography and computer analysis. Since December 1988, under commission of the Chief Geologic of the U.S. Geological Survey and the mandate of the National Geologic Mapping Program (with added representation from the Association of American State Geologists), a committee has been designing a comprehensive set of scientific map standards. Three primary issues were: (1) selecting scientific symbology and its digital representation; (2) creating an appropriate digital coding system that characterizes geologic features with respect to their physical properties, stratigraphic and structural relations, spatial orientation, and interpreted mode of origin; and (3) developing mechanisms for reporting levels of certainty for descriptive as well as measured properties. Approximately 650 symbols for geoscience maps, including present usage of the U.S Geological Survey, state geological surveys, industry, and academia have been identified and tentatively adopted. A proposed coding system comprises four-character groupings of major and minor codes that can identify all attributes of a geologic feature. Such a coding system allows unique identification of as many as 105 geologic names and values on a given map. The new standard will track closely the latest developments of the Proposed Standard for Digital Cartographic Data soon to be submitted to the National Institute of Standards and Technology by the Federal Interagency Coordinating Committee on Digital Cartography. This standard will adhere generally to the accepted definitions and specifications for spatial data transfer. It will require separate specifications of digital cartographic quality relating to positional accuracy and ranges of measured and interpreted values such as geologic age and rock composition. Provisional digital geologic map standards will be published for trial implementation. After approximately two years, when comments on the proposed standards have been solicited and modifications made, formal adoption of the standards will be recommended. Widespread acceptance of the new standards will depend on their applicability to the broadest range of earth-science map products and their adaptability to changing cartographic technology.

Conference Paper↗

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↗

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

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↗