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

نویسندگان

1 استادیار دانشکده مهندسی علوم‌زمین، دانشگاه صنعتی اراک، اراک، ایران

2 دانش‌آموخته کارشناسی‌ارشد، دانشکده مهندسی علوم‌زمین، دانشگاه صنعتی اراک، اراک، ایران

چکیده

در پژوهش حاضر به­منظور شناخت کامل کانسار مس دره­زرشک یزد، مدل سه­بعدی چاه­نگار لیتولوژی گمانه­ها، نیمرخ­های عرضی لیتولوژی و عیارسنجی منطبق بر پروفیل­های پوشش­دهنده کل محدوده و مدل سه­بعدی لیتولوژی کانسار به­کمک نرم­افزار Rockworks ترسیم شد. بر این اساس لیتولوژی غالب منطقه را سنگ­آهک، سنگ­های اسکارنی، ماسه­سنگ، دیوریت، گرانودیوریت، کوارتزیت و مونزونیت تشکیل می­دهند و سطح محدوده نیز بیشتر توسط روباره، شیل و آهک پوشیده شده است. بسته به شرایط کانسار برای تخمین ذخیره از روش کلاسیک عکس مجذور فاصله و روش زمین­آماری کریجینگ با کمک نرم­افزار Surpac استفاده شد. بر اساس واریوگرافی کانسار با استفاده از نرم­افزار WinGslib کانسار دارای ناهمسانگردی هندسی است و مقدار شعاع تجسس بیشینه در راستای با آزیموت صفر درجه برابر با 130 متر، شعاع تجسس متوسط در راستای با آزیموت 90 درجه برابر با 100 متر و شعاع تجسس کمینه در راستای قائم با آزیموت صفر درجه برابر با 40 متر بدست آمد. بعد از مدلسازی هندسی و ساخت مدل بلوکی کانسار، تخمین ذخیره کانسار به­ازای عیارحدهای 1/0، 2/0، 3/0، 4/0 و 5/0 درصد مس انجام شد و نتایج آنها با یکدیگر مقایسه شد. نتایج محاسبه ذخیره نشان می­دهد که اختلاف میزان ذخیره محاسبه­شده با روش­های مذکور نسبتاً زیاد و مثلاً به­ازای عیار حد 1/0% در حدود 1547600 تن یعنی 7/5% می­باشد. اگرچه روش عکس مجذور فاصله از دقیق­ترین روش­های کلاسیک تخمین ذخیره است، اما روش زمین­آماری کریجینگ با در نظر گرفتن ساختار فضایی منطقه، به­دلیل نااریب بودن و داشتن کمترین واریانس تخمین، دقیق­تر بوده و از درجه اعتبار بیشتری برخوردار است. 

کلیدواژه‌ها


عنوان مقاله [English]

3-D modeling and reserve evaluation of Yazd, Dareh-Zereshk copper deposit using geometric and geostatistical estimation methods

نویسندگان [English]

  • R. Ahmadi 1
  • S. Fazeli Nasab 2
1 Assist. Prof., Dept., of Mining Engineering, Arak University of Technology, Arak, Iran
2 M. Sc., (graduated), Dept., of Mining Engineering, Arak University of Technology, Arak,
چکیده [English]

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. 

کلیدواژه‌ها [English]

  • Yazd
  • Darreh-Zereshk copper deposit
  • 3-D modeling
  • WinGslib software
  • Inverse distance squared
  • Kriging
  • Surpac software
احمدی، ر (1399) مقایسه نتایج روش­های زمین­آماری خطی و غیرخطی در مدلسازی و ارزیابی ذخیره کانسار مس نارباغی شمالی ساوه، فصلنامه زمین­شناسی ایران، سال 14، شماره 56، ص 43-59.
احمدی، ر (1400) مدلسازی زمین­شناسی- اکتشافی کانسار مس نارباغی شمالی ساوه و تخمین ذخیره کانسار با استفاده از رویکردهای بلوک­بندی، مدل شبکه دوبعدی و انباشتگی دوبعدی، نشریه زمین شناسی اقتصادی، جلد 13، شماره 2، ص 435-462.
احمدی، ر.، بهارلوئی، ز (1400) ارزیابی میزان تطابق کیفی و کمی داده­های ژئوفیزیکی با عـیارسـنجی کانـسار مـس دره­زرشک یزد به­منظور پیشنهاد موقعیت بهینه گمانه­های اکتشافی جدید، نشریه زمین­شناسی مهندسی دانشگاه خوارزمی، جلد 15، شماره 1، ص 34-15.
احمدی، ر.، سادات کودهی، س. م (1397) تخمین و رده­بندی ذخیره کانسار سرب و روی اره­گیجه رباط خمین استان مرکزی با استفاده از روش­های زمین­آماری، یافته­های نوین زمین­شناسی کاربردی، دوره 12، شماره 24، ص 39-53.
بهارلویی، ز (1398) پیشنهاد موقعیت حفر گمانه های اکتشافی جدید براساس میزان تطابق و تلفیق اکتشافات ژئوفیزیکی با داده­های عیارسنجی کانسار مس دره زرشک یزد، پایان نامه کارشناسی­ارشد، دانشگاه صنعتی اراک، 150ص.
حاج ملاعلی، ع.، علوی­نائینی، م (1372) نقشه زمین­شناسی خضرآباد با مقیاس 1:100000، سازمان زمین­شناسی کشور.
حسنی­پاک، ع. ا (1379) مدلسازی کانسارهای فلزی و غیرفلزی و کاربرد اکتشافی آن­ها، انتشارات دانشگاه تهران، 512 ص.
حسنی­پاک، ع. ا.، شرف­الدین، م (1380) تحلیل داده­های اکتشافی، انتشارات دانشگاه تهران، 987ص.
حکمت­نژاد، ا.، حسنی­پاک، ع. ا (1390) تخمین توزیع عیار مس با استفاده از تخمینگر غیرخطی کریجینگ انفصالی، نشریه علمی-پژوهشی مهندسی معدن، دوره ششم، ص 11-20.
خسروتهرانی، خ (1372) چینه­شناسی کرتاسه زیرین در نواحی غرب و جنوب غربی یزد.
خوئی، ن.، قربانی، م.، تاج­بخش، پ (1378) کانسارهای مس در ایران، انتشارات سازمان زمین­شناسی و اکتشافات معدنی کشور، 421 ص.
سادات­کودهی، س. م (1396) ارزیابی ذخیره کانسار سرب و روی رباط خمین با اسـتفاده از روش­های کلاسـیک و زمین­آماری، پایان­نامه کارشناسی­ارشد، دانشگاه صنعتی اراک، 128ص.
شرکت کنـگاوران (1384) گـزارش زمیـن­شناسـی مـعـدن دره­زرشک.
شرکت ملی صنایع مس ایران (1381) گزارشات تهیه شده در مقاطع مختلف.
صمان کاو، مهندسین مشاور (1385) گزارش مطالعات تکمیلی ژئوفیزیک به روش پلاریزاسیون القایی (IP/RS) در محدوده معدن مس دره­زرشک استان یزد.
فاضلی­نسب، ص (1399) مدلسازی و ارزیابی ذخیره کانسار مس دره­زرشک یزد با استفاده از روش­هـای زمـین­آمـاری، پایان­نامه کارشناسی­ارشد، دانشگاه صنعتی اراک، 115ص.
مدنی، ح (1371) مبانی زمین­آمار، انتشارات دانشگاه صنعتی امیرکبیر، واحد تفرش، 659ص.
مدنی، ح (1376) اصول پی­جویی، اکتشاف و ارزیابی ذخایر معدنی، انتشارات خانه فرهنگ، 816ص.
نبوی، م. ح (1355) دیباچه ای بر زمین­شناسی ایران، سازمان زمین­شناسی و اکتشافات معدنی کشور.
یزدانی، م.، علی­نیا، ف (1397) بکارگیری مطالعات زمین­آماری برای شناسایی آنومالی­های Au-Cu-Fe-W در دره­ سه هزار تنکابن، علوم زمین، شماره 109، ص 83-96.
Annels, A. E (2012) Mineral deposit evaluation: A practical approach, Springer Science & Business Media, 436 p.
Bohling, G (2007) SGeMS Tutorial Notes, 26 p.
Choudhury, S (2015) Comparative study on linear and non-linear geostatistical estimation methods, A Case Study on Iron Deposit, Procedia Earth and Planetary Science, 11: 131–139.
Daya, A (2012) Reserve estimation of central part of Choghart north anomaly iron ore deposit through ordinary kriging method. International Journal of Mining Science and Technology, 22 (4): 573-577.
Erickson, Jr. A. J (1992) Geological interpretation, modeling and representation, SME Mining Engineering Handbook (2nd edition), Hartman, H. (ed.), Chapter 5. SME-AIME, New York, 333-343.
France Coffimine (1989) Co-ali abad dareh zereshk report.
Remy, N., Boucher, A., Wu, J (2006) SGeMS User’s Guide, 129 p.
Rendu, J. M (1981) An introduction to geostatistical methods of mineral evaluation, South frican institute of mining and metallurgy, 84 p.  
Taghavi, A., Maanijou., M., Lentz., D., Sepahi, A. A (2019) Partial sub-pixel and pixel-based alteration mapping of porphyry system using ASTER data: regional case study in western Yazd, Iran. International Journal of Image and Data Fusion, 10 (4): 300-326.