تخمین پارامترهای جنبش نیرومند زمین با استفاده از روش حل معکوس تعمیم یافته در منطقه قم

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

Evaluation of Strong-Motion Parameters in Qom Province in Iran Using Generalized Inversion Method

Being located on the seismic belt, Iran has experienced frequent and sometimes intense seismic activities. Many researchers discussed about seismotectonic characteristics of Iran, and based on similarities in seismotectonic parameters, they divided Iran into different tectonic provinces. Most acceptable seismotectonic divisions for Iran is provided by Mirzaie which introduces five tectonic provinces; AlborzAzerbaijan (North and Northwest of Iran), Kopet Dagh (Northeast of Iran), Zagros (Southwest of Iran), Makran (Southeast of Iran) and CentralEast of Iran. In this study, Qom region located in CentralEast province is discussed, which is a seismic active region and faced large earthquakes such as Avaj (22/6/2002) with magnitude 6.1 epicentered in 35.71˚N and 49.02˚E, Silakhor(30/3/2006) with magnitude 6.2 epicentered in 33.69˚N and 48.95˚E and Kahak (18/6/2007) with magnitude 5.7 epicentered in 34.54˚N and 50.94˚E, hence it is very important to have estimates of earthquake parameters such as anelastic and geometrical attenuation, average stress parameter, and site amplification for this region, which has high seismicity rates and especially has populated areas with high seismic risk. In this study, generalized inversion method is used which is common approach for deriving the modeling parameters from the empirically observed data. Generalized inversion method introduced by Andrews. Andrews recast the spectral ratios approach into an inverse problem for earthquakes and stations with multiple records. This method has been later improved by other researchers. Groundmotion amplitudes are affected by source, path and site effects. Source term includes the effect of magnitude, stress parameter and the type of faulting on the observed groundmotion amplitude. Path term includes both anelastic and geometrical attenuation. Site term includes the effect of the near surface features like geology and topography on the observe data that consists of site amplification and high frequency attenuation. High frequency attenuation is calculated by estimating kappa using HoughAnderson method. Convolution of these effects is what stations are recording as earthquake records. Generalized inversion method introduces a way to divide recorded spectrum into mentioned effects. At first, by transforming timedomain data into frequencydomain and then with the aid of natural logarithm operator, a linear equation is achieved. Afterward, we can derive modeling parameters for each of the effective terms. In this study, for source term, the Brune source model (omega square) is used. To remove the tradeoff between anelastic attenuation and geometrical spreading in path effect, geometrical spreading model from previous studies is used. In this formulation, all the site spectra may be multiplied by an arbitrary function of frequency and all the source spectra divided by the same function without changing the fit to the data that adds a degree of freedom to our system, which in this study is resolved by using predefined site term for reference site (Arak station) as it is located on hard rock (Vs30 >700 m/s) and can be evaluated using BooreJoyner generic rock site model. Database is constructed from 39 number of strongmotion records, recorded in 14 stations from 9 earthquakes with magnitude greater than 4 within target region. At first, horizontal and vertical components are extracted (horizontal component is considered as geometric mean of two horizontal components) and corrected using wavelet denoising method. Shear wave extracted from all records using Husid method for beginning of Swave and Kinoshita method for end of Swave and Fourier spectrum is sampled in frequencies between 0.4 and 15 Hz. A system of overdetermined equations is resulted. Stress drop, corner frequency, moment magnitude and seismic radius is calculated for each earthquake using Gridsearch method. The relation between Swave energy and moment magnitude is expressed by , regional stress drop is 81 bar and frequency dependent quality factor is  that aligns well with previous studies in this region. Calculated quality factor indicates that the region of study has high seismicity rates. Also site amplification is calculated for all the stations in the selected frequency range of 0.415 Hz. Results of this study can further be used for a better evaluation of seismic hazard in a region of interest.