Numerical simulation using a Hamiltonian particle method for effective elastic properties in cracked media
Junichi Takekawa 1 2 Hitoshi Mikada 1 Tada-nori Goto 11 Department of Civil and Earth Resources Engineering, Kyoto University, Kyoto, 615-8540, Japan.
2 Corresponding author. Email: takekawa@tansa.kumst.kyoto-u.ac.jp
Exploration Geophysics 45(2) 116-124 https://doi.org/10.1071/EG13098
Submitted: 3 December 2013 Accepted: 3 December 2013 Published: 21 January 2014
Abstract
We apply a Hamiltonian particle method, one of the particle methods, to simulate seismic wave propagation in a cracked medium. In the particle method, traction free boundaries can be readily implemented and the spatial resolution can be chosen in an arbitrary manner. Utilisation of the method enables us to simulate seismic wave propagation in a cracked medium and to estimate effective elastic properties derived from the wave phenomena. These features of the particle method bring some advantages of numerical efficiencies (e.g. calculation time, computational memory) and the reduction of time for pre-processing.
We describe first our strategy for the introduction of free surfaces inside a rock mass, i.e. cracks, and to refine the spatial resolution in an efficient way. We then model a 2D cracked medium which contains randomly distributed, randomly oriented, rectilinear, dry and non-intersecting cracks, and simulate the seismic wave propagation of P- and SV-plane waves through the region. We change the crack density in the cracked region and determine the effective velocity in the region. Our results show good agreement with the modified self-consistent theory, one of the effective medium theories. Finally, we investigate the influence of the ratio of crack length to particle spacing on the calculated effective velocities. The effective velocity obtained becomes almost constant when the ratio of crack length to particle spacing is more than ~20. Based on this result, we propose to use more than 20 particles per crack length.
Key words: cracked media, digital rock physics, effective elastic property, particle method.
References
Aoi, S., and Fujiwara, H., 1999, 3D finite-difference method using discontinuous grids: Bulletin of the Seismological Society of America, 89, 918–930Budiansky, B., and O’Connell, R., 1976, Elastic moduli of a cracked solid: International Journal of Solids and Structures, 12, 81–97
| Elastic moduli of a cracked solid:Crossref | GoogleScholarGoogle Scholar |
Dahm, T., and Becker, T. H., 1998, On the elastic and viscous properties of media containing strongly interacting in-plane cracks: Pure and Applied Geophysics, 151, 1–16
| On the elastic and viscous properties of media containing strongly interacting in-plane cracks:Crossref | GoogleScholarGoogle Scholar |
Davis, P. M., and Knopoff, L., 1995, The elastic modulus of media containing strongly interacting antiplane cracks: Journal of Geophysical Research, 100, 18253–18258
| The elastic modulus of media containing strongly interacting antiplane cracks:Crossref | GoogleScholarGoogle Scholar |
Del Valle-Garcia, R., and Sanchez-Sesma, F. J., 2003, Rayleigh waves modeling using an elastic lattice model: Geophysical Research Letters, 30, 12-1–12-4
| Rayleigh waves modeling using an elastic lattice model:Crossref | GoogleScholarGoogle Scholar |
Grechka, V., 2007a, Comparison of the non-interaction and differential schemes in predicting the effective elastic properties of fractured media: International Journal of Fracture, 144, 181–188
| Comparison of the non-interaction and differential schemes in predicting the effective elastic properties of fractured media:Crossref | GoogleScholarGoogle Scholar |
Grechka, V., 2007b, Reply to comment by E. H. Saenger: International Journal of Fracture, 146, 293–294
| Reply to comment by E. H. Saenger:Crossref | GoogleScholarGoogle Scholar |
Grechka, V., and Kachanov, M., 2006, Effective elasticity of rocks with closely spaced and intersecting cracks: Geophysics, 71, D85–D91
| Effective elasticity of rocks with closely spaced and intersecting cracks:Crossref | GoogleScholarGoogle Scholar |
Henyey, F. S., and Pomphrey, N., 1982, Self-consistent elastic moduli of a cracked solid: Geophysical Research Letters, 9, 903–906
| Self-consistent elastic moduli of a cracked solid:Crossref | GoogleScholarGoogle Scholar |
Kondo, M., Suzuki, Y., and Koshizuka, S., 2010, Suppressing local particle oscillations in the Hamiltonian particle method for elasticity: International Journal for Numerical Methods in Engineering, 81, 1514–1528
| Suppressing local particle oscillations in the Hamiltonian particle method for elasticity:Crossref | GoogleScholarGoogle Scholar |
Kruger, O. S., Saenger, E. H., and Shapiro, S. A., 2005, Scattering and diffraction by a single crack: an accuracy analysis of the rotated staggered grid: Geophysical Journal International, 162, 25–31
| Scattering and diffraction by a single crack: an accuracy analysis of the rotated staggered grid:Crossref | GoogleScholarGoogle Scholar |
Mariotti, C., 2007, Lamb’s problem with the lattice model Mka3D: Geophysical Journal International, 171, 857–864
| Lamb’s problem with the lattice model Mka3D:Crossref | GoogleScholarGoogle Scholar |
O’Brien, G. S., and Bean, C. J., 2004, A 3D discrete numerical elastic lattice method for seismic wave propagation in heterogeneous media with topography: Geophysical Research Letters, 31, L14608
| A 3D discrete numerical elastic lattice method for seismic wave propagation in heterogeneous media with topography:Crossref | GoogleScholarGoogle Scholar |
O’Brien, G. S., Bean, C. J., and Tapamo, H., 2009, Dispersion analysis and computational efficiency of elastic lattice methods for seismic wave propagation: Computers & Geosciences, 35, 1768–1775
| Dispersion analysis and computational efficiency of elastic lattice methods for seismic wave propagation:Crossref | GoogleScholarGoogle Scholar |
O’Connell, R., and Budiansky, B., 1974, Seismic velocities in dry and saturated cracked solids: Journal of Geophysical Research, 79, 5412–5426
| Seismic velocities in dry and saturated cracked solids:Crossref | GoogleScholarGoogle Scholar |
Okamoto, K., Mikada, H., Goto, T., and Takekawa, J., 2013, Numerical analysis of the relationship between time-variant coda-Q and the variation in crustal stress: Geophysical Journal International, 195, 575–581
| Numerical analysis of the relationship between time-variant coda-Q and the variation in crustal stress:Crossref | GoogleScholarGoogle Scholar |
Saenger, E. H., and Shapiro, S. A., 2002, Effective velocities in fractured media: a numerical study using the rotated staggered finite-difference grid: Geophysical Prospecting, 50, 183–194
| Effective velocities in fractured media: a numerical study using the rotated staggered finite-difference grid:Crossref | GoogleScholarGoogle Scholar |
Saenger, E. H., Gold, N., and Shapiro, S. A., 2000, Modeling the propagation of elastic waves using a modified finite-difference grid: Wave Motion, 31, 77–92
| Modeling the propagation of elastic waves using a modified finite-difference grid:Crossref | GoogleScholarGoogle Scholar |
Saenger, E. H., Kruger, O. S., and Shapiro, S. A., 2004, Effective elastic properties of randomly fractured soils: 3D numerical experiments: Geophysical Prospecting, 52, 183–195
| Effective elastic properties of randomly fractured soils: 3D numerical experiments:Crossref | GoogleScholarGoogle Scholar |
Saenger, E. H., Uribe, D., Janicke, R., Ruiz, O., and Steeb, H., 2012, Digital material laboratory: wave propagation effects in open-cell aluminium foams: International Journal of Engineering Science, 58, 115–123
| Digital material laboratory: wave propagation effects in open-cell aluminium foams:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BC38XhtVequr3P&md5=b12324b0396adcee13959b3bd6e3b690CAS |
Suzuki, Y., and Koshizuka, S., 2008, A Hamiltonian particle method for non-linear elastodynamics: International Journal for Numerical Methods in Engineering, 74, 1344–1373
| A Hamiltonian particle method for non-linear elastodynamics:Crossref | GoogleScholarGoogle Scholar |
Suzuki, Y., Koshizuka, S., and Oka, Y., 2007, Hamiltonian moving-particle semi-implicit (HMPS) method for incompressible fluid flows: Computer Methods in Applied Mechanics and Engineering, 196, 2876–2894
| Hamiltonian moving-particle semi-implicit (HMPS) method for incompressible fluid flows:Crossref | GoogleScholarGoogle Scholar |
Takekawa, J., Madariaga, R., Mikada, H., and Goto, T., 2012, Numerical simulation of seismic wave propagation produced by earthquake by using a particle method: Geophysical Journal International, 191, 1305–1316
| Numerical simulation of seismic wave propagation produced by earthquake by using a particle method:Crossref | GoogleScholarGoogle Scholar |
Takekawa, J., Mikada, H., and Goto, T., 2013, Accuracy of a particle method for modelling of Rayleigh waves: Butsuri Tansa, 66, 85–95
Toomey, A., and Bean, C. J., 2000, Numerical simulation of seismic waves using a discrete particle scheme: Geophysical Journal International, 141, 595–604
| Numerical simulation of seismic waves using a discrete particle scheme:Crossref | GoogleScholarGoogle Scholar |