Estimativas do teor de umidade empregando o método GPR: uma avaliação comparativa em experimentos de laboratório e campo

Este trabalho faz uma análise das estimativas de teores de umidade obtidas com o método Ground Penetrating Radar (GPR) comparativamente às determinadas com os métodos Time Domain Reflectometry (TDR) e gravimétrico. Os dados foram obtidos em dois experimentos diferentes: um experimento controlado em laboratório buscando reproduzir um meio homogêneo onde foram obtidas as medidas de umidade com GPR (antenas de 1 GHz) e TDR, e outro experimento de campo onde foram obtidos dados com GPR (antenas de 200 MHz) e de amostras de solos do local. Para a obtenção das estimativas a partir do método GPR foram analisados os eventos relativos à onda de transmissão direta entre as antenas, onda refratada criticamente e onda refletida em interfaces com diferentes propriedades elétricas.O GPR mostrou-se sensível às variações de umidades presentes nos dois experimentos e apresentou boa correlação com os dados obtidos com TDR (REQM de0,007 m³m-3) e das amostras (REQM de 0,039 m³m-3).

GPR para a verificação do nível d'água subterrânea em transição floresta amazônica e cerrado

Um estudo do monitoramento do nível de água foi realizado com medidas diretas e com Radar de Penetração no Solo (GPR) em uma floresta tropical de transição para o cerrado. Três poços de monitoramento do nível de água foram instalados durante 2001/2002 em três locais diferentes: o primeiro em uma área de floresta permanente, outro em área de floresta manejada e outra em uma área de pasto. Os perfis de GPR mostram que o nível do topo do lençol freático aparece como um refletor horizontal forte em março e em agosto de 2002, e como um refletor fraco durante medidas em maio e outubro de 2001 com descontinuidades devido a diversas lentes de solo laterítico no solo. O topo do lençol de água é facilmente confundido com a presença de tais lentes. A umidade do solo teve uma influência nestes sinais da reflexão, mudando a constante dielétrica do solo. A profundidade do topo do lençol da água variou 1,8 m sob a floresta permanente, 0,9 m sob a floresta manejada e 3,7 m sob o pasto.

Comparison of time-lapse GPR data collected under natural and forced infiltration conditions to estimate vadose zone hydraulic parameters

Time-lapse crosshole ground-penetrating radar (GPR) data, collected while infiltration occurs, can provide valuable information regarding the hydraulic properties of the unsaturated zone. In particular, the stochastic inversion of such data provides estimates of parameter uncertainties, which are necessary for hydrological prediction and decision making. Here, we investigate the effect of different infiltration conditions on the stochastic inversion of time-lapse, zero-offset-profile, GPR data. Inversions are performed using a Bayesian Markov-chain-Monte-Carlo methodology. Our results clearly indicate that considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions

Examining the information content of time-lapse crosshole GPR data collected under different infiltration conditions to estimate unsaturated soil hydraulic properties

Time-lapse geophysical data acquired during transient hydrological experiments are being increasingly employed to estimate subsurface hydraulic properties at the field scale. In particular, crosshole ground-penetrating radar (GPR) data, collected while water infiltrates into the subsurface either by natural or artificial means, have been demonstrated in a number of studies to contain valuable information concerning the hydraulic properties of the unsaturated zone. Previous work in this domain has considered a variety of infiltration conditions and different amounts of time-lapse GPR data in the estimation procedure. However, the particular benefits and drawbacks of these different strategies as well as the impact of a variety of key and common assumptions remain unclear. Using a Bayesian Markov-chain-Monte-Carlo stochastic inversion methodology, we examine in this paper the information content of time-lapse zero-offset-profile (ZOP) GPR traveltime data, collected under three different infiltration conditions, for the estimation of van Genuchten-Mualem (VGM) parameters in a layered subsurface medium. Specifically, we systematically analyze synthetic and field GPR data acquired under natural loading and two rates of forced infiltration, and we consider the value of incorporating different amounts of time-lapse measurements into the estimation procedure. Our results confirm that, for all infiltration scenarios considered, the ZOP GPR traveltime data contain important information about subsurface hydraulic properties as a function of depth, with forced infiltration offering the greatest potential for VGM parameter refinement because of the higher stressing of the hydrological system. Considering greater amounts of time-lapse data in the inversion procedure is also found to help refine VGM parameter estimates. Quite importantly, however, inconsistencies observed in the field results point to the strong possibility that posterior uncertainties are being influenced by model structural errors, which in turn underlines the fundamental importance of a systematic analysis of such errors in future related studies.

Fracture imaging within a granitic rock aquifer using multiple-offset single-hole and cross-hole GPR reflection data

The sparsely spaced highly permeable fractures of the granitic rock aquifer at Stang-er-Brune (Brittany, France) form a well-connected fracture network of high permeability but unknown geometry. Previous work based on optical and acoustic logging together with single-hole and cross-hole flowmeter data acquired in 3 neighbouring boreholes (70-100 m deep) has identified the most important permeable fractures crossing the boreholes and their hydraulic connections. To constrain possible flow paths by estimating the geometries of known and previously unknown fractures, we have acquired, processed and interpreted multifold, single- and cross-hole GPR data using 100 and 250 MHz antennas. The GPR data processing scheme consisting of timezero corrections, scaling, bandpass filtering and F-X deconvolution, eigenvector filtering, muting, pre-stack Kirchhoff depth migration and stacking was used to differentiate fluid-filled fracture reflections from source generated noise. The final stacked and pre-stack depth-migrated GPR sections provide high-resolution images of individual fractures (dipping 30-90°) in the surroundings (2-20 m for the 100 MHz antennas; 2-12 m for the 250 MHz antennas) of each borehole in a 2D plane projection that are of superior quality to those obtained from single-offset sections. Most fractures previously identified from hydraulic testing can be correlated to reflections in the single-hole data. Several previously unknown major near vertical fractures have also been identified away from the boreholes.

Structural joint inversion of time-lapse crosshole ERT and GPR traveltime data

Time-lapse geophysical monitoring and inversion are valuable tools in hydrogeology for monitoring changes in the subsurface due to natural and forced (tracer) dynamics. However, the resulting models may suffer from insufficient resolution, which leads to underestimated variability and poor mass recovery. Structural joint inversion using cross-gradient constraints can provide higher-resolution models compared with individual inversions and we present the first application to time-lapse data. The results from a synthetic and field vadose zone water tracer injection experiment show that joint 3-D time-lapse inversion of crosshole electrical resistance tomography (ERT) and ground penetrating radar (GPR) traveltime data significantly improve the imaged characteristics of the point injected plume, such as lateral spreading and center of mass, as well as the overall consistency between models. The joint inversion method appears to work well for cases when one hydrological state variable (in this case moisture content) controls the time-lapse response of both geophysical methods. Citation: Doetsch, J., N. Linde, and A. Binley (2010), Structural joint inversion of time-lapse crosshole ERT and GPR traveltime data, Geophys. Res. Lett., 37, L24404, doi: 10.1029/2010GL045482.

Uso do radar penetrante no solo (GPR) na investigação dos solos dos tabuleiros costeiros no litoral norte do estado da Bahia

Foram estudados nove perfis ao longo de uma toposseqüência sobre os sedimentos do Grupo Barreiras, na Fazenda Rio Negro, município de Entre Rios (BA), utilizando a prospecção eletromagnética por meio do Radar Penetrante no Solo - "Ground-penetrating radar - GPR", objetivando analisar a utilização dessa ferramenta na aquisição de informações sobre as feições que ocorrem no solo, mediante a comparação entre os radargramas obtidos e a descrição pedológica. O equipamento utilizado foi um Geophysical Survey System modelo GPR SR system-2, com antena de 80 MHz. A análise radargramétrica confirmou o aparecimento dos fragipãs e duripãs em profundidade, que ocorrem sempre acompanhados de um processo de transformação dos solos do tipo Latossolo Amarelo e Argissolo Amarelo em Espodossolo. Os padrões de reflexão mostram claramente os domínios dos solos argilosos e dos solos arenosos, com e sem a presença dos horizontes endurecidos.

Uso de sistema GPR (Ground Penetring Radar) na avaliação de atributos de um solo sob plantio de cana-de-açúcar

O GPR utiliza técnica indireta para investigação das estruturas de subsuperfície, sendo caracterizado como método não-invasivo, o qual permite extrair informações ao longo do perfil de solo sem perfurar, sondar ou escavar. Este trabalho teve o objetivo de avaliar a utilização do GPR no estabelecimento de atributos físicos do solo em área de cultivo com cana-de-açúcar, além de estabelecer sua correlação com o método do cilindro volumétrico. A área de estudo está localizada no município de Bom Sucesso - PR, sendo utilizada com cultura da cana há mais de 20 anos. Foram obtidas amostras de solos indeformadas com cilindro volumétrico nas profundidades de 0-0,1, 0,1-0,2, 0,2-0,3, 0,3-0,4 e 0,4-0,5 m, ao longo de três transeptos paralelos, sendo o central referente à linha de plantio (transepto 2) e os outros dois correspondentes às entrelinhas laterais (transeptos 1 e 3). Os valores de densidade e porosidade total, obtidos por meio das amostras físicas, foram inicialmente interpretados para verificar o estado do solo nas linhas e nas entrelinhas de cana, nas cinco profundidades. Posteriormente, foram correlacionados com os valores das constantes dielétricas (K) adquiridos pelo georradar, porém apenas para a camada de 0,4-0,5 m, uma vez que as camadas superficiais até a profundidade de 0,4 m não puderam ser avaliadas devido às limitações técnicas do equipamento, aliadas às condições físicas e químicas do solo, que acabaram por prejudicar as reflexões. A densidade mostrou-se menor nas linhas, quando comparada às entrelinhas, para as camadas de 0-0,1 e 0,2-0,3 m. A porosidade total foi menor nas entrelinhas do que nas linhas na camada de 0-0,1 m. Os valores de K foram maiores nas linhas, onde a umidade e a porosidade total foram maiores, e menores nas entrelinhas, onde as densidades encontravam se próximas ao limite crítico. A correlação de K com a densidade e porosidade total não foi significativa. Pequenas modificações na porosidade do solo em função do aumento da densidade foram suficientes para diminuir o valor de K nas entrelinhas, devido à redução na quantidade de água nos poros. Esse fato pode ser confirmado pela correlação de K com a umidade. Assim, pode-se dizer que o georradar foi extremamente sensível à variação de umidade no solo, apresentando melhores resultados para esta variável do que para a densidade e porosidade total determinada por meio do cilindro volumétrico.

Single-hole GPR reflection imaging of solute transport in a granitic aquifer

Identifying transport pathways in fractured rock is extremely challenging as flow is often organized in a few fractures that occupy a very small portion of the rock volume. We demonstrate that saline tracer experiments combined with single-hole ground penetrating radar (GPR) reflection imaging can be used to monitor saline tracer movement within mm-aperture fractures. A dipole tracer test was performed in a granitic aquifer by injecting a saline solution in a known fracture, while repeatedly acquiring single-hole GPR sections in the pumping borehole located 6 m away. The final depth-migrated difference sections make it possible to identify consistent temporal changes over a 30 m depth interval at locations corresponding to fractures previously imaged in GPR sections acquired under natural flow and tracer-free conditions. The experiment allows determining the dominant flow paths of the injected tracer and the velocity (0.4-0.7 m/min) of the tracer front. Citation: Dorn, C., N. Linde, T. Le Borgne, O. Bour, and L. Baron (2011), Single-hole GPR reflection imaging of solute transport in a granitic aquifer, Geophys. Res. Lett., 38, L08401, doi: 10.1029/2011GL047152.

Geostatistical inversion of seismic and GPR reflection images: what can we actually resolve?

Estimation of the spatial statistics of subsurface velocity heterogeneity from surface-based geophysical reflection survey data is a problem of significant interest in seismic and ground-penetrating radar (GPR) research. A method to effectively address this problem has been recently presented, but our knowledge regarding the resolution of the estimated parameters is still inadequate. Here we examine this issue using an analytical approach that is based on the realistic assumption that the subsurface velocity structure can be characterized as a band-limited scale-invariant medium. Our work importantly confirms recent numerical findings that the inversion of seismic or GPR reflection data for the geostatistical properties of the probed subsurface region is sensitive to the aspect ratio of the velocity heterogeneity and to the decay of its power spectrum, but not to the individual values of the horizontal and vertical correlation lengths.

Joint inversion of crosshole GPR and seismic traveltime data

Joint inversion of crosshole ground-penetrating radar and seismic data can improve model resolution and fidelity of the resultant individual models. Model coupling obtained by minimizing or penalizing some measure of structural dissimilarity between models appears to be the most versatile approach because only weak assumptions about petrophysical relationships are required. Nevertheless, experimental results and petrophysical arguments suggest that when porosity variations are weak in saturated unconsolidated environments, then radar wave speed is approximately linearly related to seismic wave speed. Under such circumstances, model coupling also can be achieved by incorporating cross-covariances in the model regularization. In two case studies, structural similarity is imposed by penalizing models for which the model cross-gradients are nonzero. A first case study demonstrates improvements in model resolution by comparing the resulting models with borehole information, whereas a second case study uses point-spread functions. Although radar seismic wavespeed crossplots are very similar for the two case studies, the models plot in different portions of the graph, suggesting variances in porosity. Both examples display a close, quasilinear relationship between radar seismic wave speed in unconsolidated environments that is described rather well by the corresponding lower Hashin-Shtrikman (HS) bounds. Combining crossplots of the joint inversion models with HS bounds can constrain porosity and pore structure better than individual inversion results can.

Constraining 3-D electrical resistance tomography with GPR reflection data forimproved aquifer characterization

Surface-based ground penetrating radar (GPR) and electrical resistance tomography (ERT) are common tools for aquifer characterization, because both methods provide data that are sensitive to hydrogeologically relevant quantities. To retrieve bulk subsurface properties at high resolution, we suggest incorporating structural information derived from GPR reflection data when inverting surface ERT data. This reduces resolution limitations, which might hinder quantitative interpretations. Surface-based GPR reflection and ERT data have been recorded on an exposed gravel bar within a restored section of a previously channelized river in northeastern Switzerland to characterize an underlying gravel aquifer. The GPR reflection data acquired over an area of 240×40 m map the aquifer's thickness and two internal sub-horizontal regions with different depositional patterns. The interface between these two regions and the boundary of the aquifer with then underlying clay are incorporated in an unstructured ERT mesh. Subsequent inversions are performed without applying smoothness constraints across these boundaries. Inversion models obtained by using these structural constraints contain subtle resistivity variations within the aquifer that are hardly visible in standard inversion models as a result of strong vertical smearing in the latter. In the upper aquifer region, with high GPR coherency and horizontal layering, the resistivity is moderately high (N300 Ωm). We suggest that this region consists of sediments that were rearranged during more than a century of channelized flow. In the lower low coherency region, the GPR image reveals fluvial features (e.g., foresets) and generally more heterogeneous deposits. In this region, the resistivity is lower (~200 Ωm), which we attribute to increased amounts of fines in some of the well-sorted fluvial deposits. We also find elongated conductive anomalies that correspond to the location of river embankments that were removed in 2002.

Inversion of multiple intersecting high-resolution crosshole GPR profiles for hydrological characterization at the Boise Hydrogeophysical Research Site

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. In particular, crosshole ground-penetrating radar (GPR) tomography has shown much promise in hydrology because of its ability to provide highly detailed images of subsurface radar wave velocity, which is strongly linked to soil water content. Here, we develop and demonstrate a procedure for inverting together multiple crosshole GPR data sets in order to characterize the spatial distribution of radar wave velocity below the water table at the Boise Hydrogeophysical Research Site (BHRS) near Boise, Idaho, USA. Specifically, we jointly invert 31 intersecting crosshole GPR profiles to obtain a highly resolved and consistent radar velocity model along the various profile directions. The model is found to be strongly correlated with complementary neutron porosity-log data and is further corroborated by larger-scale structural information at the BHRS. This work is an important prerequisite to using crosshole GPR data together with existing hydrological measurements for improved groundwater flow and contaminant transport modeling.

Resultados preliminares en la caracterización geométrica 3D de un litosoma arenoso de progradación deltaica mediante Georadar (GPR) (Fm. Arenisca de Roda, Eoceno, Cuenca de Graus-Tremp).

A set of GPR profiles have been recorded in order to determine the 3D geometry of a prograding delta-front sandbody (Roda sandstone formation, Eocene, Graus-Tremp basin). Common Mid Points (CMP) also have been recorded to obtain the velocity of the electromagnetic wave in ground. In order to build the topsurface of a 3D prism a set of topographic points have been acquired. Most of the GPR profiles are oriented parallel to the progradation direction (NNE-SSW) and show the expected geometries. The 3D prism has been built from the individual profiles, which shows the three dimensional geometry of the sandy lithosome.