柳建平

(山西省煤炭地質(zhì)115勘查院，山西大同 037002)

煤田地質(zhì)勘查需要制作大量的圖件，其中不乏有各種各樣的等值線圖，如煤層底板等高線圖、煤層厚度等值線圖、巖溶等水位線圖、突水性分區(qū)圖等。這些圖件大都需要利用鉆孔測(cè)量的各種數(shù)據(jù)勾繪等值線，以達(dá)到反映某種趨勢(shì)的目的。在沒有掌握Mapgis此項(xiàng)功能時(shí)大都是手工勾繪，計(jì)算機(jī)掃描之后利用Mapgis描圖，工作繁瑣且效率低下。利用Mapgis的數(shù)字地面模型子系統(tǒng)即可輕松實(shí)現(xiàn)等值線的自動(dòng)勾繪，且能反映出它們的自然趨勢(shì)，方法科學(xué)簡(jiǎn)單，效率高，尤其是在有大量數(shù)據(jù)的情況下更顯現(xiàn)出該方法的優(yōu)勢(shì)。該方法已成為地質(zhì)工作者快速完成數(shù)據(jù)圖形化工作的法寶。以下就對(duì)該方法的實(shí)際操作做詳細(xì)的介紹。

1 投影變換

將鉆孔的編號(hào)、坐標(biāo)和測(cè)量數(shù)據(jù)分別輸入Excel中，由于Mapgis投影變換系統(tǒng)不能識(shí)別Excel文件，故需另存為文本文檔［1］。在沒有坐標(biāo)文件時(shí)，也可利用Mapgis文件轉(zhuǎn)換或SECTION將點(diǎn)文件屬性輸出到表格中。文章以某井田勘探區(qū)的鉆孔突水系數(shù)為例，繪制突水系數(shù)等值線圖。1)投影變換:Mapgis6.7平臺(tái)下利用實(shí)用服務(wù)—投影變換—用戶文件投影變換打開原始數(shù)據(jù)文本文檔，見圖1。點(diǎn)擊指定數(shù)據(jù)起始位置的第二行數(shù)據(jù)，使得Mapgis能夠自動(dòng)讀取原始數(shù)據(jù)的坐標(biāo)［2］;設(shè)置用戶文件選項(xiàng)按指定分隔符、保留單列、生成點(diǎn);由于Mapgis中X，Y坐標(biāo)與實(shí)際坐標(biāo)相反，故X為第三列數(shù)據(jù)，Y為第二列數(shù)據(jù)。2)設(shè)置分隔符:點(diǎn)擊設(shè)置分隔符，將分隔符號(hào)中Tab鍵打勾，屬性名稱所在行選擇第二行，確定，如圖2所示。3)設(shè)置點(diǎn)圖元參數(shù):自由選擇子圖號(hào)及寬度大小，與數(shù)據(jù)坐標(biāo)的尺度要協(xié)調(diào)，以方便觀察為準(zhǔn)。4)設(shè)置完成后點(diǎn)擊不需要投影，數(shù)據(jù)生成，得到了鉆孔點(diǎn)文件。在Mapgis輸入編輯中編輯點(diǎn)文件，根據(jù)屬性標(biāo)注將突水系數(shù)標(biāo)注到點(diǎn)文件中，得到了含有鉆孔位置及相應(yīng)點(diǎn)的突水系數(shù)值的點(diǎn)文件，以備成圖使用。

圖1 投影變換操作界面（一）

圖2 投影變換操作界面（二）

2 空間分析

在Mapgis空間分析的數(shù)字地面模型(DTM)子系統(tǒng)中，進(jìn)入數(shù)字高程模型(Grd)，首先將離散數(shù)據(jù)網(wǎng)格化，然后進(jìn)行平面等值線繪制。下面做詳細(xì)介紹。

1)離散數(shù)據(jù)網(wǎng)格化:Grd模型(數(shù)字高程模型)分析是建立在網(wǎng)格化的數(shù)據(jù)的基礎(chǔ)上的，所以對(duì)于非網(wǎng)格化數(shù)據(jù)必須先網(wǎng)格化才能進(jìn)行有關(guān)分析［3］。進(jìn)入空間分析—DTM分析—Grd模型—離散數(shù)據(jù)網(wǎng)格化，選擇原始數(shù)據(jù)文本文檔打開后得到如圖3所示對(duì)話框。將X，Y，Z分別與文件中的數(shù)據(jù)相對(duì)應(yīng):X選擇文件中的第三列，Y選擇第二列，Z選擇第四列。系統(tǒng)提供了四種網(wǎng)格化方法，分別為距離冪函數(shù)反比加權(quán)函數(shù)、泛克立格法網(wǎng)格化、稠密數(shù)據(jù)中值選取網(wǎng)格化和稠密數(shù)據(jù)高斯距離權(quán)網(wǎng)格化，不同的方法適合不同的數(shù)據(jù)模型［4］。本文采用泛克立格網(wǎng)格化法。設(shè)置后對(duì)文件換名保存到指定的文件中，對(duì)話框如圖4所示。設(shè)置完成后點(diǎn)擊確定，得到一個(gè)網(wǎng)格化的GRD文件。

2)平面等值線圖繪制:Grd模型可以繪制三種圖件，分別是網(wǎng)格化立體圖、平面等值線圖和彩色等值立體圖。本文采用平面等值線圖，點(diǎn)擊Grd模型中的平面等值線圖繪制，選擇網(wǎng)格化的GRD文件，得到如圖5所示對(duì)話框。

勾選等值線套區(qū)、繪制色階、保留邊界線、等值線光滑處理，光滑度依個(gè)人習(xí)慣而定，高程度繪制的圖更加圓滑美觀，但有時(shí)可能會(huì)與部分原始數(shù)據(jù)有所偏離，所以推薦使用中程度。色階(顏色)也可以依個(gè)人習(xí)慣喜好自由設(shè)置，以方便觀察為準(zhǔn)。制圖幅面選擇原始數(shù)據(jù)范圍，若采用自動(dòng)檢測(cè)設(shè)置，則所得的圖不能與原始數(shù)據(jù)坐標(biāo)網(wǎng)格相吻合，也就沒有了實(shí)際意義。注意參數(shù)中頻度代表標(biāo)注的稀密程度，此處應(yīng)改為零，得到如圖6所示對(duì)話框。等值線層值的設(shè)置以實(shí)際情況而定，一般設(shè)置成整數(shù)以方便制作后圖件的觀察使用。點(diǎn)擊等值線層值，將起始Z值、終止Z值和步長(zhǎng)增量全部改為與原始數(shù)據(jù)小數(shù)點(diǎn)位數(shù)一致的值(這里原始數(shù)據(jù)精確到小數(shù)點(diǎn)后兩位)，如圖7所示。設(shè)置完成后點(diǎn)擊確定，得到了等值線草圖(如圖8所示)，將點(diǎn)線區(qū)文件保存于指定文件夾中。

圖3 網(wǎng)格化設(shè)置對(duì)話框（一）

圖4 網(wǎng)格化設(shè)置對(duì)話框（二）

圖5 參數(shù)設(shè)置對(duì)話框（一）

圖6 參數(shù)設(shè)置對(duì)話框（二）

圖7 等值線層設(shè)置對(duì)話框

圖8 等值線草圖

3 圖形處理

上述步驟所得的等值線草圖中，等值線值在圖例色階中，還需要手動(dòng)將等值線值輸入圖中，并將前面的鉆孔點(diǎn)文件添加到該圖中，得到突水系數(shù)等值線圖。根據(jù)突水系數(shù)值進(jìn)行分區(qū)，文章用0.06和0.10兩個(gè)分界線進(jìn)行分區(qū)，用以反映煤層頂?shù)装灏l(fā)生充水的危險(xiǎn)程度。

1)輸入編輯添加等值線圖點(diǎn)線區(qū)文件和鉆孔點(diǎn)文件。根據(jù)色階值編輯等值線值，將相應(yīng)的值編輯到對(duì)應(yīng)的等值線上;而后將色階、邊框等點(diǎn)線區(qū)文件刪除，保存后將點(diǎn)線文件添加背景圖框，即得到突水系數(shù)等值線圖，如圖9所示。

2)利用等值線繪制突水分區(qū)圖，以0.06和0.10兩條等值線為界進(jìn)行分區(qū)，輸入弧段造區(qū)或者重新描線利用拓?fù)渲亟ㄔ靺^(qū)，得到突水性分區(qū)圖，其中突水系數(shù)大于0.10及與其聯(lián)通斷層附近區(qū)域?yàn)榘l(fā)生充水的危險(xiǎn)地段，如圖10所示。

圖9 突水系數(shù)等值線圖

圖10 突水性分區(qū)圖

4 結(jié)語(yǔ)

應(yīng)用Mapgis數(shù)字地面模型子系統(tǒng)繪制的圖件可以很好的與原始數(shù)據(jù)套合且比例尺相同，應(yīng)用起來非常方便快捷，而且圖件準(zhǔn)確、整潔，是等值線圖件繪制的不二方法。該方法不僅能夠應(yīng)用于煤田地質(zhì)勘查中，在化探工作中也得到非常廣闊的應(yīng)用，社會(huì)經(jīng)濟(jì)效益顯著。

［1］ 蘇新平.MAPGIS6.7環(huán)境下化探數(shù)據(jù)預(yù)處理方法［J］.煤炭技術(shù)，2012，31(7):128-129.

［2］ 中地?cái)?shù)碼.MAPGIS實(shí)用教程［Z］.2003:187-202.

［3］ 陸可喜，石文學(xué)，劉繼朝.簡(jiǎn)述一種基于MAPGIS和Office工具的繪圖方法［J］.地下水，2008，30(3):113-115.

［4］ 王聲喜.MAPGIS地質(zhì)制圖基礎(chǔ)講義［Z］.2010:191-196.