مدلسازی آماری ناپیوستگی های توده سنگ در سه بعد، بر اساس پارامترهای هندسی درزه‎ها اعم از موقعیت، اندازه و چگالی آنها در واحد حجم توده سنگ صورت می‎گیرد. تعیین اندازه و چگالی ناپیوستگی ها در واحد حجم توده سنگ به این دلیل که صفحات درزه ها، درون توده سنگ بوده و نمی‎توان آن‎ها‎ را مستقیماً اندازه‎گیری کرد بسیار دشوار است. در این تحقیق ناپیوستگی های طبیعی اطراف تونل انحراف آب سد قردانلو در سازند تیرگان، شبیه‎سازی آماری گردیدند. در این مدلسازی، ناپیوستگی‎ها بصورت دیسک فرض شدند. توزیع‎های تئوری قطر درزه‎ها و همچنین چگالی آنها در واحد حجم توده سنگ با توجه به توزیع‎های تجربی بدست آمده از برداشت ناپیوستگی‎های ساختگاه، تخمین زده شدند. پس از مدلسازی ناپیوستگی ها، ارتباط آنها در دو بعد و تشخیص چند ضلعی های ناشی از تقاطع طول های اثر، با استفاده از یک برنامه کامپیوتری توسعه داده شده به نام rockblock که در محیط matlab اجرا می شود مورد بررسی قرار گرفت. حجم و طول سطح لغزش بلوک‎های اطراف تونل محاسبه و نتایج بدست آمده با نتایج حاصل از نرم‎افزار udec مقایسه گردید. هر دو روش مدلسازی عددی و آماری ناپیوستگی‎ها، نواحی ناپایدار تونل را یکسان تشخیص دادند. در نهایت سیستم نگهداری مناسب، برای تونل طراحی شد.