Coulibaly Yaya, Li Zhen Cheng*, Tikou Belem and Sylvain Lortie
University of Québec in Abitibi-Témiscamingue, Canada
Received Date: October 20, 2016; Accepted Date: December 05, 2016; Published Date: December 08, 2016
Citation: Yaya C, Cheng LZ, Belem T, Lortie S (2016) Auscultation of a Dam for Mining Restoration in Quebec, Canada. J Civil Environ Eng 6:259. doi: 10.4172/2165-784X.1000259
Copyright: © 2016 Yaya C, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Civil & Environmental Engineering
The tailing management is a long term concern for environment and social security. Following the discovery of a new gold mineralization, a mining company aims at reusing one of tailings storage facility at the mining site. However, there is a small treated water reservoir (polishing pond) downstream of the Northwest dam. It is therefore important to know the present state of the Northwest dam’s internal structure. Geophysical methods of electrical resistivity and georadar were used for the auscultation of the Northwest dam. Numerical models were subsequently used to assess the geotechnical behavior of the dam in different deposition situations by simulating the restoration.
The image processing results of georadar data show a layered structure near the surface. Changes in electrical resistivity offer deeper information than georadar. The lateral variability of the electrical resistivity corresponds to heterogeneity within each layer. As we cannot collect samples for the characterization of materials, the geophysical interpretation results help to estimate the compositional structure of the dam; ultimately it helps in numerical modeling on the safety factor estimation.
Dam; Embankment; Electrical resistivity; Factor of safety; Georadar; Internal structure; Tailings pond
The Abitibi greenstone belt of Canada is a world well-known region rich in mineral resources. By-products from exploitation of those natural resources could cause acid mine drainage (AMD) due to the interaction between water discharges, atmospheric oxygen and sulphide minerals. In some cases, metals as nickel (Ni), zinc (Zn), cobalt (Co), copper (Cu) and metalloids as arsenic (As) and antimony (Sb) may generate contaminated neutral mine drainage (DNC). The AMD and the DNC are considered the most important environmental issues, as they may affect the population’s quality of life by contamination of surface water and groundwater. To minimize environmental risks, mining companies have developed some waste management strategies: either 1) make cemented paste for backfilling underground; or 2) by stacking them to the surface with impermeable cover; or 3) dropping in the lake bottom around the operation site (tailings). Tailings dams play an important role to separating contaminated and uncontaminated areas. The stability of dams is thus a key element to be monitoring, especially in intensive mining activities region as the Abitibi greenstone belt. Tailings are a mixture of water and fine material from rocks after concentrated minerals being removed and they are piped into a storage facility (e.g. lakes) belted by dams. However, dams can sometimes yield to different types of hazards and dumping of pulp residues and cause extensive damage and major environmental impacts. Figure 1 shows the breakdown of the dam around tailings from Mount Polley in Canada with over 24 million m3 spilled . Damage with the death of 19 people is also recorded due to the bursting of tailings dam in Mariana, Brazil in 2015 . Figure 2 shows the socioeconomic impact of disruptions in the world with the loss of life and on the environment that was increased from 2000 to 2010. 76% of incidents worldwide were related to the upstream construction method . 51% of incidents (breakdowns and accidents) identified by ICOLD (1996, 2001) are due to slope instability of the dam following excessive stress in the foundation soil, the embankment of the dam, inadequate control of water pressure .
The causes of breaking dams are multiple. Apart from construction problems, poor maintenance or unusual weather as cited by Azam and Li , other vulnerability may also be the cause such as the excess pore pressures accumulation due to rapid rising of the dam. The pore pressure (especially in muddy areas) will reduce the effective stresses and resistance to shear residues [6,7]. The static and seismic liquefaction are also considered a common cause of levees broke particularly those enhanced by the upstream method [8,9]. As the break is a physical process (mechanical, hydraulic), in general, the breakdown for embankments is made according to four classical mechanisms [10-12]: external erosion, internal erosion, external instability and liquefaction (Figure 3).
Runoff of rain water can be the cause of the external erosion. The settlements on the crest cause cracks that promote water infiltration into the dam and this could cause internal erosion or slippage in an area of weakness.
Issue of study
A Canadian mining company discovered a new deposit of gold mineralization around an old mining site in the Abitibi greenstone belt of Québec, Canada. The restoration of a tailings lake to store waste rocks is then planned. According to information gathered at the site, a collapsing zone on the Northwest dam was observed in spring 2004 (Figure 4). Following this incident, the dam was repaired by cement grout injection. The objective of this study is to investigate the internal structure of the dam with two geophysical techniques (resistivity and ground penetrating radar), in order to evaluate its present status.
Internal Structure of the Dam from Geophysical Observation
The principle of the GPR (ground penetrating radar) method is:
electr oma gnetic waves at relative high frequency (40 MHz to 1500 MHz) emitted by a transmitter antenna propagate in the soil, and then reflected while they encounter an interface between materials having different dielectric permittivity. The reflected waves are recorded by another antenna in term of time of arrival and amplitude. In assessing changes in travel time two-way (between the emitter and reflector to the receiver) and the magnitude of the reflection (absolute magnitude) alon g a profile, the characteristics of the reflector can be identified (equations below for non-magnetic materials). The speed of the wave prop a g a tion mainly depends on electromagnetic properties of the material, and these are governed by the water content [13-16].
c is the wave speed in free space (0.3 m/ns), V is the wave speed of material, is the relative dielectric permittivity, and D is estimated depth o f the interface where there is a contrast of the dielectric permittivity.
The G PR method is widely used in civil engineering to detect c r a c ks between concrete and steel [17-29]; in hydrogeology and sedimentology [13,30-32]; in archelogy [33,34]; in agriculture for soil types . Despite of the small depth of investigation of GPR, this nondestructive technique has become more popular in recent years for the evaluation of near surface structures through its simplicity and speed of implementation. In general, the frequency assigned to the antennas affects the depth of investigation and resolution of the radar image. The higher the frequency, the depth of the penetration is smaller but the resolution is higher . The depth of penetration also varies depending on the water content [14,16,36]. Neal  did an excellent review on those complexities in sedimentology. The efficiency of the GPR method to distinguishing geological units depends on the contrast in physical properties such as difference in relative dielectric permittivity and in conductivity, also depends on difference in reflection patterns. Below, Table 1 cites those parameters of several soils and rocks commonly seen in dam environment.
|Materials||Dielectric permittivityεr||Conductivityσ (S/m)||Velocity V (m/ns)|
Table 1: Dielectric permittivity ( εr ), conductivity ( σ ) and the velocity (V ) of materials [14,16].
Field work and data processing: Two georadar systems: RIS MF Hi-Mod (IDS) and PulseEKKO (Sensors and Software) were used to collect data in present study. The RIS MF Hi-Mod uses two blinded antennas of 200 MHz and 600 MHz frequencies. The PulseEKKO has three antennas of 50 MHz, 100 MHz and 200 MHz individually. Figure 5 illustrates the weather condition when those surveys conducted. The GPR data were collected along 24 profiles, but only PulseEKKO data were interpreted for the subsurface structure of the dam because GPS data weren’t successfully recorded by RIS HI-Mod. As to PulseEKKO measurements Differential GPS was used for localizing georadar trace location, the base station of GPS is near the northeast end of the dam.
Using REFLEXW, the data processing flow is as follow:
• Since two GPR surveys were conducted in winter time (March 7 and March 17, 2014), we did the time-zero drift correction.
• Elimination of the low frequency noise due to near-surface reflection using the 1D filter function-Subtract mean (Dewow).
• Apply 2D filter-”Background removal” in order to highlight high-frequency signals looking for diffraction hyperboles that could be generated by cavities.
• AGC-Gain or energy decay to enhance low amplitude signals.
• Apply the air velocity (0.3 m/ns) to remove surface reflectors’ diffraction effect and then use Time-cut to visualize the processing result for the first 10 m of the depth.
Interpretation: The Northwest dam is oriented along the direction NE-SW (Figure 6). Three GPR profiles have been chosen to present reflection images for the northwest border (Line 11), the central (Line 9 ) and the southeast edge (West 04-69) of the dam. A consistent subparallel reflection pattern through those GPR images demonstrates a shallow layered subsurface structure. Obviously, the first parallel r eflectors on the top of GPR image are generated by airwave and ground wave including thin snow layer reflection. The reflector located at 2 m to 3 m of depth (indicated by red arrow) is clearly seen on the profile (Line 11) at the northwest border of the dam, but it is becoming blurred towards the southeast edge of the dam. It seems there is another reflector deeper than 4m (indicated by yellow arrows) on line 11, but difficult to see from other two profiles.
B orehole logs from geotechnical drilling F-95-1 are available f or justifying GPR data interpretation. According to the borehole o bservation, the dam is composed of backfill rocks about 5m in thickness on the top, and they overly on tailings. Figure 6 shows that the GPR method has very limited depth of penetration, if the reflector i ndicated by yellow arrows on line 11 corresponds to the interface b etween the rockfill layer and tailings, the reflection at 2 m to 3 m depth may be due to the water table. If it is true, increasing the water content in the lower part of rockfill layer made GPR signal attenuate. For deep structure features, the electric resistivity method may help to go further.
There is a difference in reflection pattern between June (line 09 and line 11) and March (line West04-69) surveys. This could be explained by a high water content of the dam due to the snow melt in June.
The electric resistivity of water is lower than that of rock, variations in the apparent resistivity measured by the electrical method are often correlated with the change in lithology also with the degree of water saturation in rocks (e.g. porous or fractured zones). This technique is widely used to detect voids [37-41], underground water [42-44], fractural zone recognition [45-57], and mineral exploration [58,59]. The difference in resistivity between the target and its surrounding materials is a key factor to ensure this technique is viable. The resistivity of materials used in the construction of dams varies over a range of more than 5 orders of magnitude (Figure 7). This is a major asset for the characterization of materials by the electrical method . Although, the electrical resistivity can be influenced by the porosity and the degree of water saturation of the medium, the clay content, the type of clay, the medium temperature, dry density, pressure and content ions dissolved in the pore water [43,56,61-64], however each substance has a limited range of variation of the resistivity as Figure 8 shows. This is reason why it is possible to estimate the thickness of the materials based on the distribution of the apparent resistivity. Anyhow, the importance of validating the resistivity data interpretation by geotechnical drilling must be considered.
Field work and data interpretation: A maximum of seventy-two electrodes of the electrical imaging system ‘Syscal R1 Plus Switch 72’ (IRIS instrument) was used to auscultate the dam (Figure 9). The 14th of July 2015, the first pseudo-section was performed at the middle of the dam using both the Wenner-Schlumberger and Wenner configurations. The spacing between 72 electrodes is 5 m. The second pseudo-section is on the shelf between a tailings lake (Figure 10) and the Northwest dam, with 54 electrodes and the Wenner-Schlumberger and Wenner configurations, a spacing of 5 m between electrodes.
The survey of July 17 made three pseudo sections using 72 electrodes; one is at a spa c ing of 5 m, and another at both spacing of 5 m and 2.5 m, with the W enner-Schlumberger and Wenner configurations. The pseudo-section with 5 m spacing is superimposed on that of July 14. They show a consistent result. Since the surveys of July 14 involve the dam also the tailings, we show those results in the interpretation section.
The resistivity data are then processed by the RES2DINV inversion software (GEOTOMO Software). It estimates apparent resistivity values for each rectangular block of the subsurface by fitting the observation based on the least-square principal. The software determines the error criterion, the inverted resistivity model is then modified to reduce the degree of error between the apparent resistivity measured and calculated. The inversion operation is repeated iteratively until the error reaches to tolerant value . As a layered structure was preliminary interpreted from georadar survey on the dam, the default smoothnessconstrained inversion method is used for the electric resistivity data inversion.
By comparing the inversion result of the pseudo-section in the middle of the dam (Figure 11) with that on the shelf, they show consistent structure with the observation in a geotechnical drilling F-95-1 (Golder Associates in 1996). According to geotechnical logs, from surface to the end of drilling the dam consists of rockfill, tailings slits, the silty clays, the silt and moraine layer. An outcrop was observed at the North-east end of the dam. This outcrop corresponds to schistose sandstone; its resistivity value was estimated between 100 Ω • m to 1000 Ω • m by Savoie et al. Our inversion results show that there is a resistivity contrast between different materials; however, it is difficult to distinguish the tailings silts and silty clays from the resistivity. With the help of drilling data, we interpret the distribution of the electrical resistivity in term of changes in materials as cited in Table 2.
|Geotechnicaldriling F-95-1||Electric resistivity|
|Rockfill||0-5 m||73-300 Ω • m (0-5 m)|
|Tailings||5-14 m||18-73 Ω • m (5-27 m)|
|Morain layer||>27 (?)||73-500 Ω • m (27-50 m)|
|Bedrock||?||> 570 Ω • m (> 65 m ?)|
Table 2: Comparison the stratigraphic profile of geotechnical drilling and lithology estimated with the resistivity interpretation.
Taking into account the variation of water content which affects the electrical resistivity it is sometimes difficult to interpret lithologic boundaries according to the 2D model of the resistivity inversion. A greater electrical conductivity could be due to high water content as shown by Maqsud et al. , where low resistivity values are coincide with high values of measured water content. Based on the amended soil and rock resistivity values from Palacky  the internal structure of the dam is interpreted as follow:
• The rockfill at the surface is highly heterogeneous in grain size also in water content. Probably this explains the large range of the resistivity variations from 73 Ω.m up to 300 Ω.m. The thickness of this layer is less than 5 m.
• As the shelf is covered by tailings that are comprised of very fine particles and water-saturated, the inversion result shows this top layer a low resistivity less than 73.5 Ω.m. The comparison between this inversion result with that in the middle of the dam, the low resistivity (<100 Ω.m) layer between 5 m to 25 m is interpreted as the tailings and silty clay layer (Figure 12).
• Since the geotechnical drilling stopped at 27.5 m, from an outcrop next to the Northeast end of the dam we guess that the bedrock may have the same composition as the outcrop. The resistivity of the outcrop is more than 500 Ω.m, the bedrock could be at more than 50 m of the depth (Figure 13). Therefore, the moraine between the tailings/silty clay layer and bedrock is about 20 m in the thickness.
In general, the main objectives associated with the stability analysis of a dam are to: 1) determine the conditions of the rock mass or soil slope stability; 2) to study potential failure mechanisms; 3) determine the sensitivity of slopes to different rupture trigger mechanisms; 4) test and compare the different support options and stabilization; 5) design slopes optimally for security and saving cost. The purpose of our numerical modeling is to determine the safety factor, in a perspective of assess slope stability of the embankment northwest, using SEEP/W and SLOPE/W.
The geometry of the dam
According to the observations in issue, the slope of the Northwest dam varies between 26.6° (upstream side) to 33.7° (downstream side). It reaches 15 meters high. For numerical simulations, it is the geometry of the dam at the end of the construction that will be reproduced. The succession, the thickness and extent of the layers in this initial model are shown on Figure 14. The lower limit at the depth would be the very dense moraine layer overly on the bedrock. These two layers are not included in the digital model of stability analysis because they can be already considered as stable and remain in balance even after the construction of the dam.
Estimation on properties of materials
The parameter values used in SLOPE/W (Table 3) mainly come from a private geotechnical report for the mining company. The value of the specific weigh t of the rockfill is estimated based on previous works at similar situation. For example, the value of 19 kN/m3 from Massiera et al.  and 20 kN / m3 from Ormann et al. . Therefore, the average value of 19.5 kN /m3 is used in the present study.
|Soil type||Weight (kN / m³)||Friction angle (degrees)||Cohésion (kPa)||Undrainedshearstrength (kPa)|
Table 3: Parameters values used in modeling (SLOPE/W).
Determination the level of the water table
Before the simulation with SLOPE/W, we have to determine the level of groundwater that plays an important role in the safety factor. The numerical code SEEP/W has been used repeatedly to find the water level in function of variation of the height of water in the upstream lake. Following saturated hydraulic conductivity values (Table 4) collected from the publications of Chapuis and Aubertin , Mbonimpa et al. , Ormann et al. . The private geotechnical report also shows a vertical hydraulic conductivity estimated in the range of 5 × 10-8 cm/sec for the silty clay and between 5 × 10-3 cm/sec to 5 × 10-6 cm/sec for silts.
|Ormann et al. ||Kx (cm/sec)||10-3||10-8||-||10-9|
|Mamert et al. (2002)||K (cm/sec)||10-8||10-6-10-11|
|Chapuis ||K (cm/sec)||-||-||10-8||10-9|
|Values used in SEEP/W||Kx=Ky (cm/sec)||10-6||10-7||10-8||10-9|
Table 4: Reference of saturated hydraulic conductivity values (ksat).
We estimated the volumetric water content to 0.5 for each of the materials, since we know that the dam is not entirely saturated (above the water table).
Stability analysis with SLOPE/W
By convention, the design the geometry of dams would ensure a minimum safety factor of 1.4 to 1.5 for each phase of construction. The minimum safety factor of 1.5 was also used in the dam stability analysis in Sweden by Ormann et al.  and the same value is recommended by the ACB for the long-term analysis (Table 3). We notice that values of the safety factor prescribed in the “Guide and how to prepare for restoration of mining sites in Quebec” by the Ministry of Energy and Natural Resources of Qu ébec (1997), are also between 1.3 and 1.5. Therefore, we use the same geometry as the initial model for varying the level of water in the upstream lake while filling it with sterile in the simulation.
The boundary conditions are set to zero pressure at the downstream end of the dam in SEEP/W. An average value of 2.10-8 m/sec is used for potential infiltration in the dam. This value was estimated from the rainfall data of Mont-Brun station (Rouyn-Noranda) according to the Institut de la Statistique du Québec. Figure 15 illustrates: a) the water table level only due to the precipitation, without rockfill layer in the downstream lake; b) the water table level only due to the precipitation, there is a rockfill layer in the downstream lake; c) iso-values of hydraulic load in relation to the water table; d) iso-interstitial pressure values.
Once basic parameters are being set up, using Morgenstern- Price analysis method and Mohr-Coulomb failure criterion; defining the specific gravity for each material, the effective cohesion, and the effective angle of internal friction, then we start the analysis process and end with the calculation of a minimum safety factor.
Influence of the water level on the upstream side of the dam: Figure 16 show that when the water level rises in the upstream lake, the minimum safety factor decreases. There is no storage of materials and also considering that the depth of water in the downstream lake is fixed.
When the water level (upstream side) is 11 m to zero elevation the safety factor is 2.15. When the water level reaches almost 13 m, the safety factor decreases to 2 while at 14 m of water level the safety factor is only 1.88.
Effect of filling the upstream lake of the dam: In a more realistic case where the tailings in the upstream lake form a slightly inclined layer, the safety factor increases compared to the case where the lake contains only water. When more fills in the upstream lake with waste rock, assuming there is only water infiltration from rainfall, the safety factor goes from 1.75 to 2.16. If the water level in the downstream lake rises to 8 m and then in stationary regime, there is an even higher safety factor of 2.32 for the dam.
By filling part of the downstream lake with waste rock and the infiltration conditions of potential water either side of the dam, the safety factor i s 2.99. Adding another layer of waste rock over the existing one in downstream side and under the same conditions where there is only the infiltration of water from rainfall on either side of the dam, the safety factor increases significantly to 5.27 as that is expected [69,70].
The two geophysical t echniques have identified the internal structure of the dam which is consistent with the observed lithology in geotechnical drilling F-95-1.
The method of electrical resistivity demonstrates some heterogeneity of materials of th e dam, particularly in the upper layer of riprap by a strong variation of resistivity, probably caused by the variation in particle size of ma terials. This heterogeneity could potentially cause differential settlement of the dam in time. The saturated zone in water as the layer of t a ilings is characterized by low electrical resistivity values (<100 Ω⋅m ). This study demonstrates the great potential of geophysical methods for monitoring for embankments.
By combi n ing the s tudy of the internal structure and stability analyzes, the results will be useful not only for the present condition assessment of the dam; thus they can be considered when planning new developments of the mine.
An analysis of deformation and effective stress would be useful in future work. If the effective stress decreases, the shear strength of the material decreases and this may generate rupture. Wherever, possible, it would be necessary to carry out geotechnical drilling to collect samples in order to determine real geomechanical parameters, untimely to perform numerical simulations of more updated stability analysis.