Cellular automata for seismic travel time estimation

• Autorzy: Szostek K. , Dwornik M.
• Języki publikacji: EN

Abstrakty

EN

In this work, the application of the cellular automata for the seismic first breaks time estimation is presented. Cellular automata (CA) algorithms are usually described by set of simple rules applied onto the grid of cells, which can represent one of the few different discrete states. In spite of this simplicity, these algorithms can still simulate variety of complex physical processes, i.e. model the fluids or gases behavior, but also simulate and predict forests fire or the spread of diseases (Turcotte 1997). What is more, these algorithms are easy to accelerate using one of the most advanced parallelization techniques, like graphics card programming - the calculation can be then performed even 1000 times faster. In this work one of the popular CA, presented by Hardy et al. (1973), was modified, according to Rothman (1987), to perform the simulation of the seismic wave propagation in geological medium. This CA was initially designed to simulate gas behavior in a reservoir, but it was modified in the following way: each cell of the grid contains set of values, which represents particles - positive and negative - that describes expansion and compression, respectively, of that cell. Then, both particle sets are calculated separately in two steps: the first is the propagation, when each particle is moved to next cell according to its movement; the second is collision, when rules are applied on each particle, if there is more than one particle in the cell. The collision rules can be divided into categories that represent their complexity, i.e. more complex behavior can be achieved if static particles are implemented, but this requires additional rules of collisions with moving particles. All the rules have to preserve the principle of conservation of mass and momentum. The major modification of the algorithm described by Rothman (1987) is the implementation of the Boltzmann lattice method (Huang 2007), which changes the regular rectangular cellular grid for the triangular network. This modification provides more option when implementing algorithm, especially the collisions rules. On the one hand, triangular grid requires slightly different approach to represent it in computer memory, just like the implementations of the rules. On the other hand, triangular grid can produce better simulation results compared to the classical rectangular grid. The results of this simulation can be used in further calculations, i.e. to solve the inverse problem. To make this CA useful, it has to provide either better than similar methods errorless results or significantly shorter computation time with comparable errors. In this work the CA is compared with mathematical wave propagation solution and the Shortest Path Method (SPM). SPM uses graph theory to reconstruct seismic ray trajectory in a real, geological medium. This method was first described in Moser (1991). As it is more time-consuming than linear method, for large model the parallel approach is necessary (Pięta et al. 2013). One of the major advantages is the possibility of optimization of the modeling quality and calculation time by adjusting number of accessory nodes. This method is easily extendable to cover the anisotropy case. This work describes, verifies and validates the CA algorithms to checks the potential of use in further calculations for solving the forward problem in seismic modeling. The research shows that the presented CA method, in spite of its speed and parallelization options, can not be applied if high precision is required.

Słowa kluczowe

EN

CA algorithm; Shortest Path Method; SPM

• Wydawca: Wydawnictwa AGH
• Czasopismo: Geology, Geophysics and Environment
• Rocznik: 2012
• Tom: Vol. 38, no. 4
• Strony: 541--542
• Opis fizyczny: Bibliogr. 6 poz.

Twórcy

autor

Szostek K.

• AGH University of Science and Technology, Faculty of Geology, Geophysics and Environmental Protection, Department of Geoinformatics and Applied Computer Science; al. Mickiewicza 30, 30-059 Krakow, Poland

autor

Dwornik M.

• AGH University of Science and Technology, Faculty of Geology, Geophysics and Environmental Protection, Department of Geoinformatics and Applied Computer Science; al. Mickiewicza 30, 30-059 Krakow, Poland

Bibliografia

• 1. Hardy J., de Pazzis O. & Pomeau Y., 1973. Time evolution of a two-dimensional model system. I. Invariant states and time correlation functions. Journal of Mathematical Physics, 14, 1746-1579.
• 2. Huang L., 2007. A Lattice Boltzmann Approach to Acoustic-Wave Propagation. Advances in Geophysics, 48, 517—559.
• 3. Moser T.J., 1991. Shortest path calculation of seismic rays. Geophysics, 56, 59-67.
• 4. Pieta A. & Dwornik M., 2013. Parallel implementation of ray tracing procedure in anisotropic medium. TASK Quarterly, 16, 135-143.
• 5. Rothman D.H., 1987. Modeling Seismic P-waves with Cellular Automata. Geophysical Research Letters, 14, 1, 17-20.
• 6. Turcotte D.L., 1997. Fractals and Chaos in Geology and Geophysics. 2n edition. Cambridge University Press, Cambridge.