子域估值法在甲瑪銅多金屬礦資源量估算中的應用

在傳統品位插值過程中，一個礦化域一般只使用單個搜索橢球體。然而對于某些礦體形態復雜或者局部礦化變化強烈的礦化域時，單一的搜索橢球體并不能滿足對該礦化域進行合理估值，從而影響了資源儲量估算可信度和資源分類，以及后續的采礦設計和排產等技術工作。以甲瑪銅多金屬礦為研究對象，利用Surpac三維軟件子域估值法對甲瑪銅多金屬礦進行估值和評價，尋求解決這種局部形態或礦化變化強烈的礦體資源儲量估算方法，子域估值法是一種很好的方法，更能確保資源儲量可信度。

1 地質數據庫

地質數據庫是數據庫技術在礦山地質工作中的實際應用，是礦山礦產資源儲量估算的基礎，是采礦設計的依據

，是礦山生產管理的有效手段。地質數據的完整性和可靠性直接影響礦山生產經營和決策

。因此，建立地質數據庫是實現地質資源信息數字化管理的最基礎工作

。

通過收集甲瑪多金屬礦歷年以來地質鉆孔、探槽、坑內鉆、巷道刻線等數據，獲得此次地質資源數據庫地質屬性信息，并建立相應的數據表，以存儲相應的地質資源信息

。

2 礦化域三維模型

2.1 礦化域范圍圈定

在充分研究甲瑪銅多金屬礦床地質特征、成礦理論和控礦因素的基礎上，根據礦床地質特征圈定礦化域。在同一礦化域內，以Cu≥0.3%、Mo≥0.03%、Pb≥0.3%、Zn≥0.5%作為邊界品位進行混圈。若一樣品中任意一種金屬元素品位大于或等于上述邊界品位時，該樣品則圈入礦化域內。

在圖書館創客空間運動如火如荼開展的熱潮中，學者們對創客空間建設的必要性進行了探討，主要分為兩種觀點，一種觀點認為建設創客空間非常必要，但也有學者從不同角度提出冷思考。

2.2 建立礦化域三維實體模型

以勘探線剖面進行礦化域地質解譯，按地質規范進行外推。完成整個礦化域解譯后，再以三角網化工具手動連結三角網，最終形成驗證通過可用于塊體模型估算的礦化域實體模型。

3 數據處理

3.1 樣長統計分析

提取礦內樣品經統計分析得到樣本數107777個，樣長均值1.67m、中值為2.00m，結合甲瑪礦山地質勘探及生產取樣實際，得到工程組合樣長為2m。

式中：H為測壓管水頭，m；c為水錘波速度，m/s；x，t分別為水錘波沿直線傳播的距離和時間；g為當地重力加速度，m/s2；f(v)是速度的非線性函數。

3.2 特異值處理

根據對勘探工程組合數據分析得到Cu品位最小值0.0001%、最大值32.694%、Cu均值0.80%，標準差0.99，Sichel-t 0.94，Cu元素的品位累積頻率曲線總體上呈一條直線,說明Cu元素分布基本上符合正態分布,同時，Cu元素品位累積頻率曲線(圖1)上部又都存在少量離散點，且偏離直線方向,這說明樣品中存在特高品位

。通過數據轉換得到Cu(log10)累計頻率直方圖（圖2），偏度-0.75，峰度4.40，基本符合正態分布條件，使用數理統計的方法進行特高品位下限值處理方法，95%的置信區間

(CI)的上限值=Cu均值+1.96標準差=0.8+1.96×0.99，得到Cu特高品位下限值2.74%。

由計算結果可知，滲透系數的變化與排水量呈正比關系，其中有限元計算結果的比例系數為3.66，方法一計算結果的比例系數為3.13。而滲透系數的變化對剩余水頭高度則沒有影響。

通過線串運算的方法截取Cu品位的特異值，統計分析Cu品位最小值0.0001%、最大值2.74%、均值0.74%、Sichel-t值為0.89%，特高品位處理合理，變化系數95%，Cu元素在礦化域分布均勻。

4 子域劃分

（1）陡Cu1子域變異函數及橢球體

子域變異函數研究：

應用實體內外區分點功能把組合樣線文件按礦化域產狀分割成三個獨立的礦化域，以每個子礦化域為獨立礦化域進行變異函數分析研究，通過變異函數球狀模型擬合，得到各子域變異函數的塊金值、基臺值、變程及橢球體參數。

2017年6月21日的國務院會議中，李克強總理提出“分享經濟”對于新常態下的中國經濟發展具有重要意義。在當今的中國“分享經濟”充分利用互聯網技術，不斷的創造出新的產品和服務模式，為供給側改革、化解產能、提供工作崗位做出了重要貢獻。與“分享經濟”相匹配，人才市場也出現了許多“分享經濟”的特征，出現了全新的雇傭模式。李克強總理提到，他在南方省市考察時發現當地出現了“零工市場”，工人們沒有固定的雇主，工作按小時收費，一天可以服務多個雇主。

在以往甲瑪銅多金屬礦資源儲量估算時，一個礦化域使用一個橢球體參數，而礦體具有明顯的陡→緩→陡的變化特征(圖3)。因此，使用單一橢球體(圖4)不能代表整個礦化域產狀，不能合理對整個礦化域插值評價，從而影響資源儲量估算精度。為保證資源儲量估算可信度，在研究以往地質資料的基礎上

，結合礦化域產狀將其分為三個子域。以Cu元素為例，分為陡Cu1、緩Cu、陡Cu2三個子域（圖5）。

當前，解決人民日益增長的美好生活需求和不平衡不充分的發展之間的矛盾，改善企業社區居民生活條件已是勢在必行。改造企業社區，政府發揮主導作用的同時，積極鼓勵多元力量參與企業創新型社區的轉型和建設。利用、改造、發展、引導，多管齊下，實現共享共治，打造全新的泰煤家園！

（2）緩Cu子域變異函數及橢球體

分別對三個獨立子礦化域插值,組合樣與估值結果驗證，得到整體插值Cu(綠色線)及各子域插值Cu(黃色線)Swath曲線(圖9)。

1.1 臨床資料 選取2016年1月至2018年1月在我院呼吸內科門診就診及住院的長期吸煙哮喘患者。入選條件：符合中華醫學會制定的支氣管哮喘防治指南哮喘診斷標準；年齡18-60歲之間；吸煙史≥10包年。兩組患者治療前一般資料比較差異無統計學意義（P＞0.05）。

（3）陡Cu2子域變異函數及橢球體

（4）模型驗證

兩組干預前腕背伸、肘屈曲、肩外展、膝屈伸活動度組間無顯著差異(P>0.05);干預后觀察組各項指標均明顯優于對照組(P<0.05)。見表1。

以上三個子域變異函數通過交叉驗證得到滿足理想結果的三個條件：一是，真值與估計值誤差的均值趨近于“0”；二是，誤差方差與克里格估計方差比值趨近于“1”

；三是，兩個標準差之間的置信區間大于等于95%。

陡Cu2子域樣本數2448個，Cu品位最小值0.0001%、最大值2.74%、均值0.46%、Sichel-t0.65%，變化系數124%，偏度2.21，峰度8.30，通過log10對數轉換后統計分析偏度0.70、峰度3.04,基本符合正太分布。通過對緩Cu2變異函數球狀模型擬合（圖8）和交叉驗證(表3)，得到陡Cu2橢球體參數，方位75°，傾伏角0°，傾角45°，主/次軸比1.11，主/最小軸比5.77。

緩Cu子域樣本數56972個，Cu品位最小值0.0001%、最大值2.74%、均值0.78%、Sichel-t0.95%，變化系數88.8%，偏度1.22，峰度3.85，通過log10對數轉換后統計分析偏度1、峰度5，基本符合正太分布。通過對緩Cu變異函數球狀模型擬合（圖7）和交叉驗證(表2)，得到緩Cu子域橢球體參數，方位120°，傾伏角0°，傾角45°，主/次軸比1.58，主/最小軸比5.07。

陡Cu1子域樣本數33782個，Cu品位最小值0.0001%、最大值2.74%、均值0.70%、Sichel-t0.73%，變化系數104%，偏度1.50，峰度4.4，通過log10數據轉換后統計分析偏度0.54、峰度3.79,符合正太分布。通過對陡Cu1變異函數球狀模型擬合（圖6）和交叉驗證(表1)，得到了Cu橢球體參數，方位105°，傾伏角0°，傾角45°，主/次軸比1.28，主/最小軸比1.74。

從Swath曲線得出，使用子域估值法估值，估值結果更接近原始組合樣，有利于資源儲量估算可信度和資源分類，以及后續的采礦設計和排產等技術工作。

（5）資源儲量對比。

由表5得出，用子域估值法估值，以Cu 0.3%為邊界品位，從礦石體積、礦石量、金屬量都相應減少，Cu資源儲量體積減少0.023%，噸位減少0.024%，金屬量減少0.005%。

5 結論

（1）在甲瑪銅多金屬礦根據礦體產狀分割成若干個子域，分別進行子域變異函數研究及插值得到較好的應用。

（2）對于解決這種局部形態和礦化變化劇烈的礦體，子域估值法是一種很好的方法，更能確保資源儲量可信度。

（3）子域估值法有效剔除了廢石量，提高了資源儲量估算的精度，能夠更好的指導后續的采礦生產活動，增加經濟效益。

（4）子域界限的劃分難度和工作量比較大，對于特別復雜礦體的估算耗時且容易出錯不適用。

[1]李雪峰,宋志偉,呂海棟,等.一種礦山地質資源儲量管理方法:,CN106529755A[P].2017.

[2]張國亮,張海峰,周勇.基于Surpac的某珍珠巖礦三維建模與露天境界圈定[J].采礦技術,2012,012(001):11-13.

[3]車海龍.Micromine軟件在某金礦床資源儲量估算中的應用[J].吉林地質,2020,39(3):7.

[4]鄭文寶,黎楓佶,唐菊興,等.基于Micromine軟件下地質統計學在甲瑪矽卡巖型銅多金屬礦儲量計算中的應用[J].地質與勘探,2011,47(4):726-736.

[5]尚瓊.安徽省沙溪銅金礦床鳳臺山礦段三維建模及成礦元素空間分析[D].合肥工業大學,2014.

[6]周娣.滑繩溜車保護在多繩摩擦提機中的應用[J].中國金屬通報,2020(4):2.

[7]唐攀,唐菊興,唐曉倩,等.傳統方法和地質統計學在礦產資源/儲量分類中的對比分析[J].金屬礦山,2013.