淺析工程地質(zhì)三維建模和可視化

1、淺析工程地質(zhì)三維建模和可視化            摘要：現(xiàn)有的地理信息系統(tǒng)(GIS)都主要表達二維的地表地物的圖形和屬性信息，要擴展到真三維包含地下地質(zhì)結(jié)構(gòu)的地質(zhì)信息系統(tǒng)還有差距。本文分析了工程地質(zhì)體三維建模與可視化研究的關(guān)鍵技術(shù)問題,并結(jié)合金沙江某水電工程壩址區(qū)的三維地質(zhì)建模進行了初步開發(fā)和應用。 關(guān)鍵詞：工程地質(zhì) 三維建模與可視化 地質(zhì)信息  1 前言 現(xiàn)有的地理信息系統(tǒng)(GIS)都主要表達二維的地表地物的圖形和屬性信息，要擴展到真三維包含地下地質(zhì)結(jié)構(gòu)的地質(zhì)信

2、息系統(tǒng)還有差距。一個大型地質(zhì)工程項目從可行性研究階段、初步設(shè)計階段到詳細設(shè)計階段，乃至到工程施工與運行階段，往往積累了大量的地質(zhì)資料，用三維模型圖形圖像來表達和解釋如此龐大的資料，比光靠數(shù)據(jù)庫和圖表圖紙等傳統(tǒng)手段來得有效的多。建立工程地質(zhì)體的三維模型，處理巖層界面與結(jié)構(gòu)面組合關(guān)系，逼真反映地下主要地質(zhì)結(jié)構(gòu)全貌，將為工程地質(zhì)工作者分析研究工程地質(zhì)現(xiàn)象和發(fā)現(xiàn)掌握巖土體結(jié)構(gòu)規(guī)律，提供一種嶄新的研究手段和研究方法。國外三維地質(zhì)建模和可視化研究發(fā)展較快。加拿大阿波羅科技集團公司推出的三維建模與分析軟件MicroLYNX，通過對離散點采樣、鉆探采樣和探槽采樣等空間數(shù)據(jù)的處理，產(chǎn)生剖面、塊和面等模型，確定礦

3、藏分布和等級變化并計算礦藏儲量。加拿大Gemcom Software International Inc.公司開發(fā)的Gemcom軟件通過鉆孔、點、多邊形等數(shù)據(jù)，利用實用的圖形編輯和生成工具，顯示鉆孔孔位分布，運用不規(guī)則三角網(wǎng)建立表面和實體模型，運用多義線圈閉巖層和礦體邊界進行儲量和品位分析，提供了交互操作功能并允許用戶根據(jù)自己的經(jīng)驗和專家知識勾畫地質(zhì)模型，實現(xiàn)任意剖面切割任意角度觀察和實體與實體或?qū)嶓w與表面的交切與布爾運算等。國外軟件主要是瞄準采礦工程，能夠較好地滿足采礦工程活動中的礦產(chǎn)資源勘探和評價、地下礦井和露天礦坑設(shè)計和規(guī)劃、礦產(chǎn)資源管理和采礦生產(chǎn)管理等需求。美國Kinetix公司開發(fā)的

4、3D Studio MAX，Alias/Wavefront公司開發(fā)的Maya和微軟公司開發(fā)的Softimage等大眾化的三維建模軟件，在構(gòu)建工業(yè)和建筑模型與動畫制作方面有其獨到之處，但交互查詢的功能較弱，與工程勘測數(shù)據(jù)庫結(jié)合并應用于工程地質(zhì)三維建模方面還有較大距離。張菊明等對風化帶分布、多層地層等地質(zhì)信息的可視化和斷層錯斷巖層的表達和顯示的算法1,2進行了較為深入的研究，為工程地質(zhì)三維可視化軟件的開發(fā)準備了數(shù)學基礎(chǔ)，并借助AutoCAD平臺實現(xiàn)了復雜三維地質(zhì)圖形的顯示。國內(nèi)的靈圖VRMap地理信息系統(tǒng)軟件有較強的地形模擬和地表地物的查詢功能，但不是真三維的地質(zhì)建模工具。北京東方泰坦科技有限公司

5、開發(fā)TITAN三維建模軟件，基于框架建模的思想，利用平行或基本平行的剖面數(shù)據(jù)，建立起三維空間復雜形狀物體的真三維實體模型，但目前只是初步的三維建模與圖形處理的引擎，在面向具體專業(yè)時，需要添加或擴充專業(yè)模塊，比如工程地質(zhì)專業(yè)模塊等。縱觀國內(nèi)外幾種軟件的研究與開發(fā)現(xiàn)狀，它們?yōu)楣こ痰刭|(zhì)三維建模與可視化打下了很好的技術(shù)基礎(chǔ)，提供了很寶貴的開發(fā)經(jīng)驗。但是，對于工程地質(zhì)專業(yè)的地質(zhì)體建模與可視化分析的針對性不強，不能夠很好地滿足工程地質(zhì)生產(chǎn)與研究的專業(yè)功能需要。因此本文將從分析工程地質(zhì)的三維建模和可視化的關(guān)鍵技術(shù)問題入手，簡單描述作者在工程地質(zhì)三維建模和可視化方面的初步開發(fā)研究成果。2 關(guān)鍵技術(shù)問題分析 2

6、.1離散數(shù)據(jù)的插值與擬合 工程地質(zhì)復雜地質(zhì)體中的各種地質(zhì)信息，包括地表地形、地下水位、地層界面、斷層、節(jié)理、風化帶分布、侵入體及各種地球物理、地球化學、巖土體的物理力學參數(shù)或數(shù)據(jù)的等值面(線)等，都可以看作是三維空間中的函數(shù)，它們的擬合函數(shù)要根據(jù)實際勘測數(shù)據(jù)建立，實測數(shù)據(jù)越豐富，越能夠真實描繪出這些信息的空間分布規(guī)律。地表地形測量數(shù)據(jù)、地下水位埋深測量信息等的單值曲面圖形生成可歸結(jié)為雙自變量離散數(shù)據(jù)的插值和擬合，多值曲面如倒轉(zhuǎn)褶皺和空間等值面等，則應采用多參變量插值等其他一些較復雜的方法。空間曲面插值函數(shù)有以下構(gòu)造方法，如與距離成反比的加權(quán)方法(Shepard 方法)，徑向基函數(shù)插值法

7、0;   2007-04-22        (Multiquadric方法)3，平面彈性理論插值法1,2等，它們同樣適用于單個連續(xù)地層界面、地球物理勘探數(shù)據(jù)、地球化學勘探數(shù)據(jù)以及巖土體物理力學參數(shù)在地質(zhì)體空間的分布。2.2 三維數(shù)據(jù)結(jié)構(gòu)        工程地質(zhì)體一般是不規(guī)則形體，在計算機圖形學中曲線和曲面總是分別通過很多微小直線段和微小三角面逼近來模擬地層巖性界線和巖層曲面，即巖層界面（和地表曲線、地下水位面等地質(zhì)層面界線）

8、和巖層曲面都分別是許多微小直線段和微小三角面的集合。地質(zhì)體三維空間數(shù)據(jù)結(jié)構(gòu)是工程地質(zhì)三維建模和可視化的基礎(chǔ)，這就要求必須具備有效的分層的三維數(shù)據(jù)結(jié)構(gòu)，能夠確保人機交互和查詢的實現(xiàn)。2.3 曲面求交 地質(zhì)體中存在大量各種層面，當出現(xiàn)地層不整合、斷層錯斷巖層、地層尖滅和地下水出露于河谷地表等情形時，就自然會遇到曲面間求交的問題；地質(zhì)體三維模型的上部邊界是地表曲面，通過數(shù)學方法擬合出的巖層面或地下水位面不應超出地表曲面，即超出部分不應顯示。同樣的，當顯示多層地層時，下面的每一巖層應以其上一巖層為邊界。因此，為了可視化地層界面必須要解決地層面與地表、斷層面和其他地層面的求交問題。另一方面，在剖面圖成圖

9、時，地質(zhì)界線的繪制是通過顯示剖面(平面)與各種地質(zhì)界面(曲面)求交所得出的交線。因此曲面求交包括地質(zhì)界面（層面）之間的相交，和地質(zhì)界面與剖面的相交兩類問題。2.4 三維拓撲結(jié)構(gòu)分析 從地質(zhì)學角度看，拓撲是地質(zhì)對象間關(guān)系的表格，拓撲表存儲層位間上覆、下伏和交切(被斷層切割后地層的拓撲表達)等的地層學關(guān)系及地質(zhì)空間位置關(guān)系。拓撲也可視為允許這些地質(zhì)關(guān)系合理儲存的數(shù)據(jù)結(jié)構(gòu)。例如，考慮多層地層，上一個巖層的底面和與其相鄰的下一個巖層的頂面是上下巖層這兩個實體的公共部分或共享邊界，它們之間的拓撲關(guān)系就是相鄰和同一的關(guān)系，在存儲數(shù)據(jù)時只存儲上一個巖層的底面或其相鄰的下一個巖層的頂面，即相鄰巖層的邊界曲面可

10、以存為一個地層曲面，大大減少數(shù)據(jù)存儲量。評價地質(zhì)模型系統(tǒng)的優(yōu)缺點往往決定于描述地質(zhì)對象所用的拓撲結(jié)構(gòu)4。2.5 可視化技術(shù)        工程地質(zhì)復雜地質(zhì)體可視化，是利用計算機技術(shù)將工程勘測獲得的數(shù)據(jù)，轉(zhuǎn)換為形象直觀的便于進行交互分析的地下地質(zhì)結(jié)構(gòu)空間形態(tài)的立體圖和剖面圖形，其基礎(chǔ)是工程數(shù)據(jù)和測量數(shù)據(jù)的可視化5。利用可視化技術(shù)可以從龐大的地質(zhì)勘測數(shù)據(jù)中構(gòu)造出地質(zhì)工程中對于邊破穩(wěn)定性和地下硐室變形破壞等起關(guān)鍵作用的巖層和結(jié)構(gòu)面，并顯示其范圍、走向和相互交切關(guān)系，幫助工程地質(zhì)人員對原始數(shù)據(jù)做出正確解釋，繼而為工程地質(zhì)分析具體問題提供決策支

11、持。    3 工程地質(zhì)三維可視化技術(shù)的初步開發(fā)與應用 3.1 研究框圖 工程地質(zhì)復雜地質(zhì)體三維建模與可視化的研究框圖如圖    2007-04-22        1所示?；陔x散采樣數(shù)據(jù)的插值與擬合的思想，即將離散數(shù)據(jù)轉(zhuǎn)化為連續(xù)曲線曲面, 工程地質(zhì)復雜地質(zhì)體三維建模與可視化的過程是，從勘探數(shù)據(jù)庫中提取各種地質(zhì)信息的坐標位置及巖土體的物理力學參數(shù)，通過不同的擬合與插值函數(shù)得到地質(zhì)層面(曲面)和地質(zhì)實體的三維計算機圖形顯示，表達地質(zhì)信息

12、在研究區(qū)域內(nèi)的分布規(guī)律。生成地質(zhì)巖層面和地質(zhì)實體后，實現(xiàn)從任意角度觀察建立的模型，實現(xiàn)根據(jù)指定的剖面走向、傾向和傾角生成垂直剖面。3.2 初步開發(fā)與應用 在收集整理現(xiàn)場勘測數(shù)據(jù)后錄入金沙江某水電工程勘測空間數(shù)據(jù)庫各分項數(shù)據(jù)表，這些數(shù)據(jù)表不僅包括地質(zhì)信息的位置數(shù)據(jù)，更重要的是提供屬性數(shù)據(jù)。 以地層巖性數(shù)據(jù)表為例，要求錄入鉆孔編號、巖層起始深度、巖層終止深度、層厚、巖性（地層名稱）、地層代碼（地層年代）、巖層走向、巖層傾向、巖層傾角、接觸關(guān)系、地質(zhì)描述等數(shù)據(jù)。隨著工程勘測的進展，能夠方便地修改補充和管理勘測數(shù)據(jù)。圖2是工程勘測數(shù)據(jù)庫中鉆孔地層系統(tǒng)數(shù)據(jù)表的管理界面。 通過孔口坐標和測量數(shù)據(jù)等的離散數(shù)據(jù)的擬合和插值法繪制壩址區(qū)的右岸地表曲面網(wǎng)格（圖3），進而可在三維圖形環(huán)境中進行虛擬現(xiàn)實瀏覽觀察（圖4）。 利用工程勘測數(shù)據(jù)，建立了壩址區(qū)右岸三維立體地質(zhì)圖。該壩址區(qū)自上而下地層巖性組合為：第四系崩坡堆積物，侏羅系泥巖、粉砂質(zhì)泥巖、泥質(zhì)粉砂巖，三疊系上統(tǒng)厚至巨厚層狀細至中粒砂巖，三疊系上統(tǒng)薄至中厚層狀粉細紗巖、粉砂巖，三疊系上統(tǒng)中厚至厚層狀中粗砂巖。通過有