Geophysical methods have been applied to a wide range of hydrogeological problems. With improvement in geophysical inversion algorithms and measurement tools, significant achievements have been made in the characteriz...Geophysical methods have been applied to a wide range of hydrogeological problems. With improvement in geophysical inversion algorithms and measurement tools, significant achievements have been made in the characterization of subsurface architecture, time-lapse monitoring of hydrogeological process and contaminant plumes delineation. In this paper, we summarize the geophysical methods that are most widely used in hydrogeology including Electrical Resistivity Tomography(ERT), Induced Polarization(IP), Ground Penetrating Radar(GPR) and Electromagnetic Induction(EMI). Three examples including lab and field works are used to demonstrate current application of geophysical methods for characterizing subsurface architecture and contaminant plumes. Though great progress has been made in hydrogeohysics over the last few decades at home and abroad, challenges still remain in practical applications. More recently, hydrogeophysics continues to develop in the areas of establishment of hydrogeophysical models, large-scale architecture characterization, uncertainty analysis, biogeochemical process monitoring and ecosystem science.展开更多
The dynamic monitoring of landslides in engineering geology has focused on the correlation among landslide stability,rainwater infiltration,and subsurface hydrogeology.However,the understanding of this complicated cor...The dynamic monitoring of landslides in engineering geology has focused on the correlation among landslide stability,rainwater infiltration,and subsurface hydrogeology.However,the understanding of this complicated correlation is still poor and inadequate.Thus,in this study,we investigated a typical landslide in southwestern China via time-lapse electrical resistivity tomography(TLERT) in November 2013 and August 2014.We studied landslide mechanisms based on the spatiotemporal characteristics of surface water infiltration and flow within the landslide body.Combined with borehole data,inverted resistivity models accurately defined the interface between Quaternary sediments and bedrock.Preferential flow pathways attributed to fracture zones and fissures were also delineated.In addition,we found that surface water permeates through these pathways into the slipping mass and drains away as fissure water in the fractured bedrock,probably causing the weakly weathered layer to gradually soften and erode,eventually leading to a landslide.Clearly,TLERT dynamic monitoring can provide precursory information of critical sliding and can be used in landslide stability analysis and prediction.展开更多
Geophysical techniques can help to bridge the inherent gap that exists with regard to spatial resolution and coverage for classical hydrological methods. This has led to the emergence of a new and rapidly growing rese...Geophysical techniques can help to bridge the inherent gap that exists with regard to spatial resolution and coverage for classical hydrological methods. This has led to the emergence of a new and rapidly growing research domain generally referred to as hydrogeophysics. Given the differing sensitivities of various geophysical techniques to hydrologically relevant parameters, their inherent trade-off between resolution and range, as well as the notoriously site-specific nature of petrophysical parameter relations, the fundamental usefulness of multi-method surveys for reducing uncertainties in data analysis and interpretation is widely accepted. A major challenge arising from such endeavors is the quantitative integration of the resulting vast and diverse database into a unified model of the probed subsurface region that is consistent with all available measurements. To this end, we present a novel approach toward hydrogeophysical data integration based on a Monte-Carlo-type conditional stochastic simulation method that we consider to be particularly suitable for high-resolution local-scale studies. Monte Carlo techniques are flexible and versatile, allowing for accounting for a wide variety of data and constraints of differing resolution and hardness, and thus have the potential of providing, in a geostatistical sense, realistic models of the pertinent target parameter distributions. Compared to more conventional approaches, such as co-kriging or cluster analysis, our approach provides significant ad- vancements in the way that larger-scale structural information eontained in the hydrogeophysieal data can be accounted for. After outlining the methodological background of our algorithm, we present the results of its application to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the detailed local-scale porosity structure. Our procedure is first tested on pertinent synthetic data and then applied to a field dataset collected at the Boise Hydrogeophysical Research Site. Finally, we compare the performance of our data integration approach to that of more conventional methods with regard to the prediction of flow and transport phenomena in highly heterogeneous media and discuss the implications arising.展开更多
This research applies Electrical Resistivity Tomography (ERT) and Time Domain Electromagnetic Method (TDEM) to study the hydrogeology of the Taubaté basin, which is characterized by half-grabens with about 850 m ...This research applies Electrical Resistivity Tomography (ERT) and Time Domain Electromagnetic Method (TDEM) to study the hydrogeology of the Taubaté basin, which is characterized by half-grabens with about 850 m of maximum sediments thickness. The study area is in Taubaté city, São Paulo State, Brazil, where the Taubaté aquifer is an important water source. The Taubaté Group is the main sedimentary package of the basin;it is formed mainly by shales that form aquicludes, and thin layers of sandstones that form the aquifer. There are 40 groundwater exploration wells in Taubaté city that provide important information. The study purpose is to characterize the geoelectrical stratigraphy of the subsurface to locate the contact between the Quaternary and Tertiary sediments and to identify the Taubaté aquifer. The ERT is used for shallow investigations (tens of meters) and the TDEM can reach a great investigation depth (hundreds of meters). Therefore, these geophysical methods are complementary. The ERT data were acquired with the pole-dipole array with 20 m of electrodes spacing and 400 m length, and the TDEM data with the central-loop array with a 200 × 200 m transmitter loop. The results permit to define the contact between the Quaternary and Tertiary sediments around 15 m depth, the Pindamonhangaba Formation between 15 m and 30 m depth and the Taubate Group between 30 m and 300 m depth. The TDEM method defined the Taubaté Group as a single geoelectric layer because the shale and the sandstone layers are all very conductive. The basement is formed by gneiss, which is a very resistive rock. The TDEM method is not able to identify a high conductor/resistor contrast. Overall, the results are consistent with the known geology and the wells information.展开更多
This work uses 2D TEM (Transient Electromagnetic) modeling for a hydrogeological study in the Paraná sedimentary basin. The study area is located at the northern region of the state of S?o Paulo, Brazil, where gr...This work uses 2D TEM (Transient Electromagnetic) modeling for a hydrogeological study in the Paraná sedimentary basin. The study area is located at the northern region of the state of S?o Paulo, Brazil, where groundwater is exploited from two aquifer systems: one sedimentary, shallow, and the other crystalline, deep. The interest in applying the TEM method in this area owes to the high exploitation rates of groundwater from the crystalline aquifer system for irrigation, which is triggering considerable seismic activity locally. This aquifer system is composed of fractured basalt within the Serra Geral Formation and is about 120 m deep. Eighty-six TEM soundings were acquired at this location, but in nine cases the data did not fit the modelled curve for 1D geoelectrical models due to the geological complexity of the area. This paper shows 2D geoelectrical modeling results based on the FDTD (Finite Differences in Time Domain) method to explain the lateral resistivity variation within the geological setting. A 2D model was generated for each sounding and compared with 1D inversion models as well as with direct information from wells. The results show some vertical variations of about 10 to 30 meters on the upper interface of the basalt layer from Serra Geral Formation. They are located at approximately 60 meters from the center of the soundings. The existence of these 2D structures in the subsurface can be related to the drainage system in the study area. The presence of these structures may indicate a connection between the shallow and deep aquifer systems, acting like a conduit that may contribute to the seismic activity reported.展开更多
This study investigates the groundwater aquifer located in Fayuim oasis. In this study, two of the electromagnetic measurement methods have been used in determining the hydrological situation in the Fayoum oasis. The ...This study investigates the groundwater aquifer located in Fayuim oasis. In this study, two of the electromagnetic measurement methods have been used in determining the hydrological situation in the Fayoum oasis. The first is airborne electromagnetic (AEM) which, sometimes is referred to as Helicopter electromagnetic (HEM) and the second is ground Time-domain Electromagnetic method (TEM). The subsurface consists of four geoelectrical layers with a rough slope towards the center. The third and the fourth layers in the succession are suggested to be the two-groundwater aquifers. The third layer saturates with fresh water overlying saline water which exists in the bottom of the second one. It is worth mentioning that the depth of the fresh water surface undulates between the surface level in two lakes in the study area and 57 meters below the ground, whereas the thickness of the fresh water aquifer varies from 13 to 36 meters. The depth of the saline water surface undulates between 59 and 81 meters below the ground. In general, airborne electromagnetic surveying has the advantage of fast resistivity mapping with high lateral resolution. Groundbased geophysical surveys are often more accurate, but they are definitely slower than airborne surveys. It depends on targets of interest, time, budget, and manpower available by the method or the combination of methods that will be chosen. A combination of different methods is useful to obtain a detailed understanding of the subsurface resistivity distribution.展开更多
Geophysical study and watershed hydrological delineation have been integrated at downstream of Alasra dam site Norh Azraq area to investigate their potential for artificial groundwater recharge. The total surface area...Geophysical study and watershed hydrological delineation have been integrated at downstream of Alasra dam site Norh Azraq area to investigate their potential for artificial groundwater recharge. The total surface area of the watershed was found to be about 195 square kilometers. The estimated annual runoff volumes for the Alasra watershed ranged between 1.2 and 1.8 MCM. Moreover, the interpretation of Ten Time Domain Electromagnetic (TDEM) soundings suggested three principal subsurface layers. The top surface layer has an intermediate resistivity (90 - 110 Ohm·m) with a thickness ranging from a few meters to around 50 m. This layer was interpreted as superficial deposits. The second subsurface layer with variably high resistivity values is composed of unsaturated massive basalt layer and probably belongs to Madhala Olivine Phyric Basalt Formation (MOB). The large variations in resistivity could be ascribed to the degree of water saturation (as a result of groundwater recharge from the nearby harvested water dam), or lithological variations (clay content) and/or due to structural control. The third subsurface layer has low resistivity values (·m to 40 Ω·m) and was found at a depth ranging from 120 to 150 m. This layer could represent a saturated basalt layer with high clay contents. The subsurface structures and major faults have been identified. Based on the results of this study, a combination of surface and subsurface artificial groundwater recharge techniques is highly recommended.展开更多
Detailed local geological and hydrogeophysical investigations were carried out for the aquifer in Yaoundé, Cameroon to delineate the architecture of different subsurface geological horizons using lithologs and ge...Detailed local geological and hydrogeophysical investigations were carried out for the aquifer in Yaoundé, Cameroon to delineate the architecture of different subsurface geological horizons using lithologs and generated vertical electrical sounding (VES) data. An attempt has also been made to estimate aquifer transmissivity from resistivity data. The transmissivity of the uncon?ned aquifer was computed by determining the Dar-Zarrouk parameters (longitudinal unit conductance and transverse unit resistance) and were compared with the actual field transmissivity. The results showed a direct relation between aquifer transmissivity and transverse resistance. The relationship established has therefore, been generalized in the study area in order to evaluate hydraulic conductivity and transmissivity at all the points where geoelectrical measurements have been carried out. This generalization allows one to derive maps of the product Kσ and transmissivity in the study area based on geoelectrical measurements. These maps are important in future modelling processes oriented towards better exploitation of the aquifers.展开更多
基金funded by the National Natural Science Fund of China (NSFC)-Xinjiang No.U1503282the NSFC No.41030746,41672229 and 41172206
文摘Geophysical methods have been applied to a wide range of hydrogeological problems. With improvement in geophysical inversion algorithms and measurement tools, significant achievements have been made in the characterization of subsurface architecture, time-lapse monitoring of hydrogeological process and contaminant plumes delineation. In this paper, we summarize the geophysical methods that are most widely used in hydrogeology including Electrical Resistivity Tomography(ERT), Induced Polarization(IP), Ground Penetrating Radar(GPR) and Electromagnetic Induction(EMI). Three examples including lab and field works are used to demonstrate current application of geophysical methods for characterizing subsurface architecture and contaminant plumes. Though great progress has been made in hydrogeohysics over the last few decades at home and abroad, challenges still remain in practical applications. More recently, hydrogeophysics continues to develop in the areas of establishment of hydrogeophysical models, large-scale architecture characterization, uncertainty analysis, biogeochemical process monitoring and ecosystem science.
基金funded by the National Basic Research Program of China(973 Program)(No.2013CB733203)the National Natural Science Foundation of China(No.41474055)
文摘The dynamic monitoring of landslides in engineering geology has focused on the correlation among landslide stability,rainwater infiltration,and subsurface hydrogeology.However,the understanding of this complicated correlation is still poor and inadequate.Thus,in this study,we investigated a typical landslide in southwestern China via time-lapse electrical resistivity tomography(TLERT) in November 2013 and August 2014.We studied landslide mechanisms based on the spatiotemporal characteristics of surface water infiltration and flow within the landslide body.Combined with borehole data,inverted resistivity models accurately defined the interface between Quaternary sediments and bedrock.Preferential flow pathways attributed to fracture zones and fissures were also delineated.In addition,we found that surface water permeates through these pathways into the slipping mass and drains away as fissure water in the fractured bedrock,probably causing the weakly weathered layer to gradually soften and erode,eventually leading to a landslide.Clearly,TLERT dynamic monitoring can provide precursory information of critical sliding and can be used in landslide stability analysis and prediction.
基金supported by the Swiss National Science Foundation
文摘Geophysical techniques can help to bridge the inherent gap that exists with regard to spatial resolution and coverage for classical hydrological methods. This has led to the emergence of a new and rapidly growing research domain generally referred to as hydrogeophysics. Given the differing sensitivities of various geophysical techniques to hydrologically relevant parameters, their inherent trade-off between resolution and range, as well as the notoriously site-specific nature of petrophysical parameter relations, the fundamental usefulness of multi-method surveys for reducing uncertainties in data analysis and interpretation is widely accepted. A major challenge arising from such endeavors is the quantitative integration of the resulting vast and diverse database into a unified model of the probed subsurface region that is consistent with all available measurements. To this end, we present a novel approach toward hydrogeophysical data integration based on a Monte-Carlo-type conditional stochastic simulation method that we consider to be particularly suitable for high-resolution local-scale studies. Monte Carlo techniques are flexible and versatile, allowing for accounting for a wide variety of data and constraints of differing resolution and hardness, and thus have the potential of providing, in a geostatistical sense, realistic models of the pertinent target parameter distributions. Compared to more conventional approaches, such as co-kriging or cluster analysis, our approach provides significant ad- vancements in the way that larger-scale structural information eontained in the hydrogeophysieal data can be accounted for. After outlining the methodological background of our algorithm, we present the results of its application to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the detailed local-scale porosity structure. Our procedure is first tested on pertinent synthetic data and then applied to a field dataset collected at the Boise Hydrogeophysical Research Site. Finally, we compare the performance of our data integration approach to that of more conventional methods with regard to the prediction of flow and transport phenomena in highly heterogeneous media and discuss the implications arising.
文摘This research applies Electrical Resistivity Tomography (ERT) and Time Domain Electromagnetic Method (TDEM) to study the hydrogeology of the Taubaté basin, which is characterized by half-grabens with about 850 m of maximum sediments thickness. The study area is in Taubaté city, São Paulo State, Brazil, where the Taubaté aquifer is an important water source. The Taubaté Group is the main sedimentary package of the basin;it is formed mainly by shales that form aquicludes, and thin layers of sandstones that form the aquifer. There are 40 groundwater exploration wells in Taubaté city that provide important information. The study purpose is to characterize the geoelectrical stratigraphy of the subsurface to locate the contact between the Quaternary and Tertiary sediments and to identify the Taubaté aquifer. The ERT is used for shallow investigations (tens of meters) and the TDEM can reach a great investigation depth (hundreds of meters). Therefore, these geophysical methods are complementary. The ERT data were acquired with the pole-dipole array with 20 m of electrodes spacing and 400 m length, and the TDEM data with the central-loop array with a 200 × 200 m transmitter loop. The results permit to define the contact between the Quaternary and Tertiary sediments around 15 m depth, the Pindamonhangaba Formation between 15 m and 30 m depth and the Taubate Group between 30 m and 300 m depth. The TDEM method defined the Taubaté Group as a single geoelectric layer because the shale and the sandstone layers are all very conductive. The basement is formed by gneiss, which is a very resistive rock. The TDEM method is not able to identify a high conductor/resistor contrast. Overall, the results are consistent with the known geology and the wells information.
文摘This work uses 2D TEM (Transient Electromagnetic) modeling for a hydrogeological study in the Paraná sedimentary basin. The study area is located at the northern region of the state of S?o Paulo, Brazil, where groundwater is exploited from two aquifer systems: one sedimentary, shallow, and the other crystalline, deep. The interest in applying the TEM method in this area owes to the high exploitation rates of groundwater from the crystalline aquifer system for irrigation, which is triggering considerable seismic activity locally. This aquifer system is composed of fractured basalt within the Serra Geral Formation and is about 120 m deep. Eighty-six TEM soundings were acquired at this location, but in nine cases the data did not fit the modelled curve for 1D geoelectrical models due to the geological complexity of the area. This paper shows 2D geoelectrical modeling results based on the FDTD (Finite Differences in Time Domain) method to explain the lateral resistivity variation within the geological setting. A 2D model was generated for each sounding and compared with 1D inversion models as well as with direct information from wells. The results show some vertical variations of about 10 to 30 meters on the upper interface of the basalt layer from Serra Geral Formation. They are located at approximately 60 meters from the center of the soundings. The existence of these 2D structures in the subsurface can be related to the drainage system in the study area. The presence of these structures may indicate a connection between the shallow and deep aquifer systems, acting like a conduit that may contribute to the seismic activity reported.
文摘This study investigates the groundwater aquifer located in Fayuim oasis. In this study, two of the electromagnetic measurement methods have been used in determining the hydrological situation in the Fayoum oasis. The first is airborne electromagnetic (AEM) which, sometimes is referred to as Helicopter electromagnetic (HEM) and the second is ground Time-domain Electromagnetic method (TEM). The subsurface consists of four geoelectrical layers with a rough slope towards the center. The third and the fourth layers in the succession are suggested to be the two-groundwater aquifers. The third layer saturates with fresh water overlying saline water which exists in the bottom of the second one. It is worth mentioning that the depth of the fresh water surface undulates between the surface level in two lakes in the study area and 57 meters below the ground, whereas the thickness of the fresh water aquifer varies from 13 to 36 meters. The depth of the saline water surface undulates between 59 and 81 meters below the ground. In general, airborne electromagnetic surveying has the advantage of fast resistivity mapping with high lateral resolution. Groundbased geophysical surveys are often more accurate, but they are definitely slower than airborne surveys. It depends on targets of interest, time, budget, and manpower available by the method or the combination of methods that will be chosen. A combination of different methods is useful to obtain a detailed understanding of the subsurface resistivity distribution.
文摘Geophysical study and watershed hydrological delineation have been integrated at downstream of Alasra dam site Norh Azraq area to investigate their potential for artificial groundwater recharge. The total surface area of the watershed was found to be about 195 square kilometers. The estimated annual runoff volumes for the Alasra watershed ranged between 1.2 and 1.8 MCM. Moreover, the interpretation of Ten Time Domain Electromagnetic (TDEM) soundings suggested three principal subsurface layers. The top surface layer has an intermediate resistivity (90 - 110 Ohm·m) with a thickness ranging from a few meters to around 50 m. This layer was interpreted as superficial deposits. The second subsurface layer with variably high resistivity values is composed of unsaturated massive basalt layer and probably belongs to Madhala Olivine Phyric Basalt Formation (MOB). The large variations in resistivity could be ascribed to the degree of water saturation (as a result of groundwater recharge from the nearby harvested water dam), or lithological variations (clay content) and/or due to structural control. The third subsurface layer has low resistivity values (·m to 40 Ω·m) and was found at a depth ranging from 120 to 150 m. This layer could represent a saturated basalt layer with high clay contents. The subsurface structures and major faults have been identified. Based on the results of this study, a combination of surface and subsurface artificial groundwater recharge techniques is highly recommended.
文摘Detailed local geological and hydrogeophysical investigations were carried out for the aquifer in Yaoundé, Cameroon to delineate the architecture of different subsurface geological horizons using lithologs and generated vertical electrical sounding (VES) data. An attempt has also been made to estimate aquifer transmissivity from resistivity data. The transmissivity of the uncon?ned aquifer was computed by determining the Dar-Zarrouk parameters (longitudinal unit conductance and transverse unit resistance) and were compared with the actual field transmissivity. The results showed a direct relation between aquifer transmissivity and transverse resistance. The relationship established has therefore, been generalized in the study area in order to evaluate hydraulic conductivity and transmissivity at all the points where geoelectrical measurements have been carried out. This generalization allows one to derive maps of the product Kσ and transmissivity in the study area based on geoelectrical measurements. These maps are important in future modelling processes oriented towards better exploitation of the aquifers.