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

  現有的地理信息系統(GIS)都主要表達二維的地表地物的圖形和屬性信息, 要擴展到真三維包含地下地質(zhì)結構的地質(zhì)信息系統還有差距。一個(gè)大型地質(zhì)工程項目從可行性研究階段、初步設計階段到詳細設計階段, 乃至到工程施工與運行階段, 往往積累了大量的地質(zhì)資料, 用三維模型圖形圖像來(lái)表達和解釋如此龐大的資料, 比光靠數據庫和圖表圖紙等傳統手段來(lái)得有效的多。建立工程地質(zhì)體的三維模型, 處理巖層界面與結構面組合關(guān)系, 逼真反映地下主要地質(zhì)結構全貌, 將為工程地質(zhì)工作者分析研究工程地質(zhì)現象和發(fā)現掌握巖土體結構規律, 提供一種嶄新的研究手段和研究方法。

 

  國外三維地質(zhì)建模和可視化研究發(fā)展較快。加拿大阿波羅科技集團公司推出的三維建模與分析軟件Micro LYNX, 通過(guò)對離散點(diǎn)采樣、鉆探采樣和探槽采樣等空間數據的處理, 產(chǎn)生剖面、塊和面等模型, 確定礦藏分布和等級變化并計算礦藏儲量。加拿大Gemcom Software I nterna tio nal Inc .公司開(kāi)發(fā)的Gemcom 軟件通過(guò)鉆孔、點(diǎn)、多邊形等數據, 利用實(shí)用的圖形編輯和生成工具,顯示鉆孔孔位分布, 運用不規則三角網(wǎng)建立表面和實(shí)體模型, 運用多義線(xiàn)圈閉巖層和礦體邊界進(jìn)行儲量和品位分析, 提供了交互操作功能并允許用戶(hù)根據自己的經(jīng)驗和專(zhuān)家知識勾畫(huà)地質(zhì)模型, 實(shí)現任意剖面切割任意角度觀(guān)察和實(shí)體與實(shí)體或實(shí)體與表面的交切與布爾運算等。國外軟件主要是瞄準采礦工程, 能夠較好地滿(mǎn)足采礦工程活動(dòng)中的礦產(chǎn)資源勘探和評價(jià)、地下礦井和露天礦坑設計和規劃、礦產(chǎn)資源管理和采礦生產(chǎn)管理等需求。美國Kinetix 公司開(kāi)發(fā)的3D Studio MAX, Alias/Wav efront 公司開(kāi)發(fā)的May a 和微軟公司開(kāi)發(fā)的Softimage 等大眾化的三維建模軟件, 在構建工業(yè)和建筑模型與動(dòng)畫(huà)制作方面有其獨到之處, 但交互查詢(xún)的功能較弱, 與工程勘測數據庫結合并應用于工程地質(zhì)三維建模方面還有較大距離。

 

  張菊明等對風(fēng)化帶分布、多層地層等地質(zhì)信息的可視化和斷層錯斷巖層的表達和顯示的算法[ 1 , 2] 進(jìn)行了較為深入的研究,為工程地質(zhì)三維可視化軟件的開(kāi)發(fā)準備了數學(xué)基礎, 并借助AutoCAD平臺實(shí)現了復雜三維地質(zhì)圖形的顯示。國內的靈圖VRMap 地理信息系統軟件有較強的地形模擬和地表地物的查詢(xún)功能, 但不是真三維的地質(zhì)建模工具。北京東方泰坦科技有限公司開(kāi)發(fā)TITAN 三維建模軟件, 基于框架建模的思想, 利用平行或基本平行的剖面數據, 建立起三維空間復雜形狀物體的真三維實(shí)體模型, 但目前只是初步的三維建模與圖形處理的引擎, 在面向具體專(zhuān)業(yè)時(shí), 需要添加或擴充專(zhuān)業(yè)模塊, 比如工程地質(zhì)專(zhuān)業(yè)模塊等。

 

  縱觀(guān)國內外幾種軟件的研究與開(kāi)發(fā)現狀, 它們?yōu)楣こ痰刭|(zhì)三維建模與可視化打下了很好的技術(shù)基礎, 提供了很寶貴的開(kāi)發(fā)經(jīng)驗。但是, 對于工程地質(zhì)專(zhuān)業(yè)的地質(zhì)體建模與可視化分析的針對性不強, 不能夠很好地滿(mǎn)足工程地質(zhì)生產(chǎn)與研究的專(zhuān)業(yè)功能需要。因此本文將從分析工程地質(zhì)的三維建模和可視化的關(guān)鍵技術(shù)問(wèn)題入手, 簡(jiǎn)單描述作者在工程地質(zhì)三維建模和可視化方面的初步開(kāi)發(fā)研究成果。

 

  1  關(guān)鍵技術(shù)問(wèn)題分析

 

  1 .1  離散數據的插值與擬合

 

  工程地質(zhì)復雜地質(zhì)體中的各種地質(zhì)信息, 包括地表地形、地下水位、地層界面、斷層、節理、風(fēng)化帶分布、侵入體及各種地球物理、地球化學(xué)、巖土體的物理力學(xué)參數或數據的等值面(線(xiàn))等, 都可以看作是三維空間中的函數, 它們的擬合函數要根據實(shí)際勘測數據建立, 實(shí)測數據越豐富, 越能夠真實(shí)描繪出這些信息的空間分布規律。地表地形測量數據、地下水位埋深測量信息等的單值曲面圖形生成可歸結為雙自變量離散數據的插值和擬合, 多值曲面如倒轉褶皺和空間等值面等, 則應采用多參變量插值等其他一些較復雜的方法?？臻g曲面插值函數有以下構造方法, 如與距離成反比的加權方法(Shepard 方法), 徑向基函數插值法(Multiquadric 方法)[ 3] , 平面彈性理論插值法[ 1 , 2] 等, 它們同樣適用于單個(gè)連續地層界面、地球物理勘探數據、地球化學(xué)勘探數據以及巖土體物理力學(xué)參數在地質(zhì)體空間的分布。

 

  1 .2  三維數據結構

 

  工程地質(zhì)體一般是不規則形體, 在計算機圖形學(xué)中曲線(xiàn)和曲面總是分別通過(guò)很多微小直線(xiàn)段和微小三角面逼近來(lái)模擬地層巖性界線(xiàn)和巖層曲面, 即巖層界面(和地表曲線(xiàn)、地下水位面等地質(zhì)層面界線(xiàn))和巖層曲面都分別是許多微小直線(xiàn)段和微小三角面的集合。地質(zhì)體三維空間數據結構是工程地質(zhì)三維建模和可視化的基礎, 這就要求必須具備有效的分層的三維數據結構,能夠確保人機交互和查詢(xún)的實(shí)現。

 

  1 .3  曲面求交

 

  地質(zhì)體中存在大量各種層面, 當出現地層不整合、斷層錯斷巖層、地層尖滅和地下水出露于河谷地表等情形時(shí), 就自然會(huì )遇到曲面間求交的問(wèn)題;地質(zhì)體三維模型的上部邊界是地表曲面,通過(guò)數學(xué)方法擬合出的巖層面或地下水位面不應超出地表曲面, 即超出部分不應顯示。同樣的, 當顯示多層地層時(shí), 下面的每一巖層應以其上一巖層為邊界。因此, 為了可視化地層界面必須要解決地層面與地表、斷層面和其他地層面的求交問(wèn)題。另一方面, 在剖面圖成圖時(shí), 地質(zhì)界線(xiàn)的繪制是通過(guò)顯示剖面(平面)與各種地質(zhì)界面(曲面)求交所得出的交線(xiàn)。因此曲面求交包括地質(zhì)界面(層面)之間的相交, 和地質(zhì)界面與剖面的相交兩類(lèi)問(wèn)題。

 

  1 .4  三維拓撲結構分析

 

  從地質(zhì)學(xué)角度看, 拓撲是地質(zhì)對象間關(guān)系的表格, 拓撲表存儲層位間上覆、下伏和交切(被斷層切割后地層的拓撲表達)等的地層學(xué)關(guān)系及地質(zhì)空間位置關(guān)系。拓撲也可視為允許這些地質(zhì)關(guān)系合理儲存的數據結構。例如, 考慮多層地層, 上一個(gè)巖層的底面和與其相鄰的下一個(gè)巖層的頂面是上下巖層這兩個(gè)實(shí)體的公共部分或共享邊界, 它們之間的拓撲關(guān)系就是相鄰和同一的關(guān)系, 在存儲數據時(shí)只存儲上一個(gè)巖層的底面或其相鄰的下一個(gè)巖層的頂面, 即相鄰巖層的邊界曲面可以成為一個(gè)地層曲面,大大減少數據存儲量。評價(jià)地質(zhì)模型系統的優(yōu)缺點(diǎn)往往決定于描述地質(zhì)對象所用的拓撲結構[ 4] 。

 

  1 .5  可視化技術(shù)

 

  工程地質(zhì)復雜地質(zhì)體可視化, 是利用計算機技術(shù)將工程勘測獲得的數據, 轉換為形象直觀(guān)的便于進(jìn)行交互分析的地下地質(zhì)結構空間形態(tài)的立體圖和剖面圖形, 其基礎是工程數據和測量數據的可視化[ 5] 。利用可視化技術(shù)可以從龐大的地質(zhì)勘測數據中構造出地質(zhì)工程中對于邊破穩定性和地下硐室變形破壞等起關(guān)鍵作用的巖層和結構面, 并顯示其范圍、走向和相互交切關(guān)系, 幫助工程地質(zhì)人員對原始數據做出正確解釋, 繼而為工程地質(zhì)分析具體問(wèn)題提供決策支持。

 

  2  工程地質(zhì)三維可視化技術(shù)的初步開(kāi)發(fā)與應用2 .1  研究框圖

 

  工程地質(zhì)復雜地質(zhì)體三維建模與可視化的研究框圖如圖1所示。

 

  基于離散采樣數據的插值與擬合的思想, 即將離散數據轉化為連續曲線(xiàn)曲面, 工程地質(zhì)復雜地質(zhì)體三維建模與可視化的過(guò)程, 是從勘探數據庫中提取各種地質(zhì)信息的坐標位置及巖土體的物理力學(xué)參數, 通過(guò)不同的擬合與插值函數得到地質(zhì)層面(曲面)和地質(zhì)實(shí)體的三維計算機圖形顯示, 表達地質(zhì)信息在研究區域內的分布規律。生成地質(zhì)巖層面和地質(zhì)實(shí)體后, 實(shí)現從任意角度觀(guān)察建立的模型, 實(shí)現根據指定的剖面走向、傾向和傾角生成垂直剖面。

 

  2 .2  初步開(kāi)發(fā)與應用

 

  2 .2 .1  工程勘測空間數據庫管理

 

  在收集整理現場(chǎng)勘測數據后錄入金沙江某水電工程勘測空間數據庫各分項數據表, 這些數據表不僅包括地質(zhì)信息的位置數據, 更重要的是提供屬性數據。

 

  以地層巖性數據表為例, 要求錄入鉆孔編號、巖層起始深度、巖層終止深度、層厚、巖性(地層名稱(chēng))、地層代碼(地層年代)、巖層走向、巖層傾向、巖層傾角、接觸關(guān)系、地質(zhì)描述等數據。隨著(zhù)工程勘測的進(jìn)展, 能夠方便地修改補充和管理勘測數據。

 

  2 .2 .2  三維瀏覽

 

  通過(guò)孔口坐標和測量數據等的離散數據的擬合和插值法繪制壩址區的右岸地表曲面網(wǎng)格, 進(jìn)而可在三維圖形環(huán)境中進(jìn)行虛擬現實(shí)瀏覽觀(guān)察。

 

  2 .2 .3  三維地質(zhì)立體圖

 

  利用工程勘測數據, 建立了壩址區右岸三維立體地質(zhì)圖。該壩址區自上而下地層巖性組合為:第四系崩坡堆積物, 侏羅系泥巖、粉砂質(zhì)泥巖、泥質(zhì)粉砂巖, 三疊系上統厚至巨厚層狀細至中粒砂巖, 三疊系上統薄至中厚層狀粉細紗巖、粉砂巖, 三疊系上統中厚至厚層狀中粗砂巖。通過(guò)有限的工程勘測數據得出的立體圖,能夠較好地滿(mǎn)足工程地質(zhì)的精度。圖3 表達了該壩址區右岸三維地質(zhì)圖。

 

  2 .2 .4  三維可視化查詢(xún)

 

  通過(guò)圖形與工程勘測數據庫中的屬性數據的鏈接, 實(shí)現可視化查詢(xún)地層巖性和其他工程地質(zhì)信息, 最終完成向三維地質(zhì)信息系統的轉變。圖4 是一簡(jiǎn)單的被斷層錯斷的水平多層地層模型, 通過(guò)模型的每個(gè)地層實(shí)體名稱(chēng)與數據表中的巖石名稱(chēng)字段對應鏈接, 能夠查詢(xún)地層的巖性, 地質(zhì)年代, 起止深度和地質(zhì)描述等工程地質(zhì)人員關(guān)心的地質(zhì)信息。

 

  3  結論

 

  (1)運用先進(jìn)的可視化技術(shù)與交互圖形技術(shù)建立數據庫, 存儲和管理現場(chǎng)勘探實(shí)測和試驗數據, 建立工程地質(zhì)體的三維模型, 工程地質(zhì)工作者可隨著(zhù)勘察或研究工作的不斷深入細致, 對研究(工作)區域隨時(shí)補充信息來(lái)自動(dòng)顯示地質(zhì)信息在研究(工作)區域內的分布, 從而不斷提高模型精度, 并且利用模型反饋回來(lái)的信息及時(shí)發(fā)現已有勘察工作中的不足, 從而及時(shí)修改勘察或研究工作方案, 指導下一步勘探或研究工作的實(shí)施。

 

  (2)工程地質(zhì)三維建模與可視化的深入研究, 可以充分利用已有現場(chǎng)勘探實(shí)測或試驗數據, 達到節約投資減少勘察或研究成本的目的。當現場(chǎng)勘探和試驗數據資料不足情況下, 通過(guò)對已有數據的插值與擬合到建立三維模型, 可以推斷和預測未知區域或研究較少區域的地質(zhì)信息或巖土體物理力學(xué)參數的分布趨勢, 從而為減少勘探工作量提供科學(xué)的可靠的依據, 達到節約花費, 為生產(chǎn)或研究部門(mén)產(chǎn)生直接經(jīng)濟效益的目的。

 

  (3)工程地質(zhì)巖土體是復雜的不規則形體, 存在各種地質(zhì)巖性層面、結構面以及各種空間分布的地質(zhì)與力學(xué)信息, 完全表達地質(zhì)信息的空間分布及巖層和結構面間的位置關(guān)系, 工程地質(zhì)三維建模與可視化研究是大有作為的。