Tight shales have two pore networks. The matrix — nanometer-scale pores in the rock itself — stores fluid but conducts it slowly. The fracture network — cracks, fissures, induced or natural — conducts fluid rapidly but stores little. Between them, mass transfers at a rate that depends on pressure difference, contact area, and permeability contrast. The two networks are coupled but operate at different scales.
Classical finite element methods handle single-porosity flow well. Dual-porosity flow introduces discontinuities at network interfaces that generate spurious oscillations on any fixed mesh. The mesh can't resolve both the millimeter-scale fractures and the centimeter-scale matrix blocks without becoming computationally prohibitive.
An adaptive physics-informed neural network framework (arXiv:2603.19561) encodes the governing equations — conservation, Darcy flow, mass exchange — directly into the loss function. The key innovations: adaptive weight tuning balances the competing physics terms, dynamic collocation point selection concentrates computational effort where the solution changes fastest, and a shared neural architecture captures the coupled behavior across both networks.
The structural insight: the framework is inherently mesh-free. It doesn't discretize space — it learns a continuous function that satisfies the physics everywhere. This eliminates the mesh-dependent oscillations by construction. The same architecture handles inverse problems: inferring the mass transfer coefficient (which can't be measured directly) from observable pressure data. The physics doesn't just constrain the forward model. It enables the inverse one.