賈艷輝，閻偉，高巍，馬春芽，李金山

基于地下水?dāng)?shù)值模擬的寶山農(nóng)場(chǎng)引江灌溉模式研究

賈艷輝1，閻偉2，高巍3，馬春芽1，李金山1*

（1.中國(guó)農(nóng)業(yè)科學(xué)院 農(nóng)田灌溉研究所，河南 新鄉(xiāng) 453002；2.河北省水利水電勘測(cè)設(shè)計(jì)研究院，天津 300250；3.河北省水利水電第二勘測(cè)設(shè)計(jì)研究院，石家莊 050021）

【】確定寶山農(nóng)場(chǎng)合理的農(nóng)業(yè)引江灌溉模式，為區(qū)域水資源優(yōu)化配置提供科學(xué)依據(jù)。針對(duì)松花江流域引江灌溉產(chǎn)生的地表水和地下水的綜合開(kāi)發(fā)利用問(wèn)題，利用三江平原水文地質(zhì)資料和寶山農(nóng)場(chǎng)地下水水位長(zhǎng)系列實(shí)測(cè)資料對(duì)MODFLOW-2000模型進(jìn)行了參數(shù)識(shí)別、驗(yàn)證及模型敏感性分析，并模擬分析了4種情景對(duì)寶山農(nóng)場(chǎng)地下水水位的影響。構(gòu)建的地下水模擬模型能夠反映灌區(qū)井群地下水位的實(shí)際變化過(guò)程。敏感性分析表明，長(zhǎng)序列地下水變化過(guò)程及地下水位分布情況受初始地下水位影響很?。徽f(shuō)明當(dāng)研究區(qū)域內(nèi)利用地下水灌溉及部分地區(qū)（沿江10 km）采用引江水灌溉模式時(shí)，地下水資源基本上可以維持供采平衡，在保證地下水合理開(kāi)采的前提下不會(huì)發(fā)生內(nèi)澇；而研究區(qū)域內(nèi)均采用引江灌溉模式則會(huì)使地下水位-抬高導(dǎo)致內(nèi)澇，內(nèi)澇情況隨著與河道距離的增加而愈發(fā)嚴(yán)重。引江水灌溉時(shí)，需加強(qiáng)對(duì)地下水位的觀測(cè)，當(dāng)?shù)叵滤贿^(guò)高時(shí)，應(yīng)減少引江灌溉用水，適當(dāng)增加地下水抽取量進(jìn)行灌溉，以保障研究區(qū)地下水和地表水的高效利用。

寶山；地下水模型；地表水；地下水；灌溉

0 引言

【研究意義】淺層地下水是可利用地下水資源的重要組成部分，也是陸地水資源的重要組成部分[1]，更是農(nóng)業(yè)灌溉用水及工業(yè)用水的重要來(lái)源[2]。不合理超采地下水或超量引水灌溉均會(huì)影響區(qū)域水資源的高效利用[3-4]。如何高效可持續(xù)開(kāi)發(fā)利用區(qū)域地表及地下水資源已成為水文、地質(zhì)及農(nóng)田水利等領(lǐng)域的研究熱點(diǎn)問(wèn)題。

合理地利用地表水和地下水資源，首要問(wèn)題是要對(duì)區(qū)域地下水進(jìn)行符合實(shí)際的模擬，進(jìn)而摸清區(qū)域地下水的運(yùn)動(dòng)規(guī)律?！狙芯窟M(jìn)展】目前針對(duì)地下水位動(dòng)態(tài)變化問(wèn)題的研究方法可分為2類(lèi)，一類(lèi)是基于長(zhǎng)期觀測(cè)資料的分析方法如Mann-Kendall檢驗(yàn)[5-6]、Yamamoto檢驗(yàn)[7]、雙向回歸分析方法[8-9]、R/S分析法[10-11]；另一類(lèi)是基于地下水運(yùn)動(dòng)模型的模擬方法如SWAT[12-14]、MODFLOW[15-17]等。前一類(lèi)方法一般需要大量的長(zhǎng)期觀測(cè)資料，而且一般只能預(yù)測(cè)某一點(diǎn)或幾點(diǎn)的地下水位而MODFLOW則可為地表水文過(guò)程與地下水系統(tǒng)之間的關(guān)系進(jìn)行耦合定量評(píng)價(jià)[18]、并可精確計(jì)算地下水不給量，為農(nóng)業(yè)灌溉地區(qū)的地下水資源合理規(guī)劃利用提供科學(xué)依據(jù)[19]。三江平原在“旱改水”過(guò)程中，水稻面積從最初的26.67萬(wàn)hm2發(fā)展到173.33萬(wàn)hm2，農(nóng)業(yè)灌溉對(duì)地下水的開(kāi)采量和引江灌溉需求與日俱增，致使很多地區(qū)出現(xiàn)了地下水位降落漏斗和內(nèi)澇?！厩腥朦c(diǎn)】因此，有必要利用較成熟的地下水運(yùn)動(dòng)模擬模型MODFLOW對(duì)寶山農(nóng)場(chǎng)地下水位動(dòng)態(tài)變化過(guò)程進(jìn)行建模，了解該區(qū)地下水的水位變動(dòng)情況【擬解決的關(guān)鍵問(wèn)題】利用MODFLOW軟件建模模擬分析不同情景（不同引江灌溉模式）條件下研究區(qū)域地下水位的動(dòng)態(tài)變化，以期促進(jìn)區(qū)域水資源的優(yōu)化配置和農(nóng)業(yè)可持續(xù)發(fā)展。

1 研究區(qū)概況

研究區(qū)位于三江平原上游的寶山農(nóng)場(chǎng)（130°47'E—130°57'E，46°55'N—47°01'N），研究區(qū)屬松花江流域，總面積2 491.7 km2，土壤為草甸土類(lèi)，黑土層30～50 cm，有機(jī)質(zhì)5.16%。研究區(qū)北部以松花江為邊界，東部以安邦河為邊界，南部以山腳線為邊界（圖1）。研究區(qū)屬大陸性季風(fēng)氣候。年日照時(shí)間2 300～2 700 d，無(wú)霜期130～150 d。年降水量500～650 mm。夏季受東南亞季風(fēng)影響，雨量充沛且集中，冬季干燥而少雨。

第四系含水層厚度一般100～150 m，靠近山前臺(tái)地區(qū)50～100 m。低平原大部分地區(qū)單井涌水量1 000～5 000 m3/d，局部地區(qū)大于5 000 m3/d。水文地質(zhì)條件分布情況如圖2所示。

圖1 研究區(qū)位置

圖2 三江平原水文地質(zhì)

2 研究區(qū)地下水動(dòng)態(tài)模型構(gòu)建

2.1 概念模型

三江平原屬于一個(gè)大型地下水匯水盆地，區(qū)域地下水流場(chǎng)主要受控于當(dāng)?shù)氐匦巍⑺暮蜌夂虻纫蛩豙26]。重力勢(shì)是驅(qū)動(dòng)盆地地下水運(yùn)動(dòng)的主要?jiǎng)恿Γ?dāng)降水或地表水轉(zhuǎn)變?yōu)榈叵滤畷r(shí)，便將相應(yīng)的重力勢(shì)能加諸地下水。在地形較高的山前地帶，地下水接受補(bǔ)給水位持續(xù)抬升，重力勢(shì)能積累，構(gòu)成勢(shì)源，為流場(chǎng)的高勢(shì)能區(qū)。地形低洼的河谷區(qū)，通常是流場(chǎng)的低勢(shì)區(qū)。地下水在重力勢(shì)能驅(qū)動(dòng)下，由山前地帶高勢(shì)能區(qū)向河谷的低勢(shì)能區(qū)運(yùn)動(dòng)，進(jìn)而形成三江平原上游區(qū)域地下水流動(dòng)系統(tǒng)。

研究區(qū)北部以松花江為界，東部以安邦河為界，南部以完達(dá)山山腳線為界。這樣北部和東部可以處理為第一類(lèi)邊界（水頭邊界）；南部山區(qū)的降水匯流后，一部分形成地表徑流，一部分形成地下徑流，徑流流量通過(guò)查閱當(dāng)?shù)氐乃牡刭|(zhì)資料并帶入模型進(jìn)行率定后確定，所以南部山區(qū)處理成第二類(lèi)邊界條件（流量邊界），水文地質(zhì)概念模型如圖3所示。根據(jù)地質(zhì)調(diào)查資料得知，含水層由南向北逐漸增厚。研究區(qū)的系統(tǒng)輸入可概化為4項(xiàng)：降水入滲、河流入滲、灌溉回補(bǔ)和上游地下水側(cè)向徑流補(bǔ)給；系統(tǒng)輸出也可以概化為4項(xiàng)：潛水蒸發(fā)、向河流排泄、人工開(kāi)采和向下游地下水側(cè)向徑流。

圖3 研究區(qū)水文地質(zhì)概念模型

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

MODFLOW模型，又稱(chēng)模塊化三維有限差分地下水流動(dòng)模型，廣泛應(yīng)用于水利學(xué)、農(nóng)業(yè)、生態(tài)環(huán)境以及城鄉(xiāng)發(fā)展規(guī)劃等領(lǐng)域。其數(shù)學(xué)模型為：

式中：K、K和K為滲透系數(shù)在、、方向上的分量；為水頭；為源匯項(xiàng)；s貯水率；為時(shí)間。

本文主要研究第四系潛水的動(dòng)態(tài)情況，建立研究區(qū)潛水運(yùn)動(dòng)方程：

式中：為滲透系數(shù)（m/d）；為地下水位（m）；為源匯項(xiàng)，地下水抽水量（d-1）；為給水度（m-1）；為時(shí)間（d）。

上邊界條件為灌溉、降水補(bǔ)給，計(jì)算式為：

式中：為縱坐標(biāo)（m）；(,)為地表；為補(bǔ)給系數(shù)；為降水及灌溉水量（m）。

研究區(qū)北部和東部邊界為第一類(lèi)邊界條件其數(shù)學(xué)表達(dá)式為：

式中：1、2分別為松花江及安邦河；1(,,)、2(,,)分別為松花江及安邦河水位函數(shù)。

式中：為垂直于邊界方向；3為南部山腳；(,,)為補(bǔ)給流量。

分別選用1999—2002年和2003—2007年較為完整的地下水位數(shù)據(jù)做模型的率定和驗(yàn)證，因此，使用1999年1月初地下水位作為初始條件：

式中：0(,)為初始地下水位分布。

2.3 模型離散

為了使用數(shù)值方法準(zhǔn)確求解研究區(qū)地下水位動(dòng)態(tài)變化過(guò)程，需對(duì)構(gòu)建的數(shù)學(xué)模型進(jìn)行空間和時(shí)間上的離散，首先對(duì)研究區(qū)進(jìn)行空間離散（圖4），把研究區(qū)剖分成1 km×1 km，黑色網(wǎng)格為活動(dòng)網(wǎng)格，綠色網(wǎng)格為非活動(dòng)網(wǎng)格。

圖4 研究區(qū)網(wǎng)格剖分

對(duì)研究區(qū)域64 a（1951—2014年）歷史氣象數(shù)據(jù)在時(shí)間上進(jìn)行離散，每年分為3個(gè)應(yīng)力期（5—8月為第1應(yīng)力期，應(yīng)力期內(nèi)有降水、灌溉提水；9—10月為第2應(yīng)力期，應(yīng)力期內(nèi)有降水；11月—次年4月為第3應(yīng)力期，應(yīng)力期內(nèi)有降水）共192個(gè)應(yīng)力期。

2.4 研究區(qū)水田分布

通過(guò)調(diào)查得知，2011年黑龍江省水稻種植分布[20]如圖5所示，對(duì)研究區(qū)的水田和旱田進(jìn)行概化處理后，約2/3的耕地面積是水田，其他為旱田，三江平原上游水田灌溉定額為716 mm，旱田灌溉定額為150 mm[21]。分布情況如圖6所示。

圖5 黑龍江省水稻種植分布

圖6 研究區(qū)水田旱田分布

2.5 參數(shù)率定及驗(yàn)證

把調(diào)查得到的水文地質(zhì)參數(shù)輸入模型中，采用1999—2002年寶山農(nóng)場(chǎng)地下水位監(jiān)測(cè)數(shù)據(jù)進(jìn)行模型參數(shù)率定，結(jié)果見(jiàn)圖7。

圖7 地下水位模擬值與實(shí)測(cè)值相關(guān)關(guān)系（1999—2003年）

將率定好的參數(shù)輸入模型，與2003—2007年的地下水位監(jiān)測(cè)數(shù)據(jù)進(jìn)行模型參數(shù)驗(yàn)證。模擬值與實(shí)測(cè)值分析結(jié)果（圖8）表明，模擬值與觀測(cè)值不同時(shí)段地下水位的變化趨勢(shì)基本一致，模擬值與觀測(cè)值擬合2=0.696 4，說(shuō)明構(gòu)建的地下水模擬模型能夠反映研究區(qū)地下水位的實(shí)際變化趨勢(shì)。

2.6 模型對(duì)初始地下水位的敏感性分析

地下水動(dòng)力學(xué)模型的定解條件包括邊界條件和初始條件，但是由于監(jiān)測(cè)點(diǎn)較少，初始條件中的地下水位分布信息與實(shí)際地下水位分布會(huì)有差別Δ0(,)，這樣的差別會(huì)對(duì)模擬的地下水位變化過(guò)程產(chǎn)生影響。因此，需分析0(,)對(duì)(,,)的影響。其他條件不變，只變化初始地下水位分布（全部初始水位+1 m和-1 m）后運(yùn)行程序，得到變化地下水位初始值后的水位分布（第64年末），模擬結(jié)果如圖9所示。

從圖9可以看出，初始地下水位+1 m和-1 m后，水位分布（第64年末）模擬結(jié)果與原始結(jié)果差別不大。另外，針對(duì)其中的1個(gè)觀測(cè)點(diǎn)（3號(hào)觀測(cè)點(diǎn)）模擬全過(guò)程的地下水位變化進(jìn)行對(duì)比，結(jié)果如圖10所示，同一觀測(cè)位置，變化初始地下水位后，前6年的地下水位有差別，1957—2015年地下水位變化基本一致。從觀測(cè)點(diǎn)地下水位及最終地下水位分布的比較可以看出，初始水位對(duì)長(zhǎng)期的地下水變化過(guò)程及地下水位分布情況影響很小。

圖8 地下水位模擬值與實(shí)測(cè)值相關(guān)關(guān)系（2003—2007年）

圖10 變化初始水位后觀測(cè)點(diǎn)地下水位變化對(duì)比

3 地下水位動(dòng)態(tài)模擬結(jié)果及討論

3.1 模擬結(jié)果

目前，研究區(qū)的水田及旱田均為井灌，引江水干支等骨干渠道都已建成，引松花江水灌溉正逐步推進(jìn)，所以，本文以三江平原地下水與地表水聯(lián)合利用為背景，設(shè)置3種聯(lián)合調(diào)度情景，與目前的全部使用地下水灌溉情況對(duì)比討論：情景1，旱田與水田都是井灌，與目前的水資源利用情況相同；情景2，沿松花江10 km的范圍是引江水灌溉，其余使用井灌；情景3，全部水田使用引江水灌溉，旱田使用井水灌溉；情景4，水田和旱田均為引江水灌溉。

分別對(duì)以上4種情景進(jìn)行了模擬，得到各時(shí)段地下水位和固定觀測(cè)點(diǎn)的地下水位過(guò)程。首先列出模擬期末（64年末）的地下水分布情況如圖11所示，從圖11中可以看出，從灌溉水利用情景1—情景4，隨著引江水量及引水灌溉范圍的增加，各級(jí)地下水位等值線都有向北移動(dòng)的趨勢(shì)；70 m地下水位等值線的形狀有向北凹變?yōu)橄虮蓖?，位置變化不明顯。75 m地下水位等值線的形狀和位置都有較大變化。100 m地下水位等值線的變化不大；地下水位對(duì)引江灌溉水量的響應(yīng)程度與地下水補(bǔ)給邊界和排泄邊界的距離有關(guān)，與補(bǔ)給排泄邊界距離越近響應(yīng)程度就越低，距離越遠(yuǎn)響應(yīng)程度越高。

圖11 研究區(qū)4種灌溉情景地下水分布情況

對(duì)上面4種情景的地下水位分布在圖中位置做剖面，可得到4條地下水位剖面線，如圖12所示。

圖12 研究區(qū)4種情景模擬期末地下水剖面線

從圖12中4條剖面線和地表線的比較可以看出，同一位置，隨著引江水量及引水灌溉范圍的增加，地下水位會(huì)升高，河邊的升高值較?。磺榫?、情景2、情景3，基本上實(shí)現(xiàn)了地下水資源的可持續(xù)開(kāi)采利用，不出現(xiàn)地下水超采（地下水位持續(xù)下降）或內(nèi)澇（地下水位上升至作物根區(qū)范圍）；情景4會(huì)使地下水排泄不及，從而導(dǎo)致內(nèi)澇。離河道越遠(yuǎn)，內(nèi)澇越嚴(yán)重。

對(duì)研究區(qū)內(nèi)4個(gè)觀測(cè)點(diǎn)（圖6）的地下水位變化過(guò)程進(jìn)行比較，如圖13所示，從選擇的4個(gè)觀測(cè)點(diǎn)64 a模擬地下水位變化過(guò)程可以看出：隨著與松花江距離的增加，年內(nèi)及年際地下水位變化幅度都會(huì)增加。并且比較同一觀測(cè)點(diǎn)的4種情景還可看出，隨著引江水量及引水灌溉范圍的增加，各觀測(cè)點(diǎn)的地下水位都有所增加，情景4的地下水位增加較多，說(shuō)明全部使用引江水灌溉水田有可能破壞地下水資源的采補(bǔ)平衡（灌溉對(duì)地下水補(bǔ)給量過(guò)多，存在地下水水位持續(xù)升高引起農(nóng)田漬害風(fēng)險(xiǎn)）。

圖13 研究區(qū)4種情景地下水位變化過(guò)程

3.2 討論

隨著三江平原水田灌溉模式和范圍不斷擴(kuò)大，對(duì)水資源的需求量也越來(lái)越大，使得水資源供需矛盾不斷加劇，淺層地下水位有逐年下降的趨勢(shì)。而充分利用江水資源，補(bǔ)充本區(qū)的不足水量，開(kāi)發(fā)流域農(nóng)業(yè)資源，通過(guò)地表水、地下水的聯(lián)合運(yùn)用，可從根本上解決研究區(qū)域內(nèi)水資源供需矛盾。牛松濤等[20]研究表明區(qū)域地下水動(dòng)態(tài)與農(nóng)業(yè)發(fā)展密切相關(guān)，水位與農(nóng)業(yè)灌溉需水量變化趨勢(shì)相反，主要原因是引灌水量過(guò)大，且有灌無(wú)排，地下水位大幅度上升會(huì)導(dǎo)致一系列的農(nóng)業(yè)生態(tài)環(huán)境問(wèn)題，其中最為突出的是漬澇鹽堿災(zāi)害[22]。本研究表明，研究區(qū)內(nèi)在目前開(kāi)采地下水灌溉的模式下，小范圍的使用引江水灌溉，不會(huì)造成區(qū)域地下水水位持續(xù)下降或持續(xù)上升。主要是因?yàn)榻邓凸喔仁怯绊憣毶睫r(nóng)場(chǎng)地下水位動(dòng)態(tài)變化的主要因子，但是提水灌溉不會(huì)導(dǎo)致地下水位持續(xù)下降[3]。李偉業(yè)等[23]在三江平原1997—2005年地下水的補(bǔ)水量和排泄量的比較中，虧損量累計(jì)達(dá) 6.89×108 m3，說(shuō)明該區(qū)域的補(bǔ)水能力在現(xiàn)有的開(kāi)采條件下在不斷地下降，地下水資源有超采的趨勢(shì)；另有研究表明，在1997—2017年，三江平原區(qū)域內(nèi)水稻田種植面積由 14.67萬(wàn)hm2增加到65.66萬(wàn)hm2，灌溉用水對(duì)地下水的需求量不斷增加，地下水位呈持續(xù)下降趨勢(shì)，下降幅度為1.60～9.29 m[24]。劉偉坡等[25]認(rèn)為三江平原地下水過(guò)量開(kāi)采引起的地下水流場(chǎng)變化，激發(fā)了河水入滲補(bǔ)給地下水能力，增加了地下水側(cè)向徑流補(bǔ)給量，改善了地下水徑流條件。而本文采用的引水灌溉模式，是灌區(qū)地下水埋深呈穩(wěn)步增大趨勢(shì)的主要影響因子，而在農(nóng)業(yè)灌溉對(duì)地下水的開(kāi)采量需求與日俱增的情況下，引江灌溉也是研究區(qū)維持地下水采補(bǔ)平衡的有效措施。

4 結(jié)論

1）MODFLOW模型可以較好模擬研究區(qū)地下水位時(shí)空動(dòng)態(tài)變化過(guò)程。

2）當(dāng)研究區(qū)域采用現(xiàn)有引江灌溉模式時(shí)，可維持當(dāng)?shù)氐叵滤裳a(bǔ)平衡；小范圍（沿江10 km范圍內(nèi)）使用引江水灌溉，不會(huì)造成區(qū)域地下水水位持續(xù)下降或持續(xù)上升；而大范圍（全部水田范圍）采用引江水灌溉時(shí)，導(dǎo)致研究區(qū)域內(nèi)部分排水不暢地區(qū)地下水水位持續(xù)上升，誘發(fā)農(nóng)田澇漬災(zāi)害。

3）引江水灌溉時(shí)，需加強(qiáng)對(duì)地下水位的觀測(cè)，如地下水位過(guò)高應(yīng)及時(shí)減少引江灌溉用水，適當(dāng)增加地下水抽取量進(jìn)行灌溉，以保障研究區(qū)水資源的平衡利用。

[1] 陳夢(mèng)熊. 試論地下水資源的基本概念與評(píng)價(jià)原則[J]. 自然資源, 1983, 5(4): 11-18.

CHEN Mengxiong. On the basic concept and evaluation principle of groundwater resources. Resources Science, 1983, 5(4): 11-18.

[2] 李佩成. 論實(shí)行最嚴(yán)格地下水資源管理[J]. 地下水, 2012, 34(3): 1-5.

LI Peicheng. Discussion on carrying out the most strict groundwater resources management[J]. Ground Water, 2012, 34(3): 1-5.

[3] 賈艷輝, 武玉剛, 朱文江, 等. 寶山農(nóng)場(chǎng)地下水動(dòng)態(tài)分析[J]. 灌溉排水學(xué)報(bào), 2018, 37(9): 73-78.

JIA Yanhui, WU Yugang, ZHU Wenjiang, et al. Dynamic of groundwater table in Baoshanfarm[J]. Journal of Irrigation and Drainage, 2018, 37(9): 73-78.

[4] 蓋兆梅, 劉仁濤, 付強(qiáng), 等. 地下水脆弱性風(fēng)險(xiǎn)預(yù)警和評(píng)估的區(qū)間三角模糊多屬性決策模型[J]. 數(shù)學(xué)的實(shí)踐與認(rèn)識(shí), 2016, 46(13): 208-214.

GAI Zhaomei, LIU Rentao, FU Qiang, et al. The multiple attribute decision making model and its application in the groundwater vulnerability risk early warning and assessment[J]. Mathematics in Practice and Theory, 2016, 46(13): 208-214.

[5] 婁寧. 基于時(shí)空數(shù)據(jù)挖掘技術(shù)的華東區(qū)域暴雨時(shí)空特征[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2017, 28(12): 4 043-4 050.

LOU Ning. Spatial-temporal dynamics of rainstorm in East China based on spatial-temporal data mining technology[J]. Chinese Journal of Applied Ecology, 2017, 28(12): 4 043-4 050.

[6] NOURANI V, DANANDEH MEHR A, AZAD N. Trend analysis of hydroclimatological variables in Urmia lake basin using hybrid wavelet Mann-Kendall and ?entests[J]. Environmental Earth Sciences, 2018, 77(5): 1-18.

[7] 殷紅, 張美玲, 辛明月, 等. 近50年沈陽(yáng)氣溫變化與城市化發(fā)展的關(guān)系[J]. 生態(tài)環(huán)境學(xué)報(bào), 2011, 20(3): 544-548.

YIN Hong, ZHANG Meiling, XIN Mingyue, et al. The changes of air temperature in Shenyang in recent 50 years and its relation to urbanization[J]. Ecology and Environmental Sciences, 2011, 20(3): 544-548.

[8] 馬天儒, 夏達(dá)忠, 張行南, 等. 基于雙向交互式回歸模型的鎮(zhèn)江站高潮位預(yù)報(bào)[J]. 水電能源科學(xué), 2015, 33(2): 50-52.

MA Tianru, XIA Dazhong, ZHANG Xingnan, et al. Peak tide level forecasting for Zhenjiang hydrological station based on bidirectional interactive regression model[J]. Water Resources and Power, 2015, 33(2): 50-52.

[9] 陳斌, 包為民, 瞿思敏, 等. 雙向線性回歸法在椒江臨海站水位預(yù)報(bào)中的應(yīng)用[J]. 水文, 2008, 28(3): 45-48.

CHEN Bin, BAO Weimin, QU Simin, et al. Application of bidirectional propagation and multivariate linear regression method in stage forecasting at Linhai station on Jiaojiang river[J]. Journal of China Hydrology, 2008, 28(3): 45-48.

[10] 李建林, 昝明軍, 李寶玲. 基于R/S分析的黑河出山年徑流量灰色預(yù)測(cè)[J]. 地域研究與開(kāi)發(fā), 2014, 33(5): 127-131.

LI Jianlin, ZAN Mingjun, LI Baoling. Grey prediction of out-mountainous annual runoff of Heihe river based on R/S analysis[J]. Areal Research and Development, 2014, 33(5): 127-131.

[11] ZHAO G S, ZHOU G Q, ZHAO X D, et al. R/S analysis for stress evolution in shaft lining and fracture prediction method[J]. Advanced Materials Research, 2011, 374-377.

[12] SUN H, CORNISH P S. Estimating shallow groundwater recharge in the headwaters of the Liverpool Plains using SWAT[J]. Hydrological Processes, 2005, 19(3): 795-807.

[13] NASIR M J, KHAN S, ZAHID H, et al. Delineation of groundwater potential zones using GIS and multi influence factor (MIF) techniques: A study of district Swat, Khyber Pakhtunkhwa, Pakistan[J]. Environmental Earth Sciences, 2018, 77(10): 1-11.

[14] AHN S R, KIM S J. Analysis of water balance by surface-groundwater interaction using the SWAT model for the Han River basin, South Korea[J]. Paddy and Water Environment, 2018, 16(3): 543-560.

[15] BRUNNER P, SIMMONS C T, COOK P G, et al. Modeling surface water-groundwater interaction with MODFLOW: Some considerations [J]. Groundwater, 2010, 48(2): 174-180.

[16] BUSHIRA K M, HERNANDEZ J R, SHENG Z P. Surface and groundwater flow modeling for calibrating steady state using MODFLOW in Colorado River Delta, Baja California, Mexico[J]. Modeling Earth Systems and Environment, 2017, 3(2): 815-824.

[17] SULE B F, AYENIGBA S E. Application of GMS - MODFLOW to investigate groundwater development potential in river meme catchment, kogi state, Nigeria[J]. International Journal of Sciences, 2017, 3(9): 39-51.

[18] TESTONI R, LEVIZZARI R, SALVE M D. Coupling of unsaturated zone and saturated zone in radionuclide transport simulations[J]. Progress in Nuclear Energy, 2017, 95: 84-95.

[19] TWARAKAVI N K C, SIMUNEK J, SEO S. Evaluating interactions between groundwater and vadose zone using the HYDRUS-based flow package for MODFLOW[J]. Vadose Zone Journal, 2008, 7(7): 757-768.

[20] 牛松濤, 鄭迎春, 韓淑敏. 農(nóng)業(yè)發(fā)展對(duì)區(qū)域地下水資源的影響[J]. 南水北調(diào)與水利科技, 2012, 10(6): 104-109.

NIU Songtao, ZHENG Yingchun, HAN Shumin. Effects of agricultural development on regional groundwater resources[J]. South-to-North Water Diversion and Water Science & Technology, 2012, 10(6): 104-109.

[21] 劉克寶, 陸忠軍, 劉述彬, 等. 基于RS的黑龍江省水稻種植空間分布格局研究[J]. 黑龍江農(nóng)業(yè)科學(xué), 2015(8): 136-141.

LIU Kebao, LU Zhongjun, LIU Shubin, et al. Study on spatial distribution of rice planting in Heilongjiang Province based on RS[J]. Heilongjiang Agricultural Sciences, 2015(8): 136-141.

[22] 中國(guó)科學(xué)院黃土高原綜合科學(xué)考察隊(duì). 黃土高原地區(qū)水資源問(wèn)題及其對(duì)策[M]. 北京: 中國(guó)科學(xué)技術(shù)出版社, 1991.

[23] 李偉業(yè), 付強(qiáng). 基于GIS的三江平原地下水資源評(píng)價(jià)[J]. 水土保持研究, 2007, 14(4): 92-95.

LI Weiye, FU Qiang. Evaluation of groundwater resources using GIS in the Sanjiangplain[J]. Research of Soil and Water Conservation, 2007, 14(4): 92-95.

[24] 席文娟．寧夏石嘴山地區(qū)地下水水質(zhì)模擬與變化趨勢(shì)分析[D]．西安: 長(zhǎng)安大學(xué)，2013.

XI Wenjuan. Simulation of Groundwater Quality and Ariation Tendency in Shuizuishan area, Ningxia, Northwest China[D]. Xi’an: Chang’an University, 2013.

[25] 劉偉坡, 沙娜, 程旭學(xué). 三江平原建三江地下水動(dòng)態(tài)變化特征[J]. 灌溉排水學(xué)報(bào), 2020, 39(5): 96-101.

LIU Weipo, SHA Na, CHENG Xuxue. Groundwater dynamics in catchment of Jiansanjiang in Sanjiangplain[J]. Journal of Irrigation and Drainage, 2020, 39(5): 96-101.

Using Numerical Simulation to Aid Surface and Subsurface Water Irrigation in Baoshan of Heilongjiang Province

JIA Yanhui1, YAN Wei2, GAO Wei3, MA Chunya1, LI Jinshan1*

(1.Farmland Irrigation Research Institute, Chinese Academy of Agricultural Sciences, Xinxiang 453002, China;2.Investigation Design and Research Institute of Water Conservancy and Hydropower of Hebei Province, Tianjin 300250, China; 3.The Second Design and Research Institute of Water Conservancy and Hydropower of Hebei Province, Shijiazhuang 050021, China)

【】The project of transforming drylands to paddy fields in Sanjiang Plain in the northeast China has increased its rice cultivation areas from 266 700 hm2prior to the project to 1 733 300 hm2after implementing the project. As a result, its demand for water also increased. Baoshan Farm is in the upper reach of Sanjiang plain on a piedmont, and its agricultural production had relied on groundwater since its establishment. Since 1992, it has gradually converted drylands to paddy field, and such a land conversion has resulted in a conflict between increased demand for water and limited water availability.【】The purpose of this paper is to analyze how to optimize the use of surface water and groundwater to resolve the conflict between limited water supply and the increased demand for water due to the land use change, so as to avoid waterlogging while in the meantime alleviating continued groundwater table dropping.【】The calculation was based on groundwater flow with the groundwater dynamics under different managements simulated using the software Modflow-2000. The aquifer parameters, model validation and sensitivity analysis of the parameters were carried out using long-term hydrogeological data measured from the area. We then compared groundwater dynamics under irrigations using different combinations of surface water and groundwater.【】The calibrated groundwater dynamics model correctly reproduced the measured changes in groundwater level in the studied area; sensitivity analysis revealed that the initial groundwater level did not have a noticeable influence on long-term groundwater variation. Using the river water to irrigate areas 10 km in the proximity of the river while irrigating other areas using the groundwater can keep groundwater resource balance while in the meantime alleviating waterlogging. In contrast, irrigating too many areas using transferred river water could raise the groundwater table and result in waterlogging in a wide range of areas.【】Rational use of groundwater and surface water in irrigation can maintain groundwater balance and avoid waterlogging in the farmlands of Baoshai; a monitoring system should be in place to timely measure the change in groundwater table.

Baoshai; groundwater model; surface water; groundwater; irrigation

賈艷輝, 閻偉, 高巍, 等. 基于地下水?dāng)?shù)值模擬的寶山農(nóng)場(chǎng)引江灌溉模式研究[J]. 灌溉排水學(xué)報(bào), 2021, 40(9): 111-117.

JIA Yanhui, YAN Wei, GAO Wei, et al. Using Numerical Simulation to Aid Surface and Subsurface Water Irrigation in Baoshan of Heilongjiang Province[J]. Journal of Irrigation and Drainage, 2021, 40(9): 111-117.

S24

A

10.13522/j.cnki.ggps.2020528

1672 – 3317（2021）09 - 0111 - 07

2020-09-20

中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)田灌溉研究所基本科研業(yè)務(wù)費(fèi)項(xiàng)目（FIRI20210203）；中國(guó)農(nóng)業(yè)科學(xué)院創(chuàng)新工程黃淮海冬小麥-夏玉米集約高效用水技術(shù)創(chuàng)新與應(yīng)用

賈艷輝（1982-），男。副研究員，博士，主要從事灌區(qū)水資源研究。E-mail: jyh_5151@126.com

李金山（1972-），男。研究員，主要從事節(jié)水灌溉研究。E-mail: lijinshan72@126.com

責(zé)任編輯：趙宇龍