ناپیوستگی ها (به طور کلی درزه ها) عامل اصلی ناهمگونی و ناپایداری در مناطق توده سنگی مورد مطالعه مهندسین عمران، زمین شناسی، معدن و نفت می باشند. بنابراین نیاز به یک نمایه واقعی از هندسه توده سنگ احساس می گردد که انواع روش های مدل سازی می توانند این نیاز را برای مهندسان برطرف سازند. 311 درزه به روش خط اسکن از کوه رحمت برداشت شد و درزه های برداشت شده با نرم افزار stereonet دسته بندی شده و 4 دسته درزه اصلی به همراه درزه های تصادفی بدست آمد. با استفاده از مشخصات درزه های برداشت شده (شیب و امتداد)، دقت توزیع های آماری (نرمال، نمایی، گاما و بتا) بررسی گردید، همچنین برازش درزه های برداشت شده برای مدلهای فوق تعیین گردید. مبنای مقایسه توزیع های آماری بررسی گردیده، ضریب برازش (r2)، مربعات خطا (sse) و ریشه میانگین مربعات خطا (rmse) بوده است. برای 311 درزه برداشت شده، توزیع بتا نمایش دقیق تری نسبت به توزیع های دیگر برای شیب درزه ها نشان داد. در صورتیکه توزیع امتداد درزه ها از مدلهای فوق پیروی نمی کرد. بنابراین توسط روش درون یابی اسپلاین مکعبی قطعه ای یک بعدی بهترین تابع چگالی احتمال برای آن بدست آمد، سپس توسط توابع چگالی احتمال و توابع توزیع تجمعی بکار گرفته شده بالا توزیع درزه ها در توده سنگ آهک به صورت یک بعدی شبیه سازی گردید. نتایج ارائه شده ارتباط نزدیکی بین مدل ارائه شده و درزه های واقعی نشان می دهد. برای شبیه سازی دو بعدی درزه های برداشت شده، با نرم افزار surfer المان بندی درزه ها انجام گرفت. برای بدست آورد تابع چگالی توأم امتداد و شیب درزه ها کدنویسی در matlab انجام گرفت که به ازای هر المان براساس درون یابی اسپلاین مکعبی دو بعدی، یک رویه محاسبه گردید. برای تولید توأم امتداد و شیب درزه ها کدنویسی در matlab انجام گرفت. سپس درزه های تولید شده با نرم افزار stereonet مدل سازی گردید. نتایج ارائه شده ارتباط نزدیکی بین مدل ارائه شده و درزه های واقعی نشان می دهد