Телефон: 522-81-45

Geoinformatika 2016; 4(60) : 29-38 (in Russian)  

3D FINITE-DIFFERENCE MIGRATION WITH PARALLELIZATION OF COMPUTING PROCESS

O.O. Verpakhovska, N.V. Pylypenko, O.V. Pylypenko

1Institute of Geophysics, National Academy of Science of Ukraine, 32, Palladin Ave., Kiev, 03680, Ukraine, e-mail: alversim@gmail.com, vpylypenko@gmail.com, lenasloboda@gmail.com

The purpose of the paper is to develop a three-dimensional post stack finite-difference migration with parallelization of computing process to process 3D seismic observations on a separate computer, taking into account the ever-increasing requirements for detailed results and the complexity of the structure of the study area. The effectiveness of the developed version of the post stack finite-difference migration is proved with practical examples, in particular, for the structure of medium complicated by the presence of the salt body.
Design/methodology/approach. Algorithm of a three-dimensional post stack migration is based on reverse full-wave field continuation in the medium that is carried out by solving the differential wave equation using a finite difference method. In this case the differential wave equation is approximated by the difference one  in the four-dimensional space-time grid. In this approach, the coordinate system is converted according to the seismic wave propagation from the depths to the surface, ensuring optimal finite-difference wave field continuation. Parallelization of computations process is on time slices which are defined with a queue. The calculation of the previous slice is stored in the intermediate memory cube. In parallelizing we used library pthread, which is present in all modern versions of the Linux system and allows you to control the execution of the task in the separation calculations on several processors. The software was developed according to the proposed algorithm of a three-dimensional post stack finite-difference migration with parallelization of the computing process for processing.
Findings. The developed three-dimensional post stack full-wave finite-difference migration with parallelization of computing process for processing has been tested in a model and real mode. The obtained results confirm the applicability of the developed method of migration in studding the geological environment of deep structure with different degrees of complexity. To demonstrate the effectiveness of this method we compared the results obtained with those obtained using standard software migration system ProMAX and the developed program with parallelization of computing process on the cluster, developed at the Institute of Geophysics. The comparison of the results allows us to emphasize a general similarity of the results and the difference in some of the particular parts of the image (larger extension of reflecting horizons in the direction of the salt body as a result of applying the proposed variant of post stack full-wave finite-difference migration), which may be attributed to a more accurate migration transformation in the developed program.
Practical value/implications. The developed version of a three-dimensional post stack full-wave finite-difference migration with parallelization of computing process allows on to carry out high-quality processing of large volumes of spatial seismic data in a short time, which depend on the number of processors in the cluster. The results of the testing of developed programs as in model and real data areal seismic survey, prove accuracy and efficiency even in the structure of the medium complicated by the presence of salt body, which is important for oil and gas companies, exploration activities.

Keywords: seismic survey, seismic wave propagation, wave field, modeling, full-wave field continuation, finite-difference method, three-dimensional post stack migration, space-time grid, computer technology, algorithm, cluster, data processing, parallelization of computing process.

The full text of papers

References:

  1. Verpakhovska A.O., Sydorenko G.D., Pylypenko V.N., Pylypenko E.V. AFinite-Deference Method of Pre-Stack Migration Procedures and Modeling of the Wave Field with Parallelizing of Calculation on Cluster. Geoinformatika, 2013, no. 3, pp. 47-58 (in Russian).
  2. Verpakhovska A.O., Pylypenko V.N. 3-D finite-difference modeling of wave field with parallelization of computing process. Sovremennye metody seysmorazvedki pri poiskakh mestorozhdeniy nefti i gaza v usloviyakh slozhnopostroennykh struktur (Seysmo-2011), 2 mezhdunarodnaya nauchno-prakticheskaya konferentsiya, 18-24 sentyabrya 2011, Ukraina, AR Krym, Feodosiya, Kurortnoe: tezisy dokladov. CD-ROM, Kiev, 2011, 5 p. (in Russian).
  3. Godunov S.K., Ryaben’kiy V.S. Raznostnye skhemy. Moscow: Nauka, 1977, 440 p. (in Russian)
  4. Pylypenko V.N., Verpakhovska A.O., Pylypenko O.V. The migration of common shot seismograms in the time domain using finite-difference method. Sovremennye metody seysmorazvedki pri poiskakh mestorozhdeniy nefti i gaza v usloviyakh slozhnopostroennykh struktur (Seysmo-2011), 2 mezhdunarodnaya nauchno-prakticheskaya konferentsiya, 18-24 sentyabrya 2011, Ukraina, AR Krym, Feodosiya, Kurortnoe: tezisy dokladov. CD-ROM. Kiev, 2011, 5 p. (in Russian).
  5. Pylypenko V.N., Verpakhovska A.O., Budkevych V.B., Pylypenko O.V. The formation of three-dimensional image of the medium on the CDP sum to study the geological structure of the mine fields. Geofizicheskiy zhurnal, 2015, vol. 37, no. 4, pp. 104-113 (in Russian).
  6. Wang Y., Li C., Tian Y., Yan H., Zhao C., Zhang J.A Parallel Algorithm of Kirchhoff Pre-stack Depth Migration Based on GPU. Algorithms and Architectures for Parallel Processing, 14th International Conference. Proceedings. Part II, august, 24-27, 2014. Dalian, China, 2014, pp. 207-218.
  7. Claebout J.F., Doherty S.M. Downward continuation of moveout-corrected seismograms. Geophysics, 1972, vol. 37, iss. 5, pp. 741-768.
  8. Shi X., Li C., Wang S., Wang H. Computing prestack Kirchhoff time migration on general purpose GPU. Computers & Geosciences,2011, vol. 37, iss. 10, pp.1702-1710.
  9. Huang T., Li X., Zhang T., Lu D. GPU-accelerated Direct Sampling method for multiple-point statistical simulation. Computers & Geosciences, 2013,57, рp. 13-23.
  10. Li Z. Wave-field extrapolation by the linearly transformed wave equation. Geophysics, 1986, vol. 51, iss. 8, pp. 1538-1551.
  11. Chakraborty S., Yerneni S., Phadke S., Bhardwaj D. Parallelization Strategies for Seismic Modelling Algorithms. Journal of Indian Geophysical Union, 2003, vol. 7, no.1, pp. 11-14.
  12. Geist A., Beguelin A., Dongarra J., Jiang W., Manchek R., Sunderam V. Parallel Virtual Machine: A Users’ Guide and Tutorial for Networked Parallel Computing. Cambridge: The MIT Press, England, 1994, 298p.
  13. Sun W.,Shu J.,Zheng W. Parallel Seismic Propagation Simulation in Anisotropic Media by Irregular Grids Finite Difference Method on PC Cluster. Computational Science and Its Applications – ICCSA, 2005: International Conference, Singapore, May 9-12, 2005. Proceedings. Part IV, pp. 762-771.