三維地質模型構建研究

李輝，韓曉暉，王新廣

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

?

三維地質模型構建研究

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

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

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

關鍵詞：三維；地質；模型；構建方法

1引言

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

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

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

2三維地質模型構建算法分析

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

2.1克里金插值法

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

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

2.2Delaunay三角剖分法

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

3三維地質模型構建方法的研究與實現

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

3.1構建區域范圍

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

3.2提取地層數據

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

3.3計算地層信息

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

3.4插值計算

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

3.5構建網格

以上述的地層數據為基礎，按照地層結構構建網格插入點。以網格插入點為特征點，生成適宜的空間多邊體數據結構。

3.6地層調整

在三維地層組織結構中經常會出現地層缺失和地層透鏡體情況，應根據地層實際情況調整與該地層相關聯的其他地層結構。具體步驟是：

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

(2)當出現有透鏡體地層的鉆孔數據，其關聯的地層結構會在其周圍最近距離的其他鉆孔數據之間尖滅，同時調整該透鏡體地層緊鄰的上層和下層地層的結構。

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

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

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

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

//計算向量

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)計算地層紋理坐標。即根據法線坐標計算面的軸向，再依據面頂點坐標計算紋理坐標值。主要代碼如下：

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成果展示與分析

通過上述步驟，基于自主開發的三維系統，可以生成三維地質模型。可以看到，各個地層都比較清晰，層與層之間的銜接也很準確，如圖2所示。

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

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

5結語

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

參考文獻

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

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

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

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

[7]徐能雄，段慶偉等. 三維地質建模方法及程序實現[M]. 北京：地質出版社，2011：23，67.

[8]程朋根，文紅. 三維空間數據建模及算法[M]. 北京：國防工業出版社，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—)，女，工程師，現主要從事地理信息系統開發與研究。

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