Gas Flow Behavior through Inorganic Nanopores in Shale Considering Confinement Effect and Moisture Content

Constructing a Microdiffusion-Seepage-Stress Multifield Coupling Model for Nanopore Gas

Existing research is difficult to fully capture the correlation between gas molecules and pore wall interactions, multiphase flow, and stress distribution in nanopores. Taking gas as an example, a microscopic model was constructed. At the same time, diffusion, seepage, and stress were considered to accurately predict and manage gas transport in nanopores. First, molecular dynamics (MD) simulation methods were adopted to simulate the motion trajectories and interactions of gas molecules in nanopores. Second, a multiscale model was established based on continuum mechanics to consider the interaction between pore walls and gas molecules, and a diffusion equation was established to describe the diffusion process of gas molecules in pores. Then, finite element analysis and porous media models were used to simulate the seepage behavior of gas in the nanopores. Finally, the stress distribution in the pores was analyzed, and the influence of the interaction between the pore wall and gas molecules on stress was considered. The multifield coupling model was experimentally evaluated from three aspects: diffusion coefficient, seepage behavior, and stress distribution. The root-mean-square error (RMSE) and mean absolute error (MAE) of the model in different testing directions were calculated using different simulation tools, such as COMSOL, ANSYS, OpenFOAM, and CFX. The mean values of RMSE and MAE were lower than 0.20 and 0.17, respectively. The constructed model can comprehensively describe gas transmission within nanopores, improving the management accuracy and efficiency.

Deprotonation of formic, acetic acids and bicarbonate ion in slit silica nanopores at infinite dilution and in the presence of electrolytes

Dielectric effects and the coupled electrostatics between the nanoconfined and the internal/external aqueous media contribute to the observed deviations of chemistry within the nanoconfined environment when compared with unconfined systems. A systematic understanding has remained elusive, especially with respect to background salt concentration and boundary condition effects like the nanopore surface chemistry and the reference state used to calculate free energies. We utilize molecular dynamics simulations along with thermodynamic integration to determine the free energy difference associated with acid-base chemistry in 2 nm and 4 nm slit pores open to a bulk-like reservoir. pKa increases are predicted when confining acetic acid, formic acid, and bicarbonate in the slits at infinite dilution conditions. We find that confinement weakens the acids, and the modulation of outer pore surface dipole magnitudes can tune the pKa shift values, suggesting that purely "intrinsic" electrostatic effect on confinement may not exist. At sufficiently high salt concentrations, the dielectric/electrostatic effects on pKa values diminish due to charge screening effects. These discoveries enable future modifications of nanopore chemistries to achieve desirable properties for industrial applications.

Phase Behavior of Hydrocarbon Fluids in Shale Reservoirs, Considering Pore Geometries, Adsorption, and Water Film

Phase behavior of hydrocarbon fluids in nanopores is different from that observed in a PVT cell due to the confinement effect. While scholars have established various models for studying the phase behavior in nanopores, the authors often ignore the effect of pore geometries, which can significantly affect the critical fluid properties in shale nanopores. In this study, we extend the Soave-Redlich-Kwong equation of state (SRK EOS) using potential theory and establish models of critical property shift, considering pore geometries, adsorption, and water film. Our research shows that the critical property shifts, considering fluid adsorption, begin at rp ≤ 10 nm and are seriously strengthened with nanopore radius reduction. The extended SRK EOS is applied to compute phase diagrams of the 50% C1-50% nC10 mixture at different pore sizes and find that the thickness of adsorption and water film causes a depression in the P-T diagram and that the bubble point pressure is lower in cylindrical pores. At pressures above 6 MPa, the irreducible water saturation and pore geometries greatly impact the vapor-liquid ratio. This study is significant for evaluating residual oil distribution and studying fluid flow laws in shale reservoirs.

An Apparent Permeability Model in Organic Shales: Coupling Multiple Flow Mechanisms and Factors

It is of great significance to study shale apparent permeability under the action of multiple flow mechanisms and factors because shale reservoirs possess complex pore structures and flow mechanisms. In this study, the confinement effect was considered, with the thermodynamic properties of gas being modified, and the law relating to the conservation of energy adopted to characterize bulk gas transport velocity. On this basis, the dynamic change of pore size was assessed, from which shale apparent permeability model was derived. The new model was verified by three steps: experimental and molecular simulation results of rarefied gas transport, shale laboratory data, and comparison with different models. The results revealed that, under the conditions of low pressure and small pore size, the microscale effects became obvious, which significantly improved gas permeability. Through comparisons, the effects of surface diffusion and matrix shrinkage, including the real gas effect, were obvious in the smaller pore sizes; nevertheless, the stress sensitivity effect was stronger in larger pore sizes. In addition, shale apparent permeability and pore size decreased with an increase in permeability material constant and increased with increasing porosity material constant, including internal swelling coefficient. The permeability material constant had the greatest effect on gas transport behavior in nanopores, followed by the porosity material constant; however, the internal swelling coefficient had the least effect. The results of this paper will be important for the prediction and numerical simulation of apparent permeability relating to shale reservoirs.

Modification of the Calculation Method for Dynamic Reserves in Tight Sandstone Gas Reservoirs

The determination of dynamic reserves is important for tight sandstone gas reservoirs in production. Based on the geological and gas data of the Yan'an gas field, the influence of pressure on the properties of natural gas is studied by mathematical methods. At the same time, the modified flowing material balance equation is established considering the changes in gas viscosity and compressibility. The result shows that (1) the viscosity of natural gas increases rapidly with pressure; (2) the deviation factor decreases with pressure (P < 15 MPa) and then increases (P > 15 MPa) with temperature; (3) the compressibility decreases rapidly with pressure and increases with temperature; (4) compared with the results of the material balance method, the average error of the flowing material balance method is 33.95%, and the accuracy of the modified flowing material balance method is higher with an average error of 1.25%; and (5) a large change in the production will affect the accuracy of the modified flowing material balance method, especially a shut-in for a long time before the pressure drop production is calculated at a certain time, so data points that are relatively consistent should be selected as far as possible to calculate the dynamic reserves. The findings of this study can help in the accurate evaluation of dynamic reserves of the tight gas reservoir in the Yan'an gas field and are an important guide for the formulation of a rational plan for the gas reservoir and its economic and efficient development.

Molecular Investigation of CO2/CH4 Competitive Adsorption and Confinement in Realistic Shale Kerogen

The adsorption behavior and the mechanism of a CO2/CH4 mixture in shale organic matter play significant roles to predict the carbon dioxide sequestration with enhanced gas recovery (CS-EGR) in shale reservoirs. In the present work, the adsorption performance and the mechanism of a CO2/CH4 binary mixture in realistic shale kerogen were explored by employing grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) simulations. Specifically, the effects of shale organic type and maturity, temperature, pressure, and moisture content on pure CH4 and the competitive adsorption performance of a CO2/CH4 mixture were investigated. It was found that pressure and temperature have a significant influence on both the adsorption capacity and the selectivity of CO2/CH4. The simulated results also show that the adsorption capacities of CO2/CH4 increase with the maturity level of kerogen. Type II-D kerogen exhibits an obvious superiority in the adsorption capacity of CH4 and CO2 compared with other type II kerogen. In addition, the adsorption capacities of CO2 and CH4 are significantly suppressed in moist kerogen due to the strong adsorption strength of H2O molecules on the kerogen surface. Furthermore, to characterize realistic kerogen pore structure, a slit-like kerogen nanopore was constructed. It was observed that the kerogen nanopore plays an important role in determining the potential of CO2 subsurface sequestration in shale reservoirs. With the increase in nanopore size, a transition of the dominated gas adsorption mechanism from micropore filling to monolayer adsorption on the surface due to confinement effects was found. The results obtained in this study could be helpful to estimate original gas-in-place and evaluate carbon dioxide sequestration capacity in a shale matrix.

Magnetohydrodynamic nanofluid radiative thermal behavior by means of Darcy law inside a porous media

Radiative nanomaterial thermal behavior within a permeable closed zone with elliptic hot source is simulated. Darcy law is selected for simulating permeable media in existence of magnetic forces. Contour plots for various buoyancy, Hartmann numbers and radiation parameter were illustrated. Carrier fluid is Al2O3-water with different shapes. Outputs prove that conduction mode augments with enhance of Ha. Nu augments with considering radiation source term.