Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shale”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Science-informed Machine Learning to Increase Recovery Efficiency in Unconventional Reservoirs

Production of hydrocarbons from fractured, unconventional reservoirs is inherently inefficient. But machine learning offers a pathway both to increasing recovery efficiency at a site and to improving forecasts of production, thereby improving the economics of operations in unconventional reservoirs. Los Alamos—in partnership with DOE, NETL, and WVU—has been developing a science-informed workflow and platform for optimizing pressure-drawdown at a site, which will allow an operator to make reservoir-management decisions that optimize recovery in consideration of future production. This work relies on a hybridization of physics-based prediction and machine learning, whereby accurate synthetic data (in combination with available site data) can enable the application of machine learning methods for rapid forecasting and optimization. The physics-based prediction is built upon experimental and theoretical work to determine transport characteristics in shale at various scales, with an emphasis on materials from MSEEL-I; this fundamental shale R&D was conducted in partnership with DOE, NETL, and several other national labs. This work has resulted from a coordinated leveraging of developments across several projects within DOE FE30, along with internal investments from Los Alamos via LDRD. The development has utilized data from the MSEEL–I site for calibration and demonstration; however, the workflow and platform are readily extendable to operations at other sites, plays, and basins. This machine-learning method can aid operators to improve both recovery efficiency and competitiveness; to this end, future work would quantify processes for other plays/basins and integrate production details with economics.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2020

CO₂ Storage prospeCtive Resource Estimation Excel aNalysis (CO₂-SCREEN) User’s Manual

This user’s manual guides the use of the National Energy Technology Laboratory’s (NETL) CO₂ Storage prospeCtive Resource Estimation Excel aNalysis (CO₂-SCREEN) tool, which was developed to aid users screening geologic formations for prospective CO₂ storage resources. This manual is specific to the CO₂-SCREEN 4.0 version which is based in Python. CO₂-SCREEN applies U.S. Department of Energy (DOE) methods and equations for estimating prospective CO₂ storage resources for saline formations, shale formations, and residual oil zones (ROZ). CO₂-SCREEN was developed to be substantive and user-friendly and provide a consistent method for calculating prospective CO₂ storage resources. CO₂-SCREEN uses a Java based graphical user interface for data inputs and uses Python to calculate prospective CO₂ storage resources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Evaluating Production Implications of Pressure Maintenance in Unconventional Oil and Gas Wells using a Machine Learning Modeling Approach: Case Study in the Permian Basin

This study implements the proprietary machine learning-based model (model) developed under the 2022 study titled “Evaluating the Impact of Proprietary Oil & Gas Data on Machine Learning Model Performance Using a Quasi-Experimental Analytical Approach” for forecasting unconventional oil and gas production using well data from the Permian Basin. The model, developed using an exclusive dataset that includes time series production data from an operator in the Permian Basin, is designed to jointly predict daily oil, gas, and water production for horizontal wells as a function of bottomhole pressure drawdown, spatial placement across the study domain, and well-completion attributes. In this study, the model was explicitly applied to explore its utility to evaluate the impact of varying drawdown strategies on the production forecast of a well in the Permian Basin dataset. Managing pressure drawdown has been identified as a way to improve estimated ultimate recovery (EUR) from unconventional shale wells due to the stress-dependent nature of fractures in shale reservoirs. Research has shown that applying a lower pressure drawdown helps to maintain the reservoir conductivity, resulting in higher productivity over the life of a well. Historic bottomhole pressure data from the well over time was used as a benchmark from which to set more and less aggressive pressure decline rates as bounding modeling cases. All pressure decline rates/strategies were forecasted over 5 years, and the model was used to generate oil, water, and gas prediction over the same timeframe. This report presents the results from the production forecast and discusses potential operational and economic implications, with contrasting perspectives between well productivity and profitability given typical oil and gas economics and the volatility in the oil and gas market.

02 PETROLEUM↗
DOE OSTI2025

Marcellus assessment shows continued expansion

Increasing uncertainty regarding the remaining volume of recoverable gas in Appalachian shale plays provided the original motivation for this study. Assessments of the Marcellus and Utica shale plays published by the US Geological Survey (USGS) in 2019 established resource volumes dramatically lower than ongoing well performance would indicate. Further, as production in the Marcellus increased and wells matured, operators observed that wells were producing more than the estimated gas-in-place (GIP).

03 NATURAL GAS↗
DOE OSTI2021

Deposition, Diagenesis, and Sequence Stratigraphy of the Pennsylvanian Morrowan and Atokan Intervals at Farnsworth Unit

Farnsworth Field Unit (FWU), a mature oilfield currently undergoing CO2-enhanced oil recovery (EOR) in the northeastern Texas panhandle, is the study area for an extensive project undertaken by the Southwest Regional Partnership on Carbon Sequestration (SWP). SWP is characterizing the field and monitoring and modeling injection and fluid flow processes with the intent of verifying storage of CO2 in a timeframe of 100–1000 years. Collection of a large set of data including logs, core, and 3D geophysical data has allowed us to build a detailed reservoir model that is well-grounded in observations from the field. This paper presents a geological description of the rocks comprising the reservoir that is a target for both oil production and CO2 storage, as well as the overlying units that make up the primary and secondary seals. Core descriptions and petrographic analyses were used to determine depositional setting, general lithofacies, and a diagenetic sequence for reservoir and caprock at FWU. The reservoir is in the Pennsylvanian-aged Morrow B sandstone, an incised valley fluvial deposit that is encased within marine shales. The Morrow B exhibits several lithofacies with distinct appearance as well as petrophysical characteristics. The lithofacies are typical of incised valley fluvial sequences and vary from a relatively coarse conglomerate base to an upper fine sandstone that grades into the overlying marine-dominated shales and mudstone/limestone cyclical sequences of the Thirteen Finger limestone. Observations ranging from field scale (seismic surveys, well logs) to microscopic (mercury porosimetry, petrographic microscopy, microprobe and isotope data) provide a rich set of data on which we have built our geological and reservoir models.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Pyrolysis of Oils from Unconventional Resources

In this study, oils from various sources were subjected to pyrolysis conditions; that is, without oxidizer, as the samples were heated to 500 °C, and held at that temperature. The oils studied included: (1) heavy oil from Grassy Creek, Missouri; (2) oil from tar sands of Asphalt Ridge in Utah; (3) mid-continent oil shales of three formations (two of Chattanooga formation, Pennsylvanian (age) formation, and Woodford formation); and (4) a Colorado Piceance Basin shale. Differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA) with either gas chromatography (GC) or mass spectrometry (MS) were used to quantify the produced gases evolved in the tests. Purge gases of helium, argon, and humid carbon dioxide were utilized. Larger scale pyrolysis tests were conducted in a tube furnace coupled to a MS and a GC. The results consistently showed that pyrolysis occurred between 300 °C and 500 °C, with the majority of gases being mainly hydrogen and light alkanes. This behavior was essentially consistent, regardless of the oil source.

42 ENGINEERING↗
DOE OSTI2023

Temporal patterns of fine‐root dynamics have little influence on seasonal soil CO 2 efflux in a mixed, mesic forest

Among the contributors to soil CO 2 efflux, there remains uncertainty about the contribution of root activity to the overall soil efflux. Soil water and temperature frequently have been used to predict a large portion of the variation in soil CO 2 efflux. We hypothesized that fine-root dynamics explain most of the remaining variability in soil CO 2 efflux that cannot be explained by soil temperature and water content. We anticipated that seasonal increases in root production, mortality via decomposition, and standing crop would result in corresponding increases in soil CO 2 efflux. We tested our hypotheses by collecting and analyzing two years of minirhizotron and soil chamber CO 2 flux data from plots distributed throughout the Shale Hills Catchment of the Susquehanna-Shale Hills Critical Zone Observatory in Central Pennsylvania, USA. Here we showed that: (1) seasonal fluctuations in fine-root dynamics yielded only a very small increase in the predictability of soil CO 2 efflux; (2) fine-root mortality effects on soil CO 2 efflux were strongly tied to soil temperature; (3) fluctuations in fine-root presence or standing mass independent of temperature and moisture had little effect on soil CO 2 efflux; and (4) new fine-root length and root length mortality had limited impacts on soil CO 2 efflux rates. We conclude that, at least in temperate forests on rocky soils, characterizing fine-root dynamics may provide only limited improvement in the estimation of soil CO 2 efflux.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Manganese uptake by red maples in response to mineral dissolution rates in soil

Manganese (Mn) is an essential plant micronutrient that influences photosynthesis, ecosystem productivity, and soil carbon storage. Our objective was to quantify how Mn uptake by forest vegetation relates to Mn release into soil solution through mineral dissolution. In this work, a greenhouse pot experiment was conducted to quantify Mn uptake by red maple saplings as a function of mineral solubility to test whether Mn uptake was limited by the supply of Mn to soil solution. Differences in soil microbial community composition between treatments, particularly amongst Mn cycling bacteria and fungi, were also evaluated to assess potential microbial impacts on observed Mn fluxes. Plant Mn uptake was highest in systems supplied with dissolved Mn(II) because it was not kinetically limited by mineral weathering. Mn uptake was also higher in systems supplied with a fast-weathering substrate (shale containing Mn(II)-bearing pyrite) than a slow-weathering substrate (Mn(IV)-oxide). However, vegetation accumulated a decreasing proportion of available Mn with increasing solubility, indicating that uptake was tempered relative to solubility. The presence of bacterial phyla containing known Mn-oxidizing bacteria indicates the potential for Mn-oxidizing bacteria to influence Mn solubility within Mn-oxide and dissolved Mn treatments. A relatively low abundance of these bacteria points to their limited capacity to oxidize and retain Mn in the shale treatments, consistent with substantial Mn leaching. We conclude that Mn uptake was primarily controlled by dissolution rates of Mn-bearing minerals but modified by biological processes. The quantitative framework presented here can guide understanding of how biogeochemical processes control element cycling between plants and soils.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Incorporating geological structure into sensitivity analysis of subsurface contaminant transport

Simulating subsurface contaminant transport at the kilometer-scale often entails modeling reactive flow and transport within and through complex geologic structures. These structures are typically meshed by hand and as a result geologic structure is usually represented by one or a few deterministically generated geological models for uncertainty studies of flow and transport in the subsurface. Uncertainty in geologic structure can have a significant impact on contaminant transport. In this study, the impact of geologic structure on contaminant tracer transport in a shale formation is investigated for a simplified generic deep geologic repository for permanent disposal of spent nuclear fuel. An open-source modeling framework is used to perform a sensitivity analysis study on transport of two tracers from a generic spent nuclear fuel repository with uncertain location of the interfaces between the stratum of the geologic structure. The automated workflow uses sampled realizations of the geological structural model in addition to uncertain flow parameters in a nested sensitivity analysis. Concentration of the tracers at observation points within, in line with, and downstream of the repository are used as the quantities of interest for determining model sensitivity to input parameters and geological realization. Finally, the results of the study indicate that the location of strata interfaces in the geological structure has a first-order impact on tracer transport in the example shale formation, and that this impact may be greater than that of the uncertain flow parameters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

A critical literature review of the advances in methane dehydroaromatization over multifunctional metal-promoted zeolite catalysts

We report the drastic rise in shale gas production has encouraged the quest for alternative uses of methane as a chemical feedstock in the manufacturing industry. While two-step syngas routes for methane valorization are deployed commercially, direct one-step routes for methane conversion are attracting much attention. As steam cracking installations have shifted from using oil-based naphtha to shale-based natural gas liquids, production of aromatics has dropped. Methane dehydroaromatization (MDA) is a one-step reaction capable of valorizing methane to hydrogen and benzene. Challenges with the MDA reaction are two-fold: the reaction is thermodynamically limited with low one-pass methane conversion and even the best catalytic systems, Mo/zeolites, suffer rapid deactivation from coking. A catalyst design strategy to improve stability is the use of multifunctional Mo-X/zeolite systems where X is a dopant capable of modulating the stability. In this paper we provide a complete overview of the main Mo-X/zeolite systems used in MDA and critically draw connections among the different types of dopants (X) employed, as a function of the role they play in the reaction/deactivation pathway. We have also dedicated a section to emerging trends with non-Mo based catalysts. The goal of this review article is to establish a basis that will facilitate the identification of useful multifunctional catalytic systems, and recognize gaps in the knowledge of these systems that deserve more attention. Improving MDA systems to the point to which they can be commercially deployed requires a multifaceted approach that combines optimization of the designs of both the catalyst and the reactor configuration. We therefore also provide a brief overview of the most recent advances in process intensification strategies employed with different reactor configurations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Determination of transition metal ions in fossil fuel associated wastewaters using chelation ion chromatography

Here, this study outlines the development and subsequent validation of a method using chelation ion chromatography (CIC) pretreatment followed by traditional ion chromatography (IC) and post column UV/vis detection to measure transition metals in fossil fuel wastewaters, such as oil &amp; gas (O&G) brines and coal mine drainage (CMD) waters. Measurement of transition metals is often an important characterization step in the research of environmental and energy systems. IC represents one way to measure these metals with the advantages of being versatile, simple and relatively low cost compared to other analytical methods. However, high concentrations of alkali and alkaline earth metals present in fossil fuel wastewaters will decrease IC detectability of transition metals in these waters. In this study, a CIC method was developed for the analysis of transition metal ions (Fe 3+ , Cu 2+ , Ni 2+ , Zn 2+ , Co 2+ , Mn 2+ , and Fe 2+ ) in fossil fuel associated wastewaters such as Appalachian CMD and O&G wastewaters from the Permian and Bakken shale basins in the United States. CIC system incorporated an on-line chelator column (e.g., the MetPac CC-1) with high selectivity for transition metals over alkali and alkaline earth metals for salt matrix removal prior to transition metal separation and detection. Additional method developments also included acidifying all samples to 2% v/v HCl and using gradient elution rather than isocratic. The recoverability of transition metals in simple salt solutions commonly found in CMD and brine samples (e.g. NaCl, Na 2 SO 4 , CaCl 2 ) using CIC was evaluated and compared to that using traditional IC. Our results found that the CIC system significantly improved transition metal recoveries for samples in 10,000 mg/L CaCl 2 matrix, reaching 87%-108% recovery for all analytes, as opposed to 2-323% recovery in traditional IC. The limits of detection in this study achieved 10.09 – 161.2 μg/L, comparable to reported values in similar IC studies. The developed method was also verified with certified water samples, resulting in 89% - 111% recoveries in samples with higher analyte concentrations (i.e. >4x the LoDs). The developed method achieved 87% -112% recoveries for most analytes in CMD samples and 72%-138% recoveries for Bakken shale samples, relative to ICP-MS values. Overall, the current IC method can be a very good screening tool for fast and cheap analysis for transition metals at mg/L level, to facilitate selection of samples for more detailed ICP-MS analysis.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

How to avoid multiple scattering in strongly scattering SANS and USANS samples

Small Angle Neutron Scattering (SANS) and Ultra Small Angle Neutron Scattering (USANS) are the only available experimental techniques to provide seamless non-destructive measurements of the geometry of the accessible and inaccessible pore structure of rocks from sub-nanopore size to the scale of macropores. They have therefore become the measurement of choice for tight reservoir rocks such as organic rich shales. A simplifying assumption in the analysis is, however, that during the path of neutrons through the sample each neutron is only scattered once. Shales are samples with a high scattering power and Multiple Scattering (MS) may occur which requires special modelling for deconvolution of the results. The approach to avoid MS is to simply reduce the sample thickness to <0.15–0.5 mm. Here, in this work, we present a systematic method on wavelength selection and preparation of samples to optimise extraction of microstructural data and minimise parasitic errors. Experimentally measured SAS transmission (TSAS) values are used as a practical criterion for estimation of the extent of MS. Generous beamtime allocations allowed robust testing revealing that sample thicknesses can be twice as thick as predicted using the standard protocol. Analysing thicker samples is particularly beneficial for statistically relevant characterisation of heterogeneous samples making the new protocol the method of choice for such samples.

(U)SANS↗
DOE OSTI2024

A dynamic solvent chamber propagation estimation framework using RNN for warm solvent injection in heterogeneous reservoirs

Warm solvent injection (WSI), injecting low-temperature solvent into formations to reduce the viscosity of heavy oil, is a clean technology for heavy oil production through reducing greenhouse gas emissions and water usage. The success of WSI operation depends on the uniform development and propagation of solvent chambers in reservoirs. However, reservoir heterogeneity stemming from shale barriers plays a detrimental role in the conformance of solvent chamber development and oil production rate. In this work, we developed a novel recurrent neural network (RNN)-based framework with the capability of efficiently tracking and estimating the solvent chamber positions in heterogeneous reservoirs based on only production time-series data. The developed estimation model utilizes the “sequence-to-sequence" mapping methodology to correlate observed production time-series sequence and solvent chamber edge sequence via a long short-term memory (LSTM) algorithm. The trained RNN models exhibit high accuracy, evidenced by the predicted dynamic solvent chamber locations match the corresponding true locations from numerical simulation, with a high coefficient of determination (R 2 ) and a low mean squared error. Specifically, the achieved R 2 values exceed 0.98 on both the training and testing data. The developed RNN-based workflow was tested via several cases from both regularly- and irregularly-shaped shale barriers, and the results were promising. The predicted solvent chambers showed strong agreement with those obtained from numerical simulations. The major benefits of this workflow include reducing computational time and saving overall monitoring and tracking costs for conventional techniques. In conclusion, the present work would provide a good demonstration of the capability of practical integration of machine learning methods in solving engineering problems.

58 GEOSCIENCES↗
DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗
DOE OSTI2020

Role of interaction between hydraulic and natural fractures on production

One of the main purposes of hydraulic fracturing in unconventional resources such as shale is to improve the connectivity between existing natural fractures and the production well. Since the inherent permeability of shale is extremely low, this newly formed fracture connectivity provides faster paths for the hydrocarbon stored in preexisting natural fractures. In this study, we performed numerical simulations of the free gas production at large scale fractured reservoir with different pre-existing natural fracture network intensities and with variable hydraulic fracture sizes and staging. We explored production characteristics as a function of the properties of both, natural fractures (intensity) and hydraulic fractures (size and staging). Our hypothesis is that once maximum fracture connectivity is achieved, additional stimulation does not increase production significantly. In this numerical study we use an advanced modeling tool to study fractures connectivity, Discrete Fracture Network (DFN) approach, which can represent fracture networks similar to those observed on individual reservoir sites. In our model the horizontal production well and vertical hydraulic fractures were deterministically defined and surrounded by stochastically generated natural fracture networks. The performed numerical experiments show how fracture connectivity depends on hydraulic fracture settings in fractured reservoir with different preexisting fracture densities. We numerically showed that cumulative free gas production is not affected significantly by hydraulic fractures size (25 m–150 m) or spacing (15 m–45 m) for dense natural fracture networks. However, the hydraulic fracture size is important in reservoirs with sparse natural fractures, although spacing (number of stages) does not play a significant role.

58 GEOSCIENCES↗
DOE OSTI2021

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗