The Design of a Three-Dimensional Physical Modeling System for Real-Time Groundwater Flows

In the past decades,physical modeling has been widely used in hydrogeology for teaching,studying and exhibition purposes.Most of these models are used to illustrate hydrogeological profiles,but few can depict three-dimensional groundwater flows,making it impossible to validate groundwater flows simulated by numerical methods with physical modeling.

The research of three-dimensional numerical simulation of groundwater-flow: taking the Ejina Basin, Northwest China as example

Water is a primary controlling factor for economic development and ecological environmental protection in the inland river basins of arid western China. And it is groundwater, as the most important component of total water resources, that plays a dominant role in the development of western China. In recent years, the use-ratio of surface water has been raised, the groundwater recharge rate from surface water has been reduced, and groundwater has been exploited on a large scale. This has led to the decline of ground-water levels and the degradation of eco-environments in the Heihe watershed. Therefore, the study on the change in groundwater levels in recent years, as well as simulating and predicting groundwater levels in the future, have become very significant for im-proving the ecological environment of the Heihe River Basin, to coordinate the water contradiction among upper, middle and lower reaches of Heihe River Basin and to allocate the water resources. The purpose of this study is to analyze the groundwa-ter-level variations of the Ejina region based on a large scale, to develop and evaluate a conceptual groundwater model in Ejina Basin, to establish the groundwater flow model using the experimental observation data and combining Modular Three-Dimensional Groundwater Flow Model (MODFLOW) and GIS software, to simulate the regional hydrologic regime in re-cent 10 years and compare various water-delivery scenarios from midstream, and to determine which one would be the best plan for maintaining and recovering the groundwater levels and increasing the area of Ejina oasis. Finally this paper discusses the pos-sible vegetation changes of Ejina Basin in the future.

Hydrochemical Constraints on the Flowing Paths of Groundwater in Limestone Reservoirs beneath the Pingdingshan Coalfield in North China

The karst groundwater of Cambrian limestone may become an important water source for industry and agriculture in the Pingdingshan area,and is also a potential threat to mining safety.Therefore,to find out the origin,flow paths,and hydrogeochemical processes of karst groundwater beneath the Pingdingshan coalfield,a total of 16 water samples were collected.Our findings confirmed that the karst groundwater is mainly recharged by precipitation.The precipitation can directly supply the deep aquifer of the karst water system through the southwest limestone outcrops,and this area mostly includes the southern part of mines No.11,No.9,and the hidden outcrops in the southern part of mine No.2.What is more,the areas adjacent to the synclinal axis,including mines No.10,No.12,and No.8,may be the main discharge areas.A mixing model of^(87)Sr/^(86)Sr and Sr showed that in the southwest Pingdingshan coalfield,the proportion of precipitation decreased gradually from the recharge area to the discharge area,ranging from 89.1%to 17.1%.Besides,the northeast Pingdingshan coalfield is another recharge area for the whole karst system,thus,the infiltrating groundwater will indirectly supply the deep aquifer through Quaternary deposition near the mine No.13.Our research results can provide theoretical support for the prevention and control of groundwater disasters and the development and utilization of regional groundwater resources in the coalfield in Northern China.

Impact of water table on hierarchically nested groundwater flow system

Water table configuration gives rise to hierarchically nested groundwater flow systems.However,there remains a lack of comprehensive understanding regarding the controlling factors of water table and its impact on flow systems.Moreover,it remains challenging to identify characteristics of water table space variation through limited groundwater observations at the regional scale.Based on two ideal two-dimensional cross-section analytical models,this study presents a simplified approach to preliminarily assess the nonlinear interactions between water table variation and three driving factors:Topography,geol-ogy and climate.Two criteria,C1 and C2,are utilized to address issues at different scales ranging from basin to local:(i)the influence of various factors on water table configuration;and(ii)the influence of water table on groundwater flow pattern.Then,the Ordos Plateau is taken as an example to explore the role of the water table in nested groundwater systems using the provided approach and criterion.The applica-tion of this approach in the Ordos Plateau demonstrates its appropriateness as a practical method for prelim-inarily determining the characteristics of water table configuration and its impact on flow systems.The study explores the mechanism influencing spatial variation in the water table and improves understanding of the interaction between topography,geology,and climate on groundwater flow patterns.

Introducing Fractal Dimension to Estimation of Soil Sensitivity to Preferential Flow

Food dye Brilliant Blue was introduced as the tracer in a dye-tracing experiment to obtain dye profile patterns of sandy loam soil, aeolian sandy soil, percolating paddy soil and permeable paddy soil. The dyed soil profiles were then photographed and the photos were scanned into a computer. Edited with certain software, only the dyed areas were left on the profile photos, which indicted the preferential flow paths for water and solute transport. Fractal dimensions of the dye patterns were calculated according to Arnold's function. Soil particle size distribution was analyzed by pipette method. The regression analysis showed that there was significant relationship between soil clay content and fractal dimension D of the dye pattern of soil profile. Based on the experiment results, the possibility of introducing fractal dimension to estimation of soil sensitivity to preferential flow is discussed.

Groundwater flow through fractured rocks and seepage control in geotechnical engineering: Theories and practices

Groundwater flow through fractured rocks has been recognized as an important issue in many geotechnical engineering practices.Several key aspects of fundamental mechanisms,numerical modeling and engineering applications of flow in fractured rocks are discussed.First,the microscopic mechanisms of fluid flow in fractured rocks,especially under the complex conditions of non-Darcian flow,multiphase flow,rock dissolution,and particle transport,have been revealed through a com-bined effort of visualized experiments and theoretical analysis.Then,laboratory and field methods of characterizing hydraulic properties(e.g.intrinsic permeability,inertial permeability,and unsaturated flow parameters)of fractured rocks in different flow regimes have been proposed.Subsequently,high-performance numerical simulation approaches for large-scale modeling of groundwater flow in frac-tured rocks and aquifers have been developed.Numerical procedures for optimization design of seepage control systems in various settings have also been proposed.Mechanisms of coupled hydro-mechanical processes and control of flow-induced deformation have been discussed.Finally,three case studies are presented to illustrate the applications of the improved theoretical understanding,characterization methods,modeling approaches,and seepage and deformation control strategies to geotechnical engi-neering projects.

Numerical simulation of the dimensional transformation of atomization in a supersonic aerodynamic atomization dust-removing nozzle based on transonic speed compressible flow

To simulate the transonic atomization jet process in Laval nozzles,to test the law of droplet atomization and distribution,to find a method of supersonic atomization for dust-removing nozzles,and to improve nozzle efficiency,the finite element method has been used in this study based on the COMSOL computational fluid dynamics module.The study results showed that the process cannot be realized alone under the two-dimensional axisymmetric,three-dimensional and three-dimensional symmetric models,but it can be calculated with the transformation dimension method,which uses the parameter equations generated from the two-dimensional axisymmetric flow field data of the three-dimensional model.The visualization of this complex process,which is difficult to measure and analyze experimentally,was realized in this study.The physical process,macro phenomena and particle distribution of supersonic atomization are analyzed in combination with this simulation.The rationality of the simulation was verified by experiments.A new method for the study of the atomization process and the exploration of its mechanism in a compressible transonic speed flow field based on the Laval nozzle has been provided,and a numerical platform for the study of supersonic atomization dust removal has been established.

Variation characteristics of aquifer parameters induced by groundwater source heat pump operation under variable flow

The variation characteristics of aquifer parameters,induced by groundwater source heat pump(GWSHP) operation under variable flow,were theoretically analyzed through a case study,in which the characteristics of building air conditioning load were considered.The results,compared with the constant flow operation,indicate that the influence on the variations of porosity,hydraulic conductivity and confined water head is decreased by 48%,51% and 71%,respectively,under variable flow operation.The security of variable flow operation is superior to that of constant flow.It is also concluded that the climate region and function of the buildings are primary factors which affect the suitability of variable flow operation in GWSHP.

Modified asymptotic Adomian decomposition method for solving Boussinesq equation of groundwater flow

The Adomian decomposition method （ADM） is an approximate analytic method for solving nonlinear equations. Generally, an approximate solution can be ob- tained by using only a few terms. However, in applications, we need to use it flexibly according to the real problem. In this paper, based on the ADM, we give a modified asymptotic Adomian decomposition method and use it to solve the nonlinear Boussinesq equation describing groundwater flows. The example shows effectiveness of the modified asymptotic Adomian decomposition method.

Changes of groundwater flow field of Luanhe River Delta under the human activities and its impact on the ecological environment in the past 30 years

The Luanhe River Delta is located in the center of the Circum-Bohai Sea Economic Zone.It enjoys rapid economic and social development while suffering relatively water scarcity.The overexploitation of groundwater in the Luanhe River Delta in recent years has caused the continuous drop of groundwater level and serious environmental and geological problems.This study systematically analyzes the evolution characteristics of the population,economy,and groundwater exploitation in the Luanhe River Delta and summarizes the change patterns of the groundwater flow regime in different aquifers in the Luanhe River Delta according to previous water resource assessment data as well as the latest groundwater survey results.Through comparison of major source/sink terms and groundwater resources,the study reveals the impacts of human activities on the groundwater resources and ecological environment in the study area over the past 30 years from 1990 to 2020.The results are as follows.The average annual drop rate of shallow groundwater and the deep groundwater in the centers of depression cones is 0.4 m and 1.64 m,respectively in the Luanhe River Delta in the past 30 years.The depression cones of shallow and deep groundwater in the study area cover an area of 545.32 km^(2)and 548.79 km^(2),respectively,accounting for more than 10%of the total area of the Luanhe River Delta.Overexploitation of groundwater has further aggravated land subsidence.As a result,two large-scale subsidence centers have formed,with a maximum subsidence rate of up to 120 mm/a.The drop of groundwater level has induced some ecological problems in the Luanhe River Delta area,such as the zero flow and water quality deterioration of rivers and continuous shrinkage of natural wetlands and water.Meanwhile,the proportion of natural wetland area to the total wetland area has been decreased from 99%to 8%and the water area from 1776 km^(2)to 263 km^(2).These results will provide data for groundwater overexploitation control,land subsidence prevention,and ecological restoration in plains and provide services for water resources management and national land space planning.

Numerical simulation of response of groundwater flow system in inland basin to density changes

The developmental characteristics of groundwater flow system are not only controlled by formation lithology and groundwater recharge conditions,but also influenced by the physical properties of fluids.Numerical simulation is an effective way to study groundwater flow system.In this paper,the ideal model is generalized according to the fundamental characteristics of groundwater system in inland basins of Western China.The simulation method of variable density flow on the development of groundwater system in inland basins is established by using EOS9 module in TOUGHREACT numerical simulation software.In accordance with the groundwater streamline,the groundwater flow system is divided into three levels,which are regional groundwater flow system,intermediate groundwater flow system and local groundwater flow system.Based on the calculation of the renewal rate of groundwater,the analysis shows that the increase of fluid density in the central part of the basin will restrain the development of regional groundwater flow system,resulting in a decrease of the circulation rate from 32.28% to 17.62% and a certain enhancement to the local groundwater flow system,which increased from 37.29% to 51.94%.

Groundwater Flow Modeling for Qushtapa Plain Unconfined Aquifer in Southern Erbil Basin, Kurdistan Region, Iraq

Increasing population growth and water demand for various purposes such as irrigation, domestic and industrial production in many parts of the Kurdistan Region is causing deficit in fresh water and rising groundwater dependence. Drilling many deep wells in the area unsystematically and continuously increased pumping water from groundwater reservoirs results in lowering of water table. Therefore, it is essential to assess the management of water resources. The study focuses on the groundwater modeling for the Qushtapa District plain area in particular under steady state flow conditions. The aquifer was simulated under unconfined condition and is represented by a single layer of 100 m thickness. MODPATH was used to measure contamination track lines and travel times. This approach involved the introduction of particles at sources of contaminants in the wells and the recharge area, then the identification of the path lines and the determination of the special distribution of contaminants through steady state flow conditions. The simulation of the groundwater head shows that the groundwater head starts from the northeastern part of the plain and decreases towards Lesser Zab River in the south of the plain from 420 m to 140 m above sea level. The modeled layer was calibrated under steady state conditions using hydraulic parameters obtained from observation and pumping wells. The calibrated model is effective in producing steady-state groundwater head distribution and good compliance with observed data. The standard error was estimated as 4.88 m, the normalized root mean square error is 8.3% and the residual mean is 15.79 m. The results of the forward tracking show the source of potential pollutants from the recharge area after different travel time, the particles released at the northern boundary travels to the center and the western part toward the pollution sources. The results of the backward tracking show that the particles located in the extraction wells moved toward the recharge area in the north and northeastern part of the study area.

Conventional and futuristic approaches for the computation of groundwater recharge:A comprehensive review

Groundwater recharge is a critical hydrologic component that determines groundwater availability and sustainability.Groundwater recharge estimation can be performed in a variety of ways,ranging from direct procedures to simulation models.The optimal strategy for recharge estimation depends on several factors,such as study objectives,climatic zones,hydrogeological conditions,data availability,methodology,and temporal and spatial constraints.Groundwater recharge is influenced by uncertainties in weather and hydrology.This study discusses conventional recharge estimation techniques and their application for optimal recharge calculation,and it also offers an overview of recent advances in recharge estimation methods.Most methods provide direct or indirect estimation of recharge across a small region on a point scale for a shorter time.With recent technological advancements and increased data availability,several advanced computational tools,including numerical,empirical,and artificial intelligence models,have been developed for efficient and reliable computation of groundwater recharge.This review article provides a thorough discussion of the techniques,assumptions,advantages,limitations,and selection procedures for estimating groundwater recharge.

Influence of underground space development mode on the groundwater flow field in Xiong’an new area

The degree and scale of underground space development are growing with the continuous advancement of urbanization in China.The lack of research on the change of the groundwater flow field before and after the development of underground space has led to various problems in the process of underground space development and operation.This paper took the key development zone of the Xiong’an New Area as the study area,and used the Groundwater modeling system software(GMS)to analyse the influence on the groundwater flow field under the point,line,and surface development modes.The main results showed that the underground space development would lead to the expansion and deepening of the cone of depression in the aquifer.The groundwater level on the upstream face of the underground structure would rise,while the water level on the downstream face would drop.The“line”concurrent development has the least impact on the groundwater flow field,and the maximum rise of water level on the upstream side of the underground structure is expected to be approximately 3.05 m.The“surface”development has the greatest impact on the groundwater flow field,and the maximum rise of water level is expected to be 7.17 m.

Review of the algebraic linear methods and parallel implementation in numerical simulation of groundwater flow

The desire to increase spatial and temporal resolution in modeling groundwater system has led to the requirement for intensive computational ability and large memory space. In the course of satisfying such requirement, parallel computing has played a core role over the past several decades. This paper reviews the parallel algebraic linear solution methods and the parallel implementation technologies for groundwater simulation. This work is carried out to provide guidance to enable modelers of groundwater systems to make sensible choices when developing solution methods based upon the current state of knowledge in parallel computing.

Boundary Element Analysis (Laplace Transform Solution) of Groundwater Unsteady Flow to a Multiple Well System in a Confined Aquifer

The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time varia-ble. The image function of water head H can be solved by BEM. We derived the boundary integral equation ofthe transformed variable H and the discretization form of it, so that there is no need to discretize the bounda-ries of well walls and it becomes easier to solve the groundwater head H by numerical inversion.

ON MULTI-DIMENSIONAL SONIC-SUBSONIC FLOW

In this paper, a compensated compactness framework is established for sonicsubsonic approximate solutions to the n-dimensional （n ≥ 2） Euler equations for steady irrotational flow that may contain stagnation points. This compactness framework holds provided that the approximate solutions are uniformly bounded and satisfy Hloc^-1（Ω） compactness conditions. As illustration, we show the existence of sonic-subsonic weak solution to n-dimensional （n ≥ 2） Euler equations for steady irrotational flow past obstacles or through an infinitely long nozzle. This is the first result concerning the sonic-subsonic limit for n-dimension （n ≥ 3）.

Quantitative Analysis of Groundwater Flow near a Partially Penetrating River under Riverside Pumping

According to practical geological and hydrogeological conditions of riverside water-supply well fields in northwestern China, an ideal hydrogeological model has been generalized and a three-dimensional mathematical model has been set up. A finite difference method was applied to simulating groundwater flow near a partially penetrating river under riverside pumping, and to analyzing the effects of river width, partial penetration and permeability of riverbed sediments on groundwater recharges. Results show that riverside pumping may cause groundwater to flow beneath the partially penetrating river, and that river width, penetration and riverbed permeability obviously influence flows from the partially penetrating river and constant-head boundaries. However, the pumping output is mainly from the partially penetrating river.

Groundwater remediation engineering——Study on the flow distribution of air sparging using acetylene

Air sparging(AS) is an emerging method to remove VOCs from saturated soils and groundwater. Air sparging performance highly depends on the air distribution resulting in the aquifer. In order to study gas flow characterization, a two-dimensional experimental chamber was designed and installed. In addition, the method by using acetylene as the tracer to directly image the gas distribution results of AS process has been put forward. Experiments were performed with different injected gas flow rates. The gas flow patterns were found to depend significantly on the injected gas flow rate, and the characterization of gas flow distributions in porous media was very different from the acetylene tracing study. Lower and higher gas flow rates generally yield more irregular in shape and less effective gas distributions.

Groundwater Flow Model for a Tannery Belt in Southern India

The objective of this article is to develop a groundwater flow model for a tannery belt using Visual MOD-FLOW Premium 4.4 for analyzing groundwater velocity and its response to various pumping strategies in two stages, viz., steady and transient conditions. The steady state model was calibrated for April 2001, whereas the transient model was employed to forecast groundwater flow under various pumping strategies. The results showed that the total groundwater abstraction was about 80.43% of the groundwater recharge, but 10.25% was used up by evapotranspiration. The groundwater velocity, which is important for contaminant migration, varied from 0.21 to 0.52 m/d in the tannery cluster. The model was more sensitive to recharge from rainfall, hydraulic conductivity and specific yield. Finally, the model showed that the aquifer could sustain a pumping rate of 24892 m3/day without further decline in water level.