A comprehensive model combining Laplace-transform finite-difference and boundary-element method for the flow behavior of a two-zone system with discrete fracture network

•New models are developed by use of the solutions in Laplace domain.•The system exhibits six flow periods.•Discrete fracture network conductivity dominates the early-time flow behavior.•Local solutions for inner zone linear flow are similar with that of a planar fracture.