視頻配套操作步驟地質(zhì)統(tǒng)計(jì)學(xué)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

18、模型顯示根據(jù)屬性為模型著色4.555.塊體模型塊體模型報(bào)告565710.3塊模型值的基礎(chǔ)統(tǒng)計(jì)塊模型值的基礎(chǔ)統(tǒng)計(jì)也是檢驗(yàn)塊模型的法。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 -&