scholarly journals Geostatistical inversion of electromagnetic induction data for landfill modelling

Author(s):  
João Narciso ◽  
Leonardo Azevedo ◽  
Marc Van Meirvenne ◽  
Ellen Van De Vijver

<p>The characterization and monitoring of landfills has become a major concern, not only for assessing the associated environmental impact (e.g., groundwater contamination) but also for evaluating the potential for recovery of secondary resources, in particular for the production of raw materials and energy. For both objectives, it is crucial to have knowledge of the waste composition and the current landfill conditions (e.g. water saturation level). Near-surface geophysical surveys have been proven effective for the non-invasive investigation of landfills, in which different methods have been used depending on the specific survey targets.  Because of its sensitivity to two subsurface physical properties, electrical conductivity (EC) and magnetic susceptibility (MS), frequency-domain electromagnetic (FDEM) induction has been successfully applied to the qualitative characterization of urban and industrial landfills, including mine tailings. Yet, due to the generally complex composition and strongly heterogeneous spatial distribution of waste deposits, reconstructing a reliable landfill model from surface geophysical measurements remains challenging. Geostatistical inversion emerges as powerful tool to improve the landfill modelling from geophysical data, allowing for a more detailed description of the spatial distribution of the properties of interest and the associated uncertainty. Additionally, it provides a flexible framework for integrating data from geophysical surveys and conventional sampling from drilling or trenching.</p><p>In this work, we present a new geostatistical inversion technique able for the simultaneous inversion of FDEM data for EC and MS, which optimize the landfill modelling procedure and is sensitive towards change on the physical properties of interest. This method is based on an iterative procedure where ensembles of subsurface models of EC and MS are generated with stochastic sequential simulation and co-simulation. These simulated models are conditioned locally by existing borehole data for these properties and by a spatial continuity pattern imposed by a variogram model. Synthetic instrument response data, including both the in-phase and quadrature-phase components of the FDEM response, are generated from each model using a forward model connecting the data domain (FDEM data) with the model domain (subsurface physical properties). The misfit between the observed and forward-modelled FDEM data, weighted according to the depth sensitivity of the FDEM response toward changes in EC and MS, is used to drive the generation of a new set of models in the next iteration. We illustrate the inversion procedure with synthetic landfill example data sets which were created based on real data collected at a mine tailing in Portugal and a municipal solid waste landfill in Belgium.</p>

2021 ◽  
Author(s):  
Ruzica Dadic ◽  
Martin Schneebeli ◽  
Henna-Reeta Hannula ◽  
Amy Macfarlane ◽  
Roberta Pirazzini

<p>Snow cover dominates the thermal and optical properties of sea ice and the energy fluxes between the ocean and the atmosphere, yet data on the physical properties of snow and its effects on sea ice are limited. This lack of data leads to two significant problems: 1) significant biases in model representations of the sea ice cover and the processes that drive it, and 2) large uncertainties in how sea ice influences the global energy budget and the coupling of climate feedback. The  MOSAiC research initiative enabled the most extensive data collection of snow and surface scattering layer (SSL) properties over sea ice to date. During leg 5 of the MOSAiC expedition, we collected multi-scale (microscale to 100-m scale) measurements of the surface layer (snow/SSL) over first year ice (FYI) and MYI on a daily basis. The ultimate goal of our measurements is to determine the spatial distribution of physical properties of the surface layer. During leg 5 of the MOSAiC expedition, that surface layer changed from the  surface scattering layer (SSL),   characteristic for the melt season, to an early autumn snow pack. Here,  we will present data showing both a) the physical properties and the spatial distribution of the SSL during the late melt season and b) the transition of the sea ice surface from the SSL to the fresh autumn snowpack. The structural properties of this transition period are poorly documented, and this season is critical  for the initialization of sea ice and snow models. Furthermore, these data are crucial to interpret simultaneous observations of surface energy fluxes, surface optical and remote sensing data (microwave signals in particular), near-surface biochemical activity, and to understand the sea ice  processes that occur as the sea ice transitions from melting to freezing.</p>


Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. B187-B199 ◽  
Author(s):  
Cyril Schamper ◽  
Flemming Jørgensen ◽  
Esben Auken ◽  
Flemming Effersø

A newly developed helicopter transient electromagnetic (TEM) system has the ability to measure very early times within just a few μs after the turn off of the primary current. For such a system, careful calibration and accurate modeling of the electromagnetic (EM) response is critical to get true resistivities of the very shallow geologic layers. We discovered that this leads to resolution of the same level or in some cases even better than what can be obtained from airborne frequency EM systems. This allowed a range of important applications where high and accurate resolution is mandatory, e.g., geotechnical applications such as urban planning, railroad and road investigations, landslides or distribution of raw materials, and assessing aquifer vulnerability. We evaluated the results of a pilot survey covering the Norsminde catchment south of Aarhus, Denmark, where we found that near-surface layers (top 30 m) can be mapped with an accuracy of a few meters in a complicated glacial sedimentary environment. The mapping of the geologic layers was assessed by a detailed analysis in which we developed a general methodology for crosschecking the EM and borehole data. This methodology is general and can easily be adapted to other data types and surveys. After rating the quality of the boreholes based on a list of predefined criteria, we concluded that the EM data matched with about three-quarters of the boreholes located within less than 15 m from the closest EM soundings. The remaining quarter of the boreholes fell into two groups in which half of the boreholes were of very poor quality or had inaccurate coordinates. Only eight of all the boreholes could not be reproduced by the data, and we attributed this to be caused by very strong lateral or vertical geologic variations not resolvable by the TEM technique.


Geophysics ◽  
1990 ◽  
Vol 55 (11) ◽  
pp. 1441-1446 ◽  
Author(s):  
P. N. Shive ◽  
T. Lowry ◽  
D. H. Easley ◽  
L. E. Borgman

A companion paper (this issue) describes a method for producing three‐dimensional simulations of physical properties for different geologic situations. Here we create a simulation for a particular case, which is a near‐surface (<80 ft deep) description of a karst environment. We simulate seismic velocity, density, resistivity, and the dielectric constant for this situation. We then conduct (in the computer) hypothetical geophysical surveys at the surface of the model. These surveys are seismic refraction, microgravity, dc resistivity, and ground‐probing radar. Physical properties appropriate for cavities are then entered in the model. Repeating the geophysical surveys over the model with cavities provides a convenient method of evaluating their potential for cavity detection. Anomalies produced by normal variations in physical properties may simulate or obscure anomalies from target features. More data about the correlation of physical properties, particularly in the horizontal directions, will be required to evaluate this problem properly.


2021 ◽  
Author(s):  
Leonardo Azevedo ◽  
João Narciso ◽  
Ellen Van De Vijver

&lt;p&gt;The near surface is a complex and often highly heterogeneous system as its current status results from interacting processes of both natural and anthropogenic origin. Effective sustainable management and land use planning, especially in urban environments, demands high-resolution subsurface property models enabling to capture small-scale processes of interest. The modelling methods based only on discrete direct observations from conventional invasive sampling techniques have limitations with respect to capturing the spatial variability of these systems. Near-surface geophysical surveys are emerging as powerful techniques to provide indirect measurements of subsurface properties. Their integration with direct observations has the potential for better predicting the spatial distribution of the subsurface physical properties of interest and capture the heterogeneities of the near-surface systems.&lt;/p&gt;&lt;p&gt;Within the most common geophysical techniques, frequency-domain electromagnetic (FDEM) induction methods have demonstrated their potential and efficiency to characterize heterogeneous deposits due to their simultaneous sensitivity to electrical conductivity (EC) and magnetic susceptibility (MS). The inverse modelling of FDEM data based on geostatistical techniques allows to go beyond conventional analyses of FDEM data. This geostatistical FDEM inversion method uses stochastic sequential simulation and co-simulation to perturbate the model parameter space and the corresponding FDEM forward model solutions, including both the synthetic FDEM responses and their sensitivity to changes on the physical properties of interest. A stochastic optimization driven by the misfit between true and synthetic FDEM data is applied to iterative towards a final subsurface model. This method not only improve the confidence of the obtained EC and MS inverted models but also allows to quantify the uncertainty related to them. Furthermore, taking into account spatial correlations enables more accurate prediction of the spatial distribution of subsurface properties and a more realistic reconstruction of small-scale spatial variations, even when considering highly heterogeneous near surface systems. Moreover, a main advantage of this iterative geostatistical FDEM inversion method is its ability to flexibly integrate data with different resolution in the same framework.&lt;/p&gt;&lt;p&gt;In this work, we apply this iterative geostatistical FDEM inversion technique, which has already been successfully demonstrated for one- and two-dimensional applications, to invert a real case FDEM data set in three dimensions. The FDEM survey data set was collected on a site located near Knowlton (Dorset, UK), which is geologically characterized by Cretaceous chalk overlain by Quaternary siliciclastic sand deposits. The subsurface at the site is known to contain several archaeological features, which produces strong local in-phase anomalies in the FDEM survey data. We discuss the particular challenges involved in the three-dimensional application of the inversion method to a real case data set and compare our results against previously obtained ones for one- and two-dimensional approximations.&lt;/p&gt;


2013 ◽  
Vol 2013 ◽  
pp. 1-7 ◽  
Author(s):  
Zhen Yu ◽  
Jing He

In order to predict the fluid in thin layer precisely, this paper proposed a high-resolution method for fluid prediction. The method used geostatistical inversion with lithology masks to calculate water saturation. We applied this method to theoretical model and real data. The result was compared with that of prestack AVA simultaneous inversion for fluid prediction. It showed that this method had high resolution both in vertical and lateral directions for fluid prediction and could also predict the fluid in thin layer efficiently.


Author(s):  
H. Guliyev ◽  
Kh. Aghayev ◽  
F. Mehraliyev ◽  
E. Ahmadova

In case when the upper part of the medium has complex geological structure and geodynamic processes occur in it, the necessity of these data increases in projecting of the object under construction. Purpose. Studying of acoustic, elastic and anisotropic properties of the upper part of section of complicatedly constructed geological media. Methodology. Seismic observations are conducted in shallow wells in the areas of construction objects located in various seismogeological conditions by NSCW (Near-Surface Cross Well testing) method. Field seismic records are processed. Kinematic and dynamic parameters of pressure and differently polarized shear waves are determined. Thin-layered one-dimensional models of physical properties of the medium are created and interpreted on the basis of nonlinear theory of elastodynamics. Results. It is determined that the medium with high porous, water saturated rocks and anomalous high reservoir pressure has anomalous low value of velocities and gradient of their increase with depth. When this medium was re-examined after deep piles were built there, the overestimated seismic velocities are obtained, which is explained by a decrease in the section of anomalously high reservoir pressure and, accordingly, the porosity of the rocks after piles were built. When the hollowness is increased in unsaturated pebble rocks, the negative value of Poisson's ratio is obtained on the standard method. Seismic anisotropy related with the direction of the grains packing of the rocks is revealed on velocities of shear waves. The change of property of rocks on depth is manifested clearer on frequencies of waves than on their amplitudes. Scientific novelty. The elasticity moduli of the 3rd order are determined which are more sensible to variability of nonlinear elastic properties of rocks of the medium than the moduli of the 2nd order. The values of Poisson's ratio are recalculated for one and the same rocks located in different conditions of rock pressure on the basis of nonclassical theory of deformation. Practical importance. The obtained results can be applied to study the media characterized by complex seismogeological hydrodynamic conditions with clay-sandy rocks of high porosity and water saturation.


2020 ◽  
Author(s):  
Miguel Cooper ◽  
Renata Cristina Bovi ◽  
Cesar Augusto Moreira ◽  
Raquel Stucchi Boschi ◽  
Lucas Moreira Furlan ◽  
...  

&lt;p&gt;Piping is a type of subsurface erosion caused by subsurface water and is considered one of the most difficult erosive processes to study. The nature of this erosion process makes it very difficult to study and quantify.&amp;#160; The aim of this study was to characterize the surface and subsurface distribution of the pipes and to understand the network architecture of pipe systems in tropical forested areas.The study area is situated at the Experimental Station of Tupi, state of S&amp;#227;o Paulo, Brazil. We conducted a Digital Elevation Model allied to a superficial pipe mapping, and 2D and 3D geophysical surveys. The subsurface erosion identified by surface mapping and geophysical surveys appeared at two depths: one more superficial, in the upper part of the study area, and one at greater depth, in the lower part of the study area. The higher topographical positions presented the pipes at less developed stages (closed depressions and simple sinkholes), while the lower topographical positions showed the most advanced features (multiple sinkholes and blind gullies). The method of electroresistivity showed zones where low resistivity values correspond to water saturation (~ 70 omh m) and high values (&gt; 4040 ohm m) that define the pipe; this method was efficient in detecting the presence of collapsed and non-collapsed pipes. We concluded that the use of different methods (superficial and subsuperficial) was essential for the characterization of pipe systems. The integrated analysis of the results obtained from the superficial and 2D subsurface mapping allowed us to infer the spatial continuity of the pipes. The 3D geophysical survey was efficient in mapping soil pipe and the connectivity in situ. The 3D modeling of the pipes revealed the connection and connectivity of the pipe network&amp;#8217;s complexity and morphology.&lt;/p&gt;


Alloy Digest ◽  
2013 ◽  
Vol 62 (9) ◽  

Abstract Böhler (or Boehler) W403 VMR is a tool steel with outstanding properties, based not only on a modified chemical composition, but on the selection of highly clean raw materials for melting, remelting under vacuum (VMF), optimized diffusion annealing, and a special heat treatment. This datasheet provides information on composition, physical properties, and elasticity. It also includes information on forming and heat treating. Filing Code: TS-721. Producer or source: Böhler Edelstahl GmbH.


Geophysics ◽  
2001 ◽  
Vol 66 (1) ◽  
pp. 78-89 ◽  
Author(s):  
Donat Demanet ◽  
François Renardy ◽  
Kris Vanneste ◽  
Denis Jongmans ◽  
Thierry Camelbeeck ◽  
...  

As part of a paleoseismological investigation along the Bree fault scarp (western border of the Roer Graben), various geophysical methods [electrical profiling, electromagnetic (EM) profiling, refraction seismic tests, electrical tomography, ground‐penetrating radar (GPR), and high‐resolution reflection seismic profiles] were used to locate and image an active fault zone in a depth range between a few decimeters to a few tens of meters. These geophysical investigations, in parallel with geomorphological and geological analyses, helped in the decision to locate trench excavations exposing the fault surfaces. The results could then be checked with the observations in four trenches excavated across the scarp. Geophysical methods pointed out anomalies at all sites of the fault position. The contrast of physical properties (electrical resistivity and permittivity, seismic velocity) observed between the two fault blocks is a result of a differences in the lithology of the juxtaposed soil layers and of a change in the water table depth across the fault. Extremely fast techniques like electrical and EM profiling or seismic refraction profiles localized the fault position within an accuracy of a few meters. In a second step, more detailed methods (electrical tomography and GPR) more precisely imaged the fault zone and revealed some structures that were observed in the trenches. Finally, one high‐resolution reflection seismic profile imaged the displacement of the fault at depths as large as 120 m and filled the gap between classical seismic reflection profiles and the shallow geophysical techniques. Like all geophysical surveys, the quality of the data is strongly dependent on the geologic environment and on the contrast of the physical properties between the juxtaposed formations. The combined use of various geophysical techniques is thus recommended for fault mapping, particularly for a preliminary investigation when the geological context is poorly defined.


Sign in / Sign up

Export Citation Format

Share Document