徐月琴，孫景林，王沁宇

（1.濰坊市水文局，山東 濰坊 261061；2.山東正大地理信息工程有限公司，山東 濰坊 261000）

為綜合分析濰坊市濰北區(qū)（含昌邑市、壽光市、寒亭區(qū)）地下水位降落漏斗的變化趨勢(shì)，采用2001—2012年地下水漏斗中心點(diǎn)監(jiān)測(cè)資料。經(jīng)對(duì)地下水位資料進(jìn)行整理、比較，最終采用18眼水位監(jiān)測(cè)井的實(shí)測(cè)資料進(jìn)行分析計(jì)算。

1　漏斗區(qū)地下水位的影響因素

根據(jù)2001—2012年《濰坊市水資源公報(bào)》資料，對(duì)濰北區(qū)有關(guān)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析計(jì)算。

1.1　漏斗區(qū)地下水位與開(kāi)采量的關(guān)系

以2012年勘查發(fā)現(xiàn)的壽西水源地抽水時(shí)對(duì)地下水位的影響為例，研究地下水開(kāi)采量與地下水位的關(guān)系。通過(guò)繪制壽西水源地開(kāi)采量—水位—降水量綜合過(guò)程線分析看出，壽西水源地抽水時(shí)，7口井的水位都有不同程度的下降。從整體上看，地下水位的波動(dòng)與地下水開(kāi)采量基本一致，地下水開(kāi)采量是地下水漏斗形成的主要因素。同時(shí)，運(yùn)用SPSS軟件，分析地下水位與開(kāi)采量的關(guān)系，以驗(yàn)證這一結(jié)論。選用研究區(qū)18眼監(jiān)測(cè)井地下水位平均值進(jìn)行計(jì)算，運(yùn)行結(jié)果為顯著性0.002，相關(guān)系數(shù)為-0.773。表明濰北區(qū)的平均地下水位與地下水開(kāi)采量均存在負(fù)相關(guān)，即隨著開(kāi)采量的增加，地下水位持續(xù)下降，地下水位的值越小。

由相關(guān)分析可知，皮爾遜相關(guān)系數(shù)γ=-0.773，零假設(shè)成立概率P=0.002，具有統(tǒng)計(jì)學(xué)意義。其中，0.5＜ γ =0.773＜0.8，顯著相關(guān)，且為負(fù)相關(guān)。因此，漏斗中心水位和地下水開(kāi)采量存在負(fù)相關(guān)，即地下水開(kāi)采量越大，地下水位的值越小。確定二者相關(guān)后，進(jìn)行線性回歸分析和曲線擬合，確定濰北區(qū)漏斗水位與地下水開(kāi)采量的函數(shù)關(guān)系。

由圖1知，線性、對(duì)數(shù)和二次曲線對(duì)應(yīng)顯著水平均小于0.05，可以認(rèn)為方程顯著；R2值越大反映的兩變量的共變量比率越高，模型與數(shù)據(jù)的擬合程度越好。二次曲線函數(shù)R2最大，優(yōu)于其他函數(shù)，因此認(rèn)為濰北區(qū)漏斗水位與地下水開(kāi)采量符合二次曲線函數(shù)關(guān)系（見(jiàn)圖1）。

圖1　濰北區(qū)擬合曲線

1.2　漏斗區(qū)地下水位與社會(huì)經(jīng)濟(jì)發(fā)展的關(guān)系

根據(jù)資料統(tǒng)計(jì)，濰北區(qū)10余年來(lái)，GDP由268.24億元漲至1 287.53億元，增長(zhǎng)了4.8倍。選取濰北區(qū)漏斗中心的多個(gè)監(jiān)測(cè)井取平均值，進(jìn)行相關(guān)分析。經(jīng)分析得知，漏斗中心水位資料和GDP存在負(fù)相關(guān)，并隨著區(qū)域經(jīng)濟(jì)社會(huì)迅速發(fā)展，水資源的需求量不斷增大，導(dǎo)致地下水開(kāi)采量增加，漏斗區(qū)地下水位下降。因此，社會(huì)經(jīng)濟(jì)的發(fā)展作為間接因素影響著地下水位的變化。

1.3　漏斗區(qū)地下水位與降水的關(guān)系

地下水位的波動(dòng)主要受地下水開(kāi)采量的影響，年際變化與降水量的波動(dòng)不一致。在考慮降水量的影響時(shí)，采用2001—2006年的數(shù)據(jù)進(jìn)行分析，便于較好地反映降水對(duì)地下水位的影響。計(jì)算漏斗中心平均地下水位與年降水量、豐枯情況見(jiàn)表1。

表1　降水與地下水位分析表

由表1可以看出，2002年、2006年為枯水年，地下水位明顯下降，偏豐水年地下水位小幅上升。因此，可認(rèn)為降水的多寡對(duì)地下水位有影響，降水是影響漏斗演變的因素之一。

1.4　漏斗區(qū)地下水位與蒸發(fā)的關(guān)系

資料采用2001—2006年區(qū)域內(nèi)實(shí)測(cè)的數(shù)據(jù)進(jìn)行分析，計(jì)算其干旱指數(shù)如表2。

表2　干旱指數(shù)與地下水位分析表

從表2看出，受降水、蒸發(fā)的影響，2002年、2006年r≥3，屬于半干旱，地下水位明顯下降；其余年r＜2，屬于半濕潤(rùn)，地下水位有所上升。因此，蒸發(fā)對(duì)地下水位也有影響，在降水和蒸發(fā)的作用下漏斗區(qū)地下水位有小幅波動(dòng)。

2　地下水位下降對(duì)生態(tài)環(huán)境的影響

2.1　地下水漏斗區(qū)水位不斷下降

根據(jù)地下水監(jiān)測(cè)資料，濰北區(qū)地下水大幅下降始于1983年，最大點(diǎn)發(fā)生在壽光留呂一帶，最大水位埋深由10.68 m下降到1999年的35.0 m，共下降24.32 m，年平均下降1.52 m，面積110 km2。該漏斗自1995年至2000年，最大水位埋深面積261 km2。因此，過(guò)度開(kāi)采地下水，導(dǎo)致地下水位持續(xù)下降，造成海水入侵、地下水污染等一系列的生態(tài)問(wèn)題，使得本來(lái)就很脆弱的生態(tài)環(huán)境進(jìn)一步惡化。

2.2　地下水漏斗區(qū)形成并不斷擴(kuò)大

由于地下水長(zhǎng)期開(kāi)采量的不斷增加，咸水入侵速度明顯加快，入侵面積不斷擴(kuò)大，據(jù)統(tǒng)計(jì)：1980年濰北區(qū)漏斗面積265 km2，到1985年漏斗面積擴(kuò)大到885 km2。隨著地下水位不斷下降，南側(cè)的壽光、五臺(tái)、昌邑城區(qū)漏斗區(qū)不斷擴(kuò)展連片，1988年已擴(kuò)展成為壽光、昌邑大漏斗，漏斗區(qū)面積擴(kuò)大為1 675 km2，1995年漏斗區(qū)面積達(dá)1 833 km2。由于采取了地下水開(kāi)采控制措施，漏斗區(qū)面積逐漸縮小，到2013年漏斗區(qū)面積為1 605 km2。

2.3　咸水入侵加劇

由于咸水分界線近淡水側(cè)地下水位不斷下降，破壞了咸淡水原有的平衡，導(dǎo)致研究區(qū)海咸水入侵呈加重趨勢(shì)。自1983年形成咸水入侵，到1990年咸水入侵面積272 km2；2002年咸水入侵面積401 km2，到2011年咸水入侵面積752 km2，2011年與2002年相比增長(zhǎng)了351 km2，年入侵速率為35.1 km2/a。近年來(lái)，隨著地下水開(kāi)采量的逐年減少，年咸水入侵速率有所緩慢減少。

2.4　地下水污染

由于地下水位持續(xù)下降，特別是地下水漏斗區(qū)的形成和不斷擴(kuò)展，很容易加速地表污水的下滲速度，加劇地下水的污染程度，給當(dāng)?shù)厮鷳B(tài)環(huán)境、人們生產(chǎn)生活帶來(lái)極大的危害。根據(jù)2012年濰坊市地下水質(zhì)監(jiān)測(cè)資料結(jié)果，壽光達(dá)到Ⅲ類(lèi)水標(biāo)準(zhǔn)，有害成分均未超標(biāo)，水質(zhì)較好。寒亭和昌邑為Ⅳ類(lèi)水，主要化學(xué)組成分多數(shù)有不同程度的增加，水質(zhì)變差。

3　結(jié)論

通過(guò)以上分析，結(jié)果表明：地下水開(kāi)采量是影響地下水位降落漏斗的最主要因素，平均地下水位與地下水開(kāi)采量均存在負(fù)相關(guān)。社會(huì)經(jīng)濟(jì)的迅速發(fā)展作為間接因素，也影響著地下水位的變化，濰北區(qū)地下水位變化與GDP存在負(fù)相關(guān)。此外，降水和蒸發(fā)也對(duì)超采區(qū)地下水位有一定影響，年降水量越少，蒸發(fā)量越大，地下水位越低。