段啟杉,曹振東,孟凡濤,宋小慶
(1.貴州地質(zhì)工程勘察院,貴州 貴陽550008;2.貴州省地質(zhì)礦產(chǎn)勘查開發(fā)局111地質(zhì)大隊(duì),貴州貴陽550008)
淺層地溫能是一種清潔,分布廣泛的新型環(huán)保能源,開發(fā)利用前景極為廣闊。合理開發(fā)利用淺層地溫能,可改善能源結(jié)構(gòu),達(dá)到節(jié)能減排,實(shí)現(xiàn)環(huán)境友好型社會有重要現(xiàn)實(shí)意義。但是在西南巖溶地區(qū),開發(fā)淺層地溫能要面臨巖溶塌陷地質(zhì)災(zāi)害。西南地區(qū)巖溶面積達(dá)76×104km2,涵蓋桂、滇、黔、川、渝、鄂、湘、粵等八省(市、區(qū))[1],巖溶塌陷是區(qū)內(nèi)的最重要的地質(zhì)災(zāi)害類型,具有隱伏性、突發(fā)性和群發(fā)性等特點(diǎn)。巖溶塌陷問題已經(jīng)成為西南巖溶區(qū)開發(fā)淺層地溫能的最大障礙。地源熱泵在抽水和回灌過程中地下水局部降深較大、局部回灌抬升以及兩者疊加而產(chǎn)生的復(fù)雜的地下水流場,在建設(shè)和運(yùn)行過程中可能會發(fā)生局部降深過大和水位波動,從而引發(fā)地面塌陷,因此如何預(yù)測和防范地面塌陷已成為淺層地溫能開發(fā)亟需解決的問題。
巖溶塌陷的產(chǎn)生一般需要三個(gè)條件,即溶蝕空間、上部蓋層和致塌作用力。國內(nèi)很多專家?guī)r溶塌陷形成機(jī)理[2]、成因機(jī)制及防治途徑[3]、覆蓋型巖溶塌陷[4]、水位波動誘發(fā)巖溶塌陷[5]等方面進(jìn)行了研究,但地下水地源熱泵開發(fā)利用巖溶塌陷預(yù)防措施研究較少。
碳酸鹽巖巖溶發(fā)育程度是決定塌陷發(fā)生的基礎(chǔ)條件,而巖溶洞隙的發(fā)育具有一定的規(guī)律性,且控制著巖溶塌陷的分布。巖溶洞隙主要發(fā)育在淺部,隨深度增加而減弱,給塌陷物質(zhì)提供必要的儲存空間和運(yùn)移場所。
巖溶塌陷的發(fā)生與上覆土層組分和厚度密切相關(guān),上覆土層以殘積或殘坡積的紅粘土、含碎石、礫石的粘土為主,含水率隨深度增加而增加,土巖交界處含水率一般較大,抗破壞能力較差;巖溶塌陷一般發(fā)生在土層厚度小于10m的區(qū)域,土層越薄則越易產(chǎn)生塌陷。土洞一旦形成,需要有一定厚度的土層滿足其應(yīng)力拱高度的需要,在自然條件下才能保證土洞頂板的穩(wěn)定[6],如式(1)所示。
式中H:破壞拱高度;h0:土洞的高度;b:土洞跨度之半;:土體內(nèi)摩擦角;a=1/tan。
地下水活動是誘發(fā)巖溶塌陷的主要動力因素,如水位變化、滲流對覆蓋土層有明顯的改造作用。地下水位變化改變土層的受力狀態(tài),加速土層崩解速,促進(jìn)巖溶塌陷的形成;地下水位的變化改變土層的塑形狀態(tài),致使土層抗剪強(qiáng)度有所改變。地下水的底蝕、側(cè)蝕作用加速土層塌陷;地下水的滲流侵蝕逐漸帶走細(xì)小顆粒,降低土層粘性,使土層產(chǎn)生滲透變形,甚至發(fā)生塌陷。
地下水地源熱泵工程是抽取和回灌地下水,以其為換熱介質(zhì),開發(fā)利用淺層地溫能。工程運(yùn)行中,可能使局部地下水水位上升或者下降,改變自然狀態(tài)下地下水流場,可能誘發(fā)產(chǎn)生巖溶塌陷。在工程的運(yùn)行過程中,誘發(fā)產(chǎn)生巖溶塌陷的主要因素有水位降深過大、水位波動頻繁和水位瞬間降速過大。
覆蓋型巖溶區(qū)人為導(dǎo)致水位降深過大,會產(chǎn)生真空吸蝕、增強(qiáng)地下水侵蝕和潛蝕作用[7],降低地下水浮托力等。覆蓋型巖溶區(qū)環(huán)境承載力,在自然條件下可能產(chǎn)生巖溶塌陷,若抽水使水位降低,則加強(qiáng)巖溶塌陷的誘發(fā)因素,改變原有平衡狀態(tài)。為預(yù)防誘發(fā)發(fā)生巖溶塌陷,需嚴(yán)格控制水位降深。
以貴陽某地下水地源熱泵工程為例,進(jìn)行水位降深引起巖溶塌陷分析。該工程地處貴陽峰叢溶蝕盆地北端,地形由北向南微傾斜,地勢較平坦,地面高程約1080 m。地層為三疊系下統(tǒng)安順組(T1a),巖性為白云巖、夾溶塌角礫巖,據(jù)勘探鉆孔揭露資料反映,區(qū)內(nèi)巖溶極發(fā)育,主要巖溶個(gè)體形態(tài)為溶槽、溶孔、溶隙及溶洞,發(fā)育深度17~98.5 m,規(guī)模0.4~3.1 m。第四系人工回填土厚度4.3~9.8 m。抽水井覆蓋層厚度6.0 m,覆蓋層紅粘土物理力學(xué)參數(shù)如下:密度1.77 g/cm3,內(nèi)聚力15 KP,抗拉強(qiáng)度4 KP,內(nèi)摩擦角 4°,土粒密度2.75 g/cm3,孔隙比 1.922,崩塌系數(shù) 0.0157m/a;水的密度1.0 g/cm3,重力加速度 9.8m/s2。
在地源熱泵工程建設(shè)期間,開展了抽水試驗(yàn),當(dāng)水位降深5.10米時(shí)地表發(fā)生了巖溶塌陷。根據(jù)對現(xiàn)場條件的分析,該處水文地質(zhì)條件符合負(fù)壓條件下的塌陷預(yù)測模型(式2)[8]。根據(jù)模型計(jì)算(工況條件:塌陷處無地面荷載 P=0,抽水試驗(yàn)期取t=0.019a),令K=1.0時(shí),由公式(2)計(jì)算上覆土層極限狀態(tài)下臨界致塌水位降深△h為4.7 m。抽水試驗(yàn)實(shí)際水位降深為5.1 m,大于臨界致塌水位降深4.7 m,因而抽水井附近出現(xiàn)巖溶地面塌陷。
水位波動誘發(fā)巖溶塌陷多發(fā)生于水位波動幅度較大、波動頻率較快、波動時(shí)間較長的情況。水位下降階段,水力坡度變大,水流速度加快,產(chǎn)生較大的動水壓力,可帶走松散覆蓋層細(xì)小顆粒,降低土層力學(xué)性質(zhì),增大發(fā)生巖溶塌陷的概率。
巖溶區(qū)覆蓋土層以殘積或殘坡積的紅粘土、含碎石、礫石的粘土為主,該類土層一般具有脹縮性,水位波動導(dǎo)致含水率不斷增加和減小,從而不斷改變土體結(jié)構(gòu),導(dǎo)致力學(xué)性質(zhì)變差,是產(chǎn)生巖溶塌陷的主要原因之一。
瞬間水位降速過大,會引起真空現(xiàn)象,加速土體氣體溢出和水體滲流,水位快速下降可形成巖溶空腔負(fù)壓;同時(shí)覆蓋型巖溶區(qū)含水層中水位快速下降,地下水水力坡度瞬間增大。滲透力大小與水力坡度成正比,因此水位降速過大導(dǎo)致滲透力增大。真空現(xiàn)象和滲透壓力增大過程加強(qiáng)了巖溶塌陷的誘發(fā)因素,增大巖溶塌陷的可能性。
瞬間水位上升過大,巖溶空腔中氣體形成高壓氣團(tuán),對周圍土體產(chǎn)生正壓力,可引起氣爆,直接導(dǎo)致巖溶塌陷的發(fā)生;水位急速上升過程中,也增大了滲透力,地下水的作用也減弱土層的抗剪能力,并且使土層自重增加,極易引起巖溶塌陷。
地下水地源熱泵的開發(fā)利用要進(jìn)行抽水和回灌,人為改變區(qū)域地下水的天然流場,并可能引起水位波動頻繁、瞬間水位降速過大均會引起覆蓋型巖溶區(qū)巖溶塌陷,因此,覆蓋型巖溶區(qū)地下水地源熱泵開發(fā)利用中,巖溶塌陷形成機(jī)理的研究顯得非常重要,對預(yù)防巖溶塌陷有著積極的作用。
查清巖溶發(fā)育特征,開發(fā)利用影響區(qū)域內(nèi),控制水位降深小于致塌降深。地下水地源熱泵的開發(fā)利用有別于其他地下水開采,可通過抽水井和回灌井的合理布置,避免地質(zhì)災(zāi)害承載力較差區(qū)域工程運(yùn)行水位與天然地下水位相差較大;避免抽水井和回灌井布置于巖溶強(qiáng)發(fā)育帶,避免產(chǎn)生巖溶塌陷。
地下水地源熱泵的開發(fā)利用過程初期,抽水量和回灌量逐漸增加,避免產(chǎn)生水位急速變話,減小巖溶塌陷的可能。若出現(xiàn)抽水井水泵損壞情況,則減少或者停止該抽水井影響范圍內(nèi)回灌井的回灌過程,避免水位變幅過大。
水位波動是引起巖溶塌陷的主要因素之一,地下水地源熱泵的開發(fā)利用是群孔抽水和群孔回灌的地下水利用過程,在地下水地源熱泵的運(yùn)行過程中,人為引起水位波動的可能性極大,合理控制抽水井和回灌井的運(yùn)行,減少人為水位波動,避免巖溶塌陷的發(fā)生。
(1)覆蓋型巖溶區(qū)地質(zhì)環(huán)境承載力較低,預(yù)防巖溶塌陷是地下水開發(fā)利用過程中的首要任務(wù),若無客觀、深入的巖溶塌陷機(jī)理巖溶,可能影響淺層地溫能的開發(fā),造成資源浪費(fèi)。
(2)合理布置抽水井和回灌井,減小人為改變地下水流場的程度,避免發(fā)生巖溶塌陷。
(3)建立地下水地源熱泵運(yùn)行監(jiān)測和預(yù)警系統(tǒng),避免水位下降或者上升過大,減小巖溶塌陷發(fā)生的可能性。
[1]袁丙華.中國西南巖溶石山地區(qū)地下水資源及生態(tài)環(huán)境地質(zhì)研究[M].成都:電子科技大學(xué)出版社,2007.
[2]徐衛(wèi)國,趙桂榮.論巖溶塌陷形成機(jī)理[J].煤炭學(xué)報(bào),1986,2:1-11.
[3]劉廣潤,程伯禹.巖溶塌陷的類型、成因機(jī)制及防治途徑-兼論武漢市巖溶塌陷勘查工作要點(diǎn)[J].工程地質(zhì)學(xué)報(bào),2011,9(4):414-417.
[4]劉秀敏,陳從新,沈強(qiáng)等.覆蓋型巖溶塌陷的空間預(yù)測與評價(jià)[J].巖土力學(xué),2011(9):2785-2790.
[5]蒙彥,殷坤龍,雷明堂.水位波動誘發(fā)巖溶塌陷的概率分析[J].中國巖溶,2006(9):239-242.
[6]康彥仁,項(xiàng)式均,陳健,等.中國南方巖溶塌陷[M].南寧:廣西科學(xué)技術(shù)出版社,1990.
[7]杜正民,吳光明,洪亮.潛蝕作用導(dǎo)致巖溶塌陷地質(zhì)災(zāi)害的實(shí)例分析[J].水文地質(zhì)工程地質(zhì),2007,34(3):89-92.
[8]譚鑒益.壇狀圓柱狀巖溶塌陷預(yù)測模型與分析[J].水文地質(zhì)工程地質(zhì),2001,28(4):19-22.