A physically based numerical approach is presented for modeling fracture–matrix interaction, which is a key issue for fractured reservoir simulation. Commonly used mathematical models for dealing with such interactions employ a dualor multiple-continuum concept, in which fractures and matrix are represented as overlapping, different, but interconnected continua, described by parallel sets of co...