РОЗПОДІЛЕНЕ ДЖЕРЕЛО: РЕЗУЛЬТАТИ МОДЕЛЮВАННЯ ТА ПЕРСПЕКТИВИ ВИКОРИСТАННЯ ДЛЯ ЗАДАЧ СЕЙСМОЛОГІЇ

Д. Малицький; О. Муйла; О. Грицай; А. Павлова; О. Асташкіна; О. Обідіна; Е. Козловський

doi:10.17721/1728-2713.69.15.96-101

Authors

D. Malytskyy Carpathian Branch of Subbotin Institute of Geophysics NAS of Ukraine 3-b Naukova St., Lviv, Ukraine, 79060
O. Muyla Carpathian Branch of Subbotin Institute of Geophysics NAS of Ukraine 3-b Naukova St., Lviv, Ukraine, 79060
O. Hrytsaj Carpathian Branch of Subbotin Institute of Geophysics NAS of Ukraine 3-b Naukova St., Lviv, Ukraine, 79060
A. Pavlova Carpathian Branch of Subbotin Institute of Geophysics NAS of Ukraine 3-b Naukova St., Lviv, Ukraine, 79060
O. Astashkina Carpathian Branch of Subbotin Institute of Geophysics NAS of Ukraine 3-b Naukova St., Lviv, Ukraine, 79060
O. Obidina Carpathian Branch of Subbotin Institute of Geophysics NAS of Ukraine 3-b Naukova St., Lviv, Ukraine, 79060
E. Kozlovskyy Carpathian Branch of Subbotin Institute of Geophysics NAS of Ukraine 3-b Naukova St., Lviv, Ukraine, 79060

DOI:

https://doi.org/10.17721/1728-2713.69.15.96-101

Keywords:

extended source, seismic field, seismic moment tensor of the earthquake, isotropic medium

Abstract

The solution of the direct problem is presented for the displacement field on the free surface of layered isotropic medium using the matrix method. The results of the direct problem are used to determine the seismic moment tensor. An extended source is considered as a set of point sources, each one is presented by seismic moment tensor. An important aspect is that for the solution of inverse problem an analytical value of the direct problem is used, i.e. inversion for seismic tensor is realized by using solutions for displacement fields. The solution for extended sources is based in the fact that the wave field from such a source is the superposition of displacement fields from each point source. Thus, the statement of the direct problem is to determine the wave field on the free surface of layered half-space when the earthquake's focus is represented as an extended source in space and time. A method is described which determines the displacement field on the free surface in the spectral domain using the values of the shift for elementary sources as well as rise time and rupture time. Matrix method is used in case of seismic waves in horizontal layered half-space where heterogeneous medium is simulated by homogeneous isotropic layers with parallel boundaries. The earthquake's focus as an extended source is placed in a uniform layer. We have shown the transition from a redefined system of equations for determining a slip vector to the solution for the generalized inverse problem. The results of the inverse problem for determining the rupture plane were tested on the example of the events that took place near Malta (24.04.2011: 13h02m12s, 35.92N, 14.95E, Mw4.0)). For this event, the determination of the rise time and rupture time is shown. Correctness of the inverse problem is provided by determining of a functional in which the norm is minimized between the real data and parameters that are obtained using the proposed method. For a singular matrix it is suggested to use a singular decomposition.

References

Аки К., Ричардс П., (1983). Количественная сейсмология. Теория и методы. М.: Мир, 520 с. Aki K., Richards P., (1983). Quantitative seismology. Theory and methods. Moscow, Mir, 520 p. (In Russian).

Вербицкий Т.З. , Починайко Р.С. , Стародуб Ю.П. , Федоришин О.С., (1985). Математическое моделирование в сейсморазведке. Київ: Наук. думка, 275 с. Verbitsky T.Z., Pochinaiko R.S., Starodub Y.P., Fedorishin O.S., (1985). Mathematical modelling in seismic exploration. Kiev, Science Thought – Naukova Dumka, 275 p. (In Russian).

Малицький Д.В., (2010). Аналітично-числові підходи до обчислення часової залежності компонент тензора сейсмічного моменту. Геоінформатика, 1, 79–86. Malytskyy D.V., (2010). Analytical and numerical approaches to the calculation of time depending on the seismic moment tensor components. Geoinformatics, 1, 79-86. (In Ukrainian).

Малицький Д.В., (1998). Основні принципи розв'язання динамічної задачі сейсмології на основі рекурентного підходу. Геофіз. журн, 5, 96–98. Malytskyy D.V., (1998). Basic principles of solving of the seismology dynamic problems based on recurrent approach. Geophysics Journal, 5, 9698. (In Ukrainian).

Малицький Д.В., (2005). Про джерело сейсмічних хвиль. Геофіз. журнал, 27(2), 304–308. Malytskyy D.V., (2005). About the seismic waves source. Geophysics Journal, 27 (2), 304-308. (In Ukrainian).

Малицкий Д.В., (1994). Рекуррентный метод решения обратной динамической задачи сейсморазведки в вертикально-неоднородной среде. Геофиз. журнал, 16(3), 61–66. Malytskyy D.V., (1994). The recursive method for the inverse dynamic problem solving of seismic exploration in a vertically inhomogeneous medium. Geophysics Journal, 16 (3), 61-66. (In Russian).

Малицький Д.В., Муйла О.О., (2007). Про застосування матричного методу і його модифікацій для дослідження поширення сейсмічних хвиль у шаруватому середовищі. Теоретичні та прикладні аспекти геоінформатики, 124–136. Malytskyy D.V., Muyla O.O., (2007). About the application of the matrix method and its modifications for the seismic waves study in layered media. Theoretical and applied aspects of geoinformatics, 124-136. (In Ukrainian).

Chen Ji, (2005). Computer Simulation of Earth Movement that Spawned the Tsunami. California Institute of Technology.

Malytskyy D., D`Amico S., (2015). Moment tensor solutions through waveforms inversion. Messina, Mistral Servise sas, 25 p. 10. Müller G., (1985). The reflectivity method: a tutorial. Geophys. J., 58, 153–174.