برای این کار ابتدا مدل عددی اولیه ای با استفاده از نرم افزار تفاضل محدود flac2d ساخته شد. سپس، با توجه به اینکه مقادیر دو پارامتر نسبت تنشهای اولیه برجا (k0) و مدول تغییرشکل پذیری توده سنگ (em) از عدم قطعیت بیشتری برخوردار بود، با بکار گیری داده های همگرایی سنجی و روش تحلیل برگشتی مقادیر این دو پارامتر برای مدل ساخته شده بهینه شد. پس از بهینه سازی مقادیر دو پارامتر مذکور، تحلیل وابسته به زمان پدیده مچاله شوندگی برای ایستگاه هایی که در آنها همگرایی زیاد دیواره مشاهده شده بود انجام شد. بدین منظور از مدل ویسکوپلاستیک برگر (چهار پارامتری) موجود در نرم افزار flac2d استفاده شد. برای مدلسازی وابسته به زمان می بایستی که مقادیر چهار پارامتر مدل ویسکوپلاستیک برگر بدست آیند که برای اینکار از روش تحلیل برگشتی داده های همگرایی سنجی ثبت شده استفاده شد. پس از بهینه سازی مقادیر پارامترهای مدل ویسکوپلاستیک برگر، مقطع بحرانی (با بیشترین روباره) تونل انتقال آب بابلک که هنوز حفاری نشده است برای دو حالت بدون پدیده مچاله شوندگی و با حضور پدیده مچاله شوندگی (تحلیل وابسته به زمان) مدلسازی شد. نتایج مدلسازی ها نرخ همگرایی زیاد دیواره (mm/day 3/0) را برای مقطع مذکور پیش بینی می کند که پیشنهاد می شود برای جلوگیری از همگرایی زیاد دیواره تونل، کف بند نیز در مقطع مذکور بهمراه سیستم نگهداری موقت نصب شود.