مدل‌سازي عددي گسل جبهه كوهستان (MFF) در رخداد زلزله‌ 21 آبان 1396 سرپل‌ذهاب به‌روش عنصر مرزي

در اثر زلزله 21 آبان 1396 سرپل ذهاب با بزرگاي 7.3=Mw تغييرمكان هاي سطحي به‌ويژه زمين لغزش بسيار بزرگ مله كبود و پس لرزه هاي زيادي به‌صورت خوشه هاي پراكنده در منطقه رخ داد. وقوع پس لرزه ها به‌صورت خوشه اي نشاني از فعاليت گسل هاي كوچك از قبل موجود مي باشد. در مطالعه حاضر گسل مسبب زلزله اصلي (گسل جبهه كوهستان (MFF)) با استفاده از كد عنصر مرزي سه‌‌بعدي 3d-def مدل‌سازي شد و ميزان تغيير مكان هاي سطحي به‌ويژه جابه‌جايي حاصل از زمين‌لغزش مله كبود و همچنين محل گسل هاي ثانويه با استفاده از ميانگين تنشهاي اصلي بيشينه و كمينه به‌عنوان معياري براي تعيين محل گسل هاي ثانويه استخراج گرديد. به‌منظور صحت‌‌سنجي نتايج مدل‌سازي براي جابه‌جايي ها، از مقايسه نتايج با داده هاي به‌دست‌آمده از روش تداخل‌‌سنجي راداري (InSAR) و روش تركيبي رديابي افست (combined offset tracking) مشاهده شد كه نتايج دو روش تقريباً به همديگر نزديك مي باشند و اختلافاتي نيز در موقعيت جابه‌جايي ها و دامنه تغييرات آنها در دو روش ديده مي شود به گونه اي كه دامنه تغييرات جابه‌جايي حاصل از داده هاي ماهواره اي در منطقه زلزله، 35- تا 80 سانتي‌متر است در حالي‌كه اين محدوده براي نتايج حاصل از مدل‌سازي انجام شده، 15- تا 100 سانتي‌متر مي باشد. همچنين براي محل گسل هاي ثانويه، ميانگين تنشهاي اصلي بيشينه و كمينه حاصل از مدل‌سازي تقريباً در محل هاي داراي تراكم پس لرزه به‌ويژه در قسمت جنوبي گسل بيشينه مي باشد و عمق هاي با بيشترين ميانگين تنشهاي اصلي بيشينه و كمينه  با عمق‌هايي كه تراكم پس لرزه زياد است، به‌جز در عمق 13كيلومتري به پايين تقريباً يكي مي باشند. علاوه‌بر اين براي تعيين ميزان حساسيت مدل نسبت به پارامترهاي ورودي تحليل حساسيت انجام شد كه در نتيجه آن تغييرات پارامترهاي شيب و عرض گسل، بيشترين اثر را بر خروجي جابه‌جايي و همچنين تغييرات پارامترهاي مدول يانگ و طول گسل بيشترين اثر را بر خروجي ميانگين تنشهاي اصلي بيشينه و كمينه داشتند.