曲來超，許江濤

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

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

曲來超1，許江濤2

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

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

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

20世紀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]等。從建模方式上看，可以歸為基于塊狀、面、體以及混合模型建模。然而上述建模方法大多存在以下不足：①建模過程繁瑣，人工干預較多；②自適應能力差，模型之間相互轉(zhuǎn)化困難；③沒有結(jié)合礦山地質(zhì)體的特點，建立的模型難以應用到實際中；④沒有綜合考慮地層尖滅、斷層、煤層等情況，建立的模型缺乏真實感。

本文采用QTPV（quasi tri-prism volume）建模方法，并作了改進和完善，應用到GIS軟件系統(tǒng)中。傳統(tǒng)的GTP建模方法是基于鉆孔數(shù)據(jù)進行構(gòu)建模型，但對于復雜的含有斷層、尖滅、煤層的地質(zhì)體顯得略有不足，而且地質(zhì)勘探的鉆孔數(shù)量較為稀少，構(gòu)造精度比較低。本文利用地質(zhì)測量數(shù)據(jù)中鉆孔、地質(zhì)勘測剖面圖數(shù)據(jù)，結(jié)合二者之間的位置關(guān)系，形成更多的虛擬鉆孔，參與地質(zhì)建模，精度較高，更加合理，更加可靠。結(jié)合地質(zhì)學的地質(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ì)學的系統(tǒng)、組、段以及巖性對地層進行順序編號，遇到煤層作為其中的一個地層編號。

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

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

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

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

2.1 似三棱柱體元定義

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

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

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

圖1 三棱柱體

圖2 似三棱柱體特例

2.2 似三棱柱體建模算法

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

以鉆孔孔口坐標為數(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 含斷層的復雜地質(zhì)構(gòu)造的推理[1]

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

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

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

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

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ī)則，沿鉆孔向下擴展新的三角形（下三角形）。

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

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

6）重復步驟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ì)剖面圖上以該點為基準,向上向下擴展一條垂直鉆孔,使該條虛擬鉆孔與各地層界面相交，記錄交點坐標以及巖層屬性,存儲在鉆孔數(shù)據(jù)表格內(nèi)，如圖3、4所示。

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

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

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

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

3 研究應用

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

4 結(jié) 語

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

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

[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].武漢測繪科技大學學報,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] 韓國建,郭達志,金學林.礦體信息的八叉樹存儲和檢索技術(shù)[J].測繪學報,1992,21(1):13-17

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

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

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

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

P208

B

1672-4623（2015）01-0059-03

10.3969/j.issn.1672-4623.2015.01.020

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

2013-12-25。

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