چکیده:
زمینلرزهها حاصل گسلش و پویایی زمینساختی هر ناحیه هستند و نسبت به بزرگا، فضا و زمان توزیع فرکتالی دارند. در این تحقیق پارامترهای لرزهخیزی b-value و FD به عنوان رویکرد کمی فرکتالی در تحلیل لرزهخیزی زاگرس شمال غرب بکار رفته است. پارامتر b با توزیع فراوانی-بزرگای گوتنبرگ-ریشتر و پارامتر FD به روش مربعشمار محاسبه و همچنین با استفاده از 30 خوشه تمرکز زمینلرزه و عوامل مهم، پهنههای لرزهخیز با روش آنتروپی مشخص شدند. کاهش b-value رخداد زمینلرزههای با فراوانی کم و بزرگای بالا و افزایش FD عدم خوشه بندی و توزیع یکنواخت زمینلرزهها را نشان می دهد. همبستگی منفی این دو مؤید زمینساخت فعال است. نتایج نشان میدهند که افزایش FD با کاهش b-value (همبستگی منفی 60 درصد) همراه است. منطقه A (غرب کرمانشاه) ، بالاترین میزان (FD (1/02 و پایینترین میزان b (0/78) را داراست که به معنی توزیع بالای تنش در منطقه (و شاهد آن زمینلرزه اخیر کرمانشاه با بزرگای 7/3) است. منطقه E (محل اتصال گسلهای جبهه کوهستان و بالارود)، مشابه A است اما زمینلرزههای آن بزرگای کمتری دارند (5 و 6 ریشتر). منطقه F (محل پیوند گسلهای اصلی و جبهه کوهستان زاگرس)، با وجود تراکم بالای گسل و زمینلرزه، نقطه مقابل A است که بیانگر رها شدن تنش در قالب زمینلرزههایی کوچک است. منطقه C (حاشیه غربی با کمترین تمرکز گسل و فراوانی زمینلرزه)، مشابه F است. مناطق D و B (مجاور گسلهای معکوس و رورانده زاگرس، شرایطی زمینساختی متوسطی دارند. نتایج مدل آنتروپی نتایج FD را تأیید میکند و رابطه معکوسی با مقادیر b دارد.
Introduction The purpose of present study is seismicity analysis of Lorestan folded arc and its adjacent thrust belt using quantitative methods. To reach this aim we performed analysis of seismicity using quantitative methods to find possible vertical and horizontal changes in seismic activity across the main Zagros faults of the northwestern part of Zagros. Firstly, we used fractal geometry and frequency-magnitude distribution of earthquakes by using FD and b-value parameters, respectively. Here b-value is the main factor in Gutenberg-Richter empirical relation which indicates the exponential distribution of earthquake magnitudes (Godano et al, 2014; 1765). This parameter also is known as fractal dimension (Mirabedini & aghatabay,2015: 60). FD is fractal dimension of earthquake epicenters distribution which has been calculated by box-counting method (Turcotte 1997). On the other hand Entropy model has been applied to specify potential of seismicity by using effective factors and 30 points of earthquake concentration. The study area in northwestern part of Zagros was divided to the simply folded arc of Lorestan and faulted-folded belt of high Zagros. Several main faults pass through the area from NW to SE and divide its main morphotectonic units as High, folded and foredeep parts of Zagros (Berberian, 1995: 193). Material and methods Data in this research can be divided to two part: parameters of earthquakes (magnitude, depth, location of epicenter) and linear data of faults and anticline/syncline axes. These data have been changed into new layers by GIS software extensions (density of epicenter and depth of earthquakes, density of faults and anticline/syncline axes, distance of fault and epicenter of earthquakes, interpolation of epicenter of earthquakes) to be applied in Entropy model, in other hand frequency of magnitude clusters and surface distribution of earthquakes are main data in Gutenberg–Richter relation and Fractal methods respectively. Numerical results of mentioned methods have been calculated and drawn in excel software. Gutenberg–Richter relation (Gutenberg & Richter 1944) is defined as Log N(m)= a-bm, where N is the cumulative number of earthquakes with magnitude larger or equal to m, a is a constant (seismicity level) and b is the slope of frequency-magnitude (size distribution) (Godano, 2014). To calculate fractal dimension of distribution of earthquake epicenters, box counting method suggested by Turcotte (1997) were applied by using Hausdorff dimension, which in two quantity of size (side length of grids) and number (number of grid boxes containing earthquake) are used to calculate FD value (Schuller et al, 2001: 3). In the other section, earthquake epicenters are divided to several clusters with different magnitude, then kernel density of each cluster was applied and subsequently, the maximum concentration of each magnitude cluster was determined as a point layer. Followingly, by overlaying these point layer with effective layers in seismicity analysis, their characteristics was extracted. Finally, an Entropy matrix was calculated and using experts rating and computing the layer’s weight, seismic zones were identified (Zonggi, et al, 2010). Result and discussion Estimated b-value indicates approximately reciprocal values compared with FD values. Decrease in b-value reveals that stress level and probability of large magnitude earthquakes occurrence is quite high and increase in FD shows that earthquakes are not clustered and are distributed homogeneously along a line in understudy area. Calculated number-size values for earthquakes represent both partial and popular FD changes. Based on partial FD, three populations can be classified: (a) Background with FD larger than popular FD; (b) Threshold with FD lower than 0.7: and (c) Anomaly with FD more than two. Based on popular FD, distribution of earthquakes is linear and transition to chaos phase is not predicted. Comparison between maximum values of Entropy zoning and FD values for each box indicates that these two values show 93% correlation (regardless of the C box values due to incompatibility with value of other boxes). Conclusion Areas with high FD value and low b-value are more tectonically active. The box labeled A which represent western parts of Kermanshah in folded Zagros, has the highest FD value (1.02) and lowest b-value (0.78). The box labeled F in southern east part is in contrast with it (highest b value:1.02 and one of the lowest FD value: 0.89) in understudy areas. E (Balarud fault) and D (High and folded Zagros) parts have almost the same FD and b values. FD and b values in B (high Zagros) are equal and less than the aforementioned areas. C (that contains a part of mountain front fault) has the lowest value of FD and same b-value as B and the changes of Entropy max values are same as FD values.
خلاصه ماشینی:
رابط ه فراواني-بزرگاي ارائه شده توسط گوتنبرگ -ريشتر بيانگر ــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــ 1 Bachmanov 2 Hirata 3 Gutenberg-Richter 4 Aki & Richards 5 Sornette et al 6 Frequency-Magnitude 7 Hirata 8 Telesca et al 9 Omori خودهمساني ١ زمين لرزه ها است )گدانو و همکاران ٢، ٢٠١٤: ١٧٦٥(.
نامپالي و همکاران ١٨ )٢٠١٨( ناهمگني فضايي توالي زمين لرزه اي ٢٠١٥ نپال را با استفاده از مقادير فرکتالي و b-value ارزيابي ــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــ 1 Self-similarity 2 Godano et al 3 Kaga & Knopoff 4 Aki 5 Mandelbrot 6 King 7 Turcotte 8 Smalley et al 9 Huang & Turcotte 10 Wang 11 Öncel et al 12 Henderson 13 Legrand 14 Wyss et al 15 Mandal & Rastogi 16 Bayrak & Bayrak 17 Wu et al 18 Nampally et al نمودند که نتايج بيانگر ناهمگني بالا و تنش پايين و توزيع غير اتفاقي اما با خودسازماني بحراني پس لرزه ها است .
در واقع در اين تحقيق وضعيت لرزه خيزي زاگرس شمال غربي را با استفاده از روش هاي فرکتالي و مدل آنتروپي مورد ارزيابي قرار گرفته است .
ــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــ 1 Vernant 2 Blanc et al 3 Berberian & king 4 Hessami 5 Jackson & Fitch 6 Maggi et al (به تصویر صفحه مراجعه شود) شکل : ٢الف : نمودار گلسرخي راستاي بيشينه تنش در زاگرس غربي (فروهيد، ١٣٨٥: ٦٠)، ب :عمق کانوني زمين لرزه هاي زاگرس (مگي و همکاران ، ٢٠٠٠: ٦٢٩)، نقشه گسل هاي اصلي فعال و پراکنش رومرکز زمين لرزه هاي بزرگ تاريخي (پيش از ١٩٠٠) و دستگاهي (١٣٩٦-١٢٧٨)زاگرس شمال غرب (مؤسسه ژئوفيزيک دانشگاه تهران ، ٢٠٠٨: irsc.