Modflow在平原區(qū)深層地下水資源評價中的應(yīng)用

李華，顏坤，侍建國

(天津市水文水資源勘測管理中心，天津 300061)

0 引言

近年來，隨著計算機(jī)技術(shù)的日益發(fā)展，國內(nèi)的地下水模擬技術(shù)取得了較大的進(jìn)步，但缺乏大區(qū)域長時間序列的高精度地下水?dāng)?shù)值模型評價地下水資源的相關(guān)研究。該文以天津市為例，全面介紹了利用Modflow建立大區(qū)域長時間序列數(shù)值模型的過程，為相關(guān)技術(shù)研究提供技術(shù)支撐，為深層地下水資源評價提供科學(xué)依據(jù)。

在一次海洋基座式起重機(jī)回轉(zhuǎn)平臺物資質(zhì)量監(jiān)造中，勝利油田的設(shè)備監(jiān)造人員發(fā)現(xiàn)龍骨之間和封板角焊縫焊腳高度與圖紙要求存在1毫米偏差后，立即要求廠家對其焊縫打磨后重新焊接，否則不允許進(jìn)行下道工序。

天津市淺層地下水研究程度較高，而深層地下水研究程度相對較低，曾開展過2次區(qū)域水資源綜合評價，第一次為1956—1984年，評價內(nèi)容包括淺層地下水資源量、可開采量和深層地下水可開采量；第二次為1980—2000年，評價內(nèi)容包括淺層地下水資源量和可開采量。兩次評價的重點集中在淺層地下水資源，均未涉及深層地下水資源的構(gòu)成評價，影響了深層地下水資源的合理開發(fā)利用。

最后，必須要定期的對學(xué)生進(jìn)行階段性考核，準(zhǔn)確有效地了解學(xué)生的學(xué)習(xí)進(jìn)展情況，然后依照學(xué)生的學(xué)習(xí)進(jìn)展對教學(xué)目標(biāo)不斷的進(jìn)行調(diào)整，以此來保證教學(xué)目標(biāo)的適用性，進(jìn)而保證高中化學(xué)分層教學(xué)模式能夠有效地發(fā)揮出應(yīng)有的作用。

傳統(tǒng)地下水資源的評價一般只評價地下水的資源量和可開采量，難以進(jìn)行不同地下水開發(fā)利用方案下的地下水不同要素的預(yù)測和分析。采用新技術(shù)和手段進(jìn)行水資源調(diào)查評價，全面分析深層地下水資源構(gòu)成，揭示其戰(zhàn)略資源屬性特征，評價可開采資源量，提高地下水資源管理效率，提高地下水資源管理水平成為亟待解決的問題。

歷史模型初始水位采用1980年3月底水位值，深層地下水水位監(jiān)測井較少，水位監(jiān)測成果難以滿足模擬需求，同時此時間段內(nèi)深層地下水開采量比較穩(wěn)定，水位變幅相對較小，因此在1980年3月底深層地下水長觀孔監(jiān)測水位成果的基礎(chǔ)上，加入地下水統(tǒng)測孔水位數(shù)值。利用克里金插值法[1]，得到平原區(qū)第Ⅱ～Ⅵ承壓含水組的初始水位值，針對東南部區(qū)域第Ⅳ和第Ⅴ承壓含水組水位觀測孔更少的情況，利用已有的鉆孔標(biāo)高進(jìn)行修正，結(jié)合水位外推，確定Ⅳ和第Ⅴ承壓含水組初始水位。第Ⅵ承壓含水組及以下地下水無開采量，水位動態(tài)變化趨勢很小，采用1997年12月水位值作為初始水位，采用1997年12月31日水位值作為現(xiàn)狀模型初始水位。

1 深層地下水流—地面沉降耦合模型的建立

1.1 時空離散

1.1.1 空間剖分

原始洞窟壁畫中就有對動物油脂、礦物質(zhì)粉末、巖石、木炭等材料的使用，而文藝復(fù)興時期藝術(shù)家所使用的蛋清、石灰、亞麻仁油等，在其所處的時代也無疑是一種“新鮮的材料”。材料使用的變化一直在促使藝術(shù)概念的變化，因為每種材料因其與生俱來的物性不同，運(yùn)用的創(chuàng)作方法不同，必然導(dǎo)致其背后所蘊(yùn)含的信息也就不同。

紅色文化在高校校園文化中占比不足，紅色教育活動開展較少，教育形式也局限于征文比賽、演講比賽、合唱比賽、觀影活動等。受時間、經(jīng)費等因素的限制，去烈士陵園掃墓、參觀紅色教育基地等活動機(jī)會更是少之又少，導(dǎo)致學(xué)生接受紅色教育途徑少，校園紅色教育實效性差。除了線下活動，據(jù)調(diào)查，各大高校少有建立專門的紅色教育網(wǎng)站或在官網(wǎng)上有紅色專題宣傳板塊，即便已建立的紅色教育網(wǎng)站也存在內(nèi)容陳舊、更新慢、吸引力不強(qiáng)等問題，點擊量少，不能充分發(fā)揮網(wǎng)絡(luò)紅色教育作用，教育體系不夠完善[4]。

研究區(qū)域為天津市整個平原區(qū)，面積約1.07×104km2?？臻g剖分采用有限差分法，使用GMS軟件中的Modflow模塊進(jìn)行矩形剖分，單元格500m×500m，共剖分7層346行248列，活動單元格總數(shù)276111個，其中II層活動單元格42628個，Ⅲ～Ⅵ層每層活動單元格38171個(Ⅵ層以下概化為Ⅶ層，作為越流邊界層)。

1.1.2 時間離散

模擬期為1980—2013年，考慮到模擬期長達(dá)34年，模擬初始時間段和截止時間段內(nèi)資料的詳細(xì)程度有較大差異，為了便于調(diào)整驗證模型，將1980—2013年分為2個模擬期。1980年3月—1998年3月，由于歷史資料相對缺乏，以每3個月為一個應(yīng)力期，共72個應(yīng)力期，稱為歷史模型；1998年1月—2013年12月，資料較為豐富，以1個月為應(yīng)力期，共192個應(yīng)力期，稱為現(xiàn)狀模型。

1.2 數(shù)據(jù)處理

1.2.1 定解條件處理

(1)初始水位處理

秧床整地作畦要求上細(xì)下粗、上實下松。在播種前施肥并耕翻 1～2次，耕翻應(yīng)耕深15cm以上，做成畦寬1.2～1.4m、溝寬50cm、溝深20～30cm的畦，苗床四周開好圍溝，以保證排水通暢。畦整平后用2.0kg壯秧劑加干細(xì)土20kg充分混拌可撒施25m2苗床，再反復(fù)耙勻于0～5cm深的土層內(nèi)，用木板壓平，澆透水即可播種。

(2)邊界條件處理

這種理論在比較法上的支持是德國的實踐。早在1902年，德國帝國法院即在“Duotal Lagend”案中確認(rèn)了對于根據(jù)專利方法(一種生產(chǎn)制造專利方法)直接獲得的產(chǎn)品，在權(quán)利人投入市場后，在其上的專利保護(hù)的效力就用盡了[注]Entscheidungen des Reichsgerichts in Zivilsachen, Erster Band, Der ganzen Reihe einundfünfzigster Band., Verlag von Beit und Comp. 1903, 139.。

1.2.5 地下水開采量數(shù)據(jù)的處理

系統(tǒng)上邊界確定為潛水含水層自由水面，潛水通過系統(tǒng)上邊界與周邊環(huán)境產(chǎn)生垂向上的越流交換，各含水組之間也存在垂向上的越流交換。底界與Ⅵ含水組有水量交換，定義為水位邊界。

粘性土層概化的合理性直接關(guān)系地面沉降計算結(jié)果。將天津市平原區(qū)粘性土層厚度小于或者等于1.5m的定義為非滯后粘性土層，力學(xué)參數(shù)包括非彈性骨架儲水系數(shù)和彈性骨架儲水系數(shù)。粘性土層厚度大于1.5m的定義為滯后性粘性土層，力學(xué)參數(shù)包括垂向滲透系數(shù)、非彈性骨架儲水系數(shù)、彈性骨架儲水系數(shù)、粘土土層的有效厚度及比例因子。將非滯后性粘性土層和不同承壓含水層組的砂層均概化為非滯后壓縮層，在同一承壓含水層組中的滯后性粘性土層處理為滯后性壓縮層。

將初始沉降量設(shè)為0來滿足模型的擬合需要?？紤]到粘性土的彈性以及塑性變形，引入前期固結(jié)壓力概念。當(dāng)前期最低水位高于當(dāng)前地下水水位時，粘性土體將產(chǎn)生塑性變形變化，當(dāng)前期最低水位低于當(dāng)前地下水水位時，粘性土體將產(chǎn)生彈性變形變化，通過地下水水位觀測孔監(jiān)測到的水位值得到各承壓含水組的最低水位。

針對偏遠(yuǎn)、交通不便的村民辦證和投訴不方便等情況，依托在鄉(xiāng)鎮(zhèn)、村委設(shè)立“便民服務(wù)點”，發(fā)放《便民服務(wù)聯(lián)系卡》，實行上門服務(wù)。

1.2.2 源匯項處理

利用GIS軟件的空間分析，利用程序?qū)懭氲较鄬?yīng)的Modflow格式文件中，詳細(xì)情況見表1。

雜志里講的正是一個中年男人的情感故事，他邊看邊和自己做著比較，正看得入神，門忽然開了，走進(jìn)幾個身穿制服的人，夾著公文包，面色冷峻，其中一個說道：我們是專案組的，現(xiàn)在決定對你實行刑事拘留。兩個年輕力壯的人遂上前把他一夾，兩只手腕間多了一副錚亮的手銬，閃著耀眼的白光。他在白光中睜開眼，天色大亮，原來是一場夢。他渾身濕漉漉的，用手胡亂抹著汗水。妻子那邊的被子掀開，衣服和挎包都不見了，枕頭下也沒有手機(jī)。回想夢中的情節(jié)，尹愛群的脊背又是陣陣寒氣。

表1 源匯項處理對照

1980—2013年平均排泄量為6.90×108m3/a，開采量為6.28×108m3/a，占排泄量的90%～93%；側(cè)向流出量為0.31×108m3/a，占排泄量的4.44%；越流量為0.13×108m3/a，占排泄量的1.90%；壓縮釋水量為0.19×108m3/a，占排泄量的2.72%(表2)。

依據(jù)平原區(qū)各項水文地質(zhì)參數(shù)的不同分區(qū)，將降水入滲系數(shù)、各承壓含水組滲透系數(shù)、田間灌溉入滲補(bǔ)給系數(shù)、潛水含水層給水度和承壓含水組儲水率等不同水文地質(zhì)參數(shù)的分區(qū)初始值輸入到建立的模型中，通過擬合地下水不同承壓含水層組典型水位觀測孔和流場動態(tài)變化過程線，識別不同承壓含水組的水文地質(zhì)參數(shù)，確定各水文地質(zhì)參數(shù)的分區(qū)數(shù)值。

1.2.4 粘性土層概化及壓縮參數(shù)處理

(3)沉降的定解條件處理

寶坻斷裂南部潛水含水層多含咸水，基本處于不開采狀態(tài)，地下水水位年內(nèi)變化幅度較?。粩嗔驯辈繚撍畬犹幱诘叵滤到y(tǒng)的補(bǔ)給區(qū)，水位恢復(fù)較快且變幅較小，故不考慮潛水層的壓縮釋水。上層含水層的越流補(bǔ)給是各含水組的主要補(bǔ)給來源，如地下水開采量較大，則水位年內(nèi)變幅也較大，且呈現(xiàn)多年持續(xù)下降趨勢，因此將Ⅱ～Ⅵ含水組作為壓縮釋水層。

根據(jù)不同含水層的砂層及厚度、粘性土層的層數(shù)及厚度，計算出每一個地下水水位觀測孔的滯后壓縮層和非滯后壓縮層對應(yīng)的參數(shù)，利用克里金插值法，得到每一個模型網(wǎng)格的參數(shù)數(shù)據(jù)，應(yīng)用VBA生成Modflow2005中SUB模塊所需要的.sub文件。

北部山區(qū)和平原區(qū)分界線、環(huán)渤海海岸帶和寶坻斷裂為自然邊界，其余均為人為邊界。流量邊界確定為山區(qū)和平原區(qū)分界線，渤海灣海岸帶以水頭值為0的水頭邊界處理，平原區(qū)北部的寶坻斷裂上下部有水量交換關(guān)系，并非嚴(yán)格意義上的不透水邊界，和其他人為邊界統(tǒng)一處理為三類邊界。

建模采用的核心模塊為Modflow-2005，但在處理大區(qū)域模型時，存在運(yùn)行時間長、開采量只能處理成點狀等缺點。通過改進(jìn)RCH子程序包，開發(fā)出新的子程序包RAW(rechargeand well，RAW)，實現(xiàn)了垂向上對多層單元格補(bǔ)排量的表達(dá)，有效壓縮了源匯項文件，減少了大區(qū)域數(shù)值模型所需要讀取的數(shù)據(jù)量，提高了模型模擬的工作效率。

1.3 模型識別驗證

模型模擬的地面沉降[10-19]與實測的地面沉降增長趨勢基本一致，其中塘沽分層標(biāo)、西青分層標(biāo)和軍糧城分層標(biāo)監(jiān)測點地面沉降模擬值與實測值相差較小，大直沽天津堿廠分層標(biāo)模擬值與實測值相差較大，總體均表現(xiàn)為地面沉降值為增長的趨勢(圖5)。

1.3.1 地下水流場

模型擬合校正過程要考慮2個重要過程，一是校正不同水文地質(zhì)參數(shù)，二是推算壓縮釋水量和反演地下水開采量。通過不斷調(diào)整水文地質(zhì)參數(shù)，進(jìn)行地下水流場[7-9]、觀測孔過程線與實測數(shù)據(jù)的擬合對比。識別后的流場擬合整體效果較好，各含水組實測水位和模擬水位總體變化趨勢一致，單觀測井實測水位和模擬水位變化趨勢也基本一致(圖1—圖4)。

圖1 1997年12月Ⅳ含水組流場擬合

圖2 1997年12月Ⅴ含水組流場擬合

圖3 1980-1998年寶坻區(qū)BD18觀測井(Ⅱ含水組)水位過程線擬合

圖4 1998-2013年靜海區(qū)JH36觀測井(Ⅲ含水組)水位過程線擬合

1.3.2 沉降量的擬合

對建立的大區(qū)域、多層結(jié)構(gòu)、長時間序列地下水流-地面沉降數(shù)值模型，應(yīng)用試估-校正法進(jìn)行識別和檢驗[2-6]，總體原則如下：①模型識別的水文地質(zhì)參數(shù)均要符合實際的水文地質(zhì)條件；②模擬的地下水流場與實際的地下水流場要保持基本一致；③模擬的地下水均衡變化情況要基本符合研究區(qū)域?qū)嶋H的水均衡變化情況；④模擬的地下水水位動態(tài)變化趨勢與實測的水位動態(tài)變化趨勢基本相似。

圖5 2011—2012年地面沉降量擬合

2 深層地下水資源評價

2.1 地下水均衡項

補(bǔ)給項[20]包括通用水頭邊界流入量、給定水頭邊界流入量、山前側(cè)向補(bǔ)給量、越流補(bǔ)給量及壓縮釋水量，排泄項包括通用水頭邊界流出量、給定水頭邊界流出量、人工開采量、越流排泄量及壓縮釋水量。

2.2 地下水均衡分析

1980—2013年平均補(bǔ)給量為6.46×108m3/a，其中Ⅱ和Ⅶ含水組越流補(bǔ)給量為2.72×108m3/a，占補(bǔ)給量的33.82%；粘性土壓縮釋水量為1.48×108m3/a，占補(bǔ)給量的18.36%；彈性釋水量為0.09×108m3/a，占補(bǔ)給量的1.09%；北部邊界及邊界的側(cè)向補(bǔ)給量為2.17×108m3/a，占補(bǔ)給量的26.97%。層間越流、側(cè)向流入、壓縮釋水已經(jīng)成為深層地下水的主要補(bǔ)給來源。

1.2.3 水文地質(zhì)參數(shù)處理

總之，運(yùn)用多媒體輔助語文教學(xué)，可以說是有利有弊，關(guān)鍵在于教師的把握，在于教師的思想態(tài)度。用好多媒體可以促進(jìn)教學(xué)，提高教學(xué)效果，既有利于師，也有利于生；如果使用不當(dāng)多媒體，可能會導(dǎo)致教學(xué)效果差，學(xué)生反感，教師尷尬，先進(jìn)的設(shè)備變成害生的毒藥。

深層地下水均衡差為-0.44×108m3/a，深層地下水系統(tǒng)處于負(fù)均衡狀態(tài)。

表2 1980—2013年平均深層地下水資源組成(104m3/a)

2.3 可開采資源量評價

(1)評價原則

南水北調(diào)通水之前，天津深層地下水開采量遠(yuǎn)超過其獲得的補(bǔ)給量，從水均衡的角度將外界補(bǔ)給量作為深層地下水的可開采資源量，結(jié)果相對偏小，也不符合實際情況。因此，深層地下水可開采資源量的評價原則是：一方面要不產(chǎn)生嚴(yán)重環(huán)境地質(zhì)問題，另一方面要結(jié)合天津深層地下水開發(fā)利用現(xiàn)狀及歷史，人為給定環(huán)境約束條件。

(2)評價方法

主要考慮不同區(qū)域深層地下水的不同開發(fā)利用程度[21]。開發(fā)利用程度相對較高的區(qū)域，產(chǎn)生地質(zhì)災(zāi)害的可能性就越大，禁止利用深層地下水釋水量，可開采資源量為側(cè)向徑流補(bǔ)給量和越流補(bǔ)給量兩項之和。相反，開發(fā)利用程度相對較低的區(qū)域，產(chǎn)生地質(zhì)災(zāi)害的可能性就越小，可開采資源量在考慮側(cè)向徑流補(bǔ)給量和越流補(bǔ)給量的基礎(chǔ)上，結(jié)合現(xiàn)狀開采條件，要考慮可動用的釋水量。

深層地下水可利用總量為3.08×108m3/a，其中Ⅱ含水組為2.20×108m3/a，Ⅲ含水組為0.52×108m3/a，Ⅳ含水組為0.26×108m3/a，Ⅴ含水組為0.10×108m3/a。

3 南水北調(diào)通水后地下水流場和地面沉降變化趨勢研究

3.1 地下水關(guān)鍵控制水位

利用數(shù)值模擬進(jìn)行地下水資源評價，通常分2個階段進(jìn)行，一是解逆問題(也稱反演)，二是解正問題(也稱正演)。解逆問題是根據(jù)水位、開采量等資料進(jìn)行數(shù)值計算，求出水文地質(zhì)參數(shù)，驗證數(shù)值模型。解正問題是利用解逆問題的參數(shù)，計算不同開采方案時的水位值，進(jìn)行方案對比和水位預(yù)報。

通過建立研究區(qū)地下水流-地面沉降耦合數(shù)值模型，將地下水開采量作為約束條件代入模型，確定地下水關(guān)鍵性控制水位。以天津市為例，根據(jù)天津市地下水中長期規(guī)劃，在壓采條件下求解出不同目標(biāo)年下的關(guān)鍵控制水位，計算出2020年和2030年平均紅線水位分別為-42.44m和-37.86m，用數(shù)值模擬法確定關(guān)鍵控制水位，是一種全面考慮水文地質(zhì)因素，計算精度較高的方法。

3.2 地下水流場變化趨勢

截至2017年12月27日，天津市累計受水23.7億m3。隨著最嚴(yán)格水資源管理制度的實施，地下水壓采進(jìn)程加快，2016年，深層地下水開采量1.76億m3，比深層地下水可開采資源量少1.32億m3。由于大幅度壓采深層地下水，深層地下水各承壓含水組水位將呈現(xiàn)整體上升趨勢。深層地下水系統(tǒng)將處于正均衡狀態(tài)，不同含水組儲存變化量為正值，彈性水在含水組中儲存了起來；壓縮釋水凈補(bǔ)給量為負(fù)值，多余的地下水儲存到土骨架中，儲存量處于逐漸恢復(fù)過程中。

ZHANG Ming-ming, ZHAO Pei, YU Yue-qing, ZHANG Cui-gai, GAO Wei, LI Fang

東麗區(qū)和津南區(qū)地下水水位有先下降后上升的變化趨勢，由于深層地下水開采量仍稍大于各區(qū)可開采量，在初始時間段內(nèi)，地下水水位有小幅下降，隨著周邊地區(qū)水位的持續(xù)上升，各區(qū)之間的內(nèi)部地下水交換量逐漸增大，水位也將逐漸上升。由于下層含水組水量主要由上層含水組越流補(bǔ)給，因此與上層相比，下層含水組得到的越流補(bǔ)給量相對越小，下層含水組水位的恢復(fù)速率比上層含水組也要小。

6.聘請張掖市螢火蟲公益的成員為代理家長。此公益組織主要是包含著個體老板和記者等多個主體，經(jīng)過合理的協(xié)調(diào)組合而成，他們一般是運(yùn)用周末或者是空閑時間，給村里的孩子們輔導(dǎo)課程，為其開設(shè)講座等，讓孩子們擁有一個愉快的周末時光。對一些特別貧困的留守兒童進(jìn)行一對一資助，目前已經(jīng)有140多位留守兒童受到資助。

研究發(fā)現(xiàn)抑郁與服藥依從性密切相關(guān)[12]。腫瘤病人經(jīng)常出現(xiàn)情緒低落、悲觀、失望，特別是病情變化會加重病人恐懼、絕望心理，嚴(yán)重者甚至拒絕治療[13]。病人1：“我不看病了，看也白看，又治不好?！辈∪?：“我怎么會得這個病，治療這么長時間，現(xiàn)在還是渾身沒勁。”病人8：“有這個病1年多了，這一年不是吃這藥就是吃那藥，都成藥罐子了，一說吃藥就郁悶，就想發(fā)火”。

模型預(yù)測2018年各承壓含水組水位恢復(fù)較快，最大處水位上升4m左右，其余大部分地區(qū)水位上升2～3m。大部分區(qū)域地下水水位前期恢復(fù)較快，導(dǎo)致后期越流量減少，水位恢復(fù)呈現(xiàn)前快后慢，然后逐漸穩(wěn)定的趨勢。

第二，機(jī)械故障，導(dǎo)致這一故障出現(xiàn)的原因有很多，其中車輛的車門的相關(guān)部件由于擠壓等出現(xiàn)了變形，或者是在生產(chǎn)的過程中，尺寸的誤差比較大而使得硬件條件不符合相關(guān)要求是最重要的因素，同時在軌道車輛運(yùn)行的過程中，由于車輛車門的開關(guān)頻率比較高，從而使得車門硬件條件發(fā)生變化也是導(dǎo)致這一類型故障出現(xiàn)的主要原因。在客流量比較大的時候，乘客在上下車的時候都會擠擠攘攘，很多乘客會被擠在車門的位置，長此以往就會導(dǎo)致車門發(fā)生變形，這樣一來就會使得車門出現(xiàn)開關(guān)遲鈍等常見故障，不利于軌道車輛的正常運(yùn)行。

3.3 地面沉降變化趨勢

當(dāng)大量開發(fā)利用地下水時，孔隙水壓力降低，直接導(dǎo)致各含水組砂層壓密，當(dāng)孔隙水壓力恢復(fù)時，各含水組砂層基本能恢復(fù)原狀，砂礫類巖土呈現(xiàn)彈性變形的狀態(tài)。當(dāng)?shù)叵滤笠?guī)模超采，粘性土壓密后，其結(jié)構(gòu)將發(fā)生不可逆的變化，即使孔隙水壓力恢復(fù)，粘性土層也基本上繼續(xù)保持壓密狀態(tài)，以塑性變形為主。故而當(dāng)?shù)叵滤畨翰蓵r，各含水組砂層有回彈趨勢，粘性土層仍然處于壓密狀態(tài)。

在現(xiàn)有的壓采方案下模型預(yù)測平原區(qū)地面沉降沒有繼續(xù)加大的趨勢，但也沒有出現(xiàn)全面回彈的現(xiàn)象，截至2018年，除薊州和寧河區(qū)地面沉降速率在2mm/a之外，其余各區(qū)地面沉降處于涵養(yǎng)回彈態(tài)勢，地面沉降整體開始回彈，回彈平均速率為3.42mm/a。

4 結(jié)論

(1)在建立平原區(qū)長時間序列地下水流-地面沉降數(shù)值模型時，可根據(jù)數(shù)據(jù)精度不同，建立不同的連貫?zāi)Ｐ停瑪?shù)據(jù)精度低的時間段內(nèi)，模擬期可適當(dāng)延長，數(shù)據(jù)精度高的時間段內(nèi)，模擬期適當(dāng)縮短，前后模型的模擬期可重復(fù)，但不可缺失。

(2)用試估-校正法進(jìn)行識別和檢驗地下水?dāng)?shù)值模型時，水文地質(zhì)參數(shù)、地下水流場、地下水均衡變化和地下水水位動態(tài)變化的擬合情況要符合研究區(qū)實際的水文地質(zhì)條件。模型的擬合校正過程要考慮2個重要過程，一是校正不同水文地質(zhì)參數(shù)，二是推算壓縮釋水量和反演地下水開采量。

(3)用地下水?dāng)?shù)值模型進(jìn)行平原區(qū)深層地下水可開采資源量評價時，基本原則是既要不產(chǎn)生嚴(yán)重環(huán)境地質(zhì)問題，又要結(jié)合研究區(qū)深層地下水開發(fā)利用現(xiàn)狀及歷史，人為給定環(huán)境約束條件，評價方法要考慮研究區(qū)不同區(qū)域深層地下水不同開發(fā)利用程度，開發(fā)利用程度低的區(qū)域要考慮可動用的深層地下水釋水量。

(4)模型預(yù)測2018年平原區(qū)地下水水位恢復(fù)較快，后期呈現(xiàn)前快后慢，然后逐漸穩(wěn)定的趨勢。在現(xiàn)有的壓采方案下模型預(yù)測平原區(qū)地面沉降沒有繼續(xù)加大的趨勢，但也沒有出現(xiàn)全面回彈的現(xiàn)象。