ARTICLE

Numerical simulation in a wave tank filled with sand

MAMDOH ALAJMI1 JOSÉ M. CARCIONE2,3 AYMAN N. QADROUH1 JING BA3*
Show Less
1 SAC - KACST, P.O. Box 6086, Riyadh 11442, Saudi Arabia.,
2 Istituto Nazionale di Oceanografia e di Geofisica Sperimentale (OGS), Borgo Grotta Gigante 42c, 34010 Sgonico, Trieste, Italy.,
3 School of Earth Sciences and Engineering, Hohai University, Nanjing, 211100, P.R. China.,
JSE 2020, 29(3), 247–260;
Submitted: 7 January 2019 | Accepted: 30 January 2020 | Published: 1 June 2020
© 2020 by the Authors. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution -Noncommercial 4.0 International License (CC-by the license) ( https://creativecommons.org/licenses/by-nc/4.0/ )
Abstract

Alajmi, M., Carcione, J.M., Qadrouh, A.N. and Ba, J., 2020. Numerical simulation in a wave tank filled with sand. Journal of Seismic Exploration, 29: 247-260. We develop a pseudospectral modeling algorithm for wave propagation in anelastic media with Dirichlet and Neumman boundary conditions. The method also allows to set non-reflecting boundaries. The modeling can be adapted to laboratory experiments, namely the implementation of free-surface, rigid and non-reflecting boundary conditions at the model boundaries, as for instance, a tank to perform physical modeling. The time-domain equations for propagation in a viscoelastic medium are described by the Zener mechanical model, that gives relaxation and creep functions in agreement with experimental results. The algorithm is based on a two-dimensional Chebyshev differential operator for solving the viscoelastic wave equation. The technique allows the implementation of non-periodic boundary conditions at the four boundaries of the numerical mesh, which requires a special treatment of these conditions based on one-dimensional characteristics. In addition, spatial grid adaptation by appropriate one-dimensional coordinate mappings allows a more accurate modeling of complex media, and reduction of the computational cost by controlling the minimum grid spacing. An example is shown, where we compute microseismograms in a tank filled with lossy sand.

Keywords
wave tank
Dirichlet conditions
Neumann conditions
anelasticity
full-wave modeling
References
  1. Al-Shuhail, A.A., Alsaleh, M.H. and Sanuade, O.A., 2018. Analysis of time-depth data insand dunes from the empty quarter desert of Southeastern Saudi Arabia. Arab. J. Sci.Engineer. doi.org/10.1007/s13369-018-3178-2
  2. Barriére, J., Bordes, C., Brito, D., Sénéchal, P. and Perroud, H., 2012. Laboratorymonitoring of P waves in partially saturated sand. Geophys. J. Internat., 191:1152-1170.
  3. Berg, P., If, F., Nielsen, P. and Skovgaard, O., 1994. Analytical reference solutions. In:
  4. Helbig, K. (Ed.), Modeling the Earth for Oil Exploration. Pergamon Press, Oxford:421-427.
  5. Bergamo, P., Bodet, L., Socco, L.V., Mourgues, R. and Tournat, V., 2014. Physicalmodelling of a surface-wave survey over a laterally varying granular medium withproperty contrasts and velocity gradients, Geophys. J. Internat., 197: 233-247.
  6. Carcione, J.M., 1996. A 2-D Chebyshev differential operator for the elastic wave equation,Comput. Meth. Appl. Mech. Engineer., 130: 33-45.
  7. Carcione, J.M., Herman, G. and ten Kroode, F.P.E., 2002. Seismic modeling, Geophysics,67: 1304-1325.
  8. Carcione, J.M., 2014. Wave fields in real media: Wave propagation in anisotropic,anelastic, porous and electromagnetic media. Handbook of Geophysical Exploration,
  9. Vol. 38, 3rd revised, extended ed., Elsevier Science Publishers, Amsterdam.
  10. Carcione, J.M., Bagaini, C., Ba, J. Wang, E. and Vesnaver, A., 2018. Waves at fluid-solidinterfaces: Explicit versus implicit formulation of the boundary condition. Geophys. J.Internat., 215: 37-48,
  11. De Angelis, M., Giannini, R. and Paolacci, F., 2010. Experimental investigation on theseismic response of a steel liquid storage tank equipped with floating roof by shakingtable tests. Earthq. Engineer. Struct. Dyn., 39: 377-396.
  12. Fagin, S.W., 1992. Seismic Modeling of Geological Structures: Applications to
  13. Exploration Problems. Geophys. Develop. Series, Vol. 2. SEG, Tulsa, OK.
  14. Farina, B., Poletto, F. and Carcione, J.M., 2017. Synthetic waveforms of axial motion in aborehole with drill string, J. Acoust. Soc. Am., 141: 828-839.
  15. Kosloff, D, and,Tal-Ezer, H., 1993. A modified Chebyshev pseudospectral method withan o(N*) time step restriction. J. Comput. Phys., 104: 457-469.
  16. Krawczyk, C.M., Buddensiek, M.-L., Oncken, O. and Kukowski, N., 2013. Seismicimaging of sandbox experiments-laboratory hardware setup and first reflectionseismic sections, Solid Earth, 4: 93-104.
  17. Lie, I., 1991. Ocean/bottom acoustic interaction with arbitrary bottom profile.FFI/Report-9 1/7009, NDRE.
  18. Prasad, M. and Meissner, R., 1992. Attenuation mechanisms in sands: Laboratory versustheoretical (Biot) data. Geophysics, 57: 710-719.
  19. Sherlock, D.H., 1999. Seismic Imaging of Sandbox Models. Ph.D. thesis, CurtinUniversity of Technology, Perth.
  20. Sherlock, D.H. and Evans, B.J., 2001. The development of seismic reflection sandboxmodelling, AAPG Bull., 85: 1645-1659.
  21. Smolkin, D., 2011. Laboratory scale seismic analysis of a spatially variable hydrologicalsurface in unconfined, unconsolidated sand. M.Sc.thesis, Louisiana State University,Baton Rouge, LA.
  22. Solymosi, B., Favretto-Cristini, N., Monteiller, V., Komatitsch, D., Cristini, P., Arntsen, B.and Ursin, B., 2018. How to adapt numerical simulation of wave propagation andultrasonic laboratory experiments to be comparable - A case study for a complextopographic model. Geophysics, 83(4): T195-T207.
Share
Back to top
Journal of Seismic Exploration, Electronic ISSN: 0963-0651 Print ISSN: 0963-0651, Published by AccScience Publishing