An evaluation of five petrophysical models for hydraulic permeability predication was completed for a site-specific set of alluvial sediments. The models considered were the Kozeny-Carmen (1927–1937), Börner (1996), Revil-Cathles (1999), and two Revil-Florsch (2010) models. The river deposited sediments represent a relatively narrow grain-size distribution, and were acquired adjacent to the Kansas River, in Olathe, Kansas. Using measured physical, electrical, and hydraulic data from Slater and Glaser (2003) , a comparison of the performance of these five petrophysical models for hydraulic permeability estimation of soils was completed. For models where the key parameter is effective grain-size, three model variations were considered using d10, d50, and d90, resulting in an evaluation of a total of eleven individual models. Parameters included in the models can be classified as three different types: a) physical parameters only, b) electrical parameters only, and c) physical and electrical parameters together. The performance of each model was rated in terms of linear regression R squared, slope, and y-intercept values when plotted against the measured hydraulic conductivities. The top three models were the Kozeny-Carmen, a modified Revil-Cathles, and the Börner model. The Kozeny-Carmen model performed with the highest rating, followed by the modified Revil-Cathles model, and the Börner model rounding out the top three. There was a significant disparity between the rating associated with the top three and the fourth best performing model suggested by Revil and Florsch. However, it should be noted that the Börner model and Revil-Florsch model are based entirely on electrical measurements. The Revil-Cathles model was greatly improved when d10 was substituted for d50 for this limited sediment grain-size range.
{"title":"A Site-specific Comparison of Permeability Prediction Models in Alluvial Sediments from Physical and Geoelectrical Measurements","authors":"D. Glaser","doi":"10.32389/jeeg21-025","DOIUrl":"https://doi.org/10.32389/jeeg21-025","url":null,"abstract":"An evaluation of five petrophysical models for hydraulic permeability predication was completed for a site-specific set of alluvial sediments. The models considered were the Kozeny-Carmen (1927–1937), Börner (1996), Revil-Cathles (1999), and two Revil-Florsch (2010) models. The river deposited sediments represent a relatively narrow grain-size distribution, and were acquired adjacent to the Kansas River, in Olathe, Kansas. Using measured physical, electrical, and hydraulic data from Slater and Glaser (2003) , a comparison of the performance of these five petrophysical models for hydraulic permeability estimation of soils was completed. For models where the key parameter is effective grain-size, three model variations were considered using d10, d50, and d90, resulting in an evaluation of a total of eleven individual models. Parameters included in the models can be classified as three different types: a) physical parameters only, b) electrical parameters only, and c) physical and electrical parameters together. The performance of each model was rated in terms of linear regression R squared, slope, and y-intercept values when plotted against the measured hydraulic conductivities. The top three models were the Kozeny-Carmen, a modified Revil-Cathles, and the Börner model. The Kozeny-Carmen model performed with the highest rating, followed by the modified Revil-Cathles model, and the Börner model rounding out the top three. There was a significant disparity between the rating associated with the top three and the fourth best performing model suggested by Revil and Florsch. However, it should be noted that the Börner model and Revil-Florsch model are based entirely on electrical measurements. The Revil-Cathles model was greatly improved when d10 was substituted for d50 for this limited sediment grain-size range.","PeriodicalId":15748,"journal":{"name":"Journal of Environmental and Engineering Geophysics","volume":"39 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2021-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"73666701","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In underground tunnel constructions, belt conveyors or subway trains can be used as passive sources for time-lapse seismic assessments of surrounding formations. The key technology of passive seismic processes is seismic interferometry (SI), which has the ability to extract the Green's function between any two points in a medium of interest. However, since the process results are typically required to meet certain premise assumptions, the current SI methods cannot be directly used for belt conveyors. In the current study, this type of source is referred to as a near-field thread source and the proposed new SI method is adapted to this type of source. There are mainly two types of SI methods: SI by cross-correlation and SI by multi-dimensional deconvolution. The former requires that the passive source be uniformly distributed in the far field, while the latter relaxes the requirements for uniform distribution yet still requires far field sources. When the first assumptions are violated, the correlation function is proportional to a Green's function with a blurred source. The source blurring is then quantified by what is referred to as an interferometric point-spread function, which can be derived from the observed data. Therefore, using SI by MDD can effectively deblur the source of the Green's function by deconvolving the PSF. In this study, one set of geophones is located in one tunnel as the passive source, and another set is in an adjacent tunnel. However, since the setup almost coincided with the distribution of the passive sources, the spatial (wavenumber) spectrum of the PSF is too broad to be inverted. Therefore, a conventional MDD method is not suitable for this study's experiments. This study found that another PSF could be obtained from the data of the adjacent tunnel. Subsequently, a new MDD formula based on the new PSF is derived and compared with the derivation processes of the conventional MDD formula. Then, the feasibility of the proposed method is verified using numerical simulations.
{"title":"Multi-Dimensional Deconvolution for Near-Field Thread Seismic Sources in Tunnels","authors":"Lu Bin","doi":"10.32389/jeeg20-057","DOIUrl":"https://doi.org/10.32389/jeeg20-057","url":null,"abstract":"In underground tunnel constructions, belt conveyors or subway trains can be used as passive sources for time-lapse seismic assessments of surrounding formations. The key technology of passive seismic processes is seismic interferometry (SI), which has the ability to extract the Green's function between any two points in a medium of interest. However, since the process results are typically required to meet certain premise assumptions, the current SI methods cannot be directly used for belt conveyors. In the current study, this type of source is referred to as a near-field thread source and the proposed new SI method is adapted to this type of source. There are mainly two types of SI methods: SI by cross-correlation and SI by multi-dimensional deconvolution. The former requires that the passive source be uniformly distributed in the far field, while the latter relaxes the requirements for uniform distribution yet still requires far field sources. When the first assumptions are violated, the correlation function is proportional to a Green's function with a blurred source. The source blurring is then quantified by what is referred to as an interferometric point-spread function, which can be derived from the observed data. Therefore, using SI by MDD can effectively deblur the source of the Green's function by deconvolving the PSF. In this study, one set of geophones is located in one tunnel as the passive source, and another set is in an adjacent tunnel. However, since the setup almost coincided with the distribution of the passive sources, the spatial (wavenumber) spectrum of the PSF is too broad to be inverted. Therefore, a conventional MDD method is not suitable for this study's experiments. This study found that another PSF could be obtained from the data of the adjacent tunnel. Subsequently, a new MDD formula based on the new PSF is derived and compared with the derivation processes of the conventional MDD formula. Then, the feasibility of the proposed method is verified using numerical simulations.","PeriodicalId":15748,"journal":{"name":"Journal of Environmental and Engineering Geophysics","volume":"14 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2021-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"75215533","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
R. Adams, B. Miller, W. Kress, S. Ikard, Jason D. Payne, Walter H. Killion
Within the State of Arkansas, there is an increasing number of aging dams and levees that have little to no documentation concerning their construction or composition. Surface geophysical surveys offer a non-intrusive method for investigating these structures to describe their lithologic makeup, evaluate the materials constructed upon, and identify potential flow paths through them. Techniques such as electrical resistivity tomography, seismic refraction, and electromagnetic induction have been used to image dams and levees. They require additional information from geologic outcrops, geotechnical borings, or drill cores to make informed geologic interpretations of the geophysical models. These geologic models then allow the owners of these structures to make more informed decisions about their operation and maintenance. Between 2011 and 2018, the U.S. Geological Survey conducted geophysical and geotechnical investigations of three earthen structures in Arkansas. Electrical and electromagnetic geophysical data were used to develop lithologic models of these structures and characterize the underlying geology. Self-potential surveys were utilized to detect the movement of water through these structures and identify any possible seepage pathways. Geotechnical methods such as electric and hydraulic direct-push well logs and cores acted as a control on the geophysical interpretations and a confirmation of anomalies. This integrated approach detected the lack of an impermeable core within a levee, imaged a change in lithology of the bedrock forming the seal beneath a gravity dam, and identified a potential seepage feature within the core of an earthen dam. These results further support that this method of extending known lithologic features via surface and borehole geophysics is a useful approach for characterizing earthen water-control structures.
{"title":"Evaluation of Electrical and Electromagnetic Geophysical Techniques to Inspect Earthen Dam and Levee Structures in Arkansas","authors":"R. Adams, B. Miller, W. Kress, S. Ikard, Jason D. Payne, Walter H. Killion","doi":"10.32389/jeeg20-063","DOIUrl":"https://doi.org/10.32389/jeeg20-063","url":null,"abstract":"Within the State of Arkansas, there is an increasing number of aging dams and levees that have little to no documentation concerning their construction or composition. Surface geophysical surveys offer a non-intrusive method for investigating these structures to describe their lithologic makeup, evaluate the materials constructed upon, and identify potential flow paths through them. Techniques such as electrical resistivity tomography, seismic refraction, and electromagnetic induction have been used to image dams and levees. They require additional information from geologic outcrops, geotechnical borings, or drill cores to make informed geologic interpretations of the geophysical models. These geologic models then allow the owners of these structures to make more informed decisions about their operation and maintenance. Between 2011 and 2018, the U.S. Geological Survey conducted geophysical and geotechnical investigations of three earthen structures in Arkansas. Electrical and electromagnetic geophysical data were used to develop lithologic models of these structures and characterize the underlying geology. Self-potential surveys were utilized to detect the movement of water through these structures and identify any possible seepage pathways. Geotechnical methods such as electric and hydraulic direct-push well logs and cores acted as a control on the geophysical interpretations and a confirmation of anomalies. This integrated approach detected the lack of an impermeable core within a levee, imaged a change in lithology of the bedrock forming the seal beneath a gravity dam, and identified a potential seepage feature within the core of an earthen dam. These results further support that this method of extending known lithologic features via surface and borehole geophysics is a useful approach for characterizing earthen water-control structures.","PeriodicalId":15748,"journal":{"name":"Journal of Environmental and Engineering Geophysics","volume":"46 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2021-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"79750072","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2021-12-01DOI: 10.32389/1083-1363-26.4.323
{"title":"Best Reviewers and Authors' Biographies","authors":"","doi":"10.32389/1083-1363-26.4.323","DOIUrl":"https://doi.org/10.32389/1083-1363-26.4.323","url":null,"abstract":"","PeriodicalId":15748,"journal":{"name":"Journal of Environmental and Engineering Geophysics","volume":"19 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2021-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"82857987","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
André C.G. Kowalski, C. Mendonça, U. Ofterdinger, H. Rocha
Groundwater transport in crystalline rocks follows pathways along fractured zones because of low primary porosity and permeability in such formations. Fractured systems encompass an imbricated set of joints and fractures with different lengths, apertures and orientations resulting in complex permeable systems with heterogeneous groundwater transport properties. Geophysical well logging has proved effectiveness in detecting depth levels with denser fracture distributions as well as the apparent aperture of fractures contributing to groundwater flow. In many cases, the extension spanned by a fracture network cannot be directly inferred because it may extend beyond the radius of investigation of common well logging probes, thus preventing quantitative estimation of critical length for lateral extension a connected fractured system may have. Here we apply a percolation theory model to estimate the critical length as inferred from the linear density of fracture distribution observed at the borehole wall with an optical imaging probe. Our results are analyzed with electrical well logging data (normal resistivity and single-point resistance) cross borehole slug tests using a set of three boreholes. A critical length of 3.9 m was inferred with a percolation model which revealed consistency with the cross borehole slug tests from two wells situated 10 m and 30 m in the vicinity of the monitored borehole. Our results suggest the utility of inferring critical percolation lengths from fracture parameters obtained using standard well logging imaging techniques with potential applications to evaluate groundwater resources, characterize contaminated sites and provide geotechnical information for works in fractured formations.
{"title":"Fracture Critical Length Estimative Using Percolation Theory and Well Logging Data","authors":"André C.G. Kowalski, C. Mendonça, U. Ofterdinger, H. Rocha","doi":"10.32389/jeeg21-019","DOIUrl":"https://doi.org/10.32389/jeeg21-019","url":null,"abstract":"Groundwater transport in crystalline rocks follows pathways along fractured zones because of low primary porosity and permeability in such formations. Fractured systems encompass an imbricated set of joints and fractures with different lengths, apertures and orientations resulting in complex permeable systems with heterogeneous groundwater transport properties. Geophysical well logging has proved effectiveness in detecting depth levels with denser fracture distributions as well as the apparent aperture of fractures contributing to groundwater flow. In many cases, the extension spanned by a fracture network cannot be directly inferred because it may extend beyond the radius of investigation of common well logging probes, thus preventing quantitative estimation of critical length for lateral extension a connected fractured system may have. Here we apply a percolation theory model to estimate the critical length as inferred from the linear density of fracture distribution observed at the borehole wall with an optical imaging probe. Our results are analyzed with electrical well logging data (normal resistivity and single-point resistance) cross borehole slug tests using a set of three boreholes. A critical length of 3.9 m was inferred with a percolation model which revealed consistency with the cross borehole slug tests from two wells situated 10 m and 30 m in the vicinity of the monitored borehole. Our results suggest the utility of inferring critical percolation lengths from fracture parameters obtained using standard well logging imaging techniques with potential applications to evaluate groundwater resources, characterize contaminated sites and provide geotechnical information for works in fractured formations.","PeriodicalId":15748,"journal":{"name":"Journal of Environmental and Engineering Geophysics","volume":"75 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2021-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87007720","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Surface wave methods have increased in popularity as a means to acquire rapid and accurate shear wave velocity ( VS) profiles for engineering applications. Use of the multichannel analysis of surface waves (MASW) method, in particular, has proliferated due to multiple factors, including the ease with which strong signal-to-noise ratio can be achieved and the increased dispersion resolution offered by multichannel acquisitions. However, typical MASW processing to extract dispersion information assumes the surface waves propagate through a layered model. Errors can arise when significant lateral variability is present in the underlying stratigraphy as encountered in certain geologic settings such as residual deposits. This study investigated the effects of such variability on the dispersion information acquired with MASW. In particular, a spatially-correlated Gaussian random field was used to model the natural fluctuations in stiffness introduced by depositional processes, which differs from the approach in other studies where more specific anomalous features have been explored. Numerical modeling was subsequently performed to simulate surface wave propagation in the representative geotechnical site condition. The recovered surface waves were used to develop a subsurface stiffness profile using a dispersion-based pseudo-2D MASW approach and a tomographic approach using full waveform inversion (FWI). The results demonstrate that considerable natural spatial variability significantly complicates interpretation of dispersion information in two primary ways: (1) uncertainty can arise regarding what the dispersion curve exactly quantifies since none of the underlying VS profiles nor the average VS profile are obtained; and (2) the dispersion images exhibit evidence of false depth-related dispersion information indicative of multiple “fundamental” modes from the superposition of multiple stratigraphic units. The FWI procedure that bypasses extraction of dispersion information was found to better recover the underlying subsurface conditions when compared to the pseudo-2D MASW results at the cost of additional computational efforts.
{"title":"Full Waveform Tomography to Address Challenges with Surface Wave Dispersion Information Caused by Significant Stochastic Variability of Subsurface Stiffness","authors":"J. Coe, Siavash Mahvelati","doi":"10.32389/jeeg21-013","DOIUrl":"https://doi.org/10.32389/jeeg21-013","url":null,"abstract":"Surface wave methods have increased in popularity as a means to acquire rapid and accurate shear wave velocity ( VS) profiles for engineering applications. Use of the multichannel analysis of surface waves (MASW) method, in particular, has proliferated due to multiple factors, including the ease with which strong signal-to-noise ratio can be achieved and the increased dispersion resolution offered by multichannel acquisitions. However, typical MASW processing to extract dispersion information assumes the surface waves propagate through a layered model. Errors can arise when significant lateral variability is present in the underlying stratigraphy as encountered in certain geologic settings such as residual deposits. This study investigated the effects of such variability on the dispersion information acquired with MASW. In particular, a spatially-correlated Gaussian random field was used to model the natural fluctuations in stiffness introduced by depositional processes, which differs from the approach in other studies where more specific anomalous features have been explored. Numerical modeling was subsequently performed to simulate surface wave propagation in the representative geotechnical site condition. The recovered surface waves were used to develop a subsurface stiffness profile using a dispersion-based pseudo-2D MASW approach and a tomographic approach using full waveform inversion (FWI). The results demonstrate that considerable natural spatial variability significantly complicates interpretation of dispersion information in two primary ways: (1) uncertainty can arise regarding what the dispersion curve exactly quantifies since none of the underlying VS profiles nor the average VS profile are obtained; and (2) the dispersion images exhibit evidence of false depth-related dispersion information indicative of multiple “fundamental” modes from the superposition of multiple stratigraphic units. The FWI procedure that bypasses extraction of dispersion information was found to better recover the underlying subsurface conditions when compared to the pseudo-2D MASW results at the cost of additional computational efforts.","PeriodicalId":15748,"journal":{"name":"Journal of Environmental and Engineering Geophysics","volume":"20 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2021-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90459524","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Raytracing is a fast and effective numerical simulation method of the seismic wavefield. It plays an important role in field data acquisition design, wavefield analysis, identification, and tomography. In raytracing, pseudo-bending (PB) is a fast and efficient method, but it is unsuitable for complex media with sudden velocity changes. An improved pseudo-bending raytracing method is presented in this paper, which can be applied to any complex medium. The proposed method first decomposes complex medium into multi-scale velocity components and then applies the pseudo-bending approach to the velocity components of different scales. The numerical simulation of seismic wavefield from models shows that the improved multi-scale pseudo-bending (MSPB) method can be applied to a medium with continuous velocity variation and any complex medium with abrupt velocity change.
{"title":"Multi-Scale Pseudo-Bending Raytracing for Arbitrary Complex Media","authors":"Zhang Huan-lan, Wang Bao-Li","doi":"10.32389/jeeg19-007","DOIUrl":"https://doi.org/10.32389/jeeg19-007","url":null,"abstract":"Raytracing is a fast and effective numerical simulation method of the seismic wavefield. It plays an important role in field data acquisition design, wavefield analysis, identification, and tomography. In raytracing, pseudo-bending (PB) is a fast and efficient method, but it is unsuitable for complex media with sudden velocity changes. An improved pseudo-bending raytracing method is presented in this paper, which can be applied to any complex medium. The proposed method first decomposes complex medium into multi-scale velocity components and then applies the pseudo-bending approach to the velocity components of different scales. The numerical simulation of seismic wavefield from models shows that the improved multi-scale pseudo-bending (MSPB) method can be applied to a medium with continuous velocity variation and any complex medium with abrupt velocity change.","PeriodicalId":15748,"journal":{"name":"Journal of Environmental and Engineering Geophysics","volume":"71 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2021-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"85879602","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Zhengyu Liu, Yongheng Zhang, Xinxin Zhang, Huaihong Wang, L. Nie, Xinji Xu, Ning Wang, Ningbo Li
In recent decades, the DC resistivity method has been applied to geophysical monitoring because of its sensitivity to hydrogeological properties. However, existing inversion algorithms cannot give a reasonable image if fluid migration is sudden and unpredictable. Additionally, systematic or measurement errors can severely interfere with accurate object location. To address these issues, we propose an improved time series inversion method for cross-hole electrical resistivity tomography (cross-hole ERT) based on the Extended Kalman Filter (EKF). Traditional EKF includes two steps to obtain the current model state: prediction and correction. We improved the prediction step by introducing the grey time series prediction method to create a new regular model sequence that can infer the potential trend of underground resistivity changes and provide a prior estimation state for reference during the next moment. To include more current information in the prior estimation state and decrease the non-uniqueness, the prediction model needs to be further updated by the least-squares method. For the correction step, we used single time-step multiple filtering to better deal with the case of sudden and rapid changes. We designed three different numerical tests simulating rapid changes in a fluid to validate the proposed method. The proposed method can capture rapid changes in the groundwater transport rate and direction of the groundwater movement for real-time imaging. Model and field experiments were performed. The inversion results of the model experiment were generally consistent with the results of dye tracing, and the groundwater behavior in the field experiment was consistent with the predicted groundwater evolution process.
{"title":"Time Series Data Inversion and Monitoring Method for Cross-Hole ERT Based on an Improved Extended Kalman Filter","authors":"Zhengyu Liu, Yongheng Zhang, Xinxin Zhang, Huaihong Wang, L. Nie, Xinji Xu, Ning Wang, Ningbo Li","doi":"10.32389/jeeg20-051","DOIUrl":"https://doi.org/10.32389/jeeg20-051","url":null,"abstract":"In recent decades, the DC resistivity method has been applied to geophysical monitoring because of its sensitivity to hydrogeological properties. However, existing inversion algorithms cannot give a reasonable image if fluid migration is sudden and unpredictable. Additionally, systematic or measurement errors can severely interfere with accurate object location. To address these issues, we propose an improved time series inversion method for cross-hole electrical resistivity tomography (cross-hole ERT) based on the Extended Kalman Filter (EKF). Traditional EKF includes two steps to obtain the current model state: prediction and correction. We improved the prediction step by introducing the grey time series prediction method to create a new regular model sequence that can infer the potential trend of underground resistivity changes and provide a prior estimation state for reference during the next moment. To include more current information in the prior estimation state and decrease the non-uniqueness, the prediction model needs to be further updated by the least-squares method. For the correction step, we used single time-step multiple filtering to better deal with the case of sudden and rapid changes. We designed three different numerical tests simulating rapid changes in a fluid to validate the proposed method. The proposed method can capture rapid changes in the groundwater transport rate and direction of the groundwater movement for real-time imaging. Model and field experiments were performed. The inversion results of the model experiment were generally consistent with the results of dye tracing, and the groundwater behavior in the field experiment was consistent with the predicted groundwater evolution process.","PeriodicalId":15748,"journal":{"name":"Journal of Environmental and Engineering Geophysics","volume":"73 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2021-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81012474","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
N. Yurkevich, S. Bortnikova, V. Olenchenko, Tanya Fedorova, Y. Karin, A. Edelev, P. Osipova, O. Saeva
Mine tailings are a very active system in which the processes of oxidation, dissolution, and the re-deposition of substances occur in real-time. Time-lapse electrical resistivity tomography and soil-gas measurements have been used on abandoned mine tailings under a highly continental climate, Western Siberia, Russia. The electrical resistivity tomography method allows the structure of the tailings to be determined, namely, its electrophysical parameters, which are related to the chemical composition and geochemical characteristics of the subsurface substance. The aim of this work is to determine the variations in the geoelectrical zoning of sulfide-bearing mine tailings depending on fluctuations in environmental conditions, i.e., ground and air temperature, in conjunction with volatile compounds of environmental concern emanating from the tailings (SO2, CS2, C2H6S). The hourly observations revealed that the configuration of the geoelectrical section varies during the day. The concentration of gases in the surface air layer varied in accordance with the ambient temperature conditions. In general, the minimum gas concentrations were determined at night, and the increase in gas concentrations began when the temperature increased. The dependence of gas formation on temperature conditions differed during the daytime and nighttime. In warmer hours, gas concentrations are highest. At night, when there was a decrease in the temperature of air and then in the ground temperature, a local increase in the concentration of all measured gases occurred at the maximum temperature difference in the air (14.1 °C), and the ground remained relatively warm (20.8 °C). There is a close relationship between ground temperature, electrical resistivity, and the rate of gas production. Local anomalies with the greatest variation in electrical resistivity are associated with the zones that have the most active gas emanations.
{"title":"Time-Lapse Electrical Resistivity Tomography and Soil-Gas Measurements on Abandoned Mine Tailings Under a Highly Continental Climate, Western Siberia, Russia","authors":"N. Yurkevich, S. Bortnikova, V. Olenchenko, Tanya Fedorova, Y. Karin, A. Edelev, P. Osipova, O. Saeva","doi":"10.32389/jeeg21-004","DOIUrl":"https://doi.org/10.32389/jeeg21-004","url":null,"abstract":"Mine tailings are a very active system in which the processes of oxidation, dissolution, and the re-deposition of substances occur in real-time. Time-lapse electrical resistivity tomography and soil-gas measurements have been used on abandoned mine tailings under a highly continental climate, Western Siberia, Russia. The electrical resistivity tomography method allows the structure of the tailings to be determined, namely, its electrophysical parameters, which are related to the chemical composition and geochemical characteristics of the subsurface substance. The aim of this work is to determine the variations in the geoelectrical zoning of sulfide-bearing mine tailings depending on fluctuations in environmental conditions, i.e., ground and air temperature, in conjunction with volatile compounds of environmental concern emanating from the tailings (SO2, CS2, C2H6S). The hourly observations revealed that the configuration of the geoelectrical section varies during the day. The concentration of gases in the surface air layer varied in accordance with the ambient temperature conditions. In general, the minimum gas concentrations were determined at night, and the increase in gas concentrations began when the temperature increased. The dependence of gas formation on temperature conditions differed during the daytime and nighttime. In warmer hours, gas concentrations are highest. At night, when there was a decrease in the temperature of air and then in the ground temperature, a local increase in the concentration of all measured gases occurred at the maximum temperature difference in the air (14.1 °C), and the ground remained relatively warm (20.8 °C). There is a close relationship between ground temperature, electrical resistivity, and the rate of gas production. Local anomalies with the greatest variation in electrical resistivity are associated with the zones that have the most active gas emanations.","PeriodicalId":15748,"journal":{"name":"Journal of Environmental and Engineering Geophysics","volume":"45 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2021-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"74188248","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Groundwater is a nearly exclusive water resource, specifically for the communities which are part of the Chicago metropolitan area. However, water shortage is predicted for many communities in this region, and demand for locating and delineating groundwater is increasing to fulfill the water supply. Shallow sand and gravel aquifers within the glacial deposits of the area specifically are high volume aquifer and less stressed compare to deeper bedrock aquifer. Yet, these aquifers are poorly understood in terms of their extent and lateral variability. This study applied the shear-wave seismic reflection method to delineate the thickness, lateral extent, and internal variability of these aquifers. We acquired horizontally polarized shear-wave (SH-waves) reflection data along five profiles of a total length of 11 km using the land streamer technology in McHenry County in northern Illinois to delineate sand and gravel aquifers. As shear waves propagate through the rock matrix and less sensitive to the presence of water, information from nearby borings and water wells aided the interpretation of the acquired SH-wave seismic profiles. We delineated multiple sand and gravel units of potential aquifers of different thicknesses and lateral extent along with the seismic profiles. The relatively higher vertical and lateral resolution of the shear-waves reflection method and its insensitivity to water saturation or chemistry made it an ideal method to resolve sand and gravel units of potential aquifers within the complex geological environment if aided by water-well information.
{"title":"Seismic Shear-Wave Characterization of Sand and Gravel Groundwater Aquifers in Northern Illinois","authors":"Zonaed Sazal, A. Ismail, J. Thomason","doi":"10.32389/jeeg21-015","DOIUrl":"https://doi.org/10.32389/jeeg21-015","url":null,"abstract":"Groundwater is a nearly exclusive water resource, specifically for the communities which are part of the Chicago metropolitan area. However, water shortage is predicted for many communities in this region, and demand for locating and delineating groundwater is increasing to fulfill the water supply. Shallow sand and gravel aquifers within the glacial deposits of the area specifically are high volume aquifer and less stressed compare to deeper bedrock aquifer. Yet, these aquifers are poorly understood in terms of their extent and lateral variability. This study applied the shear-wave seismic reflection method to delineate the thickness, lateral extent, and internal variability of these aquifers. We acquired horizontally polarized shear-wave (SH-waves) reflection data along five profiles of a total length of 11 km using the land streamer technology in McHenry County in northern Illinois to delineate sand and gravel aquifers. As shear waves propagate through the rock matrix and less sensitive to the presence of water, information from nearby borings and water wells aided the interpretation of the acquired SH-wave seismic profiles. We delineated multiple sand and gravel units of potential aquifers of different thicknesses and lateral extent along with the seismic profiles. The relatively higher vertical and lateral resolution of the shear-waves reflection method and its insensitivity to water saturation or chemistry made it an ideal method to resolve sand and gravel units of potential aquifers within the complex geological environment if aided by water-well information.","PeriodicalId":15748,"journal":{"name":"Journal of Environmental and Engineering Geophysics","volume":"20 1","pages":""},"PeriodicalIF":1.0,"publicationDate":"2021-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"79860245","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}