孫祖河橡膠壩喀斯特塌陷特征與演化機(jī)制

劉猛 陸海玉 朱永峰 王帥

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

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

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

中圖分類號(hào)： P642

文獻(xiàn)標(biāo)志碼： A

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

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

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

基金項(xiàng)目： 國(guó)家自然科學(xué)基金項(xiàng)目（42002257）

第一作者簡(jiǎn)介： 劉猛（1982—），男，山東菏澤人。高級(jí)工程師，研究方向?yàn)閹r土工程勘測(cè)。E-mail： 809402205@qq.com。

通信作者簡(jiǎn)介： 陸海玉（1978—），男，山東汶上人。高級(jí)工程師，碩士，研究方向?yàn)楣こ痰刭|(zhì)。E-mail： 631700378@qq.com。

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

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

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

1? 研究區(qū)概況

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

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

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

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

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

2? 喀斯特塌陷特征分析

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

對(duì)于喀斯特塌陷的特征需要從多方面進(jìn)行分

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

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

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

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

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

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

2.3? 水動(dòng)力條件

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

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

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

3? 喀斯特塌陷成因及演化機(jī)制

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

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

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

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

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

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

4.1? 荷載分析

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

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

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

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

塌陷土體的自重力計(jì)算公式為

Ｇ=πD2 hγ4 ，（1）

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

由降水入滲引起的滲透力計(jì)算公式為

Fs=πD2v2γ08g ，（2）

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

真空負(fù)壓力計(jì)算公式為

Fp=pA=πpD24 ，（3）

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

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

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

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

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

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

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

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

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

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

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

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

塌落拱直徑D計(jì)算公式為

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

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

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

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

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

抗剪性能

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

16.820.8030

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

4.3? 塌陷預(yù)測(cè)

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

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

h0=0.828D2fk ，（8）

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

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

5? 結(jié)論

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

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

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

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

好地下水開(kāi)采最大降深。

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

參考文獻(xiàn)：

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

（責(zé)任編輯：于海琴）