连云港徐圩地面沉降BOTDR监测与评价

连云港市位于苏北沿海地区,地面沉降灾害面积较大,多地沉降速率超过20 mm/a,徐圩的沉降现状尤为严重。为了能够对徐圩地区的地面沉降进行精细化观测,文章采用BOTDR分布式光纤感测技术,对徐圩镇127 m深的钻孔地层进行了两年多的全断面精细化监测。结果表明:徐圩镇共有四个承压含水层组,I-1隔水层和I-2隔水层土体沉降量分别占总沉降量的70.29 %、24.59 %,抽水层的土体最大沉降量仅占比1.38 %。I-1隔水层和I-2隔水层的地层岩性包括淤泥质黏土（L2）、亚黏土（L3）、亚砂土夹粉砂（L4）,总厚度为44 m,由于抽水过程中隔水层向含水层失水,导致该隔水层土体固结压缩。同时,工程建设附加荷载对地面沉降的影响也不可忽视。徐圩地区现阶段的沉降仍在继续发生,但沉降速率有减小的趋势。BOTDR技术可有效获取地面沉降钻孔全断面的土层变形分布信息,为地面沉降评价提供了一种精细化的分布式监测手段。     

Groundwater occurrence and aquifer vulnerability assessment of Magodo Area,Lagos, Southwestern Nigeria

A combination of vertical electrical soundings (VES), 2D electrical resistivity imaging (ERI) surveys and borehole logs were conducted at Magodo, Government Reserve Area (GRA) Phase 1, Isheri, Southwestern Nigeria, with the aim of delineating the different aquifers present and assessing the groundwater safety in the area. The Schlumberger electrode array was adopted for the VES and dipole-dipole array was used for the 2D imaging. The maximum current electrode spread (AB) was 800 m and the 2D traverse range between 280 and 350 m in the east-west direction. The thickness of impermeable layer overlying the confined aquifer was used for the vulnerability ratings of the study area. Five lithological units were delineated: the topsoil, clayey sand, unconsolidated sand which is the first aquifer, a clay stratum and the sand layer that constitutes the confined aquifer horizon. The topsoil thickness varies from 0.6 to 2.6 m, while its resistivity values vary between 55.4 and 510.6 Ω/m. The clayey sand layers have resistivity values ranging from 104.2 to 143.9 Ω/m with thickness varying between 0.6 and 14.7 m. The resistivity values of the upper sandy layer range from 120.7 to 2195.2 Ω/m and thickness varies from 3.3 to 94.0 m. The resistivity of the clay layer varies from 11.3 to 96.1 Ω/m and the thickness ranges from 29.6 to 76.1 m. The resistivity value of the confined aquifer ranges between 223 and 1197.4 Ω/m. The longitudinal conductance (0.0017–0.02 mhos) assessment of the topsoil shows that the topsoil within the study area has poor overburden protective capacity, and the compacted impermeable clay layer shows that the underlying confined aquifer is well protected from contamination and can be utilized as a source of portable groundwater in the study area. This study therefore enabled the delineation of shallow aquifers, the variation of their thicknesses and presented a basis for safety assessment of groundwater potential zones in the study area.     

Experimental study on the vertical deformation of aquifer soils under conditions of withdrawing and recharging of groundwater in Tongchuan region,China

Land subsidence due to aquifer-system compaction accompanying groundwater extraction is a global hazard. Rising urban construction and groundwater demand necessitate increased awareness and better understanding of the geological problem. Motivated by the lack of laboratory-scale studies on this issue, an experimental investigation on the newly developed Tongchuan region, China, is presented. The study addresses the deformation behaviors of three soil samples, with the lithology of silty clay, silt, and fine sand, under the conditions of groundwater withdrawal and recharge using the GDS Consolidation Testing System. Results indicate that all three samples were characterized by elastic–plastic deformation under the conditions of withdrawing and recharging. The vertical deformation of the silty clay in the aquitard above the first confined aquifer was larger than those of the other two samples, and its deformation is a gradual and long process; thus, considerable attention should be paid to deformation in this aquitard due to the apparent creep effect and tiny rebound deformation. However, the settlement of the fine sand in the second confined aquifer cannot be ignored due to the great thickness of the aquifer. For the same soil, as the pore-water pressure declined, the unit rate of vertical deformation decreased gradually, whereas the creep effect of deformation in the later declining stage of pore-water pressure was more apparent than that in the former declining stage. These observations are highly important to the local government, which is developing measures to prevent and control subsidence.     

Hydrogeophysical investigations for assessing the groundwater potential in part of the Peshawar basin,Pakistan

A large number of valleys and basin systems are present in the northwestern part of the Himalayas in Pakistan which form significant aquifers in the region. Hydrogeophysical investigations in the western part of Nowshera District, a part of the intermontane Peshawar basin, were undertaken to help to determine the availability of groundwater resources in the region. Thirty vertical electrical resistivity soundings (VES) were acquired using a Schlumberger expanding array configuration with a maximum current electrode spacing (AB/2) of 150 m in delineating the groundwater potential in the study area. The results of the interpreted VES data using a combination of curve matching technique and computer iterative modeling methods suggest that the area is underlain by 3 to 5 geo-electric layers. The interpretation results showed that the geo-electrical succession consists of alluvium comprising of alternating layers of clay, silty clay, fine to coarse sands, sand with gravels and gravels of variable thickness. High subsurface resistivity values are correlated with gravel–sand units and low resistivity values with the presence of clays and silts. The modeled VES results were correlated with the pumping tests results and lithological logs of the existing wells. The pumping test suggests the transmissivity of the aquifer sediments is variable corresponding to different sediments within the area. The gravel–sand intervals having high resistivity value show high transmissivity values, whereas clay–silt sediments show low transmissivities. It is concluded that majority of the high resistive gravel–sand sediments belong to an alluvial fan environment. These gravel–sand zones are promising zones for groundwater abstraction which are concentrated in the central part of the study area.     

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

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.     

Disposition of aquifer system, geo-electric characteristics and gamma-log anomaly in the Kuttanad alluvium of Kerala

Kuttanad near the western coast of Kerala is a low lying area below mean sea level. The hydrogeological scenario in this wet land dominated tract, known as kole lands, was investigated by using the lithological and geophysical data generated under the R&D exploration programme of Central Ground Water Board (CGWB) during the year 2000, as also other related data generated till the year 2009. The hydrogeology of the tract is complex, and the presence of carbonaceous and organic matter in the soil and sub-soil influences the geochemical and geo-electric characteristics of the sub-surface geological formations as reflected by SP and resistivity logs. The lithologs and geophysical logs of bore wells were used to define the aquifer geometry. While the lateral and vertical variations of the aquifer systems identified were high, the bottom aquifer was relatively brackish over a large area. The gamma log behavior in these granular zones is anomalous and misguiding in nature. The radioactivities recorded against clay horizons were in the range of 130 to 200cps and in the granular zones 150 to 350 cps. The presence of monazite bearing sands at certain depths was responsible for anomalously high gamma counts against granular zones. The water quality in the granular zones was inferred from the electrical logs. The depth of occurrence of high radiation sand beds was also demarcated.     

Well log data analysis for lithology and fluid identification in Krishna-Godavari Basin,India

Well log analysis provides the information on petrophysical properties of reservoir rock and its fluid content. The present study depicts interpretation of well log responses such as gamma ray, resistivity, density and neutron logs from six wells, namely W-1, W-2, W-9, W-12, W-13 and W-14 under the study area of Krishna-Godavari (K-G) basin. The logs have been used primarily for identification of lithology and hydrocarbon-bearing zones. The gamma ray log trend indicates deposition of cleaning upward sediment. Coarsening upward, clayey-silty-sandy bodies have been evidenced from the gamma ray log. Gas-bearing zones are characterised by low gamma ray, high deep resistivity and crossover between neutron and density logs. Total 14 numbers of hydrocarbon-bearing zones are identified from wells W-9, W-12, W-13 and W-14 using conventional log analysis. Crossplotting techniques are adopted for identification of lithology and fluid type using log responses. Crossplots, namely P-impedance vs. S-impedance, P-impedance vs. ratio of P-wave and S-wave velocities (Vp/Vs) and lambda-mu-rho (LMR), have been analysed to discriminate between lithology and fluid types. Vp/Vs vs. P-impedance crossplot is able to detect gas sand, brine sand and shale whereas P-impedance vs. S-impedance crossplot detects shale and sand trends only. LMR technique, i.e. λρ vs. μρ crossplot is able to discriminate gas sand, brine sand, carbonate and shale. The LMR crossplot improves the detectability and sensitivity of fluid types and carbonate lithology over other crossplotting techniques. Petrophysical parameters like volume of shale, effective porosity and water saturation in the hydrocarbon-bearing zones in these wells range from 5 to 37%, from 11 to 36 and from 10 to 50% respectively. The estimated petrophysical parameters and lithology are validated with limited core samples and cutting samples from five wells under the study area.     

Groundwater potential assessment in a sedimentary terrain,southwestern Nigeria

Geoelectrical resistivity techniques are increasingly being applied in addressing a wide range of hydrological, environmental, and geotechnical problems. This is due to their effectiveness in near-surface characterization. In the present study, a suite of vertical electrical soundings (VESs) was integrated with 2D geoelectrical resistivity and time-domain induced polarization (IP) imaging to characterize the near-surface and delineate the underlying aquifer in a sedimentary terrain. The geophysical survey was conducted as part of preliminary studies to evaluate the potential of groundwater resource in Iyana-Iyesi and Canaan Land area of Ota, southwestern Nigeria. A high-yield confined sandy aquifer overlain by a low-yield aquitard was delineated; overlying the aquitard is a very resistive and thick layer that is predominantly composed of kaolinitic swelling clay intercalated with phosphate mineral.     

Electrical imaging of the groundwater aquifer at Banting,Selangor, Malaysia

A geophysical study was carried out in the Banting area of Malaysia to delineate groundwater aquifer and marine clay layer of the alluvial Quaternary deposits of Beruas and Gula Formations. The Beruas Formation is formed by peat and clayey materials as well as silt and sands, whereas the Gula Formation consists of clay, silt, sand and gravels. Both Formations were deposited on top of the Carboniferous shale of the Kenny Hill Formation. A 2-D geoelectrical resistivity technique was used. Resistivity measurement was carried out using an ABEM SAS 4000 Terrameter. The 2-D resistivity data of subsurface material for each survey line was calculated through inverse modelling and then compared with borehole data. The resistivity images of all the subsurface material below the survey lines show similar pattern of continuous structure of layering or layers with some lenses with resistivity ranging from 0.1 to 50 Ωm. The upper layer shows resistivity values ranging from 0.1 to 10 Ωm, representing a clay horizon with a thickness up to 45 m. The second layer with depth varies from 45 to 70 m below surface and has resistivity values ranging from 10 to 30 Ωm. Borehole data indicate coarse sand with some gravels for this layer, which is also the groundwater aquifer in the study area. The lowermost layer at a depth of 70 m below ground level shows resistivity values ranging from 30–50 Ωm and can be correlated with metasedimentary rocks consisting of shale and metaquartzite.     

Vertical electrical sounding (VES) and multi-electrode resistivity in environmental impact assessment studies over some selected lakes: a case study

A combined geophysical investigation consisting of vertical electrical sounding (VES) and multielectrode system was carried out to map the subsurface resistivity in all major lakes which are highly polluted by the discharge of sewage and other chemical effluents in greater Hyderabad, India. The structural features identified in the study area play a major role in groundwater flow and storage. The interpretation of geophysical data and lithologs indicates that a silt/clay zone (predominantly silt) has a thickness of 5–10 m all along the drainage from Patelcheruvu to the Musi River. The silt/clay zone inferred close to the lakes is a mixture of clay, silt and sand with more silt content as indicated from the lithologs during drilling. The low resistivity values obtained can be attributed to the pollutant accumulated in the silt which can reduce the resistivity values. Further, the TDS of the water samples in these wells are more than 1,000 mg/l which further confirms the above scenario. The pollution spread is less in the upstream areas whereas it is more in the downstream which can be attributed to the shallow water table conditions and also due to the interaction of surface water and groundwater.     

Application of electrical resistivity to map the stratigraphy and salinity of fluvio-deltaic aquifers: case studies from Bangladesh that reveal benefits and pitfalls

Fluvio-deltaic aquifers are the primary source of drinking water for the people of Bangladesh. Such aquifers, which comprise the Ganges-Brahmaputra-Meghna Delta, are hydrogeologically heterogeneous. Because of widespread groundwater quality issues in Bangladesh, it is crucial to know the hydrostratigraphic architecture and hydrochemistry, as some aquifer units are contaminated, whereas others are safe. Geophysical methods provide a potentially effective and noninvasive method for extensive characterization of these aquifers. This study applies and investigates the limitations of using electrical resistivity imaging (ERI) for mapping the hydrostratigraphy and salinity of an aquifer-aquitard system adjacent to the Meghna River. Some electrical resistivity (ER) sections showed excellent correlation between resistivity and grain size. These suggest that ERI is a powerful tool for mapping internal aquifer architecture and their boundaries with finer-grained aquitards which clearly appear as low-ER zones. However, in parts of some ER sections, variations in electrical properties were determined by porewater resistivity. In these cases, low ER was indicative of brine and did not indicate the presence of finer-grained materials such as silt or clay. Accordingly, the following hydrostratigraphic zones with different resistivities were detected: (1) aquifers saturated with fresh groundwater, (2) a regional silt/clay aquitard, and (3) a deeper brine-saturated formation. In addition, shallow silt/clay pockets were detected close to the river and below the vadose zone. ERI is thus a promising technique for mapping aquifers versus aquitards; however, the observations are easily confounded by porewater salinity. In such cases, borehole information and groundwater salinity measurements are necessary for ground-truthing.

     

Mapping an aquitard breach using shear-wave seismic reflection

In multi-layered hydrostratigraphic systems, aquitard breaches caused by faulting or paleo-erosion can allow substantial quantities of water of differing quality to be exchanged between aquifers. Seismic reflection technology was used to map the extent and orientation of an aquitard breach connecting a shallow alluvial aquifer to the deeper semi-confined Memphis aquifer in southwestern Tennessee, USA. Geophysical well logs indicate the presence of the aquitard at borehole locations that define the beginning and end points on two seismic survey lines, which intersect at a borehole where the aquitard is absent. A SE–NW-oriented paleochannel, 350 m wide and approximately 35–40 m deep, is interpreted from the seismic reflection surveys. The paleochannel cuts through the aquitard and into the upper part of the Memphis aquifer, thus creating a hydraulic connection between the shallow unconfined and deeper, semi-confined aquifers. The results indicate the potential of the shear-wave seismic reflection methods to resolve shallow breaches through fine-grained aquitards given availability of sufficient well control.     

Groundwater investigation of the artesian wells on the palaeochannels in parts of the Great Rann of Kachchh,Gujarat, India,using remote sensing and geophysical techniques

The Great Rann of Kachchh (GRK) in Gujarat, India, is the largest salt desert in the world, which is usually filled with seawater ingression during high tide from the Arabian Sea. As a result, the soil gets saturated with saline water that has percolated down for several meters. Groundwater exploration in Rann area is a challenging task due to the prevailing hostile environment. For this purpose, multisensor satellite data have been used to delineate the palaeochannels in search of an alternate source of drinking water. In GRK, palaeochannels represent the zone of elevated fluvial sediments with respect to the surroundings. Evolutionary history of the palaeochannels indicates upliftment of GRK area during Allah Bund faulting. For assessing the groundwater potential of the palaeochannels, high-resolution electrical resistivity tomography (HERT) surveys have been carried out with the pole-dipole method. Electrical resistivity tomograms along 710 m traverses to a depth of 250 m in Dharmsala and Gainda area show higher-resistivity zones (medium to coarse sand with brackish water) below a thick low-resistivity layer (clay with saline water). A few exploratory drillings in the area confirm the existence of the palaeochannels, which act as a confined aquifer below 100 m depth. The artesian condition of the two drilled wells at Gainda and Khardoi along the northern boundary of GRK may be attributed to hydraulic gradient along the confined layers from the Tharparkar region in Pakistan. Thus, HERT is found to be a faster and more cost-effective geophysical survey technique for study of the deep aquifer.

     

Aquifer system response to intensive pumping in urban areas of the Gangetic plains,India: the case study of Patna

A significant component of domestic demand for water of urban areas located in the Gangetic plains is met by heavy pumping of groundwater. The present study is focused on the Patna municipal area, inhabited by 17 million people and spanning over 134 km², where entire urban water demand is catered from pumping by wells of various capacities and designs. The present study examines the nature of the aquifer system within the urban area, the temporal changes in the water/piezometric level and the recharge mechanism of the deeper aquifers. The aquifer system is made up of medium-to-coarse unconsolidated sand, lying under a ~40-m-thick predominantly argillaceous unit holding 8- to 13-m-thick localised sand layers and continues up to 220 m below ground. Groundwater occurs under semi-confined condition, with transmissivity of aquifers in 5,500–9,200 m² day^?1 range. Hydraulic head of the deeper aquifer remains in 9–19 m range below ground, in contrast to 1–9 m range of that of the upper aquitard zone. The estimated annual groundwater extraction from the deeper aquifer is ~212.0 million m³, which has created a decline of 3.9 m in the piezometric level of the deeper aquifer during the past 30 years. Unregulated construction of deep tube wells with mushrooming of apartment culture may further exacerbate the problem. The sand layers within the aquitard zone are experiencing an annual extraction of 14.5 million m³ and have exhibited stable water level trend for past one and half decades. This unit is recharged from monsoon rainfall, besides contribution from water supply pipe line leakage and seepage from unlined storm water drains.     

砂泥岩薄互层测井资料解释方法

从分析薄互层的特征入手,采用反褶积方法对单条测井曲线做提高分辨率的处理,并调整电阻率测井曲线,使其与提高了分辨率的泥质含量及孔隙度测井曲线的纵向分辨率协调、在此基础上,按一定的解释模型计算出孔隙度和饱和度等参数。实际计算表明,本方法改善了储集参数的评价,降低了含水饱和度。     

基于多源数据的弱透水层水文地质参数反演研究

弱透水层是含水层系统的重要组成部分,其水文地质参数的准确获取一直以来都是研究难点,传统室内试验难以克服应力状态及环境变化对土层参数的影响,因此如何在原位状态下评价其水文地质参数有着重要意义。选取呼和浩特盆地某淤泥层为研究对象,基于抽水试验和溶质运移试验获得的多源数据,采用基于控制体积有限差分法的非结构化网格（Unstructured Grid）方法,建立三维地下水流-溶质运移耦合模型,对弱透水层水文地质参数进行反演。结果显示:模型计算的抽水井中水位和溶质浓度变化规律与实际观测数据能够较好拟合,模型识别出的抽水含水层渗透参数（4.8 m/d）与通过解析法计算得出的结果（4.17 m/d）相近,建立的数值模型符合实际水文地质条件;反演得到淤泥质弱透水层垂向渗透系数为1.2×10?4 m/d、储水率为1.0×10?5 m?1。基于多源数据的参数反演方法可为弱透水层参数研究提供借鉴。     

Numerical groundwater-flow modeling to evaluate potential effects of pumping and recharge: implications for sustainable groundwater management in the Mahanadi delta region,India

Process-based groundwater models are useful to understand complex aquifer systems and make predictions about their response to hydrological changes. A conceptual model for evaluating responses to environmental changes is presented, considering the hydrogeologic framework, flow processes, aquifer hydraulic properties, boundary conditions, and sources and sinks of the groundwater system. Based on this conceptual model, a quasi-three-dimensional transient groundwater flow model was designed using MODFLOW to simulate the groundwater system of Mahanadi River delta, eastern India. The model was constructed in the context of an upper unconfined aquifer and lower confined aquifer, separated by an aquitard. Hydraulic heads of 13 shallow wells and 11 deep wells were used to calibrate transient groundwater conditions during 1997–2006, followed by validation (2007–2011). The aquifer and aquitard hydraulic properties were obtained by pumping tests and were calibrated along with the rainfall recharge. The statistical and graphical performance indicators suggested a reasonably good simulation of groundwater flow over the study area. Sensitivity analysis revealed that groundwater level is most sensitive to the hydraulic conductivities of both the aquifers, followed by vertical hydraulic conductivity of the confining layer. The calibrated model was then employed to explore groundwater-flow dynamics in response to changes in pumping and recharge conditions. The simulation results indicate that pumping has a substantial effect on the confined aquifer flow regime as compared to the unconfined aquifer. The results and insights from this study have important implications for other regional groundwater modeling studies, especially in multi-layered aquifer systems.     

Assessment of groundwater occurrence in Olomoro,Nigeria using borehole logging and electrical resistivity methods

The objective of this study was to assess the subsurface strata and groundwater situation of Olomoro, Nigeria using borehole logging and electrical resistivity techniques. The borehole logging consisting of resistivity and spontaneous potential logs were conducted by using the Johnson Keck logger on a drilled well in the study area. The electrical resistivity survey involving 17 vertical electrical soundings (VES) with a maximum current electrode spacing of 100 to 150 m was conducted using the Schlumberger electrode configuration. Analysis of the well cuttings revealed that the lithology of the subsurface consist of topsoil, clay, very fine sand, medium grain sand, coarse sand and very coarse sand. Results of the downhole logging also revealed that the mean electrical conductivity and the total dissolved solid of the groundwater was obtained as 390 μS/cm and 245 mg/cm³ respectively. These values are within the acceptable limit set by the Standard Organization of Nigeria (SON) for drinking water. The result of the vertical electrical sounding interpreted using the computer iterative modeling revealed the presence of four to five geoelectric layers which showed a close correlation with result from the lithology and downhole logging. Results further showed that the resistivity of the subsurface aquifer ranged between 1584 and 5420 Ωm while the aquifer depths varied between 27.8 and 39.3 m. Groundwater development of the area is suggested using the depth and resistivity maps provided in this study.     

Geoelectrical resistivity survey for the evaluation of hydrogeological condition of Bagerhat Sadar and the adjacent areas, Bangladesh

This study has been conducted to find out the aquifer characteristics of Bagerhat Sadar and adjacent areas in Bangladesh using geoelectrical resistivity method and borehole logs. The interpretation of resistivity soundings (14 nos.) shows that the sub-surface lithological sequence can be divided into four geoelectric units. The deepest geoelectric unit (with resistivity from 8 Ohm-m to 18 Ohm-m) represents the deep aquifer with usually fresh water. Resistivity values of 12 ?m or more for this unit may indicate formation water as acceptable for coastal people. In the study area, the shallow aquifer inferred is not suitable for groundwater development.     

Application of mathematical modelling techniques to the alluvial aquifer system near Wagga Wagga,New South Wales

The central part of the Murrumbidgee River meanders across an alluvial plain which is 1.5 km to 5 km wide, and confined between low hills composed of granite and folded lower Palaeozoic rocks. Deposition commenced under humid and reducing conditions during the Miocene, and continued during Pliocene time. Oligomict quartz‐sand and gravel beds are dominant in the sequence, and are interlayered with grey to black clay and silty clay, which contain a prolific pollen fauna. During Pleistocene time, a marked reduction in rainfall and consequent change in vegetation resulted in the deposition of a widespread coarse polymict poorly‐sorted gravel associated with yellow and brown clay and silt.

The quartz‐sand and gravel beds within the Pliocene deposits constitute a major aquifer system and contain low salinity water suitable for irrigation and town water‐supply purposes. Since 1968, increasing demand has been made on this aquifer system, and wells sunk for various State and Local Government bodies have a total installed capacity of approximately 1000 litres/sec.

The most concentrated exploration and development has been effected east of Wagga Wagga, where a well field was constructed to provide water for a regional water supply scheme. Following a refraction seismic survey, test drilling, and mathematical simulation, four production wells were constructed.

Mathematical simulation models developed on a mini‐computer were designed primarily to study the behaviour of the aquifer system in response to pumping from the Pliocene aquifers. An analytical model was used to assist in the location of production wells, and a numerical model was used, after construction of the wells, to determine what additional losses from the Murrumbidgee River would be induced by pumping them. These losses were found to correspond to approximately 90% of the pumping rate during long term pumping.