李 丹， 李志威， 張明珠， 龐 園

(廣州市水務(wù)科學(xué)研究所， 廣東 廣州 510220)

1 研究背景

廣花盆地是廣州市重要的地下水應(yīng)急備用水源地[1-2]，同時(shí)也是珠江三角洲巖溶塌陷及地面沉降地質(zhì)災(zāi)害重點(diǎn)防治區(qū)[3]。盆地內(nèi)自來水暫未全面供給，分散性開采巖溶水的工業(yè)及商業(yè)機(jī)井仍較多，由于不注意合理布局和科學(xué)開采，歷史上廣花盆地局部地區(qū)誘發(fā)了地面沉降、塌陷、地裂縫、房裂等地質(zhì)問題[4-6]。如何讓這一重要的應(yīng)急備用水源地在水量和水質(zhì)上滿足人類社會(huì)和生態(tài)系統(tǒng)長(zhǎng)期穩(wěn)定發(fā)展的需要，且能有效避免地質(zhì)-生態(tài)災(zāi)害，則需要對(duì)地下水資源可持續(xù)開發(fā)利用問題進(jìn)行探討。

國(guó)內(nèi)學(xué)者對(duì)廣花盆地開展過地面塌陷成因、特征、形成機(jī)理、風(fēng)險(xiǎn)評(píng)估、監(jiān)測(cè)設(shè)計(jì)與防治方面[3-8]以及水質(zhì)評(píng)價(jià)及趨勢(shì)[1-2]的研究，然而尚未對(duì)廣花盆地地下水可開采量進(jìn)行計(jì)算及預(yù)測(cè)，不利于實(shí)現(xiàn)地下水資源可持續(xù)開發(fā)利用。

地下水系統(tǒng)的狀態(tài)隨著水文氣象、開采強(qiáng)度、開采地段和開采層位等因素的變化而變化，因此確定地下水合理開采量，關(guān)鍵是要選擇能夠真實(shí)客觀地反映近期地下水資源狀況的開采量。國(guó)內(nèi)諸多學(xué)者使用數(shù)值模擬的方法計(jì)算不同研究區(qū)的地下水可開采量，并通過模擬研究區(qū)按照設(shè)計(jì)開采量進(jìn)行開采時(shí)地下水水位的變動(dòng)規(guī)律、降落漏斗的發(fā)展，論證設(shè)計(jì)開采量的可行性[9-13]。本文使用數(shù)值模擬及開采系數(shù)法計(jì)算廣花盆地的可開采量，兩種方法相互驗(yàn)證，論證設(shè)計(jì)開采量的合理性。本文可為實(shí)施地表水-地下水聯(lián)合開發(fā)提供基礎(chǔ)技術(shù)支持，對(duì)實(shí)現(xiàn)地下水高效利用與可持續(xù)發(fā)展、保障廣州市供水安全、改善環(huán)境地質(zhì)問題具有重要意義。

2 研究區(qū)概況

廣花盆地的地理范圍為北緯23°00′～23°30′，東經(jīng) 112°57′～113°26′，北高南低，面積約858 km2，以平原為主，臺(tái)地和丘陵次之。地下水可分為松散巖類孔隙水、巖溶水、基巖裂隙水3大類。

大氣降雨是研究區(qū)地下水補(bǔ)給的主要來源。平原區(qū)除主要接受大氣降雨的入滲補(bǔ)給外，還接受地表水入滲、灌溉入滲、潮水頂托反補(bǔ)給、基巖山區(qū)及丘陵臺(tái)地地下水的側(cè)向補(bǔ)給?；鶐r山區(qū)及丘陵臺(tái)地的地下水主要接受大氣降雨的入滲補(bǔ)給，山塘、水庫(kù)及渠道等對(duì)其入滲補(bǔ)給的規(guī)模較小。

總體上研究區(qū)地下水從北向南流動(dòng)。外圍丘陵山區(qū)一般地勢(shì)較高，淺部基巖裂隙水獲得補(bǔ)給后經(jīng)過短暫的徑流，便以泉或滲流形式排入附近的溪流；深層地下水則通過斷層和裂隙向廣花盆地匯流。

丘陵山區(qū)基巖裂隙水多以泉或滲流形式向溝谷排泄。江村-龍歸一線以北由于地下水位高于河水位，孔隙水和巖溶水主要向河流排泄。江村-龍歸一線以南地區(qū)地勢(shì)平坦，河涌較發(fā)育，豐水期由于地下水位埋深淺，地下水受潮水頂托影響，使地表水與地下水之間呈一定互補(bǔ)關(guān)系；枯水期地下水位埋深較深，地下水接受河涌水的入滲補(bǔ)給。此外，地下水的排泄方式還有開采和地表蒸發(fā)等。

21世紀(jì)60至80年代期間，廣花盆地城鎮(zhèn)供水地下水所占比例很高，農(nóng)村用水及灌溉用水大量使用地下水，采石場(chǎng)抽取地下水特別突出，建立了江村、肖崗、石井、新市、新華、雅瑤等多個(gè)水源地。

3 概念模型

3.1 模型邊界

第四系含水層北部和東部邊界為平原和山區(qū)分界線，山前第四系的孔隙介質(zhì)富水性好，容易接受山區(qū)降水潛流形式的側(cè)向補(bǔ)給，模型中將該邊界設(shè)為給定流量邊界；研究區(qū)西部邊界是蘆苞涌和白坭河，河流與潛水有較好的水力聯(lián)系，模型中將該邊界設(shè)為給定水頭邊界；地勢(shì)北高南低，通過地下水流場(chǎng)判斷南部為流量流出邊界，其他實(shí)測(cè)資料缺乏的邊界概化為通用水頭邊界。

巖溶含水層是模擬區(qū)主要的集中開采層，補(bǔ)給主要通過第一含水層越流補(bǔ)給，小塊巖溶出露地區(qū)接受降水補(bǔ)給，南部邊界為地下水流出邊界，不同時(shí)期地下水等水位線表明，巖溶層在不同年份和季節(jié)，存在著地下水的流出，由于廣花盆地三面被斷裂包圍，而且兩側(cè)多為基巖，或兩側(cè)大部分地層不銜接，設(shè)定為隔水邊界。模擬范圍及地下水水位觀測(cè)井分布情況見圖1。

潛水的自由水面為系統(tǒng)的頂部邊界，地下水主要以垂直運(yùn)動(dòng)為主，第四系和淺層巖溶地下水水力聯(lián)系比較密切。鉆孔揭露巖溶層底板以下巖溶水流動(dòng)緩慢，幾乎沒有地下水的開采，因此將鉆孔打到的位置作為本次數(shù)值模型的底部邊界。

3.2 模型結(jié)構(gòu)

依據(jù)研究區(qū)內(nèi)水文地質(zhì)勘察成果(主要為開采井及觀測(cè)井的鉆探成果)將整個(gè)含水系統(tǒng)自上而下分為3層，第1層是潛水含水層，主要的巖性為素填土及細(xì)砂等，其底板標(biāo)高在-30～35 m之間；第2層為相對(duì)隔水層，主要巖性為粉質(zhì)黏土，其底板標(biāo)高在-50～25 m之間；第3層為巖溶含水層，主要巖性為碳酸鹽巖，其底板標(biāo)高在-100～- 50 m之間。

圖1 模擬范圍及觀測(cè)井分布圖

3.3 水文地質(zhì)參數(shù)

潛水層水文地質(zhì)參數(shù)按照巖性的不同共劃分為8個(gè)參數(shù)區(qū)，在河流及其支流的兩側(cè)以及沖洪積扇中上部由于巖性主要是粗砂和礫石，故滲透系數(shù)K取值范圍為4～50 m/d，河道中下游逐漸過渡到細(xì)砂，中下部地區(qū)K取值較小，給水度取值范圍為0.02～0.2。第四系弱透水層水文地質(zhì)參數(shù)按照巖性的不同共劃分為6個(gè)參數(shù)分區(qū)，滲透系數(shù)取值范圍為0.014～3 m/d，彈性釋水系數(shù)取值范圍為1×10-5～3×10-4。巖溶含水層水文地質(zhì)參數(shù)根據(jù)抽水試驗(yàn)的結(jié)果，用GIS建立泰森多邊形，其滲透系數(shù)取值范圍為0.03～123 m/d，彈性釋水系數(shù)取值范圍為1×10-5～3×10-4。

3.4 源匯項(xiàng)

研究區(qū)的補(bǔ)給項(xiàng)主要包括降水入滲補(bǔ)給、河流與水庫(kù)入滲補(bǔ)給和灌溉入滲補(bǔ)給。降水入滲補(bǔ)給是主要的補(bǔ)給來源，河流補(bǔ)給以及灌溉入滲補(bǔ)給較小。排泄項(xiàng)主要有泉、地下水開采、蒸發(fā)、側(cè)向排泄，其中地下水開采是主要的排泄方式，根據(jù)研究區(qū)內(nèi)的工廠和農(nóng)田布局情況，結(jié)合計(jì)算需要，做如下概化：工業(yè)區(qū)主要集中在盆地南部區(qū)域，主要開采巖溶承壓水，因此均勻布設(shè)部分承壓水開采井和一定數(shù)量的潛水開采量。

4 數(shù)學(xué)模型

4.1 數(shù)學(xué)模型的建立

巖溶發(fā)育的不連續(xù)性、非均勻性和各向異性，加之實(shí)際勘探資料的制約，使得巖溶地下水建模成為巖溶水研究的難點(diǎn)和熱點(diǎn)。廣花盆地的含水巖組主要是可溶性的碳酸鹽巖，水流運(yùn)動(dòng)理論上是不符合達(dá)西定律的，但是盆地的巖溶發(fā)育相對(duì)比較均勻，因此仍應(yīng)用等效多孔介質(zhì)數(shù)值模擬模型分析廣花盆地的地下水流運(yùn)動(dòng)。該方法基于達(dá)西流假設(shè)和典型單元體的概念，用等效水力參數(shù)表征含水介質(zhì)滲透特性。大量的實(shí)例表明，以巖溶地下水系統(tǒng)水資源調(diào)查評(píng)價(jià)為目的，利用等效介質(zhì)方法模擬巖溶裂隙含水介質(zhì)是有效可行的。MODFLOW、FEFLOW、GMS等基于連續(xù)孔隙介質(zhì)的模擬軟件，已有很多成功應(yīng)用案例。

從空間上看，地下水流整體上以淺層垂向運(yùn)動(dòng)與深層水平運(yùn)動(dòng)相交替為特征流動(dòng)，為了用模型分析不同開采條件下含水層之間的互動(dòng)，在水文地質(zhì)概念模型基礎(chǔ)上，建立潛水和承壓水水流數(shù)值模型。本研究采用地下水模擬系統(tǒng)GMS。對(duì)研究區(qū)進(jìn)行水平剖分，研究區(qū)被剖分為160行、160列，共計(jì)有4 462個(gè)有效單元。模擬期為2011年4月至2013年3月。根據(jù)觀測(cè)井的水位資料，確定本次數(shù)值模擬的初始流場(chǎng)。由于第四系和下伏巖溶層有密切的水力聯(lián)系，監(jiān)測(cè)資料也表明潛水水位和巖溶層的水位在前汛期的4月基本一致，故主要觀察巖溶地下水初始流場(chǎng)(圖2)。

圖2 巖溶含水層初始流場(chǎng)圖

4.2 模型的校正與驗(yàn)證

采用等效多孔介質(zhì)單一連續(xù)方法的關(guān)鍵是確定等效的水力參數(shù)。本文使用手工校準(zhǔn)(正向校準(zhǔn))進(jìn)行水流模型的校準(zhǔn)，即手工調(diào)節(jié)各參數(shù)，直到模擬結(jié)果(地下水流場(chǎng)、地下水水位的變化趨勢(shì)、邊界條件、水文地質(zhì)參數(shù)和含水層結(jié)構(gòu))與實(shí)地觀測(cè)值基本一致[14-16]。反復(fù)調(diào)參后得到比較理想的模型校正結(jié)果。模擬時(shí)間總共730 d，分別計(jì)算第60、180、540和730 d的標(biāo)準(zhǔn)化殘差均方根，分別為4.996%、4.043%、6.517%、6.787%。將第60、180、540和730 d的不同觀測(cè)井的地下水水位計(jì)算值與觀測(cè)值進(jìn)行對(duì)比(圖3)，計(jì)算值與觀測(cè)值總體分布在y=x線兩側(cè)，表明觀測(cè)井的觀測(cè)值與計(jì)算值相差較小，地下水位的擬合程度較高。圖4為地下水水位等值線擬合圖，從圖4可以看出，所建立模型總體反映了模擬區(qū)的實(shí)際情況，模型經(jīng)校準(zhǔn)得到的地下水流場(chǎng)與實(shí)際流場(chǎng)基本上一致，地下水水位模擬值與實(shí)測(cè)值的動(dòng)態(tài)變化過程基本吻合。研究區(qū)地下水流場(chǎng)分布與地形趨勢(shì)一致，自東北向西南，地下水位呈遞減趨勢(shì)。3口代表性觀測(cè)孔(A1-12、A1-47和A1-23)的水位擬合過程線分別見圖5～7，觀測(cè)值總體分布在計(jì)算值附近，擬合效果較好。因此，所建立模型總體反映了研究區(qū)的實(shí)際情況，可用于地下水資源計(jì)算與預(yù)測(cè)。

圖3地下水水位計(jì)算值與觀測(cè)值對(duì)比圖4地下水水位等值線擬合圖(單位：m)

圖5 A1-12觀測(cè)井?dāng)M合曲線圖6 A1- 47觀測(cè)井?dāng)M合曲線圖7 A1-23觀測(cè)井?dāng)M合曲線

5 可開采量計(jì)算結(jié)果

5.1 數(shù)值模型法

盆地內(nèi)巖溶含水層大部分為埋藏型，上覆厚度約5～25 m的第四系松散沉積物，屬于巖溶塌陷易發(fā)區(qū)。地下水的運(yùn)動(dòng)會(huì)使覆蓋層土體的力學(xué)性質(zhì)發(fā)生改變，一般認(rèn)為抽取地下水引起塌陷的水位約束值與抽水井所在地區(qū)的土層厚度相對(duì)，即水位降深以不超過巖溶含水層頂板為宜。廣花盆地地下水水位至少高于第一含水層底板5 m，因此取降深5 m作為巖溶含水層抽取水量的約束條件。《廣東省地下水保護(hù)與利用規(guī)劃》表明，在巖溶覆蓋區(qū)由于大幅度抽取地下水易引起地面塌陷，而控制開采水位降深則是避免大面積、大幅度塌陷的主要措施，開采水位降深宜控制在5～8 m 以內(nèi)。因此本文確定巖溶開采水位降深最大值為5 m。

可開采量是在不大于該開采強(qiáng)度條件下，預(yù)測(cè)時(shí)段內(nèi)的地下水平均水位能保持在安全開采水位降深之內(nèi)。以此為基礎(chǔ)，以降深為5m的水位約束條件對(duì)模擬區(qū)未來5年(2015年9月至2020年9月)地下水可開采量進(jìn)行模擬預(yù)測(cè)。

在規(guī)劃的開采條件下模擬末期研究區(qū)含水層流場(chǎng)結(jié)果(圖8)顯示，2020年廣花盆地中心城區(qū)出現(xiàn)明顯的降落漏斗區(qū)域，降落漏斗中心的水位均控制在給定的降深內(nèi)，并且漏斗的中心是廣花盆地開采井最密集的區(qū)域，由此可見供水水源地集中開采地下水是引起地下水水位降落漏斗空間位置變化的主要影響因素。圖9為模擬期3個(gè)觀測(cè)井水位過程線，由圖9可知，中心城區(qū)附近的觀測(cè)井(A1-22和A1-25)的水位隨著時(shí)間逐漸降低，但是降低的幅度在減小，最大降深約為4 m，偏遠(yuǎn)地區(qū)的監(jiān)測(cè)井(A1-40)的水位變化幅度不大。

圖8 研究區(qū)2020年9月模擬流場(chǎng)圖

中心城區(qū)降深最大的開采井的水位預(yù)測(cè)如圖10所示。隨著開采井運(yùn)行，水位開始下降，在3月份水位達(dá)到最小值，這時(shí)控制水位降深最小，為0.3 m，從4月開始地下水水位開始回升，控制水位降深也逐漸增大，直到9月份達(dá)到峰值。10月至次年3月，地下水水位又持續(xù)下降，之后又上升，如此循環(huán)，2016年9月水位降深達(dá)到最大值2.36 m。通過不斷的調(diào)節(jié)抽水井的抽水強(qiáng)度，使預(yù)測(cè)時(shí)段內(nèi)的地下水平均降深能保持在安全開采水位降深之內(nèi)，模型計(jì)算得到的地下水允許開采量為12 997.78×104m3。

圖9 模擬期3個(gè)觀測(cè)井水位過程線

圖10 中心城區(qū)觀測(cè)井水位降深變化過程線

5.2 可開采系數(shù)法

廣花盆地地下水資源量采用補(bǔ)給量法計(jì)算。

降水入滲補(bǔ)給量的計(jì)算公式為：

Q降=F·α·P

(1)

式中：F為可入滲面面積，m2，將研究區(qū)離散為若干1 km2的單元格，參與計(jì)算的各單元格面積是指該單元格扣除遙感解譯得到的地表水系和硬化地面部分后的剩余面積；P為降水量，m，將降水量插值到研究區(qū)每一個(gè)單元格內(nèi)；α為降水入滲補(bǔ)給系數(shù)，將每一個(gè)自動(dòng)監(jiān)測(cè)井的年內(nèi)數(shù)據(jù)，繪制成表，選取次降水強(qiáng)度較大的若干場(chǎng)次降水，劃定對(duì)應(yīng)次降水時(shí)段的地下水水位增幅，根據(jù)公式(2)求取土壤的降水入滲補(bǔ)給系數(shù)。

α=μ· (∑ΔH)/P次

(2)

式中：μ為水位變動(dòng)帶含水層給水度；ΔH為次降水造成的水位增幅，m；P次為引起該水位增幅對(duì)應(yīng)時(shí)段的次降水量，m。

灌溉入滲補(bǔ)給量的計(jì)算公式為：

Q灌=15a0·m·fg

(3)

式中：a0為灌溉入滲補(bǔ)給系數(shù)，農(nóng)作物a0值為0.157，園林作物a0值為0.12；m為廣州市灌溉定額量，m3/hm2；fg為灌溉面積，hm2。

平原區(qū)的山前側(cè)向補(bǔ)給量計(jì)算公式為：

Q側(cè)補(bǔ)=K·I·L·B

(4)

式中：Q側(cè)補(bǔ)為山前側(cè)向流出量，m3/d；K為剖面位置的滲透系數(shù)，m/d；I為垂直于剖面的水力坡度；L為山丘平原界線的長(zhǎng)度，m；B為剖面含水層厚度，m。

補(bǔ)給量采用達(dá)西公式計(jì)算：

Q河補(bǔ)=K·I·L·B

(5)

式中：Q河補(bǔ)為單側(cè)河道側(cè)向補(bǔ)給量，m3/d；K為剖面位置的滲透系數(shù)，m/d，通過鉆孔抽水試驗(yàn)成果插值確定；L為河段長(zhǎng)度，m；B為計(jì)算剖面厚度，m，通過鉆孔抽水試驗(yàn)成果插值確定；I為垂直于剖面的水力坡度。

經(jīng)計(jì)算，盆地內(nèi)白坭河、新街河、天馬河的側(cè)向補(bǔ)給量較大，分別為1 075.31×104、579.13×104、279.43×104m3。此外，盆地內(nèi)流溪河段、珠江(西航道)段側(cè)向補(bǔ)給量分別為136.82×104、6.76×104m3。區(qū)內(nèi)河道側(cè)向補(bǔ)給量共計(jì)2 077.45×104m3。

利用補(bǔ)給系數(shù)法計(jì)算湖、庫(kù)、渠的入滲補(bǔ)給不僅需要確定入滲面積，還需判別其水位與地下水水位之間的大小關(guān)系。計(jì)算方法為：

Q湖庫(kù)渠=F·I·ψ，I>0

(6)

式中：Q湖庫(kù)渠為湖、庫(kù)、渠及微小支流對(duì)地下水的補(bǔ)給，m3/a；F為計(jì)算單元格內(nèi)納入計(jì)算的湖、庫(kù)、渠面積，m2；I為計(jì)算單元格與周邊含水層之間的平均水力坡度；ψ為單位水力坡度的滲漏補(bǔ)給系數(shù)，m/a，本次計(jì)算中將區(qū)內(nèi)水系劃分為主要河道和湖、庫(kù)、渠及微小支流，并分別計(jì)算其與地下水之間的補(bǔ)排量，事實(shí)上這種劃分并沒有嚴(yán)格的界定方法，人為主觀因素明顯。因此，對(duì)于河道側(cè)向補(bǔ)給量，依據(jù)公式(5)與公式(6)的計(jì)算結(jié)果應(yīng)較為一致。通過試算，當(dāng)ψ=1.1 m/a時(shí)，兩者結(jié)果基本一致。由此確定單位水力坡度的滲漏補(bǔ)給系數(shù)為1.1 m/a。

經(jīng)計(jì)算，廣花盆地各補(bǔ)給項(xiàng)結(jié)果見表1。表1中地下水資源量的計(jì)算公式為：

Q平=Q降+Q灌+Q河+Q湖庫(kù)渠+Q山前

(7)

式中：Q平為平原區(qū)的地下水資源量，104m3；Q降為降水入滲補(bǔ)給量，104m3；Q灌為灌溉入滲補(bǔ)給量，104m3；Q河為河道側(cè)向補(bǔ)給量，104m3；Q湖庫(kù)渠為湖、庫(kù)、渠滲漏補(bǔ)給量，104m3；Q山前為山前側(cè)向補(bǔ)給量，104m3。則廣花盆地地下水資源量為18 972.99×104m3。

表1 廣花盆地地下水資源量計(jì)算結(jié)果 104 m3

本文根據(jù)1∶100000比例尺區(qū)域水文地質(zhì)圖資料確定廣花盆地各計(jì)算單元的富水性，按各級(jí)富水性相對(duì)應(yīng)的可開采系數(shù)范圍，結(jié)合開采條件、抽水試驗(yàn)成果及實(shí)際開采量，確定可開采系數(shù)如表2所示。經(jīng)計(jì)算，廣花盆地地下水可開采量為14015.09×104m3。

表2 廣花盆地可開采系數(shù)表

6 結(jié)論與建議

6.1 結(jié) 論

(1)本文運(yùn)用水文地質(zhì)知識(shí)分析數(shù)據(jù)資料，按照規(guī)范化的建模流程(GMS)構(gòu)建了廣州市廣花盆地的地下水水流模型。結(jié)果表明，該模型可為廣花盆地提供合理的地下水水流模擬結(jié)果，能揭示地下水流動(dòng)規(guī)律，可用于地下水資源計(jì)算與預(yù)測(cè)。

(2)采用補(bǔ)給量法計(jì)算廣花盆地的地下水資源量為18 972.99×104m3，結(jié)合廣州市地下水基礎(chǔ)資料條件，考慮開采條件、實(shí)際開采現(xiàn)狀和地下水資源狀況等因素，采用可開采系數(shù)法計(jì)算廣花盆地地下水可開采量為14 015.10×104m3；以可開采水位最大降深5 m作為巖溶含水層抽取水量的約束條件，采用數(shù)值模型計(jì)算的地下水可開采量為12 997.78×104m3，兩種方法計(jì)算的結(jié)果相對(duì)誤差為7.3%，兩種方法相互驗(yàn)證，表明本文可開采量計(jì)算結(jié)果較為可靠。

6.2 建 議

(1)結(jié)合已有的地下水工作成果及廣花盆地水文地質(zhì)條件，進(jìn)行地下水超采區(qū)劃定及評(píng)價(jià)。

(2)提出廣花盆地應(yīng)急水源地的開采控制目標(biāo)及措施：

廣花盆地巖溶地下水區(qū)一般情況下嚴(yán)禁開采、嚴(yán)格保護(hù)，但在啟動(dòng)應(yīng)急預(yù)案或?qū)９┧黠@不足時(shí)，并經(jīng)論證不會(huì)產(chǎn)生環(huán)境地質(zhì)問題的前提條件下，可適量開采。該目標(biāo)與《廣東省地下水功能區(qū)劃》關(guān)于應(yīng)急水源地管理目標(biāo)相一致。

除應(yīng)急備用水、礦泉水、地?zé)崴葘?duì)水溫、水質(zhì)有特殊要求的地下水外，不再批準(zhǔn)抽取地下水取水許可申請(qǐng)。水行政主管部門對(duì)批準(zhǔn)的地?zé)崴?、礦泉水取水工程應(yīng)核定地下水開采量和年度用水計(jì)劃；對(duì)已批準(zhǔn)的地下水為自備水源的取水戶，要實(shí)行計(jì)劃用水，采取相關(guān)節(jié)水措施，逐年核減地下水開采量。

對(duì)于在城市供水管網(wǎng)覆蓋區(qū)域內(nèi)，新建、改建或擴(kuò)建工程項(xiàng)目自備取水設(shè)施取用地下水的申請(qǐng)不予批準(zhǔn)。市區(qū)各自備水源單位原持有取水許可證到期后，城市公共供水管網(wǎng)能滿足用水需求的，不再辦理續(xù)證手續(xù)，不得再開采地下水。在公共供水管網(wǎng)未覆蓋區(qū)域內(nèi)，建設(shè)單位應(yīng)經(jīng)建設(shè)項(xiàng)目水資源論證、取水許可獲得批準(zhǔn)后，方可興建地下水取水工程或者設(shè)施。

(3)加強(qiáng)地下水動(dòng)態(tài)觀測(cè)，逐步構(gòu)建地下水水位、水溫和水質(zhì)的自動(dòng)監(jiān)測(cè)、信息共享及預(yù)警預(yù)報(bào)體系。