استخراج امواج ریلی زلزله های ایران با استفاده از سه مولفه شتاب نگاشت در حوزه زمان-فرکانس

امواج ریلی نسبت به سایر امواج، اثرات مخرب تری بر روی سازه های ساخت بشر دارند. شناسایی و استخراج امواج ریلی از شتاب نگاشت ها در حوزه زمان-فرکانس از دقت بالایی برخوردار است. در تحقیق حاضر، از تبدیل موجک پیوسته و تبدیل استوکول برای انتقال سه مولفه شتاب نگاشت هر ایستگاه زلزله به حوزه زمان-فرکانس و استخراج امواج ریلی استفاده شده است. با استفاده از مشخصات حرکت بیضوی، امواج ریلی از شتاب نگاشت ها در حوزه زمان-فرکانس قابل استخراج می باشند. همچنین با انتقال فاز حرکت قائم و با استفاده از ضرب داخلی نرمالایز شده، می توان امواج ریلی پس رونده، پیش رونده و نیز زاویه انتشار این امواج را در حوزه زمان-فرکانس به‌طور مجزا استخراج نمود و درنهایت پاسخ را در فضای زمان ارائه داد. در مطالعه حاضر، الگوریتم های مذکور در زبان برنامه‌نویسی متلب کدنویسی شده است. صحت تحلیل های صورت گرفته، با استفاده از امواج ریلی استخراج شده از سیگنال های مصنوعی و نیز داده های زلزله چی چی، ارزیابی شده است. در نهایت زلزله های ایران شامل طبس، بم و منجیل مورد بررسی قرار گرفته و امواج ریلی آنها  استخراج و پارامتر های لرزه ای آنها محاسبه شده است. نتایج حاصله نشان می دهد که در شتاب نگاشت ها سهم امواج ریلی پس رونده بیشتر از سهم امواج ریلی پیش رونده و سایر امواج است.

Extracting Rayleigh Waves of Iran Earthquakes from Three-Component Signals in Time-Frequency Domain

Rayleigh waves can be more destructive in comparison with other types of waves during earthquakes. Identification and extraction of Rayleigh waves from seismic records can be more exactly performed in timefrequency domain. In the present study, Continuous Wavelet Transform and Stockwell Transform have been used to transfer threecomponents of seismic signals to timefrequency domain. Rayleigh waves can be extracted based on elliptical characteristics of these waves. On the other hand, Retrograde and Prograde Rayleigh waves and their propagation azimuth can be separately extracted with phase transferring of vertical component as well as normalized inner product technique. In the present study, the mentioned two algorithms have been programmed. Accuracy of the methods has been verified with using synthetic and ChiChi earthquake data. Finally, Rayleigh waves and their seismic parameters have been extracted from Tabas, Bam and Manjil earthquake records. Rayleigh waves of the earthquakes have been analyzed and compared with using different seismic parameters. According to the results, Retrograde Rayleigh waves are more predominant in comparison with Prograde Rayleigh waves in earthquake signals.IntroductionExtraction of Rayleigh waves is one of the important tasks in seismology and applied geophysics due to their destructive effects. A lot of study has been conducted in order to develop an effective method for the extraction of Rayleigh waves. Most of these methods have been presented based on the ellipticity characteristics of Rayleigh waves. Each of them has their advantages and disadvantages. In this study, two of them have been investigated, programmed and compared with each other.Methodology and ApproachesRayleigh waves have been extracted from threecomponent signals in timefrequency domain in both of the investigated methods. In algorithm 1, continuous wavelet transform has been used to transfer three components of signals from time to timefrequency domain. However, Stockwell transform has been used to perform this approach in algorithm 2. Stockwell transform retain the absolute phase of each localized frequency component, the characteristic of which is required in algorithm 2 for timefrequency transform. On the other hand, Rayleigh waves have been extracted based on the instantaneous reciprocal ellipticity and phase difference in algorithm 1 and 2 respectively. The most important advantages of algorithm 2 is that retrograde and prograde Rayleigh waves and their propagation azimuth can be separately extracted while it is not possible in algorithm 2.In this study, both of these algorithms have been clearly illustrated by flowcharts and programmed in Matlab. Synthetic signals and ChiChi earthquake signals have been used for verification. After verification, results of both algorithms have been compared, and finally they have been employed in extracting Rayleigh waves of Tabas, Bam and Manjil earthquakes. Finally, maximum displacements of Rayleigh waves and total earthquake signals have been investigated and compared in these three earthquakes. In addition, energy discharge of extracted Rayleigh waves during the earthquakes have been investigated. Results and ConclusionsResults of the study showed that algorithm 1 can be effectively used if there is not required to extract retrograde and prograde Rayleigh waves and their propagation azimuth separately because the algorithm is simpler and faster in comparison with algorithm 2. However, algorithm 2 should be employed on the condition that the extraction of retrograde and prograde Rayleigh waves and their propagation azimuth have been required. In addition, it has been concluded that Rayleigh waves have more energy in earthquake signals compared to other types of waves. Consequently, elliptical polarization is more than linear polarization in earthquakes. In comparison of retrograde and prograde Rayleigh waves, it has been concluded that retrograde Rayleigh waves have more energy and therefore they can be more destructive. Finally, it has been resulted that ratio of Rayleigh waves are different in earthquakes. Induced displacement by Rayleigh waves in Manjil earthquake was more than Bam and Tabas earthquakes.