視頻配套操作步驟地質統(tǒng)計學

1、SURPAC 軟件地質統(tǒng)計學（操作步驟）GEMCOM 國際軟件公司SURPAC 中國辦事處目錄第 1章地質域11.1 查看地質域1樣品組合32.1 在 SURPAC 軟件中提取地質域中的數據3基礎統(tǒng)計73.1 顯示柱狀圖7特異值104.1 特異值104.2 確定并處理特異值104.2.1 特異值下限的確定方法10第 2章第 3章第 4章4.2.2特異值的處理方法104.2.3用 SURPAC 軟件處理特異特異值114.2.4驗證處理結果11第 5章方差圖155.1 全方向向變異函數155.2 查看主方差圖的傾斜平面175.3 計算主方向方差圖195.3.1 查看已有的主方向方差圖結果195.3

2、.2計算主方向方差圖225.4 計算次主方向變差圖265.4.1 查看已有的次主方向方差圖結果265.4.2計算次主方向方差圖285.5 異向性橢球參數325.5.1 查看已有的橢球體325.5.2異向性橢球參數33第 6章交叉驗證376.16.1 交叉驗證37創(chuàng)建 3 維塊模型397.1 創(chuàng)建一個三維塊模型39通克立格法428.1 用異向性橢球在 3 維塊模型上進行普通克立格估值428.2 各向同性橢球在 3 維塊模型上進行普通克立格估值45第 7章第 8章I第9 章部分百分比計算489.1 計算塊模型中的塊體在 3DM 實體中的百分比48第10章 模型檢驗4910.1 勘探線數據與模型比較

3、4910.1.1 實體內塊模型百分比與礦體的實體剖面比較4910.1.2 塊體模型的克立格估值與組合樣品的對比5110.2 品位噸位曲線5410.2.1 創(chuàng)建品位噸位曲線5510.3 塊模型值的基礎統(tǒng)計5810.3.1 顯示塊模型值的基本統(tǒng)計5910.4 趨勢分析6010.4.1 導出塊質心6010.4.2 趨勢分析61II第 1 章地質域1.1 查看地質域1.打開 all_composites2.str。2.選擇顯示>全部隱藏。3.選擇顯示>點>標記。4. 按下面的窗體輸入信息，點擊執(zhí)行。5. 選擇顯示> 3D 網格.6. 按下面的窗體輸入信息，點擊執(zhí)行。7. 在圖形

4、工作區(qū)中進行點擊和拖動，將看到如下的旋轉的數據。1線串文件中的點所有鉆孔以 2 米為組合樣長的樣品品位。D1 描述字段包含金的組合品位，根據品位的高低為分不同的線串，見下表。8.再打開 ore_solid1.dtm 文件。實體模型表征了一個單獨的地質域，只有落在域內部的組合樣品才能用于估計位于地質域內的塊體模型。2第 2 章樣品組合2.1 在 SURPAC 軟件中提取地質域中的數據域ore_solid1.dtm表現(xiàn)的是QV1地質帶，下面將完成中僅在QV1域中提取組合樣品的操作：1.數據庫>打開/新建2.數據庫>組合>根據勘探工程點擊執(zhí)行定義地質帶：在下一個窗體，輸入“geol

5、ogy”表中的字段“rock”被定義為地質帶，“geology”表是區(qū)間型表，“rock”字段中的“QV1”標明了每個鉆孔穿過礦體“ore_solid1.dtm”的區(qū)間。33.文件>打開>線/DTM4. 顯示> 全部隱藏5. 顯示> 點> 標記46. 顯示> 點> 屬性7. 數據庫>顯示>鉆孔8. 查看> 觀察平面> 創(chuàng)建圖形剖面59.查看>觀察平面>后一個平面執(zhí)行 7 次上述可由運行宏文件_01_創(chuàng)建鉆孔組合.tcl來實現(xiàn)。6第 3 章基礎統(tǒng)計3.1顯示柱狀圖塊體模型>地質統(tǒng)計>基礎統(tǒng)計1.文件>

6、 從線文件導入數據2.點擊執(zhí)行。3.74. 統(tǒng)計> 報告5.點擊執(zhí)行。顯示輸出文件“raw_gold.not”內容。89第 4 章特異值4.1 特異值“特異值”是一個統(tǒng)計名詞，是指在數據母體中那些比全部數值的平均值或中位數高得多的數值。例如，下面的數據母體，236 將被認為是一個特異值。1 3 5 5 8 8 8 2 3 6特異值在變異函數計算時，能產生問題。另外，在估值中有特異值的存在，將產生不切合實際的結果。用于減少特異值的影響的辦法是對特異值用一個邊界品位替代，在上面的例子中，236 的數值將被“截取”， 用數值“9”替代。1 3 5 5 8 8 8 9另外一種辦法就是去掉特異值。

7、4.2 確定并處理特異值查看一下 gold_comp2.str 的柱狀圖和輸出的報告，將會發(fā)現(xiàn)，多數的品位數據都集中在 0 至 10g/t之間，同時也發(fā)現(xiàn)幾個比 10g/t 大的特異值。4.2.1特異值下限的確定方法1.傳統(tǒng)方法:用平均品位的 68 倍做為特異值的下限2.數理統(tǒng)計法:對組合的樣品進行基本統(tǒng)計分析，按 95%的置信區(qū)間確定特異值的下限.即:特異值下限=均值+1.96× 標準差3.本實例的特異值下限=3.828+1.96×6.831=17.2174.2.2特異值的處理方法1.傳統(tǒng)方法:用平均品位替代特異值102.截取法:將高于特異值的下限的部分截取掉,相當于用特

8、異值的下限替代特異值3.本實例的特異值下限為了避免特異值對估值的影響，采取了“截取”的辦法，即將超過 17g/t 的品位值替換為 17g/t。4.2.3用 SURPAC 軟件處理特異特異值1. 線文件工具>輸入下面信息：線串運算D1 描述字段將表達式的結果：iif(d1>17，17，d1)可表述為：如果D1的值大于17， 則D1的值等于17， 否則等于D1的原來值。2.點擊執(zhí)行。為檢驗線串運算的結果，你將在基礎統(tǒng)計窗口中對數據進行分析，4.2.4驗證處理結果1.塊體模型> 地質統(tǒng)計> 基礎統(tǒng)計112. 文件>從線文件導入數據3. 點擊執(zhí)行。12一個柱狀圖和一條描述

9、累計頻數曲線。注意這時數據中的最大值是 17。4. 統(tǒng)計>報告5.點擊執(zhí)行。顯示輸出文件“gold_cut17.not”的內容。13上述操作可由宏文件03_截取特異值.tcl.實現(xiàn)。14第 5 章方差圖5.1 全方向向變異函數1.塊體模型>地質統(tǒng)計> 方差模型2.變異函數>新建線文件方差圖3.點擊高級選項卡154.點擊執(zhí)行。5.變異函數>模型。把滯后距修改為7.6165.2 查看主方差圖的傾斜平面1.文件>打開>線/DTM2.顯示>為 DTM 著色173.文件>打開>線/DTM4. 顯示> 線> 帶顏色充填18上述操作可由

10、運行宏文件 05_data_geometry.tcl 來實現(xiàn)。5.3 計算主方向方差圖5.3.1查看已有的主方向方差圖結果1. 將display_major_axis1.str文件拖到圖形工作區(qū)。2. 將primary_variogram_map1.dtm文件拖到圖形工作區(qū)。193.點擊按鈕。4.顯示>為DTM著色。5. 將strike_dip_plane1.str文件拖到圖形工作區(qū)。6. 顯示> 線> 帶顏色充填207.顯示> 3D坐標網主方向方差圖顯示在圖形工作區(qū)中215.3.2計算主方向方差圖1.塊體模型> 地質統(tǒng)計>方差模型2.方差圖> 新建方

11、差圖3.點擊高級選項卡224.點擊執(zhí)行。5.移動滯后距滑塊，并點擊選項卡6. 變異函數> 模型237. 點擊并拖動基臺值和變程的標記點，進行擬合8. 方差圖> 選擇具有最大的連續(xù)性的方向249. 方差圖>保存為DTM10. 文件> 試驗方差圖11. 文件> 方差圖模型12. 差圖窗口也可以平鋪，點擊右鍵，彈出窗口，選TileWindows25上述操作可由運行宏文件_06_primary_variogram_map.tcl.tcl來實現(xiàn)5.4 計算次主方向變差圖5.4.1查看已有的次主方向方差圖結果1.選擇主程序窗口。2. 將display_semi_major_a

12、xis1.str文件拖到圖形工作區(qū)。3. 將secondary_variogram_map1.dtm文件拖到圖形工作區(qū)。264.顯示>為DTM著色。顯示> 3D 坐標網5.27主方向和次主方向方差圖顯示在圖形工作區(qū)中5.4.2計算次主方向方差圖1. 選擇方差圖窗口2. 方差圖> 次變差圖283.移動滯后距滑塊，并選第10個選項卡變異函數> 模型4.295.點擊并拖動基臺值和變程的標記點6.方差圖>保存為最有最大的連續(xù)性方向7.方差圖>保存為DTM308.文件> 保存> 試驗方差圖9.文件> 保存> 方差圖模型10. 方差圖窗口也可以平

13、鋪，點擊右鍵，彈出窗口，選 TileWindows。31上述操作可由運行宏文件07_secondary_variogram_map.tcl來實現(xiàn)。5.5 異向性橢球參數5.5.1查看已有的橢球體1.將 ellipsoid2.str 文件拖到圖形工作區(qū)中325.5.2異向性橢球參數1.方差圖>沿軸提取試驗方差圖2.文件>打開>方差模型33調節(jié)滯后距3.4.選擇 semi_major 選項卡，調整滯后距為 32，點擊并拖動基臺值和變程的標記點。345.選擇 minor 選項卡，調整滯后距為 24，點擊并拖動基臺值和變程的標記點。上述操作可由運行宏文件08_anisotropy_e

14、llipsoid.tcl.6.方差圖> 創(chuàng)建異向性橢球報告。357.方差圖>橢球體觀察儀。8.選擇文件> 關閉36第 6 章交叉驗證6.1 交叉驗證1.重置圖形工作區(qū)。2.塊體模型> 地質統(tǒng)計方差驗證3.輸入下面信息，點擊執(zhí)行。3738第 7 章創(chuàng)建 3 維塊模型7.1 創(chuàng)建一個三維塊模型1. 刪除 geostats_example.mdl 和 QV1.CON 兩個文件。2. 塊體模型>新建/打開3. 塊體模型> 顯示394. 塊體模型>約束>新建圖形約束405. 塊體模型> 屬性> 新建上述操作可由宏文件_09_create_blo

15、ck_m.tcl來實現(xiàn)41第 8 章通克立格法8.1 用異向性橢球在 3 維塊模型上進行普通克立格估值1 塊體模型>估值>普通克立格法4243塊體模型顯示如下：448.2 各向同性橢球在 3 維塊模型上進行普通克立格估值1. 塊體模型>估值>普通克立格法4546上述操作可用宏文件_10_ordinary_kriging.tcl來實現(xiàn)47第 9 章部分百分比計算9.1 計算塊模型中的塊體在 3DM 實體中的百分比重置圖形工作區(qū)。1.2.打開 geostats_example.mdl，并進行顯示.將約束文件 qv1.con 拖到圖形工作區(qū)中估值> 部分百分比3.輸入下

16、面信息，點擊執(zhí)行。4.塊體模型> 保存5.6.塊體模型>關閉上述操作可由運行宏文件_11_partial_percentage.tcl 來實現(xiàn)。48第 10 章 模型檢驗10.1勘探線數據與模型比較10.1.1 實體內塊模型百分比與礦體的實體剖面比較1.將 ore_solid.dtm 文件拖入到圖形工作區(qū)2.將 geostats_example_final.mdl 文件拖入到圖形工作區(qū)3.將 qv1.con 文件拖入到圖形工作區(qū)塊體模型塊體模型顯示4.塊體模型顯示根據屬性為模型著色5.496.查看觀察平面創(chuàng)建圖形剖面507.顯示點屬性10.1.2 塊體模型的克立格估值與組合樣品的對

17、比1.重置圖形工作區(qū)，關閉塊模型2. 將 ore_solid1.dtm 拖入圖形工作區(qū)3. 將 surpac.ddb 拖入圖形工作區(qū)4.數據庫顯示515.將 section_example.mdl 文件拖入圖形工作區(qū)中6.將 文件拖入圖形工作區(qū)中塊體模型塊體模型顯示7.塊體模型顯示根據屬性為模型著色8.529.查看觀察平面創(chuàng)建圖形剖面53顯示點屬性10.10.2品位噸位曲線另一個方法是報告品位、噸位，并構造一條品位噸位曲線。5410.2.1 創(chuàng)建品位噸位曲線1.將 geostats_example_final.mdl 文件拖入圖形工作區(qū)中2.將 文件拖入圖形工作區(qū)中塊體模型塊體模型顯示3.塊體

18、模型顯示根據屬性為模型著色4.555.塊體模型塊體模型報告565710.3塊模型值的基礎統(tǒng)計塊模型值的基礎統(tǒng)計也是檢驗塊模型的法。58Anisotr opic Kriged Gold體積噸位Anisotropic Kriged GoldAnisotropic Kriged Gold0.0 -> 1.000001.0 -> 2.02371744743481.82861450.92.0 -> 3.0129266425853282.466366763.693.0 -> 4.059767811953553.454119692.534.0 -> 5.02523075046134.432233498.265.0 -> 6.0900041800085.4971777.366.0 -> 7.026998539966.44347821.57.0 -> 8.0.44196105.38.0 -&