Analytical and numerical studies on a moving boundary problem of non-Newtonian Bingham fluid flow in fractal porous media

Abstract

Non-Darcy flow with a threshold in fractal porous media has been widely used in the development of unconventional petroleum resources such as heavy oil and tight oil. Mathematical modeling of such challenging “threshold flow” problems with strong nonlinearity has great significance in improving petroleum science and technology. Based on a fractal theory, a new non-Darcy kinematic equation with a fractal threshold pressure gradient (TPG) is mathematically deduced in order to describe the non-Darcy flow of a non-Newtonian Bingham fluid with a threshold in fractal porous media. Then mathematical modeling of planar radial non-Darcy flow in a fractal heavy oil reservoir is performed as a nonlinear moving boundary problem. In addition, a steady analytical solution method and a transient numerical solution method are developed. The analytical solution of an ordinary differential equation system for a simple steady model is derived directly, and the transient numerical solution of a partial differential equation system for an unsteady flow model is obtained based on the finite element method with good convergence. These two model solutions are validated by cross-comparisons. It is found from the calculation results that for the steady state, the extremely disturbed moving boundary and its corresponding pressure distribution are affected only by a TPG, production pressure, and a transport exponent; by contrast, for the unsteady state, the moving boundary and its corresponding pressure distribution are affected by many more factors including the fractal dimension. Furthermore, neglect of the fractal TPG and the induced moving boundary can lead to high overestimation of well productivity.

A novel dynamic boundary model of pressure propagation for volume fractured horizontal wells in shale oil reservoirs with developed natural fractures

“…Define the strip source center as xw and the width as xf. It is known that the line source function [11] has the following form:…”

Progress of Seepage Law and Development Technologies for Shale Condensate Gas Reservoirs

“…The diagram of reservoir heterogeneity at the SRV area is shown in Figure 10. The expression of non-uniform porosity at the SRV area considering stress sensitivity is as follows [104,105]:…”

“…The expression of non-uniform permeability at the SRV area considering stress sensitivity is as follows [104,106]:…”

“…At present, a series of experimental [107,108] and theoretical [107,[109][110][111][112] studies on reservoir heterogeneity characteristics have been carried out. In 2021, using high-precision high-pressure mercury intrusion (HPMI) experimental technology, the micro/nanopore The expression of non-uniform porosity at the SRV area considering stress sensitivity is as follows [104,105]:…”

A Simplified Lattice Boltzmann Boundary Conditions for Gas Transport in Self-Affine Microchannels with an Inherent Roughness of in a Tight Reservoir

“…Since Mandelbrot's [14] introduction of the concept of the fractal, fractal theory has been applied in many fields. It is used in studies of oil and gas to describe the complexity and irregularity of the microscopic pore-throat structures within sandstones to avoid the complexity of quantitative characterization [15][16][17]. Johnson et al [18] introduced the term "fractal" to describe the roughness of the pore surface and to calculate a compensation for the viscous attenuation in porous media caused by this roughness.…”