بهينه‌سازي تعيين محل زمين‌لرزه‌‌هاي رخ داده بين سال‌هاي 2006 تا 2013 در شمال غرب ايران با استفاده از پارامتر‌هاي سَمت (آزيموت) و پرتو

Optimization of earthquake locations in northwest of Iran during the period 2006-2013 using azimuth and ray parameters

يكي از كاربردهاي مهم زلزله‌شناسي تعيين كانون زمين‌لرزه‌‌‌ها است. دقت بيشتر در اين فرايند موجب كسب نتايج بهتري در ساير بررسي‌ها مانند بررسي ساختار پوسته، لرزه‌خيزي منطقه و تعيين صفحه گسل مسبب زلزله‌‌‌ها مي‌‌‌شود. روش‌هاي كلاسيك تعيين محل بر پايه معادلات گايگر (1910) بنا شده است. در اين روش با استفاده از جمله اول بسط تيلور به‌كمك يك مدل سرعتي معلوم براي منطقه، رابطه بين زمان رسيد و مختصات كانوني زلزله خطي‌سازي مي‌‌‌شود. پاسخ بهينه در روش‌هاي خطي شده، از كمينه كردن اختلاف داده‌‌هاي محاسبه شده و مشاهده شده به‌ روش كمترين مربعات گام‌به‌گام صورت مي‌‌‌پذيرد. اغلب به علت نبود پوشش سَمتي (آزيموتي) مناسب و فقط به واسطه استفاده از داده زمان سير، خطاي قابل‌توجهي در تعيين محل كانون زلزله‌‌‌ها ايجاد مي‌‌‌شود. در راستاي كاهش خطاي ناشي از تعيين محل، پژوهشگران از سَمت و پارامتر پرتو براي بهينه‌سازي نتايج استفاده مي‌‌‌كنند. در اين تحقيق سعي شده است بهينه‌سازي موقعيت رويدادهاي لرزه‌‌اي در محدوده °44 تا °50 طول شرقي و°36 تا °40 عرض شمالي واقع در شمال غرب ايران صورت گيرد. 11200 زمين‌لرزه با بزرگي برابر يا بيش از 4/1 در مقياس ناتلي، رخداده بين سال‌هاي 2006 تا 2013 ، ابتدا فقط با استفاده از زمان سير موج P و سپس با استفاده از زمان سير، پارامتر پرتو و سَمت موج P تعيين محل شد. نتايج حاصل بيانگر تاثير بسيار مهم پارامتر پرتو و سَمت به همراه زمان سير، در كاهش خطاي افقي تعيين محل زمين‌لرزه‌‌‌ها به ميزان قابل‌توجهي است.

An important aim of seismology is to determine focuses of earthquakes. The more precisely the process is carried out, the better results will be achieved in various studies, for example, in vestigation of the crustal structure and seismogenic zone or the determination of the plane of faults causing earthquakes. The classic methods of locating are based on Geigerʹs equations (1910). The relation between the arrival time and coordination of the earthquake is linearized in this method using the first term of the Taylorʹs expansion and a velocity model available for the study area. The optimum responses in the linearized methods are obtained by minimizing the differences between the observed and calculated data using iterative Least Squares (L.S.) method. A considerable error is usually observed in determination of earthquake focus because of the absence of a proper azimuthal bearing, as well as the mere use of the travel time data. With the purpose of reduction of the errors, researchers have used azimuth and ray path parameters for optimization of their performance. The study area is located within 44-50o E and northern latitude 36-40oN, which is located in northwestern Iran, forming part of the central Alborz-Azerbaijan tectonic region. In this research the earthquakes of a magnitude equal or greater than 1.4 in the body wave scale occurred in the area within the period from 2006 to 2013, included 11,200 events, first were located using only P wave traveltimes and then process was repeated by adding ray path and P wave azimuthal parameters data. To this effect, the aforesaid earthquakes were classified into two general classes, consisting of earthquakes of magnitudes 1.4-3.4 and earthquakes of magnitudes exceeding 3.4. Concerning the first class of earthquakes, processing was made only for the data obtained from the stations in the local seismic network of the Institute of Geophysics, University of Tehran. The velocity model used in this part of the work was local Velocity Model (Bayramnejad, 2008) that is specific to the study area. Because of a large number of earthquakes, each of the aforesaid classes of earthquakes were subdivided into two subclasses consisting of the earthquakes with azimathal gaps less than 180 degrees and those with azimuthal gaps greater than 180 degrees. The results obtained by the earthquakes locating are represented by histograms which indicate that utilization of ray parameter and azimuthal parameter considerably reduce the horizontal error in relocating of earthquakes, specially for the earthquakes with the azimuthal gaps exceeding 180 degrees compared to utilization of mere travel time data. For the second class of earthquakes, we have used all the data from the stations located within the whole seismic range of Iran. It is evident that locating of earthquakes is carried out with a greater precision where more data is available. The velocity model used in this part was similar to the one used for the whole. Optimization of relocating of earthquakesare indicated by histograms for the magnitude of 3.4-4.6earthquakes.Results indicated reduction in the horizontal error when the azimuthal and ray path parameters have been used. For the magnitude 4.6-6-2 earthquakes, it is proved, by drawing a certainty ellipse, that the use of azimuthal and ray parameters may optimize locating of earthquakes,and lower the dimensions of the certainty ellipse.