Analysis of karst spring recession curves,west of Iran

The recession hydrographs of karst springs provide important information about aquifer characteristics such as storage properties and drainage potential, karstification degree, and other hydrological features. The Faresban, Famaasb, Gonbad-e-Kabood, and Gian springs are the important karst springs which drain main karst aquifers in west of Iran. In this study, we evaluate hydrological characteristics of these four karst springs using different analytic equations. Almost, all the recession curves were well fitted by a function that consists of three main exponential terms. The initial portion of the recession curve represents the fast drainage of large fractures and conduits, ending with the slowly decreasing curve, where the drainage of rock matrix and small fissures or baseflow is dominant. The Gian spring has the largest storage volume but the lowest drained water (7%). Using Mangin’s equation, it is apparent that the baseflow supplies over 85% of the total drained water to the Gian spring, while a small quantity of total outflow is made up of quickflow periods. Of the water available to the Famaasb, Faresban, and Gonbad-e-Kabood springs, 85, 80, and 70% are drained during period of baseflow, respectively, whereas the quickflow periods have minor importance in some springs. Based on the Coutagne’s equation results, it can be concluded that the catchment of the Famaasb and Gonbad-e-Kabood springs have a karstified saturation zone with drainage channels that provide a rapid discharge with little storage capacity. At Gian catchment area, the initial discharge is not as rapid as karstification but is fairly homogeneous throughout the aquifer.     

Numerical simulation of crack propagation in layered formations

In the present study, a boundary element method based on the higher order displacement discontinuity formulation is presented to solve the general problem of hydraulic fracture propagation in layered formations. Displacement collocation technique is employed to model the higher order displacement variation along the crack and the special crack tip element near its ends. The hydraulic fracture propagation and its interaction with the layer interface in non-homogenous rock materials are studied by the proposed semi-analytical (hybridized boundary element-boundary collocation) method. The maximum tangential stress criterion (or σ-criterion) of fracture mechanics considering different elastic constants (Young modulus and Poisson’s ratio) is used to obtain the fracture path. The fracture propagation from stiff to soft and soft to stiff media for cracks having different inclination angles is modeled, and the effects of elastic constants on the hydraulic fracture propagation is studied. The results show that if the hydraulic fracture originates in the stiffer layer, its capability to cross the layer increases and is vice versa for the softer material. The comparison of the results gained from the numerical method with those in the literature show a good performance of the method in the case of propagation of hydraulic fracture in layered formations.