王明珠,萬軍偉,白 通,劉 毅,沈 芳
(1.山東省地勘局第二水文地質(zhì)工程地質(zhì)大隊(duì)(山東省魯北地質(zhì)工程勘察院),德州 253072;2.山東省地?zé)崆鍧嵞茉刺綔y(cè)開發(fā)與回灌工程技術(shù)研究中心,德州 253072)
地面沉降是由于自然因素或人類活動(dòng)引發(fā)的松散地層固結(jié)、壓縮導(dǎo)致地面高程降低的地質(zhì)現(xiàn)象,是一種緩變性的地質(zhì)災(zāi)害[1]。已有研究[2-5]表明,華北平原地面沉降的主要原因是過量開采地下水。山東省德州市德城區(qū)位于華北平原,地面沉降嚴(yán)重[6]。為控制地面沉降的發(fā)展,德州市發(fā)布了深層承壓水禁采通知[7]。
地?zé)豳Y源開發(fā)利用對(duì)環(huán)境的負(fù)面影響較小,被稱為“清潔能源”,廣泛應(yīng)用于洗浴、理療、供暖和室內(nèi)種植等領(lǐng)域,具有較高的經(jīng)濟(jì)和商業(yè)利用價(jià)值[8-10]。山東德城區(qū)砂巖熱儲(chǔ)地?zé)豳Y源開發(fā)利用程度較高,目前,地?zé)崴_采是否對(duì)地面沉降造成影響尚未開展相關(guān)研究,在缺乏相關(guān)技術(shù)支持的情況下將直接影響地方政府對(duì)地?zé)豳Y源開發(fā)的決策,也將影響地?zé)豳Y源主管部門對(duì)地?zé)崴_發(fā)利用的規(guī)范管理。地?zé)崴畾w為深層承壓水而禁采將阻礙地?zé)崆鍧嵸Y源的開發(fā)利用。本文從德城區(qū)分層標(biāo)監(jiān)測(cè)數(shù)據(jù)、土體固結(jié)理論和現(xiàn)有回灌技術(shù)論述砂巖熱儲(chǔ)地?zé)豳Y源開采對(duì)地面沉降的影響,為相關(guān)部門合理規(guī)劃地?zé)豳Y源利用提供依據(jù),為今后研究砂巖熱儲(chǔ)地?zé)崴_采對(duì)地面沉降的影響提供參考。
館陶組砂巖熱儲(chǔ)是德城區(qū)地?zé)豳Y源勘查及開采的目標(biāo)熱儲(chǔ)。館陶組在全區(qū)均有分布,頂板埋深1 026~ 1 195 m,底板埋深1 350~1 650 m,與下覆東營組呈不整合接觸。館陶組厚度一般為250~550 m,在研究區(qū)中南部市區(qū)—于官屯—黃河涯鎮(zhèn)一帶較厚(一般>500 m),向外圍逐漸變薄;西部及德州開發(fā)區(qū)館陶組厚度一般約400 m;東部減河斷裂以東由北向南館陶組厚度逐漸增加。根據(jù)物探電測(cè)井、鉆探巖屑錄井資料,館陶組熱儲(chǔ)與地層的厚度比為30%~40%,熱儲(chǔ)厚度一般為 160~200 m。
垂向上,館陶組熱儲(chǔ)呈上細(xì)下粗的正旋回沉積,底礫巖明顯。根據(jù)沉積環(huán)境與巖性特征,館陶組熱儲(chǔ)可分為上、下兩段:上段熱儲(chǔ)巖性較細(xì),一般為淺灰色、灰白色粉細(xì)砂巖或細(xì)砂巖,單層厚度小,與泥巖互層,熱儲(chǔ)與地層厚度比約為30%;下段熱儲(chǔ)巖性較粗,主要為細(xì)-中砂巖、灰白色含礫粗砂巖和砂礫巖,單層厚度大,與泥巖相間,熱儲(chǔ)與地層厚度比約為40%。礫石成分以石英和黑色燧石為主,直徑1~10 mm,磨圓度中等,砂礫巖成巖性差,呈疏松狀,孔隙度一般為24%~ 30%。
館陶組熱儲(chǔ)層被多層泥巖覆蓋,構(gòu)成隔水頂板。地?zé)崴c上覆中層咸水、淺層淡水及大氣降水無水力聯(lián)系。地?zé)豳Y源屬于可再生資源,前提是在其開發(fā)利用時(shí)進(jìn)行地?zé)嵛菜毓?取熱不取水。該區(qū)地?zé)崴倭砍练e物沉積時(shí)保留下來的沉積水和封存水外,絕大部分地?zé)崴疄槌练e物形成后在漫長(zhǎng)的地質(zhì)時(shí)期,由遠(yuǎn)、近山區(qū)側(cè)向徑流補(bǔ)給。地?zé)崴惶嫖⑷?水平徑流極其滯緩,補(bǔ)給條件差,再生能力弱,與上覆含水層之間無水力聯(lián)系,基本為封閉系統(tǒng)。因此,地?zé)崴倪\(yùn)動(dòng)特征不受現(xiàn)代水文、氣象及地形地貌的影響,而受深層水文地球化學(xué)環(huán)境條件影響較大。人工開采是重要的排泄方式,此外還有在水力梯度下的水平運(yùn)移。開采條件下,開采量主要來源于含水層的儲(chǔ)變量[11]。
地?zé)崴畡?dòng)態(tài)主要受人工開采動(dòng)態(tài)影響,在冬季采暖期水位下降,非采暖期水位略回升,但水位整體呈下降趨勢(shì)。開采量主要集中在冬季供暖期,用于洗浴的水量較小,對(duì)水位恢復(fù)無明顯影響。降雨量主要集中在每年7—8月,該時(shí)間段為地?zé)崴换謴?fù)期,但并不能說明降雨量能影響地?zé)峋乃蛔兓?/p>
1997年3月,德熱1井成井時(shí)地?zé)崴裆罡叱龅孛?.33 m。2007年后,德城區(qū)地?zé)峋杆僭龆?地?zé)崴_采量急劇增加,至2018年10月產(chǎn)能測(cè)試的地?zé)崴宦裆顬?1.30 m,19年水位下降達(dá)59.63 m,平均每年下降3.14 m。2007年供暖期前,德城區(qū)地?zé)峋當(dāng)?shù)量較少,地?zé)衢_采量低,水位動(dòng)態(tài)變化較穩(wěn)定。1998年12月,地?zé)崴畹退宦裆顬?6.15 m,2005年11月供暖期之前地?zé)崴畹退宦裆顬?4.64 m。2006年3月,地?zé)崴畹退宦裆顬?3.90 m,開采量總體呈上升趨勢(shì)。自2013年后,受地方政策及回灌影響,地?zé)崴_采量逐漸降低(圖1)。截止2019年3月,該區(qū)地?zé)崴淖畹退宦裆顬?6.16 m,平均年水位下降率為5.56 m/a??傊?德城區(qū)地?zé)崴拈_采量總體上升,開采時(shí)間主要集中在供暖季。因德州市城區(qū)有少量洗浴用水,供暖季結(jié)束后,地?zé)崴杂休^小量的開采,但對(duì)水位恢復(fù)基本無影響[12]。
圖1 1998—2017年德熱1井開采量、水位動(dòng)態(tài)變化圖Fig. 1 Dynamic variation diagram of exploitationquantity and water level of Dere well 1 from 1998 to 2017
1997年3月,華北石油康海實(shí)業(yè)公司水井工程大隊(duì)在山東省魯北地質(zhì)工程勘察院施工完成了“德熱1井”,該井為德州第一眼地?zé)峋?999年12月,山東省地勘局第二水文地質(zhì)工程地質(zhì)大隊(duì)編寫了“山東省德州市城區(qū)地?zé)豳Y源詳查報(bào)告”[13],在德州市經(jīng)濟(jì)開發(fā)區(qū)施工完成了第二眼同深度的地?zé)峋?開展了地?zé)豳Y源綜合利用研究,協(xié)助建立了凱元溫泉度假村的地?zé)峁┡?、洗浴、室?nèi)外游泳池,開創(chuàng)了該區(qū)地?zé)豳Y源綜合利用的先河[14]。德城區(qū)主要開采館陶組熱儲(chǔ),地?zé)崴疁?0~70 ℃,地?zé)豳Y源開發(fā)主要用于供暖,部分用于洗浴、理療等。2019年,德城區(qū)共有地?zé)峋?3眼,多分布于運(yùn)河以東居民小區(qū),地?zé)峋觊_采總量為950萬m3/a,近些年來德城區(qū)地?zé)崴_采量統(tǒng)計(jì)見表1。
表1 德城區(qū)地?zé)崴_采量統(tǒng)計(jì)結(jié)果
根據(jù)德城區(qū)地?zé)峋删Y料,該區(qū)地?zé)峋诮瞪?0 m以內(nèi)單井涌水量為80~ 120 m3/h,水溫為54~58 ℃,水化學(xué)類型為CI-Na型水,礦化度為4 000~5 000 mg/L,水中含有對(duì)人體健康有益的微量元素。
從構(gòu)造因素看,德城區(qū)屬于華北地臺(tái),遼冀臺(tái)向斜和魯西臺(tái)背斜,是新構(gòu)造運(yùn)動(dòng)相對(duì)穩(wěn)定的區(qū)域。該區(qū)由構(gòu)造運(yùn)動(dòng)引起的地面下沉速率為2~3 mm/a,因此,區(qū)域構(gòu)造及新構(gòu)造運(yùn)動(dòng)不是地面沉降的主要因素。
從工程角度看,德城區(qū)一般工程為建筑持力層與下臥層1.5 m以下的地層或高層建筑持力層與下臥層30 m以下的地層。區(qū)內(nèi)場(chǎng)地屬中軟土、Ⅲ類一般建筑場(chǎng)地。不考慮基坑降水對(duì)地面沉降的影響,區(qū)內(nèi)工程建筑物荷載引起的地面累計(jì)沉降量為2~3 mm/a,說明該區(qū)工程因素不是地面沉降的主要原因。
從開采條件看,該區(qū)自開采深層地下水以來,地下水水位埋深不斷下降,水頭不斷降低,深層地下水一直處于超采狀態(tài)。由深層地下水降落漏斗圖(2017年)與地面沉降量等值線圖(1991—2017年)(圖2)可知,深層地下水降落漏斗中心與地面沉降中心位置一致,均處于國棉廠一帶,且漏斗擴(kuò)展范圍和形狀基本相似。
圖2 地下水降落漏斗與地面沉降分布圖Fig. 2 Distribution of groundwater funnel and land subsidence
根據(jù)沉降原因,選取沉降中心(德城區(qū)國棉廠)(圖3)和沉降邊緣(平原縣城附近)(圖4)進(jìn)行開采量、水位和累計(jì)沉降量分析。
圖3 沉降中心開采量、水位及累計(jì)沉降量變化趨勢(shì)圖Fig. 3 Variation trend of exploitation, water level and settling volume of subsiding center
圖4 沉降邊緣開采量、水位及累計(jì)沉降量變化趨勢(shì)圖Fig. 4 Variation trend of exploitation, water level and settling volume of subsiding edge
多年來,由于地下水開采,深層地下水水位不斷下降。2017年“德州市中心城區(qū)地下水壓采工作方案”正式實(shí)施,預(yù)計(jì)至2025年,深層地下水開采井計(jì)劃全部完成封填、封存,城區(qū)內(nèi)水位呈現(xiàn)波動(dòng)回升狀態(tài)。然而,地面沉降是不可逆的地質(zhì)災(zāi)害,隨著地下水開采的不斷進(jìn)行,累計(jì)沉降量不斷增加,但沉降中心區(qū)沉降速率有減小的趨勢(shì)。
德州市地面沉降發(fā)生于20世紀(jì)80年代,包括德城區(qū)、武城、平原、陵縣、寧津和河北省吳橋、景縣、故城、諸縣部分地區(qū),與周邊滄州、衡水地面沉降區(qū)連為一片,齊河縣城也存在沉降現(xiàn)象[15-16]。
結(jié)合二等水準(zhǔn)測(cè)量、GPS監(jiān)測(cè)、InSAR監(jiān)測(cè)資料,1991—2017年,德城區(qū)多年累計(jì)沉降量最大達(dá)1 300 mm,主要分布在運(yùn)河街道—廣川街道—天衢街道,面積約18 km2(圖5)。廣川街道—宋官屯街道一線以西,多年累計(jì)沉降量>1 000 mm,面積約104 km2。新華街道—袁橋鎮(zhèn)一線以西,多年累計(jì)沉降量>800 mm,面積約81 km2。德城區(qū)其他區(qū)域多年累計(jì)沉降量為600~800 mm[17]。
圖5 1991—2017年德城區(qū)累計(jì)沉降量示意圖Fig. 5 Schematic diagram of accumulated settling volume in Decheng district from 1991 to 2017
2013年,山東省地勘局第二水文地質(zhì)工程地質(zhì)大隊(duì)在新駐址建立了1組地面沉降監(jiān)測(cè)分層標(biāo)[18],埋設(shè)深度為500 m、300 m、60 m、2 m的4個(gè)分層標(biāo)孔,該組分層標(biāo)自2015年開始按照一月一次的頻率進(jìn)行相對(duì)沉降監(jiān)測(cè)。
因該組分層標(biāo)周圍附近地區(qū)沒有基巖標(biāo)輔助測(cè)量,僅能對(duì)地面沉降進(jìn)行相對(duì)沉降測(cè)量[19-22],測(cè)出60 m以淺、60~300 m、300~500 m地層的沉降數(shù)據(jù)。2015年1月—2019年10月,60 m以淺地層累計(jì)沉降量為6 mm,沉降速率為1.3 mm/a;60~300 m地層累計(jì)沉降量為92 mm,沉降速率為19.4 mm/a;300~500 m地層累計(jì)沉降量為77 mm,沉降速率為16.2 mm/a;500m以淺地層累計(jì)沉降量為175 mm(圖6)。
圖6 分層標(biāo)各層累計(jì)沉降量示意圖Fig. 6 Schematic diagram of accumulated settling volume of each layer
山東省地質(zhì)環(huán)境監(jiān)測(cè)總站分別于2014年10月和2017年11月對(duì)德州市進(jìn)行了水準(zhǔn)測(cè)量。根據(jù)測(cè)量結(jié)果,該組分層標(biāo)地面標(biāo)在此期間累計(jì)沉降量為154.5 mm,年平均沉降速率為51.5 mm/a。德城區(qū)因構(gòu)造因素產(chǎn)生的垂直變形沉降速率為2~3 mm/a,分層標(biāo)測(cè)量結(jié)果顯示500 m以淺地層年均沉降速率為37.9 mm/a。因此,德城區(qū)500 m以深地層因抽取地下水引起的地面沉降度率為分層標(biāo)地面標(biāo)測(cè)量沉降量減去500 m以淺地層沉降速率和因構(gòu)造因素引起的垂直變形速率[23-24],可知500 m以深地層沉降速率為11.1 mm/a。
根據(jù)德城區(qū)地下水開采井分布狀況統(tǒng)計(jì)結(jié)果[25],德城區(qū)深層地下水的主要開采層位為300~500 m含水層,開采井總數(shù)為170眼,占總井?dāng)?shù)的65%,實(shí)際開采總量為2 231.2萬m3/a。500~800 m含水層,開采井總數(shù)為93眼,占總井?dāng)?shù)的35%,實(shí)際開采量為951.1萬m3/a,占總開采量的30%(表2)。
表2 深層地下水開發(fā)利用現(xiàn)狀統(tǒng)計(jì)結(jié)果
根據(jù)地層資料,300~500 m及500~800 m壓縮層的天然孔隙比分別為0.25 和0.10。500~800 m層段開采量是300~500 m層段開采量的54%,500~800 m層段天然孔隙比是300~500 m層段的40%。由于現(xiàn)實(shí)壓縮因素涉及開采量、壓縮層厚度、地層巖性組成,本文分別將300~500 m和500~800 m概化為2個(gè)天然孔隙比不同的相同介質(zhì)地層,故推斷500~800 m層段土體壓縮量約為8.2 mm/a。
根據(jù)水準(zhǔn)及分層標(biāo)監(jiān)測(cè)數(shù)據(jù),500 m以深地層沉降速率為11.1 mm/a,其中500~800 m層段土體壓縮量約為8.2 mm/a,800 m以深地層壓縮量約為2.9 mm/a。
相對(duì)于51.5 mm/a的沉降速率,800 m以深地層(2.9 mm/a)的壓縮量,與構(gòu)造引發(fā)的地層壓縮量基本相同,故地?zé)衢_采對(duì)地面沉降的影響不明顯。
以德城區(qū)為研究對(duì)象,地層埋深為200~800 m,基本表現(xiàn)為超固結(jié)或微超固結(jié)狀態(tài)。目前,已知地面沉降主要是超量開采中深層地下水引起的。
將800 m以淺地層和1 000 m以深熱儲(chǔ)層視為相同固結(jié)狀態(tài),其孔隙度視為一致,即地層飽和重度視為一致。對(duì)比相同水位降幅下增加的土層有效應(yīng)力占各自地層自重應(yīng)力的比值,判斷地?zé)衢_采是否會(huì)引發(fā)地面沉降。
3.2.1 有效應(yīng)力原理
假設(shè)在抽水過程中土層總應(yīng)力不變,孔隙水壓力的減少必然導(dǎo)致土中有效應(yīng)力的等量增加,引起土層壓縮。潛水含水層、弱透水層及承壓含水層的孔隙水壓力分布如圖7所示。
圖7 含水層孔隙水壓力分布圖Fig. 7 Distribion of pore water pressure of aquifer
當(dāng)承壓含水層抽水后,產(chǎn)生Δh的水位降,則孔隙水壓力分布為圖中的AD′C′線,即在抽水前后土層中的孔隙水壓力產(chǎn)生了差值,計(jì)算公式為
Δuw=γw·Δh,
(1)
式中:Δuw為孔隙水壓力增量,kPa;γw為水的重度,kN/m3。
根據(jù)太沙基有效應(yīng)力原理,在地表以下一定深度地層的總應(yīng)力σ為該處有效應(yīng)力σ′與孔隙水壓力uw之和,即
σ=σ′+uw,
(2)
式中:σ為總應(yīng)力,kPa;σ′為有效應(yīng)力,kPa;uw為孔隙水壓力,kPa。
在抽水過程中孔隙水壓力下,總壓力基本保持不變。抽水后,隨著孔隙水壓力消散,孔隙水壓力Δuw的下降值如轉(zhuǎn)化為有效應(yīng)力的增量Δσ′,即獲得以下公式
σ=(σ′+Δσ′)+(u-Δuw)
Δσ′=Δuw,
(3)
3.2.2 地層應(yīng)力狀態(tài)
地層中的應(yīng)力可分為自重應(yīng)力、構(gòu)造應(yīng)力和來自地表的附加應(yīng)力[26]。自重應(yīng)力是由地層重量引起的,是地應(yīng)力的主要存在形式。
σz=γsat·z,
(4)
式中:σz為自重應(yīng)力,kPa;γsat為地層土的飽和容重,kN/m3;z為垂向深度,m。
3.2.3a值計(jì)算
將因水位降低引起的有效應(yīng)力增量Δσ′與該處地層自重應(yīng)力σz的比值定義為a值(無特征意義,方便論述),該值可間接反應(yīng)水位降幅對(duì)地層壓縮性的影響,a值越大,有效應(yīng)力變幅占地層自重應(yīng)力比值越大,則水位降幅對(duì)應(yīng)壓縮量越大,反之亦然,公式為
a=Δσ′/σz。
(5)
(1)深層地下水開采層a值計(jì)算。該層巖性主要為粉土、粉質(zhì)黏土、黏土和含水砂層。各類土的飽和重度y深約18 kN/m3,德城區(qū)深層地下水開采深度多為300~800 m,根據(jù)公式(4)計(jì)算的深層地下水開采層平均自重應(yīng)力約9.9 MPa。
(2)熱儲(chǔ)層a值計(jì)算。該區(qū)地?zé)衢_采深度多為1 000 m以深,開采層段為1 000~1 600 m。1 000 m以上地層飽和重度取值18 kN/m3,館陶組以砂巖和泥巖為主,飽和重度y熱取值為26 kN/m3。根據(jù)公式(4)計(jì)算的地層總平均自重應(yīng)力約33.6 MPa。
假設(shè)水位降幅均為5 m/a,根據(jù)公式(1)計(jì)算獲得的有效應(yīng)力Δσ′為50 kPa;根據(jù)公式(5)計(jì)算獲得的熱儲(chǔ)層a深為0.15%,a熱為0.15%。
計(jì)算結(jié)果表明,在相同固結(jié)條件下,同等水位降幅,熱儲(chǔ)層有效應(yīng)力增量占該層自重應(yīng)力的比值遠(yuǎn)低于深層地下水開采層有效應(yīng)力增量占該層自重應(yīng)力的比值。
根據(jù)鉆孔資料,300~500 m及500~800 m地層為正常固結(jié),屬松散地層。熱儲(chǔ)層地層為超固結(jié),為泥質(zhì)膠結(jié),屬半成巖地層,可壓縮性遠(yuǎn)遠(yuǎn)低于深層地下水開采層,進(jìn)一步佐證地?zé)衢_采對(duì)地面沉降的影響不顯著。
地?zé)崴烊粭l件下補(bǔ)給微弱,長(zhǎng)期集中區(qū)域、集中時(shí)段開采將導(dǎo)致儲(chǔ)層壓力下降,水位持續(xù)降低,地?zé)峁┡菜疁囟容^高,長(zhǎng)期排放將造成資源浪費(fèi),對(duì)環(huán)境造成污染?;毓嗍墙鉀Q上述問題最有效的措施[27-29],不僅可以維持儲(chǔ)層壓力,減緩水位下降,避免直接排放產(chǎn)生的環(huán)境污染,而且可以通過冷水被熱儲(chǔ)層加熱汲取更多的熱能,從而保證地?zé)豳Y源的可持續(xù)開采。
2018年,山東省國土資源廳、山東省水利廳發(fā)布的“關(guān)于切實(shí)加強(qiáng)地?zé)豳Y源保護(hù)和開發(fā)利用管理的通知”要求開采孔隙熱儲(chǔ)型地?zé)豳Y源的回灌率不低于80%,形成了“同層回灌,取熱不取水,以灌定采、采灌均衡”的德州地?zé)衢_采模式。在清潔高效的地?zé)衢_采模式下,地?zé)豳Y源開采的實(shí)際開采量大大降低,削減了地?zé)衢_采對(duì)地面沉降潛在的影響[30-31]。
德城區(qū)熱儲(chǔ)開采及回灌層為砂巖熱儲(chǔ)。二十多年的地?zé)衢_采導(dǎo)致地?zé)崴怀掷m(xù)下降,水位埋深較大,具有良好的回灌空間,為地?zé)峄毓嗵峁┝嘶A(chǔ)條件。自2006年開始,該區(qū)進(jìn)行了地?zé)峄毓嗉夹g(shù)研究,對(duì)地?zé)峋删に?、合理井距、回灌工藝流程等開展專題研究,地?zé)崴畬?shí)現(xiàn)了梯級(jí)、循環(huán)、綠色及可持續(xù)利用模式。地?zé)峁┡捎脽岜脵C(jī)組,實(shí)現(xiàn)了“采灌均衡、取熱還水”,供暖尾水進(jìn)行完全回灌。2016—2019年,利用德州市砂巖熱儲(chǔ)回灌示范工程連續(xù)開展3個(gè)供暖季的生產(chǎn)性回灌,以供暖尾水為回灌水源,回灌溫度30℃,采用除砂、過濾、排氣等尾水凈化工藝,運(yùn)用同層自然無壓回灌模式,一采一灌,穩(wěn)定回灌量為50~70 m3/h,供暖尾水100%回灌,回灌效果良好。
(1)1991—2017年,德城區(qū)多年累積沉降量最大達(dá)到了1 300 mm,面積約18 km2。多年累積沉降量>1 000 mm的面積約10 418 km2,多年累積沉降量>800 mm的面積約81 km2,其他區(qū)域多年累計(jì)沉降量為600 ~800 mm。
(2)德城區(qū)地?zé)衢_采對(duì)地面沉降的影響不顯著。在相同固結(jié)及同等水位降幅條件下,熱儲(chǔ)層有效應(yīng)力增量占該層自重應(yīng)力的比值遠(yuǎn)低于深層地下水開采層有效應(yīng)力增量占該層自重應(yīng)力的比值。熱儲(chǔ)層泥質(zhì)膠結(jié),可壓縮性遠(yuǎn)遠(yuǎn)低于深層地下水開采層。