Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation          下载免费PDF全文

Tongbi Tu  Ali Ercan  M. Levent Kavvas 《水文研究》2018,32(10):1406-1419

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.  相似文献   

Spatial measurements of stream baseflow,a relevant method for aquifer characterization and permeability evaluation. Application to a hard‐rock aquifer,the Oman ophiolite     

B. Dewandel  P. Lachassagne  A. Qatan 《水文研究》2004,18(17):3391-3400

The structure, functioning and hydrodynamic properties of aquifers can be determined from an analysis of the spatial variability of baseflow in the streams with which they are associated. Such analyses are based on simple low‐cost measurements. Through interpreting the hydrological profiles (Q = f(A)) it is possible to locate the aquifer(s) linked to the stream network and to determine the type of interrelated flow, i.e. whether the stream drains or feeds the aquifer. Using an analytical solution developed for situations with a positive linear relationship, i.e. where the baseflow increases linearly with increasing catchment size, it is also possible to estimate the permeability of the aquifer(s) concerned at catchment scale. Applied to the hard‐rock aquifers of the Oman ophiolite, this method shows that the ‘gabbro’ aquifer is more permeable than the ‘peridotite’ aquifer. As a consequence the streams drain the peridotites and ‘leak’ into the gabbro. The hydrological profiles within the peridotite are linear and positive, and indicate homogeneity in the hydrodynamic properties of these formations at the kilometre scale. The permeability of the peridotite is estimated at 5 · 10^?7 to 5 · 10^?8 m/s. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

定水头注水引起的含水层水平运动和应变   总被引：2，自引：2，他引：0  

王庆良  冉兴龙 《大地测量与地球动力学》2001,21(4):26-32

基于含水层固体颗粒与孔隙水不可压缩的假设 ,本文导出了单井注水情况下泰斯承压含水层水平运动速度与水头之间的基本关系式。然后利用注水井壁处的应力、应变边界条件 ,进一步导出了单井定水头注水引起的泰斯承压含水层水平运动速度、位移和应变解析表达式。该水平位移与应变由两部分组成 :一部分为由注水压力本身引起的经典弹性力学解项 ,它仅随半径而变化 ,与注水时间无关 ;另一部分为由地下水头变化引起的水动力学位移和应变解项。其中 ,含水层水动力学水平位移随时间加长呈指数增长特征 ,水动力学径向应变则表现为近井处拉张、远井处挤压的分区特征 ,且近井拉张区随时间加长逐渐向外扩展。单井注水含水层水动力学水平位移、应变解的导出 ,完善和发展了单孔内压经典弹性平面力学问题解  相似文献   

Slumping and shallow faulting related to the presence of salt on the Continental Slope and Rise off North Carolina     

K.V. Cashman  P. Popenoe 《Marine and Petroleum Geology》1985,2(3):260-271

Seismic reflection profiles and long- and medium-range sidescan sonar were used to investigate a salt diapir complex and area of slope instability near the base of the Continental Slope off North Carolina. Within the area of investigation three diapirs are bounded on their upslope side by a scarp 60 m high and 50 km long. The slope above the scarp is characterized by a series of shallow rotational normal faults. The bottom below the scarp is furrowed by slide tracks, which were probably carved by large blocks that broke off the scarp face and slid downslope leaving rubble and scree lobes.Extensive slumping in this area appears to be a result of uplift and faulting associated with salt intrusion, which has fractured and oversteepened the slope leading to instability and failure. Sharply defined slide tracks suggest that slope failure above the breached diapir complex is a continuing process, in contrast to much of the surrounding slope area where few instability features were observed.  相似文献   

Implications of Kali-Hindon inter-stream aquifer water balance for groundwater management in western Uttar Pradesh     

Rashid Umar  M. Muqtada A. Khan  Izrar Ahmed  Shakeel Ahmed 《Journal of Earth System Science》2008,117(1):69-78

The Kali-Hindon inter-stream region extends over an area of 395 km² within the Ganga-Yamuna interfluve. It is a fertile tract for sugarcane cultivation. Groundwater is a primary resource for irrigation and industrial purposes. In recent years, over-exploitation has resulted in an adverse impact on the groundwater regime. In this study, an attempt has been made to calculate a water balance for the Kali-Hindon inter-stream region. Various inflows and outflows to and from the aquifer have been calculated. The recharge due to rainfall and other recharge parameters such as horizontal inflow, irrigation return flow and canal seepage were also evaluated. Groundwater withdrawals, evaporation from the water table, discharge from the aquifer to rivers and horizontal subsurface outflows were also estimated. The results show that total recharge into the system is 148.72 million cubic metres (Mcum), whereas the total discharge is 161.06 Mcum, leaving a deficit balance of −12.34 Mcum. Similarly, the groundwater balance was evaluated for the successive four years. The result shows that the groundwater balance is highly sensitive to variation in rainfall followed by draft through pumpage. The depths to water level are shallow in the canal-irrigated northern part of the basin and deeper in the southern part. The pre-monsoon and post-monsoon water levels range from 4.6 to 17.7 m below ground level (bgl) and from 3.5 to 16.5 m bgl respectively. It is concluded that the groundwater may be pumped in the canal-irrigated northern part, while withdrawals may be restricted to the southern portion of the basin, where intense abstraction has led to rapidly falling water table levels.  相似文献   

A comparative study of four vulnerability mapping methods in a detritic aquifer under mediterranean climatic conditions     

M. Draoui  J. Vias  B. Andreo  K. Targuisti  J. Stitou El Messari 《Environmental Geology》2008,54(3):455-463

This paper presents the results of a comparative study relating to the application of four vulnerability mapping methods, GOD, AVI, DRASTIC and SINTACS, in a pilot detritic aquifer situated in NW Morocco, known as the Martil–Alila aquifer. The principal objective of this work is to determine the most suitable such methods for this aquifer type within a Mediterranean context, and to show the effect of the rainfall variations that are characteristic of the Mediterranean climate on the degree of vulnerability. The methods applied distinguish five classes of vulnerability, these being irregularly divided up in space, with the division varying according to the method in question. The vulnerability maps obtained by the different methods strongly suggest that the eastern half of the aquifer is more vulnerable to contamination than the western half, for all hydrological situations. The effect of climatic conditions on the degree of vulnerability is well represented by the DRASTIC, according to which the aquifer is moderately to strongly vulnerable during humid hydrological years and weakly to moderately vulnerable during dry ones. For the other methods, this climatic effect is limited to the area occupied by the two predominant classes (“High” and “Low” for GOD and “High” and “Moderate” for SINTACS) while it is null for AVI. In conclusion, DRASTIC appears the most suitable for mapping the vulnerability to contamination of Mediterranean coastal detritic aquifers such as the Martil–Alila aquifer.  相似文献   

Phosphorus transport in shallow groundwater in peri-urban Kampala, Uganda: results from field and laboratory measurements   总被引：1，自引：0，他引：1  

N. R. Kulabako  M. Nalubega  R. Thunvik 《Environmental Geology》2008,53(7):1535-1551

To understand Phosphorus (P) sources and transport processes in the subsurface in Bwaise III Parish, Kampala, P attenuation and adsorption capacities of soils were studied in situ and from laboratory measurements. Relationships between sorption parameters and soil matrix properties, rates and mechanism of the adsorption process and soil P fractions were also investigated. P was generally higher in the wet than the dry season, but for both seasons, the maximum was 5 mgP/l. P transport mechanisms appeared to be a combination of adsorption, precipitation, leaching from the soil media and by colloids with the latter two playing an important role in the wet season. The sorption process comprised two phases with the first stage rate constants being about fourfold those of the second stage. The Langmuir isotherm described the sorption data well (R ² ≥ 0.95) with the second soil layer exhibiting the highest sorption maximum (C _max) (average value 0.6 ± 0.17 mgP/gDW). The best prediction of C _max had organic carbon, Ca, available P and soil pH. Residual P consisting mostly of organics was the main fraction in all the layers followed by inorganic HCl-P and NaOH-P in the top and middle layers, respectively. Loosely bound P (NH₄Cl-P) was the least fraction (<0.4% of total P) in all layers indicating the high binding capacity of P by the soils. The study results suggest that P dynamics is related to Ca, Fe and organic carbon content of the soils.  相似文献   

Deciphering potential groundwater zone in hard rock through the application of GIS   总被引：6，自引：1，他引：5  

R. K. Prasad  N. C. Mondal  Pallavi Banerjee  M. V. Nandakumar  V. S. Singh 《Environmental Geology》2008,55(3):467-475

Remote Sensing and Geographic Information System has become one of the leading tools in the field of hydrogeological science, which helps in assessing, monitoring and conserving groundwater resources. It allows manipulation and analysis of individual layer of spatial data. It is used for analysing and modelling the interrelationship between the layers. This paper mainly deals with the integrated approach of Remote Sensing and geographical information system (GIS) to delineate groundwater potential zones in hard rock terrain. The remotely sensed data at the scale of 1:50,000 and topographical information from available maps, have been used for the preparation of ground water prospective map by integrating geology, geomorphology, slope, drainage-density and lineaments map of the study area. Further, the data on yield of aquifer, as observed from existing bore wells in the area, has been used to validate the groundwater potential map. The final result depicts the favourable prospective zones in the study area and can be helpful in better planning and management of groundwater resources especially in hard rock terrains.  相似文献   

Numerical modeling of rainstorm-induced shallow landslides in saturated and unsaturated soils     

Tung-Lin Tsai  Hung-En Chen  Jinn-Chuang Yang 《Environmental Geology》2008,55(6):1269-1277

For the assessment of shallow landslides triggered by rainfall, the physically based model coupling the infinite slope stability analysis with the hydrological modeling in nearly saturated soil has commonly been used due to its simplicity. However, in that model the rainfall infiltration in unsaturated soil could not be reliably simulated because a linear diffusion-type Richards’ equation rather than the complete Richards’ equation was used. In addition, the effect of matric suction on the shear strength of soil was not actually considered. Therefore, except the shallow landslide in saturated soil due to groundwater table rise, the shallow landslide induced by the loss in unsaturated shear strength due to the dissipation of matric suction could not be reliably assessed. In this study, a physically based model capable of assessing shallow landslides in variably saturated soils is developed by adopting the complete Richards’ equation with the effect of slope angle in the rainfall infiltration modeling and using the extended Mohr–Coulomb failure criterion to describe the unsaturated shear strength in the soil failure modeling. The influence of rainfall intensity and duration on shallow landslide is investigated using the developed model. The result shows that the rainfall intensity and duration seem to have similar influence on shallow landslides respectively triggered by the increase of positive pore water pressure in saturated soil and induced by the dissipation of matric suction in unsaturated soil. The rainfall duration threshold decreases with the increase in rainfall intensity, but remains constant for large rainfall intensity.  相似文献   

Modeling of water–rock interactions in an aquitard of sandy clay in the coastal area near Beihai,China     

Zhou?XunEmail author  Yao?Jinmei  Zhang?Hua  Li?Rui  Xu?Fang 《Environmental Geology》2008,56(1):183-188

Unconsolidated sand, gravel and clay deposits near Beihai and in the Leizhou Peninsula in southern China form an unconfined aquifer, aquitard and a confined aquifer. Water and soil samples were collected from the two aquifers in the coastal Beihai area for the determination of chemical compositions, minerals and soluble ions. Hydrogeochemical modeling of three flow paths through the aquitard are carried out using PHREEQC to determine water–rock interactions along the flow paths. The results indicate that the dissolution of anorthite, fluorite, halite, rhodochrosite and CO₂, and precipitation of potash feldspar and kaolinite may be occurring when groundwater leaks through the aquitard from the unconfined aquifer to the confined aquifer. Cation exchanges between Na and Ca can also happen along the flow paths.  相似文献