立體地質填圖三維建模技術方法與應用研究

朱　威， 王大勇， 王書民， 姚大為， 王　剛

(中國地質科學院地球物理地球化學勘查研究所，廊坊　065000)

?

立體地質填圖三維建模技術方法與應用研究

朱威， 王大勇， 王書民， 姚大為， 王剛

(中國地質科學院地球物理地球化學勘查研究所，廊坊065000)

針對綜合地球物理立體地質填圖后期進行的三維地質模型構建、圖示展示方法與技術進行了研究。提出了不同屬性體單獨建模后進行切割、融合構建三維地質模型的復雜地質結構建模方法與技術；在MAPGIS K9平臺上應用該技術構建了九瑞試驗區(qū)三維地質模型，并通過不同的圖示方法分析與展示構建的九瑞試驗區(qū)地質模型；結合該區(qū)域的成礦特征利用三維地質模型從深部地質結構出發(fā)圈定了該區(qū)的成礦有利部位。

立體地質填圖； 三維圖示； 深部找礦； 九瑞試驗區(qū)

0　引言

傳統(tǒng)的地質填圖方法只能描述地表或淺部的地質信息，對深部地質信息的描述有著先天不足。立體地質填圖工作在國外早已開展。前蘇聯(lián)早在上世紀60年代～70年代便開展了立體地質填圖的研究工作，出版了有關1∶50 000地質測量的《稀有金屬礦區(qū)立體地質填圖》[1]等參考著作，并曾組織專門會議，研討距地表深2 000 m以內的礦產預測問題，提出在典型礦區(qū)開展立體模型研究，并強調在立體填圖的基礎上進行立體找礦。我國在油氣和煤炭勘查中，采用地震等物探方法與鉆探結合的勘查摸式，獲得深部立體地質構造圖，已是常規(guī)的地質工作。近年來也開始了城市三維地質調查的淺層地質結構的立體地質填圖試點工作，但甚少涉及深部金屬礦產的立體填圖工作。為提高我國資源的保障能力，配合“加強深部找礦，開辟第二找礦空間”的戰(zhàn)略舉措，許多專家、學者提出的思路和意見都涉及到立體地質填圖問題。

基于綜合地球物理方法的立體地質填圖，主要是利用獲取的綜合地球物理資料對地下一定深度范圍內具有明顯或較明顯的物性統(tǒng)計差異的物性-地質體(地質單元、目標地質體、填圖地質單元等)的分布與結構特征做出立體(三維)推測[1]，并在MAPGIS K9平臺上運用先進的建模方法和技術繪制成推斷立體地質圖。通過一系列的展示和分析方法對繪制的模型進行展示和分析，解決傳統(tǒng)填圖方法不能對深部地質信息取得更多準確認識的問題。

1　三維地質模型的構建

1.1構建約束剖面

約束剖面是由綜合地球物理解釋推斷解譯的地質剖面而來，在三維地質建模中是作為對填圖單元的空間分布進行約束的已知信息[2]，約束剖面中通常含有物性單元層、構造、地形等信息。依據三維地質建模技術，在構建約束剖面時，需要依據三維建模平臺所要求的剖面數據類型，利用矢量化技術，構建不同屬性(地層屬性剖面、巖體屬性剖面、構造屬性剖面等約束剖面，圖1)。矢量化過程中采取簡化剖面矢量化的方法(圖1中綜合地球物理解譯剖面中包含地層、巖體及構造等多種屬性)，矢量化地層屬性剖面時，將巖體屬性部分合并到相應地層屬性(圖1中綜合地球物理解譯剖面、地層屬性剖面圈出部位)，讓同一屬性地層更為連續(xù)，便于三維模型的構建，然后，將不同屬性的剖面分別導入三維建模平臺，構建不同屬性的立體地質結構模型，圖2為導入后的地層屬性剖面和巖體屬性剖面。

圖1　矢量化后地層、巖體、構造屬性剖面圖Fig.1　Vectorisation of stratigraphic, rock and tectonic attributes(a) 綜合地球物理解釋剖面；(b) 地層屬性剖面；(c) 巖體屬性剖面；(d) 構造屬性剖面

圖2　三維建模各屬性剖面數據導入流程示意圖Fig.2　The import process sketch of attribute profile used for 3D modeling(a)地層屬性剖面導入后陣列顯示；(b)巖體屬性剖面導入后陣列顯示

1.2同屬性剖面構建同一屬性體

屬性剖面到屬性體采用基于輪廓線拼接的交互建模方法，其基本思路是從兩個鄰近的剖面中選擇同一地層的邊界線進行拼接來構建地質體的面模型，再由面構成地質塊體，最后由地質塊體組合構成完整的地質實體模型。在上述過程中，關鍵步驟是基于輪廓線數據拼接，即找到一個“最佳”的與所給定的輪廓對應的曲面，選擇一個與實際地質體對象最為接近的匹配。輪廓線拼接處理的核心是用將相鄰兩斷面的輪廓線連接成由多個三角面組成的輔助面(圖3(b))，再在兩個平面輪廓線間用多個輔助面及兩斷面構造出一個多面體。由于輪廓線拼接無法由計算機自動、準確地完成，其本質還是由控制點構建輔助線，再由封閉的線成面，最后由一圈封閉的面成體的三維建模。三維建模思路：圖4(a)中暗紅色點為輪廓線上的控制點，通過地質人員的判斷選擇相鄰兩個或者多個斷面上的輪廓線控制點通過點造線的方法構建輔助線(多個斷面之間的輔助線圓滑通過多點插值圓滑實現)；輔助線構建完畢后，圖4(b)、圖4(c)中利用斷面上的輪廓線弧段及構建的兩斷面間的輔助線可以形成一個封閉線圈，由這封閉的線圈即可構建斷面間的輔助面；當一個封閉地質體所需要的輔助地質面構建完成后，選中這個地質體所包含的輔助地質面及斷面本身即可構建一個封閉的地質體(圖4(e)藍色部分為選中的構建地質體所需的部分地質面示意圖)；圖4(f)中藍色部分地質體為構建完成的一個地質體，圖中其他部分地質體也是使用類似方法構建而成的。圖5分別為地層、巖體等屬性體分別進行三維建模的流程示意圖。

圖3　輔助面不同顯示模式示意圖Fig.3　The sketch of auxiliary plane in different display modes(a) 輔助面面模式示意圖；(b) 輔助面線模式示意圖

圖4　單屬性體三維建模流程示意圖Fig.4　The process sketch of single attribute body 3D modeling(a) 構建輔助線；(b) 選中輪廓線及輔助線(c) 構建輔助面；(d) 類似方法構建其他輔助面；(e) 選中構建地質體所需的地質面；(f) 類似方法構建全部地質體

圖5　不同屬性體三維建模流程示意圖Fig.5　The process sketch of different attribute body 3D modeling

1.3不同屬性體切割、融合構建三維地質模型

不同類型屬性體分別完成三維建模后，利用數據轉存和體切割的功能，將不同屬性體的數據模型進行切割合并融合。融合過程就是將不同屬性體模型進行三維布爾求交運算的過程，采用這種方法可以保證融合后模型的完整性，避免各種屬性體同時建模遇到的公共面不重合帶來的諸多問題。圖6(d)是九瑞填圖試驗區(qū)只涉及到地層、巖體、構造屬性體通過切割、融合后構建的三維地質模型，其切割融合方法為構造屬性體作為切割屬性體融合巖體被切割屬性體得到第一次融合模型，再使用這個融合模型作為去融合切割地層屬性體，最后得到整個試驗區(qū)的三維地質模型(圖6(a)、圖6(b)、圖6(c))展示了巖體切割地層的過程)。

圖6　不同屬性體模型相交、融合示意圖Fig.6　The sketch of different attribute body intersect and fusion(a) 被切割的地層屬性體模型；(b) 用于切割的巖體屬性體模型；(c)巖體屬性體切割地層屬性體示意圖；(d) 融合后的三維模型

2　三維地質模型的展示與應用

三維地質模型構建完成后，基于MAPGIS K9的三維地質建模平臺，模型具備了縮放、平移、剝離、疊加、拾取、紋理貼圖、透明顯示以及虛擬漫游、查詢等功能，從而可以對構建的三維地質模型進行全方位的剖析[3]。這里僅介紹以下幾種展示分析方式：

1)點、線、實體模式切換。為方便從不同角度分析模型，系統(tǒng)提供了點模式、線模式和實體模式的互切換功能。

2)交互操作。提供了對三維模型的放大、縮小，適時旋轉、平移，前后移動等窗口操作功能，并可應用鼠標和鍵盤對三維地質模型進行交互操作，實現從任意方位、任意角度觀查三維模型。

3)爆炸顯示。通過爆炸顯示操作，可更加清晰地觀察各種屬性地質模型體的形態(tài)，及分布情況。在爆炸顯示時可選擇任意爆炸、整體爆炸和延軸向爆炸等爆炸方式，還可以設置爆炸顯示的坐標軸、爆炸距離、爆炸幀數和爆炸范圍比例等參數。

4)三維交互定位。利用鼠標操作，將仿真定位點在三維空間的適時移動，可獲取定位點處的三維空間坐標值。

5)組合顯示(圖7)。任意方向切面組合顯示；切面單屬性體、多屬性體組合顯示。

6)三維模型的任意切割技術(圖7)。根據三維模型分析需要，并基于三維模型的幾何和拓撲信息，結合三角網相交處理算法TriCut，實現對三維模型的任意切割，包括：平面斜切、水平剖切，折線垂直剖切，組合切割，以及隧道開挖模擬等多種切割方式。

7)三維模型的局部全屏顯示技術。全屏技術是當前應用相當廣泛的技術。應用該技術，在進行系統(tǒng)演示時，使遠距離的觀看者也能夠清晰地了解系統(tǒng)演示的內容。例如景觀系統(tǒng)在進行漫游操作時，全屏技術可起到了很好的播放效果。

結合九瑞試驗區(qū)三維地質模型，依據該試驗區(qū)成礦規(guī)律，首先尋找熔融有震旦系-寒武系地層，且深部存在規(guī)模較大的燕山期侵入巖體；其次尋找該類侵入巖體與碳酸鹽巖、碎屑巖的接觸帶；再次依據侵入巖體與圍巖的接觸關系及其圍巖的巖性，結合區(qū)內主要控礦因素以及典型控礦層(泥盆系五通組)的分布情況，綜合分析并圈定成礦的有利部位[4-8]。

依據上述條件，依據三維地質模型圈定成礦有利部位共5處，編號為I～V，其大體位置見圖8、9。

I號位于東雷灣以西至鄧家山一線的西南部地區(qū)，該地區(qū)地表出露有受NW向構造控制的小巖珠，深部賦存的侵入巖體為封山洞-東雷灣主巖體及受NW向斷層控制的巖脈；地表出露地層為三疊系灰?guī)r地層，依據沉積層序，深部應存在典型控礦層位即泥盆系五通組地層，立體地質填圖結果顯示泥盆系五通組在該處埋深為1 000 m～1 200 m左右，礦體賦存深度從500 m～1 200 m范圍內，圍巖為三疊系-二疊系灰?guī)r時，可能存在的成礦類型以夕卡巖型為主，灰?guī)r地層層間可能存在塊狀硫化物礦體，深度1 200 m左右與五通組地層接觸部位可能存在層狀塊狀硫化物礦體。

II號位于東雷灣侵入巖體的東側貓母山以西，該處西部有東雷灣巖漿巖體出露，地表出露有三疊系和二疊系灰?guī)r，并存在蝕變，并有NWW向斷層穿越該地區(qū)，結合立體地質填圖結果，從巖漿巖、圍巖、構造以及典型控礦層等因素分析，該處為成礦有利部位。該處典型控礦層(泥盆系五通組)埋深在1 000 m～1 200 m左右，巖漿巖體與三疊系-二疊系灰?guī)r接觸，可能存在的成礦類型以夕卡巖型為主，礦體賦存深度從500 m～1 200 m范圍內；巖漿巖體與五通組地層接觸部位在深度1200 m左右，是有利的礦體賦存部位，可能存在層狀塊狀硫化物礦體。

III號位于封山洞——東雷灣巖體北側，鄧家山東北方向，該處地表出露三疊系-二疊系灰?guī)r，有NE向斷層和NW向斷層交叉存在與該處，立體地質填圖結果顯示該處深部存在隱伏巖漿巖體，且該處地表出露少量受NW向控制的巖脈，其控礦因素的分布規(guī)律符合試驗區(qū)典型礦床的規(guī)律，因此認為該處為成礦有利部位。該處侵入巖體為封山洞-東雷灣主巖體及靠近主巖體的巖脈，典型控礦層(泥盆系五通組)埋深在800 m～1 200 m左右，巖漿巖體與三疊系-二疊系灰?guī)r接觸，可能存在的成礦類型以夕卡巖型為主，礦體賦存深度從500 m～1 200 m范圍內；巖漿巖體與五通組地層接觸部位在深度1 200 m左右，是有利的礦體賦存部位，可能存在層狀塊狀硫化物礦體。

IV號位于武山地區(qū)，涵蓋武山西南部地區(qū)和武山東部地區(qū)，屬武山礦外圍地區(qū)，從立體地質填圖結果看，該處侵入巖體為赤湖主巖體的赤湖-大橋支脈，典型控礦層(泥盆系五通組)在武山西部埋深在800 m～1 100 m左右，在武山南部埋深為1 000 m～1 200 m左右，在武山東部位于復向斜北翼。依據主要控礦因素分析，在武山西南部和南部，由于侵入巖體屬于侵入巖脈的尾部，雖然侵入巖體的圍巖為三疊系-二疊系灰?guī)r，但形成大規(guī)模的夕卡巖型礦體的可能性不大，與五通組地層接觸部位形成層狀塊狀硫化物礦體的可能性較大；在武山東部，由于該處侵入巖體更靠近赤湖主巖體，且五通組地層由于受褶皺影響，主要沿復向斜北翼傾斜分布，故此處礦體賦存深度可從近地表至1 000 m范圍，受地層產狀影響，該處形成規(guī)模較大的夕卡巖型礦體的可能性不大(或規(guī)模不大，主要礦體應為五通組控制的層狀塊狀硫化物礦體)。

圖7　三維地質模型不同方式顯示示意圖Fig.7　The sketch of 3D geological model displayed in different ways(a) 三方向切片組合顯示；(b) 任意切片、屬性體組合顯示；(c) 任意折面、水平切片組合顯示；(d) 模型地表地質圖貼圖；(e) 任意折面、水平切片、單屬性體組合顯示；(f) 單屬性體顯示；(g) 單屬性體水平切片顯示；(h) 模型爆炸顯示；(i) 南北向切片剖面組合顯示；(j) 任意折面切片組合顯示

圖8　五通組地層與巖漿巖體組合俯視圖Fig.8　The combination planform Wutong 　formation and igneous

V號位于太平山東南部至寶山，該處為試驗區(qū)內主要的NW向與NE向構造交匯處，該處地表北部出露志留系地層，南部出露二疊系地層，立體地質填圖結果顯示深部有隱伏巖漿巖體，且該巖漿巖體應為寶山巖體，位于復向斜的北翼且與五通組地層接觸，該部位處于NW和NE向構造結合點范圍內，與區(qū)內成礦規(guī)律完全符合，然而，此處的侵入巖體規(guī)模較小，不利于形成大型礦體。受褶皺構造影響，五通組地層位于復向斜北翼，因此礦體賦存深度自淺部沿層間滑脫斷層分布，深度可達1 000 m，成礦類型可能存在夕卡巖型、斑巖型和塊狀硫化物型。

圖9　九瑞試驗區(qū)推斷深部找礦有利部位分布圖Fig.9　The distribution of inferred favorable parts in Jiu-Rui deep prospecting(a) 500 m深度水平切面與巖漿巖組合；(b) 1 000 m深度水平切面與巖漿巖組合

3　結論與建議

1)利用構建的九瑞試驗區(qū)三維地質模型，從深部地質結構出發(fā)，結合該區(qū)的成礦規(guī)律，依據侵入巖體與圍巖的接觸關系及其圍巖的巖性，結合區(qū)內主要控礦因素以及典型控礦層(泥盆系五通組)的分布情況，綜合分析并圈定了一些成礦的有利部位，可以為該區(qū)深部找礦提供有力的依據。

2)智能化三維地質建模技術是綜合地球物理立體地質填圖工作中的一個非常重要的環(huán)節(jié)，是將地球物理解譯的地質信息轉化為三維地質結構模型的重要表達手段。目前，國內尚沒有適用于復雜地質模型的、成熟完善的三維地質建模技術。這里利用以MAPGIS K9為平臺研發(fā)的、基于二維地質剖面構建三維地質結構模型的單屬性體手工拼接、多屬性體自動融合的三維復雜模型建模技術雖然適合于構建復雜的三維地質模型，但是智能化程度不夠，操作繁瑣且需要專業(yè)的地質人員來進行此項工作。因此，在此基礎上進一步研發(fā)智能化的三維地質建模技術，提高立體地質填圖的工作效率，是以后研究的方向。

[1]杜霍夫斯基.稀有金屬礦區(qū)立體地質填圖[M].夏衛(wèi)華譯.北京:地質出版社,1984.

DUHUO (original work).Rare metals mining area in stereoscopic geological mapping[M].XIA W H (translated).Beijing:Geological Publishing House,1984.(In Chinese)

[2]張生德,張時忠,門吉華.可視化技術及其在地質勘探中的應用淺析[J].地質勘探安全,2000(04):32-36.

ZHANG S D,ZHANG S Z,MEN J H.Visualization technology and its application in geological exploration[J].Geological Exploration and Security,2000(04):32-36.(In Chinese)

[3]楊東來,張永波,王新春,等.地質體三維建模方法與技術指南[M].北京:地質出版社,2007.

YANG D L,ZHANG Y B,WANG X C,et al.3D geologicalmodeling method and technical guideline[M].Beijing:Geological Publishing House,2007.(In Chinese)

[4]翟裕生,姚書振,林新多,等.長江中下游地區(qū)鐵、銅等成礦規(guī)律研究[J].礦床地質,1992,11(1):1-11．

ZHAI Y S,YAO S Z, X D,et al.Metal logenic Regularity of Iron and Copper Deposits in the Middle and Lower Valley of the Yangtze River[J].Mineral Deposits,1992,11(1):1-11.(In Chinese)

[5]王文斌,李文達,董平,等.論長江中下游地區(qū)含銅黃鐵礦型礦床成因[J].火山地質與礦產,1994,15(2):26-33．

WANG W B,LI W D,DONG P,et al.The Genesis of Cuprifer-ous Pyrite Deposits, Middle-Lower Yangtze Area[J]. Volcanology & Mineral Resources[J].1994,15(2):26-33.(In Chinese)

[6]李啟全,羅健雄.九瑞地區(qū)區(qū)域地質構造環(huán)境新述[J].江西地質,1997,5(3):67-73．

LI Q Q,LUO J X.New above of Jiu-rui regional geological ttectonic environment[J].Jiangxi Geology,1997,5(3):67-73.(In Chinese)

[7]董樹文.長江中下游鐵銅礦帶成因之構造分析[J].中國地質科學院院報,1991,6(2):43-56．

DONG S W.Tectonic Analysis on genesis of metallogenetic belt of middle-lower Yangtze river[J].Acta Geoscientica Si-nica, 1991,6(2):43-56. (In Chinese)

[8]王大勇.長江中下游礦集區(qū)綜合地質地球物理研究以九瑞、銅陵礦集區(qū)為例[D].長春:吉林大學,2010.

WANG D Y.The integrated geophysical and geological study in the ore belt of the middle and lower reach of the Yangtze river [D].Changchun:Jilin University,2010.(In Chinese)

Methods and application research of 3D modelling technology in spatial geological mapping

ZHU Wei, WANG Da-yong, WANG Shu-min, YAO Da-wei, WANG Gang

(Institute of Geophysical and Geochemical Exploration, Langfang065000, China)

Firstly, this paper conducts some studies on the construction of 3-D geological model and the method ﹠technology of graphic display based upon integrated geophysical 3-D geological mapping, and puts forward a complex geological structure modeling method ﹠ technology, i.e. different attributes modeling separately, then intersecting and fusing. Secondly, a 3-D geological model of Jiujiang-Ruichang test area is constructed by applying the above complex method on MAPGIS K9 platform, and then comprehensive analysis and display are carried out on this 3-D model through different graphic methods. Finally, by using the 3-D geological model, favorable metallogenetic areas of the test area are selected in terms of deep geological structure.

spatial geological mapping; three-dimensional graph; deep ore prospecting; Jiujiang-Ruichang test area

2015-05-05改回日期：2015-11-04

朱威(1987-)，男，碩士，工程師，主要從事綜合地球物理勘探研究， E-mail：zhuwei@igge.cn。

1001-1749(2016)04-0571-08

TP 302.4

A

10.3969/j.issn.1001-1749.2016.04.21