趙磊磊 崔振東 張建勇

摘? ?要：對溶洞地質(zhì)進(jìn)行精確探測和三維地質(zhì)模型構(gòu)建具有重大意義?？偨Y(jié)了溶洞發(fā)育區(qū)內(nèi)不同類型復(fù)雜地質(zhì)結(jié)構(gòu)探測技術(shù)、建模思路和方法、數(shù)據(jù)插值擬合方法，并提出相關(guān)思考和建議。結(jié)論如下：①溶洞發(fā)育區(qū)內(nèi)地質(zhì)結(jié)構(gòu)體大致分為大型溶洞、小型溶洞、溶蝕孔洞、多尺度裂縫4種類型，分別采用確定性和隨機(jī)模擬方法建模，最終融合構(gòu)建成復(fù)雜三維地質(zhì)模型;②溶洞探測和建模應(yīng)注重多源數(shù)據(jù)整合和多方法集成，以實現(xiàn)地質(zhì)數(shù)據(jù)資料的相互彌補(bǔ)和建模方法的相互融合;③溶洞建模工作應(yīng)結(jié)合相關(guān)工程背景，以專業(yè)知識為基礎(chǔ)，以應(yīng)用目的為導(dǎo)向，調(diào)整建模思路，建立能反映歷史演化過程和未來演化趨勢的溶洞模型。

關(guān)鍵詞：溶洞;探測;三維地質(zhì)建模;插值擬合;建模軟件

溶洞是以石灰?guī)r和白云巖為主的可溶性巖石因地下水長期溶蝕所形成的地下空間，在我國較為發(fā)育，分布廣泛，我國巖溶面積約占國土面積的三分之一[1]。準(zhǔn)確認(rèn)識地下溶洞分布、模式、結(jié)構(gòu)、發(fā)育規(guī)律及對工程的影響，是關(guān)乎我國各類基礎(chǔ)設(shè)施工程及旅游開發(fā)項目設(shè)計、施工及運(yùn)營安全的大局。三維地質(zhì)建模已廣泛應(yīng)用于各專業(yè)領(lǐng)域，指導(dǎo)工程建設(shè)的設(shè)計和評估、有助于儲量的精確計算、實現(xiàn)平、剖面構(gòu)造形態(tài)相容并聯(lián)動修改、進(jìn)行三維空間分析與過程模擬、介紹和展示復(fù)雜地質(zhì)條件等[2-5]。目前國內(nèi)學(xué)者用簡化方式對含溶洞的地質(zhì)體進(jìn)行三維建模，如戴自航等將高速公路路堤下的溶洞簡化為橢球型，建立三維有限元模型，實現(xiàn)對溶洞頂板穩(wěn)定性的數(shù)值模擬分析[6];蘇濤等直接將溶洞簡化為球體進(jìn)行建模，計算分析隧道圍巖穩(wěn)定性[7]。本文基于大量文獻(xiàn)調(diào)研，歸納現(xiàn)階段典型溶洞探測技術(shù)，梳理溶洞發(fā)育區(qū)內(nèi)不同類型結(jié)構(gòu)體的建模方法，對比分析常用數(shù)據(jù)插值擬合方法和建模軟件，并結(jié)合溶洞建模的工程應(yīng)用需求提出思考和建議。

1? 溶洞探測技術(shù)

建立溶洞三維地質(zhì)模型，首先需采用探測技術(shù)獲取多種原始數(shù)據(jù)和地質(zhì)信息，目前常用探測方法見表1。

2? 溶洞三維地質(zhì)模型構(gòu)建方法

2.1? 現(xiàn)有三維地質(zhì)建模方法

常用的三維地質(zhì)建模方法大致分為確定性建模和隨機(jī)模擬建模方法兩大類[13]，具體分類見表2。

2.2? 溶洞三維地質(zhì)建模方法

傳統(tǒng)建模思路和流程見圖1[14]。

有學(xué)者提出“二步法”的建模思路[15]，第一步采用不同建模方法進(jìn)行單一類型模型構(gòu)建，第二步基于一定組合規(guī)律進(jìn)行融合。①大型溶洞（高度大于5 m）。一般采用地震屬性直接截斷方法，在成因控制約束下，采用確定性建模方法進(jìn)行建模[16];②小型溶洞（高0.5～5 m）。以同一巖溶帶內(nèi)相同成因的大型溶洞為訓(xùn)練圖像，采用多點地質(zhì)統(tǒng)計學(xué)隨機(jī)模擬方法進(jìn)行建模[17];③溶蝕孔洞（直徑2～500 mm）。采用序貫高斯模擬的隨機(jī)模擬方法對溶蝕孔洞儲集層進(jìn)行建模和表征[18];④多尺度裂縫型溶洞。采用分層次建模思想，按尺度大小進(jìn)行裂縫建模[19];⑤融合建模。建立大型溶洞、小型溶洞、溶蝕孔洞、裂縫4類結(jié)構(gòu)的離散模型后，按地質(zhì)分布和發(fā)育規(guī)律組合，形成地質(zhì)模型（圖2）[20]。

2.3? 數(shù)據(jù)插值擬合方法

建立三維地質(zhì)模型的關(guān)鍵是如何根據(jù)已知控制點數(shù)據(jù)內(nèi)插、外推已知資料點之間及以外的參數(shù)估計值。常用數(shù)據(jù)插值擬合方法包括離散光滑插值、克里金插值、反距離加權(quán)插值、自然鄰點插值。常用三維地質(zhì)建模軟件插值算法見表3。

3? 溶洞建模的思考和展望

注重多源數(shù)據(jù)的整合和應(yīng)用? 鉆井?dāng)?shù)據(jù)被稱為“硬數(shù)據(jù)”，地震資料被稱為“軟數(shù)據(jù)”，用來輔助判斷井間盲區(qū)部分的地質(zhì)特征。實際工作過程中，鉆孔數(shù)據(jù)失真，探測精度提高，使“硬數(shù)據(jù)變軟，軟數(shù)據(jù)變硬”。因此，注重多源數(shù)據(jù)整合應(yīng)用可起到相互彌補(bǔ)和印證作用，有利于提高建模精度。

重視多方法集成? 這里的“多方法”同時兼具探測方法和建模方法。探測方面據(jù)實際地質(zhì)條件，盡可能使用多種探測方法，實現(xiàn)數(shù)據(jù)資料的相互彌補(bǔ)和印證。建模方面根據(jù)不同地質(zhì)特征，采用不同建模方法構(gòu)建模型?？傊?，充分結(jié)合溶洞發(fā)育規(guī)律和屬性特征，不拘泥于尋求一種普適的建模方法，往往能起到事半功倍的效果。

充分結(jié)合應(yīng)用導(dǎo)向，調(diào)整建模思路? 基于溶洞建模的應(yīng)用目的，評估油氣藏儲存能力、旅游開發(fā)展示、評估溶洞影響等[28]。因此，充分結(jié)合應(yīng)用導(dǎo)向，調(diào)整建模思路，做到有的放矢對建設(shè)項目中溶洞建模及影響評估工作尤為重要。

精進(jìn)三維建模的屬性表征? 現(xiàn)階段溶洞建模多屬結(jié)構(gòu)建模，主要注重于描述溶洞的大小、位置、幾何形態(tài)等，對地質(zhì)體按統(tǒng)一屬性賦值。但越來越多的應(yīng)用要求賦予與實際相當(dāng)?shù)膶傩詤?shù)，因此，溶洞建模的屬性表征問題需亟待解決。

結(jié)合相關(guān)力學(xué)理論，實現(xiàn)動態(tài)模型構(gòu)建? 三維地質(zhì)建模已在各專業(yè)領(lǐng)域廣泛應(yīng)用?？山Y(jié)合相應(yīng)專業(yè)知識或工程需求，在建模過程中實現(xiàn)某一種或某幾種關(guān)鍵因素的動態(tài)影響，建立能反映地質(zhì)演化過程或演化趨勢的溶洞動態(tài)模型，推動專業(yè)應(yīng)用中的進(jìn)步與升級。

4? 結(jié)論

（1） 相對傳統(tǒng)鉆探技術(shù)，探測技術(shù)的應(yīng)用已使探測精度大大提高。但各種探測技術(shù)均有其應(yīng)用局限，需堅持多種探測技術(shù)的相互補(bǔ)充和印證。

（2） 巖溶發(fā)育結(jié)構(gòu)多樣、特征復(fù)雜，分為大型溶洞、小型溶洞、溶蝕孔洞、多尺度裂縫4種類型，分別采用確定性建模、隨機(jī)模擬建模等方法建立各類型模型，融合形成溶洞發(fā)育地質(zhì)體結(jié)構(gòu)模型。

（3） 對溶洞建模技術(shù)，應(yīng)結(jié)合工程應(yīng)用背景，以專業(yè)知識為基礎(chǔ)，應(yīng)用目的為導(dǎo)向，加強(qiáng)多源數(shù)據(jù)整合應(yīng)用、多方法集成，調(diào)整建模思路，建立能反映歷史演化過程和未來演化趨勢的溶洞模型。

參考文獻(xiàn)

[1]? ? 張炳暉，張研，王偉，等.基于高斯過程二元分類模型的溶洞規(guī)模預(yù)測方法[J].中國巖溶，2020，39（2）：259-263.

[2]? ? Kaufmann O， Martin T.3D geological modelling from boreholes， cross-sections sad geological maps， application over former natural gas storages in coal mines[J]. Computers & Geosciences， 2008，34（3）：278-290.

[3]? ? 楊輝廷. 縫洞型碳酸鹽巖儲層建模研究—以塔里木輪南油田奧陶系儲層為例[D].西南石油學(xué)院，2004.

[4]? ? 陳國旭，吳沖龍，張夏林，等.三維地質(zhì)建模與地礦勘查圖件編制一體化方法研究[J].地質(zhì)與勘探，2010，46（3）：542-546.

[5]? ? 姜作勤.從2D紙質(zhì)地質(zhì)圖到3D地質(zhì)模型-地質(zhì)信息表達(dá)的重大突破[C].中國地質(zhì)學(xué)會云南省地質(zhì)學(xué)會，全國地質(zhì)制圖與GIS學(xué)術(shù)論壇論文集.2011：71-78.

[6]? ? 戴自航，范夏玲，盧才金.巖溶區(qū)高速公路路堤及溶洞頂板穩(wěn)定性數(shù)值分析[J].巖土力學(xué)，2014，35（S1）：382-390.

[7]? ? 蘇濤，廖秀宇，趙樺，等.多隱伏溶洞對隧道圍巖穩(wěn)定性影響的數(shù)值模擬分析[J].內(nèi)蒙古科技與經(jīng)濟(jì)，2018（5）：68-69.

[8]? ? Kaufmann， Y.， Quinif， O.， Geohazard map of cover-collapse sinkholes in the “Tournaisis” area southern Belgium. Eng. Geol. 2002， 65， 117-124.

[9]? ? Carbonel， D.， Rodríguez-Tribaldos， V.， Gutiérrez， F.， Galve， J.P.， Guerrero， J.， Zarroca， M.，Roqué， C.， Linares， R.， McCalpin， J.P.， Acosta， E.， Investigating a damaging buried sinkhole cluster in an urban area （Zaragoza city， NE Spain） integrating multiple techniques： geomorphological surveys， DInSAR， DEMs， GPR， ERT， and trenching. Geomorphology. 2014，229， 3-16.

[10]? Martel， R.， Castellazzi， P.， Gloaguen， E.， ERT， GPR， InSAR， and tracer tests to characterize karst aquifer systems under urban areas： the case of Quebec City. Geomorphology. 2018，310， 45-56.

[11]? Mahmut， G.D.， Identification of gypsum karstification using an electrical resistivity tomography technique： The case-study of the Sivas gypsum karst area （Turkey）. Eng. Geol. 2019，252， 78-98.

[12]? Houlding S. 3D geoscience modeling： computer techniques for geological characterization[M]. Springer Science & Business Media， 1994.

[13]? 楊輝廷，江同文，顏其彬，等.縫洞型碳酸鹽巖儲層三維地質(zhì)建模方法初探[J].大慶石油地質(zhì)與開發(fā)，2004（4）：11-12+16-89.

[14]? 劉文祥. 昆明某建筑場地地基巖溶發(fā)育規(guī)律及三維可視化研究[D].昆明理工大學(xué)，2019.

[15]? 胡向陽，權(quán)蓮順，侯加根.巖溶相控多尺度碳酸鹽巖縫洞儲集體二步法建模新技術(shù)[J].科技導(dǎo)報，2013，31（24）：42-46.

[16]? 王光付. 碳酸鹽巖溶洞型儲層綜合識別及預(yù)測方法[J]. 石油學(xué)報，2008， 29（1）： 47-51.

[17]? 侯加根， 馬曉強(qiáng)， 胡向陽， 等. 碳酸鹽巖溶洞型儲集體地質(zhì)建模的幾個關(guān)鍵問題[J]. 高校地質(zhì)學(xué)報， 2013， 19（1）： 64-69.

[18]? 胡向陽， 袁向春， 侯加根， 等. 多尺度巖溶相控碳酸鹽巖縫洞型油藏儲集體建模方法[J]. 石油學(xué)報， 2014， 35（2）： 340-346.

[19]? 趙彬， 侯加根， 劉鈺銘，等. 基于示性點過程模擬的碳酸鹽巖裂隙型儲層建模方法[J]. 科技導(dǎo)報， 2011， 29（3）： 39-43.

[20]? 李陽，侯加根，李永強(qiáng).碳酸鹽巖縫洞型儲集體特征及分類分級地質(zhì)建模[J].石油勘探與開發(fā)，2016，43（4）：600-606.

Karst Cave Exploration and 3D Geological Model Construction Method

Zhao Leilei 1， 2， Cui Zhendong 1， 3， 4， Zhang Jianyong 1， 3， 4

（1. Institute of Geology and Geophysics， Chinese Academy of Sciences;Key Laboratory of shale gas and geological engineering， Chinese Academy of Sciences，Beijing， 100029，China;2. School of mechanics and architectural engineering， China University of mining and Technology （Beijing）， Beijing， 100083，China;3. Institute of Earth Sciences，

Chinese Academy of Sciences， Beijing， 100029，China;4. School of earth and Planetary Sciences，

University of Chinese Academy of Sciences， Beijing， 100049，China）

Abstract：It is of great significance to carry out accurate exploration of karst cave geology and construct three-dimensional geological model. The detection technologies， modeling ideas and methods， data interpolation fitting methods of different types of complex geological structures of karst caves are summarized， and some related thoughts and suggestions are put forward. The conclusions are as follows： （1） Geological structures in the karst cave development area can be roughly divided into four types： large karst cave， small karst cave， dissolution cave and multi-scale fracture. The modeling methods of deterministic and stochastic simulation are adopted to model the geological structures， respectively. Finally， the complex geological model of the karst cave development area is constructed by fusing them.（2） Multi-source data integration and multi-method integration should be emphasized in karst cave exploration and modeling to realize the mutual complement of geological data and the mutual integration of modeling methods.（3） The karst cave modeling should be combined with the corresponding engineering application background， based on professional knowledge and guided by application purpose， adjust the modeling ideas， establish a karst cave model that can reflect the historical evolution process and future evolution trend.

Key words：Karst cave; Exploration; 3D geological modeling; Interpolation fitting; Modeling software