摘要
A physically based numerical approach is presented for modeling fractureaaamatrix interaction, which is a key issue for fractured reservoir simulation. Commonly used mathematical models for dealing with such interactions employ a dual- or multiple-continuum concept, in which fractures and matrix are represented as overlapping, different, but interconnected continua, described by parallel sets of conservation equations. The conventional single-point upstream weighting scheme, in which the fracture relative permeability is used to represent the counterpart at the fractureaaamatrix interface, is the most common scheme by which to estimate flow mobility for fractureaaamatrix flow terms. However, such a scheme has a serious flaw, which may lead to unphysical solutions or significant numerical errors. To overcome the limitation of the conventional upstream weighting scheme, this paper presents a physically based modeling approach for estimating physically correct relative permeability in calculating multiphase flow between fractures and the matrix, using continuity of capillary pressure at the fractureaaamatrix interface. The proposed approach has been implemented into two multiphase reservoir simulators and verified using analytical solutions and laboratory experimental data. The new method is demonstrated to be accurate, numerically efficient, and easy to implement in dual- or multiple-continuum models.