## A double-porosity model for a fractured aquifer with non-Darcian flow in fractures

HYDROLOGICAL SCIENCES JOURNAL-JOURNAL DES SCIENCES HYDROLOGIQUES, vol.53, no.4, pp.868-882, 2008 (SCI-Expanded)

• Publication Type: Article / Article
• Volume: 53 Issue: 4
• Publication Date: 2008
• Doi Number: 10.1623/hysj.53.4.868
• Journal Name:
• Journal Indexes: Science Citation Index Expanded (SCI-EXPANDED), Scopus
• Page Numbers: pp.868-882
• Keywords: double porosity, fractured aquifer, Forcheimer's equation, non-Darcian flow, POROUS-MEDIA, GROUNDWATER-FLOW, RESERVOIRS, SIMULATION, PATTERN, ROCKS, WELL
• Middle East Technical University Affiliated: Yes

#### Abstract

Non-Darcian flow in a finite fractured confined aquifer is Studied. A stream bounds the aquifer at one side and all impervious Stratum at the other. The aquifer consists of fractures capable of transmitting water rapidly, and Porous blocks which mainly store water. Unsteady flow in the aquifer due to a Sudden rise in the stream level is analysed by the double-porosity conceptual model. Governing equations for the flow in fractures and blocks are developed using the continuity equation. The fluid velocity in fractures is often too high for the linear Darcian flow so that the governing equation for fracture flow is modified by Forcheimer's equation, which incorporates a nonlinear term. Governing equations are coupled by all interaction term that controls the quasi-steady-state fracture-block interflow. Governing equations are solved numerically by the Crank-Nicolson implicit scheme. The numerical results are compared to the analytical results for the Same problem which assumes Darcian flow in both fractures and blocks. Numerical and analytical Solutions give the same results When the Reynolds number is less than 0.1. The effect of nonlinearity oil the flow appears when the Reynolds number is greater than 0.1. The effect the rate of flow From the stream to the aquifer, the higher the degree of nonlinearity. The effect of aquifer parameters oil the flow is also investigated. The proposed model and its numerical solution provide a useful application of nonlinear flow models to fractured aquifers. It is possible to extend the model to different types of aquifer, as well as boundary conditions at the stream side. Time-dependent flow rates in the analysis of recession hydrographs could also be evaluated by this model.