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

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