A distributed parallel multiple-relaxation-time lattice Boltzmann method on general-purpose graphics processing units for the rapid and scalable computation of absolute permeability from high-resolution 3D micro-CT images     

F. O. Alpak  F. Gray  N. Saxena  J. Dietderich  R. Hofmann  S. Berg 《Computational Geosciences》2018,22(3):815-832

Digital rock physics (DRP) is a rapidly evolving technology targeting fast turnaround times for repeatable core analysis and multi-physics simulation of rock properties. We develop and validate a rapid and scalable distributed-parallel single-phase pore-scale flow simulator for permeability estimation on real 3D pore-scale micro-CT images using a novel variant of the lattice Boltzmann method (LBM). The LBM code implementation is designed to take maximum advantage of distributed computing on multiple general-purpose graphics processing units (GPGPUs). We describe and extensively test the distributed parallel implementation of an innovative LBM algorithm for simulating flow in pore-scale media based on the multiple-relaxation-time (MRT) model that utilizes a precise treatment of body force. While the individual components of the resulting simulator can be separately found in various references, our novel contributions are (1) the integration of all of the mathematical and high-performance computing components together with a highly optimized code implementation and (2) the delivery of quantitative results with the simulator in terms of robustness, accuracy, and computational efficiency for a variety of flow geometries including various types of real rock images. We report on extensive validations of the simulator in terms of accuracy and provide near-ideal distributed parallel scalability results on large pore-scale image volumes that were largely computationally inaccessible prior to our implementation. We validate the accuracy of the MRT-LBM simulator on model geometries with analytical solutions. Permeability estimation results are then provided on large 3D binary microstructures including a sphere pack and rocks from various sandstone and carbonate formations. We quantify the scalability behavior of the distributed parallel implementation of MRT-LBM as a function of model type/size and the number of utilized GPGPUs for a panoply of permeability estimation problems.  相似文献   

A non-hydrostatic model for the numerical study of landslide-generated waves     

Pengfei Si  Jordan Aaron  Scott McDougall  Ji Lu  Xiping Yu  Nicholas J. Roberts  John J. Clague 《Landslides》2018,15(4):711-726

This paper proposes and demonstrates a two-layer depth-averaged model with non-hydrostatic pressure correction to simulate landslide-generated waves. Landslide (lower layer) and water (upper layer) motions are governed by the general shallow water equations derived from mass and momentum conservation laws. The landslide motion and wave generation/propagation are separately formulated, but they form a coupled system. Our model combines some features of the landslide analysis model DAN3D and the tsunami analysis model COMCOT and adds a non-hydrostatic pressure correction. We use the new model to simulate a 2007 rock avalanche-generated wave event at Chehalis Lake, British Columbia, Canada. The model results match both the observed distribution of the rock avalanche deposit in the lake and the wave run-up trimline along the shoreline. Sensitivity analyses demonstrate the importance of accounting for the non-hydrostatic dynamic pressure at the landslide-water interface, as well as the influence of the internal strength of the landslide on the size of the generated waves. Finally, we compare the numerical results of landslide-generated waves simulated with frictional and Voellmy rheologies. Similar maximum wave run-ups can be obtained using the two different rheologies, but the frictional model better reproduces the known limit of the rock avalanche deposit and is thus considered to yield the best overall results in this particular case.  相似文献   

Depositional setting,provenance, and tectonic-volcanic setting of Eocene–Recent deep-sea sediments of the oceanic Izu–Bonin forearc,northwest Pacific (IODP Expedition 352)     

Alastair H. F. Robertson  Steffen Kutterolf  Aaron Avery  Alan T. Baxter  Katerina Petronotis  Gary D. Acton 《International Geology Review》2018,60(15):1816-1854

New biostratigraphical, geochemical, and magnetic evidence is synthesized with IODP Expedition 352 shipboard results to understand the sedimentary and tectono-magmatic development of the Izu–Bonin outer forearc region. The oceanic basement of the Izu–Bonin forearc was created by supra-subduction zone seafloor spreading during early Eocene (c. 50–51 Ma). Seafloor spreading created an irregular seafloor topography on which talus locally accumulated. Oxide-rich sediments accumulated above the igneous basement by mixing of hydrothermal and pelagic sediment. Basaltic volcanism was followed by a hiatus of up to 15 million years as a result of topographic isolation or sediment bypassing. Variably tuffaceous deep-sea sediments were deposited during Oligocene to early Miocene and from mid-Miocene to Pleistocene. The sediments ponded into extensional fault-controlled basins, whereas condensed sediments accumulated on a local basement high. Oligocene nannofossil ooze accumulated together with felsic tuff that was mainly derived from the nearby Izu–Bonin arc. Accumulation of radiolarian-bearing mud, silty clay, and hydrogenous metal oxides beneath the carbonate compensation depth (CCD) characterized the early Miocene, followed by middle Miocene–Pleistocene increased carbonate preservation, deepened CCD and tephra input from both the oceanic Izu–Bonin arc and the continental margin Honshu arc. The Izu–Bonin forearc basement formed in a near-equatorial setting, with late Mesozoic arc remnants to the west. Subduction-initiation magmatism is likely to have taken place near a pre-existing continent–oceanic crust boundary. The Izu–Bonin arc migrated northward and clockwise to collide with Honshu by early Miocene, strongly influencing regional sedimentation.  相似文献   

Correction: Hydraulic tomography analysis of municipal-well operation data with geology-based groundwater models     

Tong  Xin  Illman  Walter A.  Berg  Steven J.  Luo  Ning 《Hydrogeology Journal》2021,29(5):1999-2000

Hydrogeology Journal - A Correction to this paper has been published: https://doi.org/10.1007/s10040-021-02351-x  相似文献   

Recharge response and kinematics of an unusual earthflow in Liechtenstein     

Aaron  Jordan  Loew  Simon  Forrer  Markus 《Landslides》2021,18(7):2383-2401

Landslides - Understanding landslide behavior over medium and long timescales is crucial for predicting landslide hazard and constructing accurate landscape evolution models. The behavior of...  相似文献   

An Automatic Well Planner for Complex Well Trajectories     

Kristoffersen  Brage S.  Bellout  Mathias C.  Silva  Thiago L.  Berg  Carl F. 《Mathematical Geosciences》2021,53(8):1881-1905

Mathematical Geosciences - A data-driven automatic well planner procedure is implemented to develop complex well trajectories by efficiently adapting to near-well reservoir properties and geometry....  相似文献   

A diverse Pleistocene marsupial trackway assemblage from the Victorian Volcanic Plains,Australia     

Stephen P. Carey  Aaron B. Camens  Matthew L. Cupper  Rainer Grün  John C. Hellstrom  Stafford W. McKnight  Iain Mclennan  David A. Pickering  Peter Trusler  Maxime Aubert 《Quaternary Science Reviews》2011,30(5-6):591-610

A diverse assemblage of late Pleistocene marsupial trackways on a lake bed in south-western Victoria provides the first information relating to the gaits and morphology of several megafaunal species, and represents the most speciose and best preserved megafaunal footprint site in Australia. The 60–110 ka volcaniclastic lacustrine sedimentary rocks preserve trackways of the diprotodontid Diprotodon optatum, a macropodid (probably Protemnodon sp.) and a large vombatid (perhaps Ramsayia magna or ‘Phascolomys’ medius) and possible prints of the marsupial lion, Thylacoleo carnifex. The footprints were imprinted within a short time period, demonstrating the association of the taxa present, rather than the time-averaged accumulations usually observed in skeletal fossil deposits. Individual manus and pes prints are distinguishable in some trackways, and in many cases some digital pad morphology is also present. Several parameters traditionally used to differentiate ichnotaxa, including trackway gauge and the degree of print in-turning relative to the midline, are shown to be subject to significant intraspecific variation in marsupials. Sexual dimorphism in the trackway proportions of Diprotodon, and its potential for occurrence in all large bodied, quadrupedal marsupials, is identified here for the first time.  相似文献   

Travel time controls the magnitude of nitrate discharge in groundwater bypassing the riparian zone to a stream on Virginia's coastal plain     

Samuel A. Flewelling  Janet S. Herman  George M. Hornberger  Aaron L. Mills 《水文研究》2012,26(8):1242-1253

Groundwater that bypasses the riparian zone by travelling along deep flow paths may deliver high concentrations of fertilizer‐derived NO₃^? to streams, or it may be impacted by the NO₃^? removal process of denitrification in streambed sediments. In a study of a small agricultural catchment on the Atlantic coastal plain of Virginia's eastern shore, we used seepage meters deployed in the streambed to measure specific discharge of groundwater and its solute concentrations for various locations and dates. We used values of Cl^? concentration to discriminate between bypass water recharged distal to the stream and that contained high NO₃^? but low Cl^? concentrations and riparian‐influenced water recharged proximal to the stream that contained low NO₃^? and high Cl^? concentrations. The travel time required for bypass water to transit the 30‐cm‐thick, microbially active denitrifying zone in the streambed determined the extent of NO₃^? removal, and hydraulic conductivity determined travel time through the streambed sediments. At all travel times greater than 2 days, NO₃^? removal was virtually complete. Comparison of the timescales for reaction and transport through the streambed sediments in this system confirmed that the predominant control on nitrate flux was travel time rather than denitrification rate coefficients. We conclude that extensive denitrification can occur in groundwater that bypasses the riparian zone, but a residence time in biologically active streambed sediments sufficient to remove a large fraction of the NO₃^? is only achieved in relatively low‐conductivity porous media. Instead of viewing them as separate, the streambed and riparian zone should be considered an integrated NO₃^? removal unit. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Tungsten skarn mineralizations in a regional metamorphic terrain in northern Norway: a possible metamorphic ore deposit   总被引：2，自引：0，他引：2  

R. Berg Larsen 《Mineralium Deposita》1991,26(4):281-289

Scheelite- and molybdenite-bearing scapolite skarn and plagioclase skarn occur as stratabound mineralizations in a terrain which has suffered regional metamorphism. Scapolite skarn formed as replacement skarns in plagioclase-bearing hornblendic gneiss whereas plagioclase skarn precipitated in fractures as vein skarn mineralizations. The genesis of these skarn deposits are closely connected to an episode of boudinage on the limbs of a large-scale dome-shaped fold. Fluid inclusion studies revealed that the ore-forming fluids were composed of CO₂ with 2 to 17 mole% CH₄, which formed due to decarbonation reactions during peak metamorphism at the end of the Caledonian orogenesis. Thermodynamic considerations indicate that scapolite and plagioclase formed above 550 °C and 2 Kbar pressure and were followed by precipitation of scheelite between 400 and 570 °C.  相似文献   

Correlation of foreshocks and aftershocks and asperities     

Vindell Hsu  Charles E. Helsley  Eduard Berg  David A. Novelo-Casanova 《Pure and Applied Geophysics》1984,122(6):878-893

A close correlation in spatial distribution of local seismic activity and energy release patterns before and after the 1979 Petatlan, Mexico earthquake suggests heterogeneity within the fault plane of this major low-angle thrust event associated with subduction along the Middle America Trench. A simple two-asperity model is proposed to account for the complexity. Foreshocks and aftershocks of the neighboring 1981 Playa Azul earthquake showed a similar pattern. As both events occurred at the junction of the Orozco Fracture Zone and the Middle America Trench, we speculate that the observed complex fault plane is caused by subduction of the rugged ocean floor of the Orozco Fracture Zone. Short-term precursory seismicity prior to the Petatlan earthquake can be explained by using the asperity model and migration of a slip front from the south-east to the north-west across the main shock source region.  相似文献