曲來超，許江濤

（1.信息工程大學(xué) 地理空間信息學(xué)院，河南 鄭州 451400；2.河南工業(yè)和信息化職業(yè)學(xué)院，河南 焦作 454000）

用于數(shù)字礦山的復(fù)雜地質(zhì)體三維建模方法

曲來超1，許江濤2

（1.信息工程大學(xué) 地理空間信息學(xué)院，河南 鄭州 451400；2.河南工業(yè)和信息化職業(yè)學(xué)院，河南 焦作 454000）

利用三維地質(zhì)體的優(yōu)點提出了一種建模新方法。該方法以似三棱柱作為建模的基本體元，根據(jù)似三棱柱的特點和知識推理規(guī)則，進(jìn)行復(fù)雜地質(zhì)構(gòu)造的自動建模,避免了不必要的人為干預(yù)，擴(kuò)展了似三棱柱建模的適用范圍和表現(xiàn)能力。最后通過研究實例展示該方法的實際建模效果。

數(shù)字礦山；似三棱柱；地質(zhì)體

20世紀(jì)70年代以來，三維地質(zhì)建模技術(shù)不斷發(fā)展，如CSG（constructive solid geometry）模型[1]、B-Reps（boundary representation） 模 型[2]、TIN（triangulation irregular network）與GRID的集成[3]、TEN（tetrahedron network）[4,5]與八叉樹的混合模型[6]等。從建模方式上看，可以歸為基于塊狀、面、體以及混合模型建模。然而上述建模方法大多存在以下不足：①建模過程繁瑣，人工干預(yù)較多；②自適應(yīng)能力差，模型之間相互轉(zhuǎn)化困難；③沒有結(jié)合礦山地質(zhì)體的特點，建立的模型難以應(yīng)用到實際中；④沒有綜合考慮地層尖滅、斷層、煤層等情況，建立的模型缺乏真實感。

本文采用QTPV（quasi tri-prism volume）建模方法，并作了改進(jìn)和完善，應(yīng)用到GIS軟件系統(tǒng)中。傳統(tǒng)的GTP建模方法是基于鉆孔數(shù)據(jù)進(jìn)行構(gòu)建模型，但對于復(fù)雜的含有斷層、尖滅、煤層的地質(zhì)體顯得略有不足，而且地質(zhì)勘探的鉆孔數(shù)量較為稀少，構(gòu)造精度比較低。本文利用地質(zhì)測量數(shù)據(jù)中鉆孔、地質(zhì)勘測剖面圖數(shù)據(jù)，結(jié)合二者之間的位置關(guān)系，形成更多的虛擬鉆孔，參與地質(zhì)建模，精度較高，更加合理，更加可靠。結(jié)合地質(zhì)學(xué)的地質(zhì)構(gòu)造知識，通過對底層的合理編號，構(gòu)造出含煤層，能夠了解到煤層上下2部分的巖層信息。

1 數(shù)據(jù)來源與數(shù)據(jù)處理

1.1 數(shù)據(jù)來源

地質(zhì)信息數(shù)據(jù)多樣，來源廣泛，本文主要采用鉆孔柱狀圖和地質(zhì)勘測剖面圖數(shù)據(jù)（表1）。

表1 地質(zhì)信息來源與分類

1.2 數(shù)據(jù)處理

首先，根據(jù)地質(zhì)勘探剖面圖和地質(zhì)學(xué)的系統(tǒng)、組、段以及巖性對地層進(jìn)行順序編號，遇到煤層作為其中的一個地層編號。

其次，插入虛擬鉆孔，在地質(zhì)勘測剖面圖上,對于斷層、夾層的臨界點進(jìn)行插入虛擬鉆孔處理。

最后，把鉆孔數(shù)據(jù)和虛擬鉆孔數(shù)據(jù)存入數(shù)據(jù)表格中，作為基礎(chǔ)數(shù)據(jù)存放起來（表2）。

表2 地質(zhì)信息基礎(chǔ)數(shù)據(jù)

2 基于似三棱柱三維數(shù)據(jù)模型及其相關(guān)算法

2.1 似三棱柱體元定義

標(biāo)準(zhǔn)的三棱柱（orthogonal tri-prism，OTP）[7,8]是由一個三角形沿某一方向拉伸一定距離之后，所形成的體元。它包含6個頂點、相互平行的3條棱邊和3個側(cè)面四邊形，如圖1a所示。

似三棱柱體是體元的上下三角形不互相平行、側(cè)棱長度不等且不相互平行、側(cè)面四邊形的4個頂點不在一個平面內(nèi)的三棱柱體，如圖1b。

實際應(yīng)用中，上述的似三棱柱體仍難以描述復(fù)雜的地質(zhì)體，通常是似三棱柱體特例的組合。所謂特例是似三棱柱體的3個頂點部分重合引起的，根據(jù)不同頂點重合的組合，概括起來可以有4種特例，如圖2?？梢钥闯鯽是四棱椎體，b、d是四面體，它們都是似三棱柱的特例。事實上，一個標(biāo)準(zhǔn)的似三棱柱可以拆分為2個四面體，這也為體元模型之間相互轉(zhuǎn)換提供了理論依據(jù)。

圖1 三棱柱體

圖2 似三棱柱體特例

2.2 似三棱柱體建模算法

2.2.1 地表面TIN的構(gòu)建

以鉆孔孔口坐標(biāo)為數(shù)據(jù)點，按照如下步驟構(gòu)建Delaunay三角網(wǎng)。

1）如果沒有露出地面的斷層數(shù)據(jù)點，采用逐點插入法建模，基本步驟：①定義一個包含所有鉆孔數(shù)據(jù)點的初始多邊形。②在初始多邊形中建立初始三角網(wǎng)，然后迭代一下步驟，直到所有的鉆孔點都被處理：插入數(shù)據(jù)點P，在三角網(wǎng)中找出包含P的三角形T，把P與T的3個頂點相連，生成3個新的三角形。③用局部最優(yōu)方法(LOP)優(yōu)化三角網(wǎng)，使其符合Delaunay法則。

2）如果有露出地表面的斷層，以斷層線為約束，采用約束Delaunay三角剖分算法[9]。

2.2.2 含斷層的復(fù)雜地質(zhì)構(gòu)造的推理[1]

對上述Delaunay構(gòu)建的三角形，遍歷每個三角形，對其中的一個三角形作擴(kuò)展知識邏輯推理處理：

1）如果當(dāng)前帶擴(kuò)展的三角形3個頂點均無斷層點，則檢查3個點的編號。若3個頂點的編號相同，那么，沿鉆孔向下擴(kuò)展新的三角形，其3個頂點為當(dāng)前3個頂點所在相應(yīng)鉆孔的下一個點。若3個頂點的編號不同，那么，編號小的頂點沿相應(yīng)鉆孔擴(kuò)展下一個點為新的三角形頂點，而編號大的頂點保持不變。

2）如果當(dāng)前待擴(kuò)展的三角形3個頂點中有1個或2個斷層點 ，則檢查其余頂點所在相應(yīng)鉆孔中自該點向下的屬性，若含有與上述1個或2個斷層點屬于同一斷層面的點，則該斷層點無論地層編號大小均不向下擴(kuò)展，直到其余頂點按編號大小向下擴(kuò)展到與該斷層點的屬性相同為止。

3）如果當(dāng)前帶擴(kuò)展三角形的3個頂點全為斷層點，按照規(guī)則2）檢查各自地層編號，然后根據(jù)編號大小向下擴(kuò)展新三角形。

2.2.3 含斷層的地質(zhì)體三維建模方法[10]

1）根據(jù)§2.2.1構(gòu)建表面Delaunay三角形。

2）從地表面Delaunay三角網(wǎng)中提取1個三角形，將這個三角形設(shè)置為第一個廣義三棱柱的上三角形。

3）根據(jù)上述三角形3個頂點的編號，按照§2.2.2中三維構(gòu)造知識推理規(guī)則，沿鉆孔向下擴(kuò)展新的三角形（下三角形）。

4）根據(jù)上下三角形對應(yīng)關(guān)系及鉆孔列表構(gòu)建廣義三棱柱，記錄其描述信息，并將下三角形置為上三角形。

5）重復(fù)步驟3）、4），直到上三角形的頂點均為鉆孔底部點為止。

6）重復(fù)步驟1）～5），直到地表表面的TIN的三角形遍歷。

2.3 似三棱柱體建模特殊問題處理

2.3.1 地層相交點處、斷層點、煤層點建模處理

一般來說，野外鉆探的鉆孔切穿點不會與地層界面的相交處恰好相切，如果這樣建模的話，就會漏掉地層中尖滅現(xiàn)象，因而，要對此做特殊處理，根據(jù)勘探剖面圖，找出鉆孔與地層的相交點。

斷層是改變了地層走向形成的地質(zhì)錯斷走向，需要在勘探剖面圖上，找出地層與斷層線的交點。

煤層是礦山地質(zhì)體關(guān)注的重點，需要單獨建立。由于煤層有上、下表面，就要在勘探剖面圖找出煤層上、下表面尖滅點。

由于上述3種點改變了地質(zhì)體的構(gòu)造，我們把這幾種點廣義地稱為斷點。

2.3.2 建模步驟

1）插入虛擬鉆孔。上述的3種情況相交點確定后,在地質(zhì)剖面圖上以該點為基準(zhǔn),向上向下擴(kuò)展一條垂直鉆孔,使該條虛擬鉆孔與各地層界面相交，記錄交點坐標(biāo)以及巖層屬性,存儲在鉆孔數(shù)據(jù)表格內(nèi)，如圖3、4所示。

2）鉆孔三角化。將上述虛擬鉆孔并入真實鉆孔中，然后根據(jù)虛擬鉆孔與真實鉆孔中心進(jìn)行三角化。

3）采用向下擴(kuò)展三角形方法建立似三棱柱模型。

圖3 地層相交處、煤層表面點插入虛擬鉆孔示意圖

圖4 斷層點插入虛擬鉆孔示意圖

3 研究應(yīng)用

根據(jù)本文的建模方法，以SQL Server 2010為數(shù)據(jù)庫支撐，依據(jù)Direct 9.0c和C#語言，采用河北某礦區(qū)為實驗數(shù)據(jù)，建立了應(yīng)用于礦山的地質(zhì)體模型。該礦區(qū)大約10 km2，鉆孔數(shù)據(jù)位539個，生成單個似三棱柱體元11 035個，耗時53 s。地層共分15層構(gòu)建，有2個煤層，分別是2號煤和9號煤，2號煤層平均厚度5 m，9號煤平均厚度10 m，都是適合開采的煤層。研究表明，該建模方法效率高、精度好，能真實有效地還原地下三維場景，能應(yīng)用到實際的三維軟件系統(tǒng)中，并能支撐如剖切、防治水壓力測試等各種功能開發(fā)。

4 結(jié) 語

結(jié)合礦山真實數(shù)據(jù)，通過QTPV建模方式構(gòu)建了三維地質(zhì)模型，并根據(jù)特殊地質(zhì)構(gòu)造進(jìn)行了處理,對QTPV建模方法做了進(jìn)一步的改進(jìn)和完善。改進(jìn)后建模方法具有以下優(yōu)點：節(jié)省了存儲空間；能夠精確地描述實體內(nèi)部結(jié)構(gòu)；建模結(jié)果唯一；能夠?qū)崿F(xiàn)面模型與體模型的相互轉(zhuǎn)換；可以實現(xiàn)地表與地下統(tǒng)一建模；方便模型的可視化；方便建立空間拓?fù)潢P(guān)系，有利于空間分析。

礦山地質(zhì)現(xiàn)象的高度復(fù)雜性和不確定性，以及局部區(qū)域的嚴(yán)重變形、扭曲，需要專業(yè)知識和工程經(jīng)驗融入到實際建模過程中。建立模型內(nèi)部點、線、面相應(yīng)的拓?fù)潢P(guān)系，從而對地質(zhì)體進(jìn)行剖面分析，儲量分析，這是本論文進(jìn)一步研究的方向。

[1] Li R X.Date Structure and Application Issues in 3D Geographic Information System[J].Geometrics,1994,48(3):209-244

[2] 龔健雅,夏宗國.矢量與柵格集成的三維數(shù)據(jù)模型[J].武漢測繪科技大學(xué)學(xué)報,1997,22(1):7-15

[3] Holding S W.3D Geoscience Modeling:Computer Techniques for Geological Characterization[M].Springer-Verlag,1994

[4] Pilout M,Tempfli K,Molenaar M.A Tetrahedron-Based on 3D Vector Data Model for Geo information Advanced Geographic Data Modeling[J].Delft Netherlands,1994(40):129-140

[5] Chen X Y.A Workstation for Three-Dimensional Spatial Data Research[C].The Fourth International Symposium of LIESMARS,Wuhan,1995

[6] 韓國建,郭達(dá)志,金學(xué)林.礦體信息的八叉樹存儲和檢索技術(shù)[J].測繪學(xué)報,1992,21(1):13-17

[7] 武強(qiáng),徐華.虛擬礦山系統(tǒng)中三維斷層模擬技術(shù)[J].遼寧工程技術(shù)大學(xué)學(xué)報,2005,24(3):316-319

[8] 李立新,譚建榮.約束Delaunay三角剖分中強(qiáng)行嵌入約束邊的多對角線交換算法[J].計算機(jī)學(xué)報,1999,22(10):1 115-1 118

[9] 陳學(xué)習(xí),吳立新,車德福，等.基于鉆孔數(shù)據(jù)的含斷層地質(zhì)體三維建模方法[J].煤田地質(zhì)與勘探,2005,33(5):5-8

[10] 程朋根,龔健雅,史文中,等.基于似三棱柱體的地質(zhì)體三維建模與應(yīng)用研究[J].武漢大學(xué)學(xué)報:信息科學(xué)版,2007,29(7):603-607

P208

B

1672-4623（2015）01-0059-03

10.3969/j.issn.1672-4623.2015.01.020

曲來超，碩士，研究方向為地理信息系統(tǒng)的應(yīng)用。

2013-12-25。

項目來源：國家高技術(shù)研究發(fā)展計劃資助項目（2001AA131021）；教育部科學(xué)技術(shù)研究重點資助項目（03008）；三維空間信息獲取與地學(xué)應(yīng)用教育部重點實驗室開放基金資助項目。