岳瑜素，王宏偉，江恩慧，彭緒庶，時(shí)芳欣

（1.黃河水利科學(xué)研究院，河南 鄭州 450003；2.水利部堤防安全與病害防治工程技術(shù)研究中心，河南 鄭州 450003；3.中國(guó)社會(huì)科學(xué)院 數(shù)量經(jīng)濟(jì)與技術(shù)經(jīng)濟(jì)研究所，北京 100732）

1 研究背景

黃河下游河道內(nèi)分布有廣闊的灘地，占了下游河道總面積的三分之二，涉及到豫魯兩省15個(gè)市45個(gè)縣（區(qū)），其中居住人口近190萬(wàn)，現(xiàn)有耕地22萬(wàn)hm2。黃河下游灘區(qū)是黃河防洪工程體系的重要組成部分，不僅要發(fā)揮行洪、滯洪、沉沙功能，同時(shí)也是廣大灘區(qū)民眾賴以生存和生活的場(chǎng)所。因此，灘區(qū)的防洪運(yùn)用和社會(huì)經(jīng)濟(jì)發(fā)展矛盾長(zhǎng)期存在，民眾生產(chǎn)生活水平低，是受流域和地方政府共同關(guān)注的特殊區(qū)域。習(xí)近平總書(shū)記在2019年9月18日召開(kāi)的黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展座談會(huì)上明確提出，要實(shí)施黃河下游河道和灘區(qū)綜合提升治理工程。

面對(duì)新的挑戰(zhàn)，國(guó)內(nèi)眾多學(xué)者做了諸多有益的探索，圍繞黃河下游生產(chǎn)堤留存與廢除、灘區(qū)治理模式及補(bǔ)償政策、下游河道灘槽協(xié)同治理等取得了一批重要成果，并在水利部黃河水利委員會(huì)的組織下促成了“黃河下游治理方略專家研討會(huì)”[1]。胡春宏[2]院士在對(duì)黃河水沙變化及發(fā)展趨勢(shì)分析的基礎(chǔ)上指出，面對(duì)未來(lái)水沙持續(xù)減少的新形勢(shì)，可考慮適度利用控導(dǎo)工程、生產(chǎn)堤等建設(shè)防洪堤，束窄現(xiàn)狀河道寬度，使大部分灘區(qū)成為永久安全區(qū)。何予川等[3]指出，束窄河道的治理方略可以更好地解決灘區(qū)安全與發(fā)展的問(wèn)題，但考慮到尚不完善的黃河下游防洪減淤體系，其可能會(huì)帶來(lái)較大的防洪風(fēng)險(xiǎn)和壓力。江恩慧等[4]基于黃河下游中常高含沙洪水實(shí)體模型試驗(yàn)系列成果，結(jié)合黃河下游不同河段的河勢(shì)特征，提出了相對(duì)完善的治理運(yùn)行建議，即在高村以上灘區(qū)修建高標(biāo)準(zhǔn)生產(chǎn)堤以保障灘區(qū)安全和經(jīng)濟(jì)建設(shè)，高村以下灘區(qū)實(shí)施灘區(qū)淹沒(méi)補(bǔ)償政策。此外，牛玉國(guó)等[5]提出了灘區(qū)分區(qū)治理的思路與實(shí)施方式，張金良[6]也提出了灘區(qū)再造與生態(tài)治理的方案等等。然而，目前研究大多站位于防洪保安全，從水利工程建設(shè)和管理角度開(kāi)展黃河下游灘區(qū)治理研究，未將防洪與自然、經(jīng)濟(jì)發(fā)展統(tǒng)籌考慮，得到的結(jié)果也具有一定的局限性。

黃河下游灘區(qū)治理是一個(gè)復(fù)雜的系統(tǒng)工程問(wèn)題，涉及洪水治理、土地資源承載和區(qū)域經(jīng)濟(jì)發(fā)展需求等多個(gè)方面，存在多重直接或間接的聯(lián)動(dòng)關(guān)系[7]。本文在深入分析黃河下游寬灘區(qū)滯洪沉沙作用和社會(huì)經(jīng)濟(jì)發(fā)展需求的基礎(chǔ)上，綜合考慮防洪安全、產(chǎn)業(yè)發(fā)展、農(nóng)村建設(shè)、資源配置等影響因素及其聯(lián)動(dòng)關(guān)系，以自然-經(jīng)濟(jì)-社會(huì)協(xié)調(diào)發(fā)展為目標(biāo)，構(gòu)建適應(yīng)新時(shí)代的灘區(qū)系統(tǒng)動(dòng)力學(xué)模型（SD模型），并通過(guò)仿真模擬研究探索灘區(qū)可持續(xù)發(fā)展的優(yōu)化治理模式，以期為黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展提供科學(xué)依據(jù)。

2 研究方法

2.1 SD模型復(fù)合系統(tǒng)理論認(rèn)為，復(fù)合系統(tǒng)包含多個(gè)子系統(tǒng)，這些子系統(tǒng)之間既相互獨(dú)立，又相互作用。復(fù)合系統(tǒng)不是由各子系統(tǒng)簡(jiǎn)單迭加而成，而是子系統(tǒng)有機(jī)復(fù)合后的總體優(yōu)化與協(xié)調(diào)結(jié)果[8]。黃河下游灘區(qū)自然-經(jīng)濟(jì)-社會(huì)復(fù)合系統(tǒng)屬于典型的農(nóng)業(yè)復(fù)合系統(tǒng)[9-10]，其中自然子系統(tǒng)主要包括灘區(qū)土地、水、光照等自然資源及其之間的相互作用關(guān)系；經(jīng)濟(jì)子系統(tǒng)主要包括灘區(qū)百姓為生存和發(fā)展而發(fā)生的農(nóng)業(yè)及配套服務(wù)等活動(dòng)；社會(huì)子系統(tǒng)主要包括人口，體制、文化等。三個(gè)子系統(tǒng)在時(shí)間、空間、數(shù)量、結(jié)構(gòu)、秩序方面生態(tài)耦合和相互作用，共同決定了復(fù)合生態(tài)系統(tǒng)的發(fā)展與演替方向[11]。

本文考慮從行洪輸沙和經(jīng)濟(jì)社會(huì)協(xié)調(diào)發(fā)展的角度對(duì)灘區(qū)自然－經(jīng)濟(jì)－社會(huì)復(fù)合系統(tǒng)進(jìn)行分析，構(gòu)建適用于黃河下游灘區(qū)的系統(tǒng)動(dòng)力學(xué)模型[12-13]，并利用因果關(guān)系圖、流圖和系統(tǒng)動(dòng)力學(xué)方程刻畫(huà)洪水治理、土地資源配置和灘區(qū)經(jīng)濟(jì)社會(huì)發(fā)展之間的聯(lián)動(dòng)關(guān)系，結(jié)果見(jiàn)圖1。可以看出，灘區(qū)行洪輸沙-經(jīng)濟(jì)-社會(huì)協(xié)調(diào)發(fā)展因果關(guān)系網(wǎng)絡(luò)主要包括三個(gè)回路，分別為灘區(qū)耕地面積反饋回路、灘區(qū)農(nóng)業(yè)生產(chǎn)率反饋回路和灘區(qū)人口變動(dòng)反饋回路。

圖1 黃河下游灘區(qū)行洪輸沙-經(jīng)濟(jì)-社會(huì)協(xié)調(diào)發(fā)展系統(tǒng)動(dòng)力學(xué)流圖

其中，模型將黃河來(lái)水來(lái)沙和治理模式選擇等外生影響，抽象簡(jiǎn)化為對(duì)應(yīng)仿真情形下淹沒(méi)的灘區(qū)面積，農(nóng)村居民當(dāng)季的播種面積相應(yīng)減少，并進(jìn)一步造成一定程度的農(nóng)作物損失和農(nóng)村居民經(jīng)營(yíng)性收入減少。黃河來(lái)水來(lái)沙和不同治理模式選擇等外生沖擊通過(guò)灘區(qū)的農(nóng)業(yè)生產(chǎn)活動(dòng)，逐步傳導(dǎo)向整個(gè)灘區(qū)經(jīng)濟(jì)社會(huì)系統(tǒng)，引起灘區(qū)土地資源配置、農(nóng)村居民經(jīng)營(yíng)性收入和農(nóng)村人口等因素的互動(dòng)影響和循環(huán)反饋。

2.2 模型構(gòu)建在構(gòu)建黃河下游灘區(qū)自然－經(jīng)濟(jì)－社會(huì)復(fù)合系統(tǒng)動(dòng)力學(xué)模型時(shí)，本研究基于系統(tǒng)動(dòng)力學(xué)流圖中的系統(tǒng)運(yùn)動(dòng)過(guò)程因果關(guān)系，構(gòu)建不同的系統(tǒng)因子計(jì)算公式。

（1）灘區(qū)耕地面積。耕地面積增加和減少受到收入差距、灘區(qū)總面積、建設(shè)面積等因素的影響。具體函數(shù)關(guān)系如下所示：

式中：GDAR為耕地面積狀態(tài)變量；GDAIC為耕地面積增加面積；GDADC為耕地面積減少面積；TQAR為灘區(qū)總面積；JZAR為村莊建設(shè)用地面積。灘區(qū)總面積與村莊建設(shè)用地的差值作為耕地面積模擬的初始值，耕地面積增加與減少的差值作為各年凈增加的耕地面積。式（2）中，SFL為設(shè)定的耕地安全線[14]，采用灘區(qū)總面積與農(nóng)村建設(shè)面積底線的差額，農(nóng)村建設(shè)面積底線為農(nóng)村最小人均住房面積與常住人口的乘積；JZARIC為各年建設(shè)用地復(fù)耕面積；INCMD為農(nóng)村居民收入差距；α1和α2分別為耕地調(diào)整系數(shù)和調(diào)整常量，二者根據(jù)各年各縣農(nóng)村居民收入差距和耕地面積增加量擬合回歸得到。式（3）為耕地面積增加函數(shù)，其中PLE為遷出居民人數(shù)；PL為常住人口數(shù)量。以上涉及到的面積單位均為km2。

（2）灘區(qū)農(nóng)業(yè)生產(chǎn)。在特定標(biāo)準(zhǔn)防護(hù)堤或其他灘區(qū)治理方式下，洪水可能造成部分灘區(qū)淹沒(méi)，減少當(dāng)季農(nóng)村居民的農(nóng)作物播種面積、造成部分農(nóng)作物損失。具體函數(shù)關(guān)系如下式所示：

式中：BZAR為播種面積；mltc為復(fù)種指數(shù)；CYPD為農(nóng)作物產(chǎn)量；CPYPA為單位面積農(nóng)作物產(chǎn)量；CPLS為洪水淹沒(méi)農(nóng)作物損失。式（6）（7）分別為各年洪水淹沒(méi)面積和淹沒(méi)農(nóng)作物損失的估算方程，由1958年和1977年來(lái)水來(lái)量、淹沒(méi)面積和農(nóng)作物損失量回歸擬合，得到參數(shù)β1、β2、β3和β5。以上涉及到的面積單位均為km2。

根據(jù)灘區(qū)農(nóng)業(yè)生產(chǎn)率反饋回路反映的因果互動(dòng)關(guān)系，農(nóng)業(yè)機(jī)械動(dòng)力和農(nóng)業(yè)勞動(dòng)力兩類(lèi)生產(chǎn)要素投入增加是影響單位面積農(nóng)作物產(chǎn)量，即農(nóng)業(yè)生產(chǎn)效率的重要因素，且兩類(lèi)生產(chǎn)要素之間能夠互相替代。具體的函數(shù)關(guān)系如下所示：

式（8）為農(nóng)業(yè)生產(chǎn)率函數(shù)，其中CPYPA為單位面積農(nóng)作物產(chǎn)量；γ0+γ1t代表技術(shù)進(jìn)步，γ1為技術(shù)進(jìn)步系數(shù)；LIPA為單位面積勞動(dòng)力投入；δ1為勞動(dòng)投入彈性系數(shù)；MCPA為單位面積機(jī)械動(dòng)力投入；δm為機(jī)械投入彈性系數(shù)。式（9）為單位面積勞動(dòng)投入函數(shù)，其中EMPLR為農(nóng)村居民就業(yè)率。式（10）為單位面積機(jī)械動(dòng)力投入函數(shù)，其中MC為農(nóng)業(yè)機(jī)械動(dòng)力。式（11）為農(nóng)業(yè)機(jī)械動(dòng)力函數(shù)，其中MC0為期初農(nóng)業(yè)機(jī)械存量；NMC為當(dāng)年新增農(nóng)業(yè)機(jī)械；MCDC為農(nóng)業(yè)機(jī)械折舊。式（12）為新增農(nóng)業(yè)機(jī)械的決定函數(shù)，本文假定新增農(nóng)業(yè)機(jī)械與農(nóng)村居民收入（INCM）水平線性相關(guān)，參數(shù)θ1、θ2根據(jù)歷年農(nóng)村居民收入水平和各年新增農(nóng)業(yè)機(jī)械回歸擬合得到。

灘區(qū)農(nóng)村居民收入中，較大部分來(lái)自農(nóng)業(yè)經(jīng)營(yíng)性收入。洪水淹沒(méi)灘區(qū)耕地面積，導(dǎo)致當(dāng)季播種面積減少和農(nóng)作物損失，農(nóng)業(yè)收入可能相應(yīng)下降，農(nóng)村居民收入水平可能低于周?chē)貐^(qū)平均水平。具體函數(shù)關(guān)系如下所示：

式（13）為農(nóng)業(yè)收入函數(shù)，其中AGINCM為農(nóng)業(yè)經(jīng)營(yíng)性收入；CRPUP為農(nóng)作物單價(jià)；CPYD為農(nóng)作物產(chǎn)量。式（14）為農(nóng)村居民人均收入函數(shù)，其中INCMPC為人均收入水平；為人均農(nóng)業(yè)收入水平；PAGINCM為農(nóng)業(yè)經(jīng)營(yíng)收入占總收入的比重。式（15）為收入差距函數(shù)，其中INCMGAP為收入差距水平；INCMOT為周?chē)貐^(qū)農(nóng)村居民人均收入水平；INCMPC為當(dāng)?shù)剞r(nóng)村居民人均收入水平。

（3）灘區(qū)農(nóng)業(yè)人口。農(nóng)村人口變動(dòng)子系統(tǒng)的結(jié)構(gòu)相對(duì)簡(jiǎn)單，僅反映農(nóng)村常住人口變動(dòng)的主要影響因素。具體函數(shù)關(guān)系如下所示：

式（16）為農(nóng)村常住人口函數(shù)，其中PL為農(nóng)村常住人口，是期初人口數(shù)與各期人口增減的累加；PL0為期初農(nóng)村人口數(shù)量；PLI為人口增加量；PLD為人口減少量。式（17）為人口增長(zhǎng)函數(shù)，主要影響因素為人口出生，其中bthr為人口出生率。式（18）為人口減少函數(shù)，主要影響因素為人口死亡和人口遷出，其中dthr為人口死亡率；PLE為人口遷出量。式（19）為人口遷出函數(shù)，其中參數(shù)ρ1、ρ2分別為遷出系數(shù)和常數(shù)，可根據(jù)各年凈遷出人口和收入差距回歸擬合得到。

2.3 模擬情景設(shè)置本文在采用系統(tǒng)動(dòng)力學(xué)模型進(jìn)行模擬前，先按灘區(qū)不同運(yùn)用方案為設(shè)置情景。主要包括兩類(lèi)：①無(wú)防護(hù)堤方案，全面破除生產(chǎn)提，洪水自然漫灘；②防護(hù)堤方案，基本以現(xiàn)狀防護(hù)堤為基礎(chǔ)，通過(guò)調(diào)整、改造、加固，建設(shè)成為保護(hù)堤，以抵御一定標(biāo)準(zhǔn)的洪水，在防護(hù)堤內(nèi)生產(chǎn)生活，防護(hù)堤防洪標(biāo)準(zhǔn)分別為6000、8000和10 000 m3/s，均采用有閘門(mén)的防護(hù)堤方案[15-16]。河道地形以2013年河道邊界地形為基礎(chǔ)條件。不同的灘區(qū)治理方案和水沙過(guò)程，產(chǎn)生的淹沒(méi)面積也不同，對(duì)灘區(qū)經(jīng)濟(jì)社會(huì)系統(tǒng)運(yùn)行的影響也會(huì)有差異。本文采用這兩種模擬情景，模擬2016—2030年黃河下游來(lái)水來(lái)沙對(duì)灘區(qū)經(jīng)濟(jì)社會(huì)系統(tǒng)發(fā)展產(chǎn)生的影響。

3 模型檢驗(yàn)

為評(píng)估構(gòu)建的系統(tǒng)動(dòng)力學(xué)模型對(duì)現(xiàn)狀的模擬能力，采用歷史統(tǒng)計(jì)數(shù)據(jù)與仿真結(jié)果對(duì)模型進(jìn)行誤差檢驗(yàn)[16]。選取黃河下游灘區(qū)2006—2015年的耕地面積、播種面積、人均收入和常住人口等社會(huì)經(jīng)濟(jì)因子作為模型仿真的歷史數(shù)據(jù)，仿真步長(zhǎng)為1年。由于黃河下游灘區(qū)并非獨(dú)立的行政建制區(qū)，缺乏直接的社會(huì)經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)，本文選擇灘區(qū)所在的行政建制鄉(xiāng)鎮(zhèn)和縣（市）資料進(jìn)行研究，共調(diào)研收集了河南14個(gè)縣（市）97個(gè)鄉(xiāng)鎮(zhèn)以及山東菏澤東明縣、鄄城、鄆城三個(gè)縣15個(gè)鄉(xiāng)鎮(zhèn)的農(nóng)業(yè)數(shù)據(jù)，模型歷史檢驗(yàn)結(jié)果見(jiàn)表1?？梢钥闯觯瑑H有部分變量在部分年份的仿真值出現(xiàn)了超過(guò)10%的誤差，誤差主要集中在人均收入方面。造成人均收入預(yù)測(cè)誤差的原因主要是由于現(xiàn)在農(nóng)村發(fā)展需求導(dǎo)致大部分農(nóng)村勞動(dòng)力進(jìn)城務(wù)工，現(xiàn)階段灘區(qū)的經(jīng)濟(jì)效益主要來(lái)源于耕地，勞動(dòng)力的流失導(dǎo)致人均收入遠(yuǎn)低于預(yù)期值。但整體而言，模型預(yù)估灘區(qū)社會(huì)經(jīng)濟(jì)系統(tǒng)因子的變化趨勢(shì)基本與統(tǒng)計(jì)值趨勢(shì)相符，該系統(tǒng)動(dòng)力模型基本可以代表黃河下游灘區(qū)行洪輸沙-經(jīng)濟(jì)-社會(huì)協(xié)調(diào)發(fā)展的現(xiàn)狀，可以用于預(yù)測(cè)未來(lái)灘區(qū)治理情形下的發(fā)展趨勢(shì)。

4 模擬仿真結(jié)果及分析

采用構(gòu)建的系統(tǒng)動(dòng)力學(xué)模型，基于不同治理方案對(duì)灘區(qū)經(jīng)濟(jì)社會(huì)發(fā)展情景進(jìn)行仿真模擬，從耕地面積、播種面積、單位面積農(nóng)作物產(chǎn)量、農(nóng)業(yè)收入水平、農(nóng)村居民人均收入水平、常住人口等方面進(jìn)行比較，分析黃河下游來(lái)水來(lái)沙對(duì)灘區(qū)經(jīng)濟(jì)社會(huì)系統(tǒng)發(fā)展產(chǎn)生的影響，結(jié)果見(jiàn)圖2—圖7。其中，現(xiàn)狀指現(xiàn)有無(wú)防護(hù)堤的狀態(tài)；6000、8000和10 000 m3/s指具有該流量防護(hù)堤的狀態(tài)。

表2 模型歷史檢驗(yàn)結(jié)果統(tǒng)計(jì)

圖2 不同治理方案下灘區(qū)耕地面積情景模擬

圖3 不同治理方案下灘區(qū)播種面積情景模擬

圖4 不同治理方案下灘區(qū)單位面積農(nóng)作物產(chǎn)量情景模擬

圖5 不同治理方案下灘區(qū)農(nóng)業(yè)收入情景模擬

圖6 不同治理方案下灘區(qū)農(nóng)村居民人均收入情景模擬

圖7 不同治理方案下灘區(qū)農(nóng)村常住人口情景模擬

結(jié)合圖1中灘區(qū)耕地面積反饋回路可知，不同防護(hù)堤方案下，雖然保障條件提高有利于增加耕地面積，但耕地面積變化同時(shí)還受灘區(qū)總面積和建設(shè)面積的約束，后者則受灘區(qū)人口規(guī)模變動(dòng)影響。相應(yīng)地，堤防工程的防洪標(biāo)準(zhǔn)越高，灘區(qū)保護(hù)的農(nóng)村居民人口越多，居民農(nóng)業(yè)生產(chǎn)用地與生活居住用地此消彼長(zhǎng)的關(guān)系也更明顯，灘區(qū)耕地面積呈現(xiàn)出一定的周期性波動(dòng)特征。相應(yīng)地，無(wú)堤方案和6000 m3/s標(biāo)準(zhǔn)防護(hù)堤方案下，灘區(qū)農(nóng)村居民人口在長(zhǎng)期將減少或保持穩(wěn)定，灘區(qū)耕地面積變動(dòng)幅度相對(duì)較小。8000和10 000 m3/s標(biāo)準(zhǔn)防護(hù)堤方案下，灘區(qū)農(nóng)村人口在長(zhǎng)期將呈增長(zhǎng)趨勢(shì)，耕地面積周期性波動(dòng)幅度更大，其中8000 m3/s 標(biāo)準(zhǔn)防護(hù)堤方案下的人口增幅較小，耕地面積波動(dòng)的周期也較短；10 000 m3/s標(biāo)準(zhǔn)防護(hù)堤方案下的人口增幅更大，耕地面積波動(dòng)的周期相對(duì)較長(zhǎng)。同理，播種面積變化也具有類(lèi)似特征和趨勢(shì)。不同的是，播種面積還受復(fù)種指數(shù)影響。在無(wú)堤方案和6000 m3/s標(biāo)準(zhǔn)防護(hù)堤方案中，灘區(qū)農(nóng)作物復(fù)種指數(shù)分別為1.1和1.2，而在8000和10 000 m3/s 標(biāo)準(zhǔn)防護(hù)堤方案中，由于對(duì)農(nóng)業(yè)生產(chǎn)保護(hù)程度更高，復(fù)種指數(shù)達(dá)到1.4。因此，農(nóng)作物播種面積雖然仍隨耕地面積變化呈現(xiàn)周期性波動(dòng)，但波動(dòng)幅度明顯減小。

從經(jīng)濟(jì)社會(huì)發(fā)展情況來(lái)看：（1）無(wú)堤方案和6000 m3/s防護(hù)堤方案對(duì)灘區(qū)農(nóng)業(yè)生產(chǎn)和居民生活的保障較弱，洪水淹沒(méi)面積較大，造成農(nóng)作物損失也較大，灘區(qū)農(nóng)村居民人口長(zhǎng)遠(yuǎn)將減少，農(nóng)業(yè)收入水平也相應(yīng)較低；（2）8000和10 000 m3/s標(biāo)準(zhǔn)防護(hù)堤方案能夠較好保障灘區(qū)居民的生產(chǎn)生活，農(nóng)作物復(fù)種指數(shù)基本相當(dāng)于周邊地區(qū)平均水平，農(nóng)業(yè)生產(chǎn)受洪水影響較小，農(nóng)業(yè)收入水平相對(duì)較高；（3）8000和10 000 m3/s 標(biāo)準(zhǔn)防護(hù)堤方案均為高標(biāo)準(zhǔn)防護(hù)堤方案，但二者的長(zhǎng)期人口變動(dòng)趨勢(shì)存在差異。其中，10 000 m3/s標(biāo)準(zhǔn)防護(hù)堤方案的人口增幅較大，農(nóng)村居民生產(chǎn)用地和生活居住用地之間的沖突難以通過(guò)促進(jìn)技術(shù)進(jìn)步、提高單位面積農(nóng)作物產(chǎn)量等方式進(jìn)行緩解，農(nóng)業(yè)收入水平將相對(duì)低于8000 m3/s防護(hù)堤方案。

綜上發(fā)現(xiàn)，建設(shè)可以防御一定標(biāo)準(zhǔn)洪水的防護(hù)堤工程，保障部分農(nóng)業(yè)生產(chǎn)用地，對(duì)于穩(wěn)定提高灘區(qū)農(nóng)業(yè)收入水平具有重要作用。但同時(shí)需要注意的是，防洪標(biāo)準(zhǔn)的提高，將會(huì)導(dǎo)致人口加快增長(zhǎng)。因此，穩(wěn)定灘區(qū)人口增速、適當(dāng)控制灘區(qū)人口規(guī)模是實(shí)現(xiàn)灘區(qū)自然-經(jīng)濟(jì)-社會(huì)協(xié)同可持續(xù)發(fā)展的另一個(gè)必要措施。

5 結(jié)論和建議

5.1 結(jié)論本文在分析黃河下游寬灘區(qū)滯洪沉沙和社會(huì)經(jīng)濟(jì)發(fā)展的基礎(chǔ)上，以自然-經(jīng)濟(jì)-社會(huì)協(xié)調(diào)發(fā)展為目標(biāo)，構(gòu)建了適用于黃河灘區(qū)的系統(tǒng)動(dòng)力學(xué)模型。在此基礎(chǔ)上，厘清了黃河下游灘區(qū)洪水治理、土地資源配置和灘區(qū)經(jīng)濟(jì)社會(huì)發(fā)展之間的聯(lián)動(dòng)關(guān)系，模擬了黃河下游灘區(qū)治理方案下未來(lái)經(jīng)濟(jì)社會(huì)系統(tǒng)發(fā)展變化情況。2016—2030年黃河下游灘區(qū)經(jīng)濟(jì)社會(huì)系統(tǒng)發(fā)展?fàn)顩r預(yù)估分析結(jié)果發(fā)現(xiàn)，無(wú)堤方案和6000 m3/s防護(hù)堤方案對(duì)灘區(qū)農(nóng)業(yè)生產(chǎn)和居民生活的保護(hù)有限，8000和10 000 m3/s標(biāo)準(zhǔn)防護(hù)堤方案能夠較好保障灘區(qū)居民的生產(chǎn)生活。但考慮到10 000 m3/s標(biāo)準(zhǔn)防護(hù)堤方案下可能存在的農(nóng)村居民生產(chǎn)用地和生活居住用地之間的矛盾，因此認(rèn)為8000 m3/s防護(hù)堤方案整體較優(yōu)，建議未來(lái)黃河下游灘區(qū)主要以8000 m3/s防護(hù)堤方案開(kāi)展灘區(qū)自然-經(jīng)濟(jì)-社會(huì)協(xié)同的可持續(xù)發(fā)展模式研究。

5.2 建議本文通過(guò)構(gòu)建的灘區(qū)自然-經(jīng)濟(jì)-社會(huì)協(xié)同發(fā)展的系統(tǒng)動(dòng)力學(xué)模型仿真模擬結(jié)果表明，對(duì)于防護(hù)堤方案，防護(hù)標(biāo)準(zhǔn)越高，在不施加其他外來(lái)政策干預(yù)情況下，灘區(qū)人口增長(zhǎng)越快，人口峰值更高，且人口峰值到來(lái)時(shí)間越晚，延續(xù)時(shí)間越長(zhǎng)。這意味著，高標(biāo)準(zhǔn)防護(hù)方案雖然提高了黃河下游灘區(qū)的防洪標(biāo)準(zhǔn)和灘區(qū)居民生產(chǎn)生活保障水平，但宜農(nóng)宜居的環(huán)境對(duì)人口的吸引力也逐漸增大，形成了對(duì)人口增加的一種變相鼓勵(lì)，這又將加劇灘區(qū)土地利用與灘區(qū)滯沙沉沙功能之間的矛盾。為解決這一矛盾，必須從黃河下游灘區(qū)所在區(qū)域全局協(xié)調(diào)發(fā)展的角度，實(shí)施更加積極的人口管理，尤其把更積極的人口政策與扶貧政策結(jié)合起來(lái)，通過(guò)灘區(qū)外新建城鎮(zhèn)、整體搬遷、人口疏解等多種方式，有序加快灘區(qū)人口外遷，通過(guò)搬遷扶貧、產(chǎn)業(yè)扶貧和轉(zhuǎn)業(yè)扶貧等多種方式，增強(qiáng)區(qū)外生活吸引力，固化灘區(qū)人口外遷成效，實(shí)現(xiàn)在防洪標(biāo)準(zhǔn)提高和防洪功能強(qiáng)化的同時(shí)，最大限度減少人口因素的干擾。

此外，本研究對(duì)不同防護(hù)方案下的仿真模擬僅到2030年，時(shí)間長(zhǎng)度仍有些不足，要實(shí)現(xiàn)流域管理與區(qū)域社會(huì)協(xié)同發(fā)展，還需要進(jìn)一步加深研究，將黃河下游灘區(qū)的土地利用與防汛安全、社會(huì)經(jīng)濟(jì)發(fā)展及生態(tài)環(huán)境保護(hù)結(jié)合起來(lái)，通過(guò)對(duì)高風(fēng)險(xiǎn)區(qū)的土地資源利用進(jìn)行合理地規(guī)劃和管理，提高灘區(qū)發(fā)展的可持續(xù)性，同時(shí)降低、避免災(zāi)害造成的破壞。通過(guò)立法形式科學(xué)劃定灘區(qū)管理邊界，明晰管理職責(zé)和范圍，通過(guò)分區(qū)分類(lèi)管理和制定土地用途管制制度，建立灘區(qū)土地利用保障機(jī)制等，實(shí)現(xiàn)灘區(qū)自然-經(jīng)濟(jì)-社會(huì)協(xié)調(diào)可持續(xù)發(fā)展。