Разработка и реализация математических алгоритмов построения миграционных изображений упругих трещиноватых геологических сред (1187420), страница 5
Текст из файла (страница 5)
70, no. 1. — S30–S37.5530.Yan J., Sava P. Isotropic angle-domain elastic reverse-time migration // Geophysics. — 2008. — Vol. 73, no. 6. — S229–S239.31.Tsvankin I., Lynn H. B. Special section on azimuthal dependence of P-waveseismic signatures—Introduction // Geophysics. — 1999. — Vol. 64, no.4. — Pp. 1139–1142.32.Nihei K.
T. Natural fracture characterization using passive seismic illumination // Lawrence Berkeley National Laboratory. — 2003.33.Zhang M., Gao K., Huang L. Anisotropic revese-time migration for imagingfracture zones at Eleven-Mile canyon // Workshop on Geothermal ReservoirEngineering. — Stanford. 2017.34.Numerical computation of wave propagation in fractured media by applyingthe grid-characteristic method on hexahedral meshes / V. I. Golubev [et al.] //Computational Mathematics and Mathematical Physics. — 2015.
— Vol. 55,no. 3. — Pp. 509–518.35.Zhdanov M. S. Geophysical Inverse Theory and Regularization Problems.Vol. 36. — First. — Elsevier Science, 2002. — (Methods in Geochemistryand Geophysics).36.Claerbout J. F. Earth Soundings Analisys: Processing versus Inversion. —2004. — URL: sep.stanford.edu/sep/prof/pvi.pdf ; AccessedAugust 2015.37.Tromp J., Tape C., Liu Q. Seismic tomography, adjoint methods, time reversal and banana-doughnut kernels // Geophysical Journal International. —2005.
— Vol. 160, no. 1. — Pp. 195–216.38.Irons T. Marmousi2 model. — 2014. — URL: www.reproducibility.org/RSF/book/data/marmousi2/paper_html/paper.html ;Accessed June 2017.39.URL: www.reproducibility.org.5640.Beylkin G. On applications of unequally spaced fast Fourier transform //Mathematical Geophysics Summer School, Stanford. — 1998. — Vol. 127.41.Beylkin G., Kurcz C., Monzón L. Grids and transforms for band-limited functions in a disk // Inverse Problems. — 2007.
— Vol. 23, no. 5. — Pp. 2059–2088.57ПриложенияA. Программная реализация миграции БорнаПри наивной реализации описанный в разд. 1.1 алгоритм миграции обладает вычислительной сложностью O (T R S N ), где T, R, S, N – числа записанных отсчетов, приемников, источников и узлов сетки. Если же входныеданные являются сейсморазрезами постоянного удаления, приемники располагаются на плоскости постоянного значения одной из координат, а периодих расположения по двум другим кратен шагам сетки по соответствующейкоординате (или шаг кратен периоду), то выражение для условия визуализации приобретает вид свертки, и его можно эффективно вычислять с помощьюбыстрого преобразования Фурье или неравномерного быстрого преобразования Фурье (см.
напр. [40; 41]). Например, когда приемники расположенына горизонтальной плоскости, вычислительная сложность алгоритма можетбыть уменьшена до O (T Nz Nx Ny log Nx log Ny ), где Ni – число узлов сеткипо соответствующей координате.Также, в рамках поставленной задачи вычисление миграционных параметров среды на разных глубинах (или при разных значениях координаты, при постоянном значении которой расставлены приемники) можно осуществлять независимо, что позволяет параллелизовать алгоритм с практически идеальным ускорением (табл. 1 на с. 59).58Число ядер Время расчета, мин Эффективность12911.21451.3970.9974730.9995590.9896490.9957420.9978370.9849320.99410300.97311270.98512250.973Таблица 1: Характеристики расчета миграционного изображения длятестовой модели при разном числе ядер.59.