涇惠渠灌區(qū)地下水位埋深的動(dòng)態(tài)特征及影響因素分析

張曼菲　周維博

摘要：地下水作為陜西涇惠渠灌區(qū)的重要水源，對(duì)其動(dòng)態(tài)特征的研究在灌區(qū)水資源開(kāi)發(fā)利用方面具有重要的指導(dǎo)意義。根據(jù)涇惠渠灌區(qū)1981-2010年的地下水位平均埋深資料，采用線(xiàn)性?xún)A向估計(jì)法、累積距平法、Sen斜率估計(jì)法、MannKendall趨勢(shì)檢驗(yàn)法對(duì)涇惠渠灌區(qū)地下水位埋深動(dòng)態(tài)序列進(jìn)行年際變化分析；采用滑動(dòng)T檢驗(yàn)法和雙向回歸突變分析法分析序列突變性；采用主成分分析法和灰色關(guān)聯(lián)分析法確定影響灌區(qū)地下水埋深動(dòng)態(tài)的主要因素。結(jié)果表明：灌區(qū)1981-2010年地下水位平均埋深基本呈增加趨勢(shì)，突變點(diǎn)出現(xiàn)在1994年，1994年之前影響灌區(qū)地下水埋深動(dòng)態(tài)的主要因素為蒸發(fā)量，1994年之后主要影響因素為蒸發(fā)量、地下水開(kāi)采量和灌渠引水量。

關(guān)鍵詞：涇惠渠灌區(qū)；地下水位埋深；趨勢(shì)分析；突變檢驗(yàn)；影響因素

中圖分類(lèi)號(hào)：P641文獻(xiàn)標(biāo)志碼：A文章編號(hào)：

16721683（2018）02016407

Abstract：

Groundwater is an important water source of Jinghuiqu Irrigation District.The analysis of its dynamic characteristics provides important guidance to the development and utilization of water resources.Based on the groundwater depth data of Jinghuiqu Irrigation District during 19812010，we analyzed its interannual variation using linear tendency estimation method，cumulative departure，Sen slope estimation method，and MannKendall trend test.We analyzed the mutation of groundwater depth using smoothing T test and two-dimensional regression mutation analysis.We determined the most important influence factors of groundwater depth using principal component analysis and grey relationship analysis.The results showed that the groundwater depth tended to increase during 19812010.The mutation happened in 1994.Before 1994，evaporation was the most important influence factor of groundwater depth.After 1994，the most important influence factors became evaporation，groundwater extraction，and canal water diversion.

Key words：

Jinghuiqu Irrigation District；groundwater depth；trend analysis；mutation test；influence factors

自20世紀(jì)80年代以來(lái)，涇惠渠灌區(qū)地下水位持續(xù)下降，對(duì)灌區(qū)的生態(tài)環(huán)境帶來(lái)了許多不利的影響[1]，不利于灌區(qū)的可持續(xù)發(fā)展。因此對(duì)灌區(qū)地下水動(dòng)態(tài)特征和影響因子的分析研究顯得尤為重要。

已有大量學(xué)者對(duì)涇惠渠灌區(qū)的地下水動(dòng)態(tài)開(kāi)展了大量的研究工作。劉燕[1]通過(guò)分析灌區(qū)地下水位下降的自然和人為原因，指出地表引水灌溉的減少，造成了地下水采補(bǔ)失衡，從而使得地下水位大幅下降；李娜等[23]利用ARIMA模型、ANFIS模型對(duì)灌區(qū)地下水動(dòng)態(tài)進(jìn)行預(yù)測(cè)，指出灌區(qū)地下水位將持續(xù)下降；鄧康婕等[4]利用趨勢(shì)分析法和統(tǒng)計(jì)學(xué)方法，研究了灌區(qū)地下水埋深的時(shí)空變異性；宋寶德等[5]分析了灌區(qū)典型干旱時(shí)段的地下水動(dòng)態(tài)特征；王建瑩等[6]利用ArcGIS軟件分析了灌區(qū)地下水動(dòng)態(tài)及其影響因素。但以上研究大多是基于地下水動(dòng)態(tài)在時(shí)間和空間上的變化趨勢(shì)，鮮少涉及地下水動(dòng)態(tài)變化的突變特征及其影響因素分析。對(duì)于水文序列的突變研究大多數(shù)集中在降水、徑流、氣溫等研究領(lǐng)域，采用的研究方法主要有MannKendall檢驗(yàn)[710]、Yamamoto檢驗(yàn)[1113]、滑動(dòng)T檢驗(yàn)[1416]等，以上方法均適用于均值突變、反向趨勢(shì)突變的檢驗(yàn)，但在解決地下水動(dòng)態(tài)這類(lèi)同向趨勢(shì)突變的問(wèn)題中有一定的局限性，而雙向回歸分析法在此類(lèi)問(wèn)題的分析中十分有效[1718]。

因此本文收集了涇惠渠灌區(qū)1981-2010年的地下水位平均埋深資料，采用趨勢(shì)分析法和突變檢驗(yàn)法，確定地下水位平均埋深序列的突變年份，并探討突變年份前后影響灌區(qū)地下水動(dòng)態(tài)的主要因素。

1研究區(qū)概況

涇惠渠灌區(qū)總面積1 180 km2，位于陜西省關(guān)中平原中北部，地處東經(jīng)108°34′34″-109°21′35″，北緯34°25′20″-34°41′40″，是一個(gè)自涇河自流引水、灌排結(jié)合，并以農(nóng)業(yè)灌溉為主的井渠雙灌灌區(qū)。地勢(shì)自西北向東南傾斜，海拔高程350～450 m，地面坡降1/300～1/600。灌區(qū)位置及范圍詳見(jiàn)圖1。

涇惠渠灌區(qū)屬中緯度暖溫帶半干旱氣候，夏季高溫，蒸發(fā)大降雨多且集中；冬季寒冷干燥，降水稀少。多年平均降水量5389 mm，時(shí)空分布不均，多年平均蒸發(fā)量1 212 mm，多年平均氣溫134 ℃。

灌區(qū)內(nèi)第四系地層發(fā)育，自西向東、由北向南沉積厚度遞增，顆粒由粗變細(xì)，根據(jù)地質(zhì)構(gòu)造和地形地貌將灌區(qū)劃分為4個(gè)水文地質(zhì)分區(qū)和7個(gè)亞區(qū)，如圖1所示。區(qū)內(nèi)地下水類(lèi)型主要有潛水和承壓水，潛水是當(dāng)?shù)刂饕喔人?，埋深較淺；承壓水埋深大于100 m。[HJ1〗

近年來(lái)，隨著經(jīng)濟(jì)的快速發(fā)展，井渠灌溉用水比由20世紀(jì)70年代的0.3左右上升至目前的1.3左右，大大增加了地下水開(kāi)采量，使得灌區(qū)地下水超采情況嚴(yán)重，引起灌區(qū)地下水資源量減少、降落漏斗不斷擴(kuò)大等問(wèn)題[19]。因此研究灌區(qū)地下水動(dòng)態(tài)特征及其影響因素，對(duì)今后灌區(qū)的水資源開(kāi)發(fā)利用有重要的指導(dǎo)意義。

2研究方法

對(duì)地下水位平均埋深序列的突變檢驗(yàn)選擇線(xiàn)性?xún)A向估計(jì)法[20]、累積距平法[20]、Sen斜率估計(jì)法[21]和MannKendall趨勢(shì)檢驗(yàn)法[20]。前兩種方法的優(yōu)點(diǎn)在于直觀；Sen斜率估計(jì)法與基于最小二乘法的線(xiàn)性?xún)A向估計(jì)法相比，優(yōu)點(diǎn)在于可避免序列數(shù)據(jù)缺失和分布形態(tài)、以及異常值對(duì)分析結(jié)果的影響[22]；而對(duì)于地下水位平均埋深此類(lèi)非正態(tài)分布的數(shù)據(jù)序列，MannKendall趨勢(shì)檢驗(yàn)法的應(yīng)用更為廣泛，該法不要求序列滿(mǎn)足一定的分布，且能夠剔除異常值的干擾[23]。

對(duì)地下水位平均埋深序列的突變檢驗(yàn)選擇滑動(dòng)T檢驗(yàn)法[15]和雙向回歸突變分析法[1718]?；瑒?dòng)T檢驗(yàn)法是常見(jiàn)的突變檢驗(yàn)方法，優(yōu)點(diǎn)是計(jì)算簡(jiǎn)單，通過(guò)考察兩組樣本均值的差異顯著性檢驗(yàn)變異點(diǎn)，但缺點(diǎn)是子序列長(zhǎng)度人為控制性強(qiáng)，只能得到可能的突變點(diǎn)[23]。雙向回歸突變分析法除了前述所說(shuō)對(duì)同向趨勢(shì)突變有效外，其優(yōu)點(diǎn)還在于對(duì)序列的逐個(gè)點(diǎn)進(jìn)行分析檢驗(yàn)，從而確定序列的突變點(diǎn)。由于雙向回歸突變分析法在地下水領(lǐng)域的應(yīng)用較少，現(xiàn)將該方法介紹如下：

假設(shè)時(shí)間序列的突變點(diǎn)為C點(diǎn)，則此序列的雙向回歸模型可表示為

3.2突變分析

對(duì)灌區(qū)1981-2010年地下水位平均埋深序列進(jìn)行滑動(dòng)T檢驗(yàn)，n=30，選取步長(zhǎng)為6，顯著性水平α=0001，自由度[WTB1X]v[WTBX]=10，查表得知t0.001=±4587?；瑒?dòng)T統(tǒng)計(jì)量曲線(xiàn)見(jiàn)圖4（a），T的統(tǒng)計(jì)量共有7次超過(guò)0001顯著性水平，即存在7個(gè)可能的突變點(diǎn)，分別為1989、1991、1992、1993、1994、1999和2004年。

為進(jìn)一步確定突變的年份，采用雙向回歸突變分析法對(duì)灌區(qū)1981-2010年地下水位平均埋深序列再次進(jìn)行趨勢(shì)突變分析。由式12可知，U值與假定突變點(diǎn)C對(duì)應(yīng)的殘差平方和S呈負(fù)相關(guān)，因此利用U值與假定突變點(diǎn)C對(duì)應(yīng)時(shí)間的負(fù)相關(guān)關(guān)系來(lái)確定突變點(diǎn)，突變點(diǎn)即為通過(guò)0001顯著性水平檢驗(yàn)的最大值，查表得F（3，26）=736，結(jié)果見(jiàn)圖4（b），U的最大值點(diǎn)出現(xiàn)在1994年。綜合上述兩種方法，得出灌區(qū)地下水位平均埋深序列的突變點(diǎn)出現(xiàn)在1994年。

4影響因子分析

4.1影響因子的選取

根據(jù)地下水均衡方程式，含水層中地下水量的變化取決于兩項(xiàng)：補(bǔ)給項(xiàng)和排泄項(xiàng)。

根據(jù)現(xiàn)有的涇惠渠灌區(qū)資料，地下水補(bǔ)給主要來(lái)源有大氣降水入滲補(bǔ)給、河流水入滲補(bǔ)給以及灌溉入滲補(bǔ)給，因此選擇年降水量、年徑流量、渠灌引水量這3個(gè)指標(biāo)，此外涇惠渠灌區(qū)屬于井渠雙灌灌區(qū)，因此表示渠灌水量與井灌水量用水比例的指標(biāo)—渠井用水比也是影響灌區(qū)地下水動(dòng)態(tài)的因素。

灌區(qū)地下水排泄項(xiàng)主要有蒸發(fā)排泄和人為開(kāi)采，因此選擇蒸發(fā)量和地下水開(kāi)采量這2個(gè)指標(biāo)。

4.2影響因子的年際變化

[JP+1]上述7個(gè)指標(biāo)隨時(shí)間的變化曲線(xiàn)見(jiàn)圖5。圖5（a）為4個(gè)自然因素指標(biāo)的變化曲線(xiàn)，可以看出，年降水量和年徑流量隨時(shí)間整體呈減小趨勢(shì)，蒸發(fā)量整體呈增大趨勢(shì)，年平均氣溫的變化并不明顯；圖5（b）為3個(gè)人為因素指標(biāo)的變化曲線(xiàn)，渠灌引水量和渠井用水比整體呈減小趨勢(shì)，而地下水開(kāi)采量整體呈增大趨勢(shì)。

為分析各因素在突變年1994年前后的變化規(guī)律，利用線(xiàn)性?xún)A向回歸估計(jì)法計(jì)算突變年前后各因素變化曲線(xiàn)的回歸系數(shù)，計(jì)算結(jié)果見(jiàn)表3。在突變年前后，年徑流量和渠灌引水量先減小后增大，渠井用水比一直呈現(xiàn)減小趨勢(shì)，其余因素均呈現(xiàn)增大趨勢(shì)。回歸系數(shù)絕對(duì)值的大小可以反映序列變化速率的大小。

由表3可知，年降水量、渠灌引水量和渠井用水比在1994年之前的變化率大于1994年之后，尤其是渠灌引水量；其余4個(gè)指標(biāo)在1994年之前的變化率均小于1994年之后，且蒸發(fā)量的變化最為明顯

4.3基于主成分分析法的影響因子提取

利用SPSS軟件得到突變年前后7個(gè)變量指標(biāo)的相關(guān)系數(shù)矩陣，由表4可知，7個(gè)指標(biāo)間存在不同程度的相關(guān)性。在突變前，相關(guān)系數(shù)達(dá)到07以上的有年降水量（I1）和年徑流量（I2）、年徑流量（I2）和地下水開(kāi)采量（I6）、蒸發(fā)量（I3）和地下水開(kāi)采量（I6），其中年徑流量（I2）和地下水開(kāi)采量（I6）呈負(fù)相關(guān)，其余兩組呈正相關(guān)。在突變后，相關(guān)系數(shù)最高的是地下水開(kāi)采量（I6）和渠井用水比（I7），但相關(guān)系數(shù)僅為0.647。

利用SPSS軟件對(duì)上述7個(gè)指標(biāo)提取主成份。經(jīng)檢驗(yàn)，突變前后地下水埋深序列的樣本充足度分別為0541和0655，均大于05，可以進(jìn)行主成分分析。以變量指標(biāo)特征值大于1且累計(jì)貢獻(xiàn)率達(dá)到80%以上為準(zhǔn)則，分別提取1981-1994年和1995-2010年7個(gè)指標(biāo)的公共因子，結(jié)果見(jiàn)表5。1994年前后的兩個(gè)序列均提取到3個(gè)主因子，1981-1994年3個(gè)主因子的累計(jì)貢獻(xiàn)率達(dá)到85664%，1995-2010年3個(gè)主因子的累計(jì)貢獻(xiàn)率達(dá)到83579%。

在旋轉(zhuǎn)因子荷載矩陣（表6）中剔除絕對(duì)值小于0.8的因子，提取出突變前后影響地下水位平均埋深的主要因素。1994年之前影響地下水位平均埋深的第一主成分主要由蒸發(fā)量決定，相關(guān)系數(shù)達(dá)到

綜上所述，1994年之前，灌區(qū)地下水位平均埋深主要受年降水量、年徑流量、蒸發(fā)量、渠灌引水量這4個(gè)因素的影響，可以看出自然因素是影響1994年之前地下水位平均埋深的主要因素。1994年之后，灌區(qū)地下水位平均埋深主要受年降水量、蒸發(fā)量、年平均氣溫、渠灌引水量、地下水開(kāi)采量、渠井用水比這6個(gè)因素的影響，自然因素和人為因素均影響1994年后的地下水位平均埋深。

4.4基于灰色關(guān)聯(lián)度法的主導(dǎo)因子確定

采用灰色關(guān)聯(lián)度法分別確定灌區(qū)突變年前后，各影響因子對(duì)灌區(qū)地下水位平均埋深的影響程度。根據(jù)上面對(duì)灰色關(guān)聯(lián)度法的介紹，分別計(jì)算1981-1994年年降水量、年徑流量、蒸發(fā)量、渠灌引水量這4個(gè)因素的對(duì)地下水位平均埋深影響程度，1995-2010年年降水量、蒸發(fā)量、年平均氣溫、渠灌引水量、地下水開(kāi)采量、渠井用水比這6個(gè)因素對(duì)地下水位平均埋深影響程度，計(jì)算結(jié)果見(jiàn)表7。將這些因素與地下水位平均埋深的關(guān)聯(lián)度從大到小排序，1981-1994年關(guān)聯(lián)度排序?yàn)檎舭l(fā)量>年徑流量>年降水量>灌渠引水量，1995-2010年關(guān)聯(lián)度排序?yàn)檎舭l(fā)量>地下水開(kāi)采量>灌渠引水量>氣溫>年降水量>渠井用水比。

將關(guān)聯(lián)度大于0.7的因子視為主導(dǎo)因子，則1981-1994年影響地下水位平均埋深的主導(dǎo)因子為蒸發(fā)量，1995-2010年影響地下水位平均埋深的主導(dǎo)因子依次為蒸發(fā)量、地下水開(kāi)采量和灌渠引水量，表明灌區(qū)地下水動(dòng)態(tài)的影響因素由自然因素轉(zhuǎn)變?yōu)槿藶橐蛩睾妥匀灰蛩亍?/p>

5結(jié)論

（1）灌區(qū)1981-2010年地下水位平均埋深基本呈增加趨勢(shì)，平均增加速率為0.34 m/a，1981-1994年灌區(qū)水位平均埋深值呈現(xiàn)減小趨勢(shì)，且變化率較快，1994年以后水位平均埋深值呈現(xiàn)增大趨勢(shì)，變化率變慢。

（2）對(duì)灌區(qū)1981-2010年地下水位平均埋深序列采用滑動(dòng)T檢驗(yàn)法和雙向回歸突變分析法，得出灌區(qū)地下水位平均埋深序列的突變點(diǎn)出現(xiàn)在1994年。

（3）利用主成分分析法提取突變年前后影響灌區(qū)地下水位平均埋深的主要影響因子，結(jié)果表明灌區(qū)地下水位平均埋深在1994年之前，主要受年降水量、年徑流量、蒸發(fā)量、渠灌引水量這4個(gè)因素的影響；1994年之后，主要受年降水量、蒸發(fā)量、年平均氣溫、渠灌引水量、地下水開(kāi)采量、渠井用水比這6個(gè)因素的影響。采用灰色關(guān)聯(lián)度法確定影響突變年前后灌區(qū)地下水位平均埋深的主導(dǎo)因子，結(jié)果表明，1981-1994年主導(dǎo)因子為蒸發(fā)量，1995-2010年主導(dǎo)因子為蒸發(fā)量、地下水開(kāi)采量和灌渠引水量，灌區(qū)地下水動(dòng)態(tài)的影響因素由自然因素轉(zhuǎn)變?yōu)槿藶橐蛩睾妥匀灰蛩亍?/p>

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

[1]劉燕.涇惠渠灌區(qū)地下水位動(dòng)態(tài)變化特征及成因分析[J].人民長(zhǎng)江，2010，41（8）：100103，107.（LIU Y.Dynamics of groundwater table in Jinghui Canal Irrigation District of China[J].Yangtze River，2010，41（8）：100103，107.（in Chinese））.DOI：10.3969/ j.issn.10014179.2010.08.027.

[2]李娜，周維博，董起廣，等.基于ARIMA模型的灌區(qū)地下水動(dòng)態(tài)預(yù)報(bào)[J].人民長(zhǎng)江，2014（9）：4347.（LI N，ZHOU W B，DONG G Q，et al.Groundwater dynamic prediction in irrigation area based on ARIMA Model[J].Yangtze River，2014（9）：4347.（in Chinese））.DOI：10.3969/ j.issn.10014179.2014.09.011.

[3]李娜，周維博，王金鳳，等.基于主成分分析法的ANFIS模型及其應(yīng)用[J].人民黃河，2014（6）：8083.（LI N，ZHOU W B，WAGN J F，et al.ANFIS Model and its application in groundwater prediction in Irrigation District based on principal component[J].Yellow River，2014（6）：8083.（in Chinese））.DOI：10.3969/ J.issn.10001379.2014.06.026.

[4]鄧康婕，魏曉妹，降亞楠，等.基于地統(tǒng)計(jì)學(xué)的涇惠渠灌區(qū)地下水位時(shí)空變異性研究[J].灌溉排水學(xué)報(bào)，2015，34（3）：7580.（DENG K J，WEI X M，JIANG Y N，et al.Spatiotemporal variability of groundwater table in Jinghui Irrigation District based on geostatistics[J].Journal of Irrigation and Drainage，2015，34（3）：7580.（in Chinese））.DOI：10.13522/ j.cnki.ggps.2015.03.016.

[5]宋寶德，鄭新崳，吳建華，等.涇惠渠灌區(qū)干旱時(shí)段內(nèi)的地下水動(dòng)態(tài)[J].水利科技與經(jīng)濟(jì)，2010，16（4）：386388.（SONG B D，ZHEN X Y，WU J H，et al.Groundwater dynamics of Jinghuiqu Irrigation Area in arid period[J].Water Conservancy Science and Technology and Economy，2010，16（4）：386388.（in Chinese））.DOI：10.3969/ j.issn.1006 7175.2010.04.011.

[6]王建瑩，劉燕，姚阿漫，等.涇惠渠灌區(qū)地下水位動(dòng)態(tài)分析[J].灌溉排水學(xué)報(bào)，2015，34（2）：6770.（WNAG J Y，LIU Y，Y A M，et al.Dynamics of the groundwater table in Jinghui Canal Irrigation District of China[J].Journal of Irrigation and Drainage，2015，34（2）：6770.（in Chinese））.DOI：10.13522/ j.cnki.ggps.2015.02.016.

[7]鄭江禹，張強(qiáng)，史培軍，等.珠江流域多尺度極端降水時(shí)空特征及影響因子研究[J].地理科學(xué)，2017，37（2）：283291.（ZHENG J Y，ZHANG Q，SHI P J，et al.Spatiotemporal Characteristics of extreme precipitation regimes and related driving factors in the Pearl River Basin[J].Scientia Geographica Sinica，2017，37（2）：283291.（in Chinese））.DOL：10.13249/j.cnki.sgs.2017.02.015.

[8]孫寧，程伍群.河北燕山山區(qū)徑流變化規(guī)律研究[J].河北農(nóng)業(yè)大學(xué)學(xué)報(bào)，2017，40（2）：123128.（SUN N，CHENG W Q.Study on the runoff variation rules of Hebei Yanshan Mountains[J].Journal of Agricultural University of Hebei，2017，40（2）：123128.（in Chinese））.DOI：10.13320 /j.cnki.jauh.2017.0048.

[9]蘇向明，劉志輝，魏天峰，等.新疆精河徑流變化及其原因分析[J].水文，2016，36（5）：9296，91.（SU X M，LIU Z H，WEI T F，et al.Analysis of runoff variation and its causes in Jinghe Basin，Xinjiang[J]，Journal of China Hydrology，2016，36（5）：9296，91.（in Chinese））.DOI：10.3969/ j.issn.10000852.2016.05.018.

[10][ZK（#]張范平，方少文，周祖昊，等.鄱陽(yáng)湖水位多時(shí)間尺度動(dòng)態(tài)變化特性分析[J].長(zhǎng)江流域資源與環(huán)境，2017，26（1）：126133.（ZHANG F P，F(xiàn)ANG S W，ZHOU Z H，et al.Research on multitimescale dynamic characteristics of waterlevel fluctuation of the Poyang Lake in China[J].Resources and Environment in the Yangtze Basin，2017，26（1）：126133.（in Chinese））.DOI：10.11870 /cjlyzyyhj201701015.

[11]李艷玲，暢建霞.基于Morlet小波的徑流突變檢測(cè)[J].西安理工大學(xué)學(xué)報(bào)，2012，28（3）：322325.（LI Y L，CHANG J X.Detection of the abrupt changes in runoff based on the Morlet Wavelet[J].Journal of Xi′an University of Technology，2012，28（3）：322325.（in Chinese））.DOI：10.3969/j.issn.10064710.2012.03.014.

[12]史建橋，白淑英，周寅，等.近50年甘肅日照時(shí)數(shù)時(shí)空變化特征及突變分析[J].南京信息工程大學(xué)學(xué)報(bào)，2016，8（1）：6470.（SHI J Q，BAI S Y，ZHOU Y，et al.Spatiotemporal variation and abrupt change analysis of sunshine duration in Gansu in the last 50 years[J].Journal of Nanjing University of Information Science & Technology，2016，8（1）：6470.（in Chinese））.DOI：10.13878/j.cnki.jnuist.2016.01.007.

[13]王秀萍，金巍.19642013年大連地區(qū)暴雨氣候特征及變化規(guī)律[J].氣象與環(huán)境學(xué)報(bào)，2015（3）：7580.（WANG X P，JIN W.Climate characteristics of rainstorm and its variation in Dalian region from 1964 to 2013[J].Journal of Meteorology and Environment，2015 （3）：7580.（in Chinese））.DOI：10.3969/j.issn.1673503X.2015.03.012.

[14]李斌，解建倉(cāng)，胡彥華，等.1961年2013年渭河流域降水與徑流變化特征[J].南水北調(diào)與水利科技，2017，15（2）：2936.（LI B，XIE J C，HU Y H，et al.Characteristics of precipitation and runoffs during 19612013 in Weihe River basin[J].SouthtoNorth Water Transfers and Water Science & Technology，2017，15（2）：2936.（in Chinese））.DOI：10.13476/j.cnki.nsbdqk.2017.02.005.

[15]吳凱，王曉琳，王高旭，等.1961年2014年西遼河流域降水時(shí)空變異性診斷[J].南水北調(diào)與水利科技，2017，15（2）：2228.（WU K，WANG X L，WANG G X，et al.Spatialtemporal variation of precipitation in West Liao River Basin during 19612014[J].SouthtoNorth Water Transfers and Water Science & Technology，2017，15（2）：2228.（in Chinese））.DOI：10.13476/j.cnki.nsbdqk.2017.02.004.

[16]李榮波，魏鵬，紀(jì)昌明，等.雅礱江流域近60 a徑流趨勢(shì)特征分析[J].人民長(zhǎng)江，2017，48（5）：3842.（LI R B，WEI P，JI C M，et al.Analysis of runoff tendency characteristics of Yalong River for past 6 decades[J].Yangtze River，2017，48（5）：3842.（in Chinese））.DOI：10.16232/j.cnki.10014179.2017.05.009.

[17]危潤(rùn)初，肖長(zhǎng)來(lái)，方樟.黑龍江建三江地區(qū)地下水動(dòng)態(tài)趨勢(shì)突變點(diǎn)分析[J].吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2016，（1）：202210.（WEI R C，XIAO C L，F(xiàn)ANG Z.Trends mutation nodes of groundwater dynamic in Jiansanjiang Area of Heilongjiang Province[J].Journal of Jilin University（Earth Science Edition），2016，（1）：202210.（in Chinese））.DOI：10.13278/j.cnki.jjuese.2016.01.204.

[18]林卓宏，梁敏妍，梁軍，等.江門(mén)市雷暴多發(fā)期氣候特征研究[J].氣象，2012（12）：15021507.（LIN Z H，LIANG M Y，LIANG J，et al.A study on climatic characters in high frequency period of thunderstorm in Jiangmen of Guangdong province[J].Meteorological Monthly，2012 （12）：15021507.（in Chinese））

[19]李慧，周維博，馬聰，等.涇惠渠灌區(qū)地下水調(diào)蓄研究[J].干旱區(qū)研究，2015，32（4）：828833.（LI H，ZHOU W B，MA C，et al.Research on regulating capacity of groundwater in Jinghui Irrigation District[J].Arid Zone Research，2015，32（4）：828833（（in Chinese））DOI：10.13866/j.azr.2015.04.29.