基于GOCAD的寧蕪盆地云臺山地區三維地質建模

，， ，

(1.江蘇省地質調查研究院，江蘇南京210018； 2.北京大學地球與空間科學學院，北京100871)

基于GOCAD的寧蕪盆地云臺山地區三維地質建模

趙增玉1,2，陳火根1，潘懋2，賈根1

(1.江蘇省地質調查研究院，江蘇南京210018； 2.北京大學地球與空間科學學院，北京100871)

闡述了利用GOCAD軟件進行三維地質模型構建的思路及主要步驟，重點介紹其模擬地質界面的3種主要方法。以寧蕪盆地云臺山地區為例，利用GOCAD模擬了該區地層、巖體及斷裂。在構建地形、地層及構造面時，利用剖面建模法以及離散光滑插值技術，快速構建研究區的表面模型。由面模型向體模型轉化時，選用GOCAD的SGrid代替Solid功能，更加方便準確地反應了真實地質情況。

GOCAD軟件；三維地質模型；SGrid；云臺山；寧蕪盆地;江蘇；安徽

0 引 言

從20世紀80年代開始，各國學者(Carlson，1987；Houlding，1994；Wu et al，2005；Turner，2006)相繼進行三維地質建模理論與方法的研究，其中，三維空間數據模型是三維地質建模的重要研究內容，到目前為止，三維空間數據模型的種類達數十種，其中代表性的有30余種。一般地， 將三維空間數據模型分為基于面表示的數據模型、基于體表示的數據模型和混合數據模型(Gong et al, 2004)。如最為常見的規則網格(Grid)、不規則三角網(TIN)就是基于面表面的數據模型。

隨著三維地質建模技術的發展和廣泛應用，產生了很多非常實用的建模軟件，比較有影響的有GOCAD、Geomodeller3D、GSI3D、Earthvision、Vulcan、Micromine、GSIS、3Dgrid，其中，法國Nancy理工大學研發的GOCAD(Geological Object Computer Aided Design)發展較為成熟且應用廣泛。近年來，我國很多專家學者將GOCAD應用于石油、地質災害等領域，取得了良好效果。以寧蕪云臺山地區為例，利用GOCAD進行研究區的三維地質模型構建。

1 GOCAD軟件概述

GOCAD地質建模軟件是主要應用于地質領域的三維可視化建模軟件，在地質工程、地球物理勘探、礦業開發和水利工程中均有廣泛的應用，是國際上公認的主流建模軟件，該軟件充分考慮了地質資料的多源性，提供了多種建模方法，既可建立表面模型，表達空間幾何對象，也可建立實體模型，展現空間屬性分布，空間分析功能強大。三維地質建模軟件實現的關鍵技術如下。

(1) 開發了針對地質建模特點的空間插值算法——離散光滑插值算法(Discrete Smooth Interpolation, DSI)。DSI用一系列具有物體幾何與物理特性的相互連接的節點來模擬地質體，用DSI方法模擬幾何與物理特性時，已知節點和地質學中的典型信息將被轉化為線形約束，引入到模型生成的過程中。該方法基于圖形拓撲，它適用于構建復雜模型和處理模型表面不連續的情況。DSI算法可自由選擇和自動調整格網模型，實時交互操作，能夠處理一些不確定的數據等特點，決定其在地質建模和可視化中占有重要的位置(董梅等，2008；魏群等，2008；孫濤等，2011)。

(2) 實現Workflow半自動化建模，提供了構造模型(Structural Modeling)建模和三維儲層柵格模型(3D Reservoir Grid Construction)建模。構造建模可建立地質體構造模型，模擬地層面或斷層面的空間展布形態、位置和相互關系，同時，構造建模還是地震勘探過程中地震反演的重要手段；三維儲層柵格建模可根據建立的構造模型，將其轉換為實體模型，并可對地質體內部屬性(如品位)進行計算和綜合分析，形成地質體的屬性模型。

2 GOCAD模擬地質界面的主要方法

地層、斷層建模是三維地質建模最核心的部分，這些面存在于地質體內部，無法直接觀察到，GOCAD提供多種方法進行地層、斷層的模擬。

2.1 鉆孔分層點建模

若研究區通過鉆孔等手段獲取了大量地質信息，并且地層構造簡單，無復雜地質體，此時建立層面、構造面模型相對簡單，只需將同一層面或結構面上的數據以X、Y、Z的形式整理成文本文件并導入GOCAD，然后利用Surface—New—From PointSets命令生成層面或結構面。這種建模方法自動化程度高，適用于簡單的工程地質層建模，缺點是無法處理斷層。然而，在大多數情況下建模所需要的地質信息是通過野外踏勘和少量工程獲得的，并且地質背景復雜，無法直接通過鉆孔分層點建立全區的三維模型。

2.2 產狀建模

根據地層和構造的地表行跡和產狀模擬地表以下較淺范圍內的地質體。GOCAD具體實現時，可將地層或構造面產狀轉換為面的切向量。若X軸正向為東，Y軸正向為北，向上為Z方向，某層面傾向為α，傾角為β，則該層面切向量為X=cosβ·sinα，Y=cosβ·cosα，Z=-sinβ。根據層面的地表形跡以及切向量，可基本建立該層面的三維空間形態(圖1)。

圖1 根據層面產狀建模示例

圖1a為某層面或構造面的地表形跡及其產狀，圖b為根據產狀建模后在三維視圖中的正北側視圖。首先，將地表Curve對象沿面切向量拉伸一定的距離得到一個面對象，即為該地層面的宏觀基本形態，該面對象與地表出露線及產狀一致。然后，根據鉆孔揭露的該層面或構造面的位置，設置控制點，將面對象擬合到該控制點，得到目標曲面，若其空間展布與地質規律不矛盾，則可認為該曲面可以模擬該地層面，既符合地表測量的產狀數據，又經過鉆孔所控制的空間位置。重復以上過程，對每一個地層面、斷層面分別建模，就得到整個建模區域的巖性、構造分布情況。

根據地質圖產狀建模原理，可在層面產狀變化小、地質背景簡單、建模精度要求不大的區域采用該方法建模較為快捷，并且該方法還經常用于斷層三維模型的構建。

2.3 剖面建模

利用GOCAD的CrossSection對象可以很方便地將鉆孔投影到自定義剖面上，基于剖面上的鉆孔分層信息，并結合地質背景和地表出露情況，對剖面圖進行人工解譯。首先確定1組用于生成剖面的鉆孔，然后執行CrossSection—New—From Digitalized Polyline命令，在視圖中編輯自定義勘探線，使得勘探線盡可能地通過鉆孔，然后輸入要創建剖面的名稱，設置剖面的基本信息，如剖面延伸深度，即可生成自定義鉆孔剖面圖，且同時生成地表地質界線與剖面的交點，根據鉆孔與地表出露情況，可解譯各地層、巖體、構造在剖面上的走向，在剖面地質解譯完成后，對象視圖中可以看到剖面在三維模型中的空間展布。

此外，剖面控制程度需要達到地質建模所要求的精度，當解譯完每個剖面后，通過GOCAD 中Surface—New—From Several Curves命令，選擇同一層面或構造面的不同地質界線，生成對應層面或構造面，即生成了各地層、構造面的三維模型。在實際建模工作中，地層、構造的交切關系相當復雜，這就需要對交切部位進行面的局部調整。GOCAD具有強大的編輯功能，可以通過增加約束、局部DSI等命令，通過調整三角網的控制點來對曲面進行局部編輯。同時，由于生成的地層或構造面的范圍不一定與研究區范圍一致，因此，為了將所有地層、構造約束在研究范圍內，生成規則的三維模型，必須建立統一的約束，使地層、構造三維模型限制在同一范圍內。在GOCAD中選擇研究區域的邊界線，通過Surface—New—Built In Forms—Tube命令建立研究區邊界線的垂向切面，然后利用Tools—Border—Extend命令擴展沒有靠近邊界面的層面或構造面，利用Cut By Surfaces命令刪除超出邊界面的層面或構造面三角網。

3 應用實例

云臺山地區位于寧蕪盆地中段的東部區域，區內主要礦種為黃鐵礦，包括鶯子山、云臺山、天臺山、富而崗、公雞山、母雞山等多個礦段，建模范圍面積約35 km2，深度為850 m。對區內的主要地層、構造、次火山巖體進行三維空間展布分析，以進一步探索三維立體填圖技術方法、成果表達方式等。

3.1 地形面模擬

云臺山地區目前有相應的1∶1萬地理底圖資料，包括MapGIS等高線數據，但沒有高程屬性值，需要根據高程標注進行人工賦值，利用Section軟件(張燕飛等，2011；東海宇，2012)對云臺山等高線MapGIS數據進行了高程賦值。然后將MapGIS等高線數據(*.wl)轉換為CAD格式(*.DXF)，接著導入GOCAD。

通過File—Important Objects—Cultural Data—DXF加載地形文件，GOCAD會自動分辨DXF中各圖層，并以不同文件名保存在Curve單元(PL文件)中，文件名一般與圖層名相同。勾選等高線所在的PL文件，在右側視圖窗口中查看三維等高線是否正確。然后在PointSet選項下，通過New—From PointsSet Curve or Surface命令將等高線離散，并生成PointSet的VS文件(圖2)。也可以通過MapGIS先將等高線離散，以X、Y、Z的坐標形式保存在文本文檔或Excel等文件中，然后直接導入GOCAD生成VS文件。

圖2 云臺山地區地形離散點分布圖(Z×3)

GOCAD常用3種方法生成地形曲面：(1) 根據離散點分布形式在Curve中擬合地形邊界線，擬合的精度根據需要調整。然后在Surface中通過邊界線和離散點以相似三角(Homogeneous triangles)網格構建曲面，即Surface—New—From PointsSet and Curve。這種方法直接通過原始點進行地形三角網的構建，沒有插值運算，得出的地形模型較為粗糙。(2) 根據Wizard—Surface Creation—From Data(Without Internal Borders)向導生成面，對于地形面與原始控制點吻合較差的區域可以通過增加相應的控制點，再進行局部光滑平順處理(DSI)，多次反復直到符合要求。具體方法為：選擇控制點及邊界線，通過Constrain All Borders On Straight Line將邊界線轉換為模糊控制點，引入插值方向，GOCAD會自動進行DSI插值。(3) 根據離散點生成中間平面，再以離散點為該平面控制點，進行DSI插值。

GOCAD既可以模擬三維地形面，也可以將地表出露的地層、構造等信息賦予地形面，生成三維地形地質面。云臺山地區地形面采用Wizard向導中插值法模擬，并將地表出露地層及構造行跡投影到地形面(圖3)。

圖3 云臺山地區三維地形地質圖(Z×3)

3.2 鉆孔的生成

鉆孔是三維地質模型構建的重要原始資料。在GOCAD中，垂直或傾斜鉆孔的加載方式與直井或斜井相同。但不同文件格式有不同的加載方法，在文本文件中，可以直接加載鉆孔名稱、鉆孔位置(X、Y、Z)、補芯海拔KB、孔深等信息。鉆孔生成后可以在其對應的maker選項中加載或錄入分層、斷裂信息，如地層名稱、厚度、產狀等。

結合研究區不同礦段探礦工程的資料，加載了CK202，CK206，ZK1001等46個鉆孔，并且每個鉆孔都生成了相應的maker，鉆孔三維空間分布見圖4。

圖4 鉆孔三維空間分布圖

同時，由于地形是等高線插值生成的，因此，鉆孔孔口位置與地形不一定吻合，這就需要根據鉆孔的孔口位置校正地形，GOCAD提供了點校正面的功能，即Constraints—Control Points—Set Control Points，可使鉆孔的孔口嚴格地與地形重合。

3.3 地層、斷層模擬

云臺山地區受方山—小丹陽、方山—南山等深大斷裂帶多期活動影響，地形地質條件復雜。區內地層由老到新為早中三疊世黃馬青組、早侏羅世鐘山組、中侏羅世北象山組沉積基底，晚侏羅世龍王山組和晚侏羅—早白堊世大王山組火山碎屑巖，早白堊世赤山組沉積碎屑巖。區內基底褶皺為寧蕪向斜，其向斜東南翼發育次級褶皺。寧蕪斷陷盆地由于受多期構造活動的影響，基底斷裂發育，主要有NNE、NE、NW 向3組，這些斷裂構成了寧蕪斷陷盆地北段網格狀構造輪廓，云臺山地區主要表現為北東向的2個大斷裂。

利用剖面建模方法建立了研究區的地層、斷層三維模型。首先建立11條貫穿研究區的北西向剖面，為彌補區調鉆孔深度的不足，結合研究區地質背景，對各剖面進行人機交互解譯，在解譯過程中，不僅要考慮到剖面上各鉆孔的實際分層及斷裂點情況，更要結合區域地質背景，使每條剖面上的斷層、地層均符合區域構造、地層的展布，只有這樣，建立的模型才更能直觀地反映地層構造格架。

在正確解譯各剖面之后，先根據各剖面上的斷層軌跡，由線生成面，建立斷層面的三維模型，再根據各剖面解譯的地層界線，生成三維地層面，最后處理地層與斷層之間的交切關系，在地層面與斷面相交處，進行面相交計算，切除多余地層，而在未相交處，對地層面進行趨勢延伸計算，直至相交到斷層面上。對每個地層與斷層進行交切處理后，建立圖5所示的云臺山結構面三維模型。

圖5 云臺山結構面三維模型

3.4 三維體模型創建

GOCAD可以通過SGrid(網格模型對象)和Solid(實體模型對象)2種方式建立三維體模型。其中，Solid地質實體模型對象可通過點集、封閉曲面或單個面及厚度來創建，這就決定了若要通過Solid方式生成較為精確的模型，則需要足夠詳細的數據支持，方法也較為復雜。SGrid是由六面體格子組成的單元集，確定地層分界面及網格密度后，GOCAD可以通過六面體單元“充填”兩層面間的地層，所有這些六面體單元組成SGrid。劃分的六面體單元越小，擬合效果越好，鋸齒現象越不明顯。此外，六面體單元的節點信息可以導出，有利于進一步對地質體進行空間分析及數值模擬。因此，選用SGrid進行三維體模型的構建。

在Workflow中新建3D Reservoir Grid Builder，GOCAD將以向導模式半自動化建立模型。主要步驟如下：(1) Select Horizons，確定模擬地層的范圍，單擊已經建立好的最新地層頂面和最老地層底面。(2) 在Specify Gridding中，可以選擇自動、交互或已有模式計算網格(Gridding)和地層厚度(Pillars)，無特殊要求可以選擇自動計算。(3) Define Units，定義各中間地層，并按新老順序排列。(4) 最后在Build Grid中將三維網格細化，減少鋸齒，生成體模型。通過該方法生成的體模型每一個單元都有獨立的節點坐標和屬性點，可賦予多個屬性值。云臺山地區的SGrid模型見圖6。

圖6 云臺山地區三維地質模型

4 結 論

GOCAD是一款強大的地質建模軟件。一方面可以識別多種數據來源，極大地簡化建模，而且在建模時能通過DSI插值技術建立表面模型，建立的三維模型更加逼近真實地質體。另一方面，隨著研究的深入，常需要在地表、鉆孔或坑道中取樣進行各種化學元素分析，把試驗所得數據以取樣點屬性的方式賦予三維模型中對應鉆孔上，然后將有限的實驗數據通過GOCAD軟件提供的克里格插值法，即可得到對應的屬性模型或變量圖，進而進行專項分析。

董梅，慎乃齊，胡輝，等. 2008.基于GOCAD的三維地質模型構建方法[J].桂林工學院學報，28(2):188-192.

東海宇. 2012. SECTION在地質數字化制圖中的應用[J].礦山測量，(4): 35-36.

孫濤，劉亮明，趙義來，等. 2011.基于GOCAD平臺的復雜地質體系的動力學建模研究[J].礦產與地質，25(2):163-167.

魏群，黨麗娟，張俊紅，等. 2008.GOCAD在巖體三維可視化模擬中的應用[J].煤田地質與勘探，36(5):15-19.

張燕飛，朱杰勇，張威.2011.基于GOCAD的三維地質模型構建[J].河北工程大學學報：自然科學版，28(4):69-73.

CARLSON E.1987.Three-dimensional conceptual modeling of subsurface structures[C]// ASPRS, ACSM Annual Convertion. Baltimore,Maryland, USA: ASPRS, 1(4):188-200.

GONG JIANYA, CHENG PENGGEN, WANG YANDONG. 2004.Three-dimensional modeling and application in geological exploration engineering[J].Computers & Geosciences,30(4):391-404.

HOULDING S W.1994.3D Geoscience Modeling: Computer Techniques for Geological Characterization[M]. London,UK: Springer-Verlag.

TURNER A K.2006.Challenges and trends for geological modelling and visualisation[J].Bulletin of Engineering Geology and the Environment,65(2):109-127.

WU QIANG, XU HUA, ZOU XUKAI.2005.An effective method for 3D geological modeling with multi-source data[J].Computers & Geosciences,31(1):35-43.

3D geological modeling of Yuntaishan area in Ningwu Basin based on GOCAD software

ZHAOZeng-yu1,2,CHENHuo-gen1,PANMao2,JIAGen1

(1. Geological Survey of Jiangsu Province, Nanjing 210018, China; 2. School of Earth and Space Sciences, Peking University, Beijing 100871, China)

The authors illustrated the ideas and main steps of constructing 3D geological models based on GOCAD with focuses on threes main methods simulating the geological surfaces. Taking Yuntaishan area in Nanjing-Wuhu Basin for example, the strata, the rock mass and the fault of the area were simulated by using GOCAD software. Surface model was finished based on cross-section and Discrete Smooth Interpolation technology (DSI). By choosing the function of SGrid instead of the function of Solid when changing from surface model to volume model, it showed the real geological condition more conveniently and accurately.

GOCAD software; 3D geological model; SGrid; Yuntaishan; Ningwu (Nanjing-Wuhu) Basin; Jiangsu; Anhui

10.3969/j.issn.1674-3636.2014.04.652

2014-03-19；

：2014-05-07；編輯：陸李萍

趙增玉(1984— )，女，工程師，博士，構造地質學專業，主要從事地質礦產調查評價工作，E-mail：zengyu1206@163.com

P628

：A

：1674-3636(2014)04-0652-05