Characterization of a Landfill Using
Geophysical Methods
L. A. Konstantaki
*
, R. Ghose, D. Draganov, G. Diaferia (Deltares), T. J. Heimovaara
*
Email: l.a.konstantaki@tudelft.nl
We acquired seismic reflection and electrical resistivity (ER) data in the Wieringermeer landfill in the summer of 2013. We used shear-wave (S-wave) sources and receivers, with a source spacing of 1 m and a receiver spacing of 0.5 m over a line of 23.5 m length. The ER measurements had a 1 m electrode spacing covering a 60 m line; both Wenner and dipole-dipole arrays were used.
The goal of this study is to identify higher-density zones in the landfill. The zones act as a barrier for the leachate flow and thus affect the treatment procedure of the landfill. The landfill is a very heterogeneous subsurface and high-resolution measurements and precise processing steps are required for the identification of these possible stagnant zones.
1. Acquisition
The results shown in Figure 1 and 2 provide a good understanding of the subsurface; however, it is also important for the landfill operators to obtain values of the mechanical properties, especially density. Choudhury and Savoikar (2009) performed a literature study of over 30 published works that provided unit weight values and S-wave-velocity values for landfills and so derived an empirical relationship between unit weight and S-wave velocities specifically for landfills:
𝑉𝑠 = 1 (0.0174 − 0.000978𝛾), where Vs is the S-wave
velocity and 𝛾 the unit weight. Using the relation between unit weight and density (𝜌) 𝜌 = 𝛾/𝑔, where g
is the acceleration due to gravity (9.81 m/s2), we can
calculate the density values of the subsurface, see Figure 3 (the different colors for the values are for visualization purposes).
Figure 4 shows the comparison of the study of Choudhury and Savoikar, and our results. Besides seismic reflection analysis we performed also multichannel analysis of surface waves (MASW). Both results show a good fit with the theory, giving confidence to our interpretation.
3. Quantitative Characterization
Seismics
Electrical
resistivity
The processing of the data for obtaining a subsurface image required special steps to account for the heterogeneous subsurface and the noisy environment. The recorded seismic data were first cleaned of the unwanted arrivals (noise from nearby gas pipe, refraction, direct arrivals and surface waves) by muting, leaving primarily the desired reflections.
Next, we performed velocity analysis iteratively to extract a velocity model of the subsurface. These velocities were then used in a prestack depth migration to better handle the heterogeneity of the subsurface. In Figure 1 the prestack depth-migrated result is depicted showing the main features of the landfill: the landfill bottom (yellow line) and several scatterers that correspond to higher-density areas inside the landfill (red ellipses).
4. Next Steps
Besides S-wave and ER data, compressional-wave (P-wave) reflection data were acquired as well. Using the P- and S- wave data, the Poisson’s ratio can be calculated. In addition, the P/S ratio can provide another indication for the location of possible wet pockets. Using mixing rules, we will investigate the possibility to calculate the density of the subsurface to confirm our results. Finally, seismic interferometry will be applied to our data, which has shown to be promising for such a heterogeneous subsurface (Konstantaki et al., 2013).
Figure 2
The electrical resistivity tomography (ERT) image is the result of a joint inversion of the Wenner and dipole-dipole measurements. Figure 2 shows a comparison between the ERT and the interval S-wave velocities extracted from the reflection seismics. The letters indicate locations where we see features in both results. ERT shows areas where dry or wet pockets are located (higher or lower resistive areas), whereas the velocity image shows the higher-density zones (higher S-wave velocities) which is the reason for leachate accumulation (or not) in this area. For example, at event A, a higher-density zone might not allow the leachate to penetrate, resulting in a wet pocket.
Figure 4
Reflection(R) MASW Choudhury & Savoikar(C&S)
Acknowledgements
2. Subsurface Images
Interval velocity field
Joint inversion ERT Distance (m)
LogResistivity (Ohm*m)
Density field (kg/m^3)
CDP CDP
•
D. Choudhury, P. Savoikar, 2009. Simplified method to characterize municipal solid waste properties under seismic conditions. Waste management, v.29, p. 924-933.•
L.A. Konstantaki, D. Draganov, T. Heimovaara, R. Ghose, 2013. Imaging scatterers in landfills using seismic interferometry. Geophysics, v. 78, no.6, p. EN107-EN116.References
This research is supported by the STW under the project number 11035. D.D is supported by CATO2 and by the Division for Earth and Life Sciences (ALW) with financial aid from the Netherlands Organization for Scientific Research (NWO).
100 200 300 400 0 100 200 300 400 0 0 5 10 15 2 5 12 17 25 2 5 12 17 25 4.5 9.5 14.5 19.5 24.5 27 18 28 38 48 58 68 78 88 98 108 18 28 38 48 58 68 78 88 98 108 Figure 3 70 90 150 300 m/s Distance (m) 0 10 20 Landfill Bottom 5 10 15 20 25 Figure 1 1 2 1.9 1.7 1.5 1.3 0.5 1.8 1.4 1.1 0.9 0.7 1.2 1.6 0.8 0.6 R(R)² = 0.944 R(MASW)² = 0.974 R(C&S)² = 0.942 0 50 100 150 200 250 300 0 5 10 15 20 Vs (m/ s) Unit Weight (kN/m^3)