劉猛 陸海玉 朱永峰 王帥

文章編號：1671-3559（2024）03-0306-06DOI：10.13349/j.cnki.jdxbn.20230607.001

摘要： 為了避免喀斯特塌陷對孫祖河橡膠壩工程項目實施造成威脅，通過野外調(diào)查、 工程地質(zhì)測繪、 鉆探、 現(xiàn)場試驗、 穩(wěn)定性分析等手段，系統(tǒng)分析研究區(qū)喀斯特塌陷的特征，探討喀斯特塌陷的成因演化機制，查明喀斯特塌陷的誘發(fā)因素，并對研究區(qū)喀斯特塌陷的穩(wěn)定性及其演化趨勢進行預(yù)測。結(jié)果表明，研究區(qū)目前處于土洞形成階段，外界水動力條件的改變是研究區(qū)喀斯特塌陷穩(wěn)定性的主要影響因素，具備進一步發(fā)展形成喀斯特塌陷的空間條件、 物質(zhì)條件和水動力條件，在進行地下水開采時應(yīng)控制地下水位最大降深。

關(guān)鍵詞： 工程地質(zhì)學(xué)； 喀斯特塌陷； 土洞； 水動力條件； 孫祖河

中圖分類號： P642

文獻標志碼： A

開放科學(xué)識別碼（OSID碼）：

Karst Collapse Characteristics and Evolution Mechanism of

the Sunzu River Rubber Dam

LIU Meng， LU Haiyu， ZHU Yongfeng， WANG Shuai

（Shandong Survey and Design Institute of Water Conservancy Co.， Ltd.， Jinan 250013， Shandong， China）

Abstract： To avoid threat of karst collapse to implementation of the Sunzu River rubber dam project， characteristics of karst collapse in the research area were systematically analyzed by means of field survey， engineering geological mapping， drilling， field test， and stability analysis， so as to explore genetic evolution mechanism of karst collapse and identify inducing factors of karst collapse. Stability and evolution trend of karst collapse in the research area were predicted. The results show that the research area is in the stage of soil cave formation at present， and the change of external hydrodynamic conditions is the main influencing factor of the stability of karst collapse in the research area. There are space， material， and hydrodynamic conditions for further development of karst ground collapses， and the maximum groundwater level depth should be controlled during groundwater exploitation.

Keywords： engineering geology; karst collapse; soil cave; hydrodynamic condition; the Sunzu River

喀斯特塌陷（土洞）是指在有覆蓋土層的喀斯特發(fā)育區(qū)， 在自重應(yīng)力和地下水等動力因素影響下， 覆蓋層土體遭到流失遷移而形成土中洞穴并不斷向地面擴張， 最后土洞頂板失穩(wěn)產(chǎn)生塌落或引發(fā)地面變形沉陷的一種不良地質(zhì)現(xiàn)象［1］。 水利工程建設(shè)與地質(zhì)環(huán)境的聯(lián)系非常密切， 特別是在喀斯特塌陷區(qū)， 喀斯特地質(zhì)環(huán)境問題已經(jīng)成為制約水利工程建設(shè)的關(guān)鍵， 嚴重影響水利工程建設(shè)進度， 制約社會經(jīng)濟發(fā)展。

國內(nèi)外學(xué)者對于喀斯特塌陷成因演化機制進行了廣泛的研究，提出了多種喀斯特塌陷理論機制。20世紀90年代以前主要是系統(tǒng)理論研究階段，主要

收稿日期： 2023-02-28????????? 網(wǎng)絡(luò)首發(fā)時間：2023-06-08T15：33：26

基金項目： 國家自然科學(xué)基金項目（42002257）

第一作者簡介： 劉猛（1982—），男，山東菏澤人。高級工程師，研究方向為巖土工程勘測。E-mail： 809402205@qq.com。

通信作者簡介： 陸海玉（1978—），男，山東汶上人。高級工程師，碩士，研究方向為工程地質(zhì)。E-mail： 631700378@qq.com。

網(wǎng)絡(luò)首發(fā)地址： https：//kns.cnki.net/kcms2/detail/37.1378.N.20230607.1737.002.html

是對喀斯特塌陷的特征及其形成機制進行分析，對不同喀斯特塌陷間聯(lián)系與區(qū)別進行甄別歸納分類，從而形成系統(tǒng)理論。徐衛(wèi)國等［2］針對采礦排水、 礦井突水、 水資源開發(fā)等導(dǎo)致的喀斯特塌陷問題提出了真空吸蝕論?？祻┤剩?］按塌陷產(chǎn)生的主導(dǎo)因素和受力狀態(tài)，將喀斯特地面塌陷劃分為重力塌陷、 潛蝕塌陷、 沖爆塌陷、 真空吸蝕塌陷、 振動塌陷、 荷載塌陷、 溶蝕塌陷和根蝕塌陷8種類型。程星等［4］以工程實例為基礎(chǔ)將喀斯特塌陷概化為單一阻水型蓋層地質(zhì)概化模型、 單一透水型蓋層地質(zhì)概化模型、 無蓋層地質(zhì)概化模型、 阻-透型蓋層地質(zhì)概化模型、透-阻型蓋層地質(zhì)概化模型、 阻-透-阻型蓋層地質(zhì)概化模型與透-阻-透蓋層地質(zhì)概化模型等6種概化模型。羅小杰等［5-6］從土體塌陷方式著手提出了土洞型塌陷、 沙漏型塌陷和泥流型塌陷的3種喀斯特塌陷理論。上述研究極大地促進了喀斯特塌陷的理論的發(fā)展。20世紀90年代以后，國內(nèi)學(xué)者對于喀斯特塌陷的特征、 形成機制及發(fā)育規(guī)律等已經(jīng)有了較系統(tǒng)的認識，加上國家的大力支持，各種新技術(shù)、 新方法得到應(yīng)用，對于喀斯特塌陷的研究轉(zhuǎn)入試驗、 模擬、 監(jiān)測及防治等方向。陶小虎等［7］從一維地下水運動和滲透力學(xué)的角度，分析比較潛水位上升與承壓水位下降對喀斯特地區(qū)透-阻型蓋層中阻水層滲透穩(wěn)定性的影響，認為潛水位上升或者承壓水位下降均不利于喀斯特的穩(wěn)定。熊啟華等［8］通過物理模型試驗探索負壓作用下喀斯特塌陷幾何形貌演化過程，認為負壓增大或土體強度參數(shù)弱化，都將降低土拱穩(wěn)定性。張藝凡等［9］對鉆孔誘發(fā)二元結(jié)構(gòu)覆蓋型喀斯特塌陷穩(wěn)定性進行預(yù)測，結(jié)果表明，上覆蓋層越薄、 覆蓋層中砂土層厚度與黏性土層厚度的比值（砂黏比）越大、 黏土抗剪強度越小，喀斯特的穩(wěn)定性越差。高宗軍等［10］提出基于喀斯特水動態(tài)監(jiān)測的設(shè)想，通過對喀斯特水位、 水量、 濁度以及主要化學(xué)組分等指標進行動態(tài)監(jiān)測以實現(xiàn)對喀斯特塌陷的預(yù)測預(yù)報工作。孫金輝［11］通過模擬地下水波動、 降雨入滲以及真空吸蝕等不同影響因素的組合作用下黏土路基的變形響應(yīng)情況，得出覆蓋型喀斯特塌陷產(chǎn)生的先決因素是喀斯特底部裂隙寬度達到了臨界裂隙寬度。

學(xué)者們從多個方面對喀斯特塌陷進行研究，在致塌模式、 致塌作用力及塌陷演化規(guī)律等方面取得了一系列成果，提出了喀斯特塌陷的形成機制和影響方式，極大地豐富了喀斯特塌陷理論研究。由于喀斯特塌陷地質(zhì)條件的復(fù)雜性和喀斯特塌陷的突發(fā)性，在喀斯特塌陷發(fā)育地段，如何評判當前喀斯特塌陷所處的階段，如何對其誘發(fā)因素進行控制，并對其下一步的發(fā)展趨勢進行預(yù)測、 預(yù)報及防治，相關(guān)研究涉獵較少。本文中以山東省沂南縣孫祖河喀斯特塌陷為例，通過野外調(diào)查、 工程地質(zhì)測繪、 鉆探、 現(xiàn)場試驗、 穩(wěn)定性分析等手段，對孫祖河橡膠壩工程區(qū)可能形成的喀斯特塌陷的空間條件、 物質(zhì)條件和水動力條件3個方面進行聯(lián)合系統(tǒng)分析，判斷了孫祖橡膠壩工程區(qū)喀斯特塌陷的當前所處的發(fā)展階段，得出了其誘發(fā)因素，預(yù)測了其下一步演化趨勢，進一步完善喀斯特塌陷研究理論體系，為喀斯特塌陷工程實踐及監(jiān)測預(yù)警工作提供參考。

1? 研究區(qū)概況

孫祖河是東汶河流域下游的支流河道， 發(fā)源于山東省沂南縣孫祖鎮(zhèn)西北丘陵區(qū)， 在張莊鎮(zhèn)前漢沿村匯入東汶河， 河道全長20 km，流域面積為120 km2，河道平均比降為1.7‰，屬典型山區(qū)型河流。

研究區(qū)地貌屬于魯中南強—弱切割構(gòu)造剝蝕中低山丘陵亞區(qū)，主要包括丘陵地貌以及河谷地貌，總體地勢為西北高、 東南低，兩側(cè)山脊高、 中間河谷低。河道兩岸沿線地貌形態(tài)主要為侵蝕剝蝕型丘陵，以構(gòu)造剝蝕為主，長期受強—弱度剝蝕切割，切割深度為30～80 m。丘陵大多山勢低緩，丘頂渾圓呈饅頭狀，巖石直接出露。孫祖河河谷呈不對稱的U型，河漫灘較發(fā)育，河谷及漫灘地層巖性主要以淤泥、 礫質(zhì)粗砂為主，局部沉積有少量壤土，第四系覆蓋層厚度為7～11 m。

流域內(nèi)出露的主要地層單元為寒武系中統(tǒng)饅頭組C-2-3 m、 寒武系中統(tǒng)朱砂洞組C-2z^，地層巖性以灰?guī)r為主，局部夾砂巖及頁巖，巖層產(chǎn)狀為59°∠10°。

研究區(qū)在地質(zhì)構(gòu)造上位于魯西隆起區(qū)（Ⅱ）—魯中隆起（Ⅱa）—馬牧池—沂源斷?。á騛7）—馬牧池凸起（Ⅱ3a7），主要斷裂有新泰—蒙陰斷裂、 銅冶店—孫祖斷裂、 鄌郚—葛溝斷裂。

區(qū)內(nèi)年均降水差異較大，年內(nèi)分布極不均勻，隨季度變化明顯，74.7%的降水集中在每年汛期6—9月份。主要含水層為礫質(zhì)粗砂及灰?guī)r，區(qū)內(nèi)喀斯特水徑流通暢，受降水和開采影響，動態(tài)變化較大。

2? 喀斯特塌陷特征分析

在對山東省沂南縣孫祖水源工程孫祖河橡膠壩勘察時， 在距橡膠壩上游5.0 km處發(fā)現(xiàn)有一喀斯特地面塌陷坑（見圖1）， 該塌陷坑直徑12 m， 坑深4 m，坑底見水。 鉆探揭示壩址區(qū)下部基巖喀斯特強發(fā)育， 存在半徑5.0～7.0 m的喀斯特溶洞。

對于喀斯特塌陷的特征需要從多方面進行分

析，一般認為空間條件、 物質(zhì)條件和水動力條件是喀斯特地面塌陷發(fā)生的3個先決條件［12］?？λ固厮?/p>

陷發(fā)育的空間條件主要是指可溶巖裂隙發(fā)育， 即基巖面存在溶洞、? 溶溝、 裂隙等， 是喀斯特塌陷產(chǎn)生的必要條件。 喀斯特塌陷的物質(zhì)條件主要是指蓋層結(jié)構(gòu)， 決定了喀斯特塌陷的類型， 對喀斯特塌陷蓋層結(jié)構(gòu)的深入了解是正確認識喀斯特塌陷災(zāi)害成因及機制的前提。 喀斯特塌陷的水動力條件主要是指喀斯特水對蓋層的滲透、 沖刷和軟化作用等，是喀斯特塌陷產(chǎn)生的原動力。要想分析喀斯特塌陷的形成、 演化機制對塌陷區(qū)基巖喀斯特及裂隙的發(fā)育程度的影響，蓋層結(jié)構(gòu)及其水動力條件的研究必不可少。

2.1? 下伏可溶巖發(fā)育特征

研究區(qū)下伏基巖主要為寒武系中統(tǒng)饅頭組C-2-3 m以及寒武系中統(tǒng)朱砂洞組C-2z^地層， 地層巖性以灰?guī)r為主， 局部夾砂巖及頁巖， 鉆探、 物探及現(xiàn)場工程地質(zhì)測繪均顯示， 研究區(qū)喀斯特現(xiàn)象發(fā)育， 能夠為喀斯特塌陷提供空間條件。 研究區(qū)可溶巖喀斯特發(fā)育特征見表1。

2.2? 蓋層結(jié)構(gòu)

研究區(qū)地層巖性特征如圖2所示。 鉆探資料顯示， 孫祖河橡膠壩處河床及漫灘內(nèi)第四系地層主要為沖積堆積的礫質(zhì)粗砂， 局部表層為薄層黏性土， 整體厚度為7～11 m， 在礫質(zhì)粗砂和可溶巖間無隔水層， 第四系孔隙水與喀斯特水聯(lián)系密切。 砂層與下伏可溶巖直接接觸， 可溶巖頂板較薄， 鉆探揭示溶洞發(fā)育部位， 頂板厚度為0.2 m， 可見研究區(qū)喀斯特頂板極薄， 為單一透水性蓋層結(jié)構(gòu)， 該蓋層結(jié)構(gòu)易發(fā)沙漏型喀斯特地面塌陷。

2.3? 水動力條件

研究區(qū)水文地質(zhì)特征如圖3所示。孫祖河為山溪雨源型河流，位于孫祖斷塊水文地質(zhì)單元內(nèi)，汛期水量大，非汛期水量小，主要的地下水類型有松散巖類孔隙水、 碎屑巖類孔隙裂隙水、 碳酸鹽巖巖溶裂隙水、 塊狀巖類裂隙水4類［13］。碳酸鹽巖巖溶裂隙水與松散巖類孔隙水間聯(lián)系密切，無隔水層，受降水影響明顯，地下水由兩岸向?qū)O祖河河谷匯流。調(diào)查發(fā)現(xiàn)，研究區(qū)內(nèi)喀斯特水長期被周邊采砂場的過量開采，喀斯特水的開采不利于喀斯特地區(qū)的穩(wěn)定，利于喀斯特塌陷的形成。在枯水期當補給量達不到開采量時，地下水位急劇下降，此時承壓水會接受松散巖類孔隙水的補給，松散巖類孔隙水沿裂隙或喀斯特通道下滲，對灰?guī)r產(chǎn)生溶蝕、 水解、 軟化等多種作用，破壞喀斯特頂板，加劇喀斯特現(xiàn)象的發(fā)育。汛期降水量增大，喀斯特地下水位上升，當水位上升至喀斯特頂板時，喀斯特水會對上覆巖土產(chǎn)生滲壓效應(yīng)。滲壓效應(yīng)是降雨或地表水等上部水體補給喀

斯特水時對上覆巖土產(chǎn)生的綜合作用［14］，可為喀斯特地面塌陷提供充分的水動力條件。

伴隨著汛期、 枯水期的交替， 研究區(qū)巖土體常年受干濕循環(huán)條件的影響。 有研究表明， 隨著干濕循環(huán)次數(shù)的增加， 土體強度逐漸減小， 且在第1次干濕循環(huán)過程中土體強度減小尤為顯著［15］。 若喀斯特頂板被擊穿， 上覆松散細顆粒即會在地下水的滲流、 軟化等作用下不斷流失， 當上覆土體所受致塌力大于抗塌力時， 蓋層失穩(wěn)， 喀斯特塌陷發(fā)生。

3? 喀斯特塌陷成因及演化機制

根據(jù)研究區(qū)的空間條件、物質(zhì)條件和水動力條件，引起土體坍塌的原因主要有水動力作用效應(yīng)、 土力學(xué)作用效應(yīng)和氣動作用效應(yīng)［16］。研究區(qū)喀斯特塌陷主要經(jīng)歷潛蝕軟化、 土拱形成、 地面塌陷3個階段，如圖4。

1）潛蝕軟化階段。研究區(qū)下部基巖喀斯特發(fā)育程度高，溶孔、 溶洞發(fā)育，為地下水及砂性土提供了良好的空間條件和運移通道。區(qū)內(nèi)孔隙潛水和喀斯特水水力聯(lián)系密切，不存在相對隔水層，喀斯特水與孔隙潛水間補給、排泄順暢?？λ固厮艿竭^量開采及降水的影響，喀斯特水上下波動及水平流動，不斷沖刷上覆巖土體，巖土體強度逐漸降低，此時土體受到喀斯特水的浮力和土體自身抗滑力的作用，土體處于穩(wěn)定狀態(tài)。

2）土洞形成階段。當上覆喀斯特頂板被擊穿時，砂性土直接鄰水，砂性土中的松散細顆粒隨著地下水的快速流動而不斷流失， 隨著喀斯特水水位下

降，潛水滲透產(chǎn)生破壞作用，土洞逐漸形成。

3）地面塌陷階段。 隨著土洞的逐漸發(fā)展， 當上覆土層逐漸失去支撐時， 在土體自重、 真空吸力、 潛水滲透力和喀斯特水浮力減小等因素的作用下， 致塌力大于抗塌力， 土體失穩(wěn)， 喀斯特地面塌陷發(fā)生。

4? 喀斯特塌陷穩(wěn)定性分析

4.1? 荷載分析

根據(jù)極限平衡理論，喀斯特塌陷土體主要受到向下的致塌力和向上的抗塌力作用。其中致塌力包括塌陷土體的自重G、 雨水下滲滲透力Fs、 真空負壓力Fp、 抽取地下水引起地下水浮力增量Fw以及地表外力等； 抗塌力主要是指塌落體的側(cè)摩阻力f。當土體的抗塌力小于作用于土洞的致塌力時，溶洞出現(xiàn)塌陷，具體土體受力分析如圖5所示。

圖中箭頭表示受力方向； f為側(cè)摩阻力；

G為土體自重力； Fs為滲透力； Fp真空負壓力；

Fw為抽取地下水引起地下水浮力增量。

塌陷土體的自重力計算公式為

Ｇ=πD2 hγ4 ，（1）

式中： γ為土體的容重，地下水位以上使用飽和容重，地下水位以下使用天然容重； h為覆蓋土體厚度； D為塌落拱或漏斗直徑。

由降水入滲引起的滲透力計算公式為

Fs=πD2v2γ08g ，（2）

式中： γ0為水的容重； v為地下水流動速度； g為重力加速度。

真空負壓力計算公式為

Fp=pA=πpD24 ，（3）

式中： p為落水洞與漏斗內(nèi)的真空負壓差，取值為0～50 kPa［17］； A為真空負壓作用于塌陷土體的面積。

抽取地下水引起地下水浮力增量計算公式為

Fw=πD2γ0Δhw4 ，（4）

式中Δhw為地下水下降幅度。

地表外力主要為新建橡膠壩壩體荷載，當前地表外力為0。

土洞上方圓柱形塌落體的側(cè)摩阻力計算公式［18］為

f=πDk0γh22 tan φ+ch ，（5）

式中： k0為靜止側(cè)壓力系數(shù)， 取為0.6； φ為礫質(zhì)粗砂內(nèi)摩擦角； c為礫質(zhì)粗砂的黏聚力。

4.2? 穩(wěn)定性分析

將塌陷土體視為剛體，則在極限平衡條件下，塌陷土體的穩(wěn)定系數(shù)計算公式為

k=2k0 γh2 tan φ+4chD（hγ+v2γ0/2g+Δhwγ0+p）+Fl ，（6）

式中Fl為塌陷土體承受的外部荷載。

塌落拱直徑D計算公式為

D=d+2R tan（45°-φr/2），（7）

式中： d為鉆孔直徑； R為溶洞頂、 底板高度差； φr為灰?guī)r內(nèi)摩擦角。

礫質(zhì)粗砂物理力學(xué)參數(shù)見表2?？辈炱陂g正直汛期，潛水位距離現(xiàn)狀地面以下1.8 m，礫質(zhì)粗砂的加權(quán)容重γ取20.0 kN/m3，根據(jù)實際鉆探及物探資料，研究區(qū)礫質(zhì)粗砂覆蓋層厚度為8.8 m，區(qū)內(nèi)喀斯特覆蓋層頂板已被鉆孔揭穿，鉆孔直徑為0.091 m。鉆孔揭示灰喀斯特洞頂、 底板高度差為5.0～7.0 m，灰?guī)r內(nèi)摩擦角φr取60°，地下水流動速度為3 m/s。

表2? 礫質(zhì)粗砂物理力學(xué)參數(shù)

容重/（kN/m3）天然容重飽和容重

抗剪性能

天然黏聚力/kPa天然內(nèi)摩擦角/（°）

16.820.8030

由此計算得塌落拱最大直徑D為3.84 m。 若不考慮外界因素， 僅考慮土體自重條件下， 穩(wěn)定性系數(shù)為1.58， 此時研究區(qū)處于穩(wěn)定狀態(tài)， 因此， 外界水動力條件的改變是影響該喀斯特塌陷穩(wěn)定性的主要影響因素。 外界水動力條件的改變主要受地下水開采的影響， 因此， 真空負壓分別取0、 25、 50 kPa時，該喀斯特塌陷穩(wěn)定性系數(shù)隨地下水位降深關(guān)系見圖6。 由圖可知， 在真空負壓取50 kPa的最不利條件下， 此時地下水位降深為4.88 m， 該喀斯特塌陷達到極限平衡狀態(tài)， 當?shù)叵滤_采降深超過4.88 m時， 塌陷土體即會失穩(wěn)， 因此在抽取地下水時， 為了減緩或避免覆蓋型塌陷的發(fā)生，應(yīng)控制好地下水位最大降深。

4.3? 塌陷預(yù)測

研究區(qū)現(xiàn)狀地下水位降深小于4.88 m，經(jīng)喀斯特塌陷穩(wěn)定性計算，該塌陷點尚處于穩(wěn)定狀態(tài)，但現(xiàn)場調(diào)查發(fā)現(xiàn)，研究區(qū)周邊村民自家打井抽取的承壓水每年都會一段時間出現(xiàn)渾濁的現(xiàn)象，推測為在喀斯特水的作用下砂性土產(chǎn)生的小規(guī)模塌陷。由此可推斷，孫祖河橡膠壩工程區(qū)當前處于土洞形成階段。

根據(jù)普氏地壓理論，在任何深度的巖層中形成土洞，土洞上部會形成壓力拱，以承受拱上部的覆蓋巖層重量而不致于坍塌，土拱高度計算公式［19］為

h0=0.828D2fk ，（8）

式中fk為覆蓋層土體的堅固性系數(shù)， 礫質(zhì)粗砂取0.7。

經(jīng)推算，研究區(qū)可形成的塌陷土拱最大高度為2.27 m。

5? 結(jié)論

本文中對孫祖河橡膠壩工程區(qū)可能形成的喀斯特塌陷的空間條件、 物質(zhì)條件和水動力條件等3個方面進行聯(lián)合系統(tǒng)分析，得到以下主要結(jié)論：

1）根據(jù)研究區(qū)喀斯特塌陷的空間條件、 物質(zhì)條件和水動力條件，認為研究區(qū)喀斯特塌陷主要經(jīng)歷潛蝕軟化、 土洞形成、 地面塌陷3個階段，且該塌陷點當前處于土洞正在形成階段，并具備進一步發(fā)展的趨勢。

2）研究區(qū)土體穩(wěn)定性分析結(jié)果表明，外界水動力條件的改變是影響研究區(qū)喀斯特塌陷的誘發(fā)因

素，在最不利條件下，研究區(qū)地下水最大開采降深為4.88 m，為了減緩或避免喀斯特塌陷的發(fā)生，應(yīng)控制

好地下水開采最大降深。

3）研究區(qū)若發(fā)生喀斯特塌陷，其塌陷類型為沙漏型，形成的塌陷土拱最大高度為2.27 m。

參考文獻：

［1］? 賈龍， 蒙彥， 管振德. 巖溶土洞演化及其數(shù)值模擬分析［J］. 中國巖溶， 2014， 33（3）： 294.

［2］? 徐衛(wèi)國， 趙桂榮. 試論巖溶礦區(qū)地面塌陷的真空吸蝕作用［J］. 地質(zhì)論評， 1981， 27（2）： 174.

［3］? 康彥仁. 試論巖溶地面塌陷的類型劃分［J］. 中國巖溶， 1984， 3（2）： 146.

［4］? 程星， 黃潤秋. 巖溶塌陷的地質(zhì)概化模型［J］.水文地質(zhì)工程地質(zhì)， 2002（6）： 30.

［5］? 羅小杰， 羅程. 沙漏型巖溶地面塌陷物理模型［J］. 中國巖溶， 2017， 36（1）： 88.

［6］? 羅小杰， 羅程. 巖溶地面塌陷三機理理論及其應(yīng)用［J］. 中國巖溶， 2021， 40（2）： 171.

［7］? 陶小虎， 趙堅， WANG Xiaoming， 等. 地下水位變化對透-阻型巖溶塌陷影響的分析［J］. 中國巖溶， 2017， 36（6）： 777.

［8］? 熊啟華， 高旭， 涂婧， 等. 負壓作用下土洞型巖溶塌陷機理及力學(xué)模型研究［J］. 人民長江， 2022， 53（9）： 163.

［9］? 張藝凡， 徐光黎， 胡志新， 等. 鉆孔誘發(fā)二元結(jié)構(gòu)覆蓋型巖溶塌陷預(yù)測［J］. 工程勘察， 2021， 49（2）： 7.

［10］? 高宗軍， 魯統(tǒng)民， 王敏， 等. 基于巖溶水動態(tài)的巖溶地面塌陷預(yù)測預(yù)報方法［J］. 中國巖溶， 2019， 38（5）： 739.

［11］? 孫金輝. 覆蓋型巖溶塌陷臨界參數(shù)模型試驗與數(shù)值模擬研究［D］. 成都： 西南交通大學(xué)， 2011： 49.

［12］? 葉青. 干濕循環(huán)作用下巖溶區(qū)土洞塌陷機理研究［D］. 淮南： 安徽理工大學(xué)， 2021： 21.

［13］? 李常鎖， 尚宇寧， 張海林， 等. 山東省1∶5萬區(qū)域水文地質(zhì)調(diào)查報告（沂南縣幅）［R］. 濟南： 山東省地礦工程勘察院， 2016.

［14］? 黃啟云， 朱川曲， 李青鋒， 等. 山腳樹礦淺表巖溶塌陷誘因及機理分析［J］. 湖南科技大學(xué)學(xué)報（自然科學(xué)版）， 2021， 36（3）： 6.

［15］? 楊和平， 肖奪. 干濕循環(huán)效應(yīng)對膨脹土抗剪強度的影響［J］. 長沙理工大學(xué)學(xué)報（自然科學(xué)版）， 2005， 2（2）： 1.

［16］? 高培德， 王林峰. 覆蓋型巖溶塌陷的塌陷機制分析［J］. 中國巖溶， 2017， 36（6）： 771.

［17］? 賀可強， 王濱， 萬繼濤. 棗莊巖溶塌陷形成機理與致塌模式研究［J］.巖土力學(xué)， 2002， 23（5）： 569.

［18］? 魯亞楠， 李炳權(quán). 韶關(guān)港覆蓋型巖溶工程地質(zhì)特征及巖溶穩(wěn)定性評價［J］. 中國勘察設(shè)計， 2021（8）： 101.

［19］? 王濱， 李治廣， 董昕， 等. 巖溶塌陷的致塌力學(xué)模型研究：以泰安市東羊婁巖溶塌陷為例［J］. 自然災(zāi)害學(xué)報， 2011， 20（4）： 122.

（責任編輯：于海琴）