اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدNon-equilibrium transport of inhomogeneous shale gas under ultra-tight confinement
106
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
The non-equilibrium transport of inhomogeneous and dense gases highly confined by surface is encountered in many engineering applications. For example, in the shale gas production process, methane is extracted from ultra-tight pores under high pressure so the gas is inhomogeneous and dense. Currently, the complex non-equilibrium transport of inhomogeneous and dense gases where gas surface interactions play a key role is commonly investigated by molecular dynamics or on a continuum-assumption basis. Here, a tractable kinetic model based on the generalized Enskog equation and the mean-field theory is employed to couple the effects of the volume exclusion and the long-range intermolecular attraction forces. The interactions between gas molecules and confined surface are modelled by a 10-4-3 Lennard-Jones potential, which can capture gas surface adsorption. The cross-sectional density profiles of methane under different confinements are in good agreement with the molecular dynamics results reported in the literature, and the transport behaviors are validated by the non-equilibrium molecular dynamics. The velocity of methane flow in shale matrix is plug-like due to its dense characteristics in nanopores. The influence of pressure, temperature, pore size and shale composition on density and velocity profiles is analyzed quantitatively. Our results show that the Klinkenberg correction is not applicable to model shale gas flow in the production process; the Navier-Stokes model using the second-order slip boundary condition cannot produce the proper velocity profiles, and consequently fails to predict the accurate flow rate in nanopores. This study sheds new light on understanding the physics of non-equilibrium dense gas flows in shale strata.
قيم البحث
اقرأ أيضاً
While the hydraulic fracturing technology, aka fracking (or fraccing, frac), has become highly developed and astonishingly successful, a consistent formulation of the associated fracture mechanics that would not conflict with some observations is sti
ll unavailable. It is attempted here. Classical fracture mechanics, as well as the current commercial softwares, predict vertical cracks to propagate without branching from the perforations of the horizontal well casing, which are typically spaced at 10 m or more. However, to explain the gas production rate at the wellhead, the crack spacing would have to be only about 0.1 m, which would increase the overall gas permeability of shale mass about 10,000$times$. This permeability increase has generally been attributed to a preexisting system of orthogonal natural cracks, whose spacing is about 0.1 m. But their average age is about 100 million years, and a recent analysis indicated that these cracks must have been completely closed by secondary creep of shale in less than a million years. Here it is considered that the tectonic events that produced the natural cracks in shale must have also created weak layers with nano- or micro-cracking damage. It is numerically demonstrated that a greatly enhanced permeability along the weak layers, with a greatly increased transverse Biot coefficient, must cause the fracking to engender lateral branching and the opening of hydraulic cracks along the weak layers, even if these cracks are initially almost closed. A finite element crack band model, based on recently developed anisotropic spherocylindrical microplane constitutive law, demonstrates these findings.
Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic---specifically, transversely isotropic---and closely linked to the layered composite
structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semi-analytical homogenization of a bi-layer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic Modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.
We report on a novel and flexible experiment to investigate the non-equilibrium melting behaviour of model crystals made from charged colloidal spheres. In a slit geometry polycrystalline material formed in a low salt region is driven by hydrostatic
pressure up an evolving gradient in salt concentration and melts at large salt concentration. Depending on particle and initial salt concentration, driving velocity and the local salt concentration complex morphologic evolution is observed. Crystal-melt interface positions and the melting velocity are obtained quantitatively from time resolved Bragg- and polarization microscopic measurements. A simple theoretical model predicts the interface to first advance, then for balanced drift and melting velocities to become stationary at a salt concentration larger than the equilibrium melting concentration. It also describes the relaxation of the interface to its equilibrium position in a stationary gradient after stopping the drive in different manners. We further discuss the influence of the gradient strength on the resulting interface morphology and a shear induced morphologic transition from polycrystalline to oriented single crystalline material before melting.
The effects of volume exclusion and long-range intermolecular attraction are investigated by the simplified kinetic model for surface-confined inhomogeneous fluids. Gas dynamics of the ideal gas, the hard-sphere fluid and the real gas are simulated b
y the Boltzmann equation, the Enskog equation and the simple kinetic equation, respectively. Only the Knudsen minimum appears for the ideal gas, while both the Knudsen minimum and the Knudsen maximum occur for the hard-sphere fluid and the real gas under certain confinements, beyond which the maximum and minimum may disappear. The Boltzmann equation and the Enskog equation overestimates and underestimates the mass flow rate of the real gas dynamics under confinement, respectively, where the volume exclusion and the long-range intermolecular attractive potential among molecules are not ignorable. With the increase of the channel width, gas dynamics of the hard-sphere fluid and the real gas tends to the Boltzmann prediction gradually. The density inhomogeneity, which hinders the flow under confinement, is more obvious when the solid fraction is larger. The anomalous slip occurs for real gas under constant confinement. The flow at a smaller Knudsen number (larger solid fraction or channel width) contributes more practical amount of mass transfer, although the rarefaction effects is more prominent at larger Knudsen numbers. The temperature has no effect on density and velocity profiles of the ideal gas and the hard-sphere fluid, but the energy parameter among the real gas molecules decreases with the increasing temperature and the real gas dynamics tends to the hard-sphere ones consequently.
Rate- and state-dependent friction law for velocity-step and healing are analysed from a thermodynamic point of view. Assuming a logarithmic deviation from steady-state a unification of the classical Dieterich and Ruina models of rock friction is proposed.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
