姜靈芝 應曉菲

（廣東省佛山地質局，廣東 佛山 528000）

1 引言

隨著城市化進程的加快和城市經濟的快速發(fā)展，一方面，城市建設導致了地面沉降、地面塌陷、地下水資源緊缺、土壤和地下水污染等各種城市地質問題，成為制約和影響城市可持續(xù)與安全發(fā)展的主要因素之一；另一方面，當城市建筑越來越多，位置要求越來越精確，地下地質信息的完備和準確，對于城市的發(fā)展規(guī)劃、設計、施工以及地下空間開發(fā)利用等方面的作用至關重要［1］。

完成城市地質的數(shù)據庫建設和三維建模，提供豐富的三維地質模型和二三維空間分析功能，二維、三維地質成果有機結合，地下、地表、空中一體化展示，建成面向專業(yè)研究的基礎工作平臺、面向政府規(guī)劃管理的立體可視化決策支撐平臺、面向社會公眾的地質信息共享服務平臺［2］，能很好地解決城市發(fā)展帶來的地質方面問題，能有力促進城市的發(fā)展。

2 三維地質建模方法研究

2.1 建模數(shù)據的處理

三維建模前期數(shù)據包括各類地質調查數(shù)據、鉆孔數(shù)據、試驗測試數(shù)據、地下空間調查數(shù)據等。數(shù)據整理工作包括鉆孔數(shù)據整理入庫、屬性數(shù)據入庫、鉆孔時代、分層、巖性等標準化、鉆孔坐標的標準化等工作。

（1）鉆孔數(shù)據錄入 ：研發(fā)相應軟件對鉆孔數(shù)據（工程信息、鉆孔信息、分層信息、實驗測試數(shù)據）進行錄入，整理并錄入工作區(qū)內鉆孔（參與三維建模 ）。

（2）統(tǒng)一坐標系 ：由于收集的鉆孔資料在坐標系、高程系、數(shù)據類型方面各不相同，需要將不同種類的坐標系統(tǒng)一到同一個坐標系。

（3）將鉆孔時代成因、巖土名稱、巖土描述標準化。

2.2 多源數(shù)據耦合層狀地質體建模

多源數(shù)據耦合層狀地質體快速建模的主要思想是集鉆孔、剖面、平面圖等各種多源數(shù)據于一體，充分利用已知的數(shù)據源，結合地質學的規(guī)律和數(shù)學算法，模擬位置區(qū)域的地質數(shù)據，從而建立整個區(qū)域的地質模型。該方法不僅大大降低了半自動建模的繁瑣程度，節(jié)省了半自動建模時間，還兼顧了模型的準確性，可大大提高建模效率。

以佛山南海區(qū)為例，該地區(qū)地層結構較復雜，存在很多山區(qū)溶洞等多種復雜地質體。根據這個結構特點，選用北京航空大學組織開發(fā)的深探地學建模軟件。該軟件具有原始剖面導入、鉆孔管理、斷面層建模、地層建模、透鏡體建模、斷塊建模等功能，優(yōu)勢在于可以處理任意復雜的斷層、地層不整合面，以及嵌入體等復雜地質情況，使用戶避免因建模軟件技術能力不足，而簡化價格不菲的地質資料來建三維地質模型的弊端。此外，所建三維地質模型可以無簡化地進行網格剖分，使屬性插值結果更加準確。

結合該軟件特點，針對工程地質、水文地質等簡單層狀地質體，采用“鉆孔-剖面/等值線-地層實體”構模的整體建模思路，利用所有地層界面共享的網格模板來構建各個地層面，再根據建模范圍和精度（網格間距）要求生成地形網格，從基礎數(shù)據庫中提取鉆孔點位和分層信息，疊加等值線數(shù)據，生成地層面強約束點；從剖面中提取有關地層邊界線信息，基于地形網格應用這兩類數(shù)據進行插值計算，構造各地層面模型，最后根據地層之間的疊覆關系等地質信息生成地層實體模型［3］。同時，對于地表模型可添加地形約束，構建出真實地形地貌單元的地質模型。地質模型建立后可以不斷添加各種約束數(shù)據，指定約束數(shù)據的影響范圍，對地質模 型進行反復重構更新，從而更精確地表達真實的地質形態(tài)。多源數(shù)據耦合建模思路如圖1 所示。

圖1 多源數(shù)據耦合建模思路

多源數(shù)據耦合建模是指用各種來源的數(shù)據進行地質建模，包括地球物理數(shù)據、地球化學數(shù)據、鉆孔數(shù)據、工程地質剖面、高程點數(shù)據、基巖埋深等值線數(shù)據、基巖出露區(qū)、等高線數(shù)據等。

2.3 斷層建模

一般情況下，斷層數(shù)據的來源有兩類，一類來源于地震解釋系統(tǒng)，斷層數(shù)據為整個斷層面；另一類來源于其他建模軟件，如：Petrel、RMS 等，數(shù)據并非全部斷層，而是目的層范圍內的斷層數(shù)據。

在建立斷層模型之初，先根據斷層數(shù)據情況判斷一個最小二乘面，定位斷面位置。所謂的最小二乘面是滿足所有數(shù)據點到該面的距離之和最短。

對于整體斷層數(shù)據來講，最小二乘面和實際斷面傾向基本吻合，而目的層斷層數(shù)據的分布往往呈條帶狀，這種情況計算的最小二乘面往往是水平面（如圖2a 所示），而非斷層傾向方向的最小二乘面（如圖2b所示），需要人機交互進行編輯。

圖2 整體斷層數(shù)據和局部斷層數(shù)據

2.4 地層建模

以斷層模型為約束，建立地層模型用于描述地質體界面形態(tài)、相互關系。

1)船舶減速會增加船舶占用航道的時間，后續(xù)船舶受船舶隊列順序的制約，造成船舶通航時間的延誤，航道的利用率下降。

（1）地層的生成

地層的生成有四種方法 ；最小曲率法、巖丘圈閉法、反距離加權法和克里金方法，各種方法的適用情況和特點有所不同，建模時應根據實際情況合理選擇。

（2）處理地層不整合接觸（尖滅、超覆）

地層之間如果存在不整合接觸，往往需要操作者根據地質認識處理地層形態(tài)，得到正確的地層體模型。軟件是通過人機交互指定不整合地層之間的主輔關系，用主地層裁剪輔地層網格，達到表達不整合層面的目的。因此，正確指定地層之間的主輔關系，對地層劃分空間和正確的網格裁剪起決定作用。

地層劃分空間的原理：工作區(qū)中包含H1、H2 和H3 三套地層，構造特征如圖3 所示，H2 尖滅于H1，H3 與H2 整合接觸，H2 為H1 的輔地層。

圖3 工作區(qū)中的地層及地層體

先判斷是否存在地層主輔關系，主地層優(yōu)先劃分空間。H1 地層先將工區(qū)上下一分為二，將a 空間與b1、b2、c 和d 分隔開來；由于H2 屬于H1 輔地層，暫不考慮，H3在H1地層空間劃分的基礎上進一步劃分，將b1、b2、c 和d 空間上下一分為二，將d 空間與b1、b2 和c 空間分隔開；最后考慮H2 地層，對b1、b2 和c 空間劃分，將c 空間與b1 和b2 空間分隔開。根據這樣的空間分隔原理，生成地層體，得到a 空間為地表的地層體，b1 和b2 空間為H1 的地層體，c 空間為H2 的地層體，d 空間為H3 與基底的地層體。

指定地層主輔關系是地層存在尖滅、超覆、侵入巖等不整合情況時需要處理的，通過指定地層主輔關系，用軟件正確裁剪輔地層網格面并進行空間劃分，如圖4 所示。

圖4 指定地層主輔關系

3 三維地質模型構建

3.1 結構建模

先導入所有建模需要的數(shù)據，根據數(shù)據設置好工區(qū)范圍、網格大小等，創(chuàng)建好初始模型，然后根據導入的數(shù)據生成斷層模型，再在斷層模型的約束下依據地層數(shù)據生成地層模型，最后生成地層體完成初始模型的建立。初始模型生成之后需要進行檢驗和校正，檢驗其是否與原始數(shù)據相吻合，是否符合地質規(guī)律。如果檢驗時發(fā)現(xiàn)問題，需要回到斷層和地層建模環(huán)節(jié)，進行迭代更新。

通過鉆孔數(shù)據、剖面數(shù)據、地質圖數(shù)據、DEM 數(shù)據、遙感影像數(shù)據進行三維地質結構模型的建設，具體工作流程如圖5 所示。

結構模型結合剖面數(shù)據和地質圖數(shù)據進行建模工作，模型按地質時代劃分地層，并按新舊順序排列。斷層建模前需要先導入斷層數(shù)據，一般情況導入的斷層數(shù)據為斷層格式。如果斷棱數(shù)據沒有進行拆分和命名，則需要先對斷層進行拆分和命名，然后軟件根據斷層數(shù)據范圍自動計算斷層的邊界，一鍵式生成全部斷面。初始生成的斷面由于異常數(shù)據的影響，一般會形態(tài)異常，需要人機交互通過編輯斷棱、編輯邊界等操作進行修正，再通過指定相交斷層之間的關系，做斷層削截處理。通過剖面數(shù)據對地層和斷層進行約束，剖面數(shù)據經過處理后，在三維模型里顯示。

圖5 城市地質三維結構模型建設流程

地層建模前需要先導入地層數(shù)據，地層數(shù)據往往是離散點格式，每個地層會在地層集下掛出相應的節(jié)點。三維建模軟件會在斷層的約束下生成地層面，層面的范圍為建模初設的模型邊界范圍。初始生成的地層面可能會有局部異常，多是地層數(shù)據問題或斷層形態(tài)異常引起，需要使用恰當?shù)膶Σ呓鉀Q相應的問題。完成層面異常修正后，需要對地層節(jié)點進行排序，使地層節(jié)點的順序與實際地層的沉積次序相符，否則和沉積次序不符的地層體將無法生成。如果存在不整合接觸的地層，還需要定義地層之間的主輔關系，做地層削截處理，以上工作全部完成后，可以一鍵式生成地層體。

3.2 網格剖分

在任意指定的空間分辨率條件下，基于矩形截斷網格對矢量結構模型執(zhí)行空間剖分，得到精確表達結構模型邊界形態(tài)的網格模型。超高的空間分辨率可能得到超大規(guī)模的網格數(shù)據。網格生成的基本原理是首先生成二維平面網格，然后用二維網格對三維地層體進行縱向剖切。

構造模型的和網格模型的模型網格化主要是基于應用的需求。在模型從頂?shù)降讛?shù)據精度不同時，可以采用多分辨率的網格；在結構模型和無結構模型下，采用不同類型的網格剖分；同一模型在使用不同的屬性進行插值時，可根據數(shù)據的需要剖分成不同大小的網格單元（LOD 網格）；模型工區(qū)非常大，模型精度要求非常高時，海量的網格單機無法承載，可分區(qū)塊剖分，邏輯上是整體，物理存儲分塊，如圖6 所示。

圖6 大工區(qū)構造模型及大工區(qū)網格模型

數(shù)據或認識發(fā)生改變后，結構模型可以快速更新，在更新后的結構模型基礎上，重新生成網格，網格參數(shù)可直接引用原數(shù)據，也可根據更新后的結構模型進行適當更改，也就是說網格化是基于數(shù)據驅動，而不直接編輯網格，如圖7 所示。

圖7 更新模型

3.3 屬性建模

屬性參數(shù)建模是將模擬區(qū)域三維網格化，然后在網格中填充屬性值。對于任意一種儲層物性，每一個網格有且只能有一個值。由于網格精度通常要遠遠低于測井曲線的采樣精度，所以測井曲線數(shù)據必須要進行粗化才能采樣到網格中。數(shù)據粗化是屬性建模的前提和數(shù)據基礎。

圖8 利用鉆孔最小曲率算法生成地表面及插值后的巖性屬性模型

屬性建模是基于構造模型通過地質統(tǒng)計學方法計算生成的屬性地質模型，包括離散屬性（沉積相和巖相）或者連續(xù)屬性（含水量、重度、滲透率和飽和度等）。地質領域比較常用的插值算法是克里金方法和序貫模擬等方法。巖性屬性模型利用多點地質統(tǒng)計學，采用序貫指示的插值算法生成屬性模型，插值后的巖性屬性模型如圖8 所示。

4 應用實例

通過上述方法和步驟，運用鉆孔數(shù)據、剖面數(shù)據、地質圖數(shù)據、DEM 數(shù)據、遙感影像數(shù)據對佛山市地質結構進行建模，較好地完成了佛山市城市地質三維模型的構建，成果如圖9 所示。

圖9 佛山市地質結構模型與地質圖對比效果

5 總結

本文結合佛山市城市地質項目工作的安排，通過分析現(xiàn)有資料情況，參考國家行業(yè)地方等標準，并結合佛山地質數(shù)據特點，對接城市地質綜合服務系統(tǒng)，通過對各種多源數(shù)據的整理和標準化，建成以基礎地質、礦產、水文地質、工程地質、環(huán)境地質、地球物理、地球化學、遙感數(shù)據為核心數(shù)據的三維空間數(shù)據庫和元數(shù)據庫。此外，還以該數(shù)據庫為基礎數(shù)據，建立了佛山市和三龍灣城市地質三維結構和屬性一體化的模型，并能加載融合地下管線、隧道、基礎等構筑物的三維模型，成功完成了多源數(shù)據融合的三維地質模型構建的研究，同時為更加安全合理開發(fā)利用城市地下空間提供了科學依據。