基于GIM 技術的城市地質體模型構建及應用

陶二永，門陽陽，董云奇

（1.山東省煤田地質局第四勘探隊，山東 濰坊 261200；2.山東交通職業學院，山東 濰坊 261206）

引言

近年來，建筑信息模型BIM（Building Information Modeling）技術被廣泛應用在建筑、結構、施工等工程建設領域［1-4］。由于地質體的不規則性和不確定性，BIM 技術難以在巖土工程領域應用，地質信息模型GIM（Geological Information Modeling）技術由此誕生。GIM 的核心內容是描述已有地質對象的幾何形態和工程特性，主要用于構建地下三維可視化勘察信息模型。

隨著我國社會經濟的迅速發展和城市化進程的全面推進，城市地下空間作為城市空間資源極具開發潛力的部分，政府及相關專家對地下空間給予高度重視［5-7］。目前，有關地下三維可視化研究的方法尚處于摸索階段。張照杰等［8］運用三維激光掃描技術獲取點云數據對地下停車場進行空間三維建模；鄭國棟［9］研究了巖土工程勘察信息模型系統對樁基礎進行了三維可視化和巖土參數化分析；施炎等［10］研究了基于地質信息模型（GIM）的灌漿工程三維可視化分析方法和應用。然而目前針對城市地層分布模型的可視化研究較少。

1 GIM 技術

GIM 技術面向巖土工程地質，運用計算機圖形學技術，將工程勘察數據信息通過光滑離散插值技術（Discrete Smooth Interpolation，DSI）［11］構建三維可視化地質模型，其表達方式更直觀，可完整展示地質信息，并且可以獲取任意剖切面的地質情況，也可以選擇任意一層地層查看地質信息。

GIM 軟件為BM_GeoModeler，采用模塊化方式系統整合勘察、建模、設計與分析等巖土工程全流程應用功能的獨立平臺，其中的建模方法包括針對土體和平緩巖層獨創的一鍵建模、土層分類建模、巖性建模等批量自動建模方法，允許人工干預的流程建模、物探建模、針對復雜透鏡體的透鏡體合并等。

2 城市地質體GIM 建模方法

地質建模是描繪出一個已經存在的自然對象，且可供建模的資料往往很少，采用離散光滑插值（DSI）理論和約束技術模擬其形態特征的中心任務是利用這些少數已知數據、合理地 “推測” 出其余部分的形態，其中的其余部分具有不確定性，是地質體不確定性的典型代表。基于上述特點，地質體建模基于離散數學理論，BM_GeoModeler 采用DSI 理論進行建模和數據（包括物探指標、礦山資源化驗結果、土體靜探或試驗值等）處理，可處理任意不規則面的局部修正、彼此相交（針對地層尖滅、斷層錯動、透鏡體）等。

離散光滑插值算法理論是基于離散數學的一種插值理論。與傳統的建模算法相比，它可以根據各種約束條件及其組合擬合構建復雜如多Z 值（褶皺、透鏡體）、非連續性（斷層、覆蓋層）地質體模型，并且能進行局部修改，不必因地質勘探資料的變化而重新建模。（如果將地質界面視為離散化的不連續界面，地質點及地質勘探揭示的鉆孔平硐數據等作為約束條件，DSI 實際上就是通過在這些約束條件下求解目標函數—全局粗糙度函數的最優解來得到符合約束條件的最優化地質界面）。

根據實際約束情況可以得到不同條件下的約束系數，進而通過上述公式迭代求解最優化的Φ值，最終擬合得到符合約束條件的幾何模型。DSI 的約束可以分為軟約束和硬約束，軟約束通常指條件放松的（約等于）通過全局最小平方和進行擬合的約束；硬約束則是必須完全通過等式或不等式擬合的約束。非線性約束可以通過泰勒公式轉化為近似的線性約束進行擬合。由于DSI 方法考慮了節點與鄰域節點之間的關系，可以比較好地擬合非連續性幾何模型（比如地質上的斷層上下兩盤）。另外，DSI 可以根據實際情況的約束條件擬合非常復雜的模型，因此，特別適合復雜的三維地質建模。

城市地質體三維模型構建流程見圖1。主要步驟：（1）通過地質勘測獲取鉆孔數據，整理使用鉆孔共計328 個，其中318 個鉆孔深度達30 m，10 個鉆孔深度達到50 m，并針對各鉆孔的分層情況整合出所有土層類型，得到鉆孔數據錄入的巖性定義和地層定義。（2）鉆孔數據在導入BM_GeoModeler 數據庫后，可重復調用建立多個項目，便于項目各階段的修改。為確保導入的鉆孔數據在地層建模分層時不出現錯誤，可在鉆孔數據庫建立后進行層序檢查，根據檢查結果修改錯誤地層數據，直至檢查無任何錯誤方可使用。（3）地表模型的構建以鉆孔孔口標高為主、高程點數據為輔經樣條插值擬合而成。鉆孔孔口標高可從鉆孔數據庫中提取得到孔口標高點集，修改其屬性為地表，地質代號選為top；高程點數據需將高程圖轉化為.dxf 格式導入項目。（4）建立地層模型，地層模型主要參考的數據是鉆孔的豎向分層數據，是對各鉆孔間的地層根據相鄰鉆孔的分層情況做出的推測。地層模型與實際地層分布有一定誤差，因此地層模型分層的精度取決于鉆孔數據的準確性及鉆孔的密度。

圖1 城市地質體三維模型構建流程

3 城市地質體GIM 建模結果及應用

3.1 城市地質體GIM 建模結果

將孔口標高點集和高程點全局顯示后確定地表建模范圍，該范圍應將所有參考點位置包含在內。為確保在劃分網格時出現過小線段單元從而影響網格劃分進度，BM_GeoModeler 提供線集工具對高程點數據進行過濾，在地表建模范圍內進行網格劃分，網格的大小綜合考慮工成需求及擬合時間確定，地表建模方法采用單一界面建模，在網格全局加密的同時進行擬合計算，地表模型擬合結果見圖2。

圖2 地表擬合結果

地層模型建模步驟：建立地層模型邊界→劃分三角網格→檢查與調整地層→分層計算。分層計算方法采用離散光滑插值算法理論（DSI 技術）計算。

為準確定位地層模型上各位置，可對模型進行表面貼材質處理，即通過表面貼材質功能將該地區的地圖貼于模型表面。表面貼材質貼完材質后的地層模型見圖3。該地區地下管線錯綜復雜，包括各類電信管線、電力管線、燃氣管線、排水管線、熱力管線，地下管線成果以《城市地下管線探測技術規程》（CJJ 61—2017）為主要標準，設置管線數據代號及顏色。

圖3 地層模型

BM_GeoModeler 支持將管線數據以.dxf 格式保留三維位置信息導入至地層模型中，并自動進行布爾操作嵌入地層模型，地層模型加入地下管線后見圖4。

圖4 地下管線導入地層模型

3.2 城市地質體模型應用

運用BM_GeoModeler 進行剖面檢查與出圖特點：（1）內置的三維地層拓撲算法，從三維建模的角度識別地層連接關系、計算地層尖滅位置、識別透鏡體。（2）剖面內允許地質人員的檢查與干預調整、編輯與調整在剖面間聯動作用，可以在智能化識別的基礎上，實現局部位置的人為認識干預，并實時更新因人工編輯造成其他相關剖面地層連線的變動，確保各二維剖面成果在三維空間中的正確性。（3）引入了不同地質成因分界面尖滅機制，實現如沖洪積地層可以平緩尖滅于殘積層或基巖的特征，使得地質剖面繪制更加符合地質成因規律。

在任意剖面下，可根據地質判斷，調整程序智能判斷給出的地層連線結果，調整過程中，相關剖面會跟著聯動更新，降低了調整修正工作量，同時也時刻保證剖面間的編輯不影響三維空間中地層拓撲關系的正確性。

選取8 條主要干道作為剖面見圖5。在BM_GeoModeler 中可通過創建非封閉折線生成任意剖面線位置，再通過生成剖面形成該剖面線位置的剖面見圖6。可以看出研究區地勢整體為西低東高，南低北高。個別剖面圖中有透鏡體出現，考慮地層分層定義時主層有多個亞層，部分透鏡體為主層的一個亞層形成，其地層巖性與其相鄰的主層相近。

圖5 研究區主干道位置

圖6 研究區主干道剖面

4 結語

（1）選取鉆孔孔口標高作為地表建模的主要參考點，再以高程點為次要參考數據，根據數據所在平面區域確定地表建模范圍，再以立方網劃分網格，從而擬合計算得出地表模型。（2）對鉆孔數據進行層序檢查，采用離散光滑插值算法理論（DSI）進行分層計算，并劃定地層建模范圍，再根據鉆孔數據構建了地層模型。（3）將構建的三維地質體模型應用于工勘項目常用的地質剖面圖的繪制，并分析了該城區地層結構和地勢走向。