三維地質(zhì)模型構建研究

李輝，韓曉暉，王新廣

(1.天津市勘察院，天津　300191；　2.星際空間(天津)科技發(fā)展有限公司，天津　300222)

?

三維地質(zhì)模型構建研究

李輝1*，韓曉暉2，王新廣1

(1.天津市勘察院，天津300191；2.星際空間(天津)科技發(fā)展有限公司，天津300222)

摘要：為改變傳統(tǒng)的城市地質(zhì)工作和城市地下空間利用的粗放形式，彌補傳統(tǒng)的地質(zhì)學研究模式、工作方法與技術手段的不足，在總結(jié)分析傳統(tǒng)地質(zhì)研究方法與研究手段的基礎上，通過提取、完善數(shù)據(jù)庫中的鉆孔信息和地層信息，利用三維平臺接口將地層信息進行三維可視化建設，構建三維地質(zhì)模型，從而為地質(zhì)工作提供科學的依據(jù)和技術支撐。

關鍵詞：三維；地質(zhì)；模型；構建方法

1引言

地質(zhì)問題已成為貫穿于數(shù)字城市、減災防災、城建工程、地下工程、水電工程、交通工程、環(huán)境工程、資源開發(fā)的一個基礎性問題。傳統(tǒng)的二維分析方法廣泛應用在資源勘探與開發(fā)、工程地質(zhì)、水文地質(zhì)和礦床地質(zhì)等地質(zhì)相關領域中，通常采用剖面圖、等值線圖等二維圖件來表達地質(zhì)體的構造形態(tài)特征和屬性特征。這些二維圖件直觀性較差，地質(zhì)人員需要通過想象才能重構出地質(zhì)體的空間形態(tài)。

如果用三維地質(zhì)模型來直觀地表達地質(zhì)體，將大大提高傳統(tǒng)地質(zhì)學與地質(zhì)工作的研究水平。給地質(zhì)人員提供可視化的、詳細的地質(zhì)體資料，使其能直接、完整、準確地理解和分析地質(zhì)體的特征，這將給地質(zhì)工作帶來極大的便利。

三維地質(zhì)模型模擬的關鍵技術是空間對象三維模型的建立，充分利用現(xiàn)有的二維GIS數(shù)據(jù)(主要是鉆孔信息、地層信息、實驗數(shù)據(jù)等)，進行三維地質(zhì)模型建模，是三維地質(zhì)模型可視化技術難點及關鍵所在。因為這些信息能夠提供地質(zhì)(地層)構造的分布信息，能完整、準確地表達復雜地質(zhì)現(xiàn)象的邊界條件及地質(zhì)體內(nèi)包含的各種地質(zhì)構造，所以利用離散的鉆孔進行地質(zhì)曲面可視化重構是地質(zhì)構造最為直觀的研究手段。

2三維地質(zhì)模型構建算法分析

地質(zhì)模型構建方法利用已知散亂的鉆孔數(shù)據(jù)對曲面上的點進行擬合插值，插值后建立數(shù)學模型，使用三角構網(wǎng)方法將鉆孔數(shù)據(jù)連接成為曲面幾何信息，利用傳統(tǒng)的計算機圖形學方法加以顯示，從而獲取對數(shù)據(jù)場內(nèi)部信息的觀察。常用的插值方法為克里金插值方法，采用的三角構網(wǎng)方法為Delannay三角剖分法。

2.1克里金插值法

克里金插值法也被稱為空間局部估計或空間局部插值法，是地質(zhì)統(tǒng)計學的主要內(nèi)容之一。克里金插值法是在空間相關范圍分析的基礎上，用相關范圍內(nèi)的采樣點來估計待插點屬性值，在有限區(qū)域內(nèi)對區(qū)域化變量的取值進行無偏最優(yōu)估計。

克里金插值法具有較好的空間形態(tài)恢復能力，誤差較小、適用性較高，插值效果較好，逐漸被人們廣泛地應用到三維地質(zhì)模型的構建中。但它是局部插值方法，適用于精度要求較高的場合。本文將以克里金算法為基礎進行研究。

2.2Delaunay三角剖分法

三角剖分對數(shù)值分析和圖形學都是極為重要的一項預處理技術。尤其是Delaunay三角網(wǎng)是“最接近于規(guī)則化”的三角網(wǎng)，具有唯一性、空圓性、最大化最小角等特性。基本思路是利用散亂數(shù)據(jù)構造插值曲面，在此基礎上可以利用線性或有限元方法進行插值來繪制曲面或等值線，也可直接采用三角網(wǎng)進行等值線或曲面繪制。Delaunay三角剖分法算法簡單易實現(xiàn)，結(jié)構良好，數(shù)據(jù)結(jié)構簡單，數(shù)據(jù)冗余度小，存儲效率高，可適應各種分布密度的數(shù)據(jù)。

3三維地質(zhì)模型構建方法的研究與實現(xiàn)

作者在考查了大量三維地質(zhì)模型算法和工具的基礎上，從復雜度、適用性、精確性等方面綜合考慮，提出了三維地質(zhì)模型構建的初步解決方案。本方案通過提取、完善數(shù)據(jù)庫中的鉆孔信息和地層信息，通過三維平臺接口將地層信息進行三維可視化建設，生成地質(zhì)地層模型。具體構建過程如圖1所示：

3.1構建區(qū)域范圍

根據(jù)鼠標輸入或者對外接口輸入的坐標值串為邊界，確定建模區(qū)域范圍。

3.2提取地層數(shù)據(jù)

利用空間疊加算法選擇待建模區(qū)域內(nèi)的鉆孔數(shù)據(jù)、地層數(shù)據(jù)及其實驗數(shù)據(jù)，將用于構建三維地質(zhì)模型必須的標高、深度、平面坐標、巖性等信息進行提取，完成基礎信息搜集與記錄。

3.3計算地層信息

提取數(shù)據(jù)庫中鉆孔數(shù)據(jù)對應的地層分層信息，利用對應鉆孔的標高記錄，計算地層的上限和下限，完成項目鉆孔對應的地層信息搜集。

3.4插值計算

基于每個地層的信息，根據(jù)事先設定的插入密度，運用克里金算法計算該層所有特征點處以及地層模型邊界點處的插入值。

3.5構建網(wǎng)格

以上述的地層數(shù)據(jù)為基礎，按照地層結(jié)構構建網(wǎng)格插入點。以網(wǎng)格插入點為特征點，生成適宜的空間多邊體數(shù)據(jù)結(jié)構。

3.6地層調(diào)整

在三維地層組織結(jié)構中經(jīng)常會出現(xiàn)地層缺失和地層透鏡體情況，應根據(jù)地層實際情況調(diào)整與該地層相關聯(lián)的其他地層結(jié)構。具體步驟是：

(1)當出現(xiàn)有缺失地層的鉆孔數(shù)據(jù)，其周圍最近距離的其他鉆孔數(shù)據(jù)關聯(lián)的地層結(jié)構會在兩者之間尖滅，同時調(diào)整該缺失地層緊鄰的上層和下層地層的結(jié)構。

(2)當出現(xiàn)有透鏡體地層的鉆孔數(shù)據(jù)，其關聯(lián)的地層結(jié)構會在其周圍最近距離的其他鉆孔數(shù)據(jù)之間尖滅，同時調(diào)整該透鏡體地層緊鄰的上層和下層地層的結(jié)構。

地層尖滅算法涉及地層比較復雜，本文采用的算法思路是：首先判斷尖滅地層所關聯(lián)的地層的類型，將相同地層進行合并；再根據(jù)尖滅地層與其關聯(lián)地層的垂直方向的關系，確定尖滅地層的尖滅方向；第三步構建尖滅地層結(jié)構；最后構建關聯(lián)地層結(jié)構，完成地層構建。3.7構建三維地質(zhì)模型

將空間多邊體數(shù)據(jù)結(jié)構進行渲染可視化，形成可展示的地層三維模型，并進行模型優(yōu)化和紋理設置。具體步驟是：

(1)將地層標示圖例統(tǒng)一化處理；

(2)依據(jù)模型各個點的坐標以及坐標所在面的軸向，計算坐標的法線坐標。面的法線可以根據(jù)該面上任意三個點的坐標值求出，主要代碼如下：

//計算向量

vector1.x=pt1.x - pt2.x；

vector1.y=pt1.y - pt2.y；

vector1.z=pt1.z - pt2.z；

vector2.x=pt3.x - pt2.x；

vector2.y=pt3.y - pt2.y；

vector2.z=pt3.z - pt2.z；

//計算法線

normal.x=vector1.y*vector2.z-vector1.z*vector2.y；

normal.y=vector1.z*vector2.x-vector1.x*vector2.z；

normal.z=vector1.x*vector2.y-vector1.y*vector2.x；

//法線單位化

len=sqrt(normal._v[0] * normal._v[0]+ normal._v[1] * normal._v[1]+ normal._v[2] * normal._v[2] )；

if(len == 0.0f)len=1.0f；

normal.x /= len；

normal.y /= len；

normal.z /= len；

(3)計算地層紋理坐標。即根據(jù)法線坐標計算面的軸向，再依據(jù)面頂點坐標計算紋理坐標值。主要代碼如下：

for (unsigned i=0；igetNums()；++i)

{

N.set(vx[i].x()，vx[i].y()，vx[i].z())；//法線坐標

int axis=0；//軸向

double dnx=abs(N.x())；

double dny=abs(N.y())；

double dnz=abs(N.z())；

if(dny>dnx && dny>dnz)

axis=1；

if(dnz>dnx && dnz>dny)

axis=2；

//計算紋理坐標

switch (axis)

{

case 0： uv.set(P.y()，P.z())；break；

case 1： uv.set(P.x()，P.z())；break；

case 2： uv.set(P.x()，P.y())；break；

}

}

(4)逐層設置地層紋理。

4成果展示與分析

通過上述步驟，基于自主開發(fā)的三維系統(tǒng)，可以生成三維地質(zhì)模型?？梢钥吹剑鱾€地層都比較清晰，層與層之間的銜接也很準確，如圖2所示。

通過網(wǎng)格模型可以看到，地層起伏的擬合比較平滑，沒有漏洞和棱角，擬合效果較好，如圖3所示。

從成果來看，該三維地質(zhì)模型構建方法能較準確地反映地質(zhì)地層的情況，達到了預期的目標，可以作為進一步研究的基礎和依據(jù)。

5結(jié)語

本文初步實現(xiàn)了運用克里金算法和Delaunay三角構網(wǎng)方法，基于鉆孔信息的三維地質(zhì)模型的構建。此方法的優(yōu)點是實現(xiàn)簡單、速度較快，精度高，對地質(zhì)地層擬合較準確。不足之處是對于斷層、褶皺等情況缺少支持。針對研究現(xiàn)狀，下一步作者將解決不足，進一步優(yōu)化此算法，實現(xiàn)三維地質(zhì)模型構建方法在三維平臺上的典型應用。

參考文獻

[1]劉娜，謝英情，楚亮等. 基于鉆孔數(shù)據(jù)的三維地質(zhì)空間插值方法對比研究[J]. 地震研究，2008(S2)：69～74.

[2]靳國棟，劉衍聰，牛文杰. 距離加權反比插值法和克里金插值法的比較[J]. 長春工業(yè)大學學報，2003(3)：54～58.[3]羅智勇. 面向地質(zhì)勘查的三維可視化系統(tǒng)研制與開發(fā)[D]. 成都：成都理工大學，2008.

[4]楊耀紅. 三維地質(zhì)模型可視化分析技術研究與應用[D]. 武漢：中國地質(zhì)大學(武漢)，2007.

[5]郝海森，吳立新. 基于強約束Delaunay-TIN的三維地學模擬與可視化[J]. 地理與地理信息科學，2003(2)：15～18.[6]朱小弟. 三維地質(zhì)模型可視化技術研究[D]. 北京：中國礦業(yè)大學(北京)，1999.

[7]徐能雄，段慶偉等. 三維地質(zhì)建模方法及程序?qū)崿F(xiàn)[M]. 北京：地質(zhì)出版社，2011：23，67.

[8]程朋根，文紅. 三維空間數(shù)據(jù)建模及算法[M]. 北京：國防工業(yè)出版社，2011：141～143.

Build 3D Geological Model

Li Hui1，Han Xiaohui2，Wang Xinguang1

(1.Tianjin Institute of Geotechnical Investigation Surveying，Tianjin 300191，China；2.StarGIS (Tianjin) Technology Delelopment Co.，Ltd. ，Tianjin 300222，China)

Key words：3d；geological；model；method of building

Abstract：In order to change the traditional urban geological work and the extensive form of urban underground space utilization，make up for the traditional geology research model，working method and the shortage of technical means，in the summary analysis on the basis of traditional geological research methods and research means，by extracting，improve the drilling information and stratigraphic information in the database，using three-dimensional platform interface will stratigraphic information construction of 3 d visualization，build 3 d geological model，so as to provide scientific basis for geological work and technical support.

文章編號：1672-8262(2016)03-158-03

中圖分類號：P628.4

文獻標識碼：B

*收稿日期：2016—03—16

作者簡介：李輝(1982—)，女，工程師，現(xiàn)主要從事地理信息系統(tǒng)開發(fā)與研究。

基金項目：國家自然科學基金(21373210)