姬廣軍 朱吉祥

摘要：三維地質(zhì)建模是定量化研究三維空間結(jié)構(gòu)特征與變化規(guī)律的重要方法，已成為模擬分析地質(zhì)現(xiàn)象與時(shí)空演變規(guī)律重要有效手段之一。本文從數(shù)據(jù)模型、建模方法、建模軟件、不確定性及應(yīng)用領(lǐng)域五個(gè)方面闡述了三維地質(zhì)建模技術(shù)的研究現(xiàn)狀，并對(duì)三維地質(zhì)建模技術(shù)發(fā)展的提出相應(yīng)建議。

關(guān)鍵詞：三維地質(zhì)建模;研究現(xiàn)狀;不確定性;數(shù)據(jù)模型;建模方法

三維地質(zhì)建模是在各種原始數(shù)據(jù)的基礎(chǔ)上，以適當(dāng)?shù)臄?shù)據(jù)模型建立地質(zhì)特征的數(shù)學(xué)模型，通過對(duì)地質(zhì)體的幾何形態(tài)、地質(zhì)體間的關(guān)系和物性等計(jì)算機(jī)模擬，最后形成復(fù)雜整體三維模型的過程。[1]三維地質(zhì)建?？梢灾庇^地表達(dá)地質(zhì)體的空間形態(tài)與特征，能夠進(jìn)行真三維的空間操作與分析，較大程度地促進(jìn)地質(zhì)信息的科學(xué)管理與共享。本文從空間數(shù)據(jù)模型、建模方法、建模軟件、不確定性和應(yīng)用領(lǐng)域五個(gè)方面入手，系統(tǒng)闡述了三維地質(zhì)建模的研究現(xiàn)狀，并對(duì)三維地質(zhì)建模發(fā)展提出相應(yīng)建議。

1 空間數(shù)據(jù)模型

數(shù)據(jù)模型是為方便數(shù)據(jù)的存儲(chǔ)、管理、交換和使用而構(gòu)建的數(shù)據(jù)結(jié)構(gòu)、數(shù)據(jù)操作和數(shù)據(jù)約束的總稱。1978年Hunter首先提出了以八叉樹為代表的數(shù)據(jù)模型，[2]Victor[3]和Pilouk[4]提出了不規(guī)則四面體的數(shù)據(jù)模型;Wu提出了廣義三棱柱的數(shù)據(jù)模型[5]等。這些數(shù)據(jù)模型可以分為三類：基于面、基于體以及混合數(shù)據(jù)模型。這三類模型中常用的模型的優(yōu)缺點(diǎn)見表1。

3 三維地質(zhì)建模軟件

三維地質(zhì)建模技術(shù)的不斷成熟和計(jì)算機(jī)技術(shù)的快速發(fā)展，促使國(guó)外開發(fā)研制了許多較為成熟的商業(yè)化軟件。其中應(yīng)用較廣泛的軟件有GoCAD、Surpac、3DMine、Vulcan、MicroLynx、GSI3D等。GoCAD軟件的三維可視化功能強(qiáng)大，能處理復(fù)雜的地質(zhì)構(gòu)造，是目前國(guó)外應(yīng)用最廣的三維地質(zhì)建模軟件之一。

國(guó)內(nèi)在軟件方面起步較晚，且大都集中在三維可視化等方面，[10]主要包括MapGIS、GSIS建模系統(tǒng)及QuantyView等。國(guó)內(nèi)軟件具有較強(qiáng)的專業(yè)針對(duì)性，如Longruan GIS面向數(shù)字礦山、DeepInsight面向地球物理與油藏等，它們基本可以滿足大多數(shù)用戶的日常需求。[22]

4 不確定性分析

由于目前人類知識(shí)的有限性及地下空間的不可見性與復(fù)雜性，通過有限的信息來(lái)推斷未知的地下信息會(huì)存在各種形式的不確定性。[23]國(guó)內(nèi)外對(duì)模型的不確定性做出了不少研究，但仍處在初始階段。[24]三維地質(zhì)建模的不確定性包括數(shù)據(jù)不確定、[24]結(jié)構(gòu)不確定、[25]可視化不確定[26]和響應(yīng)不確定，[26]其中數(shù)據(jù)不確定和結(jié)構(gòu)不確定性研究較多，這里主要討論數(shù)據(jù)不確定性和結(jié)構(gòu)不確定性。

4.1 結(jié)構(gòu)不確定性

地質(zhì)學(xué)家可以利用地球物理勘探、露頭描述及其延伸、鉆孔揭露等方式獲取有限范圍的地質(zhì)信息，這些方法通過專家解釋后獲得的地質(zhì)結(jié)構(gòu)可能并不能真正的反映實(shí)際的地質(zhì)結(jié)構(gòu)，從而產(chǎn)生一定程度的不確定性。另一方面，地質(zhì)體或地質(zhì)現(xiàn)象本身就具有復(fù)雜性，[27]通過各種常規(guī)方法建立的模型也會(huì)存在一定程度的不確定性。Thore等[25]提出了一種確定和處理結(jié)構(gòu)不確定性的方法并詳細(xì)地分析不確定的來(lái)源;Gerritsen和Caers[28]指出了不確定的來(lái)源并且認(rèn)為數(shù)據(jù)可以用來(lái)約束模型的不確定性，此外，他還給出了基于概率、變異函數(shù)、逆向建模等方法的詳細(xì)解釋和步驟。

4.2 數(shù)據(jù)不確定性

幾乎所有的地質(zhì)數(shù)據(jù)都會(huì)有不同程度的不確定性，這些不確定性是由測(cè)量誤差、數(shù)據(jù)不足、取樣限制、不完美的概念或假說和必要的簡(jiǎn)化等因素引起的，[28]這也是限制模型精度的重要原因。Wellmann等[29]分析了由于數(shù)據(jù)質(zhì)量而引起的不確定性并給出了具有詳細(xì)步驟的處理辦法，他認(rèn)為數(shù)據(jù)的不確定來(lái)源可以分為3類：

（1）精度不夠或測(cè)量誤差，是由于原始數(shù)據(jù)的不確定性引起的，如沉積邊界的位置或構(gòu)造的方向。

（2）地質(zhì)變量本身的隨機(jī)性，是由已知點(diǎn)插值得到未知點(diǎn)不確定性引起的。

（3）認(rèn)知有限性，是由于對(duì)存在的地質(zhì)體或地質(zhì)現(xiàn)象認(rèn)知不完整和不精確引起的。[30]

由采樣點(diǎn)的數(shù)據(jù)所建立的三維地質(zhì)模型一定存在不確定性，一般不確定性隨著數(shù)據(jù)量的減小而增大。Yamamoto等[30]通過在計(jì)算模型中再取樣和后處理的方式達(dá)到減小不確定性的效果。

5 應(yīng)用領(lǐng)域

經(jīng)過幾十年的發(fā)展，三維地質(zhì)建模已經(jīng)滲入到地質(zhì)的各個(gè)領(lǐng)域，包括礦產(chǎn)、[31]石油、[32]水文地質(zhì)、[33，34]工程地質(zhì)、[7，35]地?zé)岬刭|(zhì)、[36]災(zāi)害地質(zhì)等。三維地質(zhì)建模首先使用到礦產(chǎn)、石油領(lǐng)域，[23]且在這些領(lǐng)域中技術(shù)比較成熟。在礦產(chǎn)中，主要用于分析礦體的分布規(guī)律、估算礦產(chǎn)資源量以及進(jìn)行資源量評(píng)價(jià)等服務(wù);在石油領(lǐng)域，主要進(jìn)行儲(chǔ)量計(jì)算、表達(dá)儲(chǔ)層構(gòu)造形態(tài)與品質(zhì)分布等。在水文地質(zhì)學(xué)科中，三維地質(zhì)建?？梢杂糜诜治鲅芯繀^(qū)的水文地質(zhì)條件、調(diào)查污染物運(yùn)移規(guī)律等服務(wù)。在工程地質(zhì)中，三維地質(zhì)建模用于對(duì)巷道、工程建設(shè)的承載力模擬分析等服務(wù)。三維地質(zhì)建模還可以估算地?zé)醿?chǔ)量、分析地質(zhì)災(zāi)害致災(zāi)機(jī)制等服務(wù)。

6 結(jié)論

三維地質(zhì)建模突破了傳統(tǒng)二維空間表達(dá)的束縛，還地質(zhì)以本來(lái)面目，其可以直觀、形象地展示了地質(zhì)體的空間分布及其規(guī)律，使地質(zhì)研究可以在多維空間中分析地質(zhì)問題，表達(dá)研究結(jié)果，更符合地質(zhì)研究的實(shí)際情況。雖然目前有眾多的數(shù)據(jù)模型、建模方法和應(yīng)用軟件，但他們各自的適應(yīng)能力有限，在解決復(fù)雜區(qū)域性模型構(gòu)建等方面仍面臨著較多亟待解決的難題，這些難題既涉及到地質(zhì)領(lǐng)域相關(guān)知識(shí)，同時(shí)也涉及到計(jì)算機(jī)、地質(zhì)統(tǒng)計(jì)學(xué)等學(xué)科相關(guān)知識(shí)。這就要求不同專業(yè)的相關(guān)人員進(jìn)行深層次的交流，共同推進(jìn)三維地質(zhì)建模事業(yè)的發(fā)展。

另一方面，國(guó)內(nèi)有眾多的研究人員和研究機(jī)構(gòu)進(jìn)行三維地質(zhì)建模理論研究，不同的機(jī)構(gòu)或研究人員研究思路、方法可能不同，但他們都各自積累了豐富的經(jīng)驗(yàn)，促進(jìn)這些研究機(jī)構(gòu)或研究人員進(jìn)行全方位合作，注重理論與實(shí)際相結(jié)合的研究，也是三維地質(zhì)建模發(fā)展的可行途徑。

參考文獻(xiàn)：

[1]Houlding S.3D geoscience modeling：computer techniques for geological characterization[M].Springer-Verlag，1994.

[2]韓國(guó)建，郭達(dá)志，金學(xué)林.礦體信息的八叉樹存儲(chǔ)和檢索技術(shù)[J].測(cè)繪學(xué)報(bào)，1992（1）：13-17.

[3]Victor J D，Alan P.Delaunay Tetrahedral Data Modelling for 3-D GIS Applications[J].GIS/LIS'93，1993：671-678.

[4]Pilouk M，Tempfli K，Molenaar M.A tetrahedron-based 3D vector data model for geo-information.[J].Agdm94 Spatial Data Modelling & Query Languages Ford Andd Applications M，1994.

[5]Wu L.Topological relations embodied in a generalized tri-prism（GTP）model for a 3D geoscience modeling system[J].Computers & Geosciences，2004，30（4）：405-418.

[6]Li R.Data structures and application issues in 3-D geographic information systems[J].Geomatica，1994，3.

[7]Xu N X，Tian H.Wire frame：A reliable approach to build sealed engineering geological models[J].Computers & Geosciences，2009，35（8）：1582-1591.

[8]Wu L X.Topological relations embodied in a generalized tri-prism（GTP）model for a 3D geoscience modeling system[J].Computers & Geosciences，2004，30（4）：405-418.

[9]Jrgensen F，Mller R R，Nebel L，et al.A method for cognitive 3D geological voxel modelling of AEM data[J].Bulletin of Engineering Geology & the Environment，2013，72（3）：421-432.

[10]潘懋，方裕，屈紅剛.三維地質(zhì)建模若干基本問題探討[J].地理與地理信息科學(xué)，2007，23（3）：1-5.

[11]Herbert M H，Jones C B，Tudhope D S.Three-dimensional reconstruction of geoscientific objects from serial sections[J].Visual Computer，1995，11（7）：343-359.

[12]Horsman J，Bethel W.Methods of constructing a 3D geological model from scatter data：Advanced visual systems conference，Boston，MA（United States），1995[C].

[13]Lemon A M，Jones N L.Building solid models from boreholes and user-defined cross-sections[J].Computers & Geosciences，2003，29（5）：547-555.

[14]He X，Koch J，Sonnenborg T O，et al.Transition probability‐based stochastic geological modeling using airborne geophysical data and borehole data[J].Water Resources Research，2014，50（4）：3147-3169.

[15]Hyer A S，Jrgensen F，F(xiàn)oged N，et al.Three-dimensional geological modelling of AEM resistivity data — A comparison of three methods[J].Journal of Applied Geophysics，2015，115：65-78.

[16]Zhu L F，Zhang C J，Li M J，et al.Building 3D solid models of sedimentary stratigraphic systems from borehole data：An automatic method and case studies[J].Engineering Geology，2012，127（1）：1-13.

[17]Jessell M.Three-dimensional geological modelling of potential-field data[J].Computers & Geosciences，2001，27（4）：455-4653.

[18]Calcagno P，Chilès J P，Courrioux G，et al.Geological modelling from field data and geological knowledge.Part I.Modelling method coupling 3D potential-field interpolation and geological rules[J].Physics of the Earth & Planetary Interiors，2008，171（1）：147-157.

[19]Kaufmann O，Martin T.3D geological modelling from boreholes，cross-sections and geological maps，application over former natural gas storages in coal mines[M].Pergamon Press，Inc.，2008.

[20]Wu Q，Xu H，Zou X.An effective method for 3D geological modeling with multi-source data integration[J].Computers & Geosciences，2005，31（1）：35-43.

[21]Li M C，Han Y Q，Liao Z，et al.Alternative 3D Modeling Approaches Based on Complex Multi-Source Geological Data Interpretation[J].Transactions of Tianjin University，2014，20（1）：7-14.

[22]李青元，張麗云，魏占營(yíng)，等.三維地質(zhì)建模軟件發(fā)展現(xiàn)狀及問題探討[J].地質(zhì)學(xué)刊，2013，37（4）：554-561.

[23]李瑞航，張永波，劉天霸，等.三維水文地質(zhì)建模技術(shù)研究綜述[J].南水北調(diào)與水利科技，2011，09（6）：98-101.

[24]Zhu L F，Zhuang Z Y.Framework system and research flow of uncertainty in 3D geological structure models[J].International Journal of Mining Science and Technology，2010，20（2）：306-311.

[25]Thore P，Shtuka A，Lecour M，et al.Structural uncertainties：Determination，management，and applications[J].Geophysics，2002，67（3）：840.

[26]卡爾斯杰夫.地球科學(xué)中的不確定性建模[M].北京：石油工業(yè)出版社，2016.

[27]Lindsay M，Ailleres L，Betts P，et al.Integrating geological uncertainty into combined geological and potential field inversions[J].2010：1-4.

[28]Gerritsen M，Caers J.Modeling Spatial and Structural Uncertainty in the Subsurface[M].Springer New York，2013.

[29]Wellmann J F，Horowitz F G，Schill E，et al.Towards incorporating uncertainty of structural data in 3D geological inversion[J].Tectonophysics，2010，490（34）：141-151.

[30]Yamamoto J K，Koike K，Kikuda A T，et al.Post-processing for uncertainty reduction in computed 3D geological models[J].Tectonophysics，2014，633：232-245.

[31]張曉坤，章浩.三維地質(zhì)建模在礦產(chǎn)資源開發(fā)中的應(yīng)用[J].金屬礦山，2010，V39（3）：106-110.

[32]Deutsch C V，Wang L.Hierarchical object-based stochastic modeling of fluvial reservoirs[J].Mathematical Geology，1996，28（7）：857-880.

[33]Wycisk P.3D Geological and Hydrogeological Modelling-Integrated Approaches in Urban Groundwater Management[J].Springer International Publishing，2015：3-12.

[34]Delhomme J P.Spatial variability and uncertainty in groundwater flow parameters：A geostatistical approach[J].Water Resources Research，1979，15（2）：269-280.

[35]Dong M，Neukum C，Hu H，et al.Real 3D geotechnical modeling in engineering geology：a case study from the inner city of Aachen，Germany[J].Bulletin of Engineering Geology & the Environment，2015，74（2）：281-300.

[36]Guglielmetti L，Comina C，Abdelfettah Y，et al.Integration of 3D geological modeling and gravity surveys for geothermal prospection in an Alpine region[J].Tectonophysics，2013，608（6）：1025-1036.

基金項(xiàng)目：水環(huán)所基本業(yè)務(wù)費(fèi)（SK201615）;雄安新區(qū)資源環(huán)境承載力綜合監(jiān)測(cè)和透明雄安數(shù)字平臺(tái)建設(shè)（DD20189144）

作者簡(jiǎn)介：姬廣軍（1993-），男，漢族，山東濟(jì)寧人，碩士，研究方向：區(qū)域地質(zhì)結(jié)構(gòu)建模和城市地質(zhì)信息化應(yīng)用研究。