數(shù)據(jù)庫實戰(zhàn)：建庫、導(dǎo)入與避坑指南)
簡介這份PDF文檔是Surpac Vision地質(zhì)數(shù)據(jù)庫的官方培訓(xùn)教程面向礦山地質(zhì)、測量及采礦工程領(lǐng)域的技術(shù)人員與初學(xué)者幫助其系統(tǒng)掌握地質(zhì)數(shù)據(jù)的建庫、管理與三維可視化方法。資源包內(nèi)僅含1個PDF文件大小約2.32MB內(nèi)容按章節(jié)組織涵蓋數(shù)據(jù)庫介紹、結(jié)構(gòu)解析、創(chuàng)建流程、編輯瀏覽管理、三維空間顯示、地質(zhì)解譯與剖面品位計算、組合樣品、地質(zhì)統(tǒng)計與分析、線文件剖面輸出及信息提取報告等模塊并附有3D轉(zhuǎn)換實例與多邊形資源評估兩個附件。教程從表與字段分析、空庫創(chuàng)建、數(shù)據(jù)導(dǎo)入驗證到映射數(shù)據(jù)庫逐步展開同時講解鉆孔顯示風(fēng)格、隱藏刷新、勘探線定義與礦體圈定等實操要點便于讀者對照軟件界面完成練習(xí)。目前已有221人學(xué)習(xí)適合作為Surpac入門與日常功能查閱的參考手冊。1. 從一份 PDF 說起Surpac 地質(zhì)數(shù)據(jù)庫到底在管什么如果你手里正拿著一份名為“2-surpac地質(zhì)數(shù)據(jù)庫教程.pdf”的資料大概率你已經(jīng)被 Surpac 里那套“數(shù)據(jù)庫”概念繞得有點暈。鉆孔、探槽、坑道、樣品、巖性、品位、測斜——這些詞單獨看都認(rèn)識一旦塞進(jìn) Surpac 的數(shù)據(jù)庫結(jié)構(gòu)里很多人第一次打開就懵了為什么一個鉆孔要拆成 collar、survey、geology、sample 好幾張表為什么樣品長度和品位要按區(qū)間存為什么合并數(shù)據(jù)庫時總提示字段不匹配這份教程真正要解決的不是“怎么點按鈕”而是讓你理解 Surpac 地質(zhì)數(shù)據(jù)庫的底層邏輯它本質(zhì)上是一套以“孔口坐標(biāo) 測斜 區(qū)間數(shù)據(jù)”為核心的關(guān)系型數(shù)據(jù)組織方式。你把它想成 Excel 也行但比 Excel 多了空間約束和區(qū)間校驗。適合誰看地質(zhì)勘查、礦山地質(zhì)、測量、三維建模崗位的從業(yè)者尤其是需要把野外原始數(shù)據(jù)整理成能直接用于品位估算和資源量分類的人。這一章先把“數(shù)據(jù)庫在 Surpac 里到底管什么”講清楚后面再動手建庫、導(dǎo)數(shù)據(jù)、排錯。2. 建庫之前先想清楚表結(jié)構(gòu)、字段類型與區(qū)間邏輯2.1 為什么 Surpac 要把一個鉆孔拆成四張表常見做法是collar 表存孔口坐標(biāo)和孔深survey 表存測斜數(shù)據(jù)geology 表存巖性區(qū)間sample 表存樣品區(qū)間和品位。這樣拆的好處是每張表只負(fù)責(zé)一類信息更新巖性不會動到品位補測斜也不會影響樣品。壞處是新手容易在關(guān)聯(lián)時搞錯孔號字段導(dǎo)致“數(shù)據(jù)明明導(dǎo)進(jìn)去了三維視圖里卻看不到孔”。我一般會先畫一張字段對照表把每張表的必填字段列出來再開始建庫。Surpac 對字段類型很敏感孔號用字符串坐標(biāo)和品位用浮點區(qū)間起止深度用浮點巖性代碼用字符串。如果你把孔號存成數(shù)字前導(dǎo)零會丟后面合并數(shù)據(jù)庫時就會對不上。2.2 用命令行方式建一個最小可用的數(shù)據(jù)庫Surpac 的數(shù)據(jù)庫可以通過界面建但更穩(wěn)的方式是用命令腳本建方便復(fù)現(xiàn)。下面是一個最小示例建一個名為geo_demo的數(shù)據(jù)庫包含 collar、survey、geology、sample 四張表。# 進(jìn)入 Surpac 命令行環(huán)境后執(zhí)行 # 創(chuàng)建數(shù)據(jù)庫目錄 create-dir geo_demo # 建 collar 表孔號、東坐標(biāo)、北坐標(biāo)、標(biāo)高、孔深 create-table geo_demo/collar \ hole_id string 16 \ east float \ north float \ rl float \ depth float # 建 survey 表孔號、測點深度、傾角、方位角 create-table geo_demo/survey \ hole_id string 16 \ depth float \ dip float \ azimuth float # 建 geology 表孔號、起始深度、終止深度、巖性代碼 create-table geo_demo/geology \ hole_id string 16 \ from_depth float \ to_depth float \ litho string 8 # 建 sample 表孔號、起始深度、終止深度、品位 create-table geo_demo/sample \ hole_id string 16 \ from_depth float \ to_depth float \ grade float邏輯說明每張表都用hole_id作為關(guān)聯(lián)鍵深度字段統(tǒng)一用浮點巖性代碼留 8 位字符串足夠。參數(shù)說明string 16表示最大 16 個字符float默認(rèn)雙精度。建完后用list-table檢查字段是否和預(yù)期一致尤其是字段順序Surpac 導(dǎo)入時按位置匹配順序錯了會直接串列。2.3 區(qū)間數(shù)據(jù)必須遵守的三條硬規(guī)則第一條from_depth必須小于to_depth且不能為負(fù)。第二條同一孔同一類型的區(qū)間不能重疊比如樣品區(qū)間不能和另一個樣品區(qū)間交叉。第三條區(qū)間不能超出 collar 表里的孔深。這三條在導(dǎo)入時不會全部自動檢查很多錯誤是到三維顯示或品位估算時才暴露。我習(xí)慣在導(dǎo)入前用一段 Python 腳本先做區(qū)間校驗把重疊和超深的數(shù)據(jù)篩出來。# 校驗 sample 表區(qū)間是否重疊或超深 import csv holes {} with open(collar.csv) as f: for row in csv.DictReader(f): holes[row[hole_id]] float(row[depth]) with open(sample.csv) as f: samples list(csv.DictReader(f)) errors [] for s in samples: hid s[hole_id] fd float(s[from_depth]) td float(s[to_depth]) if fd td: errors.append((hid, fd, td, 起止深度倒置)) if td holes.get(hid, 0): errors.append((hid, fd, td, 超出孔深)) # 檢查同孔重疊 from collections import defaultdict by_hole defaultdict(list) for s in samples: by_hole[s[hole_id]].append((float(s[from_depth]), float(s[to_depth]))) for hid, intervals in by_hole.items(): intervals.sort() for i in range(1, len(intervals)): if intervals[i][0] intervals[i-1][1]: errors.append((hid, intervals[i-1], intervals[i], 區(qū)間重疊)) for e in errors: print(e)這段腳本不依賴 Surpac純 Python 就能跑適合在導(dǎo)入前做一輪清洗。參數(shù)說明collar.csv和sample.csv是導(dǎo)出后的文本文件字段名要和腳本里一致。跑完沒有輸出才說明區(qū)間邏輯干凈。3. 把野外數(shù)據(jù)導(dǎo)進(jìn) SurpacCSV 映射、測斜計算與三維顯示3.1 CSV 導(dǎo)入時字段映射的四個關(guān)鍵點Surpac 導(dǎo)入 CSV 時界面會讓你把 CSV 列映射到數(shù)據(jù)庫字段。這里最容易翻車的地方有四個第一分隔符選錯中文環(huán)境下的 CSV 有時用逗號有時用制表符第二孔號列被識別成數(shù)字前導(dǎo)零丟失第三深度列帶了單位“m”導(dǎo)入后變成字符串第四空行沒刪導(dǎo)入后生成空記錄。我一般會先把 CSV 用文本編輯器打開確認(rèn)第一行是字段名沒有多余空格所有深度和品位列都是純數(shù)字。然后在 Surpac 里選擇“按字段名映射”而不是“按位置映射”這樣即使列順序變了也不會串。如果孔號有前導(dǎo)零在 CSV 里就寫成001這種形式或者干脆在導(dǎo)入前把孔號列統(tǒng)一加一個字母前綴。3.2 測斜數(shù)據(jù)怎么算從傾角方位角到三維軌跡測斜表里存的是每個測點的深度、傾角和方位角。Surpac 用這些數(shù)據(jù)計算鉆孔的三維軌跡。常見做法是孔口點已知第一個測點從孔口開始按測段長度和平均傾角方位角推算下一個點。如果測斜數(shù)據(jù)缺失Surpac 會默認(rèn)垂直孔這時候傾斜孔的三維顯示就會明顯偏位。下面是一個用 Python 快速驗算測斜軌跡的片段方便你在導(dǎo)入前檢查數(shù)據(jù)是否合理。import math def compute_trace(collar, surveys): # collar: (east, north, rl) # surveys: list of (depth, dip, azimuth) e, n, rl collar trace [(0.0, e, n, rl)] prev_depth 0.0 for depth, dip, az in surveys: seg depth - prev_depth if seg 0: continue # dip 為負(fù)表示向下 rad_dip math.radians(dip) rad_az math.radians(az) de seg * math.cos(rad_dip) * math.sin(rad_az) dn seg * math.cos(rad_dip) * math.cos(rad_az) drl seg * math.sin(rad_dip) e de n dn rl drl trace.append((depth, e, n, rl)) prev_depth depth return trace # 示例孔口 (1000, 2000, 500)兩個測點 collar (1000.0, 2000.0, 500.0) surveys [(0.0, -90.0, 0.0), (50.0, -80.0, 45.0), (100.0, -75.0, 45.0)] for p in compute_trace(collar, surveys): print(p)邏輯說明dip為負(fù)表示向下鉆進(jìn)azimuth從北順時針。參數(shù)說明seg是測段長度rad_dip和rad_az是弧度轉(zhuǎn)換。算出來的軌跡點可以和你在地質(zhì)圖上看到的孔位做對比偏差大就說明測斜數(shù)據(jù)有問題。3.3 三維顯示前必須檢查的三項第一collar 表的孔口坐標(biāo)是否和測量坐標(biāo)一致尤其是標(biāo)高正負(fù)號。第二survey 表的深度是否從 0 開始如果第一個測點深度不是 0Surpac 會從第一個測點開始畫孔口到第一個測點之間會缺失。第三geology 和 sample 表的區(qū)間是否連續(xù)如果中間有空白段三維顯示時會出現(xiàn)斷孔。我一般會在 Surpac 里先用“顯示鉆孔軌跡”功能看一遍整體形態(tài)再疊加巖性區(qū)間確認(rèn)沒有明顯錯位。4. 避坑與排查Surpac 地質(zhì)數(shù)據(jù)庫最常見的五類翻車4.1 導(dǎo)入后孔號對不上三維視圖里鉆孔消失現(xiàn)象CSV 導(dǎo)入成功數(shù)據(jù)庫里也有記錄但三維視圖里看不到鉆孔。原因collar 表的孔號字段和 survey 表的孔號字段類型不一致一個是字符串一個是數(shù)字關(guān)聯(lián)失敗。解決用list-table檢查兩張表的孔號字段類型統(tǒng)一改成字符串重新導(dǎo)入 survey 表。4.2 樣品區(qū)間重疊導(dǎo)致品位估算報錯現(xiàn)象做品位估算時提示“區(qū)間重疊”或“樣品長度為零”。原因sample 表里同一孔有兩個樣品區(qū)間交叉或者from_depth等于to_depth。解決用第 2 章里的 Python 腳本先篩一遍把重疊區(qū)間拆開或合并確保每個區(qū)間起止深度嚴(yán)格遞增。4.3 測斜數(shù)據(jù)缺失傾斜孔被畫成直孔現(xiàn)象三維視圖里所有鉆孔都是垂直的但實際是斜孔。原因survey 表為空或者測點深度沒有從 0 開始。解決補測斜數(shù)據(jù)確保第一個測點深度為 0傾角和方位角單位是度不是弧度。如果確實沒有測斜至少在 survey 表里加一條深度 0、傾角 -90、方位角 0 的記錄表示垂直孔。4.4 中文巖性代碼導(dǎo)入后變成亂碼現(xiàn)象geology 表里的巖性代碼顯示為問號或亂碼。原因CSV 文件編碼是 GBKSurpac 默認(rèn)按 UTF-8 讀取。解決導(dǎo)入前把 CSV 另存為 UTF-8 無 BOM 格式或者把巖性代碼改成英文縮寫。我一般直接用英文代碼比如GRN表示花崗巖SCH表示片巖省去編碼麻煩。4.5 合并兩個數(shù)據(jù)庫時字段順序不一致現(xiàn)象合并數(shù)據(jù)庫時提示“字段不匹配”或者合并后品位列變成了巖性列。原因兩個數(shù)據(jù)庫的 sample 表字段順序不同Surpac 按位置合并。解決合并前用list-table對比兩張表的字段順序用reorder-field命令調(diào)整成一致。更穩(wěn)的做法是統(tǒng)一用同一套建表腳本從源頭保證字段順序相同。5. 進(jìn)階技巧用 Surpac 命令腳本批量更新品位與生成報告5.1 批量更新樣品品位的命令寫法當(dāng)化驗室返回一批新品位數(shù)據(jù)時手動逐條改容易出錯。常見做法是準(zhǔn)備一個 CSV包含孔號、起始深度、終止深度、新品位然后用 Surpac 命令腳本批量更新。# 批量更新 sample 表品位 # 假設(shè) new_grade.csv 格式hole_id,from_depth,to_depth,grade open-database geo_demo update-table sample from new_grade.csv \ match hole_id,from_depth,to_depth \ set grade close-database邏輯說明match后面跟的是用來定位記錄的字段組合set后面是要更新的字段。參數(shù)說明CSV 里的字段名要和數(shù)據(jù)庫字段名一致深度匹配允許微小誤差時加tolerance 0.01。更新前先備份數(shù)據(jù)庫避免匹配錯誤導(dǎo)致品位被覆蓋。5.2 用報告模板輸出品位區(qū)間統(tǒng)計Surpac 的報告功能可以按巖性、按孔、按區(qū)間統(tǒng)計品位。我一般會建一個報告模板包含孔號、巖性、樣品數(shù)、平均品位、最大品位、最小品位。模板建好后每次更新數(shù)據(jù)只需重新運行報告不用重新配置。字段含義示例hole_id孔號ZK001litho巖性代碼GRNsample_count樣品數(shù)12avg_grade平均品位1.85max_grade最大品位4.20min_grade最小品位0.30報告生成后導(dǎo)出為 CSV方便在 Excel 里做進(jìn)一步統(tǒng)計。注意平均品位要按樣品長度加權(quán)不能直接算術(shù)平均否則長樣品和短樣品權(quán)重一樣結(jié)果會偏。5.3 一個我常用的習(xí)慣每次導(dǎo)入前先跑校驗?zāi)_本這些年我養(yǎng)成了一個習(xí)慣不管數(shù)據(jù)來自哪里導(dǎo)入 Surpac 之前先跑一遍校驗?zāi)_本檢查孔號重復(fù)、區(qū)間重疊、深度超限、品位為空。腳本不長但能省掉后面大量排查時間。有一次一批樣品數(shù)據(jù)里混了三個孔號重復(fù)的記錄導(dǎo)入后三維視圖里多出三條短孔查了半天才發(fā)現(xiàn)是 CSV 里復(fù)制粘貼多出來的。從那以后校驗?zāi)_本就成了固定動作。希望幫到你。本文還有配套的精品資源點擊獲取 Surpac地質(zhì)數(shù)據(jù)庫實戰(zhàn)：建庫、導(dǎo)入與避坑指南
簡介這份PDF文檔是Surpac Vision地質(zhì)數(shù)據(jù)庫的官方培訓(xùn)教程面向礦山地質(zhì)、測量及采礦工程領(lǐng)域的技術(shù)人員與初學(xué)者幫助其系統(tǒng)掌握地質(zhì)數(shù)據(jù)的建庫、管理與三維可視化方法。資源包內(nèi)僅含1個PDF文件大小約2.32MB內(nèi)容按章節(jié)組織涵蓋數(shù)據(jù)庫介紹、結(jié)構(gòu)解析、創(chuàng)建流程、編輯瀏覽管理、三維空間顯示、地質(zhì)解譯與剖面品位計算、組合樣品、地質(zhì)統(tǒng)計與分析、線文件剖面輸出及信息提取報告等模塊并附有3D轉(zhuǎn)換實例與多邊形資源評估兩個附件。教程從表與字段分析、空庫創(chuàng)建、數(shù)據(jù)導(dǎo)入驗證到映射數(shù)據(jù)庫逐步展開同時講解鉆孔顯示風(fēng)格、隱藏刷新、勘探線定義與礦體圈定等實操要點便于讀者對照軟件界面完成練習(xí)。目前已有221人學(xué)習(xí)適合作為Surpac入門與日常功能查閱的參考手冊。1. 從一份 PDF 說起Surpac 地質(zhì)數(shù)據(jù)庫到底在管什么如果你手里正拿著一份名為“2-surpac地質(zhì)數(shù)據(jù)庫教程.pdf”的資料大概率你已經(jīng)被 Surpac 里那套“數(shù)據(jù)庫”概念繞得有點暈。鉆孔、探槽、坑道、樣品、巖性、品位、測斜——這些詞單獨看都認(rèn)識一旦塞進(jìn) Surpac 的數(shù)據(jù)庫結(jié)構(gòu)里很多人第一次打開就懵了為什么一個鉆孔要拆成 collar、survey、geology、sample 好幾張表為什么樣品長度和品位要按區(qū)間存為什么合并數(shù)據(jù)庫時總提示字段不匹配這份教程真正要解決的不是“怎么點按鈕”而是讓你理解 Surpac 地質(zhì)數(shù)據(jù)庫的底層邏輯它本質(zhì)上是一套以“孔口坐標(biāo) 測斜 區(qū)間數(shù)據(jù)”為核心的關(guān)系型數(shù)據(jù)組織方式。你把它想成 Excel 也行但比 Excel 多了空間約束和區(qū)間校驗。適合誰看地質(zhì)勘查、礦山地質(zhì)、測量、三維建模崗位的從業(yè)者尤其是需要把野外原始數(shù)據(jù)整理成能直接用于品位估算和資源量分類的人。這一章先把“數(shù)據(jù)庫在 Surpac 里到底管什么”講清楚后面再動手建庫、導(dǎo)數(shù)據(jù)、排錯。2. 建庫之前先想清楚表結(jié)構(gòu)、字段類型與區(qū)間邏輯2.1 為什么 Surpac 要把一個鉆孔拆成四張表常見做法是collar 表存孔口坐標(biāo)和孔深survey 表存測斜數(shù)據(jù)geology 表存巖性區(qū)間sample 表存樣品區(qū)間和品位。這樣拆的好處是每張表只負(fù)責(zé)一類信息更新巖性不會動到品位補測斜也不會影響樣品。壞處是新手容易在關(guān)聯(lián)時搞錯孔號字段導(dǎo)致“數(shù)據(jù)明明導(dǎo)進(jìn)去了三維視圖里卻看不到孔”。我一般會先畫一張字段對照表把每張表的必填字段列出來再開始建庫。Surpac 對字段類型很敏感孔號用字符串坐標(biāo)和品位用浮點區(qū)間起止深度用浮點巖性代碼用字符串。如果你把孔號存成數(shù)字前導(dǎo)零會丟后面合并數(shù)據(jù)庫時就會對不上。2.2 用命令行方式建一個最小可用的數(shù)據(jù)庫Surpac 的數(shù)據(jù)庫可以通過界面建但更穩(wěn)的方式是用命令腳本建方便復(fù)現(xiàn)。下面是一個最小示例建一個名為geo_demo的數(shù)據(jù)庫包含 collar、survey、geology、sample 四張表。# 進(jìn)入 Surpac 命令行環(huán)境后執(zhí)行 # 創(chuàng)建數(shù)據(jù)庫目錄 create-dir geo_demo # 建 collar 表孔號、東坐標(biāo)、北坐標(biāo)、標(biāo)高、孔深 create-table geo_demo/collar \ hole_id string 16 \ east float \ north float \ rl float \ depth float # 建 survey 表孔號、測點深度、傾角、方位角 create-table geo_demo/survey \ hole_id string 16 \ depth float \ dip float \ azimuth float # 建 geology 表孔號、起始深度、終止深度、巖性代碼 create-table geo_demo/geology \ hole_id string 16 \ from_depth float \ to_depth float \ litho string 8 # 建 sample 表孔號、起始深度、終止深度、品位 create-table geo_demo/sample \ hole_id string 16 \ from_depth float \ to_depth float \ grade float邏輯說明每張表都用hole_id作為關(guān)聯(lián)鍵深度字段統(tǒng)一用浮點巖性代碼留 8 位字符串足夠。參數(shù)說明string 16表示最大 16 個字符float默認(rèn)雙精度。建完后用list-table檢查字段是否和預(yù)期一致尤其是字段順序Surpac 導(dǎo)入時按位置匹配順序錯了會直接串列。2.3 區(qū)間數(shù)據(jù)必須遵守的三條硬規(guī)則第一條from_depth必須小于to_depth且不能為負(fù)。第二條同一孔同一類型的區(qū)間不能重疊比如樣品區(qū)間不能和另一個樣品區(qū)間交叉。第三條區(qū)間不能超出 collar 表里的孔深。這三條在導(dǎo)入時不會全部自動檢查很多錯誤是到三維顯示或品位估算時才暴露。我習(xí)慣在導(dǎo)入前用一段 Python 腳本先做區(qū)間校驗把重疊和超深的數(shù)據(jù)篩出來。# 校驗 sample 表區(qū)間是否重疊或超深 import csv holes {} with open(collar.csv) as f: for row in csv.DictReader(f): holes[row[hole_id]] float(row[depth]) with open(sample.csv) as f: samples list(csv.DictReader(f)) errors [] for s in samples: hid s[hole_id] fd float(s[from_depth]) td float(s[to_depth]) if fd td: errors.append((hid, fd, td, 起止深度倒置)) if td holes.get(hid, 0): errors.append((hid, fd, td, 超出孔深)) # 檢查同孔重疊 from collections import defaultdict by_hole defaultdict(list) for s in samples: by_hole[s[hole_id]].append((float(s[from_depth]), float(s[to_depth]))) for hid, intervals in by_hole.items(): intervals.sort() for i in range(1, len(intervals)): if intervals[i][0] intervals[i-1][1]: errors.append((hid, intervals[i-1], intervals[i], 區(qū)間重疊)) for e in errors: print(e)這段腳本不依賴 Surpac純 Python 就能跑適合在導(dǎo)入前做一輪清洗。參數(shù)說明collar.csv和sample.csv是導(dǎo)出后的文本文件字段名要和腳本里一致。跑完沒有輸出才說明區(qū)間邏輯干凈。3. 把野外數(shù)據(jù)導(dǎo)進(jìn) SurpacCSV 映射、測斜計算與三維顯示3.1 CSV 導(dǎo)入時字段映射的四個關(guān)鍵點Surpac 導(dǎo)入 CSV 時界面會讓你把 CSV 列映射到數(shù)據(jù)庫字段。這里最容易翻車的地方有四個第一分隔符選錯中文環(huán)境下的 CSV 有時用逗號有時用制表符第二孔號列被識別成數(shù)字前導(dǎo)零丟失第三深度列帶了單位“m”導(dǎo)入后變成字符串第四空行沒刪導(dǎo)入后生成空記錄。我一般會先把 CSV 用文本編輯器打開確認(rèn)第一行是字段名沒有多余空格所有深度和品位列都是純數(shù)字。然后在 Surpac 里選擇“按字段名映射”而不是“按位置映射”這樣即使列順序變了也不會串。如果孔號有前導(dǎo)零在 CSV 里就寫成001這種形式或者干脆在導(dǎo)入前把孔號列統(tǒng)一加一個字母前綴。3.2 測斜數(shù)據(jù)怎么算從傾角方位角到三維軌跡測斜表里存的是每個測點的深度、傾角和方位角。Surpac 用這些數(shù)據(jù)計算鉆孔的三維軌跡。常見做法是孔口點已知第一個測點從孔口開始按測段長度和平均傾角方位角推算下一個點。如果測斜數(shù)據(jù)缺失Surpac 會默認(rèn)垂直孔這時候傾斜孔的三維顯示就會明顯偏位。下面是一個用 Python 快速驗算測斜軌跡的片段方便你在導(dǎo)入前檢查數(shù)據(jù)是否合理。import math def compute_trace(collar, surveys): # collar: (east, north, rl) # surveys: list of (depth, dip, azimuth) e, n, rl collar trace [(0.0, e, n, rl)] prev_depth 0.0 for depth, dip, az in surveys: seg depth - prev_depth if seg 0: continue # dip 為負(fù)表示向下 rad_dip math.radians(dip) rad_az math.radians(az) de seg * math.cos(rad_dip) * math.sin(rad_az) dn seg * math.cos(rad_dip) * math.cos(rad_az) drl seg * math.sin(rad_dip) e de n dn rl drl trace.append((depth, e, n, rl)) prev_depth depth return trace # 示例孔口 (1000, 2000, 500)兩個測點 collar (1000.0, 2000.0, 500.0) surveys [(0.0, -90.0, 0.0), (50.0, -80.0, 45.0), (100.0, -75.0, 45.0)] for p in compute_trace(collar, surveys): print(p)邏輯說明dip為負(fù)表示向下鉆進(jìn)azimuth從北順時針。參數(shù)說明seg是測段長度rad_dip和rad_az是弧度轉(zhuǎn)換。算出來的軌跡點可以和你在地質(zhì)圖上看到的孔位做對比偏差大就說明測斜數(shù)據(jù)有問題。3.3 三維顯示前必須檢查的三項第一collar 表的孔口坐標(biāo)是否和測量坐標(biāo)一致尤其是標(biāo)高正負(fù)號。第二survey 表的深度是否從 0 開始如果第一個測點深度不是 0Surpac 會從第一個測點開始畫孔口到第一個測點之間會缺失。第三geology 和 sample 表的區(qū)間是否連續(xù)如果中間有空白段三維顯示時會出現(xiàn)斷孔。我一般會在 Surpac 里先用“顯示鉆孔軌跡”功能看一遍整體形態(tài)再疊加巖性區(qū)間確認(rèn)沒有明顯錯位。4. 避坑與排查Surpac 地質(zhì)數(shù)據(jù)庫最常見的五類翻車4.1 導(dǎo)入后孔號對不上三維視圖里鉆孔消失現(xiàn)象CSV 導(dǎo)入成功數(shù)據(jù)庫里也有記錄但三維視圖里看不到鉆孔。原因collar 表的孔號字段和 survey 表的孔號字段類型不一致一個是字符串一個是數(shù)字關(guān)聯(lián)失敗。解決用list-table檢查兩張表的孔號字段類型統(tǒng)一改成字符串重新導(dǎo)入 survey 表。4.2 樣品區(qū)間重疊導(dǎo)致品位估算報錯現(xiàn)象做品位估算時提示“區(qū)間重疊”或“樣品長度為零”。原因sample 表里同一孔有兩個樣品區(qū)間交叉或者from_depth等于to_depth。解決用第 2 章里的 Python 腳本先篩一遍把重疊區(qū)間拆開或合并確保每個區(qū)間起止深度嚴(yán)格遞增。4.3 測斜數(shù)據(jù)缺失傾斜孔被畫成直孔現(xiàn)象三維視圖里所有鉆孔都是垂直的但實際是斜孔。原因survey 表為空或者測點深度沒有從 0 開始。解決補測斜數(shù)據(jù)確保第一個測點深度為 0傾角和方位角單位是度不是弧度。如果確實沒有測斜至少在 survey 表里加一條深度 0、傾角 -90、方位角 0 的記錄表示垂直孔。4.4 中文巖性代碼導(dǎo)入后變成亂碼現(xiàn)象geology 表里的巖性代碼顯示為問號或亂碼。原因CSV 文件編碼是 GBKSurpac 默認(rèn)按 UTF-8 讀取。解決導(dǎo)入前把 CSV 另存為 UTF-8 無 BOM 格式或者把巖性代碼改成英文縮寫。我一般直接用英文代碼比如GRN表示花崗巖SCH表示片巖省去編碼麻煩。4.5 合并兩個數(shù)據(jù)庫時字段順序不一致現(xiàn)象合并數(shù)據(jù)庫時提示“字段不匹配”或者合并后品位列變成了巖性列。原因兩個數(shù)據(jù)庫的 sample 表字段順序不同Surpac 按位置合并。解決合并前用list-table對比兩張表的字段順序用reorder-field命令調(diào)整成一致。更穩(wěn)的做法是統(tǒng)一用同一套建表腳本從源頭保證字段順序相同。5. 進(jìn)階技巧用 Surpac 命令腳本批量更新品位與生成報告5.1 批量更新樣品品位的命令寫法當(dāng)化驗室返回一批新品位數(shù)據(jù)時手動逐條改容易出錯。常見做法是準(zhǔn)備一個 CSV包含孔號、起始深度、終止深度、新品位然后用 Surpac 命令腳本批量更新。# 批量更新 sample 表品位 # 假設(shè) new_grade.csv 格式hole_id,from_depth,to_depth,grade open-database geo_demo update-table sample from new_grade.csv \ match hole_id,from_depth,to_depth \ set grade close-database邏輯說明match后面跟的是用來定位記錄的字段組合set后面是要更新的字段。參數(shù)說明CSV 里的字段名要和數(shù)據(jù)庫字段名一致深度匹配允許微小誤差時加tolerance 0.01。更新前先備份數(shù)據(jù)庫避免匹配錯誤導(dǎo)致品位被覆蓋。5.2 用報告模板輸出品位區(qū)間統(tǒng)計Surpac 的報告功能可以按巖性、按孔、按區(qū)間統(tǒng)計品位。我一般會建一個報告模板包含孔號、巖性、樣品數(shù)、平均品位、最大品位、最小品位。模板建好后每次更新數(shù)據(jù)只需重新運行報告不用重新配置。字段含義示例hole_id孔號ZK001litho巖性代碼GRNsample_count樣品數(shù)12avg_grade平均品位1.85max_grade最大品位4.20min_grade最小品位0.30報告生成后導(dǎo)出為 CSV方便在 Excel 里做進(jìn)一步統(tǒng)計。注意平均品位要按樣品長度加權(quán)不能直接算術(shù)平均否則長樣品和短樣品權(quán)重一樣結(jié)果會偏。5.3 一個我常用的習(xí)慣每次導(dǎo)入前先跑校驗?zāi)_本這些年我養(yǎng)成了一個習(xí)慣不管數(shù)據(jù)來自哪里導(dǎo)入 Surpac 之前先跑一遍校驗?zāi)_本檢查孔號重復(fù)、區(qū)間重疊、深度超限、品位為空。腳本不長但能省掉后面大量排查時間。有一次一批樣品數(shù)據(jù)里混了三個孔號重復(fù)的記錄導(dǎo)入后三維視圖里多出三條短孔查了半天才發(fā)現(xiàn)是 CSV 里復(fù)制粘貼多出來的。從那以后校驗?zāi)_本就成了固定動作。希望幫到你。本文還有配套的精品資源點擊獲取