巖溶體數據建模及屬性查詢研究

周美霞

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

巖溶體數據建模及屬性查詢研究

周美霞

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

地質體的數據建模技術是研究地質體特征的一項重要手段，特別是地質體切片模型的展示，能更加清晰地表達地質體的內部特征、地層分布等情況。主要探討了巖溶體(多Z值)的數據建模，并且在鉆孔數據的基礎上，利用C++語言基于OpenGL編寫程序實現巖溶體的三維可視化，展現了地層的切片模型及地層的屬性查詢。

巖溶體;OpenGL;地層切片;屬性查詢

0 引言

三維地層的構造與分析已經成為地下空間開發和利用的重要手段之一。建立了三維地層模型后，就能夠對某一區域的地下結構有更精確的掌握，能清楚地知道地層的分布狀況，特別是對地層的缺失、倒轉、尖滅等特殊地質現象有很好的了解(趙宏堅等，2010)。一般正常無倒轉的地層都是按照新老地層的順序依次堆積形成的，而巖溶體結構則出現了反常規現象，地層出現了倒轉交替現象，在此地層上鉆探，鉆孔會出現多個高程值的現象。目前，對于一般正常地層的三維可視化建模研究，國內外地質專家和學者都做了充分的解釋，而類似巖溶體這樣有多個Z值(高程值)的地質體的建模情況，國內外還處于探索研究階段。三維地質建模是運用計算機技術，在三維環境下，將空間信息管理、地質解釋、空間分析和預測、地學統計、實體內容分析以及圖形可視化等工具結合起來，并用于地質分析的技術(Houlding，1994)。地層切片技術在地震解釋，石油勘查預測中應用很廣泛，切片模型能夠深入表達地質體的內部形態特征，有益于地學工作者對于地質體內部屬性的把握，在地質體的特征分析中有著重要作用。筆者以二層地層為基礎研究巖溶體的數據建模、地層切片以及地質體屬性查詢。

1 開發環境的構建

實現巖溶體的數據建模，需要搭建開發環境，選擇開發工具，開發環境選擇Visual C++6.0中的MFC可視化程序，結合OpenGL圖形庫繪制三維地層圖形，目前OpenGL是工業界公認的先進而強大的三維圖形編程接口(3D API)。它有著功能完備的二維和三維圖形處理能力，是理想的三維應用程序開發工具(吳愛蘭等，2007)。

2 巖溶體數據建模

2．1 鉆孔數據特點

鉆孔數據排列規則，一般以行為單位規則排列，每一行上的鉆孔數相同，且富含地層特定位置的真實物質，所以能準確地體現地層的巖性、地層分布情況和屬性資料。如X、Y、Z坐標值、孔深、孔口直徑和地層巖性等。圖1顯示的就是某一區域的鉆孔數據。

圖1 鉆孔數據顯示圖

2．2 鉆孔數據的整理與存儲

如何將鉆孔數據整理存儲到計算機中是數據建模的基礎，計算機識別不了實際獲得的鉆孔數據，通過一定的原理方法將原始鉆孔數據規整到計算機能夠識別的形式。在計算機中存儲數據需要考慮到數據結構。線性表是最簡單、常用的一種數據結構，其主要特點是：除了起始節點和終端節點外，其余節點都僅有1個前驅節點和后繼節點，起始節點沒有前驅節點，終端節點沒有后繼節點。線性表有2種存儲方式：順序存儲和鏈式存儲。鏈式存儲是用一組任意的存儲單元存放線性表中的數據元素，簡稱鏈表(史麗燕，2010)。單鏈表是一種基本的線性表。如果每個節點僅包含1個指向其后繼元素的指針，則稱此類節點構成的是單鏈表(史麗燕，2010)。本實例就是將鉆孔數據以單鏈表的形式存儲起來。

2．3 鉆孔插值加密

由于工程鉆探得到鉆孔數據的數量有限，而建模需要足夠的數據才能達到準確性，所以，對于工程鉆探得到的鉆孔數據要進行插值加密處理，設定好鉆孔坐標之間的ΔX、ΔY。結合單鏈表的存儲性質很容易實現源數據的插值加密，存儲加密的鉆孔數據時，要考慮到多Z值的情況。

2．4 巖溶體表面模型的構建

巖溶體的表面是由很多個面元(Surface)集合而成的，每一個面元又是由4個鉆孔點構成，利用插值加密的鉆孔數據形成一系列的面元，巖溶體的表面就形成了。面元的Z值根據鉆孔的Z值判定，一般以左下角的鉆孔Z值代表面元的Z值。用到的函數是OpenGL中的glBegin(GL_QUAD_STRIP)，因巖溶體鉆孔數據多Z值的情況，相應的面元也具有多Z值。

2．5 巖溶體體模型的構建

巖溶體三維可視化建模既是將自然狀態下的巖溶體運用一定的數據模型與方法將其展現在計算機上，那么如何在計算機上將復雜的現實模型表達——任何復雜的數據都可以由簡單的數據集合而成，現實中的地質體模型可以由體元集合展示。在VC++環境中結合OpenGL繪制體元，而體元的繪制要根據體元與對應面元的Z值大小作比較，賦予不同的顏色值，在MFC中結合OpenGL三維可視化展現。在比較的過程中要考慮面元的多Z值情況，分別設定不同的算法，體元的繪制只需通過三重“for循環”，分別在X、Y、Z 3個方向上循環1遍。體元的大小要和面元一致。遇到面元有2個或者3個Z值的設定算法確保體元繪制正確。算法的正確性直接決定體模型建立的正確與否。

2．5．1 遇到面元有2個Z值情況 依據研究區域的巖溶體特點判斷算法的設定，主要是確定面元是在研究區域邊界還是中間位置。(1)第一種情況：面元處于研究區域邊界，當體元Z值小于面元的最小Z值時，設置體元有顏色;當體元Z值大于面元的最小Z值時，設置體元呈透明顯示(或者另一種不同顏色，只要統一即可)。(2)第二種情況：位于中間，當體元Z值小于面元的最大Z值時，設置體元有顏色，當體元Z值大于面元的最大Z值時，設置體元透明顯示。

2．5．2 遇到面元有3個Z值情況 當體元Z值小于面元的最小Z值時，設置體元有顏色;當體元Z值大于面元的最小Z值小于中間Z值時，設置體元透明顯示;當體元Z值大于面元的中間Z值小于最大Z值時，設置體元有顏色;當體元Z值大于面元的最大Z值時，設置體元透明顯示。

研究實例如圖2，為巖溶體的建模效果圖。

圖2 巖溶體模型效果圖

3 巖溶體切片模型

在體模型的基礎上，進一步展現切片模型，對于地質體的研究有更深的意義，地震切片分析技術是地震資料解釋的一種常用手段，其在油氣勘探中的作用愈來愈大(張軍華等，2007;王鵬等，2008)。通過切片可以更清晰地看到地層內部的巖性，為進一步研究地質體的物理、化學屬性提供了依據。切片模型的展示既是在體模型的基礎上設置一定的算法使特定位置的地質體展現，而其他位置的地質體不顯示，那么整個模型就呈現出切片的形式。定義一個實現切片模型的函數void Sectionmodel();算法實現如下：

切片模型實現如圖3、圖4。

圖3 切片模型1

圖4 切片模型2

4 巖溶體屬性查詢的實現

本實例在編程實現地質體數據建模的同時，利用OpenGL里的glLoadName(name)函數，分別給每一層地層命名，然后再為不同的地層附上相應的屬性資料，在可視化的過程中實現點擊查詢地層的屬性。glLoadName(name)既是在渲染模型的同時賦予一個名字，當點擊到這個名字的地層模型時調用相應的函數實現功能。功能實現界面如圖5。

圖5 屬性查詢界面圖

5 結論

三維地學模擬是三維地質信息GIS可視化的核心研究內容。隨著科學計算可視化技術和三維GIS技術的發展，地質信息的三維可視化也將得到進一步的深化(朱良峰等，2006)。如今，地質問題已成為貫穿于數字城市、減災防災、城建工程、地下工程、水電工程、交通工程、環境工程、資源開發的一個基礎性問題，三維地學模擬是解決這些領域的地質問題的關鍵(周翠英等，2005)。隨著計算機技術的不斷發展，地質研究者們也在努力應用各種方法將在實地考察中得到的地質資料和地層信息移位到計算機屏幕上，使勘查得到的地質資料可視化表達及資源共享。通過不斷的努力和研究，目前國內外也取得了一定的成就，但是學術研究是無止境的，需要更多的人為此付諸努力和汗水。地層切片模型既是對地質體模型的進一步分析，通過切片又可以了解到地層內部變化情況，為各種地質現象的判斷提供了依據。

史麗燕．2010．單鏈表基本操作的實現［J］．軟件導刊，(2)：21－22．

吳愛蘭，樓建列．2007．基于OpenGL的三維圖形的繪制［J］．中國科技信息，(9)：122 －124．

王鵬，高偉，張紅斌．2008．地震切片演化技術在乍得A區塊的應用［J］．石油地球物理勘探，43(增刊1)：115－118．

周翠英，陳恒，劉祚秋，等．2005．重大工程地下環境信息系統的特點［J］．巖土力學，26(2)：216 －220．

朱良峰，潘信，吳信才．2006．三維地質建模及可視化系統的設計與開發［J］．巖土力學，27(5)：828 －832．

張軍華，周振曉，譚明友，等．2007．地震切片解釋中的幾個理論問題［J］．石油地球物理勘探，42(3)：348－352．

趙宏堅，周翠英．2010．基于實體建模的三維地層構造［J］．巖土力學，31(4)：1257 －1263．

HOULDING S W．1994．3D Geoscience Modeling：Computer Techniques for Geological Characterization［M］．New York：Springer-Verlag．

Study on karst data modeling and attribute inquiry

ZHOU Mei-xia

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

Geological body data modeling technology was an important way to study the properties of the geological body，to display geological body section model，to express the inner properties and stratigraphic distribution of the geological body．The author discussed the data modeling of karst body(multi Z value)，realized a 3D visualization of karst body based on OpenGL programming and by the use of C++language on the foundation of bore data．In addition，the stratigraphic section display and attribute inquiry were undertaken besides 3D visualization．

Karst;OpenGL;Stratigraphic section;Attribute inquiry

TP391

A

1674－3636(2012)03－0316－04

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

2012－05－20;編輯：陸李萍

周美霞(1988— )，女，碩士研究生，地圖制圖學與地理信息工程專業，主要從事地質體的數據建模與三維可視化研究，E-mail：zmx19880222@126．com