Mathematical Modeling of Pressure-Dependent Variation in the Hydrodynamic Parameters of Gas Fields

This study introduces a mathematical framework for analyzing unsteady gas filtration in porous media with pressure-dependent porosity variations. The physical process is formulated as a nonlinear parabolic boundary value problem that captures the coupled interaction between pressure evolution and porosity changes during gas production. To solve the equation, a numerical strategy is developed by integrating the Alternating Direction Implicit (ADI) scheme with quasi-linearization iterations, employing finite difference discretization on a two-dimensional spatial grid. Extensive computational experiments are performed to investigate the influence of key reservoir parameters—including porosity coefficient, permeability, gas viscosity, and well production rate—on the spatiotemporal behavior of pressure and porosity during long-term extraction. The results indicate significant porosity variations near the wellbore driven by local pressure depletion, reflecting strong sensitivity of the system to formation properties. The validated numerical model provides valuable quantitative insights for optimizing reservoir management and improving production forecasting in gas field development. Overall, the proposed methodology serves as a practical tool for oil and gas engineers to assess long-term reservoir performance under diverse operational conditions and to design efficient extraction strategies that incorporate pressure-dependent formation property changes.