王 倩 韓京增

（中國(guó)建筑材料工業(yè)地質(zhì)勘查中心山東總隊(duì)）

礦山生產(chǎn)排放的廢渣或廢石被隨意堆積在山坡或溝谷，加之其透水性降低，一旦受到暴雨便會(huì)對(duì)地下水補(bǔ)給產(chǎn)生影響［1］。為了降低該影響，并且為礦區(qū)生態(tài)環(huán)境及其周邊地下水的修復(fù)治理提供參考依據(jù)，研究非金屬礦山開(kāi)采土壤擾動(dòng)重構(gòu)對(duì)地下水補(bǔ)給的影響。本文通過(guò)構(gòu)建土壤擾動(dòng)重構(gòu)對(duì)地下水補(bǔ)給影響數(shù)值模型，以此為基礎(chǔ)，計(jì)算水位演變情況，確定入滲系數(shù)和含水率與水位，以期能有效分析非金屬礦山開(kāi)采土壤擾動(dòng)重構(gòu)對(duì)地下水補(bǔ)給的影響，為非金屬礦山開(kāi)采提供參考。

1 土壤擾動(dòng)重構(gòu)對(duì)地下水補(bǔ)給影響因素分析

從土層厚度、土壤溫度、地表徑流量等角度出發(fā)，研究分析土壤擾動(dòng)重構(gòu)對(duì)地下水補(bǔ)給影響因素，建立一種非均質(zhì)性、非均值、非穩(wěn)定的三維地下水模型，由此構(gòu)建土壤擾動(dòng)重構(gòu)對(duì)地下水補(bǔ)給影響數(shù)值模型。通過(guò)對(duì)土壤影響因素和地下水補(bǔ)充關(guān)系的歸納，得出了兩者之間的聯(lián)系［2］。土壤擾動(dòng)重構(gòu)主要影響因素：

（1）土壤密度及孔隙度。土壤密度越高，則土壤滲透率就越低；土壤孔隙度越大，土壤便越疏松。由于土壤密度太低，巖溶水便會(huì)不斷地從地面滲出，從而會(huì)及時(shí)補(bǔ)充地下水。

（2）溫度。溫度對(duì)土壤的影響比較大，土壤溫度會(huì)隨著深度的增大而逐漸上升，進(jìn)而影響地下水位的變化［3］。

（3）地表徑流量。在土壤干擾重構(gòu)過(guò)程中，地表徑流量是影響地下水補(bǔ)給最大、最直接的因素。

結(jié)合上述3個(gè)因素，著重分析其對(duì)土壤干擾重構(gòu)特征及影響因素之間變化規(guī)律和主要影響因子間交互作用機(jī)制，以便土壤擾動(dòng)重構(gòu)下地下水補(bǔ)給數(shù)值模型的構(gòu)建。

2 模型建立及初始條件設(shè)置

研究的采礦區(qū)域以大氣降水為主，為山地基巖裂縫水的橫向徑流，其流向與地貌基本相同，整體呈自東北向西南方向的趨勢(shì)，并在最后排入下游湖裙帶。由于礦山多年的開(kāi)采，導(dǎo)致部分地下水向露天坑口涌出，該區(qū)地下水以蒸發(fā)和礦排為主；結(jié)構(gòu)主要為砂石、煤層等組成的多層結(jié)構(gòu)，土壤組合結(jié)構(gòu)為腐殖土、黏土、中砂。通過(guò)GMS 建模，選擇20 個(gè)鉆孔數(shù)據(jù)導(dǎo)入GMS，對(duì)模擬區(qū)進(jìn)行不規(guī)則三角剖分，由克里金插值法進(jìn)行插值計(jì)算生成地質(zhì)結(jié)構(gòu)實(shí)體，建立研究區(qū)三維結(jié)構(gòu)模型，如圖1所示。試驗(yàn)研究擬設(shè)置了5 m 地下水位埋深，采用有限差分法，建立的非金屬礦山開(kāi)采區(qū)域地下水流場(chǎng)模型（圖2），在水流場(chǎng)下研究非金屬礦山開(kāi)采土壤擾動(dòng)重構(gòu)對(duì)地下水補(bǔ)給影響。為了降低初始條件對(duì)模擬結(jié)果影響，采用數(shù)值模擬方法，將2019—2021 年非金屬礦山采場(chǎng)開(kāi)采的土壤擾動(dòng)重構(gòu)引入到模型中。其中，模擬的時(shí)間步長(zhǎng)是1 d，而空間步長(zhǎng)是3 cm。

在進(jìn)行地下水流動(dòng)形態(tài)建模時(shí)，可以將其視為一種非均質(zhì)性、非均值、非穩(wěn)定的三維地下水模式，則：

式中，ψ1、ψ2為第一流量和第二流量邊界；H為計(jì)算期間水位，m；t為時(shí)間，s；λ為滲透系數(shù)；z為透水層的垂向滲透的厚度，m；κ為導(dǎo)水系數(shù)。

則相應(yīng)的初始條件為

式中，e2為已知第二水頭邊界；x、y為笛卡爾坐標(biāo)。

則相應(yīng)的邊界條件為

式中，e1為已知第一水頭邊界；σ0（x，y，t）為側(cè)向補(bǔ)給流量，m3。

使用多源數(shù)據(jù)整合方法，基于區(qū)域多年氣溫、降水、蒸發(fā)量等數(shù)據(jù)，地下水補(bǔ)給的最大可能受降水以及礦山水文地質(zhì)條件構(gòu)建土壤擾動(dòng)重構(gòu)的影響模型，充分考慮地質(zhì)條件和水文地質(zhì)條件的影響對(duì)地下水補(bǔ)給特性影響，對(duì)模擬結(jié)果展開(kāi)詳細(xì)分析［4-6］。

3 模擬結(jié)果分析

通過(guò)水位演變情況，分析礦山開(kāi)采區(qū)域土壤擾動(dòng)重構(gòu)水分入滲系數(shù)、含水率與水位之間關(guān)系。本文研究礦山開(kāi)采區(qū)域的年平均降雨量400～550 mm，年最大降雨量可達(dá)到1 100 mm，其中65%～85%的降雨量分布在7—9 月，而排泄量在30～32 m3/min，因?yàn)楣┧烤哂兄芷谛宰兓瑢?dǎo)致地下水位有逐年的季節(jié)性變化。由于土壤擾動(dòng)重構(gòu)對(duì)地下水的補(bǔ)充具有逆轉(zhuǎn)風(fēng)險(xiǎn)，使得地下水的有效補(bǔ)給量顯著提高。為了進(jìn)一步對(duì)土壤擾動(dòng)重構(gòu)模擬方法加以改進(jìn)和完善，研究水分入滲系數(shù)、含水率與水位之間關(guān)系，結(jié)果如圖3所示。

由圖3 可知，隨著入滲系數(shù)和含水率的升高，水位呈現(xiàn)上升的情況，并且不同入滲系數(shù)和含水率所引起的水位升高是不同的，入滲系數(shù)和含水率與水位上升呈線性關(guān)系，因此，通過(guò)分析結(jié)果，可以獲得其相關(guān)方程：

式中，C1、C2為大氣降水、土壤滲水；η、ζ為入滲系數(shù)、含水系數(shù)。

在非金屬礦山中，礦體的穩(wěn)態(tài)含水率與入滲系數(shù)隨深度的關(guān)系見(jiàn)圖4。由圖4（a）可知，土壤擾動(dòng)重構(gòu)區(qū)域土壤深度為100 cm 以內(nèi)時(shí)，含水率隨土壤擾動(dòng)重構(gòu)深度增加呈增大趨勢(shì)，但變化幅度較小。但是當(dāng)土壤深度超過(guò)100 cm 時(shí)，含水率隨土壤擾動(dòng)重構(gòu)深度增加呈增大趨勢(shì)，且變化幅度較大，大于土壤未重構(gòu)區(qū)域，說(shuō)明現(xiàn)有土壤擾動(dòng)重構(gòu)方式對(duì)超過(guò)100 cm 土壤含水率具有一定促進(jìn)作用。由圖4（b）可知，隨著土壤深度的增加，降水入滲系數(shù)也隨之增大，對(duì)于土壤擾動(dòng)重構(gòu)區(qū)域最大入滲系數(shù)為0.10、土壤未重構(gòu)區(qū)域最大入滲系數(shù)為0.135，說(shuō)明現(xiàn)有土壤擾動(dòng)重構(gòu)方式縮減了降水對(duì)地下水有效補(bǔ)給?；诖?，分析非金屬礦山開(kāi)采區(qū)域地下水位，如表1所示。

由表1可知，地下水位總體呈現(xiàn)先上升后下降趨勢(shì)，在2020 年12 月份的水位低于2019 年12 月，2021年3月水位低于2020年3月，說(shuō)明土壤擾動(dòng)重構(gòu)后使得地下徑流受到了一定影響。

為了深入分析土壤擾動(dòng)重構(gòu)對(duì)地下水補(bǔ)給的影響，將研究區(qū)域現(xiàn)場(chǎng)與草原區(qū)域現(xiàn)場(chǎng)進(jìn)行對(duì)比，草原土壤結(jié)構(gòu)為腐殖土和黃土，研究入滲系數(shù)隨土壤深度的變化規(guī)律如圖5所示。

?

根據(jù)圖5 可知，不同的土壤結(jié)構(gòu)的入滲系數(shù)不同，并且隨著土壤深度的變化，入滲系數(shù)的變化程度不同，其中草原研究區(qū)域的入滲系數(shù)高于本文研究的非金屬礦山開(kāi)采土壤區(qū)域，并且其在50 cm 后，入滲系數(shù)呈大幅度上升趨勢(shì)，表明其對(duì)50 cm 后的水分具有一定的促進(jìn)作用。綜合分析可知，由于黏土的滲透性能較差，土壤水在黏土中向下運(yùn)移受阻，導(dǎo)致黏土上部含水率較高，即黏土的“隔水效應(yīng)”對(duì)土壤水分的運(yùn)移起主控作用［7-8］。

4 結(jié) 語(yǔ)

通過(guò)對(duì)非金屬礦山開(kāi)采土壤擾動(dòng)重構(gòu)對(duì)地下水補(bǔ)給影響要素進(jìn)行分析可知，由于礦山開(kāi)采區(qū)域土地的擾動(dòng)和重建，使其地下水文要素及環(huán)境發(fā)生了較大變化，從而導(dǎo)致其地下水在原有基礎(chǔ)上得到了一定程度補(bǔ)充，從而造成了礦山地面的變形。為了降低采礦擾動(dòng)對(duì)地下水的影響，必須采取進(jìn)一步措施，以有效地抑制因礦山開(kāi)采對(duì)地下水補(bǔ)充能力的負(fù)面影響。