賀曉浪 段中會(huì) 苗霖田 張建軍

摘 要：針對(duì)沙漠區(qū)潛水含水層富水性評(píng)價(jià)問(wèn)題，基于野外調(diào)查和理論研究，總結(jié)認(rèn)為毛烏素沙漠區(qū)內(nèi)潛水含水層富水性主要由3層次5個(gè)因素綜合控制，結(jié)合層次分析法、模糊聚類，提出了該區(qū)內(nèi)潛水含水層富水性評(píng)價(jià)模型，并分析了各影響因素在一個(gè)水文年內(nèi)的變化特征，綜合評(píng)價(jià)了區(qū)內(nèi)薩拉烏蘇組潛水含水層的富水性分區(qū)特征及其動(dòng)態(tài)變化規(guī)律。結(jié)果表明：區(qū)內(nèi)潛水含水層以富水性弱區(qū)分布最廣，且明顯大于其他區(qū)域，同時(shí)該區(qū)的面積、形態(tài)也最容易受到枯水期和豐水期變化的影響。而較弱富水區(qū)的面積對(duì)季節(jié)變化不敏感，但其分布的區(qū)域發(fā)生了較大變化。相比之下，較強(qiáng)和強(qiáng)富水含水區(qū)則受到的影響較小。

關(guān)鍵詞：地質(zhì)資源與地質(zhì)工程;潛水含水層;富水性評(píng)價(jià);保水采煤

中圖分類號(hào)：TD 742.3?文獻(xiàn)標(biāo)志碼：ADOI：10.13800/j.cnki.xakjdxxb.2019.0113文章編號(hào)：1672-9315（2019）01-0088-08

Unconfined aquifer water?richness and dynamic change

features of Mu Us Desert Area in northern Shaanxi

HE Xiao?lang，DUAN Zhong?hui3，MIAO Lin?tian1，3，ZHANG Jian?jun3

（1.College of Geology and Environment，Xi’an University of Science and Technology，Xi’an 710054，China;

2.China Coal Energy Research Institute Co.，Ltd.，Xi’an 710054，China;

3.Key Laboratory of Coal Resources Exploration and Comprehensive Utilization，Ministry of Land and Resources，Xi’an 710021，China）

Abstract：Based on the field investigation and theoretical research，it is summarized that the water richness of phreatic aquifer in Mu Us desert area is mainly controlled by 5 factors which can be formed of 3 hierarchies，AHP（analytic hierarchy process）and fuzzy clustering method were adopted to build the water richness evaluation model of unconfined aquifer in the district，and the characteristics of each index in one hydrological year were analyzed.Partition characteristics and the dynamic change features ofSalawusu phreatic aquifer were comprehensively evaluated.The results show that the area with weakest water yield is most widely distributed and significantly greater than other areas，its area and shape are also the most vulnerable to the wet and dry season change，while the area of weaker water?richness district would be less influenced by this

change，but significant changes have taken place in its distribution shape.Stronger and strongest water richness zones，in contrast，were less affected.

Key words：geological resources and geological engineering;unconfined aquifer;water?richness evaluation;coal mining with water protection

0?引?言

陜北榆神府礦區(qū)屬毛烏素沙漠區(qū)，薩拉烏蘇組潛水含水層分布范圍廣闊，其存在不僅為煤炭、能源化工基地提供了直接水源，也是陜北毛烏素沙漠地區(qū)水泉的直接補(bǔ)給來(lái)源和維系區(qū)內(nèi)生態(tài)環(huán)境的重要物質(zhì)基礎(chǔ)[1-3]，是區(qū)內(nèi)最重要的含水層，因而分析研究區(qū)內(nèi)潛水的富集規(guī)律，進(jìn)行潛水含水層的富水性評(píng)價(jià)，為該含水層的合理開(kāi)發(fā)利用及保護(hù)提供依據(jù)則顯得尤為重要。

當(dāng)前國(guó)內(nèi)外富水性研究主要包含現(xiàn)場(chǎng)抽（放）水試驗(yàn)[4]、物探定性研究[5-6]、武強(qiáng)、邱士利等人研究的[7-9]多因素綜合分析預(yù)測(cè)法等。但多從礦井防治水角度來(lái)研究含水層的富水性，進(jìn)而使得研究目標(biāo)主要集中在靠近煤層的承壓含水層，而對(duì)于潛水含水層的研究較少。近年來(lái)保水開(kāi)采理論及相關(guān)技術(shù)研究快速發(fā)展[10]，但研究的方向卻主要集中在討論煤水間距與“兩帶”發(fā)育高度間的位置關(guān)系及其控制方面，對(duì)于保水開(kāi)采的目標(biāo)含水層—潛水含水層的富水性研究卻并未隨之明顯增多和深入?？紤]到潛水含水層與承壓含水層在補(bǔ)徑排及賦存條件的巨大差異以及毛烏素沙漠區(qū)干旱少水的水文地質(zhì)特征，開(kāi)展該區(qū)域內(nèi)的富水性動(dòng)態(tài)綜合評(píng)價(jià)，為加強(qiáng)生態(tài)潛水含水層的保護(hù)、合理開(kāi)發(fā)以及區(qū)內(nèi)煤礦的水害防治提供指導(dǎo)等方面都具有重要的現(xiàn)實(shí)意義。

文中意在以小保當(dāng)井田為例，選取多個(gè)指標(biāo)并分析其變化特征，再結(jié)合層次分析、模糊聚類法構(gòu)建潛水含水層富水性綜合評(píng)價(jià)模型，綜合分析井田薩拉烏蘇組潛水含水層的富水性分區(qū)特征及其動(dòng)態(tài)變化規(guī)律。

1?井田水文地質(zhì)特征及概況

1.1?井田概況

小保當(dāng)井田處于陜西省神木縣境內(nèi)，位于榆神礦區(qū)腹地。井田地貌主要為第四系風(fēng)積半固定沙丘和固定沙丘，屬毛烏素沙漠丘陵地貌（圖1）。區(qū)內(nèi)降雨集中，6～9月份降水量占了全年的70%.全年蒸發(fā)量為1 668.0 mm，超過(guò)了降雨量的4倍。作為第四代現(xiàn)代化煤礦示范項(xiàng)目，小保當(dāng)井田地質(zhì)構(gòu)造簡(jiǎn)單，開(kāi)采條件好，其中一號(hào)煤礦年產(chǎn)1 500萬(wàn)t，二號(hào)煤礦年產(chǎn)1 300萬(wàn)t，屬特大型建設(shè)煤礦。但本地區(qū)地下水資源匱乏，生態(tài)地質(zhì)環(huán)境脆弱，對(duì)潛水的依賴性強(qiáng)。

1.2?潛水含水層水文地質(zhì)特征

作為區(qū)內(nèi)最重要的含水層，潛水含水層由薩拉烏蘇組與風(fēng)積沙構(gòu)成，其主要特征如下：①巖性特征。從平面分布上來(lái)看，全區(qū)以細(xì)沙為主，中沙次之，僅在個(gè)別地方有粉沙分布（圖2）;②水文地質(zhì)參數(shù)。該含水層出露3個(gè)下降泉，單泉流量6.60～6.92 L/s.據(jù)該層鉆孔抽水資料，統(tǒng)降統(tǒng)徑單位涌水量q=

0.011 04～0.4524 L/s·m，平均滲透系數(shù)K=0.094 67～4.037 9 m/d，礦化度M=0.271～0.467 g/L（表1）;③含水層厚度。風(fēng)積沙與薩拉烏蘇組累計(jì)厚度20.0～75.3 m，沙層幾乎覆蓋全區(qū)，表現(xiàn)為低洼區(qū)厚，梁峁區(qū)薄，即井田邊界附近厚，中部?。▓D3）。④潛水含水層水位埋深。潛水含水層水位埋深0～30.50 m，整體埋藏較淺，大部分在5 m以下（圖4），灘地區(qū)內(nèi)更甚。

潛水含水層下伏保德組相對(duì)隔水層，區(qū)內(nèi)有零星出露，一般結(jié)構(gòu)致密，為半堅(jiān)硬狀。在風(fēng)沙灘地附近薄，甚至缺失。該隔水層的頂面形態(tài)對(duì)上覆沙層的厚度的控制作用明顯：頂面標(biāo)高高則上覆沙層薄，反之則較厚。

1.3?潛水含水層補(bǔ)徑排條件

潛水含水層的補(bǔ)給源主要為大氣降雨，井田地表沙層廣泛分布，其入滲系數(shù)為0.30～0.60，平均滲透系數(shù)K=1.805 4 m/d，極有利于降雨的下滲與運(yùn)移?？紤]地形地貌和補(bǔ)徑排條件的差別，將井田以地面分水嶺為界，劃分為東、西2個(gè)水文單元（圖5）。分水嶺東部潛水總體向東北或東南方向徑流，部分潛水經(jīng)保德組紅土“天窗”補(bǔ)給下伏含水層，也有部分潛水在地勢(shì)較低處出露，形成海子或泉;分水嶺西部潛水總體向西南或者西北方向徑流，在低洼地帶匯聚。補(bǔ)給區(qū)與排泄區(qū)之間的地形高差大，因此兩地內(nèi)水力坡度有一定差距：補(bǔ)給區(qū)水力坡度較大，而排泄區(qū)水力坡度變小。區(qū)內(nèi)地表和地下水分水嶺的水文地質(zhì)邊界基本一致，因而潛水的補(bǔ)徑排條件分區(qū)特性明顯。

2?潛水含水層富水性綜合分區(qū)

2.1?評(píng)價(jià)指標(biāo)體系的建立

潛水含水層富水性的影響因素和評(píng)價(jià)指標(biāo)眾多，各個(gè)相關(guān)因素關(guān)系復(fù)雜。單個(gè)指標(biāo)因素難以準(zhǔn)確描述富水性這種受控于多因素的含水層水文性質(zhì)。筆者認(rèn)為毛烏素沙漠區(qū)內(nèi)潛水含水層的富水性指標(biāo)體系可由3層次，5個(gè)指標(biāo)構(gòu)成。包含有：?jiǎn)挝挥克?、滲透系數(shù)、補(bǔ)給模數(shù)和含水層厚度和巖性粒度[11-16]（圖6）。

其中單位涌水量、滲透系數(shù)和含水層厚度可以直接測(cè)定，而巖性粒度指數(shù)可以通過(guò)以下方式確定（補(bǔ)給模數(shù)的確定方法見(jiàn)后文）：依據(jù)《巖土工程勘察規(guī)范》（GB50021?2001）可將含水層的巖性粒度分成礫砂、粗砂等5個(gè)等級(jí)?？紤]到對(duì)于潛水來(lái)說(shuō)，粒度越大則越有利于潛水的存儲(chǔ)和滲流[17]，因而以礫砂為最高級(jí)賦予其粒度指數(shù)為9，其他依次類推見(jiàn)表2.

2.2?評(píng)價(jià)指標(biāo)權(quán)重的確定

采用層次分析法（Analytic Hierarchy Process，簡(jiǎn)稱AHP法）計(jì)算指標(biāo)權(quán)重值，計(jì)算各指標(biāo)參數(shù)對(duì)于總目標(biāo)的權(quán)重值見(jiàn)表5.

2.3?評(píng)價(jià)指標(biāo)的融合及評(píng)價(jià)模型的建立

將層次分析法所確定的權(quán)重系數(shù)與各個(gè)主控因素相結(jié)合，建立小保當(dāng)井田內(nèi)薩拉烏蘇組孔隙潛水含水層富水性評(píng)價(jià)的綜合指標(biāo)——富水性指數(shù)[18]。

D=nk=1Wk·fk（x，y）

=0.036 4f1（x，y）+0.579 8f2（x，y）+0.139 4f3（x，y）+0.172 2f4（x，y）+0.072 2f5（x，y）（1）

式中D為富水性指數(shù);Wk為影響因素的權(quán)重值;fk（x，y）為單因素影響值函數(shù);x，y為坐標(biāo);n為影響因素的個(gè)數(shù)。

為使不同評(píng)價(jià)參數(shù)具有可比性，需進(jìn)行數(shù)據(jù)量綱統(tǒng)一化處理（文中采用極大值法），以消除各個(gè)相關(guān)參數(shù)不同量綱數(shù)據(jù)對(duì)評(píng)價(jià)結(jié)果的影響，并使數(shù)據(jù)相對(duì)化[19]，并針對(duì)數(shù)據(jù)模型開(kāi)展了一致性性檢驗(yàn)[20]，滿足要求。

2.4?綜合評(píng)價(jià)指標(biāo)分級(jí)標(biāo)準(zhǔn)的確定

對(duì)小保當(dāng)井田內(nèi)薩拉烏蘇組孔隙潛水含水層富水性評(píng)價(jià)綜合指標(biāo)D進(jìn)行模糊聚類，從而確定分區(qū)閾值為0.143，0.310，0.561，根據(jù)分區(qū)閾值將研究區(qū)劃分為弱、較弱、較強(qiáng)和強(qiáng)4個(gè)不同等級(jí)的區(qū)域（表6）。

3?富水性指標(biāo)的變化分析

3.1?變化時(shí)間段的劃分

通過(guò)對(duì)區(qū)內(nèi)及周邊井田水文長(zhǎng)觀井、長(zhǎng)觀鉆孔的調(diào)查發(fā)現(xiàn)：潛水位的變化與區(qū)內(nèi)降雨量的變化情況相類似，亦有豐水期與枯水期之分，雖兩者變化略有差別，但時(shí)間差值較小。因而綜合降雨和井田潛水水位變化情況，將7～10月份定義為井田豐水期，一個(gè)水文內(nèi)11月份至次年6月份定義為枯水期。

3.2?豐水期與枯水期內(nèi)含水層水文參數(shù)

在評(píng)價(jià)指標(biāo)中，部分指標(biāo)在一個(gè)水文年內(nèi)發(fā)生的變化幾乎可以忽略不計(jì)，如巖性粒度和滲透系數(shù)（滲透系數(shù)大小主要受巖性顆粒影響。在一個(gè)水文年內(nèi)巖性特征基本不會(huì)發(fā)生變化的情況下，該參數(shù)也不會(huì)發(fā)生明顯變化）。也有部分指標(biāo)可能會(huì)發(fā)生較大的變化，如含水層厚度，補(bǔ)給模數(shù)和單位涌水量。

3.2.1?含水層厚度

潛水含水層的厚度是指潛水面到隔水底板間的垂直距離。據(jù)小保當(dāng)和相鄰小壕兔井田內(nèi)水文長(zhǎng)觀資料，兩處長(zhǎng)觀井水位變化的平均值均在10.6 cm左右，即得：7～10月的水位值相較于一個(gè)水文年內(nèi)得其他月份整體上漲10.6 cm，此外還可以降雨入滲量來(lái)驗(yàn)證該水位變化值（含水層厚度變化值）的準(zhǔn)確性。

3.2.2?含水層單位涌水量q

裘布衣公式推導(dǎo)出的單位涌水量的計(jì)算式形式復(fù)雜，涉及的參數(shù)眾多，利用該公式難以直接判定豐水期與枯水期變化對(duì)單位涌水量的影響。因此引入導(dǎo)水系數(shù)，將單位涌水量的計(jì)算簡(jiǎn)化以方便判斷。而據(jù)蘭太權(quán)[21-22]等人的研究，單位涌水量q與導(dǎo)水系數(shù)T在數(shù)值上存在線性關(guān)系，即

q=aT（2）

式中?q為單位涌水量，m3/（d·m）;a為比例系數(shù);T為導(dǎo)水系數(shù)，m3/（d·m）;而導(dǎo)水系數(shù)一般用來(lái)表征含水層導(dǎo)水能力的大小，在數(shù)值上等于滲透系數(shù)（K）與含水層厚度（H）的乘積，即T=KH.為比較出一個(gè)水文年內(nèi)7～10月份潛水含水層的單位涌水量q′和其余月份的單位涌水量q的變化情況，文中采用比值法（作商法）進(jìn)行分析對(duì)比。

q′q

=aK′H′aKH（3）

式中?H′，H為同一水文年內(nèi)7～10月份和其他月份的含水層厚度，m;K′，K分別為同一水文年內(nèi)7～10月和其他月份的含水層滲透系數(shù)，m/d.已知滲透系數(shù)在該時(shí)間段內(nèi)基本不發(fā)生變化，即K=K′. 而H′=H+0.106.將其帶入式（3）化簡(jiǎn)可得

q′q

=1+0.106H（4）

區(qū)內(nèi)潛水含水層厚度平均約21 m.此時(shí)q′與q兩者間的比值約為1.005.考慮到潛水含水層的實(shí)測(cè)單位涌水量q=0.059 6～1.971 0 L/s·m，因而在一個(gè)水文年內(nèi)，單位涌水量的變化非常有限。文中研究時(shí)將單位涌水量以不變處理。

3.2.3?含水層補(bǔ)給模數(shù)

在以降雨為主要補(bǔ)給來(lái)源的情況下，本次研究將降雨的入滲補(bǔ)給量作為補(bǔ)給模數(shù)[23-25]式（5），即

M=αPF（5）

式中?M為補(bǔ)給模數(shù)，m3;α為降雨入滲系數(shù)，無(wú)量綱;P為降雨量，mm;F為補(bǔ)給區(qū)面積，km2.

考慮到區(qū)內(nèi)沙層厚度變化大，包氣帶條件不同，而分別計(jì)算分水嶺兩側(cè)的補(bǔ)給模數(shù)值[26]：分水嶺西側(cè)的豐水期補(bǔ)給模數(shù)：分水嶺西側(cè)井田面積F1=78 232 448.90 m2，入滲補(bǔ)給系數(shù)α1=

0.535，7～10月份的豐水期內(nèi)的降雨量為284.605 mm，將以上數(shù)據(jù)帶入式（5），得出：豐水嶺西側(cè)在豐水期的總補(bǔ)給模數(shù)M總=1 191.20 m3，月補(bǔ)給模數(shù)為M=297.80×104 m3/mo.分水嶺西側(cè)的枯水期補(bǔ)給量和分水嶺東側(cè)在不同時(shí)期的補(bǔ)給模數(shù)計(jì)算方法同上，計(jì)算結(jié)果見(jiàn)表7.

4?富水性動(dòng)態(tài)綜合評(píng)價(jià)

按照上文毛烏素沙漠區(qū)潛水含水層富水性指數(shù)的構(gòu)建方法，將小保當(dāng)井田豐水期和枯水期內(nèi)潛水含水層的各水文地質(zhì)參數(shù)代入模型運(yùn)算，即可得到豐水期、枯水期和全年期內(nèi)小保當(dāng)井田薩拉烏蘇組孔隙潛水含水層的富水性指數(shù)，在此基礎(chǔ)之上利用GIS強(qiáng)大的空間處理能力按照分區(qū)標(biāo)準(zhǔn)繪制分區(qū)評(píng)價(jià)圖（圖8）。

4.1?不同時(shí)期內(nèi)潛水含水層富水性分區(qū)特征

整體來(lái)看，井田內(nèi)以弱富水區(qū)分布最廣，其次為強(qiáng)富水區(qū)，較弱和較強(qiáng)區(qū)介于兩者之間。各區(qū)的面積、形態(tài)在不同時(shí)期各不相同。

分布面積方面，對(duì)季節(jié)變化最敏感的區(qū)域?yàn)槿醺凰畢^(qū)，由豐水期到枯水期沿分水嶺兩側(cè)向外擴(kuò)張，面積增加了井田總面積的19.0%，其次為較強(qiáng)富水區(qū)，減少14.3%，;再為強(qiáng)富水區(qū)，減少5.3%;而對(duì)于較弱富水區(qū)的影響很小，僅僅增加了0.6%（圖9）。

分布形態(tài)方面，較強(qiáng)富水區(qū)的整體形態(tài)變化不大，主要為在一號(hào)井田內(nèi)向東部收縮;強(qiáng)富水區(qū)在井田的北部和東部區(qū)域有不同程度的收縮，但變化也都不大;對(duì)于較弱富水區(qū)域來(lái)說(shuō)，雖面積變化小，但其分布的區(qū)域卻發(fā)生了較大的變化：分布區(qū)域更分散，由原來(lái)相互連接的弓形區(qū)域變?yōu)楸灰?、二?hào)井田分界線所截?cái)嗟拈L(zhǎng)條帶狀分布區(qū)。

綜上所述，小保當(dāng)井田潛水含水層從豐水期到枯水期，面積和形態(tài)變化最容易受到影響的是弱富水區(qū)，究其原因主要在于毛烏素沙漠區(qū)內(nèi)潛水含水層的滲透系大，弱富水區(qū)內(nèi)的含水層厚度、補(bǔ)給模數(shù)等指標(biāo)，更易受到豐、枯水季變化的影響。而較弱富水區(qū)雖然整體面積變化不大，但是其分布的范圍發(fā)生了較大變化。造成這種現(xiàn)象的主要原因在于較弱富水區(qū)屬于弱-強(qiáng)富水區(qū)的過(guò)渡區(qū)域，隨著豐、枯水季的轉(zhuǎn)換，區(qū)內(nèi)含水層厚度和補(bǔ)給模數(shù)也迅速變化，導(dǎo)致了較弱富水區(qū)的分布區(qū)域發(fā)生了較大變化，而受到分區(qū)閾值的控制，較弱富水區(qū)主要與弱富水區(qū)進(jìn)行轉(zhuǎn)化，這也導(dǎo)致了較弱富水區(qū)的面積變化較小[27-29]。

4.2?潛水含水層富水性綜合分區(qū)特征分析

通過(guò)對(duì)比井田富水性的分區(qū)特征和分水嶺的展布形態(tài)，不難發(fā)現(xiàn)井田內(nèi)分水嶺的展布形態(tài)對(duì)富水性分區(qū)具有明顯的控制作用：分水嶺附近區(qū)域的含水層富水性要明顯弱于遠(yuǎn)離分水嶺分布的區(qū)域。而由于分水嶺自身亦非處在一個(gè)平面上，即分水嶺本身也存在較大的高差（尤其是保德組頂面高差），這也決定了沿著分水嶺會(huì)出現(xiàn)不同富水性區(qū)域。

保德組的頂板起伏，決定了分水嶺的位置以及分水嶺本身的起伏特征，這也就進(jìn)一步使得弱富水性區(qū)域在在接受降雨補(bǔ)給后，主要通過(guò)地下徑流的方式補(bǔ)給了水位標(biāo)高相對(duì)較低的富水區(qū)域，甚至在這些區(qū)域形成小型湖泊或海子，進(jìn)而使得弱富水區(qū)域與相對(duì)強(qiáng)富水區(qū)域內(nèi)的富水性差距進(jìn)一步擴(kuò)大。

5?結(jié)?論

1）通過(guò)對(duì)毛烏素沙漠區(qū)水文地質(zhì)條件分析，得出潛水含水層的富水性評(píng)價(jià)體系可有由3層次，5個(gè)指標(biāo)構(gòu)建，并采用層次分析、模糊聚類等方法來(lái)綜合評(píng)價(jià)該區(qū)內(nèi)潛水含水層的富水性;

2）一個(gè)水文年內(nèi)，潛水含水層的滲透系數(shù)、巖性粒度和單位涌水量值基本不發(fā)生變化;豐水期相比枯水期，補(bǔ)給模數(shù)的變化值較大，而含水層厚度也增加0.106 m;

3）小保當(dāng)井田內(nèi)潛水含水層以富水性弱區(qū)分布最廣，且明顯大于其他區(qū)域。從豐水期到枯水期，面積和形態(tài)變化最容易受到影響的是弱富水區(qū)，而較弱富水區(qū)雖然面積變化不大，但是其分布的區(qū)域發(fā)生了較大變化。相比之下，較強(qiáng)和強(qiáng)富水含水區(qū)則受到豐水期和枯水期變化的影響較小。且分水嶺以及保德組頂面形態(tài)對(duì)富水性分區(qū)有明顯的控制作用。

參考文獻(xiàn)（References）：

[1] 李文平，段中會(huì)，華解明，等.陜北榆神府礦區(qū)地質(zhì)環(huán)境現(xiàn)狀及采煤效應(yīng)影響預(yù)測(cè)[J].工程地質(zhì)學(xué)報(bào)，2000，8（3）：324-333.LI Wen?ping，DUAN Zhong?hui，HUA Jie?ming，et al.Evaluation of present geological environment and prediction of its variation caused by mining in Yushenfu mine area of north Shaanxi[J].Journal of Engineering Geology，2000，8（3）：324-333.

[2]范立民，馬雄德，冀瑞君.西部生態(tài)脆弱礦區(qū)保水采煤研究與實(shí)踐進(jìn)展[J].煤炭學(xué)報(bào)，2015，40（8）：1711-1717.FAN Li?min，MA Xiong?de，JI Rui?jun.Progress in engineering practice of water?preserved coal mining in western eco?environment frangible area[J]. Journal of China Coal Society，2015，40（8）：1711-1717.

[3]Ghoubachi，Saad Younes.Impact of lake nasser on the groundwater of the Nubia sandstone aquifer system inTushka area，South Western Desert，Egypt[J].Journal of King Saud University Science，2012 24（2）：101-109.

[4]張?zhí)炷?，李健，施龍青，?物探技術(shù)在礦井水防治中的應(yīng)用研究[J].西安科技大學(xué)學(xué)報(bào)，2013，33（1）：56-60.ZHANG Tian?mo，LI Jian，SHI Long?qing，et al.Application of geophysical exploration technique in the mine water prevention and control[J].Journal of Xi’an University of Science and Technology，2013，33（1）：56-60.

[5]于師建.三軟煤層上覆含水層富水性瞬變電磁法探測(cè)技術(shù)[J].煤炭科學(xué)技術(shù)，2015，43（1）：104-107，73. YU Shi?jian.Transient electromagnetic method detecting technology of water?bearing occurence of aquifers in three?soft coal seam overlying rocks[J].Coal Science and Technology，2015，43（1）：104-107，73.

[6]武?強(qiáng)，張志龍，馬積福.煤層底板突水評(píng)價(jià)的新型實(shí)用方法—主控指標(biāo)體系的建設(shè)[J].煤炭學(xué)報(bào)，2007，32（1）：42-47.WU Qiang，ZHANG Zhi?long，MA Ji?fu.A new practical methodology of the coal floor water bursting evaluating Ⅱ?the master controlling index system constru?ction[J].Journal of China Coal Society，2007，32（1）：42-47.

[7]武?強(qiáng)，殷作如.評(píng)價(jià)煤層頂板涌（突）水條件的“三圖—雙預(yù)測(cè)法”[J].煤炭學(xué)報(bào)，2000，25（1）：60-65.WU Qiang，HUANG Xiao?ling，DONG Dong?lin.“Three maps?two predictions”method to evaluate water bursting conditions on roof coal[J].Journal of China Coal Society，2000，25（1）：60-65.

[8]邱士利.基于GIS的巖溶地下水富水性評(píng)價(jià)方法研究[D].北京：中國(guó)地質(zhì)科學(xué)院，2006.

QIU Shi?li. Quantitative evaluation of karst groundwater yield property based on GIS technology[D].Beijing：Chinese Academy of Geological Sciences，2006.

[9]韓?超，泮曉華，李國(guó)梁.基于GIS多源信息集成的含水層富水性模糊層次分析法[J].水文地質(zhì)工程地質(zhì)，2012，39（4）：19-25.HAN Chao，PAN Xiao?hua，LI Guo?liang，et al.The fuzzy analytic hierarchy process of water abundance of an aquifer based on GIS and multi?source information fusion techniques[J].Hydrogeology and Engineering Geology，2012，39（4）：19-25.

[10]范立民，向茂西，彭?捷，等.西部生態(tài)脆弱礦區(qū)地下水對(duì)高強(qiáng)度采煤的響應(yīng)[J].煤炭學(xué)報(bào)，2016，41（11）：2672-2678.FAN Li?min，XIANG Mao?xi，PENG Jie，et al.Groundwater response to intensive mining in ecologically fragile area[J].Journal of China Coal Society，2016，41（11）：2672-2678.

[11]彭?濤，宣良瑞，張海潮，等.臥龍湖煤礦砂巖含水層富水性預(yù)測(cè)及評(píng)價(jià)[J].煤礦安全，2014（8）：199-202.PENG Tao，XUAN Liang?rui，ZHANG Hai?chao，et al.Prediction and evaluation of water abundance of sandstone aquifer in wolonghu coal mine[J].Coal & Safety，2012，39（4）：19-25.

[12]范立民，王雙明，劉社虎，等.榆神礦區(qū)礦井涌水量特征及影響因素[J].西安科技大學(xué)學(xué)報(bào)，2009，29（1）：7-11.FAN Li?min，WANG Shuang?ming，LIU She?hu，et al.Outcome characteristics and influ?encing factor of coalmining wastewater in Yushen mining area[J].

Journal of Xi’an University of Science and Technology，2009，29（1）：7-11.

[13]邵亞紅，姚多喜，魯海峰，等.松散層底部含水層富水性評(píng)價(jià)[J].煤礦安全，2014，45（7）：127-130.SHAO Ya?hong，YAO Duo?xi，LU Hai?feng，et al.Evaluation of water abundance of loose bed bottom aquifer[J].Safty in Coal Mines，2014，45（7）：127-130.

[14]魏廣慶.板集井田新生界底部松散層富水性分布規(guī)律及其對(duì)煤層開(kāi)采影響研究[D].淮南：安徽理工大學(xué)，2005.

WEI Guang?qing.Study on aquifer parameters of the lowermost Cenozoic unconsolidated strata in the Banji Coalfield and its influence on coal mining activites[D].Huainan：Anhui University of Science and Technology， 2005.

[15]劉瑞新.松散含水層下提高開(kāi)采上限的研究與實(shí)踐[J].煤炭科學(xué)技術(shù)，2010，38（11）：56-59.LIU Rui?xin.Study and practices on improvement of mining up limit under loose aquifer[J].Coal Science and Technology，2010，38（11）：56-59.

[16]李小龍，姚多喜.孫疃礦第四含水層分布特征及富水性評(píng)價(jià)[J].安徽理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2013，33（4）：15-19.LI Xiao?long，YAO Duo?xi.Distribution characteristics and evaluation of the fourth aquifer in Suntong[J].Journal of Anhui University of Science and Technology（Natural Science），2013，33（4）：15-19.

[17]段?鵬.包氣帶巖性結(jié)構(gòu)對(duì)地下水補(bǔ)給的影響研究[D].西安：長(zhǎng)安大學(xué)，2013.

DUAN Peng.The effect of soil texture and structure of vadose zone on recharge of groundwater[D].Xi’an：Chang’an University，2013.

[18]武?強(qiáng)，樊振麗.基于GIS的信息融合型含水層富水性評(píng)價(jià)方法——富水性指數(shù)法[J].煤炭學(xué)報(bào)，2011，36（7）：1124-1128.WU Qiang，F(xiàn)AN Zhen?li.Water?richness evaluation method of water?filled aquifer based on the principle of information fusion with GIS：

water?richness index method[J].Journal of China Coal Society，2011，36（7）：1124-1128.

[19]李樹忱，馮現(xiàn)大，李術(shù)才，等.礦井頂板突水模型試驗(yàn)多場(chǎng)信息的歸一化處理方法[J].煤炭學(xué)報(bào)，2011，36（3）：447-451.LI Shu?chen，F(xiàn)ENG Xian?da，LI Shu?cai，et al.The normalization process of the multi?field information from a coal mine water?inrush model test[J].Journal of China Coal Society，2011，36（3）：447-451.

[20]Steiner B L，Antolín F，Jacomet S.Testing of the consistency of the sieving（wash?over）process of waterlogged sediments by multiple operators[J].Journal of Archaeological Science Reports，2015，2：310-320.

[21]蘭太權(quán).再論“單位涌水量就是導(dǎo)水系數(shù)”[J].水文地質(zhì)工程地質(zhì)，2016，43（5）：173-175 LAN Tai?quan.Revisit “the specific well yield is the transmissivity”[J].Hydro?Geology and Engineering Geology，2016，43（5）：173-175.

[22]葛天民.關(guān)于導(dǎo)水系數(shù)T與單位涌水量q的關(guān)系討論[J].勘察科學(xué)技術(shù)，1984（2）：17.GE Tian?min.The discussion on specific capacity and the transmissivity[J].Exploration of Science and Technology，1984（2）：17.

[23]Boucher M，F(xiàn)avreau G，Vouillamoz J M，et al.Estimating specific yield and transmissivity with magnetic resonance sounding in an unconfined sandstone aquifer（Niger）[J].Hydrogeology Journal，2009，17（7）：1805-1815.

[24]朱芮芮，鄭紅星，劉昌明.黃土高原典型流域地下水補(bǔ)給-排泄關(guān)系及其變化[J].地理科學(xué)，2010，30（1）：108-112.ZHU Rui?rui，

ZHENG Hong?xing，LIU Chang?min.Changes of ground water recharge and discharge in watershed of the loess plateau[J].Scientia Geographica Sinica，2010，30（1）：108-112.

[25]Chang C，Der Yeh H.Variability quantification of excess pressure head in heterogeneous deformable aquifers[J].Applied Mathematical Modelling，2016，40（19-20）：8580-8591.

[26]Masetti M，Pedretti D，Sorichetta A，et al.Impact of a storm?water infiltration basin on the recharge dynamics in a highly permeable aquifer[J].Water Resources Management，2016，30（1）：149-165.

[27]賀曉浪.小保當(dāng)煤礦潛水含水層富水性及受采動(dòng)影響分析[D].西安：西安科技大學(xué)，2016.

HE Xiao?lang.Water?richness evaluation of unconfined aquifer and mining effects on the aquifer in Xiaobaodang coal mine[D].Xi’an：Xi’an University of Science and Technology，2016.

[28]吳永勝，哈?斯，烏格特茉勒.毛烏素沙地南緣沙丘表面徑流特征[J].科學(xué)通報(bào)，2011，56（34）：2917-2922.

WU Yong?sheng， Hasi E， Wugetemole. Characteristics of surface runoff in a sandy area in southern Mu Us sandy land[J].Chinese Science Bulletin，2011，56（34）：2917-2922.

[29]Nishigaki T，Sugihara S，Kilasara M，et al.Surface runoff generation and soil loss under different soil and rainfall properties in the Uluguru Mountains，Tanzania[J].Land Degradation and Development，2017，28（1）：283-293.