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

نوع مقاله : Articles

نویسندگان

1 پژوهشگاه علوم و فنون هسته‌ای ایران، تهران

2 پژوهشگاه بین‌المللی زلزله‌شناسی و مهندسی زلزله، تهران، ایران

چکیده

در این پژوهش سعی شده است که پارامترهای لرزه­ای برای منطقه قم محاسبه گردد. برای نیل به این هدف از روش حل معکوس تعمیم یافته استفاده شده است که در ابتدا توسط اندروز [1] ارائه شد. برای استفاده از این روش ابتدا بانک داده­های شتاب­نگاری ایجاد و داده­های شتاب­نگاری همگی از شبکه شتاب­نگاری سازمان تحقیقات راه، مسکن و شهرسازی استخراج شد. برای تهیه بانک داده­ها، تمامی رکوردهای شتاب مربوط به منطقه مورد مطالعه از بقیه رکوردها جدا شد و سپس مؤلفه‌های قائم و افقی هر رکورد استخراج گردید. پس از تهیه بانک داده­ها، تمامی رکوردها تصحیح شد و موج برشی مستقیم از کل رکورد جداسازی گردید. طیف فوریه برای رکوردهای به‌دست‌آمده محاسبه شد و در ۲۰ فرکانس بین ۴/۰ تا ۱۵ هرتز با فواصل لگاریتمی مساوی از طیف فوریه نمونه­برداری به عمل آمد تا داده‏ها جهت استفاده در الگوریتم حل معکوس تعمیم یافته آماده شوند. با استفاده از روش حل معکوس تعمیم یافته و با استفاده از روش اندرو برای مقید کردن درجه آزادی معادله، ماتریس معلومات و مجهولات مسئله تهیه شد. معادله ماتریسی با کمک روش جداسازی مقادیر تکین حل شد. در محاسبه فاکتور کیفیت امواج برشی مدل گسترش هندسی از مدل معتضدیان [2] استفاده شده است.اثرات چشمه، مسیر حرکت امواج برشی و تأثیرات ساختگاهی که ایستگاه ثبت کننده در آن واقع شده است هرکدام به­طور مجزا محاسبه شده­اند. در این منطقه پارامتر کیفیت امواج برشی برابر   محاسبه شده است. مدل استفاده شده برای تحلیل اثر چشمه مدل  می­باشد و برای محاسبه افت تنش، فرکانس گوشه، بزرگای گشتاوری و غیره از روش جستجوی نقطه‌به‌نقطه استفاده شده است و این پارامترها برای تمامی زلزله­های بانک داده­ها محاسبه گردیده­اند. تشدید ساختگاهی نیز برای تمامی ایستگاه­های ثبت کننده محاسبه شده است.

کلیدواژه‌ها


عنوان مقاله [English]

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

نویسندگان [English]

  • Mohammad Reza Aram 1
  • Babak Khazaie Moghaddam 2
1 Nuclear Science and Technology Research Institute, Tehran, Iran
2 International Institute of Earthquake Engineering and Seismology, Tehran, Iran
چکیده [English]

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; Alborz-Azerbaijan (North and Northwest of Iran), Kopet Dagh (Northeast of Iran), Zagros (Southwest of Iran), Makran (Southeast of Iran) and Central-East of Iran. In this study, Qom region located in Central-East 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. Ground-motion 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 ground-motion 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 Hough-Anderson 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 time-domain data into frequency-domain 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 trade-off 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 Boore-Joyner generic rock site model. Database is constructed from 39 number of strong-motion 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 de-noising method. Shear wave extracted from all records using Husid method for beginning of S-wave and Kinoshita method for end of S-wave and Fourier spectrum is sampled in frequencies between 0.4 and 15 Hz. A system of over-determined equations is resulted.
Stress drop, corner frequency, moment magnitude and seismic radius is calculated for each earthquake using Grid-search method. The relation between S-wave 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.4-15 Hz. Results of this study can further be used for a better evaluation of seismic hazard in a region of interest.

کلیدواژه‌ها [English]

  • Generalized Inversion Method
  • Strong-Motion Parameters
  • Quality factor
  • Stress Parameter
  • Site Amplification
  • Qom Province
  1. Andrews, D.J. (1986) Objective determination of source parameters and similarity of earthquakes of different size. Earthquake Source Mechanics. American Geophysical Union, 37, 259-267.
  2. Motazedian, D. (2006) Region-specific key seismic parameters for earthquakes in Northern Iran. Bull. Seismol. Soc. Am., 96(4),1383-1395.
  3. Hartzell, S. (1992) Site response estimation from earthquake data. Bull. Seism. Soc. Am., 82(6), 2308-2327.
  4. Berberian, M. (2006) Contribution to the seismotectonics of Iran (part 2). Geological Survey of Iran, 39, 429-517.
  5. Nowroozi, A. (1976) Seismotectonic provinces of Iran. Bull. Seism. Soc. Am., 66(4), 1249-1276.
  6. Khademi, M.R., Nayeri, A. (1976) Seismotectonic provinces of Iran. Iranian Committee on Large Dams; Bulletin 6.
  7. Mirzaei, N., Mengtan, G., Yuntai, C. (1998) Seismic source regionalization for seismic zoning of Iran: Major seismo tectonic provinces. J. Earthquake Pred. Res., 7(4), 465-95.
  8. Babaahmadia, A., Safaei, H., Yassaghi, A., Vafa, H., Naeimi, A., Madanipour, S., Ahmadi, M. (2010) study of Quaternary structures in the Qom region, West Central Iran. Journal of Geodynamics, 50, 355-367.
  9. Ansari, A., Noorzad, A., Zare, M., (2007) Application of wavelet multi-resolution analysis for correction of seismic acceleration records, J. Geophys. Eng., 4(4), 1-16.
  10. Husid, P. (1967) Gravity Effects on the Earthquake Response of Yielding Structures. Report of Earthquake Engineering Research Laboratory, California Institute of Technology, Pasadena, California.
  11. Kinoshita, SH. (1994) Frequency-Dependent Attenuation of Shear Waves in the Crust of the Southern Kanto Area, Japan. Bull. Seism. Soc. Am., 84(5), 1387-1396.
  12. Boore, D.M. (1983) Stochastic simulation of high-frequency ground motions based on seismological models of the radiated spectra. Bull. Seism. Soc. Am., 73(6), 1865-1894.
  13. Anderson, J.G. and Hough, S.E. (1984) A Model for the Shape of the Fourier Amplitude Spectrum of Acceleration at High Frequencies. Bull. Seismol. Soc. Am., 74(5), 1969-1993.
  14. Aki, K. (1967) Scaling law of seismic spectrum. J. Geophys. Res., 72(4), 1217-1231.
  15. Yoshimoto, K., Sato, H., and Ohtake, M. (1993) Frequency-dependent attenuation of P and S waves in the Kanto area, Japan, based on the coda normalization method. Geophys. J. Int., 114, 165-174.
  16. Zafarani, H., Hassani, B. and Ansari, A. (2012) Estimation of earthquake parameters in the Alborz seismic zone, Iran using generalized inversion method. Soil Dynam. Earthq. Eng., 42, 197-218.
  17. Ghasemi, H., Kamalian, N., Hamzehloo, H., Beitollahi, A. (2006) Determination of quality factor of direct shear waves, Qs, in Alborz region using near-field strong motion records of Kojor earthquake in frequency range of 1~32 Hz. Physics of the Earth and Space, 31,103-112.
  18. Mousavi, M., Zafarani, H., Noorzad, A., Ansari, A., Bargi, KH. (2007) Analysis of Iranian strong motion data using the specific barrier model. J. Geophys. Eng., 4(4), 1-14.
  19. Raoof, M., Herrmann, R.B., and Malagnini, L. (1999) Attenuation and excitation of three-component ground motion in southern California. Bull. Seismol. Soc. Am., 89(4), 888-902.
  20. Brune, J.N. (1970) Tectonic stress and the spectra of seismic shear waves from earthquakes. J. Geophys. Res., 75(26), 4997-5009.
  21. Boore, D. (2003) Simulation of ground motion using the stochastic method. Pure Appl. Geophys., 160, 635-76.
  22. Engdahl, E.R., Jackson, J.A., Myers, S.C., Bergman, E.A., Priestley, K. (2006) Relocation and assessment of seismicity in the Iran region. Geophys. J. Int., 167(2), 761-778.