閆科偉，牛自禮

（1.山西水利職業(yè)技術(shù)學(xué) 交通工程，山西 運(yùn)城 044000；2.平頂山市公路交通勘察設(shè)計(jì)院，河南 平頂山 467000）

礦山勘探工程安全生產(chǎn)受到多方面因素的制約，其中礦井水害是最難治理的地質(zhì)災(zāi)害之一，礦井突水引發(fā)的事故會(huì)給礦區(qū)開(kāi)發(fā)造成嚴(yán)重負(fù)面影響，甚至還會(huì)危害勘探人員的生命安全。因此勘探工作前，首先需要探明礦山水文地質(zhì)結(jié)構(gòu)，但是由于多數(shù)礦山水文地質(zhì)條件復(fù)雜，對(duì)準(zhǔn)確預(yù)測(cè)礦井水害帶來(lái)巨大干擾，增大了礦山開(kāi)采的工作難度[1]。地下水模擬軟件（GMS），為三維可視化地質(zhì)建模提供了極大的幫助，利用GMS建立的地質(zhì)模型更為直接，可以有效顯示含水層和隔水層等水文地質(zhì)結(jié)構(gòu)的空間分布和狀態(tài)，對(duì)了解礦山水文地質(zhì)結(jié)構(gòu)和深入分析地下水內(nèi)在規(guī)律具有重要作用。因此本文基于GMS軟件，對(duì)礦山水文地質(zhì)結(jié)構(gòu)三維可視化系統(tǒng)進(jìn)行設(shè)計(jì)，對(duì)查明礦區(qū)水文地質(zhì)條件，準(zhǔn)確判斷礦坑涌水量提供技術(shù)支持，對(duì)提高礦山開(kāi)采的安全性和經(jīng)濟(jì)性起到一定促進(jìn)作用。

1 礦山水文地質(zhì)結(jié)構(gòu)監(jiān)測(cè)硬件

1.1 選擇芯片

由于是對(duì)野外環(huán)境中礦山水文地質(zhì)結(jié)構(gòu)的監(jiān)測(cè)，通常規(guī)模較大和節(jié)點(diǎn)數(shù)量較多，因此系統(tǒng)設(shè)計(jì)要盡量控制硬件成本。從經(jīng)濟(jì)環(huán)保的角度考慮，選擇的芯片需要選擇低能耗且依靠太陽(yáng)能供電工作，減少人工部署工作量。本文選擇SOC芯片，將必要電路集成到單一芯片上，SOC芯片集成度高，具有高性能低成本的特點(diǎn)，滿(mǎn)足本文系統(tǒng)的硬件要求。芯片在接收到信號(hào)后，對(duì)信號(hào)進(jìn)行過(guò)濾處理，將有效信息存儲(chǔ)進(jìn)寄存器中；在向外發(fā)送信號(hào)時(shí)，先將信號(hào)放入模擬數(shù)字轉(zhuǎn)換器，再經(jīng)天線發(fā)送，完成整個(gè)信號(hào)傳遞的過(guò)程。

1.2 部署傳感器

為適應(yīng)礦山環(huán)境，部署采集數(shù)據(jù)的傳感器要滿(mǎn)足大范圍無(wú)線監(jiān)測(cè)的要求，本文選擇光纖傳感器，將其部署于礦山水文地質(zhì)鉆孔中，傳感器的傾角應(yīng)與鉆孔保持一致。通過(guò)傳感器實(shí)時(shí)監(jiān)測(cè)鉆孔內(nèi)信息數(shù)據(jù)的變化，不僅可以查明含水層和滲水層情況，還確定未含水地層變形狀況，借此判斷地層裂隙發(fā)育程度，分析該地層是否與含水層產(chǎn)生聯(lián)系，總結(jié)水文地質(zhì)結(jié)構(gòu)變化規(guī)律。

2 基于GMS的三維可視化系統(tǒng)軟件

2.1 構(gòu)建礦山水文地質(zhì)結(jié)構(gòu)模型

實(shí)際礦山地下水層結(jié)構(gòu)往往十分復(fù)雜，因此需要建立水文地質(zhì)結(jié)構(gòu)模型。建立模型的主要原則是對(duì)水文條件進(jìn)行概念化，保證可以再現(xiàn)水流動(dòng)態(tài)。首先對(duì)含水層實(shí)際邊界條件進(jìn)行概念化計(jì)算，通常情況下，含水層上部為交換邊界，接收降雨和河流等的補(bǔ)給，將底部不透水基巖設(shè)定為模型隔水邊界。將模型水頭邊界設(shè)定為第一類(lèi)邊界，計(jì)算公式如下：

其中，Γ2表示二類(lèi)邊界；k表示三維空間的滲透系數(shù)；表示邊界法向量；h2(x,y,z)表示二類(lèi)邊界上已知流量函數(shù)。將模型混合邊界設(shè)定為第三類(lèi)邊界，計(jì)算公式如下：

其中，Γ3表示三類(lèi)邊界；h3(x,y,z)表示三類(lèi)邊界上流量函數(shù)；C表示邊界條件函數(shù)；h表示水頭邊界取值；h0表示天然邊界取值。

2.2 設(shè)計(jì)水文地質(zhì)結(jié)構(gòu)數(shù)據(jù)庫(kù)

為建立水文地質(zhì)結(jié)構(gòu)數(shù)據(jù)庫(kù)，首先依據(jù)數(shù)據(jù)來(lái)源劃分類(lèi)別，本文系統(tǒng)數(shù)據(jù)庫(kù)將數(shù)據(jù)劃分為含水層特征、水文地質(zhì)鉆孔和水文地質(zhì)概念模型邊界三大類(lèi)，再分為不同子類(lèi)，各類(lèi)數(shù)據(jù)具有不同信息屬性。其中水文地質(zhì)鉆孔數(shù)據(jù)表是三維可視化系統(tǒng)的主要依據(jù)，描述水文地質(zhì)結(jié)構(gòu)屬性特征，本文設(shè)計(jì)鉆孔基礎(chǔ)信息表見(jiàn)表1，地層分層信息表見(jiàn)表2，含水層信息表見(jiàn)表3。

表1 水文地質(zhì)鉆孔基礎(chǔ)信息表

表2 鉆孔地層分層信息表

表3 鉆孔含水層層段信息表

水文地質(zhì)鉆孔的三個(gè)信息表，通過(guò)鉆孔統(tǒng)一編號(hào)關(guān)聯(lián)起來(lái)，得到鉆孔空間結(jié)構(gòu)關(guān)系，為三維可視化系統(tǒng)提供信息數(shù)據(jù)支持。

2.3 選取空間插值算法

在三維可視化建模過(guò)程中，需要對(duì)采集的原始數(shù)據(jù)使用插值算法，對(duì)空間采樣不充足和缺乏的區(qū)域進(jìn)行適當(dāng)修正，來(lái)提高建模完整性和精確度。由于地質(zhì)結(jié)構(gòu)的復(fù)雜和多變性，選取的插值算法會(huì)對(duì)三維可視化精度造成一定影響，本文選擇空間插值算法對(duì)模型進(jìn)行擬合，提高系統(tǒng)準(zhǔn)確度。假定平面上分布散亂點(diǎn)p(x,y,z)，其坐標(biāo)為(xi,yi)，屬性值為Zi(i= 1,2,...,n)，根據(jù)周?chē)y點(diǎn)屬性值計(jì)算p點(diǎn)的值?？臻g插值算法可以表示為：

其中，di表示散亂點(diǎn)到p點(diǎn)的距離；k∈[0 ,2]?？臻g插值算法主要是利用未知點(diǎn)屬性進(jìn)行加權(quán)平均，使插值生成曲面較為平滑。

3 實(shí)驗(yàn)結(jié)果與分析

3.1 系統(tǒng)性能測(cè)試

為驗(yàn)證本文設(shè)計(jì)三維可視化系統(tǒng)的可靠性，以實(shí)際某礦山水文地質(zhì)勘察資料為基礎(chǔ)，對(duì)系統(tǒng)性能進(jìn)行測(cè)試。首先建立了該礦山水文地質(zhì)結(jié)構(gòu)數(shù)據(jù)庫(kù)，利用實(shí)際資料構(gòu)建三維模型，對(duì)其可視化，在計(jì)算機(jī)環(huán)境中展示水文地質(zhì)結(jié)構(gòu)及其空間分布。根據(jù)實(shí)際勘探資料中含水層介質(zhì)分布特征分析，得到水文地質(zhì)參數(shù)取值情況，如表4所示。

表4 水文地質(zhì)參數(shù)實(shí)際值

利用本文系統(tǒng)對(duì)研究區(qū)模型進(jìn)行識(shí)別，計(jì)算出滲透系數(shù)和給水度見(jiàn)表5。

表5 水文地質(zhì)參數(shù)識(shí)別值

通過(guò)參數(shù)實(shí)際值和系統(tǒng)識(shí)別值的對(duì)比，可以發(fā)現(xiàn)偏差相對(duì)較小，因此說(shuō)明含水層擬合情況較好，三維可視化結(jié)果較為準(zhǔn)確。

3.2 對(duì)比測(cè)試

為檢驗(yàn)本文設(shè)計(jì)的三維可視化系統(tǒng)實(shí)際應(yīng)用價(jià)值，將本文設(shè)計(jì)系統(tǒng)作為實(shí)驗(yàn)組，將使用GOCAD和GIS軟件的系統(tǒng)作為對(duì)照組，進(jìn)行對(duì)比測(cè)試，比較系統(tǒng)的水均衡擬合效果，對(duì)比結(jié)果見(jiàn)表6。

表6 不同系統(tǒng)對(duì)比測(cè)試結(jié)果

通過(guò)表6對(duì)比結(jié)果可知，本文系統(tǒng)的水均衡效果擬合較好，可應(yīng)用于實(shí)際礦山水文地質(zhì)勘查中。

4 結(jié)束語(yǔ)

本文基于GMS對(duì)礦山水文地質(zhì)結(jié)構(gòu)三維可視化系統(tǒng)進(jìn)行設(shè)計(jì)，為探明水文地質(zhì)結(jié)構(gòu)和含水層特征，提供可視化平臺(tái)。但本文研究是對(duì)穩(wěn)定水流狀態(tài)下的模擬預(yù)測(cè)，還存在一定局限性，未來(lái)系統(tǒng)設(shè)計(jì)還應(yīng)考慮非穩(wěn)定流模擬；水文地質(zhì)結(jié)構(gòu)的三維可視化模型的精度應(yīng)進(jìn)一步提高；本文主要對(duì)地下水的動(dòng)態(tài)信息進(jìn)行描述，后續(xù)還應(yīng)結(jié)合地層和構(gòu)造特征對(duì)系統(tǒng)進(jìn)行進(jìn)一步優(yōu)化升級(jí)，以提高三維可視化系統(tǒng)的科學(xué)性。