مدل سازی سه بعدی و ارزیابی ذخیره کانسار مس دره زرشک یزد با استفاده از روش های تخمین کلاسیک و زمین آماری

در پژوهش حاضر به منظور شناخت کامل کانسار مس دره زرشک یزد، مدل سه بعدی چاه نگار لیتولوژی گمانه ها، نیمرخ های عرضی لیتولوژی و عیارسنجی منطبق بر پروفیل های پوشش دهنده کل محدوده و مدل سه بعدی لیتولوژی کانسار به کمک نرم افزار rockworks ترسیم شد. بر این اساس لیتولوژی غالب منطقه را سنگ آهک، سنگ های اسکارنی، ماسه سنگ، دیوریت، گرانودیوریت، کوارتزیت و مونزونیت تشکیل می دهند و سطح محدوده نیز بیشتر توسط روباره، شیل و آهک پوشیده شده است. بسته به شرایط کانسار برای تخمین ذخیره از روش کلاسیک عکس مجذور فاصله و روش زمین آماری کریجینگ با کمک نرم افزار surpac استفاده شد. بر اساس واریوگرافی کانسار با استفاده از نرم افزار wingslib کانسار دارای ناهمسانگردی هندسی است و مقدار شعاع تجسس بیشینه در راستای با آزیموت صفر درجه برابر با 130 متر، شعاع تجسس متوسط در راستای با آزیموت 90 درجه برابر با 100 متر و شعاع تجسس کمینه در راستای قائم با آزیموت صفر درجه برابر با 40 متر بدست آمد. بعد از مدلسازی هندسی و ساخت مدل بلوکی کانسار، تخمین ذخیره کانسار به ازای عیارحدهای 0.1، 0.2، 0.3، 0.4 و 0.5 درصد مس انجام شد و نتایج آنها با یکدیگر مقایسه شد. نتایج محاسبه ذخیره نشان می دهد که اختلاف میزان ذخیره محاسبه شده با روش های مذکور نسبتاً زیاد و مثلاً به ازای عیار حد 0.1% در حدود 1547600 تن یعنی 5.7% می باشد. اگرچه روش عکس مجذور فاصله از دقیق ترین روش های کلاسیک تخمین ذخیره است، اما روش زمین آماری کریجینگ با در نظر گرفتن ساختار فضایی منطقه، به دلیل نااریب بودن و داشتن کمترین واریانس تخمین، دقیق تر بوده و از درجه اعتبار بیشتری برخوردار است.

3-d modeling and reserve evaluation of yazd, dareh-zereshk copper deposit using geometric and geostatistical estimation methods

in the present research, to fully recognize the yazd, darreh-zereshk copper deposit, 3-d multiple-log of the boreholes lithology, cross-sections of lithology and assay data corresponding to profiles covering the whole area and 3-d lithology model of the deposit were plotted using rockworks software. based on this, limestone, skarns, sandstone, diorite, granodiorite, quartzite and monzonite perform the dominant lithology of the region and the surface of the area is mostly covered by overburden, shale and lime. depending on the deposit conditions, to estimate the deposit reserve, inverse distance squared method as a geometric method and geostatistical kriging method were used through surpac software. according to the deposit variography using wingslib software, the deposit has geometric anisotropy and maximum range of 130 m with azimuth zero, mean range of 100 m with azimuth 90 and minimum range of 40 m on the vertical direction with azimuth zero. after geometric modeling and producing block model of the deposit, the reserve was estimated for cut-off grades of 0.1, 0.2, 0.3, 0.4 and 0.5% copper as well as their results were compared. the results of the reserve calculation show that the difference of the calculated reserves with applied methods is relatively high, i.e. for cut-off grade of 0.1% is about 1547600 tones or 5.7%. although, inverse distance squared method is the most accurate method among the geometric methods of reserve estimation, geostatistical kriging method considering spatial structure of the region, being unbiased and having minimum estimation variance is more accurate with higher degree of validity.