Computational and Experimental Investigations of Fluid Flow in Rock Materials
View this Special IssueResearch Article  Open Access
Kui Wang, Jie Liu, Ying Zou, Mingjie Zhao, Yue Qiang, "ThreeDimensional WaveField Characteristics of EarthRock Dams with Different Hidden Hazards", Advances in Civil Engineering, vol. 2018, Article ID 8941381, 10 pages, 2018. https://doi.org/10.1155/2018/8941381
ThreeDimensional WaveField Characteristics of EarthRock Dams with Different Hidden Hazards
Abstract
The threedimensional wavefield characteristics of dangerous earthrock dams are the premise and foundation to identify the hidden hazards of the earthrock dams via wave test method. Based on the equivalent elastic wave theory in isotropic mediums, the threedimensional wave field of earthrock dams with different defects was calculated using the finite element numerical simulation method. The characteristic and regularity of the threedimensional wave field was explored. The result of the analysis shows that the size and position of the defects are closely related with the wavefield information received on the dam surface. The elasticity modulus and density of the defects have great effect on the wave field, and alternately, Poisson’s ratio of the defects has little effect on the wave field. The greater the difference between the wave velocity of the hazardous materials and the dam body is, the stronger the energy of the scattered wave will be, which in turn will produce clearer lineups of the scattered wave on the timehistory profile. For a single wave movement signal, the dominant frequency amplitude is exponentially correlated with the longitudinal wave velocity of the hazardous materials. Therefore, in the process of wave detection for earthrock dam hazards, the effect of the elasticity modulus, density, and wave velocity on the test signal should be taken into account. The dominant frequency amplitude can serve as one of the control parameters for interpretation of dam hazards by the analysis of the wave test signal.
1. Introduction
Geophysical exploration techniques is an important tool for studying the Earth’s nonuniform structural and dynamic processes, and seismic tomography is one of the most widely used methods [1, 2]. According to the wave field’s characteristic parameters such as kinetic and kinematics, the range and position of any abnormalities can be obtained. The wavefield characteristics can be utilized for the inversion imaging of a particular medium or anomaly. Seismic wave detection technology has been widely used in geotechnical engineering applications, boasting fast, lossless, and simple operation [3–5]. Scholars have studied the application of seismic wave detection technology for earthrock dams. This method was initially used for the geological exploration of the dam and then was gradually utilized for detection of the hidden hazards [6–8]. Research shows that although the wave signals on the dam surface contain abundant information of the dam’s internal structure, it is difficult to improve its inversion accuracy using wave velocity parameters. Therefore, a key application of wave detection technology has become determining how to extract the information about hidden hazards from the test signal. It is necessary to analyze the wavefield distribution characteristics of the dam in order to effectively extract the information of the hidden hazard from the test signal.
The wave equation numerical simulation method is an effective means of studying the wave propagation law and analyzing the internal wavefield distribution of the structure. We can study the propagation mechanism of seismic waves and the interpretation of complex strata by simulating the wave propagation of the composite media. Therefore, more and more numerical simulation methods are used to study the propagation characteristics of seismic waves in structures. For example, Kim et al. implemented the finite element method to study the elastic wave propagation of a threedimensional solid concrete model under the influence of cracks [9]. Narayan studied the propagation characteristics of Rayleigh waves in layered media by the finite difference method [10]. Bakamjian studied the simulation of the threedimensional seismic wave propagation by the boundary integral equation method [11]. Du simulated the wave field of the anisotropic viscoelastic media via pseudospectrum method [12].
In summary, the interpretation method and wavefield characteristics of the existing nondestructive test data cannot directly serve as the basis for detecting the hidden dam’s hazard. It is necessary to describe the wavefield characteristics of an earthrock dam via wave detection technique.
2. ThreeDimensional WaveField Analysis Model for EarthRock Dams
Zhao et al. [8] have proposed a numerical simulation method for threedimensional wave field of earthrock dams. According to this method, a threedimensional finite element model is established. The model size is shown in Figures 1 and 2, and the different hazardous materials have been considered in the model. First, the threedimensional wavefield characteristics of earthrock dams with/without hazards are compared, then, the effects of hidden hazards are analyzed. The parameters of the hazardous material are shown in Tables 1 and 2.(1)Geometric division: theoretically, the wave phenomenon can be observed in a finite element model, so the unit grid size should be less than the wavelength. Usually, the unit grid size is selected less than half of a wavelength to increase accuracy [7]. The unit grid density can be increased depending on the actual model size to improve the accuracy. The data collection interval is the step length of the finite element analysis. In principle, the time step should set to a value where the wave cannot pass through the unit.
(a)
(b)
(a)
(b)


In this paper, tetrahedron unit SOLID92 with 10 nodes is selected. The grid can be divided freely. The grid size of the dam mass is 4 m, and the grid size of the hazardous mass is 2 m. The analysis time T = 0.8 s, and the time step Δt = 1 ms.(2)Boundary conditions: based on the viscoelastic boundary conditions in ANSYS software, a spring damper is utilized to mitigate the effects of boundaryreflected waves on the wavefield analysis.(3)Centrum: the centrum is excited by a delta impulse function. Its expression iswhere A is the maximum amplitude of the excitation force. f is the dominant frequency of the wave load function, where the dominant frequency magnitude directly represents the width of the frequency band. t_{0} is the main peak time of the wave load, and T is the centrum cycle.
In this paper, = −20 kN (vertical downward direction), f = 25 Hz, t_{0} = 0.02 s, and T = 0.04 s.
3. ThreeDimensional WaveField Characteristics of EarthRock Dams with/without Hazards
Wave propagation can be directly judged on the basis of the displacement field. Thus, the displacement field is used to analyze the characteristic rules of the wave field. The threedimensional wavefield characteristics obtained from numerical simulation are shown in Figure 3.
(a)
(b)
(c)
(d)
(e)
(f)
By analyzing Figure 3, it is clear that a direct wave can produce scattered waves when it interacts with a hazardous mass different from the dam material. The scattered wave spreads from the hazardous site to the surface of the dam mass over time exhibiting significant motion. The direct wave gradually spreads out to the distance until eventual absorption caused by the action of the viscouselastic boundary. At this point, the vibration of the dam particle becomes smooth.
The timehistory curves are collected at various measuring points parallel to the dam axis on the dam crest. Figure 4 shows the timehistory profile that reflects the overall wavefield characteristics.
(a)
(b)
(c)
(d)
(e)
(f)
Based on a contrastive analysis of the timehistory profile (Figure 4), it is thus clear that the wavefield information of the earthrock dam spreads out in the periphery in a very regular manner, and there is no other wave information after the direct wave over time when there is no hazardous mass. When a hazardous mass is present, a scattered wave is produced which spreads to the dam crest once the direct wave reaches the hazardous mass. Therefore, significant wave motion signals occur again on the timehistory profile in the period following the passage of the direct wave, and these signals are the scattered waves arising due to the presence of the hazardous mass. Figure 3 (wave field) and Figure 4 (timehistory profile) indicate that the scattered wave spreads to the dam crest along the surface of the dam slope and leads to significant phase changes in wave motion. The figures show that the X direction has more significant characteristics than the other two directions. Thus, the existence and characteristics of the hazardous mass within the dam mass can be judged depending on the scattered wave signals received by the surface of the dam.
4. Analysis and Discussion for the Effects of Hidden Hazards
4.1. Impact Analysis for Elasticity Modulus and Density
Wave velocity is a basic wave parameter for measuring material properties, and it can be used to measure the material properties of the hazardous mass. In cases with a fixed Poisson’s ratio, we discuss the effects of the elasticity modulus and density on the wavefield characteristics. The model is designed in the same manner as mentioned earlier. The hazard size is R = 4 m. The distance between the central position and dam crest is h = 9 m. The parameters of the hazardous material are shown in Table 1.
4.2. Characteristic Analysis for the TimeHistory Profile
The timehistory curve for all measuring points is used to establish a profile to analyze the full wavefield characteristics (Figure 5).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
The first phase of the scattering signal received in the X direction is between 0.1 s and 0.18 s and the central position in the vicinity of the 8th–10th scattering phase. The second phase of the scattering signal received in the X direction is between 0.17 s and 0.24 s and the central position in the vicinity of the 8th–10th scattering phase. It can be concluded that the time points at which the scattered wave generated by the hazardous mass in the same position received on the same measuring line of the dam crest are basically the same. However, the closer the proximity of the defects to that of the surrounding media, the less significant the scattering result. There is limited information available regarding the lineups for Y direction and Z direction. The time points for all lineups are basically the same. A greater difference in material parameters presents a clearer lineup.
4.3. Characteristic Analysis for Wave Signals
Time domain analysis: to enrich the analysis results, the timehistory data of a single measuring point is subjected to a quantitative contrastive analysis based on the analysis for the full wave field. The timehistory curve for the same measuring point on the measuring line is selected for comparison (Figure 6, No. 14 measuring point is selected). A single wave signal on the dam surface in the case of different hazard characteristics has the following regularities: the vibration amplitude of single measuring points at the same distance from the centrum decreases with the decreased difference in material properties between the defects and the surrounding media. Therefore, the wave amplitude with an initiation point of scattered signals serves as a parameter for judging the material properties of the defects.
(a)
(b)
(c)
Frequency domain analysis: a timefrequency analysis is conducted for the timehistory data of No. 8 measuring point in different models mentioned above. Accordingly, we can obtain the regularity concerning the effect of characteristics of the hazardous material on the wavefield characteristics in the frequency domain (Figure 7).
(a)
(b)
(c)
As the dam and hazard materials exhibit decreasing differences in wave parameters, the dominant frequency amplitude of the wave signal decreases significantly, and the regularities in three directions stay consistent. Obviously, the dominant frequency amplitude is highly sensitive in the wave field. Therefore, it can serve as an important parameter in an inversion study in the actual detection of wave motion, and the expected degree of damage can be estimated on the basis of the dominant frequency amplitude of the wave signal of a single measuring point in the site investigation.
4.4. Impact Analysis for Poisson’s Ratio
Based on the nature of homogeneous media, Poisson’s ratio is another parameter determining the elastic characteristics of a material. Therefore, it is necessary to study the effect of Poisson’s ratio on the wavefield information. The elasticity modulus and density of the hazardous material remain unchanged. The specific material parameters are shown in Table 2. The model size and design are the same as mentioned above.
The timehistory curves for all measuring points on the same measuring line are generated by changing Poisson’s ratio of the hazardous material. The analysis of the wave field of the hazardous earthrock dam shows that the wave field in the X direction is highly sensitive to changes in hazards. Therefore, the timehistory profiles for wave motions in the X direction are shown in Figure 8.
(a)
(b)
(c)
Figure 8 shows that the characteristic of the scattered wave is not significant by changing Poisson’s ratio alone, and there are no significant changes in phase that can be observed in the timehistory profile. A weak scattering result can still be observed when the elasticity modulus and density are changed and the wave velocity of the longitudinal wave reaches approximately 500 m/s. No scattering result can be observed when only Poisson’s ratio is changed and the wave velocity of the longitudinal wave is 484.3 m/s. Therefore, we conclude that Poisson’s ratio of the defects has little effect on the scattered wave field of the earthrock dam. In other words, the scattered wave field is not too sensitive to the changes in Poisson’s ratio of the hazardous material.
The timehistory curve for the same measuring point on the measuring line is selected for comparison (Figure 9). With increasing difference in wave velocity between the defect material and the dam mass material, the wave motion increases slightly. However, there is no significantly abnormal wave motion in the later stage. We conclude that Poisson’s ratio has little effect on the wave field of the threedimensional earthrock dam.
(a)
(b)
(c)
5. Conclusion and Suggestions
A finite element model for the analysis of the threedimensional wave field of the earthrock dam with defects was established for numerical simulation of a wave field. The wavefield characteristics of the earthrock dam with/without defects were analyzed. We concluded the following:(1)When hazards are present in the earthrock dam and the size of the defects are not significant enough to cause diffraction, it will produce a series of scattering phenomena, such as reflection and refraction, when the direct wave meets the defect. The scattered wave spreads out and reaches the surface of the dam mass and then is finally received by the detector. It indicates a significant wavefield disturbance. The timehistory curve still fluctuates after the direct wave is propagated. The lineups of the scattered wave are present in the timehistory profile reflecting the overall wave field. The lineups are an important signal that reflects the characteristics of the scattered wave field. It can be used for judging the presence or absence of hazards.(2)The result of the finite element value calculation in the wavefield analysis shows that the vertical triangular load excitation method will be used in the wave motion model of earthrock dam. The threedimensional wave field of a through leakage passage is highly sensitive to abnormal wave motions in the direction vertical to the dam’s axis.(3)The result of the numerical calculation of the finite element model of the earthrock dam with different hazardous material characteristics shows that a greater difference in wave velocity between the defect’s material and the dam’s material indicates a more intensive scattered wave field and a clearer scattering axis on the timehistory profile. In terms of a single wave motion signal, the amplitude corresponding to the dominant frequency is exponentially correlated with the longitudinal wave velocity of the hazardous material. Therefore, the amplitude corresponding to the dominant frequency in the inverse analysis for wave motion can serve as one of the control parameters for inversion.(4)When the differences in longitudinal wave velocity are the same between hazardous materials of the hazardous earthrock dam, the elasticity modulus and density of hazardous material have a significant effect on the wave field. They are highly sensitive in the wave motion analysis; Poisson’s ratio is less sensitive in the wave motion analysis. It is very difficult to distinguish between defects with different Poisson’s ratios during the process of actual exploration.
Data Availability
All the data in this paper were obtained by numerical simulation, and the numerical simulation software is ANYS. The software was purchased by Chongqing Jiaotong University.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge the financial supports from the Natural Science Foundation Project of China (Grant nos. 51609027 and 51879017), Chongqing Research Program of Basic Research and Frontier Technology (Grant nos. cstc2016jcyjA0016 and cstc2017jcyjBX0066), and Scientific and Technological Research Program of Chongqing Municipal Education Commission (Grant no. KJ1501007).
References
 D. Ma, X. Cai, Q. Li, and H. Duan, “Insitu and numerical investigation of groundwater inrush hazard from grouted karst collapse pillar in longwall mining,” Water (Switzerland), vol. 10, no. 9, p. 1187, 2018. View at: Publisher Site  Google Scholar
 L. De Giorgi and G. Leucci, “Detection of hazardous cavities below a road using combined geophysical methods,” Surveys in Geophysics, vol. 35, no. 4, pp. 1003–1021, 2014. View at: Publisher Site  Google Scholar
 H. A. Hamdan and A. Vafidis, “Joint inversion of 2D resistivity and seismic travel time data to image saltwater intrusion over karstic areas,” Environmental Earth Sciences, vol. 68, no. 7, pp. 1877–1885, 2013. View at: Publisher Site  Google Scholar
 A. M. E. Mohamed, A. S. A. Abu El Ata, F. Abdel Azim et al., “Sitespecific shear wave velocity investigation for geotechnical engineering applications using seismic refraction and 2D multichannel analysis of surface waves,” NRIAG Journal of Astronomy and Geophysics, vol. 2, no. 1, pp. 88–101, 2013. View at: Publisher Site  Google Scholar
 J. M. Crane, J. M. Lorenzo, and J. B. Harris, “A new electrical and mechanically detonatable shear wave source for near surface (0–30 m) seismic acquisition,” Journal of Applied Geophysics, vol. 91, pp. 1–8, 2013. View at: Publisher Site  Google Scholar
 S.Z. Wang, C. Tan, G. Chen et al., “The application of Rayleighwave exploration on detecting hidden dam troubles,” Progress in Geophysics, vol. 20, no. 1, pp. 262–266, 2005. View at: Google Scholar
 M.Jie Zhao, Y. U. Dong, and H.Yan Zhao, “Experimental study on velocity and resistivity combined tomography for diagnosing leakage in earth rockfill dam,” Journal of Hydraulic Engineering, vol. 43, no. 1, pp. 118–126, 2012. View at: Google Scholar
 M.J. Zhao, Y. Zou, and X. Zhang, “Threedimensional numerical simulation and characteristic analysis of wavefield for earthrock dam contain hidden troubles,” Journal of Hydraulic Engineering, vol. 47, no. 5, pp. 599–607, 2016. View at: Google Scholar
 J. Kim, J. Woo, and W. B. Na, “Finite element simulation of twopoint elastic wave excitation method for damage detection in concrete structures,” Russian Journal of Nondestructive Testing, vol. 44, no. 10, pp. 719–726, 2008. View at: Publisher Site  Google Scholar
 J. P. , S. K. Narayan, “Study of effects of focal depth on the characteristics of Rayleigh waves using finite difference method,” Acta Geophysica, vol. 58, no. 4, pp. 624–644, 2010. View at: Publisher Site  Google Scholar
 B. Bakamjian, Boundary Integrals: An Efficient Method for Modeling Seismic Wave Propagation in 3D, University of Tulsa, Tulsa, Oklahoma, 1992.
 Q.Z. Du, “Wavefield forward modelling with the pseudospectral method in viscoelastic and azimuthally anisotropic media,” ACTA Physica Sinica, vol. 53, no. 12, pp. 4428–4434, 2004. View at: Google Scholar
Copyright
Copyright © 2018 Kui Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.