劉 星，劉瑩瑩

(安徽理工大學(xué)，安徽 淮南 232001)

基于真三維GIS的地質(zhì)構(gòu)造模型與空間分析

劉 星，劉瑩瑩

(安徽理工大學(xué)，安徽 淮南 232001)

帶有斷層的地質(zhì)體建模難于層狀地層，是當(dāng)前研究的難點(diǎn)，模型完成后的三維空間分析功能目前也處于研究階段，是當(dāng)前研究的熱點(diǎn)。以某石膏礦為例，介紹了利用存在地層倒轉(zhuǎn)、地層缺失等復(fù)雜鉆孔數(shù)據(jù)進(jìn)行三維建模的流程以及模型建立后的三維空間分析方法和地學(xué)應(yīng)用，嘗試了三維GIS空間－屬性查詢、空間裁剪、縱剖面橫切面的自動(dòng)生成方法和實(shí)例研究。結(jié)果表明，地質(zhì)體和構(gòu)造的特征的三維建模和空間分析有利于地質(zhì)演化的理解和開采設(shè)計(jì)的優(yōu)化。

真三維;GIS;空間分析;地學(xué)應(yīng)用

0 引言

一個(gè)典型的三維地質(zhì)模型是由1個(gè)或多個(gè)地層、結(jié)構(gòu)面、斷層等組成的(趙曉東等，2009;明鏡等，2009)。對(duì)于平原區(qū)的地質(zhì)分析模型主要涉及第四系土層、地下水等。地層之間的關(guān)系一般比較簡(jiǎn)單，地層相互平行。在構(gòu)造復(fù)雜的情況下，存在地層相交、錯(cuò)斷、尖滅、透鏡體、夾層等地質(zhì)現(xiàn)象，使三維地質(zhì)建模過程復(fù)雜化(趙曉東等，2009;李燦輝等，2010)。

三維GIS作為一種建模工具，能夠大大提高地質(zhì)建模的效率和精度，可以滿足對(duì)復(fù)雜地質(zhì)區(qū)域的建模要求。當(dāng)前，具備三維空間分析的軟件最有代表性的是法國南希大學(xué)開發(fā)的GOCAD(李燦輝等，2010;劉光偉等，2010)，筆者試圖研究利用該軟件建立三維地質(zhì)模型一般流程：數(shù)據(jù)的分析和預(yù)處理、插值、建立三維地層界面、建立三維地層實(shí)體、剖面圖對(duì)比等步驟。

1 三維地層數(shù)據(jù)模型

目前，構(gòu)建三維地質(zhì)模型的方法可分為兩大類：基于面模型的構(gòu)建方法和基于體模型的構(gòu)建方法(向中林等，2009)。前者側(cè)重于用地質(zhì)實(shí)體分界面(如地層面、接觸面和邊界面)來形成地質(zhì)實(shí)體的空間輪廓;后者則側(cè)重于三維空間中地質(zhì)實(shí)體邊界與內(nèi)部的整體表示，通過對(duì)三維空間進(jìn)行體元分割來實(shí)現(xiàn)地質(zhì)實(shí)體的真三維表達(dá)和分析。在實(shí)際建模過程中，經(jīng)常將這兩種方法結(jié)合起來，例如，可先利用基于面模型的構(gòu)建方法形成地質(zhì)實(shí)體的分界面，然后再對(duì)其進(jìn)行體分割。因此，無論是面模型還是體模型，都需要先構(gòu)建出地質(zhì)體的分界面。

地質(zhì)實(shí)體和地質(zhì)界面一起組成三維地質(zhì)模型，它的建立過程直接影響模型的精確度。在GOCAD中，若直接利用Solid功能生成地質(zhì)實(shí)體(Solid)，是非常困難且無實(shí)際意義的。而使用GOCAD Solid功能，也可以生成滿足要求的地質(zhì)實(shí)體，該法與Solid相比有許多優(yōu)勢(shì)，主要表現(xiàn)在以下3個(gè)方面。

(1)采用SGrid(柵格)形成實(shí)體，它的實(shí)現(xiàn)方法簡(jiǎn)單易操作，利用先前生成的層面即可完成。

(2)插值過程中得到的相關(guān)數(shù)據(jù)是用于工程計(jì)算的寶貴資料。進(jìn)行地質(zhì)建模的目的，將地質(zhì)結(jié)構(gòu)更鮮明直觀地展現(xiàn)出來，是將插值過程中得到的相關(guān)數(shù)據(jù)用于進(jìn)一步的地質(zhì)評(píng)價(jià)，SGrid則將這二者完美地結(jié)合起來。

(3)柵格的節(jié)點(diǎn)信息和拓?fù)潢P(guān)系可以Excel表格的形式導(dǎo)出，SGrid導(dǎo)出的節(jié)點(diǎn)信息和拓?fù)潢P(guān)系可以通過簡(jiǎn)單的程序轉(zhuǎn)換，進(jìn)而導(dǎo)入到 ANSYS和ABAQUS等有限元程序中，實(shí)現(xiàn)交互式計(jì)算，拓展了軟件功能。

2 鉆孔數(shù)據(jù)分層與可視化

本次研究的主要數(shù)據(jù)來源于鉆孔柱狀圖，共24個(gè)鉆孔，其中的層面和標(biāo)志層劃分不夠詳細(xì)。因此，首要任務(wù)是將每個(gè)鉆孔按照建模的要求，重新厘定劃分，重新劃分的地層序列為松散層(Q+N)、白堊系(K)、銅頭尖組(T2t)和 月山組(T2y)，上面兩組為砂泥巖，東馬鞍山組第三段(T2d3)，為巖溶灰?guī)r，東馬鞍山組第二段(T2d2)，為石膏層位，東馬鞍山組第一段(T2d1)，為灰?guī)r或硅質(zhì)灰?guī)r，最下部為南陵湖組(T1n1)灰?guī)r。

根據(jù)上述分層原則，將每個(gè)鉆孔的分層界面轉(zhuǎn)換為真深度，建立各個(gè)鉆孔的分層數(shù)據(jù)庫。首先根據(jù)鉆孔開孔坐標(biāo)生成井位分布圖，然后將每個(gè)鉆孔分層數(shù)據(jù)輸入對(duì)應(yīng)的鉆孔，形成鉆孔模擬圖。在此基礎(chǔ)上，將每個(gè)鉆孔的相鄰分層組合成為區(qū)，模擬不同的地層。將存在地層次序倒轉(zhuǎn)(黃色)、缺失(綠色)的鉆孔用不同標(biāo)志表明，有利于從總體上把握構(gòu)造的分布(圖1)。

圖1 地層缺失和地層倒轉(zhuǎn)的鉆孔位置

3 非連續(xù)層面地層的地質(zhì)模型

以上方法在針對(duì)層狀連續(xù)地層的模型建立中十分方便，但在構(gòu)造發(fā)育導(dǎo)致層面錯(cuò)斷時(shí)很難應(yīng)用，究其原由，則是基于地表的數(shù)字高程模型不能適應(yīng)地下地層、構(gòu)造面等非連續(xù)、重復(fù)性地層(陳冰凌等，2009;Apel，2005)。在該石膏礦研究區(qū)，要建立地質(zhì)三維模型，必須首先劃分塊段，確定每個(gè)層面的斷塊邊界，針對(duì)單個(gè)連續(xù)的斷塊建立模型，然后在空間上進(jìn)行耦合。

該石膏礦的地層和構(gòu)造在空間上的分布特征還不完全清楚，尤其是東部地區(qū)，鉆孔數(shù)量偏少，難以控制其地層變化規(guī)律。但通過已有鉆孔揭示的信息來建模，在三維空間觀察其變化特征，比在二維狀態(tài)下分析地質(zhì)構(gòu)造的變化要直觀、有效。

研究將石膏礦劃分為4個(gè)斷塊，這在三維模型中非常明顯。當(dāng)然可能還存在小的斷裂構(gòu)造，但不影響地質(zhì)建模與分析(圖2)。

圖2 某石膏礦的塊段初步劃分

將地層建模的層序劃分為第四系、白堊系，以上兩套地層在本區(qū)域具有統(tǒng)一的界面，既可以通過鉆孔中的層面分界直接插值形成一個(gè)連續(xù)統(tǒng)一的界面。下面的地層分層比較復(fù)雜，西南角的斷塊1，白堊系之下的地層系列為：Tt3、Ty3、Td3、Td2、Td1，屬于正常層序。而斷塊2的層序序列自上而下分別為：Tt3、Ty3、Td3、Td2、Td1，也屬于正常序列，但不完整，不是全局分布的。斷塊3的地層自白堊系之下分別為Td3、Td2、Td1，在與斷塊2接觸之處形成逆斷層，部分鉆孔地層倒轉(zhuǎn)。斷塊4的地層與斷塊3類似。在4個(gè)塊體之外，分別為左上角和右下角，存在橫穿全區(qū)的大型斷裂，左上角的斷塊是一個(gè)逆斷層上盤(上升盤)，出露老的T1n1(南陵湖組)和T3d1(東馬鞍山組下段)，石膏層被剝蝕殆盡，該處的鉆孔出現(xiàn)地層順序倒轉(zhuǎn)現(xiàn)象，即石膏層之上出現(xiàn)南陵湖組老地層，說明鉆孔位于斷層之上。圖中黃色標(biāo)志說明是有地層次序倒轉(zhuǎn)的鉆孔。右下角為一個(gè)大型正斷層，在該斷層附近，出現(xiàn)地層缺失現(xiàn)象，即在白堊系之下，直接與南陵湖組地層接觸，表明該斷層下盤為上升盤，顏色為綠色的鉆孔即為地層缺失的鉆孔。按照該建模思路，需要分別界定各個(gè)斷塊每個(gè)地層在空間上的邊界，保證界線在空間上分布合理，否則建立的模型在空間上會(huì)存在交叉、分離等不合理現(xiàn)象。首先是確定斷裂切穿的層位，在本研究區(qū)，斷裂切穿了自白堊系以下的所有地層。這樣將白堊系底作為模型的最上界面，在斷塊1下，按照計(jì)算的基巖地層邊界，依次產(chǎn)生 T2t、T2y、T2d3、T2d2、T2d1、T1n 界面。

對(duì)于斷塊2，按照地質(zhì)資料分析，基巖依然存在4套地層，但本塊段只有1個(gè)鉆孔控制，無法利用數(shù)據(jù)插值方法形成層面，需要采用地質(zhì)推斷方法產(chǎn)生各個(gè)地層的層面，其依據(jù)是根據(jù)單個(gè)鉆孔將礦體頂板等高線上移一定高度，產(chǎn)生新的地層界面。本塊段產(chǎn)生的地質(zhì)層面自上而下依次為T2t、T2y、T2d3、T2d2(石膏層)、T2d1，其層面邊界的確定如圖3。

圖3 塊段2地層邊界確定和層面建立

利用上述插值方法或推算的方法得到的層面，通常會(huì)產(chǎn)生不合理的現(xiàn)象，例如由于地層的尖滅，第四系和白堊系地層在東北部產(chǎn)生重合現(xiàn)象，但插值會(huì)造成微小的K底界位于Q底界之上，這樣的層面在后面無法構(gòu)成空間體模型，需要采用腳本語言強(qiáng)令K底界位于Q底界之下，但由于兩個(gè)層面在垂向上不完全耦合，腳本語言不能夠?qū)⒂行У膬蓚€(gè)層面Z值對(duì)比，必須借助于其他方式加以解決。

同樣，塊段2的K底面有一部分位于銅頭尖組之上，這是不合理的，可以強(qiáng)令銅頭尖組地層Z值小于等于K底面的Z值，這樣銅頭尖組地層自西南向東北尖滅于白堊系之下，這與地質(zhì)演化和后期保存狀態(tài)是相符合的。

斷塊3范圍比較大，基巖地層埋深較淺，并且東北高西南低，自白堊系之下地層為東馬鞍山組第三段，東馬鞍山組第二段和東馬鞍山組第一段，缺失銅頭尖組(T2t)和月山組(T2y)，地層次序正常。

斷塊4的地層分布次序與斷塊3相似，基巖地層埋深也較淺，自白堊系之下地層為東馬鞍山組第三段、東馬鞍山組第二段、東馬鞍山組第一段，缺失銅頭尖組(T2t)和月山組(T2y)，地層次序也正常。

4 礦區(qū)真三維地層模型構(gòu)建

利用塊段劃分建立的各個(gè)地質(zhì)界面，建立該研究區(qū)的三維實(shí)體模型，三維實(shí)體模型的建立，有利于后續(xù)的地質(zhì)切面分析，地質(zhì)屬性分析，夾層分布，巖溶水分布可視化和礦體綜合評(píng)價(jià)，儲(chǔ)量估計(jì)(朱發(fā)華等，2009;Andrea et al，2009)等。具體步驟如下。

(1)根據(jù)鉆孔分層資料，提取分層的高程數(shù)據(jù)，導(dǎo)入到GOCAD中，生成原始點(diǎn)分布圖。

(2)進(jìn)行克里金插值，導(dǎo)入到GOCAD中后，生成加密后的數(shù)據(jù)點(diǎn)分布圖。

(3)以各個(gè)鉆孔分層點(diǎn)作為控制點(diǎn)，在地層尖滅、缺失處獲得交線。

(4)形成Delaunay三角網(wǎng)格化的層面。利用GOCAD中里的DSI插值算法優(yōu)化各層初始層面，用交線約束修改原尖滅、缺失地層層面，從而最終建立各地層的層面。

(5)通過上下面，生成體模型。由于上下面形狀不對(duì)稱，必須使用Workflow方法，控制上下層面的連接部位。

同一斷塊自上而下存在多個(gè)地層體，依次建立各個(gè)地層體，然后融合(圖4)。

圖4 某石膏礦三維立體模型

(6)體模型的屬性傳遞。將鉆孔中石膏層的層位，夾層厚度，作為點(diǎn)元素輸入，將各個(gè)點(diǎn)的厚度值作為屬性，傳遞給體模型，利用空間三維插值方法，將夾層數(shù)目展現(xiàn)在切片中。

(7)三維切片。將每個(gè)地層自新到老，依次賦予一定的數(shù)值，例如Q賦值為0，K賦值為1，T2t賦值為2，T2y賦值為 3，T2d3賦值為 4，T2d2賦值為 5，T2d1賦值為6。切片可以在3個(gè)方向進(jìn)行，可以顯示地層間的接觸關(guān)系、空間分布狀況和構(gòu)造的控制情況(圖 5)。

圖5 某石膏礦三維任意切片圖

通過構(gòu)造單個(gè)的任意面，生成斷面，輸入到CAD軟件進(jìn)行編輯整飾(圖6)。

圖6 通過層面推斷構(gòu)造

5 三維圖形地質(zhì)分析

5．1 地層的三維空間分布特征

由于連續(xù)層序的界面在沉積上是連續(xù)整合的，除非后期受構(gòu)造抬升發(fā)生剝蝕，其層面才會(huì)造成起伏，排除這種情況，如果沒有遭受剝蝕的層面，也會(huì)發(fā)生起伏變化，究其原因有兩種：一是沉積基底的不平，二是構(gòu)造運(yùn)動(dòng)造成的斷層。對(duì)于前者，其起伏應(yīng)當(dāng)是光滑的曲面，而斷層造成的層面起伏，則在曲率上比較大。根據(jù)該原則，利用三維空間的地質(zhì)體變化特征，利用人工解釋或者二次導(dǎo)數(shù)的方法，可以提取出未經(jīng)鉆孔驗(yàn)證的斷裂分布的特征(黃瑞婕等，2010)，供勘探布網(wǎng)參考和為采礦工程提供依據(jù)。

圖7是礦體底面，其沉積時(shí)期應(yīng)該適合下覆底層整合接觸，層面與礦體形態(tài)的不協(xié)調(diào)應(yīng)當(dāng)是后期構(gòu)造運(yùn)動(dòng)造成，現(xiàn)用紅色標(biāo)注于上圖，將來的鉆孔應(yīng)當(dāng)優(yōu)先考慮紅色區(qū)域。

圖7的地質(zhì)體難以反映層內(nèi)構(gòu)造形跡，但清楚地表明兩斷塊抬升速度有差異，上部斷塊由于出露地表剝蝕殆盡，下面的斷塊是西部抬升的較高。它們的上層面基本表明是后期的侵蝕面。

圖7 地質(zhì)體的不均衡構(gòu)造抬升

5．2 地質(zhì)構(gòu)造空間分布與解釋

通過構(gòu)造三維模型，能夠發(fā)現(xiàn)先前工作中構(gòu)造解釋的不合理性，為下一步布置勘探鉆孔提供依據(jù)。圖8為石膏礦底板三維圖，在鉆孔4—4處，底板連續(xù)性與其他各孔不協(xié)調(diào)，應(yīng)該存在1條正斷層，才造成礦體減薄。下一步鉆孔位置應(yīng)在此附近布置，以便控制礦體變化，探明構(gòu)造變化特征。

圖8 通過三維層面推斷構(gòu)造特征

通過前后兩期構(gòu)造解釋情況下的模型對(duì)比，可以發(fā)現(xiàn)在實(shí)施物探措施之后，斷裂解譯的合理性比前期解譯合理。圖9是解譯前的構(gòu)造特征，中間的斷塊起伏太大，礦體扭曲得非常嚴(yán)重，難以真實(shí)存在，而右圖則是后期解譯的構(gòu)造，中間的斷塊過渡自然，礦體平穩(wěn)過渡。

圖9 通過三維模型推斷前期構(gòu)造的合理性

5．3 物性及巖溶水空間分布與解釋

通過將巖溶水分別特征與三維模型結(jié)合，進(jìn)行切片觀看，能夠預(yù)先了解巖溶水的空間立體分布。將石膏層的灰?guī)r夾層作為立體模型的屬性，能夠利用切片的方式觀察石膏層中的夾層分布情況，有利于總體把握礦體豐度分布(圖10、圖11)。

圖10 三維切片的地質(zhì)屬性觀察

圖11 某石膏礦－300 m切面含水分布、石膏中夾層分布圖

5．4 礦石量的計(jì)算

通過三維模型的建立，能夠精確計(jì)算該研究區(qū)的石膏資源量，本次計(jì)算大約為0．84億m3，該資源量是沒有扣除灰?guī)r夾層、含水層和采礦設(shè)計(jì)的總體資源量。

6 結(jié)論

介紹了非連續(xù)斷塊地層的三維模型建立方法，利用某石膏礦的鉆孔資料，從分層厘定、數(shù)據(jù)庫建立到三維模型的展現(xiàn)以及地質(zhì)屬性空間分析，詳細(xì)敘述了建模的整個(gè)過程。

通過模型的建立和屬性分析，不僅可以方便數(shù)據(jù)管理，便于地下數(shù)據(jù)可視化，還可以從真實(shí)角度分析構(gòu)造特征，重新解譯斷裂邊界，掌握礦體真實(shí)賦存狀態(tài)，為總體分析把握礦體豐度分布提供可信的數(shù)據(jù)模型。

通過切片方式，很容易計(jì)算和分析石膏礦中夾層的分布和含水層位的空間分布，為地質(zhì)勘探和采礦工程提供依據(jù)。

陳冰凌，王曉鵬．2009．真三維地質(zhì)體建模技術(shù)及其在煤田地質(zhì)勘探的應(yīng)用［J］．中國煤炭地質(zhì)，(增刊2)：123 －126．

黃瑞婕，劉紅欣．2010．基于三維地質(zhì)建模的斷裂構(gòu)造識(shí)別技術(shù)［J］．石油天然氣學(xué)報(bào)，32(4)：217－220．

李燦輝，潘勇，吳強(qiáng)．2010．復(fù)雜礦床三維建模與地質(zhì)自動(dòng)繪圖技術(shù)研究［J］．計(jì)算機(jī)應(yīng)用研究，27(10)：3961 －3963．

劉光偉，白潤才，呂進(jìn)國，等．2010．基于三維地質(zhì)實(shí)體模型生成地質(zhì)剖面圖的應(yīng)用［J］．遼寧工程技術(shù)大學(xué)學(xué)報(bào)，29(4)：557 －559．

明鏡，潘懋，屈紅剛，等．2009．北京市新生界三維地質(zhì)結(jié)構(gòu)模型構(gòu)建［J］．北京大學(xué)學(xué)報(bào)：自然科學(xué)版，(1)：111－119．

向中林，白萬備，王妍，等．2009．基于Surpac的礦山三維地質(zhì)建模及可視化過程研究［J］．河南理工大學(xué)學(xué)報(bào)：自然科學(xué)版，28(3)：307 －311，320．

趙曉東，李利崗，彭林軍．2009．GIS和GOCAD支持下的礦山3D 地質(zhì)建模［J］．地理與地理信息科學(xué)，25(2)：34－38．

朱發(fā)華，賀懷建．2009．基于地質(zhì)雷達(dá)和鉆孔數(shù)據(jù)的三維地層建模［J］．巖土力學(xué)，30(增刊1)：267 －270．

ARNAUD BERLIOUX．1994．Building models with GOCAD［J］．Stanford Exploration Project，(80)：588 －607．

APEL M．2005．A 3D geological information system framework［J］．Geophysical Research Abstracts，7：05732．

ANDREA ZANCHI， SALVIFRANCESCAC， ZANCHETTA STEFANO，et al．2009．3D reconstruction of complex geological bodies：Examples from the Alps［J］．Computers ＆Geosciences，35(1)：49 －69．

OLIVIER KAUFMANN，THIERRY MARTI．2008．3D geological modeling from boreholes，cross-sections and geological maps，application over former natural gas storages in coal mines［J］．Computers ＆ Geosciences，34(3)：278 －290．

Geological structure model and spatial analysis based on true 3D GIS

LIU Xing，LIU Ying-ying

(Anhui University of Science and Technology，Huainan 232001，Anhui)

The geologic modeling for geologic bodies with faults was a difficulty，which was more complicated than those layered strata and its spatial analysis function after the achievement of the model was also a new research arena．Taking a gypsum mine as an example，the authors introduced a flowsheet of building a 3D model by complicated borehole data with strata loss and strata reverse，3D spatial analysis methods on those models and their geologic application．It was concluded that 3D modeling and its spatial analysis for geologic bodies and structure properties were advantageous to the understanding of geologic evolution and optimization of mining design．

True 3D;GIS;Spatial analysis;Geological application

TP391

A

1674－3636(2012)03－0274－06

10．3969/j．issn．1674-3636．2012．03．274

2012－05－20;編輯：侯鵬飛

國家自然科學(xué)基金項(xiàng)目(40872103)資助

劉星(1974— )，男，博士研究生，主要從事成礦預(yù)測(cè)研究，E-mail：liuxing0795@126．com